Identification of a promoter element within the zebrafish colXalpha1 gene responsive to runx2 isoforms Osf2/Cbfa1 and til-1 but not to pebp2alphaA2.

Type X collagen is a short chain collagen specifically expressed by hypertrophic chondrocytes during endochondral ossification. We report here the functional analysis of the zebrafish (Danio rerio) collagen Xalpha1 gene (colXalpha1) promoter with the identification of a region responsive to two isoforms of the runt domain transcription factor runx2. Furthermore, we provide evidence for the presence of dual promoter usage in zebrafish, a finding that should be important to further understanding of the regulation of its restricted tissue distribution and spatial-temporal expression during early development. The zebrafish colXalpha1 gene structure is comparable to that recently identified by comparative genomics in takifugu and shows homology with corresponding mammalian genes, indicating that its general architecture has been maintained throughout vertebrate evolution. Our data suggest that, as in mammals, runx2 plays a role in the development of the osteogenic lineage, supporting zebrafish as a model for studies of bone and cartilage development.

Can codon usage bias explain intron phase distributions and exon symmetry?

More introns exist between codons (phase 0) than between the first and the second bases (phase 1) or between the second and the third base (phase 2) within the codon. Many explanations have been suggested for this excess of phase 0. It has, for example, been argued to reflect an ancient utility for introns in separating exons that code for separate protein modules. There may, however, be a simple, alternative explanation. Introns typically require, for correct splicing, particular nucleotides immediately 5' in exons (typically a G) and immediately 3' in the following exon (also often a G). Introns therefore tend to be found between particular nucleotide pairs (e.g., G|G pairs) in the coding sequence. If, owing to bias in usage of different codons, these pairs are especially common at phase 0, then intron phase biases may have a trivial explanation. Here we take codon usage frequencies for a variety of eukaryotes and use these to generate random sequences. We then ask about the phase of putative intron insertion sites. Importantly, in all simulated data sets intron phase distribution is biased in favor of phase 0. In many cases the bias is of the magnitude observed in real data and can be attributed to codon usage bias. It is also known that exons may carry either the same phase (symmetric) or different phases (asymmetric) at the opposite ends. We simulated a distribution of different types of exons using frequencies of introns observed in real genes assuming random combination of intron phases at the opposite sides of exons. Surprisingly the simulated pattern was quite similar to that observed. In the simulants we typically observe a prevalence of symmetric exons carrying phase 0 at both ends, which is common for eukaryotic genes. However, at least in some species, the extent of the bias in favor of symmetric (0,0) exons is not as great in simulants as in real genes. These results emphasize the need to construct a biologically relevant null model of successful intron insertion.

Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection.

In numerous species, from bacteria to Drosophila, evidence suggests that selection acts even on synonymous codon usage: codon bias is greater in more abundantly expressed genes, the rate of synonymous evolution is lower in genes with greater codon bias, and there is consistency between genes in the same species in which codons are preferred. In contrast, in mammals, while nonequal use of alternative codons is observed, the bias is attributed to the background variance in nucleotide concentrations, reflected in the similar nucleotide composition of flanking noncoding and exonic third sites. However, a systematic examination of the covariants of codon usage controlling for background nucleotide content has yet to be performed. Here we present a new method to measure codon bias that corrects for background nucleotide content and apply this to 2396 human genes. Nearly all (99%) exhibit a higher amount of codon bias than expected by chance. The patterns associated with selectively driven codon bias are weakly recovered: Broadly expressed genes have a higher level of bias than do tissue-specific genes, the bias is higher for genes with lower rates of synonymous substitutions, and certain codons are repeatedly preferred. However, while these patterns are suggestive, the first two patterns appear to be methodological artifacts. The last pattern reflects in part biases in usage of nucleotide pairs. We conclude that we find no evidence for selection on codon usage in humans.

Local similarity in evolutionary rates extends over whole chromosomes in human-rodent and mouse-rat comparisons: implications for understanding the mechanistic basis of the male mutation bias.

