PubMed Health⌕ Search

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 361 records · Page 20Linked to original sources

Nucleotide bias causes a genomewide bias in the amino acid composition of proteins.

We analyzed the nucleotide contents of several completely sequenced genomes, and we show that nucleotide bias can have a dramatic effect on the amino acid composition of the encoded proteins. By surveying the genes in 21 completely sequenced eubacterial and archaeal genomes, along with the entire Saccharomyces cerevisiae genome and two Plasmodium falciparum chromosomes, we show that biased DNA encodes biased proteins on a genomewide scale. The predicted bias affects virtually all genes within the genome, and it could be clearly seen even when we limited the analysis to sets of homologous gene sequences. Parallel patterns of compositional bias were found within the archaea and the eubacteria. We also found a positive correlation between the degree of amino acid bias and the magnitude of protein sequence divergence. We conclude that mutational bias can have a major effect on the molecular evolution of proteins. These results could have important implications for the interpretation of protein-based molecular phylogenies and for the inference of functional protein adaptation from comparative sequence data.

Amino Acids↗

Cloning and sequencing of the dnaK region of Streptomyces coelicolor A3(2).

The dnaK homologue of Streptomyces coelicolor A3(2) strain M145 has been cloned and sequenced. Nucleotide sequence analysis of a 2.5-kb region revealed an open reading frame (ORF) encoding a predicted DnaK protein of 618 amino acids (M(r) = 66,274). The dnaK coding sequence displays extreme codon bias and shows a strong preference for CGY and GGY, for Arg and Gly codons, respectively. The predicted DnaK sequence has a high Lys:Arg ratio which is not typical of streptomycete proteins. The region immediately downstream from dnaK contains an ORF for a GrpE-like protein; the predicted start codon of grpE overlaps the last two codons of dnaK, indicating that the two genes are translationally coupled. This organisation differs from that reported for other prokaryotes.

Amino Acid Sequence↗

On the statistics of ultrasonic spectral parameters.

Several factors affect the accuracy and precision of ultrasonic spectrum analysis, which is used for characterization of normal and diseased tissue in a variety of organs. For example, averaging procedures and the sequence of operations affect the accuracy and precision of spectrum analysis. Averaging procedures and logarithmic conversion (i.e., conversion to dB) introduce a constant bias that affects spectral amplitudes and the values of intercept and midband fit; the bias depends on the sequencing of the log conversion and averaging as well as the number of independent spectra or spectral parameters that are averaged. We derive expressions that permit correction of such biases. Furthermore, we show that standard deviations for slope and midband-fit estimation can be minimized by averaging spectra before dB conversion and before computing spectral parameters by linear regression. Experimental results using phantoms agree remarkably with theoretical predictions for the data window functions studied in this article, Hamming and rectangular.

Humans↗

Mutational bias affects protein evolution in flowering plants.

Amino acid sequences from several thousand homologous gene pairs were compared for two plant genomes, Oryza sativa and Arabidopsis thaliana. The Arabidopsis genes all have similar G+C (guanine plus cytosine) contents, whereas their homologs in rice span a wide range of G+C levels. The results show that those rice genes that display increased divergence in their nucleotide composition (specifically, increased G+C content) showed a corresponding, predictable change in the amino acid compositions of the encoded proteins relative to their Arabidopsis homologs. This trend was not seen in a "control" set of rice genes that had nucleotide contents closer to their Arabidopsis homologs. In addition to showing an overall difference in the amino acid composition of the homologous proteins, we were also able to investigate the biased patterns of amino acid substitution since the divergence of these two species. We found that the amino acid exchange matrix was highly asymmetric when comparing the High G+C rice genes with their Arabidopsis homologs. Finally, we investigated the possible causes of this biased pattern of sequence evolution. Our results indicate that the biased pattern of protein evolution is the consequence, rather than the cause, of the corresponding changes in nucleotide content. In fact, there is an even more marked asymmetry in the patterns of substitution at synonymous nucleotide sites. Surprisingly, there is a very strong negative correlation between the level of nucleotide bias and the length of the coding sequences within the rice genome. This difference in gene length may provide important clues about the underlying mechanisms.

Amino Acids↗

Sequence evolution and phylogenetic signal in control-region and cytochrome b sequences of rainbow fishes (Melanotaeniidae).

