Comparative genomics and the study of evolution by natural selection

Genomics profoundly affects most areas of biology, including ecology and evolutionary biology. By examining genome sequences from multiple species, comparative genomics offers new insight into genome evolution and the way natural selection moulds DNA sequence evolution. Functional divergence, as manifested in the accumulation of nonsynonymous substitutions in protein-coding genes, differs among lineages in a manner seemingly related to population size. For example, the ratio of nonsynonymous to synonymous substitution (d_N/d_S) is higher in apes than in rodents, compatible with Ohta's nearly neutral theory of molecular evolution, which suggests that the fixation of slightly deleterious mutations contributes to protein evolution at an extent negatively correlated with effective population size. While this supports the idea that functional evolution is not necessarily adaptive, comparative genomics is uncovering a role for positive Darwinian selection in 10–40% of all genes in different lineages, estimates that are likely to increase when the addition of more genomes gives increased power. Again, population size seems to matter also in this context, with a higher proportion of fixed amino acid changes representing advantageous mutations in large populations. Genes that are particularly prone to be driven by positive selection include those involved with reproduction, immune response, sensory perception and apoptosis. Genetic innovations are also frequently obtained by the gain or loss of complete gene sequences. Moreover, it is increasingly realized, from comparative genomics, that purifying selection conserves much more than just the protein-coding part of the genome, and this points at an important role for regulatory elements in trait evolution. Finally, genome sequencing using outbred or multiple individuals has provided a wealth of polymorphism data that gives information on population history, demography and marker evolution.     

Sexual Selection and the Molecular Evolution of ADAM Proteins

Rapid evolution has been identified for many reproductive genes and recent studies have combined phylogenetic tests and information on species mating systems to test sexual selection. Here we examined the molecular evolution of the ADAM gene family, a diverse group of 35 proteins capable of adhesion to and cleavage of other proteins, using sequence data from 25 mammalian genes. Out of the 25 genes analyzed, all those expressed in male reproductive tissue showed evidence of positive selection. Positively selected amino acids within the protein adhesion domain were only found in sperm surface ADAM proteins (ADAMs 1, 2, 3, 4, and 32) suggesting selection driven by male × female interactions. We tested heterogeneity in rates of evolution of the adhesion domain of ADAM proteins by using sequence data from Hominidae and macaques. The use of the branch and branch-site models (PAML) showed evidence of higher d _N/d _S and/or positive selection linked to branches experiencing high postmating selective pressures (chimpanzee and macaque) for Adams 2, 18, and 23. Moreover, we found consistent higher proportion of nonsynonymous relative to synonymous and noncoding sequence substitutions in chimpanzee and/or macaque only for Adams 2, 18, and 23. Our results suggest that lineage-specific sexual selection bouts might have driven the evolution of the adhesion sperm protein surface domains of ADAMs 2 and 18 in primates. Adams 2 and 18 are localized in chromosome 8 of primates and adjacent to each other, so their evolution might have also been influenced by their common genome localization.     

Domain-Specific Positive Selection Contributes to the Evolution of Arabidopsis Leucine-Rich Repeat Receptor-Like Kinase (LRR RLK) Genes

Leucine-rich repeat receptor-like kinases (LRR RLKs) comprise the largest group within the plant receptor-like kinase (RLK) superfamily, and the Arabidopsis genome alone contains over 200 LRR RLK genes. Although there is clear evidence for diverse roles played by individual LRR RLK genes in Arabidopsis growth and development, the evolutionary mechanism for this functional diversification is currently unclear. In this study, we focused on the LRRII RLK subfamily to investigate the molecular mechanisms that might have led to the functional differentiation of Arabidopsis LRR RLK genes. Phylogenetic analysis of 14 genes in this subfamily revealed three well-supported groups (I, II, and III). RT-PCR analysis did not find many qualitative differences in expression among these 14 genes in various Arabidopsis tissues, suggesting that evolution of regulatory sequences did not play a major role in their functional divergence. We analyzed substitution patterns in the predicted ligand-binding regions of these genes to examine if positive selection has acted to produce novel ligand-binding specificities, using the nonsynonymous/synonymous rate ratio (d _N/d _S) as an indicator of selective pressure. Estimates of d _N/d _S ratios from multiple methods indicate that nonsynonymous substitutions accumulated during divergence of the three lineages. Positive selection is likely to have occurred along the lineages ancestral to groups II and III. We suggest that positive selection on the ligand-binding sites of LRRII RLKs promoted diversification of ligand-binding specificities and thus contributed to the functional differentiation of Arabidopsis LRRII RLK genes during evolution. [Reviewing Editor: Dr. Martin Kreitman]     