The sex chromosomes and autosomes spend different times in the germ line of the two sexes. If cell division is mutagenic and if the sexes differ in number of cell divisions, then we expect that sequences on the X and Y chromosomes and autosomes should mutate at different rates. Tests of this hypothesis for several mammalian species have led to conflicting results. At the same time, recent evidence suggests that the chromosomal location of genes on autosomes affects their rate of evolution at synonymous sites. This suggests a mutagenic source different from germ cell replication. To correctly interpret the previous estimates of male mutation bias, it is crucial to understand the degree and range of this local similarity. With a carefully chosen randomization protocol, local similarity in synonymous rates of evolution can be detected in human-rodent and mouse-rat comparisons. However, the synonymous-site similarity in the mouse-rat comparison remains weak. Simulations suggest that this difference between the mouse-human and the mouse-rat comparisons is not artifactual and that there is therefore a difference between humans and rodents in the local patterns of mutation or selection on synonymous sites (conversely, we show that the previously reported absence of a local similarity in nonsynonymous rates of evolution in the human-rodent comparison was a methodological artifact). We show that linkage effects have a long-range component: not one in a million random genomes shows such levels of autosomal heterogeneity. The heterogeneity is so great that more autosomes than expected by chance have rates of synonymous evolution comparable with that of the X chromosome. As autosomal heterogeneity cannot be owing to different times spent in the germ line, this demonstrates that the dominant determiner of synonymous rates of evolution is not, as has been conjectured, the time spent in the male germ line.

The evolution of isochores.

One of the most striking features of mammalian chromosomes is the variation in G+C content that occurs over scales of hundreds of kilobases to megabases, the so-called 'isochore' structure of the human genome. This variation in base composition affects both coding and non-coding sequences and seems to reflect a fundamental level of genome organization. However, although we have known about isochores for over 25 years, we still have a poor understanding of why they exist. In this article, we review the current evidence for the three main hypotheses.

A comparative test of a theory for the evolution of anisogamy.

Why are sperm small and eggs large? The dominant explanation for the evolution of gamete size dimorphism envisages two opposing selection pressures acting on gamete size: small gametes are favoured because many can be produced, whereas large gametes contribute to a large zygote with consequently increased survival chances. This model predicts disruptive selection on gamete size (i.e. selection for anisogamy) if increases in zygote size confer disproportional increases in fitness (at least over part of its size range). It therefore predicts that increases in adult size should be accompanied by stronger selection for anisogamy. Using data from the green algal order Volvocales, we provide the first phylogenetically controlled test of the model's predictions using a published phylogeny and a new phylogeny derived by a different method. The predictions that larger organisms should (i) have a greater degree of gamete dimorphism and (ii) have larger eggs are broadly upheld. However, the results are highly sensitive to the phylogeny and the mode of analysis used.

High guanine-cytosine content is not an adaptation to high temperature: a comparative analysis amongst prokaryotes.

The causes of the variation between genomes in their guanine (G) and cytosine (C) content is one of the central issues in evolutionary genomics. The thermal adaptation hypothesis conjectures that, as G:C pairs in DNA are more thermally stable than adenonine:thymine pairs, high GC content may he a selective response to high temperature. A compilation of data on genomic GC content and optimal growth temperature for numerous prokaryotes failed to demonstrate the predicted correlation. By contrast, the GC content of Structural RNAs is higher at high temperatures. The issue that we address here is whether more freely evolving sites in exons (i.e. codonic third positions) evolve in the same manner as genomic DNA as a whole, Showing no correlated response, or like structural RNAs showing a strong correlation. The latter pattern would provide strong support for the thermal adaptation hypothesis, as the variation in GC content between orthologous genes is typically most profoundly seen at codon third sites (GC3). Simple analysis of completely sequenced prokaryotic genomes shows that GC3, but not genomic GC, is higher on average in thermophilic species. This demonstrates, if nothing else, that the results from the two measures cannot be presumed to be the same. A proper analysis, however, requires phylogenetic control. Here, therefore, we report the results of a comparative analysis of GC composition and optimal growth temperature for over 100 prokaryotes. Comparative analysis fails to show, in either Archea or Eubacteria, any hint of connection between optimal growth temperature and GC content in the genome as a whole, in protein-coding regions or, more crucially at GC. Conversely, comparable analysis confirms that GC content of structural RNA is strongly correlated with optimal temperature. Against the expectations of the thermal adaptation hypothesis, within prokaryotes GC content in protein-coding genies, even at relatively freely evolving sites, cannot be considered an adaptation to the thermal environment.

