Higher Rates Of Molecular Evolution Research Articles

Transcription Factors Evolve Faster Than Their Structural Gene Targets in the Flavonoid Pigment Pathway.

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

The effect of landscape on Cucurbita pepo-pollinator interaction networks varies depending on plants\u2019 genetic diversity

In recent years, evidence has been found that plant-pollinator interactions are altered by land-use and that genetic diversity also plays a role. However, how land-use and genetic diversity influence plant–pollinator interactions, particularly in the Neotropics, where many endemic plants exist is still an open question. Cucurbita pepo is a monoecious plant and traditional crop wide distributed, with high rates of molecular evolution, landraces associated with human cultural management and a history of coevolution with bees, which makes this species a promising model for studying the effect of landscape and genetic diversity on plant-pollinator interactions. Here, we assess (1) whether female and male flowers differences have an effect on the interaction network, (2) how C. pepo genetic diversity affects flower-bee visitation network structure, and (3) what is the effect that land-use, accounting for C. pepo genetic variability, has on pumpkin-bee interaction network structure. Our results indicate that female and male flowers presented the same pollinator community composition and interaction network structure suggesting that female/male differences do not have a significant effect on network evolution. Genetic diversity has a positive effect on modularity, nestedness and number of interactions. Further, the effect of semi-natural areas on nestedness could be buffered when genetic diversity is high. Our results suggest that considering genetic diversity is relevant for a better understanding of the effect of land-use on interaction networks. Additionally, this understanding has great value in conserving biodiversity and enhancing the stability of interaction networks in a world facing great challenges of habitat and diversity loss.

Genetic insights into the evolution of genera with the eastern Asia-eastern North America floristic disjunction: a transcriptomics analysis.

Large disjunctions in species distributions provide excellent opportunities to study processes that shape biogeographic patterns. One such disjunction is the eastern Asia-eastern North America (EA-ENA) floristic disjunction. For many genera with this disjunction, species richness is greater in EA than in ENA; this pattern has been attributed, in part, to higher rates of molecular evolution and speciation in EA. Longer branch lengths have been found in some EA clades, relative to their ENA sister clades, suggesting that the EA lineages have evolved at a higher rate, possibly due to environmental heterogeneity, potentially contributing to the species richness anomaly. To evaluate whether rates of molecular evolution are elevated in EA relative to ENA, we used transcriptomes from species in 11 genera displaying this disjunction. Rates of molecular evolution were estimated for up to 385 orthologous nuclear loci per genus. No statistically significant differences were identified in pairwise comparisons between EA and ENA sister species, suggesting equal rates of molecular evolution for both species; the data also suggest similar selection pressures in both regions. For larger genera, evidence likewise argues against more species-rich clades having higher molecular evolutionary rates, regardless of region. Our results suggest that genes across multiple gene ontology categories are evolving at similar rates under purifying selection in species in both regions. Our data support the hypothesis that greater species richness in EA than ENA is due to factors other than an overall increase in rates of molecular evolution in EA.

A Cambrian-Ordovician Terrestrialization of Arachnids.

Understanding the temporal context of terrestrialization in chelicerates depends on whether terrestrial groups, the traditional Arachnida, have a single origin and whether or not horseshoe crabs are primitively or secondarily marine. Molecular dating on a phylogenomic tree that recovers arachnid monophyly, constrained by 27 rigorously vetted fossil calibrations, estimates that Arachnida originated during the Cambrian or Ordovician. After the common ancestor colonized the land, the main lineages appear to have rapidly radiated in the Cambrian–Ordovician boundary interval, coinciding with high rates of molecular evolution. The highest rates of arachnid diversification are detected between the Permian and Early Cretaceous. A pattern of ancient divergence estimates for terrestrial arthropod groups in the Cambrian while the oldest fossils are Silurian (seen in both myriapods and arachnids) is mirrored in the molecular and fossil records of land plants. We suggest the discrepancy between molecular and fossil evidence for terrestrialization is likely driven by the extreme sparseness of terrestrial sediments in the rock record before the late Silurian.

