Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Recombination, dominance and selection on amino acid polymorphism in the Drosophila genome: contrasting patterns on the X and fourth chromosomes

Surveys of nucleotide polymorphism and divergence indicate that the average selection coefficient on Drosophila proteins is weakly positive. Similar surveys in mitochondrial genomes and in the selfing plant Arabidopsis show that weak negative selection has operated. These differences have been attributed to the low recombination environment of mtDNA and Arabidopsis that has hindered adaptive evolution through the interference effects of linkage. We test this hypothesis with new sequence surveys of proteins lying in low recombination regions of the Drosophila genome. We surveyed >3800 bp across four proteins at the tip of the X chromosome and >3600 bp across four proteins on the fourth chromosome in 24 strains of D. melanogaster and 5 strains of D. simulans. This design seeks to study the interaction of selection and linkage by comparing silent and replacement variation in semihaploid (X chromosome) and diploid (fourth chromosome) environments lying in regions of low recombination. While the data do indicate very low rates of exchange, all four gametic phases were observed both at the tip of the X and across the fourth chromosome. Silent variation is very low at the tip of the X (thetaS = 0.0015) and on the fourth chromosome (thetaS = 0.0002), but the tip of the X shows a greater proportional loss of variation than the fourth shows relative to normal-recombination regions. In contrast, replacement polymorphism at the tip of the X is not reduced (thetaR = 0.00065, very close to the X chromosome average). MK and HKA tests both indicate a significant excess of amino acid polymorphism at the tip of the X relative to the fourth. Selection is significantly negative at the tip of the X (Nes = -1.53) and nonsignificantly positive on the fourth (Nes approximately 2.9), analogous to the difference between mtDNA (or Arabidopsis) and the Drosophila genome average. Our distal X data are distinct from regions of normal recombination where the X shows a deficiency of amino acid polymorphism relative to the autosomes, suggesting more efficient selection against recessive deleterious replacement mutations. We suggest that the excess amino acid polymorphism on the distal X relative to the fourth chromosome is due to (1) differences in the mutation rate for selected mutations on the distal X or (2) a greater relaxation of selection from stronger linkage-related interference effects on the distal X. This relaxation of selection is presumed to be greater in magnitude than the difference in efficiency of selection between X-linked vs. autosomal selection.

High mutation rate and predominance of insertions in the Caenorhabditis elegans nuclear genome

Mutations have pivotal functions in the onset of genetic diseases and are the fundamental substrate for evolution. However, present estimates of the spontaneous mutation rate and spectrum are derived from indirect and biased measurements. For instance, mutation rate estimates for Caenorhabditis elegans are extrapolated from observations on a few genetic loci with visible phenotypes and vary over an order of magnitude. Alternative approaches in mammals, relying on phylogenetic comparisons of pseudogene loci and fourfold degenerate codon positions, suffer from uncertainties in the actual number of generations separating the compared species and the inability to exclude biases associated with natural selection. Here we provide a direct and unbiased estimate of the nuclear mutation rate and its molecular spectrum with a set of C. elegans mutation-accumulation lines that reveal a mutation rate about tenfold higher than previous indirect estimates and an excess of insertions over deletions. Because deletions dominate patterns of C. elegans pseudogene variation, our observations indicate that natural selection might be significant in promoting small genome size, and challenge the prevalent assumption that pseudogene divergence accurately reflects the spontaneous mutation spectrum.

Recombination and base composition: the case of the highly self-fertilizing plant Arabidopsis thaliana

The effects of recombination and self-fertilization on base composition were investigated both theoretically and experimentally in the Arabidopsis genome. Levels of inbreeding modulate the effect of recombination on base composition.

Background

Rates of recombination can vary among genomic regions in eukaryotes, and this is believed to have major effects on their genome organization in terms of base composition, DNA repeat density, intron size, evolutionary rates and gene order. In highly self-fertilizing species such as Arabidopsis thaliana, however, heterozygosity is expected to be strongly reduced and recombination will be much less effective, so that its influence on genome organization should be greatly reduced.

Results

