Darwinian evolution in the light of genomics

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

Distinct evolutionary strategies in the GGPPS family from plants

Multiple geranylgeranyl diphosphate synthases (GGPPS) for biosynthesis of geranylgeranyl diphosphate (GGPP) exist in plants. GGPP is produced in the isoprenoid pathway and is a central precursor for various primary and specialized plant metabolites. Therefore, its biosynthesis is an essential regulatory point in the isoprenoid pathway. We selected 119 GGPPSs from 48 species representing all major plant lineages, based on stringent homology criteria. After the diversification of land plants, the number of GGPPS paralogs per species increases. Already in the moss Physcomitrella patens, GGPPS appears to be encoded by multiple paralogous genes. In gymnosperms, neofunctionalization of GGPPS may have enabled optimized biosynthesis of primary and specialized metabolites. Notably, lineage-specific expansion of GGPPS occurred in land plants. As a representative species we focused here on Arabidopsis thaliana, which retained the highest number of GGPPS paralogs (twelve) among the 48 species we considered in this study. Our results show that the A. thaliana GGPPS gene family is an example of evolution involving neo- and subfunctionalization as well as pseudogenization. We propose subfunctionalization as one of the main mechanisms allowing the maintenance of multiple GGPPS paralogs in A. thaliana genome. Accordingly, the changes in the expression patterns of the GGPPS paralogs occurring after gene duplication led to developmental and/or condition specific functional evolution.

Tinbergen's four questions: an appreciation and an update

This year is the 50th anniversary of Tinbergen's (1963) article 'On aims and methods of ethology', where he first outlined the four 'major problems of biology'. The classification of the four problems, or questions, is one of Tinbergen's most enduring legacies, and it remains as valuable today as 50 years ago in highlighting the value of a comprehensive, multifaceted understanding of a characteristic, with answers to each question providing complementary insights. Nonetheless, much has changed in the intervening years, and new data call for a more nuanced application of Tinbergen's framework. The anniversary would seem a suitable opportunity to reflect on the four questions and evaluate the scientific work that they encourage.

Horizontal gene transfer in eukaryotes: the weak-link model

The significance of horizontal gene transfer (HGT) in eukaryotic evolution remains controversial. Although many eukaryotic genes are of bacterial origin, they are often interpreted as being derived from mitochondria or plastids. Because of their fixed gene pool and gene loss, however, mitochondria and plastids alone cannot adequately explain the presence of all, or even the majority, of bacterial genes in eukaryotes. Available data indicate that no insurmountable barrier to HGT exists, even in complex multicellular eukaryotes. In addition, the discovery of both recent and ancient HGT events in all major eukaryotic groups suggests that HGT has been a regular occurrence throughout the history of eukaryotic evolution. A model of HGT is proposed that suggests both unicellular and early developmental stages as likely entry points for foreign genes into multicellular eukaryotes.

Replicon-dependent bacterial genome evolution: the case of Sinorhizobium meliloti

Many bacterial species, such as the alphaproteobacterium Sinorhizobium meliloti, are characterized by open pangenomes and contain multipartite genomes consisting of a chromosome and other large-sized replicons, such as chromids, megaplasmids, and plasmids. The evolutionary forces in both functional and structural aspects that shape the pangenome of species with multipartite genomes are still poorly understood. Therefore, we sequenced the genomes of 10 new S. meliloti strains, analyzed with four publicly available additional genomic sequences. Results indicated that the three main replicons present in these strains (a chromosome, a chromid, and a megaplasmid) partly show replicon-specific behaviors related to strain differentiation. In particular, the pSymB chromid was shown to be a hot spot for positively selected genes, and, unexpectedly, genes resident in the pSymB chromid were also found to be more widespread in distant taxa than those located in the other replicons. Moreover, through the exploitation of a DNA proximity network, a series of conserved “DNA backbones” were found to shape the evolution of the genome structure, with the rest of the genome experiencing rearrangements. The presented data allow depicting a scenario where the pSymB chromid has a distinctive role in intraspecies differentiation and in evolution through positive selection, whereas the pSymA megaplasmid mostly contributes to structural fluidity and to the emergence of new functions, indicating a specific evolutionary role for each replicon in the pangenome evolution.

What is evolutionary novelty? Process versus character based definitions

With the rise of EvoDevo, the topic of evolutionary novelty has received renewed attention. Indeed, it has been argued that one of the major contributions of EvoDevo to evolutionary theory is the explanation of phenotypic novelty. Despite such assertions, dispute continues over what exactly a novelty is and whether the term applies to a unique type of evolutionary phenomenon or whether it merely has informal meaning. In a recent special issue of J. Exp. Zool. (Mol. Dev. Evol.) dedicated to novelty, a new definition was introduced, linking novelty exclusively with adaptation and developmental constraint. In our commentary, we discuss how defining novelty in this process oriented manner leads to heightened difficulties with the application of the term and the identification of novelties. At the same time it conceals important implications for evolutionary studies. In contrast, we argue for a character based definition that is independent from adaptive necessities and promotes the integration of evolutionary factors not included in the standard theory. The implications of approaching novelty in this manner take the issue beyond definitional debates.

Zinc in evolution

The first part of this paper gives an account of the early work I did with Vallee on zinc enzymes. The second part is work and thoughts on zinc biochemistry subsequent to this work and following on from it. I shall then put our knowledge of zinc biochemistry into the context of the evolution of organisms. The zinc proteins are divided into two classes: (1) mainly enzymes often found in the early forms of life, prokaryotes, and continued in eukaryotes, showing virtually no exchange of zinc; (2) mainly buffers, pumps, chaperones and transcription factors all showing faster exchange and largely only found in eukaryotes. The complexity of organisms, especially in three steps, prokaryotes, single-cell and then multicell eukaryotes, coincides with the two major rises of oxygen. The rise of oxygen caused a rise of zinc in the sea as estimated from the analysis of sediments. These changes of zinc and other elements with oxygen indicate that environmental changes were inevitable. To a considerable extent then the rise of organism complexity arose from a need to use and in some cases to avoid the steps in element availability. After about 0.5 Ga there has been little change in the chemistry of the environment and little in the gross features of the chemistry of life notwithstanding the huge increase in diversity. The final paragraph acknowledges the inspiration and stimulation of the work of Vallee in the biology and biochemistry of zinc.

