Biomedical subjects

David H Ardell

Publications and source records attributed to David H Ardell.

12 recordsLinked to original sources

PubMed2006

The coevolution of genes and genetic codes: Crick's frozen accident revisited.

The standard genetic code is the nearly universal system for the translation of genes into proteins. The code exhibits two salient structural characteristics: it possesses a distinct organization that makes it extremely robust to errors in replication and translation, and it is highly redundant. The origin of these properties has intrigued researchers since the code was first discovered. One suggestion, which is the subject of this review, is that the code's organization is the outcome of the coevolution of genes and genetic codes. In 1968, Francis Crick explored the possible implications of coevolution at different stages of code evolution. Although he argues that coevolution was likely to influence the evolution of the code, he concludes that it falls short of explaining the organization of the code we see today. The recent application of mathematical modeling to study the effects of errors on the course of coevolution, suggests a different conclusion. It shows that coevolution readily generates genetic codes that are highly redundant and similar in their error-correcting organization to the standard code. We review this recent work and suggest that further affirmation of the role of coevolution can be attained by investigating the extent to which the outcome of coevolution is robust to other influences that were present during the evolution of the code.

Animals↗

PubMed2006

TFAM detects co-evolution of tRNA identity rules with lateral transfer of histidyl-tRNA synthetase.

We present TFAM, an automated, statistical method to classify the identity of tRNAs. TFAM, currently optimized for bacteria, classifies initiator tRNAs and predicts the charging identity of both typical and atypical tRNAs such as suppressors with high confidence. We show statistical evidence for extensive variation in tRNA identity determinants among bacterial genomes due to variation in overall tDNA base content. With TFAM we have detected the first case of eukaryotic-like tRNA identity rules in bacteria. An alpha-proteobacterial clade encompassing Rhizobiales, Caulobacter crescentus and Silicibacter pomeroyi, unlike a sister clade containing the Rickettsiales, Zymomonas mobilis and Gluconobacter oxydans, uses the eukaryotic identity element A73 instead of the highly conserved prokaryotic element C73. We confirm divergence of bacterial histidylation rules by demonstrating perfect covariation of alpha-proteobacterial tRNA(His) acceptor stems and residues in the motif IIb tRNA-binding pocket of their histidyl-tRNA synthetases (HisRS). Phylogenomic analysis supports lateral transfer of a eukaryotic-like HisRS into the alpha-proteobacteria followed by in situ adaptation of the bacterial tDNA(His) and identity rule divergence. Our results demonstrate that TFAM is an effective tool for the bioinformatics, comparative genomics and evolutionary study of tRNA identity.

Alphaproteobacteria↗

PubMed2006

Visualizing bacterial tRNA identity determinants and antideterminants using function logos and inverse function logos.

Sequence logos are stacked bar graphs that generalize the notion of consensus sequence. They employ entropy statistics very effectively to display variation in a structural alignment of sequences of a common function, while emphasizing its over-represented features. Yet sequence logos cannot display features that distinguish functional subclasses within a structurally related superfamily nor do they display under-represented features. We introduce two extensions to address these needs: function logos and inverse logos. Function logos display subfunctions that are over-represented among sequences carrying a specific feature. Inverse logos generalize both sequence logos and function logos by displaying under-represented, rather than over-represented, features or functions in structural alignments. To make inverse logos, a compositional inverse is applied to the feature or function frequency distributions before logo construction, where a compositional inverse is a mathematical transform that makes common features or functions rare and vice versa. We applied these methods to a database of structurally aligned bacterial tDNAs to create highly condensed, birds-eye views of potentially all so-called identity determinants and antideterminants that confer specific amino acid charging or initiator function on tRNAs in bacteria. We recovered both known and a few potentially novel identity elements. Function logos and inverse logos are useful tools for exploratory bioinformatic analysis of structure-function relationships in sequence families and superfamilies.

Base Sequence↗

PubMed2006

Influences on gene expression in vivo by a Shine-Dalgarno sequence.

The Shine-Dalgarno (SD+: 5'-AAGGAGG-3') sequence anchors the mRNA by base pairing to the 16S rRNA in the small ribosomal subunit during translation initiation. We have here compared how an SD+ sequence influences gene expression, if located upstream or downstream of an initiation codon. The positive effect of an upstream SD+ is confirmed. A downstream SD+ gives decreased gene expression. This effect is also valid for appropriately modified natural Escherichia coli genes. If an SD+ is placed between two potential initiation codons, initiation takes place predominantly at the second start site. The first start site is activated if the distance between this site and the downstream SD+ is enlarged and/or if the second start site is weakened. Upstream initiation is eliminated if a stable stem-loop structure is placed between this SD+ and the upstream start site. The results suggest that the two start sites compete for ribosomes that bind to an SD+ located between them. A minor positive contribution to upstream initiation resulting from 3' to 5' ribosomal diffusion along the mRNA is suggested. Analysis of the E. coli K12 genome suggests that the SD+ or SD-like sequences are systematically avoided in the early coding region suggesting an evolutionary significance.

