Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Biological sequence analysis”

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 109 records · Page 6Linked to original sources
PubMed1987

Propeptide of human protein C is necessary for gamma-carboxylation.

Protein C is one of a family of vitamin K dependent proteins, including blood coagulation factors and bone proteins, that contains gamma-carboxyglutamic acid. Sequence analysis of the cDNAs for these proteins has revealed the presence of a prepro leader sequence that contains a pre sequence or hydrophobic signal sequence and a propeptide containing a number of highly conserved amino acids. The pre region is removed from the growing polypeptide chain by signal peptidase, while the pro region is subsequently removed from the protein prior to secretion. In the present study, deletion mutants have been constructed in the propeptide region of the cDNA for human protein C, and the cDNAs were then expressed in mammalian cell culture. These deletions included the removal of 4, 9, 12, 15, 16, or 17 amino acids comprising the carboxyl end of the leader sequence of 42 amino acids. The mutant proteins were then examined by Western blotting, barium citrate adsorption and precipitation, amino acid sequence analysis, and biological activity and compared with the native protein present in normal plasma. These experiments have shown that protein C is readily synthesized in mammalian cell cultures, processed, and secreted as a two-chain molecule with biological activity. Furthermore, the pre portion or signal sequence in human protein C is 18 amino acids in length, and the pro portion of the leader sequence is 24 amino acids in length. Also, during biosynthesis and secretion, the amino-terminal region of the propeptide (residues from about -12 through -17) is important for gamma-carboxylation of protein C, while the present data and those of others indicate that the carboxyl-terminal portion of the propeptide (residues -1 through -4) is important for the removal of the pro leader sequence by proteolytic processing.

Amino Acid Sequence↗
PubMed2004

Pseudo-periodic partitions of biological sequences.

MOTIVATION: Algorithm development for finding typical patterns in sequences, especially multiple pseudo-repeats (pseudo-periodic regions), is at the core of many problems arising in biological sequence and structure analysis. In fact, one of the most significant features of biological sequences is their high quasi-repetitiveness. Variation in the quasi-repetitiveness of genomic and proteomic texts demonstrates the presence and density of different biologically important information. It is very important to develop sensitive automatic computational methods for the identification of pseudo-periodic regions of sequences through which we can infer, describe and understand biological properties, and seek precise molecular details of biological structures, dynamics, interactions and evolution. RESULTS: We develop a novel, powerful computational tool for partitioning a sequence to pseudo-periodic regions. The pseudo-periodic partition is defined as a partition, which intuitively has the minimal bias to some perfect-periodic partition of the sequence based on the evolutionary distance. We devise a quadratic time and space algorithm for detecting a pseudo-periodic partition for a given sequence, which actually corresponds to the shortest path in the main diagonal of the directed (acyclic) weighted graph constructed by the Smith-Waterman self-alignment of the sequence. We use several typical examples to demonstrate the utilization of our algorithm and software system in detecting functional or structural domains and regions of proteins. A big advantage of our software program is that there is a parameter, the granularity factor, associated with it and we can freely choose a biological sequence family as a training set to determine the best parameter. In general, we choose all repeats (including many pseudo-repeats) in the SWISS-PROT amino acid sequence database as a typical training set. We show that the granularity factor is 0.52 and the average agreement accuracy of pseudo-periodic partitions, detected by our software for all pseudo-repeats in the SWISS-PROT database, is as high as 97.6%.

Algorithms↗
PubMed1994

The genetic data environment an expandable GUI for multiple sequence analysis.

An X-Windows-based graphic user interface is presented which allows the seamless integration of numerous existing biomolecular programs into a single analysis environment. This environment is based on a core multiple sequence editor that is linked to external programs by a user-expandable menu system and is supported on Sun and DEC workstations. There is no limitation to the number of external functions that can be linked to the interface. The length and number of sequences that can be handled are limited only by the size of virtual memory present on the workstation. The sequence data itself is used as the reference point from which analysis is done, and scalable graphic views are supported. It is suggested that future software development utilizing this expandable, user-defined menu system and the I/O linkage of external programs will allow biologists to easily integrate expertise from disparate fields into a single environment.

