SEARCH · PubMed Health

Results for “sequence biases”

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 217 records · Page 12Linked to original sources

PubMed2002

Global features of the Pseudomonas putida KT2440 genome sequence.

The compositional bias of the G+C, di- and tetranucleotide contents in the 6 181 862 bp Pseudomonas putida KT2440 genome was analysed in sliding windows of 4000 bp in steps of 1000 bp. The genome has a low GC skew (mean 0.066) between the leading and lagging strand. The values of GC contents (mean 61.6%) and of dinucleotide relative abundance exhibit skewed Gaussian distributions. The variance of tetranucleotide frequencies, which increases linearly with increasing GC content, shows two overlapping Gaussian distributions of genome sections with low (minor fraction) or high variance (major fraction). Eighty per cent of the chromosome shares similar GC contents and oligonucleotide bias, but 105 islands of 4000 bp or more show atypical GC contents and/or oligonucleotide signature. Almost all islands provide added value to the metabolic proficiency of P. putida as a saprophytic omnivore. Major features are the uptake and degradation of organic chemicals, ion transport and the synthesis and secretion of secondary metabolites. Other islands endow P. putida with determinants of resistance and defenceor with constituents and appendages of the cell wall. A total of 29 islands carry the signature of mobile elements such as phage, transposons, insertion sequence (IS) elements and group II introns, indicating recent acquisition by horizontal gene transfer. The largest gene carries the most unusual sequence that encodes a multirepeat threonine-rich surface adhesion protein. Among the housekeeping genes, only genes of the translational apparatus were located in segments with an atypical signature, suggesting that the synthesis of ribosomal proteins is uncoupled from the rapidly changing translational demands of the cell by the separate utilization of tRNA pools.

Base Composition↗

PubMed2006

An exploration of 3'-end processing signals and their tissue distribution in Oryza sativa.

The 3' untranslated regions deeply affect many properties of eukaryotic mRNA. In plants, the polyadenine control signals contained in these regions seem to be more variable than of mammals. Three cDNA libraries derived from the leaf, endosperm and stem tissues of rice were sequenced from the 3'-end. Of the 9911 transcripts analyzed, 5723 unique transcripts were identified from the leaf sequences, 2934 from the endosperm and 1254 from the stem. The information entropy and two statistical methods were used to compile a list of rice poly(A) control signals. Based on their distribution, these signals can be roughly grouped into far-upstream element (FUE), near-upstream element (NUE), T-rich region (TRE) and downstream element (DE). The distribution of rice conserved regions is similar to the previous model from Arabidopsis and yeast, with a few differences in word constructions. Interestingly, we also found the word distributions were diverse in the cleavage site of downstream sequences of different rice tissues. The signal bias in downstream sequences may lead mRNA to be differently cleaved in different rice tissues.

3' Untranslated Regions↗

PubMed2001

A comprehensive collection of mouse zinc finger motifs compiled by molecular indexing.

The C(2)H(2) zinc finger motif found in many transcription factors is thought to be for nucleic acid binding and/or dimerization. Nearly 1% of eukaryote genes are estimated to encode this motif. We investigated the gene family encoding this motif in the Mus musculus mRNA by molecular indexing, a technique used to select a subpopulation of cDNA by ligation of adapters to cDNA fragments digested by a class IIS restriction enzyme(s). In place of oligo-dT primers in the original method, a polymerase chain reaction primer designed based on the conserved sequence of the C(2)H(2) zinc finger protein stranded cDNA was made from various mouse tissue mRNAs, digested with FokI and BsmAI, and joined with adapters. Amplification of the cDNA with an adapter primer and zinc finger-specific primer yielded products enriched in zinc finger protein genes. Fragments were separated by subsequent denaturing polyacrylamide gel electrophoresis, and characterized by DNA sequencing. Consequently, 259 C(2)H(2) zinc finger motif sequences were obtained, among which 166 were novel. Combined with the reported sequences, these mouse motif sequences were compared with those of other species such as Saccharomyces cerevisiae and Caenorhabditis elegans. Some of the amino acids in the motif sequence showed strong bias among species. Most of the novel sequences were supposed to be DNA-binding according to the surface potential of predicted tertiary structures.