Here we investigated theoretically the joint effects of recombination and self-fertilization on base composition, and tested the predictions with genomic data from the complete A. thaliana genome. We show that, in this species, both codon-usage bias and GC content do not correlate with the local rates of crossing over, in agreement with our theoretical results.

Conclusions

We conclude that levels of inbreeding modulate the effect of recombination on base composition, and possibly other genomic features (for example, transposable element dynamics). We argue that inbreeding should be considered when interpreting patterns of molecular evolution.

Meiosis-specific yeast Hop1 protein promotes synapsis of double-stranded DNA helices via the formation of guanine quartets

In most eukaryotes, genetic exchange between paired homologs occurs in the context of a tripartite proteinaceous structure called the synaptonemal complex (SC). Genetic analyses have revealed that the genes encoding SC proteins are vital for meiotic chromosome pairing and recombination. However, the number, nature and/or the mechanism used by SC proteins to align chromosomes are yet to be clearly defined. Here, we show that Saccharomyces cerevisiae Hop1, a component of SC, was able to promote pairing of double-stranded DNA helices containing arrays of mismatched G/G sequences. Significantly, pairing was rapid and robust, independent of homology in the arms flanking the central G/G region, and required four contiguous guanine residues. Furthermore, data from truncated DNA double helices showed that 20 bp on either side of the 8 bp mismatched G/G region was essential for efficient synapsis. Methylation interference indicated that pairing between the two DNA double helices involves G quartets. These results suggest that Hop1 is likely to play a direct role in meiotic chromosome pairing and recombination by its ability to promote synapsis between double-stranded DNA helices containing arrays of G residues. To our knowledge, Hop1 is the first protein shown to promote synapsis of DNA double helices from yeast or any other organism.

Selection on codon usage in Drosophila americana

Synonymous codons are not used at random, significantly influencing the base composition of the genome. The selection-mutation-drift model proposes that this bias reflects natural selection in favor of a subset of preferred codons. Previous estimates in Drosophila of the intensity of selective forces involved seem too large to be reconciled with theoretical predictions of the level of codon bias. This probably results from confounding effects of the demographic histories of the species concerned. We have studied three species of the virilis group of Drosophila, which are more likely to satisfy the assumptions of the evolutionary models. We analyzed the patterns of polymorphism and divergence in a sample of 18 genes and applied a new method for estimating the intensity of selection on synonymous mutations based on the frequencies of unpreferred mutations among polymorphic sites. This yielded estimates of selection intensities (N(e)s) of the order of 0.65, which is more compatible with the observed levels of codon bias. Our results support the action of both selection and mutational bias on codon usage bias and suggest that codon usage and genome base composition in the D. americana lineage are in approximate equilibrium. Biased gene conversion may also contribute to the observed patterns.

Recombination drives the evolution of GC-content in the human genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue in the identification of functional sequence features. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We have analyzed the pattern of neutral substitutions in 14.3 Mb of primate noncoding regions. We show that the GC-content toward which sequences are evolving is strongly correlated (r(2) = 0.61, P </= 2 10(-16)) with the rate of crossovers (notably in females). This demonstrates that recombination drives the evolution of base composition in human (probably via the process of biased gene conversion). The present substitution patterns are very different from what they had been in the past, resulting in a major modification of the isochore structure of our genome. This non-equilibrium situation suggests that changes of recombination rates occur relatively frequently during evolution, possibly as a consequence of karyotype rearrangements. These results have important implications for understanding the spatial and temporal variations of substitution processes in a broad range of sexual organisms, and for detecting the hallmarks of natural selection in DNA sequences.

Molecular correlates of genes exhibiting RNAi phenotypes in Caenorhabditis elegans