A model for continuous production of thermally induced recombinant proteins

Summary General conditions for continuous expression of heterologous genes fromP _L promoter in two fermenters connected in series have been established. The induction time of the bacterial cells is calculated as a function of the retention time in the inducing reactor. Using this model, it is possible to adapt fermentation parameters to the particular behaviour of any specific recombinant clone.Nomenclature F(t) flow at timet [ml/min] - M _T (t) culture induced, at timeT of fermentation, during a period up tot [ml] - N _T (t) culture induced, at timeT of fermentation, during a period fromt tot+dt [ml] - p(t) product yield in a discontinuous culture [units/ml] - P(t) product yield at the outlet of the fermenter [units/ml] - v(t) volume of culture entered into the inducing reactor up to timet [ml] - V volume of the inducing reactor [ml] Greek letters retention time in the inducing reactor [min] - (t) average induction time at timet [min]     

The effects of recombination,mutation and selection on the evolution of the Rp1 resistance genes in grasses

Plant immune genes, or resistance genes, are involved in a co‐evolutionary arms race with a diverse range of pathogens. In agronomically important grasses, such R genes have been extensively studied because of their role in pathogen resistance and in the breeding of resistant cultivars. In this study, we evaluate the importance of recombination, mutation and selection on the evolution of the R gene complex Rp1 of Sorghum, Triticum, Brachypodium, Oryza and Zea. Analyses show that recombination is widespread, and we detected 73 independent instances of sequence exchange, involving on average 1567 of 4692 nucleotides analysed (33.4%). We were able to date 24 interspecific recombination events and found that four occurred postspeciation, which suggests that genetic introgression took place between different grass species. Other interspecific events seemed to have been maintained over long evolutionary time, suggesting the presence of balancing selection. Significant positive selection (i.e. a relative excess of nonsynonymous substitutions (d_N/d_S>1)) was detected in 17–95 codons (0.42–2.02%). Recombination was significantly associated with areas with high levels of polymorphism but not with an elevated d_N/d_S ratio. Finally, phylogenetic analyses show that recombination results in a general overestimation of the divergence time (mean = 14.3%) and an alteration of the gene tree topology if the tree is not calibrated. Given that the statistical power to detect recombination is determined by the level of polymorphism of the amplicon as well as the number of sequences analysed, it is likely that many studies have underestimated the importance of recombination relative to the mutation rate.     

Evolutionary dynamics of Rh2 opsins in birds demonstrate an episode of accelerated evolution in the New World warblers (Setophaga)

Low rates of sequence evolution associated with purifying selection can be interrupted by episodic changes in selective regimes. Visual pigments are a unique system in which we can investigate the functional consequences of genetic changes, therefore connecting genotype to phenotype in the context of natural and sexual selection pressures. We study the RH2 and RH1 visual pigments (opsins) across 22 bird species belonging to two ecologically convergent clades, the New World warblers (Parulidae) and Old World warblers (Phylloscopidae) and evaluate rates of evolution in these clades along with data from 21 additional species. We demonstrate generally slow evolution of these opsins: both Rh1 and Rh2 are highly conserved across Old World and New World warblers. However, Rh2 underwent a burst of evolution within the New World genus Setophaga, where it accumulated substitutions at 6 amino acid sites across the species we studied. Evolutionary analyses revealed a significant increase in d_N/d_S in Setophaga, implying relatively strong selective pressures to overcome long‐standing purifying selection. We studied the effects of each substitution on spectral tuning and found they do not cause large spectral shifts. Thus, substitutions may reflect other aspects of opsin function, such as those affecting photosensitivity and/or dark–light adaptation. Although it is unclear what these alterations mean for colour perception, we suggest that rapid evolution is linked to sexual selection, given the exceptional plumage colour diversification in Setophaga.     

Genetic variation of geminiviruses: comparison between sexual and asexual host plant populations