Amino Acid Motifs↗

PubMed2002

Construction and expression of a reshaped VH domain against human CD28 molecules.

Compared with the amino acid sequence of a mouse anti-human CD28 VH domain antibody, the two most homologous sequences of human antibodies were pulled out from Genbank. One of them was used as the main template for the framework regions of the reshaped VH domain. While the original mouse antibody CDRs were inserted into the human acceptor FRs, some residues in human acceptor FRs, which were different from those of the original mouse FRs in corresponding positions, were then determined or, alternatively, mutagenized to their conservative properties in kappa classification. Based on the amino acid sequences of the designed VH domain, the nucleotide sequence was deduced by using E. coli bias codons. The sequence was split into ten 30 to 60 nucleotide fragments for synthesizing, then annealed and amplified by overlap PCR. Taq DNA polymerase was used in a buffer with high Mg2+ concentration to induce more random mutations, both in FRs and CDRs. A phage display library was constructed by cloning these PCR products. After three rounds of panning, several reshaped VH with high antigen binding activity were obtained. One of them had the same CDR amino acid sequences as that of the original mouse VH domain. Further study showed that it retained a high antigen binding affinity after being expressed in E. coli BL21 (DE3).

Antibody Affinity↗

PubMed1992

Substitutional bias confounds inference of cyanelle origins from sequence data.

Available molecular and biochemical data offer conflicting evidence for the origin of the cyanelle of Cyanophora paradoxa. We show that the similarity of cyanelle and green chloroplast sequences is probably a result of these two lineages independently developing the same pattern of directional nucleotide change (substitutional bias). This finding suggests caution should be exercised in the interpretation of nucleotide sequence analyses that appear to favor the view of a common endosymbiont for the cyanelle and chlorophyll-b-containing chloroplasts. The data and approaches needed to resolve the issue of cyanelle origins are discussed. Our findings also have general implications for phylogenetic inference under conditions where the base compositions (compositional bias) of the sequences analyzed differ.

Base Composition↗

PubMed2005

Score statistics of global sequence alignment from the energy distribution of a modified directed polymer and directed percolation problem.

Sequence alignment is one of the most important bioinformatics tools for modern molecular biology. The statistical characterization of gapped alignment scores has been a long-standing problem in sequence alignment research. Using a variant of the directed path in random media model, we investigate the score statistics of global sequence alignment taking into account, in particular, the compositional bias of the sequences compared. Such statistics are used to distinguish accidental similarity due to compositional similarity from biologically significant similarity. To accommodate the compositional bias, we introduce an extra parameter p indicating the probability for positive matching scores to occur. When p is small, a high scoring alignment obviously cannot come from compositional similarity. When p is large, the highest scoring point within a global alignment tends to be close to the end of both sequences, in which case we say the system percolates. By applying finite-size scaling theory on percolating probability functions of various sizes (sequence lengths), the critical p at infinite size is obtained. For alignment of length t, the fact that the score fluctuation grows as chi(t)1/3 is confirmed upon investigating the scaling form of the alignment score. Using the Kolmogorov-Smirnov statistics test, we show that the random variable , if properly scaled, follows the Tracy-Widom distributions: Gaussian orthogonal ensemble for p slightly larger than pc and Gaussian unitary ensemble for larger p. Although these results deepen our understanding of the distribution of alignment scores, the use of these results in practical applications remains somewhat heuristic and needs to be further developed. Nevertheless, the possibility of characterizing score statistics for modest system size (sequence lengths), via proper reparametrization of alignment scores, is illustrated.

Algorithms↗

PubMed2004

Statistical modeling of sequencing errors in SAGE libraries.