Evolution of microbes and viruses: a paradigm shift in evolutionary biology?

When Charles Darwin formulated the central principles of evolutionary biology in the Origin of Species in 1859 and the architects of the Modern Synthesis integrated these principles with population genetics almost a century later, the principal if not the sole objects of evolutionary biology were multicellular eukaryotes, primarily animals and plants. Before the advent of efficient gene sequencing, all attempts to extend evolutionary studies to bacteria have been futile. Sequencing of the rRNA genes in thousands of microbes allowed the construction of the three- domain “ribosomal Tree of Life” that was widely thought to have resolved the evolutionary relationships between the cellular life forms. However, subsequent massive sequencing of numerous, complete microbial genomes revealed novel evolutionary phenomena, the most fundamental of these being: (1) pervasive horizontal gene transfer (HGT), in large part mediated by viruses and plasmids, that shapes the genomes of archaea and bacteria and call for a radical revision (if not abandonment) of the Tree of Life concept, (2) Lamarckian-type inheritance that appears to be critical for antivirus defense and other forms of adaptation in prokaryotes, and (3) evolution of evolvability, i.e., dedicated mechanisms for evolution such as vehicles for HGT and stress-induced mutagenesis systems. In the non-cellular part of the microbial world, phylogenomics and metagenomics of viruses and related selfish genetic elements revealed enormous genetic and molecular diversity and extremely high abundance of viruses that come across as the dominant biological entities on earth. Furthermore, the perennial arms race between viruses and their hosts is one of the defining factors of evolution. Thus, microbial phylogenomics adds new dimensions to the fundamental picture of evolution even as the principle of descent with modification discovered by Darwin and the laws of population genetics remain at the core of evolutionary biology.

Rhizome of life, catastrophes, sequence exchanges, gene creations, and giant viruses: how microbial genomics challenges Darwin

Darwin's theory about the evolution of species has been the object of considerable dispute. In this review, we have described seven key principles in Darwin's book The Origin of Species and tried to present how genomics challenge each of these concepts and improve our knowledge about evolution. Darwin believed that species evolution consists on a positive directional selection ensuring the “survival of the fittest.” The most developed state of the species is characterized by increasing complexity. Darwin proposed the theory of “descent with modification” according to which all species evolve from a single common ancestor through a gradual process of small modification of their vertical inheritance. Finally, the process of evolution can be depicted in the form of a tree. However, microbial genomics showed that evolution is better described as the “biological changes over time.” The mode of change is not unidirectional and does not necessarily favors advantageous mutations to increase fitness it is rather subject to random selection as a result of catastrophic stochastic processes. Complexity is not necessarily the completion of development: several complex organisms have gone extinct and many microbes including bacteria with intracellular lifestyle have streamlined highly effective genomes. Genomes evolve through large events of gene deletions, duplications, insertions, and genomes rearrangements rather than a gradual adaptative process. Genomes are dynamic and chimeric entities with gene repertoires that result from vertical and horizontal acquisitions as well as de novo gene creation. The chimeric character of microbial genomes excludes the possibility of finding a single common ancestor for all the genes recorded currently. Genomes are collections of genes with different evolutionary histories that cannot be represented by a single tree of life (TOL). A forest, a network or a rhizome of life may be more accurate to represent evolutionary relationships among species.

Horizontal transfers of Tc1 elements between teleost fishes and their vertebrate parasites, lampreys

Horizontal gene transfer (HGT) has been recognized to be an important mechanism that shaped the evolution and genomes of prokaryotes and unicellular eukaryotes. However, HGT is regarded to be exceedingly rare among eukaryotes. We discovered massive transfers of a DNA transposon, a Tc1 element encoding a transposase, between multiple teleost fishes and lampreys that last shared a common ancestor over 500 Ma. Members of this group of Tc1 elements were found to exhibit a mosaic phylogenetic distribution, yet their sequences were highly similar even between distantly related lineages (95%-99% identity). Our molecular phylogenetic analyses suggested that horizontal transfers of this element happened repeatedly, involving multiple teleost fishes that are phylogenetically only distantly related. Interestingly, almost all the affected teleost lineages are also known to be subject to lamprey parasitism, suggesting that the horizontal transfers between vertebrates might have occurred through parasite-host interaction. The genomes of several northern hemisphere lamprey species, including that of the sea lamprey (Petromyzon marinus), were found to contain thousands of copies of the foreign elements. Impact of this event is discussed in relation to other peculiar genomic features of lampreys.

Lamarckian evolution of the giant Mimivirus in allopatric laboratory culture on amoebae

Acanthamoeba polyphaga Mimivirus has been subcultured 150 times on germ-free amoebae. This allopatric niche is very different from that found in the natural environment, where the virus is in competition with many other organisms. In this experiment, substantial gene variability and loss occurred concurrently with the emergence of phenotypically different viruses. We sought to quantify the respective roles of Lamarckian and Darwinian evolution during this experiment. We postulated that the Mimivirus genes that were down-regulated at the beginning of the allopatric laboratory culture and inactivated after 150 passages experienced Lamarckian evolution because phenotypic modifications preceded genotypic modifications, whereas we considered that genes that were highly transcribed in the new niche but were later inactivated obeyed Darwinian rules. We used the total transcript abundances and sequences described for the genes of Mimivirus at the beginning of its laboratory life and after 150 passages in allopatric culture on Acanthamoeba spp. We found a statistically significant positive correlation between the level of gene expression at the beginning of the culture and gene inactivation during the 150 passages. In particular, the mean transcript abundance at baseline was significantly lower for inactivated genes than for unchanged genes (165 ± 589 vs. 470 ± 1,625; p < 1e–3), and the mean transcript levels during the replication cycle of Mimivirus M1 were up to 8.5-fold lower for inactivated genes than for unchanged genes. In addition, proteins tended to be less frequently identified from purified virions in their early life in allopatric laboratory culture if they were encoded by variable genes than if they were encoded by conserved genes (9 vs. 15%; p = 0.062). Finally, Lamarckian evolution represented the evolutionary process encountered by 63% of the inactivated genes. Such observations may be explained by the lower level of DNA repair of useless genes.