Understanding genome-wide links between genotype and phenotype has generally been difficult due to both the complexity of phenotypes, and until recently, inaccessibility to large numbers of genes that might underlie a trait. To address this issue, we establish the association between particular RNAi phenotypes in Caenorhabditis elegans and sequence characteristics of the corresponding proteins and DNA. We find that genes showing RNAi phenotypes are long and highly expressed with little noncoding DNA and high rates of synonymous site substitution (KS). In addition, genes conferring RNAi phenotypes have significantly lower rates of nonsynonymous site substitution (KA). Collectively, these sequence features explain nearly 20% of the difference between the sets of loci that display or lack a RNAi-mediated effect, and reflect aspects both of the RNAi mechanism and the biological function of the genes. For example, the particularly low rate of evolution of genes in the sterility RNAi phenotype class suggests a role of C. elegans life history in shaping these patterns of sequence and expression characteristics on phenotypes. This approach also allows prediction of a set of heretofore-uncharacterized loci for which we expect future RNAi studies to reveal phenotypic effects (i.e., false negatives in present screens).

Factors affecting polymorphism at microsatellite loci in bread wheat [ Triticum aestivum (L.) Thell]: effects of mutation processes and physical distance from the centromere

The effects of factors known to influence the level of polymorphism at microsatellite loci were studied using 99 markers and seven lines of bread wheat. Mutational factors as well as indirect selective events shape diversity at these loci. Theory predicts that the selection of favorable alleles should reduce polymorphism at neutral neighboring loci in genomic areas with low recombination rates. In wheat, local recombination rate is positively correlated with physical distance from the centromere. Seventy four loci among the 99 used could be physically located on the chromosome. We studied how the following affected the diversity among a set of inbred lines: the length of the alleles, the motif (CA versus CT), the structure of the loci (perfect versus imperfect) and the chromosomal position of the loci. For each locus, we determined whether the polymorphism observed at a locus was compatible with the Stepwise Mutation Model (SMM) or the Two-Phase Model (TPM). Both the mutation rate and the compatibility with the SMM or the TPM were shown to be variable between loci. Wheat microsatellite loci were found to be more variable when segregating alleles were perfect and had long motifs (composed of many repetitions). Diversity observed at 19 loci was not compatible with the SMM. Loci located in distal regions, with presumably high recombination rates, had longer allele sizes and were more polymorphic than loci located in proximal regions. We conclude that both mutation factors and indirect selective events vary according to the local recombination rate and therefore jointly influence the level of polymorphism at microsatellite loci in wheat.

Recombination has little effect on the rate of sequence divergence in pseudoautosomal boundary 1 among humans and great apes

Recent studies indicated that recombination is strongly mutagenic. In particular, data from the mouse pseudoautosomal boundary (PAB) suggested that locally intensive recombination increased the nucleotide substitution rate by more than 100-fold and greatly increased the GC content. Here we study the rates of nucleotide substitution in eight introns of the human and great ape XG gene, which spans the boundary between the pseudoautosomal region 1 (PAR1) and the X-specific region. Contrary to what is expected under the above hypothesis, our sequence data from humans and great apes reveal that the PAR1 introns of XG have actually evolved slightly slower than X-specific introns. Only when a New World monkey was compared with hominoids were the rates slightly increased in the PAR1 introns. In terms of base composition, although the intergenic regions of the human PAR1 show a significant increase of G and C nucleotides, the base composition of the surveyed PAR1 introns is similar to that of the X-specific introns. Direct and indirect evidence indicates that the recombination rate is, indeed, much higher in PAR1 introns than in X-specific introns, and that the present PAB has persisted since the common ancestor of hominoids. Therefore, the mutagenic effect of recombination is far weaker than previously proposed, at least in hominoid PABs.

Intraspecific DNA variation in nuclear genes of the mosquito Aedes aegypti

Single nucleotide polymorphisms (SNPs) are an abundant source of genetic variation among individual organisms. To assess the usefulness of SNPs for genome analysis in the yellow fever mosquito, Aedes aegypti, we sequenced 25 nuclear genes in each of three strains and analysed nucleotide diversity. The average frequency of nucleotide variation was 12 SNPs per kilobase, indicating that nucleotide variation in Ae. aegypti is similar to that in other organisms, including Drosophila and the malaria vector Anopheles gambiae. Transition polymorphisms outnumbered transversion polymorphisms, at a ratio of about 2:1. We examined codon usage and confirmed that mutational bias favours G and C ending codons. Codon bias was most pronounced in highly expressed genes. Nucleotide diversity estimates indicated that substitution rates are positively correlated in coding and non-coding regions. Nucleotide diversity varied from one gene to another. The unequal distribution of SNPs among Ae. aegypti nuclear genes suggests that single base variations are non-neutral and are subject to selective constraints. Our analysis showed that ubiquitously expressed genes have lower polymorphism rates and are likely under strong purifying selection, whereas tissue specific genes and genes with a putative role in parasite defence exhibit higher levels of polymorphism that may be associated with diversifying selection.