The elevated GC content at exonic third sites is not evidence against neutralist models of isochore evolution.

The human genome is divided into isochores, large stretches (>>300 kb) of genomic DNA with more or less consistent GC content. Mutational/neutralist and selectionist models have been put forward to explain their existence. A major criticism of the mutational models is that they cannot account for the higher GC content at fourfold-redundant silent sites within exons (GC4) than in flanking introns (GCi). Indeed, it has been asserted that it is hard to envisage a mutational bias explanation, as it is difficult to see how repair enzymes might act differently in exons and their flanking introns. However, this rejection, we note, ignores the effects of transposable elements (TEs), which are a major component of introns and tend to cause them to have a GC content different from (usually lower than) that dictated by point mutational processes alone. As TEs tend not to insert at the extremities of introns, this model predicts that GC content at the extremities of introns should be more like that at GC4 than are the intronic interiors. This we show to be true. The model also correctly predicts that small introns should have a composition more like that at GC4 than large introns. We conclude that the logic of the previous rejection of neutralist models is unsafe.

Evidence for purifying selection acting on silent sites in BRCA1.

In mammals, it is usually assumed that selection cannot be strong enough to act on nucleotide mutations that do not cause a change at the protein level (i.e. 'silent' or 'synonymous' mutations). Here we report the results of a molecular evolutionary analysis of BRCA1. We find a repeatable pronounced peak in the ratio of nonsynonymous to synonymous substitutions between codons 200-300. Unusually, this peak is caused by a plummet in the silent-site rate of evolution. The most parsimonious interpretation of these data is that purifying selection is acting on silent sites.

Accelerated molecular evolution of insect orthologues of ERG28/C14orf1: a link with ecdysteroid metabolism?

We have analysed the evolution of ERG28/C14orf1, a gene coding for a protein involved in sterol biosynthesis. While primary sequence of the protein is well conserved in all organisms able to synthesize sterols de novo, strong divergence is noticed in insects, which are cholesterol auxotrophs. In spite of this virtual acceleration, our analysis suggests that the insect orthologues are evolving today at rates similar to those of the remaining members of the family. A plausible way to explain this acceleration and subsequent stabilization is that Erg28 plays a role in at least two different pathways. Discontinuation of the cholesterogenesis pathway in insects allowed the protein to evolve as much as the function in the other pathway was not compromised.

The proteins of linked genes evolve at similar rates.

Much more variation in the rate of protein evolution occurs than is expected by chance. But why some proteins evolve rapidly but others slowly is poorly resolved. It was proposed, for example, that essential genes might evolve slower than dispensable ones, but this is not the case; and despite earlier claims, rates of evolution do not correlate with amino-acid composition. A few patterns have been found: proteins involved in antagonistic co-evolution (for example, immune genes, parasite antigens and reproductive conflict genes) tend to be rapidly evolving, and there is a correlation between the rate of protein evolution and the mutation rate of the gene. Here we report a new highly statistically significant predictor of a protein's rate of evolution, and show that linked genes have similar rates of protein evolution. There is also a weaker similarity of rates of silent site evolution (see ref. 13), which appears to be, in part, a consequence of the similarity in rates of protein evolution. The similarity in rates of protein evolution is not a consequence of underlying mutational patterns. A pronounced negative correlation between the rate of protein evolution and a covariant of the recombination rate indicates that rates of protein evolution possibly reflect, in part, the local strength of stabilizing selection.

Dosage, deletions and dominance: simple models of the evolution of gene expression.

