SEARCH · PubMed Health
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.
Explore the source record for details and available documents.
Maternal phenotype and maternal environment can profoundly affect the phenotype and fitness of offspring. Yet the causes of variation in such maternal effects are rarely known. Embryos in avian eggs cannot develop without being incubated and this creates an opportunity for maternal control of duration and onset of offspring development. However, females might adjust the start of incubation (e.g., coincident with the first egg or delayed until after egg-laying) in response to environmental conditions that they experience at the time of breeding. We studied two populations of the house finch ( Carpodacus mexicanus) that breed at the climatic extremes of the species' geographical range (Montana and Alabama) and found that in both populations, the timing of incubation onset was closely associated with the bias in the sequence in which male and female eggs were laid within a clutch. When females started incubation with the first egg, they produced sons and daughters in highly biased sequence, when females delayed the onset of incubation until after the egg-laying, the sequence of sons and daughters was not biased. Because in both populations, onset of incubation was associated with the ambient temperature, these results emphasize that maternal effects on offspring can be influenced by ecological conditions experienced by parental generation.
We have sequenced the circumsporozoite protein gene from 16 isolates from China, the Philippines, Papua New Guinea and the Solomon Islands. We found very limited polymorphisms in the non-repetitive regions of the circumsporozoite gene from these isolates. All samples from China contained a 36-base insert 3' to the repeats previously seen only in a North Korean isolate. Limited variation was found in the repeat regions, which allowed these and previously sequenced isolates to be classified into groups based on repeat structure. These groupings also correlate with the geographical origin of the isolates.
Fifty-three gene sequences from E. coli containing 18,288 reading frame triplets have been characterized according to the nature and level of average codon preference. The distribution of average preferences is bimodal, with approximately half the genes using an average of only 36 codons, and the remainder just 42 codons. There is a high correlation between the level of codon bias, the tRNA population and the abundance of protein product, indicating biased patterns are exploited by the cell for the production of widely different levels of gene product. This relationship is especially striking in genes involved in the production of components for transcription and translation. Overall, the genes for these processes generate some five-fold more protein than the average in the genome, and use about five fewer codons. The very high codon bias found in the RNA polymerase gene thus provides a simple, autogenous mechanism for the coordinate synthesis of these components and RNA polymerase. A surprisingly high level of codon probability is also found in triplets of the complement of coding sequences. This is apparently due to the evolutionary dispersion of coding sequences and/or the requirement for increased levels of secondary structure in messenger RNAs.
The completion of the entire 4.2-Mb genome sequence of the gram-positive bacterium Bacillus subtilis has been a milestone for biological studies on this model organism. This paper describes bioinformatics work related to this joint European and Japanese project: methods and strategies for gene annotation and detection of sequencing errors, using an integrated cooperative computer environment (Imagene); construction of a specialized database for data management and a WWW server for data retrieval (SubtiList); DNA sequence analysis, yielding striking results on oligonucleotide bias, repeated sequences, and codon usage, all landmarks of evolutionary events shaping the B. subtilis genome.
Structural genetic alterations in cancer often involve gene loss or gene amplification. With the advent of microarray approaches for the analysis of the genome, as exemplified by array-CGH (Comparative Genomic Hybridization), scanning for gene-dosage alterations is limited only by issues of DNA microarray density. However, samples of interest to the pathologist often comprise small clusters of just a few hundred cells, which do not provide sufficient DNA for array-CGH analysis. We sought to develop a simple method that would permit amplification of the whole genome without the use of thermocycling or ligation of DNA adaptors, because such a method would lend itself to the automated processing of a large number of tissue samples. We describe a method that permits the isothermal amplification of genomic DNA with high fidelity and limited sequence representation bias. The method is based on strand displacement reactions that propagate by a hyperbranching mechanism, and generate hundreds, or even thousands, of copies of the genome in a few hours. Using whole genome isothermal amplification, in combination with comparative genomic hybridization on cDNA microarrays, we demonstrate the ability to detect gene losses in yeast and gene dosage imbalances in human breast tumor cell lines. Although sequence representation bias in the amplified DNA presents potential problems for CGH analysis, these problems have been overcome by using amplified DNA in both control and tester samples. Gene-dosage alterations of threefold or more can be observed with high reproducibility with as few as 1000 cells of starting material.
The Rat2 cell line carries 50-70 stably integrated copies per cell of a lambda/lacI shuttle vector as a target for mutagenicity testing. Rat2 cells were exposed to 1 and 10 micrograms/ml of 7,12-dimethylbenz[a]anthracene (DMBA) for 24 h at 37 degrees C in the presence of primary rat hepatocytes, and grown to confluence. The shuttle vector was rescued from untreated and mutagen-treated cells and mutant frequencies were determined. The low and high doses of DMBA induced mutant frequencies that were 7-fold (25 +/- 4.9 x 10(-5)) and 33-fold (127 +/- 19.9 x 10(-5)) higher, respectively, than the spontaneous mutant frequency (3.8 +/- 0.7 x 10(-5)). DNA sequence analysis of the DMBA-induced lacI- mutants indicated that they contained mainly basepair substitution mutations at A:T and G:C, and that A:T-->T:A and G:C-->T:A transversions were the predominant types. In addition, 23 of 28 (82%) A:T basepair substitution mutations occurred with the mutated dA, the putatively adducted base, on the coding strand. Furthermore, 20 of the 28 (71%) A:T mutations had the mutated dA flanked 5' by a dC, and 17 of these were A:T-->T:A transversions, suggesting a sequence preference for this mutation. Except for a higher proportion of G:C-->A:T transitions in the low dose data, the mutational profiles from low and high doses of DMBA were similar. These results indicate that DMBA mutagenesis in the lacI gene of Rat2 cells displays distinct DNA sequence and DNA strand preferences.
Phylogenetic relationships of 79 caniform carnivores were addressed based on four nuclear sequence-tagged sites (STS) and one nuclear exon, IRBP, using both supertree and supermatrix analyses. We recovered the three major arctoid lineages, Ursidae, Pinnipedia, and Musteloidea, as monophyletic, with Ursidae (bears) strongly supported as the basal arctoid lineage. Within Pinnipedia, Phocidae (true seals) were sister to the Otaroidea [Otariidae (fur seals and sea lions) and Odobenidae (walrus)]. Phocid subfamily and tribal designations were supported, but the otariid subfamily split between fur seals and sea lions was not. All family designations within Musteloidea were strongly supported: Mephitidae (skunks), Ailuridae (monotypic red panda), Mustelidae (weasels, badgers, otters), and Procyonidae (raccoons). A novel hypothesis for the position of the red panda was recovered, placing it as branching after Mephitidae and before Mustelidae+Procyonidae. Within Mustelidae, subfamily taxonomic changes are considered. This study represents the most comprehensive sampling to date of the Caniformia in a molecular study and contains the most complete molecular phylogeny for the Procyonidae. Our data set was also used in an empirical examination of the effect of missing data on both supertree and supermatrix analyses. Sequence for all genes in all taxa could not be obtained, so two variants of the data set with differing amounts of missing data were examined. The amount of missing data did not have a strong effect; instead, phylogenetic resolution was more dependent on the presence of sufficient informative characters. Supertree and supermatrix methods performed equivalently with incomplete data and were highly congruent; conflicts arose only in weakly supported areas, indicating that more informative characters are required to confidently resolve close species relationships.
1. There has been considerable interest in techniques recently developed for the study of arterial baroreceptor-cardiac reflex sensitivity based on analysis of spontaneous baroreflex sequences and on spectral analysis. This study examined how these newer techniques agreed with the established pharmacological methods in elderly subjects. 2. In 20 elderly subjects [10 hypertensive (clinic blood pressure 180 +/- 4/88 +/- 2 mmHg) and 10 normotensive (clinic blood pressure 136 +/- 3/73 +/- 2 mmHg)], we assessed baroreflex sensitivity from spontaneous sequences of increasing and decreasing blood pressure and pulse interval and their mean, and from spectral analysis to derive alpha, the index of overall baroreflex gain. Pharmacological baroreflex sensitivity was derived from the blood pressure and pulse interval responses to depressor (sodium nitroprusside) and pressor (phenylephrine) stimuli, and their mean. 3. Baroreflex sensitivity was significantly lower in the hypertensive group by the pharmacological, sequence and spectral methods (all P < 0.05). 4. There was acceptable agreement between pharmacological baroreflex sensitivity and sequences of the same direction, but with some systematic bias. There was also reasonable agreement between pharmacological and spectral baroreflex sensitivity and close agreement without bias between sequence and spectral methods. 5. The newer and established techniques demonstrate acceptable agreement in the elderly, albeit with some systematic bias. Pharmacological methods have enjoyed historical precedence but newer techniques give equivalent results, and are preferable in some circumstances. The newer techniques may be more descriptive of the spontaneous behaviour of the arterial baroreflex at rest rather than under artificially stimulated conditions.
It is now well-established that compositional bias in DNA sequences can adversely affect phylogenetic analysis based on those sequences. Phylogenetic analyses based on protein sequences are generally considered to be more reliable than those derived from the corresponding DNA sequences because it is believed that the use of encoded protein sequences circumvents the problems caused by nucleotide compositional biases in the DNA sequences. There exists, however, a correlation between AT/GC bias at the nucleotide level and content of AT- and GC-rich codons and their corresponding amino acids. Consequently, protein sequences can also be affected secondarily by nucleotide compositional bias. Here, we report that DNA bias not only may affect phylogenetic analysis based on DNA sequences, but also drives a protein bias which may affect analyses based on protein sequences. We present a striking example where common phylogenetic tools fail to recover the correct tree from complete animal mitochondrial protein-coding sequences. The data set is very extensive, containing several thousand sites per sequence, and the incorrect phylogenetic trees are statistically very well supported. Additionally, neither the use of the LogDet/paralinear transform nor removal of positions in the protein alignment with AT- or GC-rich codons allowed recovery of the correct tree. Two taxa with a large compositional bias continually group together in these analyses, despite a lack of close biological relatedness. We conclude that even protein-based phylogenetic trees may be misleading, and we advise caution in phylogenetic reconstruction using protein sequences, especially those that are compositionally biased.
A method has been developed for discovering patterns in DNA sequences. Loosely based on the well-known Lempel Ziv model for text compression, the model detects repeated sequences in DNA. The repeats can be forward or inverted, and they need not be exact. The method is particularly useful for detecting distantly related sequences, and for finding patterns in sequences of biased nucleotide composition, where spurious patterns are often observed because the bias leads to coincidental nucleotide matches. We show here the utility of the method by applying it to genomic sequences of Plasmodium falciparum. A single scan of chromosomes 2 and 3 of P. falciparum, using our method and no other a priori information about the sequences, reveals regions of low complexity in both telomeric and central regions, long repeats in the subtelomeric regions, and shorter repeat areas in dense coding regions. Application of the method to a recently sequenced contig of chromosome 10 that has a particularly biased base composition detects a long internal repeat more readily than does the conventional dot matrix plot. Space requirements are linear, so the method can be used on large sequences. The observed repeat patterns may be related to large-scale chromosomal organization and control of gene expression. The method has general application in detecting patterns of potential interest in newly sequenced genomic material.
MOTIVATION: Hidden Markov models (HMMs) calculate the probability that a sequence was generated by a given model. Log-odds scoring provides a context for evaluating this probability, by considering it in relation to a null hypothesis. We have found that using a reverse-sequence null model effectively removes biases owing to sequence length and composition and reduces the number of false positives in a database search. Any scoring system is an arbitrary measure of the quality of database matches. Significance estimates of scores are essential, because they eliminate model- and method-dependent scaling factors, and because they quantify the importance of each match. Accurate computation of the significance of reverse-sequence null model scores presents a problem, because the scores do not fit the extreme-value (Gumbel) distribution commonly used to estimate HMM scores' significance. RESULTS: To get a better estimate of the significance of reverse-sequence null model scores, we derive a theoretical distribution based on the assumption of a Gumbel distribution for raw HMM scores and compare estimates based on this and other distribution families. We derive estimation methods for the parameters of the distributions based on maximum likelihood and on moment matching (least-squares fit for Student's t-distribution). We evaluate the modeled distributions of scores, based on how well they fit the tail of the observed distribution for data not used in the fitting and on the effects of the improved E-values on our HMM-based fold-recognition methods. The theoretical distribution provides some improvement in fitting the tail and in providing fewer false positives in the fold-recognition test. An ad hoc distribution based on assuming a stretched exponential tail does an even better job. The use of Student's t to model the distribution fits well in the middle of the distribution, but provides too heavy a tail. The moment-matching methods fit the tails better than maximum-likelihood methods. AVAILABILITY: Information on obtaining the SAM program suite (free for academic use), as well as a server interface, is available at http://www.soe.ucsc.edu/research/compbio/sam.html and the open-source random sequence generator with varying compositional biases is available at http://www.soe.ucsc.edu/research/compbio/gen_sequence
The relationship between the codon usage bias and the sequence context surrounding the AUG translation initiation codon was examined in 1100 Drosophila melanogaster mRNA sequences. The codon usage bias measured by the "codon adaptation index" (CAI), and the effectiveness of the AUG context for translation initiation assessed by the "AUG context adaptation index" (AUGCAI), showed a significant positive relationship (correlation coefficient: r = 0.34, p <0.0001), indicating that these two factors are evolutionally under a similar natural selection constraint at the translational level. The importance of each position of the AUG context in relation to codon usage bias was examined, and the preference for the nucleotide at the -13, -12, -11, -10, -7, -6, -5, -4, -3, -2, and -1 positions showed a significant positive correlation to the codon usage bias, suggesting the action of natural selection on these very specific positions of the Drosophila genome. The relationship between AUGCAI value and gene length was also examined, and a significant negative relationship was found (r = -0.15, p <0.0001), suggesting a general tendency of higher expressivity of shorter genes, and of lower expressivity of longer genes in D. melanogaster.
The relationship between the codon usage bias and the sequence context surrounding the AUG translation initiation codon was examined in 211 Saccharomyces cerevisiae mRNA sequences. The codon usage bias and the number of matches to optimal AUG context, (A/U)A(A/C)AA(A/C)AUGUC(U/C), for translation initiation showed a positive relationship, indicating that these two factors are evolutionally under the similar natural selection constraint at the translation level. A new index (AUGCAI = AUG Context Adaptation Index) for the measure of optimal AUG context was devised, and the importance of each position of AUG context was also examined.
The complete sequence of 57 kb of the human HPRT locus has been determined using automated fluorescent DNA sequencing. The strategy employed increasingly directed sequencing methods: A randomly generated M13 library was sequenced to generate contiguous overlapping sets of sequences (contigs). M13 clones at the ends of these contigs were further sequenced using M13 (universal and reverse) and custom oligonucleotide primers to order the contigs and to complete the sequencing project. The human HPRT sequence includes 1676 bp 5' and 15,238 bp 3' to exons 1 and 9, respectively. The sequence contains 49 representatives of the Alu repeat, along with several other types of repetitive sequences. The Alu sequences exhibit a biased orientation, with those sequences in the first half of the locus oriented in the minus direction relative to transcription of the gene (3'----5' = 77%, P less than 0.005) and those sequences in the latter half of the locus oriented randomly (5'----3' = 67%, P less than 0.5). The development and performance of the sequencing strategy and the features of the human HPRT gene are presented.
The effect of different amplification and labeling methods on DNA microarray expression results has not been previously delineated. To analyze the variation associated with widely accepted T7-based RNA amplificationand labeling methods, aliquots of the Stratagene Human Universal Reference RNA were labeled using three eukaryotic target preparation methods followed by uniform replicate array hybridization (Affymetrix U95Av2). Method-dependent variability was observed in the yield and size distribution of labeled products, as well as in the gene expression results. A significant increase in short transcripts, when compared to unamplified mRNA, was observed in methods with long in vitro transcription reactions. Intramethod reproducibility showed correlation coefficients >0.99, whereas intermethod comparisons showed coefficients ranging from 0.94 to 0.98 and a nearly twofold increase in coefficient of variation. Fold amplification for each method positively correlated with the number of genes present. Our experiments uncovered two factors that introduced significant bias in gene expression data: the number of labeled nucleotides, which introduces sequence-dependent bias, and the length of the in vitro transcription reaction, which introduces transcript size-dependent bias. This study provides evidence that variability in expression data may be caused, in part, by differences in amplification and labeling protocols.
Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.
It is presently accepted that, in mammals, due to the greater number of cell divisions in the male germline than in the female germline, nucleotide substitutions occur more frequently in males. The data on mutation bias in insertions and deletions (indels) are contradictory, with some studies indicating no sex bias and others indicating either female or male bias. The sequenced rat and mouse genomes provide a unique opportunity to investigate a potential sex bias for different types of mutations. Indeed, mutation rates can be accurately estimated from a large number of orthologous loci in organisms similar in generation time and in the number of germline cell divisions. Here we compare the mutation rates between chromosome X and autosomes for likely neutral sites in eutherian ancestral interspersed repetitive elements present at orthologous locations in the rat and mouse genomes. We find that small indels are male biased: The male-to-female mutation rate ratio (alpha) for indels in rodents is approximately 2. Similarly, our whole-genome analysis in rodents indicates an approximately twofold excess of nucleotide substitutions originating in males over that in females. This is the same as the male-to-female ratio of the number of germline cell divisions in rat and mouse. Thus, this is consistent with nucleotide substitutions and small indels occurring primarily during DNA replication.