Patterns of Selection Against Transposons Inferred From the Distribution of Tc1, Tc3 and Tc5 Insertions in the mut-7 Line of the Nematode Caenorhabditis elegans

To identify the factors (selective or mutational) that affect the distribution of transposable elements (TEs) within a genome, it is necessary to compare the pattern of newly arising element insertions to the pattern of element insertions that have been fixed in a population. To do this, we analyzed the distribution of recent mutant insertions of the Tc1, Tc3, and Tc5 elements in a mut-7 background of the nematode Caenorhabditis elegans and compared it to the distribution of element insertions (presumably fixed) within the sequenced genome. Tc1 elements preferentially insert in regions with high recombination rates, whereas Tc3 and Tc5 do not. Although Tc1 and Tc3 both insert in TA dinucleotides, there is no clear relationship between the frequency of insertions and the TA dinucleotide density. There is a strong selection against TE insertions within coding regions: the probability that a TE will be fixed is at least 31 times lower in coding regions than in noncoding regions. Contrary to the prediction of theoretical models, we found that the selective pressure against TE insertions does not increase with the recombination rate. These findings indicate that the distribution of these three transposon families in the genome of C. elegans is determined essentially by just two factors: the pattern of insertions, which is a characteristic of each family, and the selection against insertions within coding regions.

Linkage Disequilibrium Patterns Across a Recombination Gradient in African Drosophila melanogaster

Previous multilocus surveys of nucleotide polymorphism have documented a genome-wide excess of intralocus linkage disequilibrium (LD) in Drosophila melanogaster and D. simulans relative to expectations based on estimated mutation and recombination rates and observed levels of diversity. These studies examined patterns of variation from predominantly non-African populations that are thought to have recently expanded their ranges from central Africa. Here, we analyze polymorphism data from a Zimbabwean population of D. melanogaster, which is likely to be closer to the standard population model assumptions of a large population with constant size. Unlike previous studies, we find that levels of LD are roughly compatible with expectations based on estimated rates of crossing over. Further, a detailed examination of genes in different recombination environments suggests that markers near the telomere of the X chromosome show considerably less linkage disequilibrium than predicted by rates of crossing over, suggesting appreciable levels of exchange due to gene conversion. Assuming that these populations are near mutation-drift equilibrium, our results are most consistent with a model that posits heterogeneity in levels of exchange due to gene conversion across the X chromosome, with gene conversion being a minor determinant of LD levels in regions of high crossing over. Alternatively, if levels of exchange due to gene conversion are not negligible in regions of high crossing over, our results suggest a marked departure from mutation-drift equilibrium (i.e., toward an excess of LD) in this Zimbabwean population. Our results also have implications for the dynamics of weakly selected mutations in regions of reduced crossing over.

Translational selection and yeast proteome evolution

The primary structures of peptides may be adapted for efficient synthesis as well as proper function. Here, the Saccharomyces cerevisiae genome sequence, DNA microarray expression data, tRNA gene numbers, and functional categorizations of proteins are employed to determine whether the amino acid composition of peptides reflects natural selection to optimize the speed and accuracy of translation. Strong relationships between synonymous codon usage bias and estimates of transcript abundance suggest that DNA array data serve as adequate predictors of translation rates. Amino acid usage also shows striking relationships with expression levels. Stronger correlations between tRNA concentrations and amino acid abundances among highly expressed proteins than among less abundant proteins support adaptation of both tRNA abundances and amino acid usage to enhance the speed and accuracy of protein synthesis. Natural selection for efficient synthesis appears to also favor shorter proteins as a function of their expression levels. Comparisons restricted to proteins within functional classes are employed to control for differences in amino acid composition and protein size that reflect differences in the functional requirements of proteins expressed at different levels.