Sexual Selection Shapes Seminal Vesicle Secretion Gene Expression in House Mice.

Reproductive proteins typically have high rates of molecular evolution, and are assumed to be under positive selection from sperm competition and cryptic female choice. However, ascribing evolutionary divergence in the genome to these processes of sexual selection from patterns of association alone is problematic. Here, we use an experimental manipulation of postmating sexual selection acting on populations of house mice and explore its consequences for the expression of seminal vesicle secreted (SVS) proteins. Following 25 generations of selection, males from populations subjected to postmating sexual selection had evolved increased expression of at least two SVS genes that exhibit the signature of positive selection at the molecular level, SVS1 and SVS2. These proteins contribute to mating plug formation and sperm survival in the female reproductive tract. Our data thereby support the view that sexual selection is responsible for the evolution of these seminal fluid proteins.

Population Genomic Analysis of a Pitviper Reveals Microevolutionary Forces Underlying Venom Chemistry.

Venoms are among the most biologically active secretions known, and are commonly believed to evolve under extreme positive selection. Many venom gene families, however, have undergone duplication, and are often deployed in doses vastly exceeding the LD50 for most prey species, which should reduce the strength of positive selection. Here, we contrast these selective regimes using snake venoms, which consist of rapidly evolving protein formulations. Though decades of extensive studies have found that snake venom proteins are subject to strong positive selection, the greater action of drift has been hypothesized, but never tested. Using a combination of de novo genome sequencing, population genomics, transcriptomics, and proteomics, we compare the two modes of evolution in the pitviper, Protobothrops mucrosquamatus. By partitioning selective constraints and adaptive evolution in a McDonald–Kreitman-type framework, we find support for both hypotheses: venom proteins indeed experience both stronger positive selection, and lower selective constraint than other genes in the genome. Furthermore, the strength of selection may be modulated by expression level, with more abundant proteins experiencing weaker selective constraint, leading to the accumulation of more deleterious mutations. These findings show that snake venoms evolve by a combination of adaptive and neutral mechanisms, both of which explain their extraordinarily high rates of molecular evolution. In addition to positive selection, which optimizes efficacy of the venom in the short term, relaxed selective constraints for deleterious mutations can lead to more rapid turnover of individual proteins, and potentially to exploration of a larger venom phenotypic space.

High Rate of Protein Coding Sequence Evolution and Species Diversification in the Lycaenids

Understanding the variations in the rate of molecular evolution among lineages may provide clues on the processes that molded extant biodiversity. Here, we report the high rate of molecular evolution in Lycaenidae and Riodinidae compared to other families of butterflies (Papilionoidea). We assembled a phylogeny of butterflies using eight molecular markers and comprising 4891 species. We found that the rate of molecular evolution is higher in Lycaenidae and Riodinidae compared to the other families, but only the nuclear gene Wingless showed a marked difference, while Elongation factor 1-alpha showed a more moderate difference. In contrast, the gene Cytochrome Oxidase subunit 1 showed no difference between lycaenids as well as riodinids, and other butterflies. In parallel, we calculated the rates of diversification in all subfamilies of the Papilionoidea using the method-of-moments estimator for stem-group ages, which does not require a fully solved phylogeny for the target clades. We found that the Nemeobiinae and Lycaeninae, from the Riodinidae and Lycaenidae families, respectively, had the highest rate of diversification among all subfamilies. Among the life-history traits that could explain differences in molecular evolution and diversification rate, lycaenids display mutualistic or antagonistic interactions with ants, a higher level of host plant specialization and reduced dispersal abilities compared to other butterfly families. Since the current study is limited by the unique event of apparition of myrmecophily analyzed, the relationship with traits cannot be evaluated statistically. Future studies should measure myrmecophily and dispersal abilities quantitatively across a more detailed phylogeny of lycaenids to test for an association between shifts in the strength of mutualism, rates of molecular evolution and the diversification of lineages.

