The correlation between synonymous and nonsynonymous substitutions in Drosophila: mutation, selection or relaxed constraints?

Kimura's Theory of Non-Adaptive Radiation and Peto's Paradox: A Missing Link?

Karyotype diversity reflects genome integrity and stability. A strong correlation between karyotype diversity and species richness, meaning the number of species in a phylogenetic clade, was first reported in mammals over forty years ago: in mammalian phylogenetic clades, the standard deviation of karyotype diversity (KD) closely corresponded to species richness (SR) at the order level. These initial studies, however, did not control for phylogenetic signal, raising the possibility that the correlation was due to phylogenetic relatedness among species in a clade. Accordingly, karyotype diversity trivially reflects species richness simply as a passive consequence of adaptive radiation. A more recent study in mammals controlled for phylogenetic signals and established the correlation as phylogenetically independent, suggesting that species richness cannot, in itself, explain the observed corresponding karyotype diversity. The correlation is, therefore, remarkable because the molecular mechanisms contributing to karyotype diversity are evolutionarily independent of the ecological mechanisms contributing to species richness. Recently, it was shown in salamanders that the two processes generating genome size diversity and species richness were indeed independent and operate in parallel, suggesting a potential non-adaptive, non-causal but biologically meaningful relationship. KD depends on mutational input generating genetic diversity and reflects genome stability, whereas species richness depends on ecological factors and reflects natural selection acting on phenotypic diversity. As mutation and selection operate independently and involve separate and unrelated evolutionary mechanisms-there is no reason a priori to expect such a strong, let alone any, correlation between KD and SR. That such a correlation exists is more consistent with Kimura's theory of non-adaptive radiation than with ecologically based adaptive theories of macro-evolution, which are not excluded in Kimura's non-adaptive theory. The following reviews recent evidence in support of Kimura's proposal, and other findings that contribute to a wider understanding of the molecular mechanisms underlying the process of non-adaptive radiation.

Decoding the effects of synonymous variants

Synonymous single nucleotide variants (sSNVs) are common in the human genome but are often overlooked. However, sSNVs can have significant biological impact and may lead to disease. Existing computational methods for evaluating the effect of sSNVs suffer from the lack of gold-standard training/evaluation data and exhibit over-reliance on sequence conservation signals. We developed synVep (synonymous Variant effect predictor), a machine learning-based method that overcomes both of these limitations. Our training data was a combination of variants reported by gnomAD (observed) and those unreported, but possible in the human genome (generated). We used positive-unlabeled learning to purify the generated variant set of any likely unobservable variants. We then trained two sequential extreme gradient boosting models to identify subsets of the remaining variants putatively enriched and depleted in effect. Our method attained 90% precision/recall on a previously unseen set of variants. Furthermore, although synVep does not explicitly use conservation, its scores correlated with evolutionary distances between orthologs in cross-species variation analysis. synVep was also able to differentiate pathogenic vs. benign variants, as well as splice-site disrupting variants (SDV) vs. non-SDVs. Thus, synVep provides an important improvement in annotation of sSNVs, allowing users to focus on variants that most likely harbor effects.

A Genome-Wide Profiling of Glioma Patients with an IDH1 Mutation Using the Catalogue of Somatic Mutations in Cancer Database

Simple Summary

Glioma patients that present a somatic mutation in the isocitrate dehydrogenase 1 (IDH1) gene have a significantly better prognosis and overall survival than patients with the wild-type genotype. An IDH1 mutation is hypothesized to occur early during cellular transformation and leads to further genetic instability. A genome-wide profiling of glioma patients in the Catalogue of Somatic Mutations in Cancer (COSMIC) database was performed to classify the genetic differences in IDH1-mutant versus IDH1-wildtype patients. This classification will aid in a better understanding of how this specific mutation influences the genetic make-up of glioma and the resulting prognosis. Key differences in co-mutation and gene expression levels were identified that correlate with an improved prognosis.

Abstract

Gliomas are differentiated into two major disease subtypes, astrocytoma or oligodendroglioma, which are then characterized as either IDH (isocitrate dehydrogenase)-wild type or IDH-mutant due to the dramatic differences in prognosis and overall survival. Here, we investigated the genetic background of IDH1-mutant gliomas using the Catalogue of Somatic Mutations in Cancer (COSMIC) database. In astrocytoma patients, we found that IDH1 is often co-mutated with TP53, ATRX, AMBRA1, PREX1, and NOTCH1, but not CHEK2, EGFR, PTEN, or the zinc finger transcription factor ZNF429. The majority of the mutations observed in these genes were further confirmed to be either drivers or pathogenic by the Cancer-Related Analysis of Variants Toolkit (CRAVAT). Gene expression analysis showed down-regulation of DRG2 and MSN expression, both of which promote cell proliferation and invasion. There was also significant over-expression of genes such as NDRG3 and KCNB1 in IDH1-mutant astrocytoma patients. We conclude that IDH1-mutant glioma is characterized by significant genetic changes that could contribute to a better prognosis in glioma patients.

Analysis of selection in protein-coding sequences accounting for common biases