Selection conflicts, gene expression, and codon usage trends in yeast

Synonymous codon usage in yeast appears to be influenced by natural selection on gene expression, as well as regional variation in compositional bias. Because of the large number of potential targets of selection (i.e., most of the codons in the genome) and presumed small selection coefficients, codon usage is an excellent model for studying factors that limit the effectiveness of selection. We use factor analysis to identify major trends in codon usage for 5836 genes in Saccharomyces cerevisiae. The primary factor is strongly correlated with gene expression, consistent with the model that a subset of codons allows for more efficient translation. The secondary factor is very strongly correlated with third codon position GC content and probably reflects regional variation in compositional bias. We find that preferred codon usage decreases in the face of three potential limitations on the effectiveness of selection: reduced recombination rate, increased gene length, and reduced intergenic spacing. All three patterns are consistent with the Hill-Robertson effect (reduced effectiveness of selection among linked targets). A reduction in gene expression in closely spaced genes may also reflect selection conflicts due to antagonistic pleiotropy.

Codon bias and frequency-dependent selection on the hemagglutinin epitopes of influenza A virus

Although the surface proteins of human influenza A virus evolve rapidly and continually produce antigenic variants, the internal viral genes acquire mutations very gradually. In this paper, we analyze the sequence evolution of three influenza A genes over the past two decades. We study codon usage as a discriminating signature of gene- and even residue-specific diversifying and purifying selection. Nonrandom codon choice can increase or decrease the effective local substitution rate. We demonstrate that the codons of hemagglutinin, particularly those in the antibody-combining regions, are significantly biased toward substitutional point mutations relative to the codons of other influenza virus genes. We discuss the evolutionary interpretation and implications of these biases for hemagglutinin's antigenic evolution. We also introduce information-theoretic methods that use sequence data to detect regions of recent positive selection and potential protein conformational changes.

Biased gene conversion: implications for genome and sex evolution

Classical genetic studies show that gene conversion can favour some alleles over others. Molecular experiments suggest that gene conversion could favour GC over AT basepairs, leading to the concept of biased gene conversion towards GC (BGC(GC)). The expected consequence of such a process is the GC-enrichment of DNA sequences under gene conversion. Recent genomic work suggests that BGC(GC) affects the base composition of yeast, invertebrate and mammalian genomes. Hypotheses for the mechanisms and evolutionary origin of such a strange phenomenon have been proposed. Most BGC(GC) events probably occur during meiosis, which has implications for our understanding of the evolution of sex and recombination.

Evidence for co-evolution of gene order and recombination rate

There is increasing evidence in eukaryotic genomes that gene order is not random, even allowing for tandem duplication. Notably, in numerous genomes, genes of similar expression tend to be clustered. Are there other reasons for clustering of functionally similar genes? If genes are linked to enable genetic, rather than physical clustering, then we also expect that clusters of certain genes might be associated with blocks of reduced recombination rates. Here we show that, in yeast, essential genes are highly clustered and this clustering is independent of clustering of co-expressed genes and of tandem duplications. Adjacent pairs of essential genes are preferentially conserved through evolution. Notably, we also find that clusters of essential genes are in regions of low recombination and that larger clusters have lower recombination rates. These results suggest that selection acts to modify both the fine-scale intragenomic variation in the recombination rate and the distribution of genes and provide evidence for co-evolution of gene order and recombination rate.

Genome evolution: recombination speeds up adaptive evolution

Mutual interference among linked genetic sites subject to selection may reduce the level of adaptation. A recent study detected this effect using data on protein sequence evolution and codon usage in Drosophila.

Context dependence of meiotic recombination hotspots in yeast: the relationship between recombination activity of a reporter construct and base composition

Borde and colleagues reported that a reporter plasmid inserted at different genomic locations in Saccharomyces cerevisiae had different levels of meiotic recombination activity. We show that the level of recombination activity is very significantly correlated with the GC content of DNA sequences flanking the insertion.