The nucleotide sequences of segments of the cytochrome b gene (351 bp), the tRNA(Pro) gene (49 bp), and the control region (approximately 313 bp) of mitochondrial DNA were obtained from 26 fish representing different populations and species of Melanotaenia and one species of Glossolepis, freshwater rainbow fishes confined to Australia and New Guinea. The purpose was to investigate relative rates and patterns of sequence evolution. Overall levels of divergence were similar for the cytochrome b and tRNA control-region sequences, both ranging from < 1% within subspecies to 15%-19% between genera. However, the patterns of sequence evolution differed. For the cytochrome b gene, transitions consistently exceeded transversions, the bias ranging from 4.2:1 to 2:1, depending on the level of sequence divergence. However, in the control-region sequence, a bias toward transitions (2:1) was observed only in comparisons between very similar sequences, and transversions outnumbered transitions in comparisons of divergent sequences. Graphic comparisons suggested that the control region was saturated for transitions at relatively low levels of sequence divergence but accumulated transversions at a greater rate than did the cytochrome b sequence. These distinct patterns of base substitution are associated with differences in A+T content, which is 70% for the tRNA control-region segment versus 50% for cytochrome b. A test for skewness in the distribution of lengths of random trees indicated that both segments contained phylogenetic signal. Parsimony analyses of the data from the two regions, with or without weighting schemes appropriate to the respective patterns of sequence evolution, identified the same five groupings of sequences, but the relationships among the groups differed. However, in most cases the branches uniting different combinations of groups were poorly supported, and the differences among topologies were insignificant. Considering the observed patterns of base substitution and the results of the phylogenetic analyses, we deduce that both the control region and cytochrome b are appropriate for population genetic studies but that the control region is less effective than cytochrome b for resolving relationships among divergent lineages of rainbow fishes.

Animals↗

Human retroelements may introduce intragenic polyadenylation signals.

In the human genome, the insertion of LINE-1 and Alu elements can affect genes by sequence disruption, and by the introduction of elements that modulate the gene's expression. One of the modulating sequences retroelements may contribute is the canonical polyadenylation signal (pA), AATAAA. L1 elements include these within their own sequence and AATAAA sequences are commonly created in the A-rich tails of both SINEs and LINEs. Computational analysis of 34 genes randomly retrieved from the human genome draft sequence reveals an orientation bias, reflected as a lower number of L1s and Alus containing the pA in the same orientation as the gene. Experimental studies of Alu-based pA sequences when placed in pol II or pol III transcripts suggest that the signal is very weak, or often not used at all. Because the pA signal is highly affected by the surrounding sequence, it is likely that the Alu constructs evaluated did not provide the required recognition signals to the polyadenylation machinery. Although the effect of pA signals contributed by Alus is individually weak, the observed reduction of "sense" oriented pA-containing L1 and Alu elements within genes reflects that even a modest influence causes a change in evolutionary pressure, sufficient to create the biased distribution.

Base Sequence↗

Annotation pattern of ESTs from Spodoptera frugiperda Sf9 cells and analysis of the ribosomal protein genes reveal insect-specific features and unexpectedly low codon usage bias.

MOTIVATION: A whole set of Expressed Sequence Tags (ESTs) from the Sf9 cell line of Spodoptera frugiperda is presented here for the first time. By this way we want to identify both conserved and specific genes of this pest species. We also expect from this analysis to find a class of protein sequences providing a tool to explore genomic features and phylogeny of Lepidoptera. RESULTS: The ESTs display both housekeeping as well as developmentally regulated genes, and a high percentage of sequences with unknown function. Among the identified ORFs, almost all ribosomal proteins (RPs) were found with high EST redundancy and hence sequence accuracy. The codon usage found among RP genes is in average surprisingly much less biased in Lepidoptera than in other organisms. Other Spodoptera genes also displayed a low bias, suggesting a general genome expression feature in this Lepidoptera. We also found that the L35A and L36 RP sequences, respectively, display 40 and 10 amino-acid insertions, both being present only in insects. Sequence analysis suggests that they are probably not subjected to a strong selective pressure and may be good phylogenetic markers for Lepidoptera. Most interestingly, the Lepidoptera sequences of 9 RP genes displayed a specific signature different from the canonical one. We conclude that the RP family allows valuable comparative genomics and phylogeny of Lepidoptera. AVAILABILITY: All EST sequence data are available from the private 'Spodo-Base' upon request.