How microbiology helps define the rhizome of life

In contrast to the tree of life (TOF) theory, species are mosaics of gene sequences with different origins. Observations of the extensive lateral sequence transfers in all organisms have demonstrated that the genomes of all life forms are collections of genes with different evolutionary histories that cannot be represented by a single TOF. Moreover, genes themselves commonly have several origins due to recombination. The human genome is not free from recombination events, so it is a mosaic like other organisms' genomes. Recent studies have demonstrated evidence for the integration of parasitic DNA into the human genome. Lateral transfer events have been accepted as major contributors of genome evolution in free-living bacteria. Furthermore, the accumulation of genomic sequence data provides evidence for extended genetic exchanges in intracellular bacteria and suggests that such events constitute an agent that promotes and maintains all bacterial species. Archaea and viruses also form chimeras containing primarily bacterial but also eukaryotic sequences. In addition to lateral transfers, orphan genes are indicative of the fact that gene creation is a permanent and unsettled phenomenon. Currently, a rhizome may more adequately represent the multiplicity and de novo creation of a genome. We wanted to confirm that the term “rhizome” in evolutionary biology applies to the entire cellular life history. This view of evolution should resemble a clump of roots representing the multiple origins of the repertoires of the genes of each species.

From environment to man: genome evolution and adaptation of human opportunistic bacterial pathogens

Environment is recognized as a huge reservoir for bacterial species and a source of human pathogens. Some environmental bacteria have an extraordinary range of activities that include promotion of plant growth or disease, breakdown of pollutants, production of original biomolecules, but also multidrug resistance and human pathogenicity. The versatility of bacterial life-style involves adaptation to various niches. Adaptation to both open environment and human specific niches is a major challenge that involves intermediate organisms allowing pre-adaptation to humans. The aim of this review is to analyze genomic features of environmental bacteria in order to explain their adaptation to human beings. The genera Pseudomonas, Aeromonas and Ochrobactrum provide valuable examples of opportunistic behavior associated to particular genomic structure and evolution. Particularly, we performed original genomic comparisons among aeromonads and between the strictly intracellular pathogens Brucella spp. and the mild opportunistic pathogens Ochrobactrum spp. We conclude that the adaptation to human could coincide with a speciation in action revealed by modifications in both genomic and population structures. This adaptation-driven speciation could be a major mechanism for the emergence of true pathogens besides the acquisition of specialized virulence factors.

An evolutionary network of genes present in the eukaryote common ancestor polls genomes on eukaryotic and mitochondrial origin

To test the predictions of competing and mutually exclusive hypotheses for the origin of eukaryotes, we identified from a sample of 27 sequenced eukaryotic and 994 sequenced prokaryotic genomes 571 genes that were present in the eukaryote common ancestor and that have homologues among eubacterial and archaebacterial genomes. Maximum-likelihood trees identified the prokaryotic genomes that most frequently contained genes branching as the sister to the eukaryotic nuclear homologues. Among the archaebacteria, euryarchaeote genomes most frequently harbored the sister to the eukaryotic nuclear gene, whereas among eubacteria, the α-proteobacteria were most frequently represented within the sister group. Only 3 genes out of 571 gave a 3-domain tree. Homologues from α-proteobacterial genomes that branched as the sister to nuclear genes were found more frequently in genomes of facultatively anaerobic members of the rhiozobiales and rhodospirilliales than in obligate intracellular ricketttsial parasites. Following α-proteobacteria, the most frequent eubacterial sister lineages were γ-proteobacteria, δ-proteobacteria, and firmicutes, which were also the prokaryote genomes least frequently found as monophyletic groups in our trees. Although all 22 higher prokaryotic taxa sampled (crenarchaeotes, γ-proteobacteria, spirochaetes, chlamydias, etc.) harbor genes that branch as the sister to homologues present in the eukaryotic common ancestor, that is not evidence of 22 different prokaryotic cells participating at eukaryote origins because prokaryotic “lineages” have laterally acquired genes for more than 1.5 billion years since eukaryote origins. The data underscore the archaebacterial (host) nature of the eukaryotic informational genes and the eubacterial (mitochondrial) nature of eukaryotic energy metabolism. The network linking genes of the eukaryote ancestor to contemporary homologues distributed across prokaryotic genomes elucidates eukaryote gene origins in a dialect cognizant of gene transfer in nature.

General trends in the utilization of structural factors contributing to biological complexity

During evolution, proteins containing newly emerged domains and the increasing proportion of multidomain proteins in the full genome-encoded proteome (GEP) have substantially contributed to increasing biological complexity. However, it is not known how these two potential structural factors are preferentially utilized at given physiological states. Here, we classified proteins according to domain number and domain age and explored the general trends across species for the utilization of proteins from GEP to various certain-state proteomes (CSPs, i.e., all the proteins expressed at certain physiological states). We found that multidomain proteins or only older domain-containing proteins are significantly overrepresented in CSPs compared with GEP, which is a trend that is stronger in multicellular organisms than in unicellular organisms. Interestingly, the strengths of overrepresentation decreased during evolution of multicellular eukaryotes. When comparing across CSPs, we found that multidomain proteins are more overrepresented in complex tissues than in simpler ones, whereas no difference among proteins with domains of different ages is evident between complex and simple tissues. Thus, biological complexity under certain conditions is more significantly realized by diverse domain organization than by the emergence of new types of domain. In addition, we found that multidomain or only older domain-containing proteins tend to evolve slowly and generally are under stronger purifying selection, which may partly result from their general overrepresentation trends in CSPs.

Evolutionary systems biology: historical and philosophical perspectives on an emerging synthesis