Base Sequence↗

PubMed2005

The length of the 5' leader of Escherichia coli tRNA precursors influences bacterial growth.

Based on a computational analysis of the 5' regions of tRNA-encoding genes, the average length of the 5' leaders in tRNA precursors in Escherichia coli appears to be 17-18 residues long. An in vivo assay based on tRNA nonsense suppression was developed and used to investigate the function of the 5' leader of the tRNA precursors on tRNA processing and bacterial growth. Our data indicate that the 5' leader influences bacterial growth but is surprisingly not absolutely necessary for growth. These findings are consistent with previous in vitro data where it was demonstrated that the 5' leader plays a role in the interaction with RNase P, the endoribonuclease responsible for removing the 5' leader in the cell. We discuss the plausible role of the 5' leader in processing and tRNA gene expression.

Base Sequence↗

PubMed2005

The genomic pattern of tDNA operon expression in E. coli.

In fast-growing microorganisms, a tRNA concentration profile enriched in major isoacceptors selects for the biased usage of cognate codons. This optimizes translational rate for the least mass invested in the translational apparatus. Such translational streamlining is thought to be growth-regulated, but its genetic basis is poorly understood. First, we found in reanalysis of the E. coli tRNA profile that the degree to which it is translationally streamlined is nearly invariant with growth rate. Then, using least squares multiple regression, we partitioned tRNA isoacceptor pools to predicted tDNA operons from the E. coli K12 genome. Co-expression of tDNAs in operons explains the tRNA profile significantly better than tDNA gene dosage alone. Also, operon expression increases significantly with proximity to the origin of replication, oriC, at all growth rates. Genome location explains about 15% of expression variation in a form, at a given growth rate, that is consistent with replication-dependent gene concentration effects. Yet the change in the tRNA profile with growth rate is less than would be expected from such effects. We estimated per-copy expression rates for all tDNA operons that were consistent with independent estimates for rDNA operons. We also found that tDNA operon location, and the location dependence of expression, were significantly different in the leading and lagging strands. The operonic organization and genomic location of tDNA operons are significant factors influencing their expression. Nonrandom patterns of location and strandedness shown by tDNA operons in E. coli suggest that their genomic architecture may be under selection to satisfy physiological demand for tRNA expression at high growth rates.

Journal Article↗

PubMed2004

The louse-borne human pathogen Bartonella quintana is a genomic derivative of the zoonotic agent Bartonella henselae.

We present the complete genomes of two human pathogens, Bartonella quintana (1,581,384 bp) and Bartonella henselae (1,931,047 bp). The two pathogens maintain several similarities in being transmitted by insect vectors, using mammalian reservoirs, infecting similar cell types (endothelial cells and erythrocytes) and causing vasculoproliferative changes in immunocompromised hosts. A primary difference between the two pathogens is their reservoir ecology. Whereas B. quintana is a specialist, using only the human as a reservoir, B. henselae is more promiscuous and is frequently isolated from both cats and humans. Genome comparison elucidated a high degree of overall similarity with major differences being B. henselae specific genomic islands coding for filamentous hemagglutinin, and evidence of extensive genome reduction in B. quintana, reminiscent of that found in Rickettsia prowazekii. Both genomes are reduced versions of chromosome I from the highly related pathogen Brucella melitensis. Flanked by two rRNA operons is a segment with similarity to genes located on chromosome II of B. melitensis, suggesting that it was acquired by integration of megareplicon DNA in a common ancestor of the two Bartonella species. Comparisons of the vector-host ecology of these organisms suggest that the utilization of host-restricted vectors is associated with accelerated rates of genome degradation and may explain why human pathogens transmitted by specialist vectors are outnumbered by zoonotic agents, which use vectors of broad host ranges.

Animals↗

PubMed2004

SCANMS: adjusting for multiple comparisons in sliding window neutrality tests.

SUMMARY: A Perl utility called SCANMS was written that operates on parametric bootstrap data generated by Hudson's simulator MS. SCANMS simulates sliding window analysis of Tajima's D or Fu and Li's D(*) or F(*) over the replicates and generates the empirical joint distribution of window minima and maxima. This distribution may then be used to calculate achieved significance values or significance cutoffs that account for multiple comparisons in sliding window analysis of real data. AVAILABILITY: SCANMS is available free of charge under the terms of the Academic Free License 2.0 and may be downloaded from http://www.lcb.uu.se/~dave/SCANMS

Algorithms↗

PubMed2003

Three neuropeptide Y receptor genes in the spiny dogfish, Squalus acanthias, support en bloc duplications in early vertebrate evolution.