Abstracting and Indexing↗

The primary structure of the Saccharomyces cerevisiae gene for 3-phosphoglycerate kinase.

The DNA sequence of the gene for the yeast glycolytic enzyme, 3-phosphoglycerate kinase (PGK), has been obtained by sequencing part of a 3.1 kbp HindIII fragment obtained from the yeast genome. The structural gene sequence corresponds to a reading frame of 1251 bp coding for 416 amino acids with no intervening DNA sequences. The amino acid sequence is approximately 65 percent homologous with human and horse PGK protein sequences and is in general agreement with the published protein sequence for yeast PGK. As for other highly expressed structural genes in yeast, the coding sequence is highly codon biased with 95 percent of the amino acids coded for by a select 25 codons (out of 61 possible). Besides structural DNA sequence, 291 bp of 5'-flanking sequence and 286 bp of 3'-flanking sequence were determined. Transcription starts 36 nucleotides upstream from the translational start and stops 86-93 nucleotides downstream from the translational stop. These results suggest a non-polyadenylated mRNA length of 1373 to 1380 nucleotides, which is consistent with the observed length of 1500 nucleotides for polyadenylated PGK mRNA. A sequence TATATATAAA is found at 145 nucleotides upstream from the translational start. This sequence resembles the TATAAA box that is possibly associated with RNA polymerase II binding.

Amino Acid Sequence↗

Increased sensitivity of tryptic peptide detection by MALDI-TOF mass spectrometry is achieved by conversion of lysine to homoarginine.

Mass spectrometric techniques for identification of proteins by "mass fingerprinting" (matching the masses of tryptic peptides from a protein digest to the theoretical peptides in a database) such as matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) are rapidly growing in popularity as the demand for high throughput analysis of the proteome increases. This is due, in part, to the ability to automate the technique and the rapid rate with which mass spectra may be acquired. An important factor in the accuracy of the technique is the number of tryptic peptides that are identified in the various searching algorithms that exist. The greater sequence coverage of the parent protein that is obtained, the higher the level of confidence in the identification that is determined. One impediment to high levels of sequence coverage is the bias of MALDI-TOF mass spectrometry to arginine-containing peptides. Increasing the sensitivity to lysine-containing peptides should increase the sequence coverage obtained. In order to achieve this result we have developed conditions to modify the epsilon-amine group of lysine in tryptic peptides with O-methylisourea. The conditions utilized result in the conversion of lysine to homoarginine with no modification of the amine terminus of the peptides. The sensitivity of MALDI-TOF mass spectrometry detection of peptides was increased dramatically following modification. The modification chemistry may be applied to tryptic peptide mixtures prior to desalting and spotting onto MALDI-TOF plates. This technique will be particularly useful for identifying proteins with a high lysine/arginine ratio.

Amino Acid Sequence↗

Sequence requirements for stabilization of a peptide reverse turn in water solution--proline is not essential for stability.

To probe the sequence requirements for stabilization of a reverse turn conformation in a short peptide in water solution, the behavior of two series of peptides was investigated by nuclear magnetic resonance (NMR) spectroscopy. The peptides have the general sequences XPGDV and AXGDV, where X is a representative subset of all 20 naturally occurring amino acids. The residues chosen at positions 3 and 4, Gly and Asp, respectively, were shown to give the greatest population of reverse turns in a previous study [Dyson, H. J., Rance, M., Houghten, R. A., Lerner, R. A. & Wright, P. E. (1988) J. Mol. Biol. 201, 161-200]. Within this framework, the identity of the first residue of the turn (X in XPGDV) does not greatly influence the turn population, although a small but significant increase is observed for residues such as Ala which have a preference for backbone conformations in the alpha region of (phi,psi) space. The series AXGDV was initially studied for completeness only, since it was expected that the turn would not be stabilized in such a small linear peptide in the absence of proline. In contrast, it appears that a significant population of type II turn conformations is to be found in peptides in the series AXGDV, although proline remains one of the most favorable residues at position 2. These results indicate that while residues at all positions within the turn can influence the turn population, the presence of Gly-Asp as the third and fourth members of the sequence gives a strong bias towards type II turn formation regardless of the residues at positions 1 and 2. Our results give a final prediction that the sequence with the highest intrinsic propensity for turn formation is APGD.

Amino Acid Sequence↗

An improved transposon for the halophilic archaeon Haloarcula hispanica.