Linkage limits the power of natural selection in Drosophila

Population genetic theory shows that the efficacy of natural selection is limited by linkage-selection at one site interferes with selection at linked sites. Such interference slows adaptation in asexual genomes and may explain the evolutionary advantage of sex. Here, we test for two signatures of constraint caused by linkage in a sexual genome, by using sequence data from 255 Drosophila melanogaster and Drosophila simulans loci. We find that (i) the rate of protein adaptation is reduced in regions of low recombination, and (ii) evolution at strongly selected amino acid sites interferes with optimal codon usage at weakly selected, tightly linked synonymous sites. Together these findings suggest that linkage limits the rate and degree of adaptation even in recombining genomes.

Hill-Robertson interference is a minor determinant of variations in codon bias across Drosophila melanogaster and Caenorhabditis elegans genomes

According to population genetics models, genomic regions with lower crossing-over rates are expected to experience less effective selection because of Hill-Robertson interference (HRi). The effect of genetic linkage is thought to be particularly important for a selection of weak intensity such as selection affecting codon usage. Consistent with this model, codon bias correlates positively with recombination rate in Drosophila melanogaster and Caenorhabditis elegans. However, in these species, the G+C content of both noncoding DNA and synonymous sites correlates positively with recombination, which suggests that mutation patterns and recombination are associated. To remove this effect of mutation patterns on codon bias, we used the synonymous sites of lowly expressed genes that are expected to be effectively neutral sites. We measured the differences between codon biases of highly expressed genes and their lowly expressed neighbors. In D. melanogaster we find that HRi weakly reduces selection on codon usage of genes located in regions of very low recombination; but these genes only comprise 4% of the total. In C. elegans we do not find any evidence for the effect of recombination on selection for codon bias. Computer simulations indicate that HRi poorly enhances codon bias if the local recombination rate is greater than the mutation rate. This prediction of the model is consistent with our data and with the current estimate of the mutation rate in D. melanogaster. The case of C. elegans, which is highly self-fertilizing, is discussed. Our results suggest that HRi is a minor determinant of variations in codon bias across the genome.

Rates and patterns of molecular evolution in inbred and outbred Arabidopsis

The evolution of self-fertilization is associated with a large reduction in the effective rate of recombination and a corresponding decline in effective population size. If many spontaneous mutations are slightly deleterious, this shift in the breeding system is expected to lead to a reduced efficacy of natural selection and genome-wide changes in the rates of molecular evolution. Here, we investigate the effects of the breeding system on molecular evolution in the highly self-fertilizing plant Arabidopsis thaliana by comparing its coding and noncoding genomic regions with those of its close outcrossing relative, the self-incompatible A. lyrata. More distantly related species in the Brassicaceae are used as outgroups to polarize the substitutions along each lineage. In contrast to expectations, no significant difference in the rates of protein evolution is observed between selfing and outcrossing Arabidopsis species. Similarly, no consistent overall difference in codon bias is observed between the species, although for low-biased genes A. lyrata shows significantly higher major codon usage. There is also evidence of intron size evolution in A. thaliana, which has consistently smaller introns than its outcrossing congener, potentially reflecting directional selection on intron size. The results are discussed in the context of heterogeneity in selection coefficients across loci and the effects of life history and population structure on rates of molecular evolution. Using estimates of substitution rates in coding regions and approximate estimates of divergence and generation times, the genomic deleterious mutation rate (U) for amino acid substitutions in Arabidopsis is estimated to be approximately 0.2-0.6 per generation.

Human SNP variability and mutation rate are higher in regions of high recombination

Understanding the co-variation of nucleotide diversity and local recombination rates is important both for the mapping of disease-associated loci and in understanding the causes of sequence evolution. It is known that single nucleotide polymorphisms (SNPs) around protein coding genes show higher diversity in regions of high recombination. Here, we find that this correlation holds for SNPs across the entire human genome, the great majority of which are not near exons or control elements. Contrasting with results from coding regions, we provide evidence that the higher nucleotide diversity in regions of high recombination is most likely due, at least in part, to a higher mutation rate. One possible explanation for this is that recombination is mutagenic.

