PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Biological sequence analysis”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 91 records · Page 5Linked to original sources
PubMed2001

Manuscript evolution.

Frequently, letters, words and sentences are used in undergraduate textbooks and the popular press as an analogy for the coding, transfer and corruption of information in DNA. We discuss here how the converse can be exploited, by using programs designed for biological analysis of sequence evolution to uncover the relationships between different manuscript versions of a text. We point out similarities between the evolution of DNA and the evolution of texts.

Evolution, Molecular↗
PubMed2001

Manuscript evolution.

Frequently, letters, words and sentences are used in undergraduate textbooks and the popular press as an analogy for the coding, transfer and corruption of information in DNA. We discuss here how the converse can be exploited, by using programs designed for biological analysis of sequence evolution to uncover the relationships between different manuscript versions of a text. We point out similarities between the evolution of DNA and the evolution of texts.

DNA↗
PubMed2000

Cloning and characterization of a Bacillus thuringiensis serovar higo gene encoding a novel class of the delta-endotoxin protein, Cry27A, specifically active on the Anopheles mosquito.

A novel gene encoding a 98-kDa mosquitocidal delta-endotoxin protein, designated Cry27A, was cloned from a Bacillus thuringiensis serovar higo strain. The Cry27A protein contained the five sequence blocks of amino acids commonly conserved in most B. thuringiensis Cry proteins. Relatively high homologies, ranging from 43.0% to 84.4%, existed between the Cry27A protein and several established classes of mosquitocidal Cry proteins (Cry4A, Cry10A, Cry19A, Cry19B, and Cry20A) in the sequence of 51 N-terminal amino acids. The complete sequence of this protein, however, showed low levels (<40%) of amino acid identity to those of the known Cry proteins. Although the expression level of the cry27A gene was low in the transformants under the control of its own promoter, the use of the cyt1A promoter resulted in high-level expression of the gene, leading to the formation of inclusions. The expressed Cry27A protein showed larvicidal activity highly specific for Anopheles stephensi, but lacked the toxicity against Culex pipiens molestus and Aedes aegypti. The results suggest that the Cry27A protein is responsible for the Anopheles-preferential toxicity of the B. thuringiensis serovar higo strain.

Amino Acid Sequence↗
PubMed2000

Whole genomes: the foundation of new biology and medicine.

Our genomic DNA sequence provides a unique glimpse of the provenance and evolution of our species, the migration of peoples, and the causation of disease. Understanding the genome may help resolve previously unanswerable questions, including perhaps which human characteristics are innate or acquired. Such an understanding will make it possible to study how genomic DNA sequence varies among populations and among individuals, including the role of such variation in the pathogenesis of important illnesses and responses to pharmaceuticals. The study of the genome and the associated proteomics of free-living organisms will eventually make it possible to localize and annotate every human gene, as well as the regulatory elements that control the timing, organ-site specificity, extent of gene expression, protein levels, and post-translational modifications. For any given physiological process, we will have a new paradigm for addressing its evolution, development, function, and mechanism.

Animals↗
PubMed1999

Detection of new DNA polymerase genes of known and potentially novel herpesviruses by PCR with degenerate and deoxyinosine-substituted primers.

A consensus primer PCR approach was used to (i) investigate the presence of herpesviruses in wild and zoo equids (zebra, wild ass, tapir) and to (ii) study the genetic relationship of the herpesvirus of pigeons (columbid herpesvirus 1) to other herpesvirus species. The PCR assay, based on degenerate primers targeting highly conserved regions of the DNA polymerase gene of herpesviruses, was modified by using a mixture of degenerate and deoxyinosine-substituted primers. The applicability of the modification was validated by amplification of published DNA polymerase genes of 16 herpesvirus species and of the previously uncharacterized DNA polymerase genes of equine herpesvirus 3 (EHV-3) and equine herpesvirus 5 (EHV-5). The modified assay was then used for partial amplification of the polymerase of columbid herpesvirus 1 which is presently classified as a beta-herpesvirus based on biological criteria. Sequence analysis of amplicons obtained from four different viral strains revealed a close relationship of columbid herpesvirus 1 to members of the subfamily Alphaherpesvirinae, especially to Marek's disease herpesvirus. This was confirmed by characterization of additional 1.6kb of the columbid herpesvirus 1 polymerase. Consensus PCR analysis of blood samples from zebras, a wild ass and a tapir revealed amplicons showing high percentages ( > 50%) of sequence identity to DNA polymerases of gamma-herpesviruses. In particular, the zebra and the wild ass sequence were closely related to each other and to the polymerases of the equine gamma-herpesviruses EHV-2 and EHV-5 with sequence identities of > 80%. This is a first indication that novel gamma-herpesviruses are present in wild and zoo equids.