Dominance of the wild-type allele over spontaneous null mutations, such as deletions, can be explained in terms of the effects of changes in enzyme dose on the flux of metabolic pathways. If ever increasing levels of enzyme activity have ever decreasing effects on the flux of the biochemical pathway, then halving of dosage will always have a lesser effect on flux than half the effect of complete removal of gene activity. Furthermore, if gene expression rates are high, then halving of dose can have a negligible effect on flux and dominance will be strong. Given that strong dominance appears to be common, this leaves open the issue of why enzyme activity levels are so high that a halving of expression rates is of minimal effect. Why produce so much surplus enzyme? One explanation, suggested by Haldane, is that selection favoured high expression levels as a defence against mutation. We model this scenario formally and show that protection from mutation is an extremely weak force determining expression levels. The selective coefficients are only of the order of the mutation rate. However, if we suppose a linear mapping of flux with fitness and a monotonic cost to increased gene expression, it follows simply that here exists an optimal level of gene expression. By contrast to the mutational model, doubling of gene expression rates when the system is distant from the optimum is associated with extremely high selective coefficients (orders of magnitude higher than the mutation rate). When the cost of gene expression is slight the optimal rate of expression is such that strong dominance will follow.

The molecular evolution of signal peptides.

Signal peptides direct mature peptides to their appropriate cellular location, after which they are cleaved off. Very many random alternatives can serve the same function. Of all coding sequences, therefore, signal peptides might come closest to being neutrally evolving. Here we consider this issue by examining the molecular evolution of 76 mouse-rat orthologues, each with defined signal peptides. Although they do evolve rapidly, they evolve about half as fast as neutral sequences. This indicates that a substantial proportion of mutations must be under stabilizing selection. A few putative signal sequences lack a hydrophobic core and these tend to be more slowly evolving than others, indicating even stronger stabilizing selection. However, closer scrutiny suggests that some of these represent mis-annotations in GenBank. It is also likely that some of the substitutions are not neutral. We find, for example, that the rate of protein evolution correlates with that of the mature peptide. This may be a result of compensatory evolution. We also find that signal peptides of immune genes tend to be faster evolving than the average, which suggests an association with antagonistic co-evolution. Previous reports also indicated that the signal peptide of the imprinted gene, Igf2r, is also unusually fast evolving. This, it was hypothesized, might also be indicative of antagonistic co-evolution. Comparison of Igf2r's signal peptide evolution shows that, although it is not an outlier, its rate of evolution is comparable to that of many of the faster evolving immune system signal sequences and 5/6 of the amino acid changes do not conserve hydrophobicity. This is at least suggestive that there is something unusual about Igf2r's signal sequence.

The evolution of gene number: are heritable and non-heritable errors equally important?

Is there a limit to the number of genes carried by an organism? Two reasons have been. First, as most mutations are deleterious, for a given per locus mutation rate there must exist an upper limit to the number of genes that is consistent with individual survival. Second, the imprecision of the mechanisms governing gene expression might also restrict genomic complexity. As gene expression errors are probably much more common than mutations, it is the latter that are more likely to impose a limit. However, these errors are not heritable and therefore cannot accumulate in populations. Which of the two sorts of effect are more likely to impose a limit? We address this issue in two ways. First, we ask about the load imposed by each sort of error. We show that the harmful effect of non-heritable failures is higher than that of heritable mutations, if (p) x (delta) > mu, where p is the rate of non-heritable failures, delta measures the harmful effect of these failures and mu is the rate of heritable mutations. Therefore, although the rate of non-heritable errors might be very high, this does not demonstrate that they are more important than mutations as their impact must be discounted by the strength of their effects. Further, we note that both theory and evidence suggest that the most common errors are of the least importance. Second, we discuss the population genetics of a new gene duplication. Previous attempts to make a connection between error rates and limits on gene number are based on group selection arguments. These fail to show a direct limitation on the spread of gene duplications. We note that empirical evidence indicates that duplication per se tends to result in expression errors that may be heritable. We therefore argue that a hybrid model, one evoking heritable expression errors, is likely to be the most realistic.