MOTIVATION: Sequencing errors may bias the gene expression measurements made by Serial Analysis of Gene Expression (SAGE). They may introduce non-existent tags at low abundance and decrease the real abundance of other tags. These effects are increased in the longer tags generated in LongSAGE libraries. Current sequencing technology generates quite accurate estimates of sequencing error rates. Here we make use of the sequence neighborhood of SAGE tags and error estimates from the base-calling software to correct for such errors. RESULTS: We introduce a statistical model for the propagation of sequencing errors in SAGE and suggest an Expectation-Maximization (EM) algorithm to correct for them given observed sequences in a library and base-calling error estimates. We tested our method using simulated and experimental SAGE libraries. When comparing SAGE libraries, we found that sequencing errors can introduce considerable bias. High abundance tags may be falsely called as significantly differentially expressed, especially when comparing libraries with different levels of sequencing errors and/or of different size. Truly, differentially expressed tags have decreased significance as 'true'-tag counts are generally underestimated. This may alter if tags near the threshold of differential expression are called significant. Moreover, the number of different transcripts present in a library is overestimated as false tags are introduced at low abundance. Our correction method adjusts the tag counts to be closer to the true counts and is able to partly correct for biases introduced by sequencing errors. AVAILABILITY: An implementation using R is distributed as an R package. An online version is available at http://tagcalling.mbgproject.org

Algorithms↗

PubMed2005

On the correlation between composition and site-specific evolutionary rate: implications for phylogenetic inference.

Model-based phylogenetic reconstruction methods traditionally assume homogeneity of nucleotide frequencies among sequence sites and lineages. Yet, heterogeneity in base composition is a characteristic shared by most biological sequences. Compositional variation in time, reflected in the compositional biases among contemporary sequences, has already been extensively studied, and its detrimental effects on phylogenetic estimates are known. However, fewer studies have focused on the effects of spatial compositional heterogeneity within genes. We show here that different sites in an alignment do not always share a unique compositional pattern, and we provide examples where nucleotide frequency trends are correlated with the site-specific rate of evolution in RNA genes. Spatial compositional heterogeneity is shown to affect the estimation of evolutionary parameters. With standard phylogenetic methods, estimates of equilibrium frequencies are found to be biased towards the composition observed at fast-evolving sites. Conversely, the ancestral composition estimates of some time-heterogeneous but spatially homogeneous methods are found to be biased towards frequencies observed at invariant and slow-evolving sites. The latter finding challenges the result of a previous study arguing against a hyperthermophilic last universal ancestor from the low apparent G + C content of its rRNA sequences. We propose a new model to account for compositional variation across sites. A Gaussian process prior is used to allow for a smooth change in composition with evolutionary rate. The model has been implemented in the phylogenetic inference software PHASE, and Bayesian methods can be used to obtain the model parameters. The results suggest that this model can accurately capture the observed trends in present-day RNA sequences.

Evolution, Molecular↗

PubMed2004

Inferring property selection pressure from positional residue conservation.

In this study, we attempt to understand and explain positional selection pressure in terms of underlying physical and chemical properties. We propose a set of constraining assumptions about how these pressures behave, then describe a procedure for analysing and explaining the distribution of residues at a particular position in a multiple sequence alignment. In contrast to previous approaches, our model takes into account both amino acid frequencies and a large number of physical-chemical properties. By analysing each property separately, it is possible to identify positions where distinct conservation patterns are present. In addition, the model can easily incorporate sequence weights that adjust for bias in the sample sequences. Finally, a test of statistical significance is provided for our conservation measure. The applicability of this method is demonstrated on two HIV-1 proteins: Nef and Env. The tools, data and results presented in this article are available at http://flan.blm.cs.cmu.edu.

Algorithms↗

PubMed1983

Homologous nucleotide sequences at the 5' termini of messenger RNAs synthesized from the yeast enolase and glyceraldehyde-3-phosphate dehydrogenase gene families. The primary structure of a third yeast glyceraldehyde-3-phosphate dehydrogenase gene.