Algorithms↗
PubMed2006

Proteomics of xenografted human breast cancer indicates novel targets related to tamoxifen resistance.

Tamoxifen is the most frequently used drug for hormone therapy of breast cancer patients, even though a high percentage of women are (or become) refractory to this treatment. The proteins involved in tamoxifen resistance of breast tumor cells as well as the mechanisms by which they interact, are still unknown. Some years ago, we established the xenograft breast tumor 3366, sensitive to tamoxifen and the 3366/TAM, resistant to tamoxifen, derived after two years of in vivo passages of the parental 3366 under tamoxifen treatment. Here, we compare the protein expression levels of both xenografts. 2-DE of proteins from total cell extracts showed very high reproducibility among tumors from each group (tamoxifen sensitive and tamoxifen resistant). The heuristic clustering analysis of these gels pooled them correctly in both groups. Twelve proteins were found up-regulated in the tamoxifen-resistant line, while nine were down-regulated. The proteins differentially expressed were identified by MS and sequence database analysis. Biological functions of these proteins are related to cell-cell adhesion and interaction, signal transduction, DNA and protein synthesis machinery, mitochondrial respiratory chain, oxidative stress processes and apoptosis. Three of the identified proteins (ALG-2 interacting protein and two GDP-dissociation inhibitors) could be directly involved in the resistance phenomenon.

Amino Acid Sequence↗
PubMed2004

Identification of a novel apolipoprotein, ApoN, in ovarian follicular fluid.

A novel apolipoprotein, designated ApoN, has been identified in bovine ovarian follicular fluid using chromatographic purification methods, amino acid sequence analysis, molecular biology, and bioinformatics. The apolipoprotein is a hydrophobic 12-kDa protein processed from the C terminus of a 29-kDa precursor expressed in a number of tissues, including the ovary, testis, the anterior chamber of the eye, skeletal muscle, uterus, and liver. Bovine, porcine, and murine ApoN display significant homology at the amino acid level across the entire precursor sequence. Surprisingly, there appears to be no orthologous protein in the human, although an APON-like pseudogene is found on chromosome 12. The N-terminal fragment of the ApoN precursor shows significant homology with the N-terminal sequence of the precursor of the cholesterol transport regulatory protein ApoF, but the corresponding C-terminal sequences of ApoN and ApoF possess no homology. ApoN is present in the high-density lipoprotein fraction of bovine serum and both the high-density lipoprotein and low-density lipoprotein fractions of bovine follicular fluid and is found in several tissues that are associated with local immunological privilege. These data suggest that ApoN may play a role in steroidogenesis and/or immunoregulation in the gonads of nonhuman species, as well as similar roles in other tissues.

Amino Acid Sequence↗
PubMed2026

Comparative essentialome analysis of six Pectobacteriaceae strains using the TNSEEK pipeline identifies conserved and strain-specific fitness determinants.

Transposon sequencing (Tn-seq) is a powerful technique for defining the essential genes required for bacterial survival. However, gene essentiality can vary significantly across taxonomic levels, and comparing large Tn-seq datasets from multiple strains presents considerable analytical challenges. To address this, we developed TNSEEK, a fully automated bioinformatics pipeline for the systematic and comparative analysis of Tn-seq experiments. We applied TNSEEK to analyse newly generated data for six soft rot Pectobacteriaceae strains, encompassing species from the Dickeya and Pectobacterium genera, grown in a rich medium. This approach identified a core essentialome of 225 genes, primarily involved in fundamental cellular maintenance, conserved across all 6 strains, a set comparable in size to that of the neighbouring Enterobacteriaceae family. Only a few genus-specific essential genes were found, highlighting interesting distinct metabolic capabilities between Dickeya and Pectobacterium genera. In striking contrast, we discovered a large variable essentialome comprising 181 strain-specific genes, many of which are of unknown function. A portion of these strain-specific essential genes are components of defence systems and prophage genomic regions. The unexpected essentiality of selected components of these modules is consistent with cellular dependency on cognate toxic, restriction or immunity functions encoded by defence-associated loci under the tested growth condition. Furthermore, a comparison with the Escherichia coli essentialome demonstrates that discrepancies in gene essentiality can often be attributed to differences in growth conditions, particularly temperature, as well as variations in genetic redundancy. In conclusion, the TNSEEK pipeline provides a reproducible framework for comparative analysis of mariner/Himar1 Tn-seq datasets across multiple strains.