Integrating genomics, bioinformatics, and classical genetics to study the effects of recombination on genome evolution

This study presents compelling evidence that recombination significantly increases the silent GC content of a genome in a selectively neutral manner, resulting in a highly significant positive correlation between recombination and "GC3s" in the yeast Saccharomyces cerevisiae. Neither selection nor mutation can explain this relationship. A highly significant GC-biased mismatch repair system is documented for the first time in any member of the Kingdom Fungi. Much of the variation in the GC3s within yeast appears to result from GC-biased gene conversion. Evidence suggests that GC-biased mismatch repair exists in numerous organisms spanning six kingdoms. This transkingdom GC mismatch repair bias may have evolved in response to a ubiquitous AT mutational bias. A significant positive correlation between recombination and GC content is found in many of these same organisms, suggesting that the processes influencing the evolution of the yeast genome may be a general phenomenon. Nonrecombining regions of the genome and nonrecombining genomes would not be subject to this type of molecular drive. It is suggested that the low GC content characteristic of many nonrecombining genomes may be the result of three processes (1) a prevailing AT mutational bias, (2) random fixation of the most common types of mutation, and (3) the absence of the GC-biased gene conversion which, in recombining organisms, permits the reversal of the most common types of mutation. A model is proposed to explain the observation that introns, intergenic regions, and pseudogenes typically have lower GC content than the silent sites of corresponding open reading frames. This model is based on the observation that the greater the heterology between two sequences, the less likely it is that recombination will occur between them. According to this "Constraint" hypothesis, the formation and propagation of heteroduplex DNA is expected to occur, on average, more frequently within conserved coding and regulatory regions of the genome. In organisms possessing GC-biased mismatch repair, this would enhance the GC content of these regions through biased gene conversion. These findings have a number of important implications for the way we view genome evolution and suggest a new model for the evolution of sex.

Recombination rate and the distribution of transposable elements in the Drosophila melanogaster genome

We analyzed the distribution of 54 families of transposable elements (TEs; transposons, LTR retrotransposons, and non-LTR retrotransposons) in the chromosomes of Drosophila melanogaster, using data from the sequenced genome. The density of LTR and non-LTR retrotransposons (RNA-based elements) was high in regions with low recombination rates, but there was no clear tendency to parallel the recombination rate. However, the density of transposons (DNA-based elements) was significantly negatively correlated with recombination rate. The accumulation of TEs in regions of reduced recombination rate is compatible with selection acting against TEs, as selection is expected to be weaker in regions with lower recombination. The differences in the relationship between recombination rate and TE density that exist between chromosome arms suggest that TE distribution depends on specific characteristics of the chromosomes (chromatin structure, distribution of other sequences), the TEs themselves (transposition mechanism), and the species (reproductive system, effective population size, etc.), that have differing influences on the effect of natural selection acting against the TE insertions.

Interactions between natural selection, recombination and gene density in the genes of Drosophila

In Drosophila, as in many organisms, natural selection leads to high levels of codon bias in genes that are highly expressed. Thus codon bias is an indicator of the intensity of one kind of selection that is experienced by genes and can be used to assess the impact of other genomic factors on natural selection. Among 13,000 genes in the Drosophila genome, codon bias has a slight positive, and strongly significant, association with recombination--as expected if recombination allows natural selection to act more efficiently when multiple linked sites segregate functional variation. The same reasoning leads to the expectation that the efficiency of selection, and thus average codon bias, should decline with gene density. However, this prediction is not confirmed. Levels of codon bias and gene expression are highest for those genes in an intermediate range of gene density, a pattern that may be the result of a tradeoff between the advantages for gene expression of close gene spacing and disadvantages arising from regulatory conflicts among tightly packed genes. These factors appear to overlay the more subtle effect of linkage among selected sites that gives rise to the association between recombination rate and codon bias.

Through a genome, darkly: comparative analysis of plant chromosomal DNA