Genomic DNA containing a third yeast glyceraldehyde-3-phosphate dehydrogenase structural gene has been isolated on a bacterial plasmid designated pgap11. The complete nucleotide sequence of this structural gene was determined. The gene contains no intervening sequences, codon usage is highly biased, and the nucleotide sequence of the coding portion of this gene is 90% homologous to the other two glyceraldehyde-3-phosphate dehydrogenase genes (Holland, J. P., and Holland, M. J. (1980) J. Biol. Chem. 255, 2596-2605). Based on the extent of nucleotide sequence divergence among the three glyceraldehyde-3-phosphate dehydrogenase genes, it is likely that they arose as a consequence of two duplication events and the gene contained on the hybrid plasmid designated pgap11 is a product of the first duplication event. All three structural genes share extensive nucleotide sequence homology in the 5'-noncoding regions adjacent to the three respective translational initiation codons. The gene contained on pgap11 is not homologous to the others downstream from the respective translational termination codon, however. The 5' termini of messenger RNAs synthesized from the three glyceraldehyde-3-phosphate dehydrogenase and two yeast enolase genes have been mapped to sites ranging from 36 to 82 nucleotides upstream from the respective translational initiation codons. In each case the 5' terminus of the mRNA maps to a region of strong nucleotide sequence homology which is shared by all five structural genes. These latter data confirm that all five structural genes are expressed during vegetative cell growth and further support the hypothesis that a portion of the 5'-noncoding flanking region of the yeast glyceraldehyde-3-phosphate dehydrogenase and enolase genes evolved from a common precursor sequence.

Base Sequence↗

PubMed1998

PCR bias toward the wild-type k-ras and p53 sequences: implications for PCR detection of mutations and cancer diagnosis.

PCR-based cancer diagnosis requires detection of rare mutations in k-ras, p53 or other genes. The assumption has been that mutant and wild-type sequences amplify with near equal efficiency, so that they are eventually present in proportions representative of the starting material. Work on factor IX suggests that this assumption is invalid for one case of near-sequence identity. To test the generality of this phenomenon and its relevance to cancer diagnosis, primers distant from point mutations in p53 and k-ras were used to amplify wild-type and mutant sequences from these genes. A substantial bias against PCR amplification of mutants was observed for two regions of the p53 gene and one region of k-ras. For k-ras and p53, bias was observed when the wild-type and mutant sequences were amplified separately or when mixed in equal proportions before PCR. Bias was present with proofreading and non-proofreading polymerases. Mutant and wild-type segments of the factor V, cystic fibrosis transmembrane conductance regulator and prothrombin genes were amplified and did not exhibit PCR bias. Therefore, the assumption of equal PCR efficiency for point mutant and wild-type sequences is invalid in several systems. Quantitative or diagnostic PCR will require validation for each locus, and enrichment strategies may be needed to optimize detection of mutants.

Base Sequence↗

PubMed1985

Sense codons are found in specific contexts.

The sequence environment of codons in structural genes has been investigated statistically, using computer methods. A set of Escherichia coli genes with abundant products was compared with a set having low gene product levels, in order to detect potential differences associated with expression. The results show striking non-randomness in the nucleotides occurring near codons. These effects are, unexpectedly, very much larger and more homogeneous among the genes with rare products. The intensity of effects in weakly expressed genes suggests that such non-random sequence environments decrease expression. In the weakly expressed set of genes, the 5' neighbor of a codon, and all positions of the 3' neighbor codon are biased. In the highly expressed genes, the first nucleotide of the next codon is a uniquely affected site. The distribution of non-randomness in weakly expressed genes suggests that sequence bias is primarily due to a constraint acting directly on the secondary or tertiary structure of the codon/anticodon. In highly expressed genes, the observed bias suggests an interaction between the codon/anticodon and a site outside the codon/anticodon. Much of the tendency to non-random near-neighbor sequences in weakly expressed genes can be ascribed to a correlation between nearby nucleotides and the wobble nucleotide of the codon, despite the fact that selection of such correlations will alter the amino acid sequence. The favored pattern, in genes expressed at low level, is R YYR or Y RRY. R indicates purine, Y indicates pyrimidine; the space is the boundary between codons. It seems likely that this preference for nearby sequences is the physical basis of the genetic context effect. Under this assumption such sequence biases will affect expression. On this basis, we predict new sites for contextual mutations which decrease expression, and suggest strategy for the design of messages having optimal translational activity.

