Sex solves Haldane's dilemma.

The cumulative reproductive cost of multi-locus selection has been considered to be a potentially limiting factor on the rate of adaptive evolution. In this paper, we show that Haldane's arguments for the accumulation of reproductive costs over multiple loci are valid only for a clonally reproducing population of asexual genotypes. We show that a sexually reproducing population avoids this accumulation of costs. Thus, sex removes a perceived reproductive constraint on the rate of adaptive evolution. The significance of our results is twofold. First, the results demonstrate that adaptation based on multiple genes-such as selection acting on the standing genetic variation-does not entail a huge reproductive cost as suggested by Haldane, provided of course that the population is reproducing sexually. Second, this reduction in the cost of natural selection provides a simple biological explanation for the advantage of sex. Specifically, Haldane's calculations illustrate the evolutionary disadvantage of asexuality; sexual reproduction frees the population from this disadvantage.

The advantage of recombination when selection is acting at many genetic Loci.

Natural selection can act at many loci across the genome. But as the number of polymorphic loci increases linearly, the number of possible genotypic combinations increases exponentially. Consequently, a finite population - even a very large population - contains only a small sample of all possible multi-locus genotypes. In this paper, we revisit the classic Fisher-Muller models of recombination, taking into account the abundant standing variation that is commonly seen in natural populations. We show that the generation of new genotypic combinations through recombination is an important component of adaptive evolution based on multi-locus selection. Specifically, high-fitness genotypes are expected to be absent from the initial population when the frequencies of favorable alleles at the selected loci are low. But as the allele frequencies rise in response to selection the missing genotypes will be generated by recombination. Given recombination, if the average frequency of the favored alleles at the various selected loci is equal to p, then the expected number of favorable alleles per chromosome will be equal to pL, where L is the number of loci. As the value of p approaches unity at the selected loci, the number of favorable alleles per chromosome will approach a value of L, i.e., at the end of the selection process a favorable allele will be found at all loci. In the absence of recombination, however, selection will be limited to the highest-fitness genotypes that are already present in the initial population. We point out that the fitness of such initial genotypes is far less than the theoretical maximum fitness because they contain a favorable allele at only a fraction of the loci. Consequently, recombination acts to unblock the adaptive response to multi-locus selection in finite populations. Using simulations, we show that the sexual population can withstand invasion by newly-arising asexual clones. These results help explain the maintenance of sexual reproduction in natural populations.

The evolution of genomic GC content undergoes a rapid reversal within the genus Plasmodium.

The genome of the malarial parasite Plasmodium falciparum is extremely AT rich. This bias toward a low GC content is a characteristic of several, but not all, species within the genus Plasmodium. We compared 4283 orthologous pairs of protein-coding sequences between Plasmodium falciparum and the less AT-biased Plasmodium vivax. Our results indicate that the common ancestor of these two species was also extremely AT rich. This means that, although there was a strong bias toward A+T during the early evolution of the ancestral Plasmodium lineage, there was a subsequent reversal of this trend during the more recent evolution of some species, such as P. vivax. Moreover, we show that not only is the P. vivax genome losing its AT richness, it is actually gaining a very significant degree of GC richness. This example illustrates the potential volatility of nucleotide content during the course of molecular evolution. Such reversible fluxes in nucleotide content within lineages could have important implications for phylogenetic reconstruction based on molecular sequence data.

Resampling the pool of genotypic possibilities: an adaptive function of sexual reproduction.

BACKGROUND: Natural populations harbor significant levels of genetic variability. Because of this standing genetic variation, the number of possible genotypic combinations is many orders of magnitude greater than the population size. This means that any given population contains only a tiny fraction of all possible genotypic combinations. RESULTS: We show that recombination allows a finite population to resample the genotype pool, i.e., the universe of all possible genotypic combinations. Recombination, in combination with natural selection, enables an evolving sexual population to replace existing genotypes with new, higher-fitness genotypic combinations that did not previously exist in the population. This process allows the sexual population to gradually increase its fitness far beyond the range of fitnesses in the initial population. In contrast to this, an asexual population is limited to selection among existing lower fitness genotypes. CONCLUSIONS: The results provide an explanation for the ubiquity of sexual reproduction in evolving natural populations, especially when natural selection is acting on the standing genetic variation.

DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics.