The evolution of animal Argonautes: evidence for the absence of antiviral AGO Argonautes in vertebrates

In addition to mediating regulation of endogenous gene expression, RNA interference (RNAi) in plants and invertebrates plays a crucial role in defense against viruses via virus-specific siRNAs. Different studies have demonstrated that the functional diversity of RNAi in animals is linked to the diversification of the Argonaute superfamily, central components of RISCs (RNA induced silencing complexes). The animal Argonaute superfamily is traditionally grouped into AGO and PIWI Argonautes. Yet, by performing phylogenetic analyses and determining the selective evolutionary pressure in the metazoan Argonaute superfamily, we provide evidence for the existence of three conserved Argonaute lineages between basal metazoans and protostomes, namely siRNA-class AGO, miRNA-class AGO and PIWI Argonautes. In addition, it shown that the siRNA-class AGO lineage is characterized by high rates of molecular evolution, suggesting a role in the arms race with viruses, while the miRNA-class AGOs display strong sequence conservation. Interestingly, we also demonstrate that vertebrates lack siRNA-class AGO proteins and that vertebrate AGOs display low rates of molecular evolution. In this way, we provide supportive evidence for the loss of the antiviral siRNA-class AGO group in vertebrates and discuss the consequence hereof on antiviral immunity and the use of RNAi as a loss of function tool in these animals.

Effects of gene choice, base composition and rate heterogeneity on inference and estimates of divergence times in cypriniform fishes

Gene choice, base compositional heterogeneity and rate heterogeneity are known to influence phylogenetic reconstruction but are often overlooked. Here, we investigate the impact of these factors in reconstructing the phylogenetic relationships and divergence times of cypriniform fishes. We compile a multilocus data set of newly sequenced and previously published nuclear protein-coding genes from species from all major lineages of Cypriniformes. We collect data on the oldest known cypriniform fossils and gather body size information. We demonstrate that body size correlates with rates of molecular evolution and contributes to rate heterogeneity not only in Cypriniformes but also in a number of other fish lineages. We find that miniaturized cypriniforms show exceptionally high rates of molecular evolution that may lead to improper phylogenetic placement of these taxa due to issues associated with branch attraction. We identify and correct for base compositional heterogeneity and find that this has a marked impact on topology and a corresponding impact on the estimation of divergence times within Cypriniformes. Using gene sequences, fossil calibrations and two different methods, we show strong evidence that Catostomidae is sister to all other cypriniform families and the major clades of cypriniform fishes diverged long before their fossil record indicates.

Molecular organization of 5S rDNA in two ukrainian populations of sycamore (Acer pseudoplatanus)

Aim. The genetic variability of tree species has still not been studied enough at the molecular level. Considering the relatively low rate of molecular evolution in perennial tree species, it is necessary to use only those regions of genome, which demonstrate a high level of variability, such as 5S rDNA. Accordingly, to estimate the potential of 5S rDNA as a molecular marker for tree species, the organization of this genomic region was compared between samples from two geographically remote Ukrainian populations of sycamore, Acer pseudoplatanus. Methods. PCR amplification, cloning and sequencing of the 5S rDNA IGS of A.pseudoplatanus. Results. It was shown that only one variant of 5S rDNA repeat with a length of 475 bp is present in the genome of A. pseudoplatanus. Also, it was found that the elements typical for the angiosperm RNA polymerase III promoter, which are localized in IGS, appear to be different from those previously described for species of other families. The level of IGS sequence similarity within the populations exceeds 99%, while the level of IGS sequence similarity between various populations is only 93,3-94,3%. Conclusions. High rate of molecular evolution of the 5S rDNA IGS makes them a convenient molecular marker for evaluation of intraspecific variation in populations of A. pseudoplatanus.Key words: 5S rDNA, molecular markers, intraspecific variability, Acer