Amino Acids↗

PubMed1994

Evidence that GHN phase bias does not constitute a framing code.

It has been suggested that a triplet repeated pattern found in coding sequences, the G-nonG-N or GHN phase bias, serves a framing function during protein synthesis. To test this idea, the framing characteristics of a highly GHN biased sequence are examined. No effects on reading frame maintenance are observed despite the use of sensitive frameshift assays. Specifically, first the GHN phase is not more accurate than the alternative overlapping phases (i.e., HNG and NGH). Second, ribosomes do not exhibit any significant tendency to slip from the alternative frames into the GHN pahse. In addition, examination of Escherichia coli programmed frameshift sites does not support roles for GHN phase bias in programmed frameshifting. Framing functions for GHN phase bias, if they occur at all, must be extremely limited.

Bacterial Proteins↗

PubMed2005

Biases in phylogenetic estimation can be caused by random sequence segments.

We consider the effects of fully or partially random sequences on the estimation of four-taxon phylogenies. Fully or partially random sequences occur when whole subsets of sequences or some sites for subsets of sequences are independent of sequence data for the other taxa. Random sequences can be a consequence of misalignment or because sites evolve at very fast rates in some portions of a tree, a situation that occurs especially in analyses involving deep divergence times. One might reasonably speculate that random sites will only add noise to the estimation of a phylogeny. We show that in the case that a random sequence is added to a three-taxa alignment, it is more likely to be a neighbor of the sequence corresponding to the longest branch in the three-taxon tree. Surprisingly, when only about half of the sites show randomness, a long-branch-repels form of small sample bias occurs, and when a minority of sites show randomness this becomes a long-branch-attraction bias again. The most serious bias, one that does not vanish with increasing sequence length, occurs when more than one sequence is partially random. If there is a large amount of overlap in the random sites for two sequences, those two sequences will be attracted to each other; otherwise, they will repel each other. Random sequences or sites can, therefore, cause complicated biases in phylogenetic inference. We suggest performing analyses with and without potentially saturated sequences and/or misaligned sites, to check that these biases are not affecting the inferred branching pattern.

Animals↗

PubMed1999

Predicting regional mutability in antibody V genes based solely on di- and trinucleotide sequence composition.

Somatic mutations are not distributed randomly throughout Ab V region genes. A sequence-specific target bias is revealed by a defined hierarchy of mutability among di- and trinucleotide sequences located within Ig intronic DNA. Here we report that the di- and trinucleotide mutability preference pattern is shared by mouse intronic JH and Jkappa clusters and by human VH genes, suggesting that a common mutation mechanism exists for all Ig V genes of both species. Using di- and trinucleotide target preferences, we performed a comprehensive analysis of human and murine germline V genes to predict regional mutabilities. Heavy chain genes of both species exhibit indistinguishable patterns in which complementarity-determining region 1 (CDR1), CDR2, and framework region 3 (FR3) are predicted to be more mutable than FR1 and FR2. This prediction is borne out by empirical mutation data from nonproductively rearranged human VH genes. Analysis of light chain genes in both species also revealed a common, but unexpected, pattern in which FR2 is predicted to be highly mutable. While our analyses of nonfunctional Ig genes accurately predicts regional mutation preferences in VH genes, observed relative mutability differences between regions are more extreme than expected. This cannot be readily accounted for by nascent mRNA secondary structure or by a supplemental gene conversion mechanism that might favor nucleotide replacements in CDR. Collectively, our data support the concept of a common mutation mechanism for heavy and light chain genes of mice and humans with regional bias that is qualitatively, but not quantitatively, accounted for by short nucleotide sequence composition.

Animals↗

PubMed2001

Computational identification of noncoding RNAs in E. coli by comparative genomics.