Systems biology (SB) is at least a decade old now and maturing rapidly. A more recent field, evolutionary systems biology (ESB), is in the process of further developing system-level approaches through the expansion of their explanatory and potentially predictive scope. This chapter will outline the varieties of ESB existing today by tracing the diverse roots and fusions that make up this integrative project. My approach is philosophical and historical. As well as examining the recent origins of ESB, I will reflect on its central features and the different clusters of research it comprises. In its broadest interpretation, ESB consists of five overlapping approaches: comparative and correlational ESB; network architecture ESB; network property ESB; population genetics ESB; and finally, standard evolutionary questions answered with SB methods. After outlining each approach with examples, I will examine some strong general claims about ESB, particularly that it can be viewed as the next step toward a fuller modern synthesis of evolutionary biology (EB), and that it is also the way forward for evolutionary and systems medicine. I will conclude with a discussion of whether the emerging field of ESB has the capacity to combine an even broader scope of research aims and efforts than it presently does.

Rho GTPases: deciphering the evolutionary history of a complex protein family

Rho GTPases constitute a significant subgroup of the eukaryotic Ras superfamily of small GTPases implicated in the regulation of diverse cellular processes, such as the dynamics of the actin cytoskeleton, establishment, and maintenance of cell polarity and membrane trafficking. Whereas a few eukaryotes lack Rho genes, a majority of species typically bear multiple Rho paralogs, raising a question about the origin of the family and the paths of its diversification in individual eukaryotic lineages. In this chapter, we ruminate on several aspects of the evolutionary history of the Rho family and methodological challenges of its reconstruction. First, we provide an updated survey of Rho GTPases in diverse eukaryotic branches, demonstrating almost ubiquitous occurrence of Rho genes across the eukaryotic phylogeny most consistent with the presence of at least one Rho gene already in the last eukaryotic common ancestor. Second, we discuss the obstacles in reconstructing the history of gene duplications giving rise to the extant diversity of Rho paralogs in different species, and point to numerous limitations posed by the current phylogenetic methodology. Third, as a case study demonstrating various issues of data collection, phylogenetic analyses and interpretations of trees, we present an analysis of the Rho family in the fungal kingdom, revealing the existence of at least four separate paralogs (Cdc42, Rac, Rho1, and Rho4) in early fungi and subsequent potentially independent expansions of the family in different fungal subgroups. We conclude with the warning that the currently dominating perception of the Rho phylogeny is biased by the metazoan (and especially vertebrate) perspective, and a new, more global view is to be worked out when a better genome sampling and more adequate methods of phylogenetic inference are employed.

Comparative interaction networks: bridging genotype to phenotype

Over the past decade, biomedical research has witnessed an exponential increase in the throughput of the characterization of biological systems. Here we review the recent progress in large-scale methods to determine protein-protein, genetic and chemical-genetic interaction networks. We discuss some of the limitations and advantages of the different methods and give examples of how these networks are being used to study the evolutionary process. Comparative studies have revealed that different types of protein-protein interactions diverge at different rates with high conservation of co-complex membership but rapid divergence of more promiscuous interactions like those that mediate post-translational modifications. These evolutionary trends have consistent genetic consequences with highly conserved epistatic interactions within complex subunits but faster divergence of epistatic interactions across complexes or pathways. Finally, we discuss how these evolutionary observations are being used to interpret cross-species chemical-genetic studies and how they might shape therapeutic strategies. Together, these interaction networks offer us an unprecedented level of detail into how genotypes are translated to phenotypes, and we envision that they will be increasingly useful in the interpretation of genetic and phenotypic variation occurring within populations as well as the rational design of combinatorial therapeutics.

Neuropsychiatric Genetics of Happiness, Friendships, and Politics: Hypothesizing Homophily ("Birds of a Feather Flock Together") as a Function of Reward Gene Polymorphisms

Mindful of the new evolutionary ideas related to an emerging scientific focus known as omics, we propose that spiritual, social, and political behaviors may be tied in part to inheritable reward gene polymorphisms, as has been demonstrated for the addictions. If so, analyses of gene polymorphisms may assist in predicting liberalism or conservatism in partisan attachments. For example, both drinking (alcohol) and obesity seem to cluster in large social networks and are influenced by friends having the same genotype, in particular the DRD2 A1 allele. Likewise, voting, voting turnout and attachment to a particular political ideology is differentially related to various reward genes (e.g., 5HTT, MOA, DRD2, and DRD4), possibly predicting liberalism or conservatism. Moreover, voters' genetic information may predict presidential outcomes more than the actual issues at hand or the presidential candidates themselves. Thus, political discussions on TV, radio, or other media may be morphed by one's reward gene polymorphisms and as such, may explain the prevalence of generations of die-hard republicans and equally entrenched democratic legacies. Indeed, even in politics, birds of a feather (homophily) flock together. We caution that our proposal should be viewed mindfully awaiting additional research before definitive statements or conclusions can be derived from the studies to date, and we encourage large scale studies to confirm these earlier reports.

Mitonuclear match: optimizing fitness and fertility over generations drives ageing within generations

Many conserved eukaryotic traits, including apoptosis, two sexes, speciation and ageing, can be causally linked to a bioenergetic requirement for mitochondrial genes. Mitochondrial genes encode proteins involved in cell respiration, which interact closely with proteins encoded by nuclear genes. Functional respiration requires the coadaptation of mitochondrial and nuclear genes, despite divergent tempi and modes of evolution. Free-radical signals emerge directly from the biophysics of mosaic respiratory chains encoded by two genomes prone to mismatch, with apoptosis being the default penalty for compromised respiration. Selection for genomic matching is facilitated by two sexes, and optimizes fitness, adaptability and fertility in youth. Mismatches cause infertility, low fitness, hybrid breakdown, and potentially speciation. The dynamics of selection for mitonuclear function optimize fitness over generations, but the same selective processes also operate within generations, driving ageing and age-related diseases. This coherent view of eukaryotic energetics offers striking insights into infertility and age-related diseases.

