Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Multiple Sequence Alignment”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources
PubMed1987

Progressive sequence alignment as a prerequisite to correct phylogenetic trees.

A progressive alignment method is described that utilizes the Needleman and Wunsch pairwise alignment algorithm iteratively to achieve the multiple alignment of a set of protein sequences and to construct an evolutionary tree depicting their relationship. The sequences are assumed a priori to share a common ancestor, and the trees are constructed from difference matrices derived directly from the multiple alignment. The thrust of the method involves putting more trust in the comparison of recently diverged sequences than in those evolved in the distant past. In particular, this rule is followed: "once a gap, always a gap." The method has been applied to three sets of protein sequences: 7 superoxide dismutases, 11 globins, and 9 tyrosine kinase-like sequences. Multiple alignments and phylogenetic trees for these sets of sequences were determined and compared with trees derived by conventional pairwise treatments. In several instances, the progressive method led to trees that appeared to be more in line with biological expectations than were trees obtained by more commonly used methods.

Algorithms↗
PubMed1991

Phylogenetic relationships among megabats, microbats, and primates.

We present 744 nucleotide base positions from the mitochondrial 12S rRNA gene and 236 base positions from the mitochondrial cytochrome oxidase subunit I gene for a microbat, Brachyphylla cavernarum, and a megabat, Pteropus capestratus, in phylogenetic analyses with homologous DNA sequences from Homo sapiens, Mus musculus (house mouse), and Gallus gallus (chicken). We use information on evolutionary rate differences for different types of sequence change to establish phylogenetic character weights, and we consider alternative rRNA alignment strategies in finding that this mtDNA data set clearly supports bat monophyly. This result is found despite variations in outgroup used, gap coding scheme, and order of input for DNA sequences in multiple alignment bouts. These findings are congruent with morphological characters including details of wing structure as well as cladistic analyses of amino acid sequences for three globin genes and indicate that neurological similarities between megabats and primates are due to either retention of primitive characters or to convergent evolution rather than to inheritance from a common ancestor. This finding also indicates a single origin for flight among mammals.

Amino Acid Sequence↗
PubMed1990

Calculating percent identity between protein or DNA sequences with a word processor.

Two macros, to calculate percentage identity between protein or DNA sequences using the Microsoft Word word processor, are described. The user prepares an alignment file of multiple sequences which is used by the macros to calculate number of matches, number of mismatches, total number of compared positions, and the percent identity. The macros are especially useful when alignment of multiple sequences is possible only by eye.

Algorithms↗
PubMedDate not supplied

Multiple alignment and hierarchical clustering of conserved amino acid sequences in the replication-associated proteins of plant RNA viruses.

We have used multiple alignment computer programs to align and hierarchically cluster the conserved amino acid "signature" sequences found in the replication-associated proteins of all plant RNA viruses sequenced so far. These regions, called "polymerase", "nucleotide-binding" and "N-terminal" are well conserved even between viruses which are only distantly related, and are thus very well suited for this type of analysis. Our results show that the clusterings obtained using these very short amino acid sequences are very robust to computing parameters and are surprisingly well matched with the taxonomic grouping of RNA plant viruses. The possibility of using this system as a new taxonomic criterion is discussed.

Amino Acid Sequence↗
PubMed1992

Analysis of conserved domains and sequence motifs in cellular regulatory proteins and locus control regions using new software tools for multiple alignment and visualization.

With the tremendous expansion of molecular sequence data in recent years, multiple alignment is arguably one of the two most important analytic techniques (the other being fast database searching). A number of useful approaches to this problem have previously been developed, but often they are limited to only a subset of multiple-alignment applications and cannot easily deal with the complex structural organization seen in an increasing number of sequences. For example, a single sequence may contain several domains of different evolutionary origins, and the multiplicities and relative ordering of these domains may be quite different among related sequences. Here we describe an integrated set of interactive Unix tools that combines several multiple-alignment techniques with traditional "dot-plot" visualization to provide a flexible environment for approaching complex sequence analysis problems. We apply these tools to the identification and characterization of "catalytic" domains in ras and rho/rac GTPase-activating proteins, to "Src homology" (SH2, SH3) domains in cytoplasmic signaling proteins, to repetitive sequence motifs in the alpha and beta subunits of protein prenyltransferases, and to regulatory DNA sequences in the locus control region of the beta-globin gene cluster.