An improved transposon (ThD73) for Haloarcula hispanica is described. Based on the halobacterial insertion sequence ISH28, it showed little target sequence specificity but was biased toward a lower G+C content. Twenty randomly selected ThD73 mutants were analyzed, and the DNA flanking their insertions revealed several recognizable sequences, including two (unrelated) ISH elements.

Base Sequence↗

cDNA cloning of rice lipoxygenase L-2 and characterization using an active enzyme expressed from the cDNA in Escherichia coli.

A full-length cDNA of rice lipoxygenase L-2 was cloned from 3-day-old seedlings. The identity of the clone was determined by amino acid sequencing of selected peptides of the purified enzyme and immunological characterization of an active enzyme that was produced from the cDNA in Escherichia coli by cultivation at 15 degrees C. The nucleotide sequence showed a strong bias toward G and C in the selection of nucleotides, especially at the third position of the codons (93% G/C). The complete amino acid sequence of the enzyme was deduced from the nucleotide sequence. The molecular mass of the enzyme was calculated to be 96,657 Da based on 865 amino acids. The amino acid sequence shares similarity with those of dicot lipoxygenases throughout the enzyme at a level of 50%. A hydropathy profile calculated from the amino acid sequence resembled those of dicot lipoxygenases, suggesting conservation of the secondary structure of these enzymes. The active enzyme, expressed in Escherichia coli, was characterized for pH dependence of the enzyme activity, intramolecular specificity, heat stability and Km. The enzyme had the same properties as the L-2 enzyme that was isolated from seedlings, but differed from the lipoxygenase L-3 isolated from mature plants.

Amino Acid Sequence↗

Precision ID mtDNA Whole Genome Panel and sequencing of telogen hairs - perspectives for validation and implementation in casework.

Shed hair is a commonly encountered type of forensic evidence. Shed telogen hairs generally contain insufficient or highly degraded nuclear DNA for STR profiling; however, mtDNA analysis of telogen hair and hair shafts remains possible. We validated whole mitochondrial genome (mtGenome) sequencing using the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific) and subsequently implemented the panel for the analysis of telogen hair, buccal, and casework samples. We analysed 90 diluted DNA samples containing 3-3,600 mtDNA copies, shed telogen hairs and their corresponding mtDNA from buccal swabs from 91 individuals, and 11 archived DNA extracts from hair samples in criminal cases. Complete mtGenome sequences were consistently recovered in 99% of samples across DNA dilution series at DNA input levels as low as 47 mtDNA copies, demonstrating the assay's robustness under low-template conditions. We obtained complete and reproducible mtGenome sequences with &#x2265;&#x2009;327 mtDNA copies/&#xb5;L from telogen hair samples. After applying ISFG recommendations and excluding low-confidence discrepancies associated with high-strand bias, heteroplasmic variants and sequencing artifacts, mtGenome sequence concordance increased from 93.4% to 100%. None of the 16 negative controls produced complete mtDNA sequences. Six negative controls showed low-level mtDNA signal (2-8 variants), consisting predominantly of common polymorphisms. These samples did not yield complete mtGenome sequences and showed no correspondence to any of the analysed samples. Finally, archived telogen hair samples from criminal cases presented complete mtGenome sequences with an average read depth of 1,037x.Our findings highlight the reliability of mtDNA analysis of telogen hairs using the Precision ID mtDNA Whole Genome Panel for implementation in forensic casework.

Forensic casework↗

Hamming chromatography.

Selection of molecules with desired properties from random pools of biopolymers has become a powerful tool in biotechnology. On designing an evolution experiment, a certain knowledge of the concomitant fitness landscape is clearly helpful to set up the optimal experimental conditions. The correlation function is a useful means of characterizing a given landscape, since it can be efficiently measured if one has a method of separating a pool of random sequences according to their Hamming distance from a moderately small number of test sequences. In this paper we describe a special type of hybridization chromatography, where a mixture of oligomers (partially) complementary to a given test sequence is hybridized to the test sequence, covalently bound to a matrix. DNA oligomers are eluted in an 'effective temperature gradient' using conditions that minimize the differences of effects of GC versus AT pairs on the melting temperatures. This method should be a means to quickly separate error classes and thus be the crucial step in characterizing fitness landscapes of biopolymers through an experimental approach. It would also be a useful tool to design sequence pools with a bias towards desired mutant spectra.