DNA barcoding aims to provide an efficient method for species-level identifications and, as such, will contribute powerfully to taxonomic and biodiversity research. As the number of DNA barcode sequences accumulates, however, these data will also provide a unique 'horizontal' genomics perspective with broad implications. For example, here we compare the goals and methods of DNA barcoding with those of molecular phylogenetics and population genetics, and suggest that DNA barcoding can complement current research in these areas by providing background information that will be helpful in the selection of taxa for further analyses.

Finite populations, finite resources, and the evolutionary maintenance of genetic recombination.

Most previous models for the evolution of sex implicitly assume infinite population sizes and limitless resources. However, because favorable mutations are very rare and eukaryotic populations are finite, it has already been shown that multiple favorable mutants virtually never occur by chance. Therefore, sex is required to combine different favorable mutations into a single lineage. Second, we show that even when multiple favorable mutations do coexist, competition between genotypes can create negative epistasis for fitness, thus favoring recombination. Competition is especially effective when selection is at the level of viability in K-selected species living in a resource-limited environment. This means that recombination is advantageous both for incorporating new favorable mutations into the gene pool and for accelerating their increase to fixation. These advantages of recombination are diminished, however, as genome sizes decrease or as the amount of competition within the species is a less important component of selection.

A universal DNA mini-barcode for biodiversity analysis.

BACKGROUND: The goal of DNA barcoding is to develop a species-specific sequence library for all eukaryotes. A 650 bp fragment of the cytochrome c oxidase 1 (CO1) gene has been used successfully for species-level identification in several animal groups. It may be difficult in practice, however, to retrieve a 650 bp fragment from archival specimens, (because of DNA degradation) or from environmental samples (where universal primers are needed). RESULTS: We used a bioinformatics analysis using all CO1 barcode sequences from GenBank and calculated the probability of having species-specific barcodes for varied size fragments. This analysis established the potential of much smaller fragments, mini-barcodes, for identifying unknown specimens. We then developed a universal primer set for the amplification of mini-barcodes. We further successfully tested the utility of this primer set on a comprehensive set of taxa from all major eukaryotic groups as well as archival specimens. CONCLUSION: In this study we address the important issue of minimum amount of sequence information required for identifying species in DNA barcoding. We establish a novel approach based on a much shorter barcode sequence and demonstrate its effectiveness in archival specimens. This approach will significantly broaden the application of DNA barcoding in biodiversity studies.

Joint stationary moments of a two-island diffusion model of population subdivision.

An expression for joint stationary moments of a diffusion approximation to a generalized Wright-Fisher model, corresponding to two finite populations of equal sizes, with migration and mutation, is derived. This gives a complete description of the stationary distribution of allele frequencies in the balance between migration, mutation and genetic drift. We derive the sampling formula in terms of the joint stationary moments, and we also prove that the diffusion process corresponding to this model of population division is not reversible.

Rapid divergence of codon usage patterns within the rice genome.

BACKGROUND: Synonymous codon usage varies widely between genomes, and also between genes within genomes. Although there is now a large body of data on variations in codon usage, it is still not clear if the observed patterns reflect the effects of positive Darwinian selection acting at the level of translational efficiency or whether these patterns are due simply to the effects of mutational bias. In this study, we have included both intra-genomic and inter-genomic comparisons of codon usage. This allows us to distinguish more efficiently between the effects of nucleotide bias and translational selection. RESULTS: We show that there is an extreme degree of heterogeneity in codon usage patterns within the rice genome, and that this heterogeneity is highly correlated with differences in nucleotide content (particularly GC content) between the genes. In contrast to the situation observed within the rice genome, Arabidopsis genes show relatively little variation in both codon usage and nucleotide content. By exploiting a combination of intra-genomic and inter-genomic comparisons, we provide evidence that the differences in codon usage among the rice genes reflect a relatively rapid evolutionary increase in the GC content of some rice genes. We also noted that the degree of codon bias was negatively correlated with gene length. CONCLUSION: Our results show that mutational bias can cause a dramatic evolutionary divergence in codon usage patterns within a period of approximately two hundred million years. The heterogeneity of codon usage patterns within the rice genome can be explained by a balance between genome-wide mutational biases and negative selection against these biased mutations. The strength of the negative selection is proportional to the length of the coding sequences. Our results indicate that the large variations in synonymous codon usage are not related to selection acting on the translational efficiency of synonymous codons.

DNA asymmetric strand bias affects the amino acid composition of mitochondrial proteins.