Pectobacterium↗
PubMed2005

Characterization of a naturally occurring recombinant isolate of Grapevine fanleaf virus.

The naturally occurring Grapevine fanleaf virus (GFLV) recombinant isolate A17b was recovered from its grapevine host by sap inoculation and serial passages onto Gomphrena globosa, a pseudo local lesion herbaceous host, and Chenopodium quinoa, a systemic herbaceous host, to characterize some of its biological properties. Sequence analysis of the CP gene, in which a recombinational event was previously detected, demonstrated the genetic stability of recombinant isolate A17b over a 5-year period in its natural host as well as in C. quinoa. Also, recombinant isolate A17b was graft transmissible, as shown by an in vitro heterologous approach, and transmitted by the nematode Xiphinema index as readily as nonrecombinant GFLV isolates. Furthermore, despite a lower pathogenicity on Chenopodium amaranticolor, recombinant isolate A17b had a similar host range and induced similar symptoms in type and severity to nonrecombinant GFLV isolates. Interestingly, the use of infectious chimeric RNA2 transcripts in combination to RNA1 transcripts of GFLV strain F13 suggested no implication of the recombination event in the CP gene of isolate A17b in the reduced pathogenicity on C. amaranticolor. Altogether, recombinant isolate A17b had similar biological properties to GFLV nonrecombinant isolates.

Base Sequence↗
PubMed2025

NanoASV: a snakemake workflow for reproducible field-based Nanopore full-length 16S metabarcoding amplicon data analysis.

SUMMARY: NanoASV is a conda environment and snakemake-based workflow using state-of-the-art bioinformatics software to process full-length SSU rRNA (16S/18S) amplicons acquired with Oxford Nanopore Sequencing technology. Its strength lies in reproducibility, portability, and the possibility to run offline, allowing in-field analysis. It can be installed on the Nanopore MK1C sequencing device and process data locally. AVAILABILITY AND IMPLEMENTATION: Source code and documentation are freely available at https://github.com/ImagoXV/NanoASV and Zenodo archive at https://doi.org/10.5281/zenodo.14730742.

Software↗
PubMed2025

Global Profiling and Analysis of 5' Monophosphorylated mRNA Decay Intermediates.

During RNA turnover, the action of endo- and exo-ribonucleases can yield RNA decay intermediates with specific 5' ends. These RNA decay intermediates have been demonstrated to be the outcome of decapping, microRNA-directed endo-cleavage, or the protected fragments of ribosomes and exon-junction complexes. Therefore, global analysis of RNA decay intermediates can facilitate studies of many RNA decay pathways. In this chapter, we describe a high-throughput sequencing protocol named parallel analysis of RNA ends (PARE), which allows genome-wide profiling of 5' monophosphorylated mRNA decay intermediates from plants or other eukaryotes. Also, we present the tools and scripts necessary for the proper analysis of RNA degradome data obtained from the PARE method. Details and modifications of library construction procedures and bioinformatic analyses to optimize sequencing quality and cope with emerging sequencing platforms and findings are highlighted.

RNA Stability↗
PubMed2026

Nested co-expression network analysis identifies compact gene clusters in a black box.