Plant nuclear genomes encompass a wide range of variation in size and nucleotide composition with diverse arrangements of chromosomal segments, repetitive sequences and distribution of genes. Comparative genomic analysis may be undertaken at different levels of organisation, which are reflected in this review, together with a focus on the genetic and functional significance of the observed variation. Patterns of genome organisation have been revealed which reflect the different underlying mechanisms and constraints driving change. Thus comparative issues of genome size, nucleotide sequence composition and genome heterogeneity are provided as a background to understanding the different levels of segmental and repetitive sequence duplication and distribution of genes. The extent of synteny and collinearity revealed by recent genetic and sequence comparisons is discussed, together with a consideration of problems associated with such analyses. The possible origins and mechanisms of variation in genome size and organisation are covered, including the prevalence of duplication at different levels of organisation. The likely genetic, functional and adaptive consequences of replicated loci are discussed with evidence from comparative studies. The scope for comparative analysis of epigenetic plant genome variation is considered. Finally, opportunities for applying comparative genomics to isolating genes and understanding complex crop genomes are addressed.

The neutral theory in the genomic era

A number of tests have been developed to detect positive selection at the molecular level. These tests are based on DNA polymorphism within and divergence between species. Applications of these tests have revealed a large collection of genes that have evolved under positive selection and some general insights into adaptive evolution. Recently, these tests have been applied on a genomic scale and have provided estimates of the frequency of adaptive substitutions and a critical test of the neutral theory.

Gene expression and molecular evolution

The combination of complete genome sequence information and estimates of mRNA abundances have begun to reveal causes of both silent and protein sequence evolution. Translational selection appears to explain patterns of synonymous codon usage in many prokaryotes as well as a number of eukaryotic model organisms (with the notable exception of vertebrates). Relationships between gene length and codon usage bias, however, remain unexplained. Intriguing correlations between expression patterns and protein divergence suggest some general mechanisms underlying protein evolution.

Evaluation of methods for detecting recombination from DNA sequences: computer simulations

Recombination is a key evolutionary process that shapes the architecture of genomes and the genetic structure of populations. Although many statistical methods are available for the detection of recombination from DNA sequences, their absolute and relative performance is still unknown. Here we evaluated the performance of 14 different recombination detection algorithms. We used the coalescent with recombination to simulate DNA sequences with different levels of recombination, genetic diversity, and rate variation among sites. Recombination detection methods were applied to these data sets, and whether they detected or not recombination was recorded. Different recombination methods showed distinct performance depending on the amount of recombination, genetic diversity, and rate variation among sites. The model of nucleotide substitution under which the data were generated did not seem to have a significant effect. Most methods increase power with more sequence divergence. In general, recombination detection methods seem to capture the presence of recombination, but they are not very powerful. Methods that use substitution patterns or incompatibility among sites were more powerful than methods based on phylogenetic incongruence. Most methods do not seem to infer more false positives than expected by chance. Especially depending on the amount of diversity in the data, different methods could be used to attain maximum power while minimizing false positives. Results shown here will provide some guidance in the selection of the most appropriate method/s for the analysis of the particular data at hand.

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.

Selection for recombination in small populations

The reasons that sex and recombination are so widespread remain elusive. One popular hypothesis is that sex and recombination promote adaptation to a changing environment. The strongest evidence that increased recombination may evolve because recombination promotes adaptation comes from artificially selected populations. Recombination rates have been found to increase as a correlated response to selection on traits unrelated to recombination in several artificial selection experiments and in a comparison of domesticated and nondomesticated mammals. There are, however, several alternative explanations for the increase in recombination in such populations, including two different evolutionary explanations. The first is that the form of selection is epistatic, generating linkage disequilibria among selected loci, which can indirectly favor modifier alleles that increase recombination. The second is that random genetic drift in selected populations tends to generate disequilibria such that beneficial alleles are often found in different individuals; modifier alleles that increase recombination can bring together such favorable alleles and thus may be found in individuals with greater fitness. In this paper, we compare the evolutionary forces acting on recombination in finite populations subject to strong selection. To our surprise, we found that drift accounted for the majority of selection for increased recombination observed in simulations of small to moderately large populations, suggesting that, unless selected populations are large, epistasis plays a secondary role in the evolution of recombination.