Alkyl and Aryl Transferases↗
PubMed1989

Evolutionary relationship between the TonB-dependent outer membrane transport proteins: nucleotide and amino acid sequences of the Escherichia coli colicin I receptor gene.

The nucleotide sequence of the Escherichia coli colicin I receptor gene (cir) has been determined. The predicted mature protein consists of 599 amino acids and has a molecular weight of 67,169. Several previously noted characteristics of other E. coli outer membrane protein sequences were also identified in the sequence of Cir. These include an overall acidic nature, the absence of long hydrophobic stretches of amino acids, and a lack of predicted alpha-helical secondary structure. Because two classes of outer membrane proteins (the TonB-dependent transport proteins and the porins) share some structural features, protein sequences from both of these groups were aligned pairwise and scored for sequence similarity. Statistical evidence suggested that the porins were not related to the proteins in the TonB-dependent group; however, there was a significant relationship between the proteins in the TonB-dependent group. On the basis of the multiple progressive sequence alignment and the similarity scores derived from it, a tree representing evolutionary distance between five TonB-dependent outer membrane transport proteins was generated.

Amino Acid Sequence↗
PubMed1995

A multiple alignment of the capsid protein sequences of nepoviruses and comoviruses suggests a common structure.

The amino acid sequences of the regions encoding the structural proteins of eleven nepoviruses and five comoviruses, two genera of the family Comoviridae, have been aligned. The properties predicted by computer analysis (three-dimensional-3D-structure, hydrophobicity) are also correlated along this alignment, and aligned to the experimentally determined 3D structure of two comoviruses. It can thus be assumed that the 3D structure of the unique nepovirus coat protein matches that of the bipartite protomer found in the comovirus particles. In this model, the spatial locations of two amino-acid motifs characteristic of nepoviruses are in close vicinity, at the external surface of the virion. The coat proteins of nepoviruses and comoviruses may thus share a common evolutionary origin. A phylogenetic analysis was made using the multiple alignment, allowing a better understanding of the molecular relationships between these two groups of viruses.

Amino Acid Sequence↗
PubMed1992

The inference of evolutionary trees from molecular data.

1. Procedures for multiple alignment of sequence data, subsequent phylogenetic inference, and testing of the trees derived are presented. 2. The assumptions underlying different approaches and the extent to which they are valid are discussed.

Amino Acid Sequence↗
PubMed1995

Evidence for an evolutionary relationship among type-II restriction endonucleases.

Type-II restriction-modification (R-M) systems comprise two enzymes, a DNA methyltransferase (MTase) and a restriction endonuclease (ENase), each of which specifically interact with the same 4-8 bp sequence. All type-II MTases share several amino acid (aa) sequence motifs, which makes an evolutionary relatedness among these enzymes probable. The type-II ENases, in contrast, except for some homologous isoschizomers, do not share significant aa sequence similarity. Therefore, ENases in general have been considered unrelated. Here we show that in addition to the analysis of the genotype (aa sequence), a comparison of the phenotype (recognition sequence) of these enzymes can provide independent information regarding evolutionary relationships, and thereby, help to analyze the significance of weak aa sequence similarities. Multistep Monte-Carlo analyses were employed to demonstrate that the recognition sequences of those ENases, which were found to be related by a progressive multiple aa sequence alignment, are more similar to each other than would be expected by chance. This analysis supports the notion that not only type-II MTases, but also type-II ENases did not arise independently in evolution, but rather evolved from one or a few primordial DNA-modifying and DNA-cleaving enzymes, respectively.