Protein-based signatures of functional evolution in Plasmodium falciparum

Background

It has been known for over a decade that Plasmodium falciparum proteins are enriched in non-globular domains of unknown function. The potential for these regions of protein sequence to undergo high levels of genetic drift provides a fundamental challenge to attempts to identify the molecular basis of adaptive change in malaria parasites.

Results

Evolutionary comparisons were undertaken using a set of forty P. falciparum metabolic enzyme genes, both within the hominid malaria clade (P. reichenowi) and across the genus (P. chabaudi). All genes contained coding elements highly conserved across the genus, but there were also a large number of regions of weakly or non-aligning coding sequence. These displayed remarkable levels of non-synonymous fixed differences within the hominid malaria clade indicating near complete release from purifying selection (dN/dS ratio at residues non-aligning across genus: 0.64, dN/dS ratio at residues identical across genus: 0.03). Regions of low conservation also possessed high levels of hydrophilicity, a marker of non-globularity. The propensity for such regions to act as potent sources of non-synonymous genetic drift within extant P. falciparum isolates was confirmed at chromosomal regions containing genes known to mediate drug resistance in field isolates, where 150 of 153 amino acid variants were located in poorly conserved regions. In contrast, all 22 amino acid variants associated with drug resistance were restricted to highly conserved regions. Additional mutations associated with laboratory-selected drug resistance, such as those in PfATPase4 selected by spiroindolone, were similarly restricted while mutations in another calcium ATPase (PfSERCA, a gene proposed to mediate artemisinin resistance) that reach significant frequencies in field isolates were located exclusively in poorly conserved regions consistent with genetic drift.

Conclusion

Coding sequences of malaria parasites contain prospectively definable domains subject to neutral or nearly neutral evolution on a scale that appears unrivalled in biology. This distinct evolutionary landscape has potential to confound analytical methods developed for other genera. Against this tide of genetic drift, polymorphisms mediating functional change stand out to such an extent that evolutionary context provides a useful signal for identifying the molecular basis of drug resistance in malaria parasites, a finding that is of relevance to both genome-wide and candidate gene studies in this genus.

Energetics and genetics across the prokaryote-eukaryote divide

BACKGROUND

All complex life on Earth is eukaryotic. All eukaryotic cells share a common ancestor that arose just once in four billion years of evolution. Prokaryotes show no tendency to evolve greater morphological complexity, despite their metabolic virtuosity. Here I argue that the eukaryotic cell originated in a unique prokaryotic endosymbiosis, a singular event that transformed the selection pressures acting on both host and endosymbiont.

RESULTS

The reductive evolution and specialisation of endosymbionts to mitochondria resulted in an extreme genomic asymmetry, in which the residual mitochondrial genomes enabled the expansion of bioenergetic membranes over several orders of magnitude, overcoming the energetic constraints on prokaryotic genome size, and permitting the host cell genome to expand (in principle) over 200,000-fold. This energetic transformation was permissive, not prescriptive; I suggest that the actual increase in early eukaryotic genome size was driven by a heavy early bombardment of genes and introns from the endosymbiont to the host cell, producing a high mutation rate. Unlike prokaryotes, with lower mutation rates and heavy selection pressure to lose genes, early eukaryotes without genome-size limitations could mask mutations by cell fusion and genome duplication, as in allopolyploidy, giving rise to a proto-sexual cell cycle. The side effect was that a large number of shared eukaryotic basal traits accumulated in the same population, a sexual eukaryotic common ancestor, radically different to any known prokaryote.

CONCLUSIONS

The combination of massive bioenergetic expansion, release from genome-size constraints, and high mutation rate favoured a protosexual cell cycle and the accumulation of eukaryotic traits. These factors explain the unique origin of eukaryotes, the absence of true evolutionary intermediates, and the evolution of sex in eukaryotes but not prokaryotes.

REVIEWERS

This article was reviewed by: Eugene Koonin, William Martin, Ford Doolittle and Mark van der Giezen. For complete reports see the Reviewers' Comments section.

How stands the Tree of Life a century and a half after The Origin?

We examine the Tree of Life (TOL) as an evolutionary hypothesis and a heuristic. The original TOL hypothesis has failed but a new "statistical TOL hypothesis" is promising. The TOL heuristic usefully organizes data without positing fundamental evolutionary truth.

Iron in evolution

Iron chemistry in the environment and in organisms is entwined. The iron surface minerals in solution for the first billion years of the planet were ferrous compounds. This ion became and has remained a major participant in organisms. The evolution of iron was due to its oxidation to insoluble ferric ions by oxygen released from organisms. The evolution of cellular iron chemistry then required uptake from this oxidised state. Use was expanded from the mainly electron transfer properties in the original reductive cell interior to employment in external oxidative chemistry. The environment/organisms evolution is that of one predictable chemical system.

Positional information resolves structural variations and uncovers an evolutionarily divergent genetic locus in accessions of Arabidopsis thaliana

Genome sequencing of closely related individuals has yielded valuable insights that link genome evolution to phenotypic variations. However, advancement in sequencing technology has also led to an escalation in the number of poor quality–drafted genomes assembled based on reference genomes that can have highly divergent or haplotypic regions. The self-fertilizing nature of Arabidopsis thaliana poses an advantage to sequencing projects because its genome is mostly homozygous. To determine the accuracy of an Arabidopsis drafted genome in less conserved regions, we performed a resequencing experiment on a ∼371-kb genomic interval in the Landsberg erecta (Ler-0) accession. We identified novel structural variations (SVs) between Ler-0 and the reference accession Col-0 using a long-range polymerase chain reaction approach to generate an Illumina data set that has positional information, that is, a data set with reads that map to a known location. Positional information is important for accurate genome assembly and the resolution of SVs particularly in highly duplicated or repetitive regions. Sixty-one regions with misassembly signatures were identified from the Ler-0 draft, suggesting the presence of novel SVs that are not represented in the draft sequence. Sixty of those were resolved by iterative mapping using our data set. Fifteen large indels (>100 bp) identified from this study were found to be located either within protein-coding regions or upstream regulatory regions, suggesting the formation of novel alleles or altered regulation of existing genes in Ler-0. We propose future genome-sequencing experiments to follow a clone-based approach that incorporates positional information to ultimately reveal haplotype-specific differences between accessions.