One of the most promising hypotheses for the evolution of sex is that sexual reproduction is advantageous because it increases the rate of adaptive evolution in response to parasites. To investigate this advantage of sex, we compared genetic variation of geminiviruses infecting sexual and asexual populations of Eupatorium (Asteraceae). The infection frequency was 37.5% in the sexual population and 87.8% in the asexual population. The lower infection frequency in the sexual population might be the result of higher genetic diversity of host plants. If geminiviruses have diverged to counter defence systems of genetically variable hosts, genetic diversity of viruses is expected to be higher in sexual host populations than in asexual host populations. To test this expectation, we used single-strand conformation polymorphism (SSCP) analysis to examine genetic diversity of the geminiviruses in a DNA region containing the open-reading frame (ORF) C4 gene, which is known to function as a host range determinant. As predicted, higher genetic diversity of viruses was observed in the sexual population: three SSCP types were found in the asexual population while six types were found in the sexual population. Sequencing of the polymerase chain reaction (PCR) products revealed further genetic diversity. Phylogenetic analysis of the sequences showed that the SSCP types belonged to four different clades. Several SSCP types from the same clade were found in the sexual population, whereas the asexual population included only one SSCP type from each clade. Amino acid replacements of ORF C4 are suggested to be accelerated in the sexual population. This evidence supports the hypothesis that sexual reproduction is advantageous as a defence against epidemic disease.     

Minimum evolution using ordinary least-squares is less robust than neighbor-joining

The method of minimum evolution reconstructs a phylogenetic tree T for n taxa given dissimilarity data d. In principle, for every tree W with these n leaves an estimate for the total length of W is made, and T is selected as the W that yields the minimum total length. Suppose that the ordinary least-squares formula S _W (d) is used to estimate the total length of W. A theorem of Rzhetsky and Nei shows that when d is positively additive on a completely resolved tree T, then for all W ≠ T it will be true that S _W (d) > S _T (d). The same will be true if d is merely sufficiently close to an additive dissimilarity function. This paper proves that as n grows large, even if the shortest branch length in the true tree T remains constant and d is additive on T, then the difference S _W (d)-S _T (d) can go to zero. It is also proved that, as n grows large, there is a tree T with n leaves, an additive distance function d _T on T with shortest edge ε, a distance function d, and a tree W with the same n leaves such that d differs from d _T by only approximately ε/4, yet minimum evolution incorrectly selects the tree W over the tree T. This result contrasts with the method of neighbor-joining, for which Atteson showed that incorrect selection of W required a deviation at least ε/2. It follows that, for large n, minimum evolution with ordinary least-squares can be only half as robust as neighbor-joining.     

A Genomic Comparison of Faster-Sex, Faster-X, and Faster-Male Evolution Between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d_N) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d_N values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d_N. Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation. [Reviewing Editor: Dr. Willie J. Swanson]     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

Genetic tests of ancient asexuality in Root Knot Nematodes reveal recent hybrid origins

Background  

The existence of "ancient asexuals", taxa that have persisted for long periods of evolutionary history without sexual recombination, is both controversial and important for our understanding of the evolution and maintenance of sexual reproduction. A lack of sex has consequences not only for the ecology of the asexual organism but also for its genome. Several genetic signatures are predicted from long-term asexual (apomictic) reproduction including (i) large "allelic" sequence divergence (ii) lack of phylogenetic clustering of "alleles" within morphological species and (iii) decay and loss of genes specific to meiosis and sexual reproduction. These genetic signatures can be hard to assess since it is difficult to demonstrate the allelic nature of very divergent sequences, divergence levels may be complicated by processes such as inter-specific hybridization, and genes may have secondary roles unrelated to sexual reproduction. Apomictic species of Meloidogyne root knot nematodes have been suggested previously to be ancient asexuals. Their relatives reproduce by meiotic parthenogenesis or facultative sexuality, which in combination with the abundance of nematode genomic sequence data, makes them a powerful system in which to study the consequences of reproductive mode on genomic divergence.     

Molecular Phylogenies and Evolution of crt Genes in Algae

ABSTRACT