The evolution of protein-coding genes is usually driven by selective processes, which favor some evolutionary trajectories over others, optimizing the subsequent protein stability and activity. The analysis of selection in this type of genetic data is broadly performed with the metric nonsynonymous/synonymous substitution rate ratio (dN/dS). However, most of the well-established methodologies to estimate this metric make crucial assumptions, such as lack of recombination or invariable codon frequencies along genes, which can bias the estimation. Here, we review the most relevant biases in the dN/dS estimation and provide a detailed guide to estimate this metric using state-of-the-art procedures that account for such biases, along with illustrative practical examples and recommendations. We also discuss the traditional interpretation of the estimated dN/dS emphasizing the importance of considering complementary biological information such as the role of the observed substitutions on the stability and function of proteins. This review is oriented to help evolutionary biologists that aim to accurately estimate selection in protein-coding sequences.

Comparative Transcriptome Profiling of the Loaches Triplophysa bleekeri and Triplophysa rosa Reveals Potential Mechanisms of Eye Degeneration

Eye degeneration is one of the most obvious characteristics of organisms restricted to subterranean habitats. In cavefish, eye degeneration has evolved independently numerous times and each process is associated with different genetic mechanisms. To gain a better understanding of these mechanisms, we compared the eyes of adult individuals of the cave loach Triplophysa rosa and surface loach Triplophysa bleekeri. Compared with the normal eyes of the surface loach, those of the cave loach were found to possess a small abnormal lens and a defective retina containing photoreceptor cells that lack outer segments. Sequencing of the transcriptomes of both species to identify differentially expressed genes (DEGs) and genes under positive selection revealed 4,802 DEGs and 50 genes under positive selection (dN/dS > 1, FDR < 0.1). For cave loaches, we identified one Gene Ontology category related to vision that was significantly enriched in downregulated genes. Specifically, we found that many of the downregulated genes, including pitx3, lim2, crx, gnat2, rx1, rho, prph2, and β|γ-crystallin are associated with lens/retinal development and maintenance. However, compared with those in the surface loach, the lower dS rates but higher dN rates of the protein-coding sequences in T. rosa indicate that changes in amino acid sequences might be involved in the adaptation and visual degeneration of cave loaches. We also found that genes associated with light perception and light-stimulated vision have evolved at higher rates (some genes dN/dS > 1 but FDR > 0.1). Collectively, the findings of this study indicate that the degradation of cavefish vision is probably associated with both gene expression and amino acid changes and provide new insights into the mechanisms underlying the degeneration of cavefish eyes.

Morphology and genome of a snailfish from the Mariana Trench provide insights into deep-sea adaptation

It is largely unknown how living organisms-especially vertebrates-survive and thrive in the coldness, darkness and high pressures of the hadal zone. Here, we describe the unique morphology and genome of Pseudoliparis swirei-a recently described snailfish species living below a depth of 6,000 m in the Mariana Trench. Unlike closely related shallow sea species, P. swirei has transparent, unpigmented skin and scales, thin and incompletely ossified bones, an inflated stomach and a non-closed skull. Phylogenetic analyses show that P. swirei diverged from a close relative living near the sea surface about 20 million years ago and has abundant genetic diversity. Genomic analyses reveal that: (1) the bone Gla protein (bglap) gene has a frameshift mutation that may cause early termination of cartilage calcification; (2) cell membrane fluidity and transport protein activity in P. swirei may have been enhanced by changes in protein sequences and gene expansion; and (3) the stability of its proteins may have been increased by critical mutations in the trimethylamine N-oxide-synthesizing enzyme and hsp90 chaperone protein. Our results provide insights into the morphological, physiological and molecular evolution of hadal vertebrates.

Evolutionary patterns of Toll-like receptor signaling pathway genes in the Suidae

BACKGROUND

The Toll-like receptor (TLR) signaling pathway constitutes an essential component of the innate immune system. Highly conserved proteins, indicative of their critical roles in host survival, characterize this pathway. Selective constraints could vary depending on the gene's position within the pathway as TLR signaling is a sequential process and that genes downstream of the TLRs may be more selectively constrained to ensure efficient immune responses given the important role of downstream genes in the signaling process. Thus, we investigated whether gene position influenced protein evolution in the TLR signaling pathway of the Suidae. The members of the Suidae examined included the European Sus scrofa (wild boar), Asian Sus scrofa (wild boar), Sus verrucosus, Sus celebensis, Sus scebifrons, Sus barbatus, Babyrousa babyrussa, Potamochoerus larvatus, Potamochoerus porcus and Phacochoerus africanus.

RESULTS

A total of 33 TLR signaling pathway genes in the Suidae were retrieved from resequencing data. The evolutionary parameter ω (dn/ds) had an overall mean of 0.1668 across genes, indicating high functional conservation within the TLR signaling pathway. A significant relationship was inferred for the network parameters gene position, number of protein-protein interactions, protein length and the evolutionary parameter dn (nonsynonymous substitutions) such that downstream genes had lower nonsynonymous substitution rates, more interactors and shorter protein length than upstream genes. Gene position was significantly correlated with the number of protein-protein interactions and protein length. Thus, the polarity in the selective constraint along the TLR signaling pathway was due to the number of molecules a protein interacted with and the protein's length.

CONCLUSION

Results indicate that the level of selective constraints on genes within the TLR signaling pathway of the Suidae is dependent on the gene's position and network parameters. In particular, downstream genes evolve more slowly as a result of being highly connected and having shorter protein lengths. These findings highlight the critical role of gene network parameters in gene evolution.

Evolutionary Rate Heterogeneity of Primary and Secondary Metabolic Pathway Genes in Arabidopsis thaliana