Amino Acid Sequence↗
PubMed1994

A novel RNA-binding motif in omnipotent suppressors of translation termination, ribosomal proteins and a ribosome modification enzyme?

Using computer methods for database search, multiple alignment, protein sequence motif analysis and secondary structure prediction, a putative new RNA-binding motif was identified. The novel motif is conserved in yeast omnipotent translation termination suppressor SUP1, the related DOM34 protein and its pseudogene homologue; three groups of eukaryotic and archaeal ribosomal proteins, namely L30e, L7Ae/S6e and S12e; an uncharacterized Bacillus subtilis protein related to the L7A/S6e group; and Escherichia coli ribosomal protein modification enzyme RimK. We hypothesize that a new type of RNA-binding domain may be utilized to deliver additional activities to the ribosome.

Amino Acid Sequence↗
PubMed1991

Eukaryotic DNA polymerase amino acid sequence required for 3'----5' exonuclease activity.

We have identified an amino-proximal sequence motif, Phe-Asp-Ile-Glu-Thr, in Saccharomyces cerevisiae DNA polymerase II that is almost identical to a sequence comprising part of the 3'----5' exonuclease active site of Escherichia coli DNA polymerase I. Similar motifs were identified by amino acid sequence alignment in related, aphidicolin-sensitive DNA polymerases possessing 3'----5' proofreading exonuclease activity. Substitution of Ala for the Asp and Glu residues in the motif reduced the exonuclease activity of partially purified DNA polymerase II at least 100-fold while preserving the polymerase activity. Yeast strains expressing the exonuclease-deficient DNA polymerase II had on average about a 22-fold increase in spontaneous mutation rate, consistent with a presumed proofreading role in vivo. In multiple amino acid sequence alignments of this and two other conserved motifs described previously, five residues of the 3'----5' exonuclease active site of E. coli DNA polymerase I appeared to be invariant in aphidicolin-sensitive DNA polymerases known to possess 3'----5' proofreading exonuclease activity. None of these residues, however, appeared to be identifiable in the catalytic subunits of human, yeast, or Drosophila alpha DNA polymerases.

Amino Acid Sequence↗
PubMed1987

Selection of circularization sites in a group I IVS RNA requires multiple alignments of an internal template-like sequence.

Circularization and reverse circularization of the Tetrahymena thermophila rRNA intervening sequence resemble the first and second steps in splicing, respectively. However, site-specific base substitutions show that different nucleotides are involved in selection of the 5' splice site and the circularization sites. Furthermore, a substitution at the major circularization site that prevents circularization can be suppressed by second substitutions at two different nucleotide positions. A model is proposed in which adjacent and overlapping sequences can function as a binding site, forming a short duplex with the sequence at the circularization site and thus directing circularization and reverse circularization. Because the 5' exon-binding site and three potential circularization binding sites fall within a contiguous eight nucleotide region, this sequence may translocate relative to the catalytic core of the ribozyme in a template-like manner.

Animals↗
PubMed2026

KCFtools: rapid alignment-free method for introgression screening and GWAS using k-mer profiles.

MOTIVATION: In the era of multiple genome references, researchers often align sequencing reads against distinct assemblies or even multiple references simultaneously. This enables applications such as the detection of introgressed segments or highly variable genomic regions, which are especially prevalent in large-genome crop species such as lettuce or wheat. However, these applications come at the cost of increased computational burden, inconsistencies in mapping methods, and reduced reproducibility across studies. To address these limitations, we developed KCFtools, a Java-based toolkit that identifies the presence and absence of k-mers in nonoverlapping genomic or transcriptomic windows by comparing query and reference genomes. This alignment-free approach enables the efficient computation of an identity score for each window, thereby facilitating robust detection of introgressed or variable regions across genomes. RESULTS: We systematically evaluated the performance and accuracy of the k-mer-based method implemented in KCFtools, benchmarking it against conventional single nucleotide variation-based introgression detection pipelines. Our results demonstrate that KCFtools effectively captures introgressed segments and structurally diverse regions, even in species with fragmented or highly divergent reference genomes. In addition, we extended KCFtools to generate genotype matrices from k-mer variation tables. These matrices are compatible with genome-wide association studies software and allow the identification of loci associated with phenotypic traits. We showcase the utility of this approach by detecting known and novel associations for downy mildew resistance in lettuce, underscoring the pipeline's potential for high-resolution, reference-agnostic population genetic analysis. AVAILABILITY AND IMPLEMENTATION: https://github.com/sivasubramanics/kcftools.