MOTIVATION: Digital analysis of biological systems requires methods capable of identifying both broad and nested gene modules reflecting complex biological processes. Existing transcriptomic methods often miss compact gene sets corresponding to subprocesses in specialized cell types, limiting insights into functional heterogeneity. RESULTS: We present Nested-WGCNA, a two-stage unsupervised network analysis algorithm designed to identify coarse-grained and fine-grained gene modules. Applied to bulk RNA-Seq data, Nested-WGCNA reveals stable modules reproducible across datasets. When validated against scRNA-Seq data, these modules correspond to both major and minor immune cell subtypes. Application to immunotherapy response datasets uncovers predictive and prognostic biomarkers, highlighting its utility in treatment stratification and biomarker discovery. AVAILABILITY: The NestedWGCNA source code and analysis pipeline are available on GitHub (https://github.com/ilyada/NestedWGCNA) and archived on Zenodo (https://doi.org/10.5281/zenodo.18959244).

Algorithms↗
PubMed2001

Toxicity analysis of N- and C-terminus-deleted vegetative insecticidal protein from Bacillus thuringiensis.

A vegetative insecticidal protein (VIP)-encoding gene from a local isolate of Bacillus thuringiensis has been cloned, sequenced, and expressed in Escherichia coli. The expressed protein shows insecticidal activity against several lepidopteran pests but is ineffective against Agrotis ipsilon. Comparison of the amino acid sequence with those of reported VIPs revealed a few differences. Analysis of insecticidal activity with N- and C-terminus deletion mutants suggests a differential mode of action of VIP against different pests.

Animals↗
PubMed2026

Likelihood-based optimization enables accurate copy number estimation for paralogous genes using exome data.

MOTIVATION: Exome sequencing is widely used for genetic studies; however, accurate detection of copy number variants (CNV) in paralogous genes is challenging due to short-read mapping ambiguity and extensive copy-number variation. The human genome contains several hundred paralogous genes, many of which are known to harbor disease-associated CNVs. Existing exome CNV callers are primarily designed for rare CNV detection in uniquely mappable regions and are not well-suited for paralogous genes. METHODS: We describe a computational method (EdgeCopy) for copy number profiling of paralogous genes using whole-exome sequence data. EdgeCopy aggregates reads mapped to all copies of paralogous genes and relates observed read depth to copy number for multiple exome samples using an approximate composite likelihood function. The likelihood function is optimized using numerical optimization to obtain gene-level fractional copy number estimates that are discretized and refined using a Hidden Markov Model to obtain exon-level copy number estimates. RESULTS: Benchmarking of Edgecopy using experimental copy number data showed high concordance (mean = 0.973) for six disease-associated paralogous genes. We evaluated performance using whole-exome data from approximately 2400 samples across five continental populations from the 1000 Genomes Project. EdgeCopy shows robust concordance with whole-genome sequencing based estimates (0.974-0.982) across populations and 130 paralogous genes spanning a wide range of copy-number variation. In comparison, copy number analysis using a state-of-the-art exome CNV caller failed to estimate copy number for paralogous genes with very high mapping ambiguity and showed much lower concordance (0.565) for CNV events compared to EdgeCopy (0.908). AVAILABILITY: EdgeCopy is freely available at https://github.com/vibansal-lab/edgecopy.

Humans↗
PubMed2003

Cloning and analysis of a DNA fragment stimulating avermectin production in various Streptomyces avermitilis strains.

To isolate a gene for stimulating avermectin production, a genomic library of Streptomyces avermitilis ATCC 31267 was constructed in Streptomyces lividans TK21 as the host strain. An 8.0-kb DNA fragment that significantly stimulated actinorhodin and undecylprodigiosin production was isolated. When wild-type S. avermitilis was transformed with the cloned fragment, avermectin production increased approximately 3.5-fold. The introduction of this fragment into high-producer (ATCC 31780) and semi-industrial (L-9) strains also resulted in an increase of avermectin production by more than 2.0- and 1.4-fold, respectively. Subclones were studied to locate the minimal region involved in stimulation of pigmented-antibiotic and avermectin production. An analysis of the nucleotide sequence of the entire DNA fragment identified eight complete and one incomplete open reading frame. All but one of the deduced proteins exhibited strong homology (68 to 84% identity) to the hypothetical proteins of Streptomyces coelicolor A3(2). The orfX gene product showed no significant similarity to any other protein in the databases, and an analysis of its sequence suggested that it was a putative membrane protein. Although the nature of the stimulatory effect is still unclear, the disruption of orfX revealed that this gene was intrinsically involved in the stimulation of avermectin production in S. avermitilis.

Bacterial Proteins↗
PubMed1987

Total synthesis and further characterization of the gamma-carboxyglutamate-containing "sleeper" peptide from Conus geographus venom.

The total synthesis of the Gla-containing "sleeper" peptide (Gly-Glu-Gla-Gla-Leu-Gln-Gla-Asn-Gln-Gla-Leu-Ile-Arg-Gla-Lys-Ser-Asn-NH2 ) from Conus geographus is described. A new strategy for the synthesis of acid-sensitive peptide amides was developed, which allowed complete deprotection and cleavage of the L-gamma-carboxyglutamate-containing peptide from the 2,4-dimethoxybenzhydrylamine resin. Synthetic sleeper peptide, after preparative high-performance liquid chromatography (HPLC) purification, was shown to be identical with the native peptide by all criteria (coelution experiments of HPLC, sequence analysis, and biological activity). In addition, a developmental switch in the behavioral symptoms induced by the peptide after intracerebral administration in mice was documented. At low doses of the peptide (4-30 pmol/g), a sleeplike state was induced in mice under 2 weeks old; in contrast, older mice became markedly hyperactive. It is proposed that, in the presence of Ca2+, the sleeper peptide assumes an alpha-helical configuration in which all the gamma-carboxyglutamate residues are located on the same side of the alpha-helix.

1-Carboxyglutamic Acid↗
PubMed2006

Topology testing of phylogenies using least squares methods.

BACKGROUND: The least squares (LS) method for constructing confidence sets of trees is closely related to LS tree building methods, in which the goodness of fit of the distances measured on the tree (patristic distances) to the observed distances between taxa is the criterion used for selecting the best topology. The generalized LS (GLS) method for topology testing is often frustrated by the computational difficulties in calculating the covariance matrix and its inverse, which in practice requires approximations. The weighted LS (WLS) allows for a more efficient albeit approximate calculation of the test statistic by ignoring the covariances between the distances. RESULTS: The goal of this paper is to assess the applicability of the LS approach for constructing confidence sets of trees. We show that the approximations inherent to the WLS method did not affect negatively the accuracy and reliability of the test both in the analysis of biological sequences and DNA-DNA hybridization data (for which character-based testing methods cannot be used). On the other hand, we report several problems for the GLS method, at least for the available implementation. For many data sets of biological sequences, the GLS statistic could not be calculated. For some data sets for which it could, the GLS method included all the possible trees in the confidence set despite a strong phylogenetic signal in the data. Finally, contrary to WLS, for simulated sequences GLS showed undercoverage (frequent non-inclusion of the true tree in the confidence set). CONCLUSION: The WLS method provides a computationally efficient approximation to the GLS useful especially in exploratory analyses of confidence sets of trees, when assessing the phylogenetic signal in the data, and when other methods are not available.

Animals↗
PubMed2025

Expanding and improving analyses of nucleotide recoding RNA-seq experiments with the EZbakR suite.

Nucleotide recoding RNA sequencing methods (NR-seq; TimeLapse-seq, SLAM-seq, TUC-seq, etc.) are powerful approaches for assaying transcript population dynamics. In addition, these methods have been extended to probe a host of regulated steps in the RNA life cycle. Current bioinformatic tools significantly constrain analyses of NR-seq data. To address this limitation, we developed EZbakR (https://github.com/isaacvock/EZbakR), an R package to facilitate a more comprehensive set of NR-seq analyses, and fastq2EZbakR (https://github.com/isaacvock/fastq2EZbakR), a Snakemake pipeline for flexible preprocessing of NR-seq datasets, collectively referred to as the EZbakR suite. Together, these tools generalize many aspects of the NR-seq analysis workflow. The fastq2EZbakR pipeline can assign reads to a diverse set of genomic features (e.g., genes, exons, splice junctions), and EZbakR can perform analyses on any combination of these features. EZbakR extends standard NR-seq mutational modeling to support multi-label analyses (e.g., s4U and s6G dual labeling), and implements an improved hierarchical model to better account for transcript-to-transcript variance in metabolic label incorporation. EZbakR also generalizes dynamical systems modeling of NR-seq data to support analyses of premature mRNA processing and flow between subcellular compartments. Finally, EZbakR implements flexible and well-powered comparative analyses of all estimated parameters via design matrix-specified generalized linear modeling. The EZbakR suite will thus allow researchers to make full, effective use of NR-seq data.

Software↗
PubMed2002

Identification and characterization of a putative transcriptional regulator controlling the expression of fouling inhibitors in Pseudoalteromonas tunicata.

The dark green pigmented marine bacterium Pseudoalteromonas tunicata colonizes living surfaces and produces a range of extracellular compounds that inhibit common fouling organisms, including marine invertebrate larvae, algae, bacteria, and fungi. We have observed a positive correlation between the antifouling activity of P. tunicata strain D2 and the expression of pigmentation. To address the hypothesis that pigmentation and antifouling may be jointly regulated in this organism and to begin to identify potential regulatory elements, we used transposon mutagenesis to generate a strain of P. tunicata deficient in antifouling activity. The data presented here describe the phenotypic and molecular characterization of a nonpigmented transposon mutant strain of P. tunicata (D2W2). Analyses of the antifouling capabilities of D2W2 demonstrate that this strain is deficient in the ability to inhibit each of the target fouling organisms. Genetic analysis of D2W2 identified a gene, designated wmpR (white mutant phenotype), with high sequence similarity to transcriptional regulators ToxR from Vibrio cholerae and CadC from Escherichia coli. Two-dimensional polyacrylamide gel electrophoresis analysis revealed that WmpR is essential for the expression of a significant subset of stationary-phase-induced proteins likely to be important for the synthesis of fouling inhibitors. The identification of a gene involved in the regulation of expression of antifouling phenotypes will contribute to the understanding of the interactions between bacteria and other surface-colonizing organisms in the marine environment.

Animals↗
PubMed2003

Toxicity of Bacillus sphaericus LP1-G against susceptible and resistant Culex quinquefasciatus and the cloning of the mosquitocidal toxin gene.

Bacillus sphaericus LP1-G, belonging to flagellar serotype H3, has been found to have moderate toxicity against two resistant Culex quinquefasciatus colonies (RLCq1 and RLCq2) and the susceptible contrast (SLCq). With an aim of screening mosquitocidal acting factor, a partial genome library was prepared from a partial HindIII digest of the total DNA from Bacillus sphaericus LP1-G. Two thousand twenty Escherichia coli clones were screened for toxicity against susceptible SLCq, and a toxic clone, designated E-UL68, was chosen for further study. The recombinant E-UL68 performed toxicity against both susceptible and two resistant colonies, having the same level of toxicity as that of wide-type strain LP1-G. Sequence analysis revealed that the inserted fragment was composed of 3876 nucleotides and contained a complete gene, whose sequence was identical to that of the mtx gene from B. sphaericus SSII-1. Because the binary toxin produced during sporulation of strain LP1-G has no activity against the target mosquitoes, this indicates that the Mtx toxin or other active factors might perhaps be responsible for the toxicity of LP1-G against different colonies of mosquito larvae.

Animals↗