Primary metabolism is essential to plants for growth and development, and secondary metabolism helps plants to interact with the environment. Many plant metabolites are industrially important. These metabolites are produced by plants through complex metabolic pathways. Lack of knowledge about these pathways is hindering the successful breeding practices for these metabolites. For a better knowledge of the metabolism in plants as a whole, evolutionary rate variation of primary and secondary metabolic pathway genes is a prerequisite. In this study, evolutionary rate variation of primary and secondary metabolic pathway genes has been analyzed in the model plant Arabidopsis thaliana. Primary metabolic pathway genes were found to be more conserved than secondary metabolic pathway genes. Several factors such as gene structure, expression level, tissue specificity, multifunctionality, and domain number are the key factors behind this evolutionary rate variation. This study will help to better understand the evolutionary dynamics of plant metabolism.

Codon and Amino Acid Usage Are Shaped by Selection Across Divergent Model Organisms of the Pancrustacea

In protein-coding genes, synonymous codon usage and amino acid composition correlate to expression in some eukaryotes, and may result from translational selection. Here, we studied large-scale RNA-seq data from three divergent arthropod models, including cricket (Gryllus bimaculatus), milkweed bug (Oncopeltus fasciatus), and the amphipod crustacean Parhyale hawaiensis, and tested for optimization of codon and amino acid usage relative to expression level. We report strong signals of AT3 optimal codons (those favored in highly expressed genes) in G. bimaculatus and O. fasciatus, whereas weaker signs of GC3 optimal codons were found in P. hawaiensis, suggesting selection on codon usage in all three organisms. Further, in G. bimaculatus and O. fasciatus, high expression was associated with lowered frequency of amino acids with large size/complexity (S/C) scores in favor of those with intermediate S/C values; thus, selection may favor smaller amino acids while retaining those of moderate size for protein stability or conformation. In P. hawaiensis, highly transcribed genes had elevated frequency of amino acids with large and small S/C scores, suggesting a complex dynamic in this crustacean. In all species, the highly transcribed genes appeared to favor short proteins, high optimal codon usage, specific amino acids, and were preferentially involved in cell-cycling and protein synthesis. Together, based on examination of 1,680,067, 1,667,783, and 1,326,896 codon sites in G. bimaculatus, O. fasciatus, and P. hawaiensis, respectively, we conclude that translational selection shapes codon and amino acid usage in these three Pancrustacean arthropods.

Effective population size does not predict codon usage bias in mammals

Synonymous codons are not used at equal frequency throughout the genome, a phenomenon termed codon usage bias (CUB). It is often assumed that interspecific variation in the intensity of CUB is related to species differences in effective population sizes (N_e), with selection on CUB operating less efficiently in species with small N_e. Here, we specifically ask whether variation in N_e predicts differences in CUB in mammals and report two main findings. First, across 41 mammalian genomes, CUB was not correlated with two indirect proxies of N_e (body mass and generation time), even though there was statistically significant evidence of selection shaping CUB across all species. Interestingly, autosomal genes showed higher codon usage bias compared to X-linked genes, and high-recombination genes showed higher codon usage bias compared to low recombination genes, suggesting intraspecific variation in N_e predicts variation in CUB. Second, across six mammalian species with genetic estimates of N_e (human, chimpanzee, rabbit, and three mouse species: Mus musculus, M. domesticus, and M. castaneus), N_e and CUB were weakly and inconsistently correlated. At least in mammals, interspecific divergence in N_e does not strongly predict variation in CUB. One hypothesis is that each species responds to a unique distribution of selection coefficients, confounding any straightforward link between N_e and CUB.

Mitochondrial-nuclear interactions: compensatory evolution or variable functional constraint among vertebrate oxidative phosphorylation genes?

Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate dS and nonsynonymous substitution rate dN) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d(N) of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d(N) between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

POSITIVE DARWINIAN SELECTION: DOES THE COMPARATIVE METHOD RULE?

To detect positive Darwinian selection it is thought essential to compare two sequences. Despite its defects, "the comparative method rules." However, genes evolving rapidly under positive selection conflict more with internal forces (the genome phenotype) than genes evolving slowly under negative selection. In particular, there is conflict with stem-loop potential. The conflict between protein-encoding potential (primary information) and stem-loop potential (secondary information) permits detection of positive selection in a single sequence. The degree to which secondary information is compromised provides a measure of the speed of transmission of primary information. Thus, the sovereignty of the comparative method is challenged not only by its own defects, but also by the availability of a single-sequence method. However, while of limited utility for positive selection, the comparative method casts new light on Darwin's great question — the origin of species. Comparison of rates of synonymous and non-synonymous mutation suggests that branching into new species begins with synonymous mutations.

Genetic architecture of metabolic rate: environment specific epistasis between mitochondrial and nuclear genes in an insect

The extent to which mitochondrial DNA (mtDNA) variation is involved in adaptive evolutionary change is currently being reevaluated. In particular, emerging evidence suggests that mtDNA genes coevolve with the nuclear genes with which they interact to form the energy producing enzyme complexes in the mitochondria. This suggests that intergenomic epistasis between mitochondrial and nuclear genes may affect whole-organism metabolic phenotypes. Here, we use crossed combinations of mitochondrial and nuclear lineages of the seed beetle Callosobruchus maculatus and assay metabolic rate under two different temperature regimes. Metabolic rate was affected by an interaction between the mitochondrial and nuclear lineages and the temperature regime. Sequence data suggests that mitochondrial genetic variation has a role in determining the outcome of this interaction. Our genetic dissection of metabolic rate reveals a high level of complexity, encompassing genetic interactions over two genomes, and genotype × genotype × environment interactions. The evolutionary implications of these results are twofold. First, because metabolic rate is at the root of life histories, our results provide insights into the complexity of life-history evolution in general, and thermal adaptation in particular. Second, our results suggest a mechanism that could contribute to the maintenance of nonneutral mtDNA polymorphism.