Amino Acid Sequence↗
PubMed1985

Activation of an N-ras gene in acute myeloblastic leukemia through somatic mutation in the first exon.

A transforming N-ras gene has been cloned from acute myeloblastic leukemia bone marrow cells, in parallel with the N-ras gene derived from fibroblasts of the same patient. N-ras derived from fibroblasts lacked focus-forming activity in NIH/3T3 cells, indicating that gene activation in the leukemia cells must have occurred by a somatic event. Construction of chimeric molecules between the transforming and the normal N-ras genes and subsequent biological and sequence analysis of these constructs revealed that the transforming gene was altered by a point mutation changing amino acid 12 of the N-ras protein from glycine to aspartic acid.

Alleles↗
PubMed2002

Conceptual data modelling for bioinformatics.

Current research in the biosciences depends heavily on the effective exploitation of huge amounts of data. These are in disparate formats, remotely dispersed, and based on the different vocabularies of various disciplines. Furthermore, data are often stored or distributed using formats that leave implicit many important features relating to the structure and semantics of the data. Conceptual data modelling involves the development of implementation-independent models that capture and make explicit the principal structural properties of data. Entities such as a biopolymer or a reaction, and their relations, eg catalyses, can be formalised using a conceptual data model. Conceptual models are implementation-independent and can be transformed in systematic ways for implementation using different platforms, eg traditional database management systems. This paper describes the basics of the most widely used conceptual modelling notations, the ER (entity-relationship) model and the class diagrams of the UML (unified modelling language), and illustrates their use through several examples from bioinformatics. In particular, models are presented for protein structures and motifs, and for genomic sequences.

Computational Biology↗
PubMed2025

vcfgl: a flexible genotype likelihood simulator for VCF/BCF files.

MOTIVATION: Accurate quantification of genotype uncertainty is pivotal in ensuring the reliability of genetic inferences drawn from NGS data. Genotype uncertainty is typically modeled using Genotype Likelihoods (GLs), which can help propagate measures of statistical uncertainty in base calls to downstream analyses. However, the effects of errors and biases in the estimation of GLs, introduced by biases in the original base call quality scores or the discretization of quality scores, as well as the choice of the GL model, remain under-explored. RESULTS: We present vcfgl, a versatile tool for simulating genotype likelihoods associated with simulated read data. It offers a framework for researchers to simulate and investigate the uncertainties and biases associated with the quantification of uncertainty, thereby facilitating a deeper understanding of their impacts on downstream analytical methods. Through simulations, we demonstrate the utility of vcfgl in benchmarking GL-based methods. The program can calculate GLs using various widely used genotype likelihood models and can simulate the errors in quality scores using a Beta distribution. It is compatible with modern simulators such as msprime and SLiM, and can output data in pileup, Variant Call Format (VCF)/BCF, and genomic VCF file formats, supporting a wide range of applications. The vcfgl program is freely available as an efficient and user-friendly software written in C/C++. AVAILABILITY AND IMPLEMENTATION: vcfgl is freely available at https://github.com/isinaltinkaya/vcfgl.

Software↗
PubMed2025

Locality-aware pooling enhances protein language model performance across varied applications.

MOTIVATION: Protein language models (PLMs) are amongst the most exciting recent advances for characterizing protein sequences, and have enabled a diverse set of applications, including structure determination, functional property prediction, and mutation impact assessment, all from single protein sequences alone. State-of-the-art PLMs leverage transformer architectures originally developed for natural language processing, and are pre-trained on large protein databases to generate contextualized representations of individual amino acids. To harness the power of these PLMs to predict protein-level properties, these per-residue embeddings are typically "pooled" to fixed-size vectors that are further utilized in downstream prediction networks. Common pooling strategies include Cls-Pooling and Avg-Pooling, but neither of these approaches can capture the local substructures and long-range interactions observed in proteins. RESULTS: We propose the use of attention pooling, which can naturally capture these important features of proteins. To make the expensive attention operator (quadratic in the length of the input protein) feasible in practice, we introduce bag-of-mer pooling, or BoM-Pooling, a locality-aware hierarchical pooling technique that combines windowed average pooling with attention pooling. We empirically demonstrate that both full attention pooling and BoM-Pooling outperform previous pooling strategies on three important, diverse tasks: (i) predicting the activities of two proteins as they are varied; (ii) detecting remote homologs; and (iii) predicting signaling protein interactions with peptides. Overall, our work highlights the advantages of biologically inspired pooling techniques in protein sequence modeling and is a step toward more effective adaptations of language models in biological settings. AVAILABILITY AND IMPLEMENTATION: https://github.com/Singh-Lab/bom-pooling.