Exploring the symbiotic pangenome of the nitrogen-fixing bacterium Sinorhizobium meliloti

Background

Sinorhizobium meliloti is a model system for the studies of symbiotic nitrogen fixation. An extensive polymorphism at the genetic and phenotypic level is present in natural populations of this species, especially in relation with symbiotic promotion of plant growth. AK83 and BL225C are two nodule-isolated strains with diverse symbiotic phenotypes; BL225C is more efficient in promoting growth of the Medicago sativa plants than strain AK83. In order to investigate the genetic determinants of the phenotypic diversification of S. meliloti strains AK83 and BL225C, we sequenced the complete genomes for these two strains.

Results

With sizes of 7.14 Mbp and 6.97 Mbp, respectively, the genomes of AK83 and BL225C are larger than the laboratory strain Rm1021. The core genome of Rm1021, AK83, BL225C strains included 5124 orthologous groups, while the accessory genome was composed by 2700 orthologous groups. While Rm1021 and BL225C have only three replicons (Chromosome, pSymA and pSymB), AK83 has also two plasmids, 260 and 70 Kbp long. We found 65 interesting orthologous groups of genes that were present only in the accessory genome, consequently responsible for phenotypic diversity and putatively involved in plant-bacterium interaction. Notably, the symbiosis inefficient AK83 lacked several genes required for microaerophilic growth inside nodules, while several genes for accessory functions related to competition, plant invasion and bacteroid tropism were identified only in AK83 and BL225C strains. Presence and extent of polymorphism in regulons of transcription factors involved in symbiotic interaction were also analyzed. Our results indicate that regulons are flexible, with a large number of accessory genes, suggesting that regulons polymorphism could also be a key determinant in the variability of symbiotic performances among the analyzed strains.

Conclusions

In conclusions, the extended comparative genomics approach revealed a variable subset of genes and regulons that may contribute to the symbiotic diversity.

Gene fragmentation: a key to mitochondrial genome evolution in Euglenozoa?

Phylum Euglenozoa comprises three groups of eukaryotic microbes (kinetoplastids, diplonemids, and euglenids), the mitochondrial (mt) genomes of which exhibit radically different modes of organization and expression. Gene fragmentation is a striking feature of both euglenid and diplonemid mtDNAs. To rationalize the emergence of these highly divergent mtDNA types and the existence of insertion/deletion RNA editing (in kinetoplastids) and trans-splicing (in diplonemids), we propose that in the mitochondrion of the common evolutionary ancestor of Euglenozoa, small expressed gene fragments promoted a rampant neutral evolutionary pathway. Interactions between small antisense transcripts of these gene fragments and full-length transcripts, assisted by RNA-processing enzymes, permitted the emergence of RNA editing and/or trans-splicing activities, allowing the system to tolerate indel mutations and further gene fragmentation, respectively, and leading to accumulation of additional mutations. In this way, dramatically different mitochondrial genome structures and RNA-processing machineries were able to evolve. The paradigm of constructive neutral evolution acting on the widely different mitochondrial genetic systems in Euglenozoa posits the accretion of initially neutral molecular interactions by genetic drift, leading inevitably to the observed 'irremediable complexity'.

Constraints and plasticity in genome and molecular-phenome evolution

Multiple constraints variously affect different parts of the genomes of diverse life forms. The selective pressures that shape the evolution of viral, archaeal, bacterial and eukaryotic genomes differ markedly, even among relatively closely related animal and bacterial lineages; by contrast, constraints affecting protein evolution seem to be more universal. The constraints that shape the evolution of genomes and phenomes are complemented by the plasticity and robustness of genome architecture, expression and regulation. Taken together, these findings are starting to reveal complex networks of evolutionary processes that must be integrated to attain a new synthesis of evolutionary biology.

The effects of low-impact mutations in digital organisms

Background

Avida is a computer program that performs evolution experiments with digital organisms. Previous work has used the program to study the evolutionary origin of complex features, namely logic operations, but has consistently used extremely large mutational fitness effects. The present study uses Avida to better understand the role of low-impact mutations in evolution.

Results

When mutational fitness effects were approximately 0.075 or less, no new logic operations evolved, and those that had previously evolved were lost. When fitness effects were approximately 0.2, only half of the operations evolved, reflecting a threshold for selection breakdown. In contrast, when Avida's default fitness effects were used, all operations routinely evolved to high frequencies and fitness increased by an average of 20 million in only 10,000 generations.

Conclusions

Avidian organisms evolve new logic operations only when mutations producing them are assigned high-impact fitness effects. Furthermore, purifying selection cannot protect operations with low-impact benefits from mutational deterioration. These results suggest that selection breaks down for low-impact mutations below a certain fitness effect, the selection threshold. Experiments using biologically relevant parameter settings show the tendency for increasing genetic load to lead to loss of biological functionality. An understanding of such genetic deterioration is relevant to human disease, and may be applicable to the control of pathogens by use of lethal mutagenesis.

A comprehensive benchmark study of multiple sequence alignment methods: current challenges and future perspectives

Multiple comparison or alignmentof protein sequences has become a fundamental tool in many different domains in modern molecular biology, from evolutionary studies to prediction of 2D/3D structure, molecular function and inter-molecular interactions etc. By placing the sequence in the framework of the overall family, multiple alignments can be used to identify conserved features and to highlight differences or specificities. In this paper, we describe a comprehensive evaluation of many of the most popular methods for multiple sequence alignment (MSA), based on a new benchmark test set. The benchmark is designed to represent typical problems encountered when aligning the large protein sequence sets that result from today's high throughput biotechnologies. We show that alignmentmethods have significantly progressed and can now identify most of the shared sequence features that determine the broad molecular function(s) of a protein family, even for divergent sequences. However,we have identified a number of important challenges. First, the locally conserved regions, that reflect functional specificities or that modulate a protein's function in a given cellular context,are less well aligned. Second, motifs in natively disordered regions are often misaligned. Third, the badly predicted or fragmentary protein sequences, which make up a large proportion of today's databases, lead to a significant number of alignment errors. Based on this study, we demonstrate that the existing MSA methods can be exploited in combination to improve alignment accuracy, although novel approaches will still be needed to fully explore the most difficult regions. We then propose knowledge-enabled, dynamic solutions that will hopefully pave the way to enhanced alignment construction and exploitation in future evolutionary systems biology studies.