Similar selective factors affect both between-gene and between-exon divergence in Drosophila

As a consequence of alternative splicing, a gene's exons will have different frequencies of inclusion into mature mRNA and different patterns of expression. These differences affect their patterns of evolutionary divergence. Using the recently reannotated genome of Drosophila melanogaster and the genome sequences of four closely related species of the melanogaster subgroup, we investigated the effect of alternative splicing, inclusion level (defined as the number of transcripts an exon is found in), and expression pattern on exon evolution across divergence times ranging from 1 to 12.5 Ma. Genes undergoing alternative splicing have a broader pattern of expression associated with a lower divergence rate in comparison with genes with a single annotated protein isoform. Within genes undergoing alternative splicing, we report a significant effect of inclusion level on exon evolution, as alternatively spliced exons are less conserved than constitutively spliced exons. More generally, there are significant negative correlations between inclusion level and exon evolutionary rates that can be associated with relaxation of selection. A significant effect of expression pattern on evolution rates is also observed. Overall, we found that similar selective factors such as the expression level and the pattern of expression are affecting both gene and exon evolution.

Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.

Population genomics, the study of genome-wide patterns of sequence variation within and between closely related species, can provide a comprehensive view of the relative importance of mutation, recombination, natural selection, and genetic drift in evolution. It can also provide fundamental insights into the biological attributes of organisms that are specifically shaped by adaptive evolution. One approach for generating population genomic datasets is to align DNA sequences from whole-genome shotgun projects to a standard reference sequence. We used this approach to carry out whole-genome analysis of polymorphism and divergence in Drosophila simulans, a close relative of the model system, D. melanogaster. We find that polymorphism and divergence fluctuate on a large scale across the genome and that these fluctuations are probably explained by natural selection rather than by variation in mutation rates. Our analysis suggests that adaptive protein evolution is common and is often related to biological processes that may be associated with gene expression, chromosome biology, and reproduction. The approaches presented here will have broad applicability to future analysis of population genomic variation in other systems, including humans.

Low-coverage genome sequences from multiple Drosophila simulans strains provide the first comprehensive view of polymorphism and divergence in the fruit fly.

Mutation rate variation in multicellular eukaryotes: causes and consequences

A basic knowledge about mutation rates is central to our understanding of a myriad of evolutionary phenomena, including the maintenance of sex and rates of molecular evolution. Although there is substantial evidence that mutation rates vary among taxa, relatively little is known about the factors that underlie this variation at an empirical level, particularly in multicellular eukaryotes. Here we integrate several disparate lines of theoretical and empirical inquiry into a unified framework to guide future studies that are aimed at understanding why and how mutation rates evolve in multicellular species.

Weak selection and recent mutational changes influence polymorphic synonymous mutations in humans

Recent large-scale genomic and evolutionary studies have revealed the small but detectable signature of weak selection on synonymous mutations during mammalian evolution, likely acting at the level of translational efficacy (i.e., translational selection). To investigate whether weak selection, and translational selection in particular, plays any role in shaping the fate of synonymous mutations that are present today in human populations, we studied genetic variation at the polymorphic level and patterns of evolution in the human lineage after human-chimpanzee separation. We find evidence that neutral mechanisms are influencing the frequency of polymorphic mutations in humans. Our results suggest a recent increase in mutational tendencies toward AT, observed in all isochores, that is responsible for AT mutations segregating at lower frequencies than GC mutations. In all, however, changes in mutational tendencies and other neutral scenarios are not sufficient to explain a difference between synonymous and noncoding mutations or a difference between synonymous mutations potentially advantageous or deleterious under a translational selection model. Furthermore, several estimates of selection intensity on synonymous mutations all suggest a detectable influence of weak selection acting at the level of translational selection. Thus, random genetic drift, recent changes in mutational tendencies, and weak selection influence the fate of synonymous mutations that are present today as polymorphisms. All of these features, neutral and selective, should be taken into account in evolutionary analyses that often assume constancy of mutational tendencies and complete neutrality of synonymous mutations.

Magic bullets and golden rules: data sampling in molecular phylogenetics

Data collection for molecular phylogenetic studies is based on samples of both genes and taxa. In an ideal world, with no limitations to resources, as many genes could be sampled as deemed necessary to address phylogenetic problems. Given limited resources in the real world, inadequate (in terms of choice of genes or number of genes) sequences or restricted taxon sampling can adversely affect the reliability or information gained in phylogenetics. Recent empirical and simulation-based studies of data sampling in molecular phylogenetics have reached differing conclusions on how to deal with these problems. Some advocated sampling more genes, others more taxa. There is certainly no 'magic bullet' that will fit all phylogenetic problems, and no specific 'golden rules' have been deduced, other than that single genes may not always contain sufficient phylogenetic information. However, several general conclusions and suggestions can be made. One suggestion is that the determination of a multiple, but moderate number (e.g., 6-10) of gene sequences might take precedence over sequencing a larger set of genes and thereby permit the sampling of more taxa for a phylogenetic study.