Software↗
PubMed1990

Herpesviral deoxythymidine kinases contain a site analogous to the phosphoryl-binding arginine-rich region of porcine adenylate kinase; comparison of secondary structure predictions and conservation.

Twelve herpesviral deoxythymidine kinases were examined for regions of sequence similarity by multiple alignment. Six highly conserved sites were observed. Site 1 corresponded to a glycine-rich loop that forms part of the ATP-binding pocket in porcine adenylate kinase (PAK), and site 5 corresponded to a region in PAK, located on one lobe of the cleft, that contains arginine residues that bind substrate phosphoryl groups. Site 3, consisting of the motif -DRH-, is thought to be involved in thymine/deoxythymidine recognition; site 4, which is nearby, probably participates in this function as well. The functions of sites 2 and 6 have not been identified. Secondary structure predictions were made by the Garnier method and averaged for each position in the multiple alignment. The structure predicted for all six sites was typically a short flexible region (turn or coil) at or adjacent to the site, flanked by rigid structures (helix or sheet) on either side.

Adenylate Kinase↗
PubMed1992

Novel GACG-hairpin pair motif in the 5' untranslated region of type C retroviruses related to murine leukemia virus.

We searched for the presence of common RNA structural motifs in mammalian type C retroviruses related to murine leukemia viruses and the closely related avian spleen necrosis virus. A novel motif consisting of a pair of hairpins, called hairpin pair motif, was detected in the 5' untranslated regions of the genomes of these retroviruses. A combination of computational analyses that included the assessment of phylogenetic sequence conservation by multiple alignment, the search for regions with unusual RNA folding properties, and the analysis of RNA secondary structure by suboptimal free-energy calculations highlighted the significance of this hairpin pair motif. The hairpin pair motif encompasses 70 to 80 nucleotides between the splice donor site and the gag translational initiation codon of these viruses. The motif is composed of two adjacent hairpins both with a perfectly conserved GACG tetraloop. We propose that the novel GACG-hairpin pair motif described here constitutes an essential component of the regulatory machinery in these type C retroviruses.

Base Sequence↗
PubMed1992

Expansion of the mammalian 3 beta-hydroxysteroid dehydrogenase/plant dihydroflavonol reductase superfamily to include a bacterial cholesterol dehydrogenase, a bacterial UDP-galactose-4-epimerase, and open reading frames in vaccinia virus and fish lymphocystis disease virus.

Mammalian 3 beta-hydroxysteroid dehydrogenase and plant dihydroflavonol reductases are descended from a common ancestor. Here we present evidence that Nocardia cholesterol dehydrogenase, E. coli UDP-galactose-4 epimerase, and open reading frames in vaccinia virus and fish lymphocystis disease virus are homologous to 3 beta-hydroxysteroid dehydrogenase and dihydroflavonol reductase. Analysis of a multiple alignment of these sequences indicates that viral ORFs are most closely related to the mammalian 3 beta-hydroxysteroid dehydrogenases. The ancestral protein of this superfamily is likely to be one that metabolized sugar nucleotides. The sequence similarity between 3 beta-hydroxysteroid dehydrogenase and the viral ORFs is sufficient to suggest that these ORFs have an activity that is similar to 3 beta-hydroxysteroid dehydrogenase or cholesterol dehydrogenase, although the putative substrates are not yet known.

3-Hydroxysteroid Dehydrogenases↗