Natural Language Processing↗
PubMed2025

Fast and flexible minimizer digestion with digest.

SUMMARY: Minimizer digestion is an increasingly common component of bioinformatics tools, including tools for de Bruijn graph assembly and sequence classification. We describe a new open source tool and library to facilitate efficient digestion of genomic sequences. It can produce digests based on the related ideas of minimizers, modimizers or syncmers. Digest uses efficient data structures, scales well to many threads, and produces digests with expected spacings between digested elements. AVAILABILITY AND IMPLEMENTATION: Digest is implemented in C++17 with a Python API, and is available open-source at https://github.com/VeryAmazed/digest. The python library is available on Bioconda. Rust bindings are available as a public crate at https://crates.io/crates/digest-rs.

Software↗
PubMed2026

seq2ribo: structure-aware integration of machine learning and simulation to predict ribosome location profiles from RNA sequences.

MOTIVATION: Ribosome dynamics are vital in the process of protein expression. Current methods rely on ribosome profiling (Ribo-seq), RNA-seq profiles, and full genomic context. This restricts their use in de novo sequence design, like messenger RNA (mRNA) vaccines. Simulation-only approaches like the Totally Asymmetric Simple Exclusion Process (TASEP) oversimplify translation by focusing solely on codon elongation times. RESULTS: We present seq2ribo, a hybrid simulation and machine learning framework that predicts ribosome A-site locations using only an mRNA sequence as input. Our method first employs a novel structure-aware TASEP (sTASEP), which models translation using a comprehensive set of fitted parameters that include codon wait times and structural features, such as local angles, base-pairing, and discrete positional buckets. The ribosome locations generated by sTASEP are then processed by a polisher model, which learns to refine the simulated ribosome distributions. seq2ribo provides high-fidelity predictions of ribosome locations across diverse cell types (iPSC, HEK293, LCL, and RPE-1), significantly outperforming baselines. seq2ribo is the first method to achieve meaningful positional correlation with observed ribosome profiles from sequence alone, reaching transcript-level Pearson correlations up to 0.920 and within-transcript shape correlations up to 0.186, where all baselines yield near-zero values on these metrics. seq2ribo also reduces elementwise error by up to 37.7% relative to the sequence-only Translatomer baseline. By adding a task-specific head, seq2ribo achieves Pearson correlations up to 0.732 with experimental translation efficiency (TE) across several cell lines, and up to 0.903 with measured protein expression. By operating from sequence alone, seq2ribo provides a new tool for synthetic biology, enabling the rational design and optimization of mRNA sequences without the need for expression-level data or genomic context. AVAILABILITY: seq2ribo is available at https://github.com/Kingsford-Group/seq2ribo.

Machine Learning↗
PubMed2026

GraphyloVar: predicting the impact of non-coding variants using a multi-species sequence model.

MOTIVATION: Understanding the functional impact of genetic variants is a key problem for precision medicine. Tools like CADD, PhyloP, and PhastCons are useful, but they often look at each position in the genome in isolation. This means they can miss important information from the evolutionary history that connects different species. In this paper, we extend our previous model, Graphylo, to predict the effects of variants. Our new model, GraphyloVar, is built to directly utilize the phylogenetic tree that relates the species. RESULTS: GraphyloVar is a deep learning model that considers both DNA sequence and evolutionary patterns from many species. It uses two main components: Graph Convolutional Networks (GCNs) to process the phylogenetic tree, and Transformer encoders to extract features from the DNA sequences. Pre-trained to predict population-level allele frequencies on the TOPMed whole-genome sequencing cohort, GraphyloVar achieves an AUROC of 0.6246 zero-shot on &#x223c;149M held-out variants, and an ensemble with CADD reaches 0.6442 (+0.020, P<10-15). Fine-tuned GraphyloVar achieves the highest AUROC across all 13 MPRA benchmark datasets. By integrating deep learning with explicit phylogenetic input, GraphyloVar offers a powerful and complementary approach to variant effect prediction that utilizes the full evolutionary history from many species to better identify and prioritize important non-coding variants. AVAILABILITY AND IMPLEMENTATION: Code and datasets are available at https://github.com/DongjoonLim/GraphyloVar under DOI: 10.5281/zenodo.20616818.

Phylogeny↗
PubMed2002

In vitro analysis of human immunodeficiency virus type 1 resistance to nevirapine and fitness determination of resistant variants.