Bacteriophage T7↗

Dual enigma of somatic hypermutation of immunoglobulin variable genes: targeting and mechanism.

The immunoglobulin loci are uniquely unstable regions of the genome which undergo as much mutation and selection in a matter of days as a species can undergo in generations of evolution. We have studied the mutational pattern and targeting of this unusual hypermutation process over the past 16 years. The pattern of somatic mutations in rearranged variable (V) genes differs from the pattern of meiotic mutations, indicating that a different mechanism generates hypermutation than generates spontaneous mutation. Hypermutations begin on the 5' end of rearranged V genes downstream of the transcription initiation site and continue through the V exon and into the 3'-flanking region before tapering off. Mutations are located randomly throughout the DNA sequence and exhibit strand bias. The targeting of mutations to the region in and around the rearranged V gene appears to require interactions between the promoter and downstream intronic DNA sequences. The same mechanism that initiates hypermutation around V genes may also produce double-strand breaks that catalyze homologous recombination between rearranged V genes on two chromosomal alleles. With this data we have built a model of hypermutation which predicts that V-region DNA is destabilized at the nuclear matrix during transcription and undergoes strand breaks.

Animals↗

A dissection of volatility in yeast.

It has been suggested that volatility, the proportion of mutations which change an amino acid, can be used to infer the level of natural selection acting upon a gene. This conjecture is supported by a correlation between volatility and the rate of nonsynonymous substitution (dN), or the ratio of nonsynonymous and synonymous substitution rates, in a variety of organisms. These organisms include yeast, in which the correlations are quite strong. Here we show that these correlations are a by-product of a correlation between synonymous codon bias toward translationally optimal codons and dN. Although this analysis suggests that volatility is not a good measure of the selection, we suggest that it might be possible to infer something about the level of natural selection, from a single genome sequence, using translational codon bias.

Base Sequence↗

Sequence of PRAT satellite DNA "frozen" in some Coleopteran species.

The intriguing diversity of highly abundant satellite repeats found even among closely related species can result from processes leading to dramatic changes in copy number of a particular sequence in the genome and not from rapid accumulation of mutations. To test this hypothesis, we investigated the distribution of the PRAT satellite DNA family, a highly abundant major satellite in the coleopteran species Palorus ratzeburgii, in eight species belonging to the related genera ( Tribolium, Tenebrio, Latheticus), the subfamily (Pimeliinae), and the family (Chrysomelidae). Dot blot analysis and PCR assay followed by Southern hybridization revealed that the PRAT satellite, in the form of low-copy number repeats, was present in all tested species. The PRAT satellite detected in the species Pimelia elevata has been sequenced, and compared with previously cloned PRAT monomers from Palorus ratzeburgii and Palorus subdepressus. Although the two Palorus species diverged at least 7 Myr ago, and the subfamily Pimeliinae separated from the genus Palorus 50-60 Myr ago, all PRAT clones exhibit high mutual homology, with average variability relative to the common consensus sequence of 1.3%. The presence of ancestral mutations found in PRAT clones from all three species as well as the absence of species diagnostic mutations illustrate extremely slow sequence evolution. This unexpectedly high conservation of PRAT satellite DNA sequence might be induced by a small bias of turnover mechanisms favoring the ancestral sequence in the process of molecular drive.

Animals↗

A high-throughput Arabidopsis reverse genetics system.

A collection of Arabidopsis lines with T-DNA insertions in known sites was generated to increase the efficiency of functional genomics. A high-throughput modified thermal asymmetric interlaced (TAIL)-PCR protocol was developed and used to amplify DNA fragments flanking the T-DNA left borders from approximately 100000 transformed lines. A total of 85108 TAIL-PCR products from 52964 T-DNA lines were sequenced and compared with the Arabidopsis genome to determine the positions of T-DNAs in each line. Predicted T-DNA insertion sites, when mapped, showed a bias against predicted coding sequences. Predicted insertion mutations in genes of interest can be identified using Arabidopsis Gene Index name searches or by BLAST (Basic Local Alignment Search Tool) search. Insertions can be confirmed by simple PCR assays on individual lines. Predicted insertions were confirmed in 257 of 340 lines tested (76%). This resource has been named SAIL (Syngenta Arabidopsis Insertion Library) and is available to the scientific community at www.tmri.org.

Agrobacterium tumefaciens↗