Variations in GC content between genomes have been extensively documented. Genomes with comparable GC contents can, however, still differ in the apportionment of the G and C nucleotides between the two DNA strands. This asymmetric strand bias is known as GC skew. Here, we have investigated the impact of differences in nucleotide skew on the amino acid composition of the encoded proteins. We compared orthologous genes between animal mitochondrial genomes that show large differences in GC and AT skews. Specifically, we compared the mitochondrial genomes of mammals, which are characterized by a negative GC skew and a positive AT skew, to those of flatworms, which show the opposite skews for both GC and AT base pairs. We found that the mammalian proteins are highly enriched in amino acids encoded by CA-rich codons (as predicted by their negative GC and positive AT skews), whereas their flatworm orthologs were enriched in amino acids encoded by GT-rich codons (also as predicted from their skews). We found that these differences in mitochondrial strand asymmetry (measured as GC and AT skews) can have very large, predictable effects on the composition of the encoded proteins.

Evidence for strong selective constraint acting on the nucleotide composition of 16S ribosomal RNA genes.

Previous studies have shown that the guanine plus cytosine (G+C) content of ribosomal RNAs (rRNAs) is highly correlated with bacterial growth temperatures. This correlation is strongest in the double-stranded stem regions of the rRNA, a fact that can be explained by selection for increased structural stability at high growth temperatures. In this study, we examined the single-stranded regions of 16S rRNAs. We reasoned that, since these regions of the molecule are subject to less structural constraint than the stem regions, their nucleotide content might simply reflect the overall nucleotide content of the genome. Contrary to this expectation, however, we found that all of the single-stranded regions are characterized by very high adenine (A) and relatively low cytosine (C) contents. Moreover, the nucleotide content of these single-stranded regions is surprisingly constant between species, despite dramatic differences in optimal growth temperatures, and despite large differences in the overall genomic G+C content. This provides compelling evidence for strong stabilizing selection acting on 16S rRNA single-stranded regions. We found that selection favors purines (A+G), and especially adenine (A), in the single-stranded regions of these rRNAs.

Synonymous codon usage is subject to selection in thermophilic bacteria.

The patterns of synonymous codon usage, both within and among genomes, have been extensively studied over the past two decades. Despite the accumulating evidence that natural selection can shape codon usage, it has not been possible to link a particular pattern of codon usage to a specific external selective force. Here, we have analyzed the patterns of synonymous codon usage in 40 completely sequenced prokaryotic genomes. By combining the genes from several genomes (more than 80 000 genes in all) into a single dataset for this analysis, we were able to investigate variations in codon usage, both within and between genomes. The results show that synonymous codon usage is affected by two major factors: (i) the overall G+C content of the genome and (ii) growth at high temperature. This study focused on the relationship between synonymous codon usage and the ability to grow at high temperature. We have been able to eliminate both phylogenetic history and lateral gene transfer as possible explanations for the characteristic pattern of codon usage among the thermophiles. Thus, these results demonstrate a clear link between a particular pattern of codon usage and an external selective force.

Structural organization and chromosomal location of the chicken alpha-amylase gene family.

Characterization of the Gallus gallus alpha-amylase gene family revealed that the chicken genome contains two distinct amy loci. One of the two loci is expressed in the chicken pancreas while cDNA clones for the second locus were detected in a library constructed from liver mRNA. Fluorescent in situ hybridization to chromosome spreads showed that the two loci are both located on chromosome 8 within the chicken genome. Moreover, each locus contains both an intact, expressed gene copy as well as a pseudogene. The expressed gene and the pseudogene are arranged in a divergent configuration in the pancreatic amy locus, while in the hepatic locus the intact gene and the pseudogene are arranged in tandem. The data suggest a complex pattern of evolution for the chicken amylase gene family which includes multiple gene duplication events, insertion/deletion events, as well as changes in spatial expression patterns.

Thermophilic prokaryotes have characteristic patterns of codon usage, amino acid composition and nucleotide content.