The soluble epoxide hydrolase gene harbors sequence variation associated with susceptibility to and protection from incident ischemic stroke

Stroke is the leading cause of severe disability and the third leading cause of death, accounting for one of every 15 deaths in the USA. We investigated the association of polymorphisms in the soluble epoxide hydrolase gene (EPHX2) with incident ischemic stroke in African-Americans and Whites. Twelve single nucleotide polymorphisms (SNPs) spanning EPHX2 were genotyped in a case-cohort sample of 1336 participants from the Atherosclerosis Risk in Communities (ARIC) study. In each racial group, Cox proportional hazard models were constructed to assess the relationship between incident ischemic stroke and EPHX2 polymorphisms. A score test method was used to investigate the association of common haplotypes of the gene with risk of ischemic stroke. In African-Americans, two common EPHX2 haplotypes with significant and opposing relationships to ischemic stroke risk were identified. In Whites, two common haplotypes showed suggestive indication of an association with ischemic stroke risk but, as in African-Americans, these relationships were in opposite direction. These findings suggest that multiple variants exist within or near the EPHX2 gene, with greatly contrasting relationships to ischemic stroke incidence; some associated with a higher incidence and others with a lower incidence.

Determination of mutation trend in proteins by means of translation probability between RNA codes and mutated amino acids

In this study, we estimate the translation probability to amino acid from RNA codon. With the determined 183 translation probabilities and amino-acid composition of eight highly mutated proteins, we construct the theoretical distributions of mutated amino acids in these proteins and then compare them with their actual distributions affected by mutations. Thereafter we trace the pattern of translation probabilities from RNA codons to mutated amino acids of 1053 point missense mutations. Finally, we statistically conclude that the natural mutation trend goes along the theoretical translation probability.

Ancient divergence in bathypelagic lake tanganyika deepwater cichlids: mitochondrial phylogeny of the tribe bathybatini

The cichlid species flock of Lake Tanganyika represents a polyphyletic assemblage of eight ancestral lineages, which colonized the emerging lake independently. Our study is focused on one of these lineages, the Bathybatini, a tribe of specialized piscivorous cichlids of the deep pelagic zone. By analyzing three mtDNA gene segments of all eight species of the tribe and two species of the closely related Trematocarini, we propose on the basis of a linearized tree analysis that the Bathybatini comprise two distinct lineages, the genera Hemibates and Bathybates, that seeded the primary lacustrine Tanganyika radiation independently. The genus Hemibates is likely to represent a distinct lineage that emerged simultaneously with the tribe Trematocarini and the genus Bathybates and should be therefore treated as a distinct tribe. Within the genus Bathybates, B. minor clearly represents the most ancestral split and is likely to have diverged from the remaining species in the course of the "primary lacustrine Tanganyika radiation" during which also the radiations of the Lamprologini and the H-lineage took place. The remaining "large" Bathybates species also diversified almost simultaneously and in step with the diversification of other Tanganyikan lineages-the Limnochromini and Cyprichromini-with B. graueri occupying the most ancestral branch, suggesting that these were induced by the same environmental changes. The lack of geographic color morphs suggests that competition and resource partitioning, rather than allopatric speciation, promoted speciation within the genus Bathybates.

Codon Bias and Noncoding GC Content Correlate Negatively with Recombination Rate on the Drosophila X Chromosome

The patterns and processes of molecular evolution may differ between the X chromosome and the autosomes in Drosophila melanogaster. This may in part be due to differences in the effective population size between the two chromosome sets and in part to the hemizygosity of the X chromosome in Drosophila males. These and other factors may lead to differences both in the gene complements of the X and the autosomes and in the properties of the genes residing on those chromosomes. Here we show that codon bias and recombination rate are correlated strongly and negatively on the X chromosome, and that this correlation cannot be explained by indirect relationships with other known determinants of codon bias. This is in dramatic contrast to the weak positive correlation found on the autosomes. We explored possible explanations for these patterns, which required a comprehensive analysis of the relationships among multiple genetic properties such as protein length and expression level. This analysis highlights conserved features of coding sequence evolution on the X and the autosomes and illuminates interesting differences between these two chromosome sets.

Patterns of synonymous codon usage in Drosophila melanogaster genes with sex-biased expression

The nonrandom use of synonymous codons (codon bias) is a well-established phenomenon in Drosophila. Recent reports suggest that levels of codon bias differ among genes that are differentially expressed between the sexes, with male-expressed genes showing less codon bias than female-expressed genes. To examine the relationship between sex-biased gene expression and level of codon bias on a genomic scale, we surveyed synonymous codon usage in 7276 D. melanogaster genes that were classified as male-, female-, or non-sex-biased in their expression in microarray experiments. We found that male-biased genes have significantly less codon bias than both female- and non-sex-biased genes. This pattern holds for both germline and somatically expressed genes. Furthermore, we find a significantly negative correlation between level of codon bias and degree of sex-biased expression for male-biased genes. In contrast, female-biased genes do not differ from non-sex-biased genes in their level of codon bias and show a significantly positive correlation between codon bias and degree of sex-biased expression. These observations cannot be explained by differences in chromosomal distribution, mutational processes, recombinational environment, gene length, or absolute expression level among genes of the different expression classes. We propose that the observed codon bias differences result from differences in selection at synonymous and/or linked nonsynonymous sites between genes with male- and female-biased expression.