Some genes produce noncoding transcripts that function directly as structural, regulatory, or even catalytic RNAs [1, 2]. Unlike protein-coding genes, which can be detected as open reading frames with distinctive statistical biases, noncoding RNA (ncRNA) gene sequences have no obvious inherent statistical biases [3]. Thus, genome sequence analyses reveal novel protein-coding genes, but any novel ncRNA genes remain invisible. Here, we describe a computational comparative genomic screen for ncRNA genes. The key idea is to distinguish conserved RNA secondary structures from a background of other conserved sequences using probabilistic models of expected mutational patterns in pairwise sequence alignments. We report the first whole-genome screen for ncRNA genes done with this method, in which we applied it to the "intergenic" spacers of Escherichia coli using comparative sequence data from four related bacteria. Starting from >23,000 conserved interspecies pairwise alignments, the screen predicted 275 candidate structural RNA loci. A sample of 49 candidate loci was assayed experimentally. At least 11 loci expressed small, apparently noncoding RNA transcripts of unknown function. Our computational approach may be used to discover structural ncRNA genes in any genome for which appropriate comparative genome sequence data are available.

Animals↗

PubMed2002

eCodonOpt: a systematic computational framework for optimizing codon usage in directed evolution experiments.

We present a systematic computational framework, eCodonOpt, for designing parental DNA sequences for directed evolution experiments through codon usage optimization. Given a set of homologous parental proteins to be recombined at the DNA level, the optimal DNA sequences encoding these proteins are sought for a given diversity objective. We find that the free energy of annealing between the recombining DNA sequences is a much better descriptor of the extent of crossover formation than sequence identity. Three different diversity targets are investigated for the DNA shuffling protocol to showcase the utility of the eCodonOpt framework: (i) maximizing the average number of crossovers per recombined sequence; (ii) minimizing bias in family DNA shuffling so that each of the parental sequence pair contributes a similar number of crossovers to the library; and (iii) maximizing the relative frequency of crossovers in specific structural regions. Each one of these design challenges is formulated as a constrained optimization problem that utilizes 0-1 binary variables as on/off switches to model the selection of different codon choices for each residue position. Computational results suggest that many-fold improvements in the crossover frequency, location and specificity are possible, providing valuable insights for the engineering of directed evolution protocols.

Aldose-Ketose Isomerases↗

PubMed2001

Evolution of transcript structure and base composition of rDNA expansion segment D3 in ticks.

Four to thirty-two copies of the rDNA 28S gene expansion segment D3 and flanking H14 stem were sequenced in six species of ticks (Ixodes: Ixodidae: Acari). Sequence match among species varied from 66% to 97%. Sequence length averaged 130 bases in I. persulcatus across eight Eurasian sites and averaged 186 bases in five other species across 19 Eurasian and North American sites. The difference in length represents one or more deletions totalling about 60 bases that correspond to stems S3 or S4 of the folded transcript. The typical transcript conformation was observed as one possible low energy structure in the five species of longer D3. The structure entails a basal loop with four stem/loop structures, S1-S4 (moving 5' to 3') atop stem H14. A secondary structure lacking S4 but possessing all other putative standard features of the D3 transcript is possible with the shorter I. persulcatus sequences. Interspecific sequence differences occur at higher frequency in loops and bulges vs. complementary pairing regions of stems. Insertion/deletion events (indels) and base substitutions accounted equally for sequence differences. Indels are flanked by similar sequences, suggesting that they occur by slippage during replication. The D3 of Ixodes species is composed of a degenerate set of subrepeats. Thus, unequal exchange among subrepeats may have caused the reduction in length of the I. persulcatus D3. Compensatory base substitution and compensatory insertion/deletion events are indicated by the failure of mutations to affect secondary structure. Transversions accounted for 64% of sequence differences and were biased toward the gain of G and U and the loss of A and C. This bias could re-establish intramolecular base pairing when disrupted by insertions or deletions that shift one side of a stem relative to the other. The distribution of sequence differences, biased substitution, and conservation of transcript conformation in D3 suggest selective constraint.

Animals↗