Do saline taxa evolve faster? Comparing relative rates of molecular evolution between freshwater and marine eukaryotes.

The major branches of life diversified in the marine realm, and numerous taxa have since transitioned between marine and freshwaters. Previous studies have demonstrated higher rates of molecular evolution in crustaceans inhabiting continental saline habitats as compared with freshwaters, but it is unclear whether this trend is pervasive or whether it applies to the marine environment. We employ the phylogenetic comparative method to investigate relative molecular evolutionary rates between 148 pairs of marine or continental saline versus freshwater lineages representing disparate eukaryote groups, including bony fish, elasmobranchs, cetaceans, crustaceans, mollusks, annelids, algae, and other eukaryotes, using available protein-coding and noncoding genes. Overall, we observed no consistent pattern in nucleotide substitution rates linked to habitat across all genes and taxa. However, we observed some trends of higher evolutionary rates within protein-coding genes in freshwater taxa-the comparisons mainly involving bony fish-compared with their marine relatives. The results suggest no systematic differences in substitution rate between marine and freshwater organisms.

Sperm Bindin Divergence under Sexual Selection and Concerted Evolution in Sea Stars.

Selection associated with competition among males or sexual conflict between mates can create positive selection for high rates of molecular evolution of gamete recognition genes and lead to reproductive isolation between species. We analyzed coding sequence and repetitive domain variation in the gene encoding the sperm acrosomal protein bindin in 13 diverse sea star species. We found that bindin has a conserved coding sequence domain structure in all 13 species, with several repeated motifs in a large central region that is similar among all sea stars in organization but highly divergent among genera in nucleotide and predicted amino acid sequence. More bindin codons and lineages showed positive selection for high relative rates of amino acid substitution in genera with gonochoric outcrossing adults (and greater expected strength of sexual selection) than in selfing hermaphrodites. That difference is consistent with the expectation that selfing (a highly derived mating system) may moderate the strength of sexual selection and limit the accumulation of bindin amino acid differences. The results implicate both positive selection on single codons and concerted evolution within the repetitive region in bindin divergence, and suggest that both single amino acid differences and repeat differences may affect sperm-egg binding and reproductive compatibility.

Patterns of molecular evolution of RNAi genes in social and socially parasitic bumblebees

The high frequency of interactions amongst closely related individuals in social insect colonies enhances pathogen transmission. Group-mediated behavior supporting immune defenses tends to decrease selection acting on immune genes. Along with low effective population sizes this might result in relaxed constraint and rapid evolution of immune system genes. Here, we show that antiviral siRNA genes show high rates of molecular evolution with argonaute 2, armitage and maelstrom evolving faster in social bumblebees compared to their socially parasitic cuckoo bumblebees that lack a worker caste. RNAi genes show frequent positive selection at the codon level additionally supported by the occurrence of parallel evolution. Their evolutionary rate is linked to their pathway specific position with genes directly interacting with viruses showing the highest rates of molecular evolution. We suggest that higher pathogen load in social insects indeed drives the molecular evolution of immune genes including antiviral siRNA, if not compensated by behavior.

Biogeography and evolution of seeder and resprouter forms of Erica coccinea (Ericaceae) in the fire-prone Cape fynbos

The genus Erica represents the epitome of plant biodiversity in the South African Cape fynbos with over 700 species. This genus is composed of seeder and resprouter species, but both species diversity and endemism are strongly linked to the seeder habit and concentrated in the southwestern Cape Floristic Region (CFR). Erica coccinea is a relatively abundant and widespread fynbos species whose most remarkable morphological feature is the existence of distinct seeder and resprouter forms, frequently—but not always—in disjunct populations. Both higher within-population genetic diversity and among-population differentiation have been found in seeders, most likely as a consequence of the shorter generation times and faster population turnovers. Resprouters, despite being less diverse, are suspected to be ancestral. However, no solid evidence has yet been provided for the ancestrality of the resprouter form, or for the demographic processes that have determined the current distribution of genetic diversity in both regeneration forms. Here, we used microsatellites and sequences of the nuclear ribosomal internal transcribed spacers to describe the phylogeographic structure of seeder and resprouter E. c occinea populations and provide good evidence for the ancestral status of the resprouter form and the comparatively high rates of molecular evolution in derived seeder populations. We also reveal that mixed populations, where both seeder and resprouter individuals co-occur, were originated by secondary contacts. This study highlights the role of fire in driving accelerated diversification in seeder lineages of highly speciose CFR fynbos taxa.