The carotenoids constitute the most widespread class of pigments in nature. Most previous work has concentrated on the identification and characterization of their chemical physical properties and bioavailability. In recent years, significant amounts of research have been conducted in an attempt to analyze the genes and the molecular regulation of the genes involved in the biosynthesis of carotenoids. However, it is important not to lose sight of the early evolution of carotenoid biosynthesis. One of the major obstacles in understanding the evolution of the respective enzymes and their patterns of selection is a lack of a well-supported phylogenic analysis. In the present research, a major long-term objective was to provide a clearer picture of the evolutionary history of genes, together with an evaluation of the patterns of selection in algae. These phylogenies will be important in studies characterizing the evolution of algae. The gene sequences of the enzymes involved in the major steps of the carotenoid biosynthetic pathway in algae (cyanobacteria, rhofophyta, chlorophyta) have been analyzed. Phylogenetic relationships among protein-coding DNA sequences were reconstructed by neighbor-joining (NJ) analysis for the respective carotenoid biosynthetic pathway genes (crt) in algae. The analysis also contains an estimation of the rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d_N), synonymous nucleotide substitution per synonymous site (d_S), and the ratio of nonsynonmous (d_N/d_S) for the test of selection patterns. The phylogenetic trees show that the taxa of some genera have a closer evolutionary relationship with other genera in some gene sequences, which suggests a common ancient origin and that lateral gene transfer has occurred among unrelated genera. The d_N values of crt genes in the early pathway are relatively low, while those of the following steps are slightly higher, while the d_N values of crt genes in chlorophyta are higher than those in cyanobacteria. Most of the d_N/d_S values exceed 1. The phylogenetic analysis revealed that lateral gene transfer may have taken place across algal genomes and the d_N values suggest that most of the early crt genes are well conserved compared to the later crt genes. Furthermore, d_N values also revealed that the crt genes of chlorophyta are more evolutionary than cyanobacteria. The amino acids' changes are mostly adaptive evolution under the influence of positive diversity selection.     

Energetic constraints on sexual activity in the male edible dormouse (Glis glis)

The aim of this study was to examine to what extent reproductive activity in male edible dormice (Glis glis) might be energetically constrained. Demographic data, morphometric data, and oral body temperature (T _or) measurements were collected in two study areas between 1993 and 2002 in southwest Germany and combined with subcutaneous body temperature (T _sc) registrations of captive dormice. T _sc measurements were collected directly after emergence from hibernation (June) until the end of the mating season (July). Wild edible dormice showed strong fluctuations in their reproductive output between years. Not all males were sexually active each year and the number of litters born was positively correlated with the number of sexually active males, which suggests that sexual activity in males is constrained and in turn limits reproductive success. A comparison of the T _or of sexually quiescent and active males revealed that sexually quiescent males had significantly lower T _or (median: 28.8°C; 25/75% quartiles: 16.4/31.0; n=31) than sexually active males (median: 34.2°C; 25/75% quartiles: 32.0/35.6; n=156). Body condition of sexually active and quiescent males was not different after emergence from hibernation. However, sexually active males showed a significant reduction in their body condition between June and July, the time of mating, while body condition of sexually quiescent males remained constant. Continuous T _sc registrations in captive sexually active male dormice showed strong circadian T _sc fluctuations. Even though daily torpor bouts with T _sc below 20°C occurred in these males, most of the time T _sc fluctuated above 30°C, which is known as the critical body temperature threshold above which testes maturation can take place in this species. These results demonstrate that male dormice incur high costs due to sexual activity and that thermoregulation is determined by a trade-off between energetic savings and reproductive activity.     

A mutant with constitutive sexual organ development in<Emphasis Type="Italic"> Marchantia polymorpha</Emphasis> L.

In lower land plants, genes controlling the transition from vegetative growth to sexual reproduction have not yet been identified. In the dioecious liverwort Marchantia polymorpha, the transition to sexual reproduction accompanied by the formation of sexual organs on the gametophytic thallus is initiated under long-day conditions. By particle bombardment-mediated mutagenesis, we generated a mutant of M. polymorpha that constitutively forms sexual organs. This mutant is fully fertile, showing that the mutation does not affect formation of male or female sexual organs per se. Genetic analysis reveals that this phenotype is caused by mutation of a single autosomal locus, suggesting that this mutation defines or controls a gene regulating the transition to sexual reproduction in M. polymorpha.     

Functional bias in molecular evolution rate of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Background  

Characteristics derived from mutation and other mechanisms that are advantageous for survival are often preserved during evolution by natural selection. Some genes are conserved in many organisms because they are responsible for fundamental biological function, others are conserved for their unique functional characteristics. Therefore one would expect the rate of molecular evolution for individual genes to be dependent on their biological function. Whether this expectation holds for genes duplicated by whole genome duplication is not known.     

Sexual selection and maintenance of sex: evidence from comparisons of rates of genomic accumulation of mutations and divergence of sex-related genes in sexual and hermaphroditic species of Caenorhabditis