Male killing can select for male mate choice: a novel solution to the paradox of the lek.

In lekking species, intense directional selection is applied to aspects of the male genotype by female choice. Under conventional quantitative genetics theory, the expectation is that this will lead to a rapid loss in additive genetic variance for the trait in question. However, despite female choice, male variation is maintained and hence it pays females to continue choosing. This has been termed the 'paradox of the lek'. Here we present a theoretical analysis of a putative sex-role-reversed lek in the butterfly Acraea encedon. Sex-role reversal appears to have come about because of infection with a male-killing Wolbachia. The bacterium is highly prevalent in some populations, such that there is a dearth of males. Receptive females form dense aggregations, and it has been suggested that males preferentially select females uninfected with the bacterium. As with more conventional systems, this presents a theoretical problem exactly analogous to the lek paradox, namely what maintains female variation and hence why do males continue to choose? We model the evolution of a male choice gene that allows discrimination between infected and uninfected females, and show that the stable maintenance of both female variation and male choice is likely, so long as males make mistakes when discriminating between females. Furthermore, our model allows the maintenance, in a panmictic population, of a male killer that is perfectly transmitted. This is the first model to allow this result, and may explain the long-term persistence of a male killer in Hypolimnas bolina.

Comparative sequence analysis of the VHL tumor suppressor gene.

Comparative genome analysis may provide novel insights into gene evolution and function. To investigate the von Hippel-Lindau (VHL) disease tumor suppressor gene, we sequenced the VHL gene in seven primate species. Comparative analysis was performed for human, primate, and rodent VHL genes and for a putative Caenorhabditis elegans VHL homologue identified by database analysis. The VHL gene has two translation initiation sites (at codons 1 and 54); however, the relative importance of the full-length translation product (pVHL30) and that translated from the second internal translation initiation site (pVHL19) is unclear. The N-terminal sequence of pVHL30 contains eight copies of a GXEEX acidic repeat motif in human and higher primates, but only three copies were present in the marmoset, and only one copy was present in rodent VHL genes. Evolutionary analysis suggested that the N-terminal repetitive sequence in pVHL30 was of less functional importance than those regions present in both pVHL30 and pVHL19. The VHL gene product is reported to form complexes with various proteins including elongin B, elongin C, VBP-1, fibronectin, Spl, CUL2, and HIF-1. Although most of the regions in pVHL that had been implicated in binding specific proteins demonstrated evolutionary conservation, the carboxy-terminal putative VBP-1 binding site was less well conserved, suggesting that VBP-1 binding may have less functional significance. Although an amino acid substitution (K171T) close to the pVHL elongin binding region was found in baboon, analysis of the structure of human pVHL suggested that this substitution would not interfere with pVHL/elongin C interaction. In general, there was a good correlation between the pVHL domains that demonstrated most evolutionary conservation and those that were most frequently mutated in tumors. Analysis of human/C. elegans conservation and human germline and somatic mutation patterns identified a highly conserved mutation cluster region between codons 74 and 90. However, this region is likely to be important for the structural integrity of pVHL rather than representing an additional protein binding domain.

The evolutionary dynamics of male-killers and their hosts.

Male-killing bacteria are cytoplasmic sex-ratio distorters that are transmitted vertically through females of their insect hosts. The killing of male hosts by their bacteria is thought to be an adaptive bacterial trait because it augments the fitness of female hosts carrying clonal relatives of those bacteria. Here we attempt to explain observations of multiple male-killers in natural host populations. First we show that such male-killer polymorphism cannot be explained by a classical model of male-killing. We then show that more complicated models incorporating the evolution of resistance in hosts can explain male-killer polymorphism. However, this is only likely if resistance genes are very costly. We also consider the long-term evolutionary dynamics of male-killers, and show that evolution towards progressively more 'efficient' male-killers can be thwarted by the appearance of host resistance. The presence of a resistance gene can allow a less efficient male-killer to outcompete its rival and hence reverse the trend towards more efficient transmission and reduced metabolic load on the host.