Correlated evolution of synonymous and nonsynonymous sites in Drosophila

Recent work has shown that Drosophila melanogaster genes with fast-evolving nonsynonymous sites have lower codon usage bias. This pattern has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. Here we have looked for this correlation in a much larger and less biased dataset, comprising 630 gene pairs from D. melanogaster and D. yakuba. We confirmed that there is a negative correlation between the rate of nonsynonymous substitutions (d(N)) and codon bias in D. melanogaster. We then tested the interference hypothesis and other alternative explanations, including one involving gene expression. We found that d(N) indeed correlates with the level of gene expression. Given that gene expression is a strong determinant of codon bias, the relationship between d(N) and codon bias might be a by-product of gene expression. However, our tests show that none of the hypotheses we consider seem to explain the data fully.

An evolutionary analysis of orphan genes in Drosophila

Orphan genes are protein-coding regions that have no recognizable homolog in distantly related species. A substantial fraction of coding regions in any genome sequenced consists of orphan genes, but the evolutionary and functional significance of orphan genes is not understood. We present a reanalysis of the Drosophila melanogaster proteome that shows that there are still between 26% and 29% of all proteins without a significant match with noninsect sequences, and that these orphans are underrepresented in genetic screens. To analyze the characteristics of orphan genes in Drosophila, we used sequence comparisons between cDNAs retrieved from two Drosophila yakuba libraries and their corresponding D. melanogaster orthologs. We find that a cDNA library from adults yields twice as many orphan genes as such a library from embryos. The orphan genes evolve on average more than three times faster than nonorphan genes, although the width of the evolutionary rate distribution is similar for the two classes. In particular, some orphan genes show very low substitution rates that are comparable to otherwise highly conserved genes. We propose a model suggesting that orphans may be involved in the evolution of adaptive traits, and that slow-evolving orphan genes may be particularly interesting candidate genes for identifying lineage-specific adaptations.

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.

Evolutionary patterns of codon usage in the chloroplast gene rbcL

In this study we reconstruct the evolution of codon usage bias in the chloroplast gene rbcL using a phylogeny of 92 green-plant taxa. We employ a measure of codon usage bias that accounts for chloroplast genomic nucleotide content, as an attempt to limit plausible explanations for patterns of codon bias evolution to selection- or drift-based processes. This measure uses maximum likelihood-ratio tests to compare the performance of two models, one in which a single codon is overrepresented and one in which two codons are overrepresented. The measure allowed us to analyze both the extent of bias in each lineage and the evolution of codon choice across the phylogeny. Despite predictions based primarily on the low G + C content of the chloroplast and the high functional importance of rbcL, we found large differences in the extent of bias, suggesting differential molecular selection that is clade specific. The seed plants and simple leafy liverworts each independently derived a low level of bias in rbcL, perhaps indicating relaxed selectional constraint on molecular changes in the gene. Overrepresentation of a single codon was typically plesiomorphic, and transitions to overrepresentation of two codons occurred commonly across the phylogeny, possibly indicating biochemical selection. The total codon bias in each taxon, when regressed against the total bias of each amino acid, suggested that twofold amino acids play a strong role in inflating the level of codon usage bias in rbcL, despite the fact that twofolds compose a minority of residues in this gene. Those amino acids that contributed most to the total codon usage bias of each taxon are known through amino acid knockout and replacement to be of high functional importance. This suggests that codon usage bias may be constrained by particular amino acids and, thus, may serve as a good predictor of what residues are most important for protein fitness.

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.

Evolutionary dynamics of mammalian mRNA untranslated regions by comparative analysis of orthologous human, artiodactyl and rodent gene pairs

Most evolutionary studies based on molecular data refer to the portion of genomes encoding for proteins. Today, however, more and more attention is paid to the so-called 'non-coding' regions, which constitute a notable portion of the metazoan nuclear genome. Among them, the untranslated regions of messenger RNAs (mRNA UTRs) are particularly important, as they are involved in the regulation of gene expression, controlling translation efficiency as well as mRNA localization and stability. Up to now, only few studies have focused on the analysis of the compositional and structural features of UTRs, or carried out to investigate quantitatively their evolutionary dynamics. For this reason we have carried out an inter-order study on the evolutionary rate of 5' and 3' UTRs with respect to the corresponding coding region in 93 triplets of orthologous genes (selected through a phylogenetic approach, for a total of 645 625 nt) belonging to Primates (Homo sapiens), Artiodactyla (Bos taurus) and Rodentia (Mus spp.). Our study, that considered only likely orthologous genes, has revealed interesting features on the evolution of these regions concerning nucleotide substitution rate and indels and repetitive element distribution. UTRs from different genes showed a remarkable heterogeneity in the evolutionary dynamics, with some homologous so highly divergent to prevent their alignment, and other rather conserved, at least in some regions, most divergent sequence pairs were excluded from our analysis. The comparison between the nucleotide substitution rates calculated for 5' and 3' UTRs with those calculated on synonymous coding position allowed us to verify and measure the existence of functional constraints acting upon the UTRs of different genes which have shown, in many cases, a positive selection driven evolutionary dynamics.

Population, evolutionary and genomic consequences of interference selection