The uncharacterized gene 1700093K21Rik and flanking regions are correlated with reproductive isolation in the house mouse, Mus musculus

Reproductive barriers exist between the house mouse subspecies, Mus musculus musculus and M. m. domesticus, members of the Mus musculus species complex, primarily as a result of hybrid male infertility, and a hybrid zone exists where their ranges intersect in Europe. Using single nucleotide polymorphisms (SNPs) diagnostic for the two taxa, the extent of introgression across the genome was previously compared in these hybrid populations. Sixty-nine of 1316 autosomal SNPs exhibited reduced introgression in two hybrid zone transects suggesting maladaptive interactions among certain loci. One of these markers is within a region on chromosome 11 that, in other studies, has been associated with hybrid male sterility of these subspecies. We assessed sequence variation in a 20Mb region on chromosome 11 flanking this marker, and observed its inclusion within a roughly 150kb stretch of DNA showing elevated sequence differentiation between the two subspecies. Four genes are associated with this genomic subregion, with two entirely encompassed. One of the two genes, the uncharacterized 1700093K21Rik gene, displays distinguishing features consistent with a potential role in reproductive isolation between these subspecies. Along with its expression specifically within spermatogenic cells, we present various sequence analyses that demonstrate a high rate of molecular evolution of this gene, as well as identify a subspecies amino acid variant resulting in a structural difference. Taken together, the data suggest a role for this gene in reproductive isolation.

Recombination is associated with the evolution of genome structure and worker behavior in honey bees

The rise of insect societies, marked by the formation of reproductive and sterile castes, represents a major unsolved mystery in evolution. Across several independent origins of sociality, the genomes of social hymenopterans share two peculiar attributes: high recombination and low but heterogeneous GC content. For example, the genome of the honey bee, Apis mellifera, represents a mosaic of GC-poor and GC-rich regions with rates of recombination an order of magnitude higher than in humans. However, it is unclear how heterogeneity in GC content arises, and how it relates to the expression and evolution of worker traits. Using population genetic analyses, we demonstrate a bias in the allele frequency and fixation rate of derived C or G mutations in high-recombination regions, consistent with recombination's causal influence on GC-content evolution via biased gene conversion. We also show that recombination and biased gene conversion actively maintain the heterogeneous GC content of the honey bee genome despite an overall A/T mutation bias. Further, we found that GC-rich genes and intergenic regions have higher levels of genetic diversity and divergence relative to GC-poor regions, also consistent with recombination's causal influence on the rate of molecular evolution. Finally, we found that genes associated with behavior and those with worker-biased expression are found in GC-rich regions of the bee genome and also experience high rates of molecular evolution. Taken together, these findings suggest that recombination acts to maintain a genetically diverse and dynamic part of the genome where genes underlying worker behavior evolve more quickly.

Genetic Variation of Highly Pathogenic H5N1 Avian Influenza Viruses in Vietnam Shows Both Species-Specific and Spatiotemporal Associations