Genomes of the most dangerous epidemic bacteria have a virulence repertoire characterized by fewer genes but more toxin-antitoxin modules

BACKGROUND

We conducted a comparative genomic study based on a neutral approach to identify genome specificities associated with the virulence capacity of pathogenic bacteria. We also determined whether virulence is dictated by rules, or if it is the result of individual evolutionary histories. We systematically compared the genomes of the 12 most dangerous pandemic bacteria for humans ("bad bugs") to their closest non-epidemic related species ("controls").

METHODOLOGY/PRINCIPAL FINDINGS

We found several significantly different features in the "bad bugs", one of which was a smaller genome that likely resulted from a degraded recombination and repair system. The 10 Cluster of Orthologous Group (COG) functional categories revealed a significantly smaller number of genes in the "bad bugs", which lacked mostly transcription, signal transduction mechanisms, cell motility, energy production and conversion, and metabolic and regulatory functions. A few genes were identified as virulence factors, including secretion system proteins. Five "bad bugs" showed a greater number of poly (A) tails compared to the controls, whereas an elevated number of poly (A) tails was found to be strongly correlated to a low GC% content. The "bad bugs" had fewer tandem repeat sequences compared to controls. Moreover, the results obtained from a principal component analysis (PCA) showed that the "bad bugs" had surprisingly more toxin-antitoxin modules than did the controls.

CONCLUSIONS/SIGNIFICANCE

We conclude that pathogenic capacity is not the result of "virulence factors" but is the outcome of a virulent gene repertoire resulting from reduced genome repertoires. Toxin-antitoxin systems could participate in the virulence repertoire, but they may have developed independently of selfish evolution.

The struggle for life of the genome's selfish architects

Transposable elements (TEs) were first discovered more than 50 years ago, but were totally ignored for a long time. Over the last few decades they have gradually attracted increasing interest from research scientists. Initially they were viewed as totally marginal and anecdotic, but TEs have been revealed as potentially harmful parasitic entities, ubiquitous in genomes, and finally as unavoidable actors in the diversity, structure, and evolution of the genome. Since Darwin's theory of evolution, and the progress of molecular biology, transposable elements may be the discovery that has most influenced our vision of (genome) evolution. In this review, we provide a synopsis of what is known about the complex interactions that exist between transposable elements and the host genome. Numerous examples of these interactions are provided, first from the standpoint of the genome, and then from that of the transposable elements. We also explore the evolutionary aspects of TEs in the light of post-Darwinian theories of evolution.

High sensitivity to aligner and high rate of false positives in the estimates of positive selection in the 12 Drosophila genomes

We investigate the effect of aligner choice on inferences of positive selection using site-specific models of molecular evolution. We find that independently of the choice of aligner, the rate of false positives is unacceptably high. Our study is a whole-genome analysis of all protein-coding genes in 12 Drosophila genomes annotated in either all 12 species (~6690 genes) or in the six melanogaster group species. We compare six popular aligners: PRANK, T-Coffee, ClustalW, ProbCons, AMAP, and MUSCLE, and find that the aligner choice strongly influences the estimates of positive selection. Differences persist when we use (1) different stringency cutoffs, (2) different selection inference models, (3) alignments with or without gaps, and/or additional masking, (4) per-site versus per-gene statistics, (5) closely related melanogaster group species versus more distant 12 Drosophila genomes. Furthermore, we find that these differences are consequential for downstream analyses such as determination of over/under-represented GO terms associated with positive selection. Visual analysis indicates that most sites inferred as positively selected are, in fact, misaligned at the codon level, resulting in false positive rates of 48%-82%. PRANK, which has been reported to outperform other aligners in simulations, performed best in our empirical study as well. Unfortunately, PRANK still had a high, and unacceptable for most applications, false positives rate of 50%-55%. We identify misannotations and indels, many of which appear to be located in disordered protein regions, as primary culprits for the high misalignment-related error levels and discuss possible workaround approaches to this apparently pervasive problem in genome-wide evolutionary analyses.

Nonadaptive processes in primate and human evolution

Evolutionary biology has tended to focus on adaptive evolution by positive selection as the primum mobile of evolutionary trajectories in species while underestimating the importance of nonadaptive evolutionary processes. In this review, I describe evidence that suggests that primate and human evolution has been strongly influenced by nonadaptive processes, particularly random genetic drift and mutation. This is evidenced by three fundamental effects: a relative relaxation of selective constraints (i.e., purifying selection), a relative increase in the fixation of slightly deleterious mutations, and a general reduction in the efficacy of positive selection. These effects are observed in protein-coding, regulatory regions, and in gene expression data, as well as in an augmentation of fixation of large-scale mutations, including duplicated genes, mobile genetic elements, and nuclear mitochondrial DNA. The evidence suggests a general population-level explanation such as a reduction in effective population size (N(e)). This would have tipped the balance between the evolutionary forces of natural selection and random genetic drift toward genetic drift for variants having small selective effects. After describing these proximate effects, I describe the potential consequences of these effects for primate and human evolution. For example, an increase in the fixation of slightly deleterious mutations could potentially have led to an increase in the fixation rate of compensatory mutations that act to suppress the effects of slightly deleterious substitutions. The potential consequences of compensatory evolution for the evolution of novel gene functions and in potentially confounding the detection of positively selected genes are explored. The consequences of the passive accumulation of large-scale genomic mutations by genetic drift are unclear, though evidence suggests that new gene copies as well as insertions of transposable elements into genes can potentially lead to adaptive phenotypes. Finally, because a decrease in selective constraint at the genetic level is expected to have effects at the morphological level, I review studies that compare rates of morphological change in various mammalian and island populations where N(e) is reduced. Furthermore, I discuss evidence that suggests that craniofacial morphology in the Homo lineage has shifted from an evolutionary rate constrained by purifying selection toward a neutral evolutionary rate.