A number of recent studies have shown that thermophilic prokaryotes have distinguishable patterns of both synonymous codon usage and amino acid composition, indicating the action of natural selection related to thermophily. On the other hand, several other studies of whole genomes have illustrated that nucleotide bias can have dramatic effects on synonymous codon usage and also on the amino acid composition of the encoded proteins. This raises the possibility that the thermophile-specific patterns observed at both the codon and protein levels are merely reflections of a single underlying effect at the level of nucleotide composition. Moreover, such an effect at the nucleotide level might be due entirely to mutational bias. In this study, we have compared the genomes of thermophiles and mesophiles at three levels: nucleotide content, codon usage and amino acid composition. Our results indicate that the genomes of thermophiles are distinguishable from mesophiles at all three levels and that the codon and amino acid frequency differences cannot be explained simply by the patterns of nucleotide composition. At the nucleotide level, we see a consistent tendency for the frequency of adenine to increase at all codon positions within the thermophiles. Thermophiles are also distinguished by their pattern of synonymous codon usage for several amino acids, particularly arginine and isoleucine. At the protein level, the most dramatic effect is a two-fold decrease in the frequency of glutamine residues among thermophiles. These results indicate that adaptation to growth at high temperature requires a coordinated set of evolutionary changes affecting (i) mRNA thermostability, (ii) stability of codon-anticodon interactions and (iii) increased thermostability of the protein products. We conclude that elevated growth temperature imposes selective constraints at all three molecular levels: nucleotide content, codon usage and amino acid composition. In addition to these multiple selective effects, however, the genomes of both thermophiles and mesophiles are often subject to superimposed large changes in composition due to mutational bias.

The DNA Barcode Linker.

DNA barcoding is based on the use of short DNA sequences to provide taxonomic tags for rapid, efficient identification of biological specimens. Currently, reference databases are being compiled. In the future, it will be important to facilitate access to these databases, especially for nonspecialist users. The method described here provides a rapid, web-based, user-friendly link between the DNA sequence from an unidentified biological specimen and various types of biological information, including the species name. Specifically, we use a customized, Google-type search algorithm to quickly match an unknown DNA sequence to a list of verified DNA barcodes in the reference database. In addition to retrieving the species name, our web tool also provides automatic links to a range of other information about that species. As the DNA barcode database becomes more populated, it will become increasingly important for the broader user community to be able to exploit it for the rapid identification of unknown specimens and to easily obtain relevant biological information about these species. The application presented here meets that need.

Highly similar noncoding genomic DNA sequences: ultraconserved, or merely widespread?

In this note, I propose an explanation for the seeming contradiction between bioinformatics-based predictions of an essential function for ultraconserved DNA sequences, and the lack of an experimental demonstration of such function.

An evolutionary footprint of age-related natural selection in mitochondrial DNA.

By comparing mtDNA sequences between different orders of mammals, we show that both longevity and generation time are significantly correlated with the nucleotide content of the mtDNA. Specifically, there is a positive correlation between generation time and mt GC content. This correlation is repeated, at a finer evolutionary scale, within the primates. Moreover, a comparison of human and chimpanzee mtDNAs shows that the effect has been very pronounced during the short evolutionary period since the divergence of these two species, with human mtDNA showing a GC-biased pattern of substitution at the variable sites. In addition to these DNA sequence patterns, comparisons between the human and the chimp mt protein sequences also revealed a surprisingly high substitution rate for threonine residues, resulting in a reduction of threonine in the human mt proteome. These patterns of both DNA and protein evolution can be explained by a balance between AT-biased mutational pressure and age-related purifying selection.

DNA barcodes provide a quick preview of mitochondrial genome composition.

DNA barcodes have achieved prominence as a tool for species-level identifications. Consequently, there is a rapidly growing database of these short sequences from a wide variety of taxa. In this study, we have analyzed the correlation between the nucleotide content of the short DNA barcode sequences and the genomes from which they are derived. Our results show that such short sequences can yield important, and surprisingly accurate, information about the composition of the entire genome. In other words, for unsequenced genomes, the DNA barcodes can provide a quick preview of the whole genome composition.

Assessing the effect of varying sequence length on DNA barcoding of fungi.

DNA barcoding shows enormous promise for the rapid identification of organisms at the species level. There has been much recent debate, however, about the need for longer barcode sequences, especially when these sequences are used to construct molecular phylogenies. Here, we have analysed a set of fungal mitochondrial sequences - of various lengths - and we have monitored the effect of reducing sequence length on the utility of the data for both species identification and phylogenetic reconstruction. Our results demonstrate that reducing sequence length has a profound effect on the accuracy of resulting phylogenetic trees, but surprisingly short sequences still yield accurate species identifications. We conclude that the standard short barcode sequences ( approximately 600 bp) are not suitable for inferring accurate phylogenetic relationships, but they are sufficient for species identification among the fungi.