Nevirapine-resistant variants were generated by serial passages in MT-2 cells in the presence of increasing drug concentrations. In passage 5, mutations V106A, Y181C and G190A were detected in the global population, associated with a 100-fold susceptibility decrease. Sequence analysis of biological clones obtained from passage 5 and subsequent passages showed that single mutants, detected in first passages, were progressively replaced in passage 15 by double mutants, correlating with a 500-fold increase in phenotypic resistance. Fitness determination of single mutants confirmed that, in the presence of nevirapine, every variant was more fit than wild-type with a fitness order Y181C>V106A>G190A>wild-type. Unexpectedly, in the absence of the drug, the Y181C resistant mutant was more fit than wild-type, with a fitness gradient Y181C>wild-type >G106A>or=V190A. Using a molecular clone in which the Y181C mutation was introduced by in vitro mutagenesis, the greater fitness of the Y181C mutant was confirmed in new competition cultures. These data exemplify the role of resistance mutations on virus phenotype but also on virus evolution leading, occasionally, to resistant variants fitter than the wild-type in the absence of the drug.

Adaptation, Physiological↗
PubMed1993

Identification of two expressed flagellin genes in the insect pathogen Bacillus thuringiensis subsp. alesti.

Flagellin from Bacillus thuringiensis subspecies alesti strain Bt75 was isolated both from the culture medium and from flagella. Two protein forms with molecular masses close to 32 kDa were obtained from flagella; one form was identical to the flagellin purified from the culture medium. The N-terminal amino acid sequences were identical for both forms. Two genes coding for flagellin have been identified in B. thuringiensis subsp. alesti. The flaB gene was cloned and sequenced in its entire length. The clone containing the flaA gene was incomplete. Both genes were expressed in the mid-exponential growth phase. The flaB gene was flanked by long (355 bp) direct repeats protruding into the coding region in both the N- and C-terminal parts of the gene. DNA sequences related to the flaB gene were found in most other B. thuringiensis subspecies, and in two of them, subsp. kurstaki and subsp. entomocidus, such sequences were present in multiple copies.

Amino Acid Sequence↗
PubMed1995

Beta 3-adrenoceptors and intestinal motility.

Early substantial evidence of the low susceptibility to beta-adrenoceptor antagonists of non alpha-adrenergic responses reducing gut motility and tone was reluctantly accepted as indicating a third beta-receptor subtype different from the beta 1 and beta 2. This applied likewise to lipolysis until new selective "lipolytic" beta-agonists poorly effective at established beta-receptors were introduced. Shortly afterwards these "lipolytic" as well as certain newer and even more selective beta-adrenoceptor agonists were shown to be potent inhibitors of intestinal motility. The latter are the "gut-specific" phenylethanolaminotetralins whose availability as pure isomers attested to the stringent stereochemical requirements for selectivity at non-beta 1, non-beta 2 beta-adrenoceptors. Acceptance of the functionally based concept of a beta 3-adrenoceptor was boosted on structural grounds by molecular biology studies. Sequence analysis indicated the existence in humans and rodents of genes coding for a third subtype of beta-receptor that, when expressed in transfected heterologous cells, had a pharmacological profile distinct from the previously established subtypes. Finally, aryloxypropanolaminotetralins have been prepared as the first selective antagonists of beta 3-adrenoceptors, thus providing unambiguous conclusive evidence of the distinctive functional features of those abundant in the rat colon. The therapeutic potential in gastroenterology of the newer compounds targetable on the beta 3-adrenoceptor is suggested by their potent intestinal action in vivo in animal models without any of the cardiovascular or other unwanted effects of conventional beta-adrenoceptor agonists and antagonists, and by the clinically confirmed importance of beta-adrenergic control of motor function throughout the alimentary canal. However, open questions include the incidence of species-related differences in beta 3-adrenoceptors, and as yet there are no data on gastrointestinal functions in humans under the influence of drugs designed to act selectively at these receptors.

Adrenergic beta-Agonists↗
PubMed2002

Geomicrobiology: how molecular-scale interactions underpin biogeochemical systems.

Microorganisms populate every habitable environment on Earth and, through their metabolic activity, affect the chemistry and physical properties of their surroundings. They have done this for billions of years. Over the past decade, genetic, biochemical, and genomic approaches have allowed us to document the diversity of microbial life in geologic systems without cultivation, as well as to begin to elucidate their function. With expansion of culture-independent analyses of microbial communities, it will be possible to quantify gene activity at the species level. Genome-enabled biogeochemical modeling may provide an opportunity to determine how communities function, and how they shape and are shaped by their environments.

Animals↗