Defining pathogenic bacterial species in the genomic era

Actual definitions of bacterial species are limited due to the current criteria of definition and the use of restrictive genetic tools. The 16S ribosomal RNA sequence, for example, has been widely used as a marker for phylogenetic analyses; however, its use often leads to misleading species definitions. According to the first genetic studies, removing a certain number of genes from pathogenic bacteria removes their capacity to infect hosts. However, more recent studies have demonstrated that the specialization of bacteria in eukaryotic cells is associated with massive gene loss, especially for allopatric endosymbionts that have been isolated for a long time in an intracellular niche. Indeed, sympatric free-living bacteria often have bigger genomes and exhibit greater resistance and plasticity and constitute species complexes rather than true species. Specialists, such as pathogenic bacteria, escape these bacterial complexes and colonize a niche, thereby gaining a species name. Their specialization allows them to become allopatric, and their gene losses eventually favor reductive genome evolution. A pathogenic species is characterized by a gene repertoire that is defined not only by genes that are present but also by those that are lacking. It is likely that current bacterial pathogens will disappear soon and be replaced by new ones that will emerge from bacterial complexes that are already in contact with humans.

Analysis of the 3' ends of tRNA as the cause of insertion sites of foreign DNA in Prochlorococcus

The purpose of this study was to investigate the characteristics of transfer RNA (tRNA) responsible for the association between tRNA genes and genes of apparently foreign origin (genomic islands) in five high-light adapted Prochlorococcus strains. Both bidirectional best BLASTP (basic local alignment search tool for proteins) search and the conservation of gene order against each other were utilized to identify genomic islands, and 7 genomic islands were found to be immediately adjacent to tRNAs in Prochlorococcus marinus AS9601, 11 in P. marinus MIT9515, 8 in P. marinus MED4, 6 in P. marinus MIT9301, and 6 in P. marinus MIT9312. Monte Carlo simulation showed that tRNA genes are hotspots for the integration of genomic islands in Prochlorococcus strains. The tRNA genes associated with genomic islands showed the following characteristics: (1) the association was biased towards a specific subset of all iso-accepting tRNA genes; (2) the codon usages of genes within genomic islands appear to be unrelated to the codons recognized by associated tRNAs; and, (3) the majority of the 3' ends of associated tRNAs lack CCA ends. These findings contradict previous hypotheses concerning the molecular basis for the frequent use of tRNA as the insertion site for foreign genetic materials. The analysis of a genomic island associated with a tRNA-Asn gene in P. marinus MIT9301 suggests that foreign genetic material is inserted into the host genomes by means of site-specific recombination, with the 3' end of the tRNA as the target, and during the process, a direct repeat of the 3' end sequence of a boundary tRNA (namely, a scar from the process of insertion) is formed elsewhere in the genomic island. Through the analysis of the sequences of these targets, it can be concluded that a region characterized by both high GC content and a palindromic structure is the preferred insertion site.

Horizontal gene transfers with or without cell fusions in all categories of the living matter

This article reviews the history of widespread exchanges of genetic segments initiated over 3 billion years ago, to be part of their life style, by sphero-protoplastic cells, the ancestors of archaea, prokaryota, and eukaryota. These primordial cells shared a hostile anaerobic and overheated environment and competed for survival. "Coexist with, or subdue and conquer, expropriate its most useful possessions, or symbiose with it, your competitor" remain cellular life's basic rules. This author emphasizes the role of viruses, both in mediating cell fusions, such as the formation of the first eukaryotic cell(s) from a united crenarchaeon and prokaryota, and the transfer of host cell genes integrated into viral (phages) genomes. After rising above the Darwinian threshold, rigid rules of speciation and vertical inheritance in the three domains of life were established, but horizontal gene transfers with or without cell fusions were never abolished. The author proves with extensive, yet highly selective documentation, that not only unicellular microorganisms, but the most complex multicellular entities of the highest ranks resort to, and practice, cell fusions, and donate and accept horizontally (laterally) transferred genes. Cell fusions and horizontally exchanged genetic materials remain the fundamental attributes and inherent characteristics of the living matter, whether occurring accidentally or sought after intentionally. These events occur to cells stagnating for some 3 milliard years at a lower yet amazingly sophisticated level of evolution, and to cells achieving the highest degree of differentiation, and thus functioning in dependence on the support of a most advanced multicellular host, like those of the human brain. No living cell is completely exempt from gene drains or gene insertions.

Network biology: a direct approach to study biological function

In this paper we discuss the dualism of gene networks and their role in systems biology. We argue that gene networks (1) can serve as a conceptual framework, forming a fundamental level of a phenomenological description, and (2) are a means to represent and analyze data. The latter point does not only allow a systems analysis but is even amenable for a direct approach to study biological function. Here we focus on the clarity of our main arguments and conceptual meaning of gene networks, rather than the causal inference of gene networks from data. WIREs Syst Biol Med 2011 3 379-391 DOI: 10.1002/wsbm.134 For further resources related to this article, please visit the WIREs website.

The energetics of genome complexity

All complex life is composed of eukaryotic (nucleated) cells. The eukaryotic cell arose from prokaryotes just once in four billion years, and otherwise prokaryotes show no tendency to evolve greater complexity. Why not? Prokaryotic genome size is constrained by bioenergetics. The endosymbiosis that gave rise to mitochondria restructured the distribution of DNA in relation to bioenergetic membranes, permitting a remarkable 200,000-fold expansion in the number of genes expressed. This vast leap in genomic capacity was strictly dependent on mitochondrial power, and prerequisite to eukaryote complexity: the key innovation en route to multicellular life.