Several hypotheses have been proposed to explain the persistenceof dioecy despite the reproductive advantages conferred to hermaphrodites,including greater efficiency at purging deleterious mutationsin the former. Dioecy can benefit from both mutation purgingand accelerated evolution by bringing together beneficial mutationsin the same individual via recombination and shuffling of genotypes.In addition, mathematical treatment has shown that sexual selectionis also capable of mitigating the cost of maintaining separatesexes by increasing the overall fitness of sexual populations,and genomic comparisons have shown that sexual selection canlead to accelerated evolution. Here, we examine the advantagesof dioecy versus hermaphroditism by comparing the rate of evolutionin sex-related genes and the rate of accumulation of deleteriousmutations using a large number of orthologs (11,493) in thedioecious Caenorhabditis remanei and the hermaphroditic Caenorhabditisbriggsae. We have used this data set to estimate the deleteriousmutation rate per generation, U, in both species and find thatalthough it is significantly higher in hermaphrodites, bothspecies are at least 2 orders of magnitude lower than the valuerequired to explain the persistence of sex by efficiency atpurging deleterious mutations alone. We also find that genesexpressed in sperm are evolving rapidly in both species; however,they show a greater increase in their rate of evolution relativeto genes expressed in other tissues in C. remanei, suggestingstronger sexual selection pressure acting on these genes indioecious species. Interestingly, the persistence of a signalof rapid evolution of sperm genes in C. briggsae suggests arecent evolutionary origin of hermaphrodism in this lineage.Our results provide empirical evidence of increased sexual selectionpressure in dioecious animals, supporting the possibility thatsexual selection may play an important role in the maintenanceof sexual reproduction.     

Patterns and impacting factors of gene evolutionary rate between wild and cultivated emmer wheat (Triticum turgidum)

Tetraploid emmer wheat (Triticum turgidum L., BBAA) is the founder progenitor of bread wheat, providing the valuable genetic resource and gene pool for wheat improvement. However, the evolutionary trajectory of tetraploid wheat, especially the evolutionary fate of different types of genes has not been well studied. In this study, the rate of non-synonymous substitution (d_N) and synonymous substitution (d_S) was calculated by comparing the orthologs between the wild emmer and cultivated durum wheat at the whole genome and subgenome levels to obtain the positively selected genes (PSGs) and negatively selected genes (NSGs). Then, mutation rate, gene length, exon number, GC content, codon bias, and expression level were comprehensively investigated and compared between the PSGs and NSGs. Within both wild emmer and cultivated durum wheat, PSGs between A and B subgenome displayed shorter gene and exon lengths as well as fewer exon numbers compared with NSGs, whereas from wild emmer to cultivated durum wheat, PSGs showed longer gene length and more exon numbers. Furthermore, PSGs displayed much higher expression levels and stronger codon usage bias, but lower genetic diversity compared with NSGs. Finally, two PSGs TdER1-6B, and TdLC7-2A, were found to play the crucial roles in regulating grain width and plant height of tetraploid wheat, respectively. This study systematically investigated the evolutionary, structural, and functional difference between PSGs and NSGs in tetraploid wheat, which will contribute to a better understanding of the selective mode and evolutionary trajectory during wheat domestication and evolution.     

Stronger selective constraint on downstream genes in the oxidative phosphorylation pathway of cetaceans

The oxidative phosphorylation (OXPHOS) pathway is an efficient way to produce energy via adenosine triphosphate (ATP), which is critical for sustaining an energy supply for cetaceans in a hypoxic environment. Several studies have shown that natural selection may shape the evolution of the genes involved in OXPHOS. However, how network architecture drives OXPHOS protein sequence evolution remains poorly explored. Here, we investigated the evolutionary patterns of genes in the OXPHOS pathway across six cetacean genomes within the framework of a functional network. Our results show a negative correlation between the strength of purifying selection and pathway position. This result indicates that downstream genes were subjected to stronger evolutionary constraints than upstream genes, which may be due to the dual function of ATP synthase in the OXPHOS pathway. Additionally, there was a positive correlation between codon usage bias and omega (ω = dN/dS) and a negative correlation with synonymous substitution rate (dS), indicating that the stronger selective constraint on genes (with less biased codon usage) along the OXPHOS pathway is attributable to an increase in the rate of synonymous substitution. Surprisingly, there was no significant correlation between protein–protein interactions and the evolutionary estimates, implying that highly connected enzymes may not always show greater evolutionary constraints. Compared with that observed for terrestrial mammals, we found that the signature of positive selection detected in five genes (ATP5J, LHPP, PPA1, UQCRC1 and UQCRQ) was cetacean‐specific, reflecting the importance of OXPHOS for survival in hypoxic, aquatic environments.