Weakly selected mutations are most likely to be physically clustered across genomes and, when sufficiently linked, they alter each others' fixation probability, a process we call interference selection (IS). Here we study population genetics and evolutionary consequences of IS on the selected mutations themselves and on adjacent selectively neutral variation. We show that IS reduces levels of polymorphism and increases low-frequency variants and linkage disequilibrium, in both selected and adjacent neutral mutations. IS can account for several well-documented patterns of variation and composition in genomic regions with low rates of crossing over in Drosophila. IS cannot be described simply as a reduction in the efficacy of selection and effective population size in standard models of selection and drift. Rather, IS can be better understood with models that incorporate a constant "traffic" of competing alleles. Our simulations also allow us to make genome-wide predictions that are specific to IS. We show that IS will be more severe at sites in the center of a region containing weakly selected mutations than at sites located close to the edge of the region. Drosophila melanogaster genomic data strongly support this prediction, with genes without introns showing significantly reduced codon bias in the center of coding regions. As expected, if introns relieve IS, genes with centrally located introns do not show reduced codon bias in the center of the coding region. We also show that reasonably small differences in the length of intermediate "neutral" sequences embedded in a region under selection increase the effectiveness of selection on the adjacent selected sequences. Hence, the presence and length of sequences such as introns or intergenic regions can be a trait subject to selection in recombining genomes. In support of this prediction, intron presence is positively correlated with a gene's codon bias in D. melanogaster. Finally, the study of temporal dynamics of IS after a change of recombination rate shows that nonequilibrium codon usage may be the norm rather than the exception.

The evolutionary analysis of "orphans" from the Drosophila genome identifies rapidly diverging and incorrectly annotated genes

In genome projects of eukaryotic model organisms, a large number of novel genes of unknown function and evolutionary history ("orphans") are being identified. Since many orphans have no known homologs in distant species, it is unclear whether they are restricted to certain taxa or evolve rapidly, either because of a lack of constraints or positive Darwinian selection. Here we use three criteria for the selection of putatively rapidly evolving genes from a single sequence of Drosophila melanogaster. Thirteen candidate genes were chosen from the Adh region on the second chromosome and 1 from the tip of the X chromosome. We succeeded in obtaining sequence from 6 of these in the closely related species D. simulans and D. yakuba. Only 1 of the 6 genes showed a large number of amino acid replacements and in-frame insertions/deletions. A population survey of this gene suggests that its rapid evolution is due to the fixation of many neutral or nearly neutral mutations. Two other genes showed "normal" levels of divergence between species. Four genes had insertions/deletions that destroy the putative reading frame within exons, suggesting that these exons have been incorrectly annotated. The evolutionary analysis of orphan genes in closely related species is useful for the identification of both rapidly evolving and incorrectly annotated genes.

A statistical analysis of nucleotide substitutions in the Drosophila Adh region reflects irregularities in molecular clocks

Substitutions rates are expected to be rather constant when a gene is compared between species. To analyze this feature, Ka/Ks ratios have been studied for Alcohol dehydrogenase (Adh) and Alcohol dehydrogenase duplication (Adh-dup) genes in Drosophila species. Adh Ka/Ks values are lower in intrasubgenus comparisons involving species of the Sophophora group than when these species are compared to the D. immigrans and S. lebanonensis, and this difference does not occur in the Adh-dup comparisons.

Comparative analysis of the nonA region in Drosophila identifies a highly diverged 5' gene that may constrain nonA promoter evolution

A genomic fragment from Drosophila virilis that contained all the no-on-transientA (nonA) coding information, plus several kilobases of upstream material, was identified. Comparisons of nonA sequences and the gene nonA-like in D. melanogaster, a processed duplication of nonA, suggest that it arose before the split between D. melanogaster and D. virilis. In both species, another gene that lies <350 bp upstream from the nonA transcription starts, and that probably corresponds to the lethal gene l(1)i19, was identified. This gene encodes a protein that shows similarities to GPI1, which is required for the biosynthesis of glycosylphosphatidylinositol (GPI), a component for anchoring eukaryotic proteins to membranes, and so we have named it dGpi1. The molecular evolution of nonA and dGpi1 sequences show remarkable differences, with the latter revealing a level of amino acid divergence that is as high as that of transformer and with extremely low levels of codon bias. Nevertheless, in D. melanogaster hosts, the D. virilis fragment rescues the lethality associated with a mutation of l(1)i19e, as well as the viability and visual defects produced by deletion of nonA(-). The presence of dGpi1 sequences so close to nonA appears to have constrained the evolution of the nonA promoter.

Substitution rates in Drosophila nuclear genes: implications for translational selection

The relationships between synonymous and nonsynonymous substitution rates and between synonymous rate and codon usage bias are important to our understanding of the roles of mutation and selection in the evolution of Drosophila genes. Previous studies used approximate estimation methods that ignore codon bias. In this study we reexamine those relationships using maximum-likelihood methods to estimate substitution rates, which accommodate the transition/transversion rate bias and codon usage bias. We compiled a sample of homologous DNA sequences at 83 nuclear loci from Drosophila melanogaster and at least one other species of Drosophila. Our analysis was consistent with previous studies in finding that synonymous rates were positively correlated with nonsynonymous rates. Our analysis differed from previous studies, however, in that synonymous rates were unrelated to codon bias. We therefore conducted a simulation study to investigate the differences between approaches. The results suggested that failure to properly account for multiple substitutions at the same site and for biased codon usage by approximate methods can lead to an artifactual correlation between synonymous rate and codon bias. Implications of the results for translational selection are discussed.