Domestic poultry act as a reservoir for persistent H5N1 endemicity in Vietnam, and the circulation of poultry flocks across farms and to market is thought to drive the spatial movement and evolution of avian influenza viruses. Using a dataset of complete or nearly full genomic sequences from highly pathogenic H5N1 avian influenza viruses collected in domestic poultry in Vietnam from 2003 to 2007, we explore potential differences in genetic characteristics according to species of isolation and the spatiotemporal characteristics of the viruses. Clustering algorithms and ANOVA indicate that H5N1 viruses in Vietnam show differences in the amount of genetic change that chicken viruses experience as compared to duck viruses, with duck viruses showing higher rates of molecular evolution on all eight of influenza's gene segments. There also exist distinct patterns of genetic differentiation according to the year in which they were isolated. These findings suggest that genetic evolution of avian influenza viruses is continuous through time but could also be mediated by the species in which the viruses occur, information that has implications for prevention efforts.

Comprehensive Phylogenetic Reconstruction of Amoebozoa Based on Concatenated Analyses of SSU-rDNA and Actin Genes

Evolutionary relationships within Amoebozoa have been the subject of controversy for two reasons: 1) paucity of morphological characters in traditional surveys and 2) haphazard taxonomic sampling in modern molecular reconstructions. These along with other factors have prevented the erection of a definitive system that resolves confidently both higher and lower-level relationships. Additionally, the recent recognition that many protosteloid amoebae are in fact scattered throughout the Amoebozoa suggests that phylogenetic reconstructions have been excluding an extensive and integral group of organisms. Here we provide a comprehensive phylogenetic reconstruction based on 139 taxa using molecular information from both SSU-rDNA and actin genes. We provide molecular data for 13 of those taxa, 12 of which had not been previously characterized. We explored the dataset extensively by generating 18 alternative reconstructions that assess the effect of missing data, long-branched taxa, unstable taxa, fast evolving sites and inclusion of environmental sequences. We compared reconstructions with each other as well as against previously published phylogenies. Our analyses show that many of the morphologically established lower-level relationships (defined here as relationships roughly equivalent to Order level or below) are congruent with molecular data. However, the data are insufficient to corroborate or reject the large majority of proposed higher-level relationships (above the Order-level), with the exception of Tubulinea, Archamoebae and Myxogastrea, which are consistently recovered. Moreover, contrary to previous expectations, the inclusion of available environmental sequences does not significantly improve the Amoebozoa reconstruction. This is probably because key amoebozoan taxa are not easily amplified by environmental sequencing methodology due to high rates of molecular evolution and regular occurrence of large indels and introns. Finally, in an effort to facilitate future sampling of key amoebozoan taxa, we provide a novel methodology for genome amplification and cDNA extraction from single or a few cells, a method that is culture-independent and allows both photodocumentation and extraction of multiple genes from natural samples.

Genetic diversity in birds is associated with body mass and habitat type

In some organisms, their habitat can broadly predict their population genetic diversity. For example, marine fishes harbor considerably more genetic diversity than do freshwater fishes because of the larger long-term evolutionary effective population sizes in the former. Body mass (BM) is another predictor of genetic variation, in that small-bodied mammals generally have higher rates of molecular evolution than large mammals. Does genetic variation in birds vary similarly? We investigated the relationships among microsatellite DNA diversity, BM and habitat type (aquatic or terrestrial) in 76 avian species. Our results show that across 1008 avian microsatellite loci, mean heterozygosity was positively correlated with the number of alleles per species. The mean level of heterozygosity and allele number in birds were similar to those of mammals and reptiles, but smaller than fishes. Terrestrial birds have greater genetic diversity (both in terms of mean heterozygosity and allelic diversity per population) than aquatic species. BM of aquatic birds was significantly larger than that of terrestrial birds and there was a negative relationship between mean heterozygosity and BM. Our results, interpreted in light of previously published data from other vertebrates, suggest that patterns of genetic diversity in birds depends on their evolutionary effective population size (determined in part by ecological and environmental features) and on the rate of molecular evolution.

Genome-wide investigation reveals high evolutionary rates in annual model plants

BackgroundRates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.ResultsAccording to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.ConclusionsThe higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the mechanism determining the molecular evolutionary rate at the genomic level.