It has been debated whether the increase in gene number during early vertebrate evolution was due to multiple independent gene duplications or synchronous duplications of many genes. We describe here the cloning of three neuropeptide Y (NPY) receptor genes belonging to the Y1 subfamily in the spiny dogfish, Squalus acanthias, a cartilaginous fish. The three genes are orthologs of the mammalian subtypes Y1, Y4, and Y6, which are located in paralogous gene regions on different chromosomes in mammals. Thus, these genes arose by duplications of a chromosome region before the radiation of gnathostomes (jawed vertebrates). Estimates of duplication times from linearized trees together with evidence from other gene families supports two rounds of chromosome duplications or tetraploidizations early in vertebrate evolution. The anatomical distribution of mRNA was determined by reverse-transcriptase PCR and was found to differ from mammals, suggesting differential functional diversification of the new gene copies during the radiation of the vertebrate classes.

Amino Acid Sequence↗

PubMed2003

Polymorphism, recombination and alternative unscrambling in the DNA polymerase alpha gene of the ciliate Stylonychia lemnae (Alveolata; class Spirotrichea).

DNA polymerase alpha is the most highly scrambled gene known in stichotrichous ciliates. In its hereditary micronuclear form, it is broken into >40 pieces on two loci at least 3 kb apart. Scrambled genes must be reassembled through developmental DNA rearrangements to yield functioning macronuclear genes, but the mechanism and accuracy of this process are unknown. We describe the first analysis of DNA polymorphism in the macronuclear version of any scrambled gene. Six functional haplotypes obtained from five Eurasian strains of Stylonychia lemnae were highly polymorphic compared to Drosophila genes. Another incompletely unscrambled haplotype was interrupted by frameshift and nonsense mutations but contained more silent mutations than expected by allelic inactivation. In our sample, nucleotide diversity and recombination signals were unexpectedly high within a region encompassing the boundary of the two micronuclear loci. From this and other evidence we infer that both members of a long repeat at the ends of the loci provide alternative substrates for unscrambling in this region. Incongruent genealogies and recombination patterns were also consistent with separation of the two loci by a large genetic distance. Our results suggest that ciliate developmental DNA rearrangements may be more probabilistic and error prone than previously appreciated and constitute a potential source of macronuclear variation. From this perspective we introduce the nonsense-suppression hypothesis for the evolution of ciliate altered genetic codes. We also introduce methods and software to calculate the likelihood of hemizygosity in ciliate haplotype samples and to correct for multiple comparisons in sliding-window analyses of Tajima's D.

Animals↗

PubMed2002

No accident: genetic codes freeze in error-correcting patterns of the standard genetic code.

The standard genetic code poses a challenge in understanding the evolution of information processing at a fundamental level of biological organization. Genetic codes are generally coadapted with, or 'frozen' by, the protein-coding genes that they translate, and so cannot easily change by natural selection. Yet the standard code has a significantly non-random pattern that corrects common errors in the transmission of information in protein-coding genes. Because of the freezing effect and for other reasons, this pattern has been proposed not to be due to selection but rather to be incidental to other evolutionary forces or even entirely accidental. We present results from a deterministic population genetic model of code-message coevolution. We explicitly represent the freezing effect of genes on genetic codes and the perturbative effect of changes in genetic codes on genes. We incorporate characteristic patterns of mutation and translational error, namely, transition bias and positional asymmetry, respectively. Repeated selection over small successive changes produces genetic codes that are substantially, but not optimally, error correcting. In particular, our model reproduces the error-correcting patterns of the standard genetic code. Aspects of our model and results may be applicable to the general problem of adaptation to error in other natural information-processing systems.

Animals↗

PubMed2002

The impact of message mutation on the fitness of a genetic code.

The Standard Genetic Code is organized such that similar codons encode similar amino acids. One explanation suggested that the Standard Code is the result of natural selection to reduce the fitness "load" that derives from the mutation and mistranslation of protein-coding genes. We review the arguments against the mutational load-minimizing hypothesis and argue that they need to be reassessed. We review recent analyses of the organization of the Standard Code and conclude that under cautious interpretation they support the mutational load-minimizing hypothesis. We then present a deterministic asexual model with which we study the mode of selection for load minimization. In this model, individual fitness is determined by a protein phenotype resulting from the translation of a mutable set of protein-coding genes. We show that an equilibrium fitness may be associated with a population with the same genetic code and that genetic codes that assign similar codons to similar amino acids have a higher fitness. We also show that the number of mutant codons in each individual at equilibrium, which determines the strength of selection for load minimization, reflects a long-term evolutionary balance between mutations in messages and selection on proteins, rather than the number of mutations that occur in a single generation, as has been assumed by previous authors. We thereby establish that selection for mutational load minimization acts at the level of an individual in a single generation. We conclude with comments on the shortcomings and advantages of load minimization over other hypotheses for the origin of the Standard Code.

Amino Acids↗