Nonrandom spatial distribution of synonymous substitutions in the GP63 gene from Leishmania

In this work we analyze the variability in substitution rates in the GP63 gene from Leishmania. By using a sliding window to estimate substitution rates along the gene, we found that the rate of synonymous substitutions along the GP63 gene is highly correlated with both the rate of amino acid substitution and codon bias. Furthermore, we show that comparisons involving genes that represent independent phylogenetic lines yield very similar divergence/conservation patterns, thus suggesting that deterministic forces (i.e., nonstochastic forces such as selection) generated these patterns. We present evidence indicating that the variability in substitution rates is unambiguously related to functionally relevant features. In particular, there is a clear relationship between rates and the tertiary structure of the encoded protein since all divergent segments are located on the surface of the molecule and facing one side (almost parallel to the cell membrane) on the exposed surface of the organism. Remarkably, the protein segments encoded by these variable regions encircle the active site in a funnel-like distribution. These results strongly suggest that the pattern of nucleotide divergence and, notably, of synonymous divergence is affected by functional constraints.

Decoupled evolution of coding region and mRNA expression patterns after gene duplication: implications for the neutralist-selectionist debate

The neutralist perspective on molecular evolution maintains that the vast majority of mutations affecting gene function are neutral or deleterious. After a gene duplication where both genes are retained, it predicts that original and duplicate genes diverge at clock-like rates. This prediction is usually tested for coding sequences, but can also be applied to another important aspect of gene function, the genes' expression pattern. Moreover, if both sequence and expression pattern diverge at clock-like rates, a correlation between divergence in sequence and divergence in expression patterns is expected. Duplicate gene pairs with more highly diverged sequences should also show more highly diverged expression patterns. This prediction is tested for a large sample of duplicated genes in the yeast Saccharomyces cerevisiae, using both genome sequence and microarray expression data. Only a weak correlation is observed, suggesting that coding sequence and mRNA expression patterns of duplicate gene pairs evolve independently and at vastly different rates. Implications of this finding for the neutralist-selectionist debate are discussed.

Evidence for a high frequency of simultaneous double-nucleotide substitutions

Point mutations are generally assumed to involve changes of single nucleotides. Nevertheless, the nature and known mechanisms of mutation do not exclude the possibility that several adjacent nucleotides may change simultaneously in a single mutational event. Two independent approaches are used here to estimate the frequency of simultaneous double-nucleotide substitutions. The first examines switches between TCN and AGY (where N is any nucleotide and Y is a pyrimidine) codons encoding absolutely conserved serine residues in a number of proteins from diverse organisms. The second reveals double-nucleotide substitutions in primate noncoding sequences. These two complementary approaches provide similar high estimates for the rate of doublet substitutions, on the order of 0.1 per site per billion years.

Synonymous rates at the RpII215 gene of Drosophila: variation among species and across the coding region

The region encompassing the RpII215 gene that encodes the largest component of the RNA polymerase II complex (1889 amino acids) has been sequenced in Drosophila subobscura, D. madeirensis, D. guanche, and D. pseudoobscura. Nonsynonymous divergence estimates (Ka) indicate that this gene has a very low rate of amino acid replacements. Given its low Ka and constitutive expression, synonymous substitution rates are, however, unexpectedly high. Sequence comparisons have allowed the molecular clock hypothesis to be tested. D. guanche is an insular species and it is therefore expected to have a reduced effective size relative to D. subobscura. The significantly higher rate of synonymous substitutions detected in the D. guanche lineage could be explained if synonymous mutations behave as nearly neutral. Significant departure from the molecular clock hypothesis for synonymous and nonsynonymous substitutions was detected when comparing the D. subobscura, D. pseudoobscura, and D. melanogaster lineages. Codon bias and synonymous divergence between D. subobscura and D. melanogaster were negatively correlated across the RpII215 coding region, which indicates that selection coefficients for synonymous mutations vary across the gene. The C-terminal domain (CTD) of the RpII215 protein is structurally and functionally differentiated from the rest of the protein. Synonymous substitution rates were significantly different in both regions, which strongly indicates that synonymous mutations in the CTD and in the non-CTD regions are under detectably different selection coefficients.

Natural selection on synonymous sites is correlated with gene length and recombination in Drosophila

Evolutionary analysis of codon bias in Drosophila indicates that synonymous mutations are not neutral, but rather are subject to weak selection at the translation level. Here we show that the effectiveness of natural selection on synonymous sites is strongly correlated with the rate of recombination, in accord with the nearly neutral hypothesis. This correlation, however, is apparent only in genes encoding short proteins. Long coding regions have both a lower codon bias and higher synonymous substitution rates, suggesting that they are affected less efficiently by selection. Therefore, both the length of the coding region and the recombination rate modulate codon bias. In addition, the data indicate that selection coefficients for synonymous mutations must vary by a minimum of one or two orders of magnitude. Two hypotheses are proposed to explain the relationship among the coding region length, the codon bias, and the synonymous divergence and polymorphism levels across the range of recombination rates in Drosophila. The first hypothesis is that selection coefficients on synonymous mutations are inversely related to the total length of the coding region. The second hypothesis proposes that interference among synonymous mutations reduces the efficacy of selection on these mutations. We investigated this second hypothesis by carrying out forward simulations of weakly selected mutations in model populations. These simulations show that even with realistic recombination rates, this interference, which we call the "small-scale" Hill-Robertson effect, can have a moderately strong influence on codon bias.