Comparison of the genome for phylogenetically related bacteriophages phiKZ and EL of Pseudomonas aeruginosa: evolutionary aspects and minimal genome size

Bacteriophages of the family Myoviridae represent one of the most widespread domains of the biosphere substantially affecting the ecological balance of microorganisms. Interestingly, sequence analysis of genomic DNAs of large bacteriophages revealed many genes coding for proteins with unknown functions. A new approach is proposed to improve the functional identification of genes. This approach is based on comparing the genome sequence for phylogenetically and morphologically related phages showing no considerable homology at the level of genomic DNA. It is assumed that gene functions essential for the development of phages of a given family are conserved and that the corresponding genes code for similar orthologous proteins even when lacking sequence homology. The genome was sequenced and compared for two Pseudomonas aeruginosa giant bacteriophages, phiKZ and EL, which belong to a group of (phiKZ-related phages. A substantial difference in genome organization was observed, suggesting specific features of phage evolution. In addition, the problem of the minimal genome of the superfamily is discussed on the basis of the difference in size and structure between the phiKZ and EL genomes.     

Evolution of living organisms as a multilevel process

There is inherent capacity to increase the degree of aggregation within each of the levels of structural organization of living matter. At the macromolecular level (MML), this is an increase in the gene number in the genomes of evolving organisms; at the cellular level (CL), an increase in cell size; and at the multicellular level (MCL), an increase in the number of cells in the multicellular aggregate. However, the increase in the degree of aggregation causes gene incompatibility in case of genome evolution and instability in case of large cells and multicellular aggregates with simple structure. Gene incompatibility may be neutralized by spacio-temporal disconnection of the products of incompatible genes at the cellular and multicellular levels. The larger cells and multicellular aggregates are stabilized by increased structural complexity which is a consequence of the origin of new genes. There is a feedback between the processes of evolution at different levels MML→CL→ MCL.The processes of evolutionary development at different levels of structural organization are also relatively independent. The coincidence of these processes gives rise to stable organisms of higher complexity, which are then subjected to natural selection and population processes to establish a new step in progressive biological evolution. In all of the normal organisms of newly evolved species there is a correspondence between the different levels of structural organization, i.e. in their degree of aggregation, their complexity and functional organization. The form of correspondence for multicellular organisms is presented.     

Evolution of genome size: multilevel selection, mutation bias or dynamical chaos?

In the past two years, new data on conceptual aspects of the evolution of eukaryotic genome size have appeared, including the adaptivity of genome enlargement, the mechanisms of genome size change and the relation of genome size to organismal complexity. New data on the hypotheses of "selfish DNA" and "mutational equilibrium" have been recently obtained. A relationship is emerging between the intragenomic distribution of noncoding DNA and differential gene expression, which suggests that noncoding DNA is involved in epigenetic organization of the genome and organismal complexity. The standpoint of dynamical chaos, which integrates multilevel selection and mutation biases, may provide a framework for studying the evolution of genome size.     

From 2R to 3R: evidence for a fish-specific genome duplication (FSGD)

An important mechanism for the evolution of phenotypic complexity, diversity and innovation, and the origin of novel gene functions is the duplication of genes and entire genomes. Recent phylogenomic studies suggest that, during the evolution of vertebrates, the entire genome was duplicated in two rounds (2R) of duplication. Later, approximately 350 mya, in the stem lineage of ray-finned (actinopterygian) fishes, but not in that of the land vertebrates, a third genome duplication occurred-the fish-specific genome duplication (FSGD or 3R), leading, at least initially, to up to eight copies of the ancestral deuterostome genome. Therefore, the sarcopterygian (lobe-finned fishes and tetrapods) genome possessed originally only half as many genes compared to the derived fishes, just like the most-basal and species-poor lineages of extant fishes that diverged from the fish stem lineage before the 3R duplication. Most duplicated genes were secondarily lost, yet some evolved new functions. The genomic complexity of the teleosts might be the reason for their evolutionary success and astounding biological diversity.     

How to Make a Rodent Giant: Genomic Basis and Tradeoffs of Gigantism in the Capybara,the World’s Largest Rodent

Gigantism results when one lineage within a clade evolves extremely large body size relative to its small-bodied ancestors, a common phenomenon in animals. Theory predicts that the evolution of giants should be constrained by two tradeoffs. First, because body size is negatively correlated with population size, purifying selection is expected to be less efficient in species of large body size, leading to increased mutational load. Second, gigantism is achieved through generating a higher number of cells along with higher rates of cell proliferation, thus increasing the likelihood of cancer. To explore the genetic basis of gigantism in rodents and uncover genomic signatures of gigantism-related tradeoffs, we assembled a draft genome of the capybara (Hydrochoerus hydrochaeris), the world’s largest living rodent. We found that the genome-wide ratio of nonsynonymous to synonymous mutations (ω) is elevated in the capybara relative to other rodents, likely caused by a generation-time effect and consistent with a nearly neutral model of molecular evolution. A genome-wide scan for adaptive protein evolution in the capybara highlighted several genes controlling postnatal bone growth regulation and musculoskeletal development, which are relevant to anatomical and developmental modifications for an increase in overall body size. Capybara-specific gene-family expansions included a putative novel anticancer adaptation that involves T-cell-mediated tumor suppression, offering a potential resolution to the increased cancer risk in this lineage. Our comparative genomic results uncovered the signature of an intragenomic conflict where the evolution of gigantism in the capybara involved selection on genes and pathways that are directly linked to cancer.     

Evolution of complexity in RNA-like replicator systems

Background  

The evolution of complexity is among the most important questions in biology. The evolution of complexity is often observed as the increase of genetic information or that of the organizational complexity of a system. It is well recognized that the formation of biological organization – be it of molecules or ecosystems – is ultimately instructed by the genetic information, whereas it is also true that the genetic information is functional only in the context of the organization. Therefore, to obtain a more complete picture of the evolution of complexity, we must study the evolution of both information and organization.     

The evolution of teleost pigmentation and the fish‐specific genome duplication

Teleost fishes have evolved a unique complexity and diversity of pigmentation and colour patterning that is unmatched among vertebrates. Teleost colouration is mediated by five different major types of neural‐crest derived pigment cells, while tetrapods have a smaller repertoire of such chromatophores. The genetic basis of teleost colouration has been mainly uncovered by the cloning of pigmentation genes in mutants of zebrafish Danio rerio and medaka Oryzias latipes. Many of these teleost pigmentation genes were already known as key players in mammalian pigmentation, suggesting partial conservation of the corresponding developmental programme among vertebrates. Strikingly, teleost fishes have additional copies of many pigmentation genes compared with tetrapods, mainly as a result of a whole‐genome duplication that occurred 320–350 million years ago at the base of the teleost lineage, the so‐called fish‐specific genome duplication. Furthermore, teleosts have retained several duplicated pigmentation genes from earlier rounds of genome duplication in the vertebrate lineage, which were lost in other vertebrate groups. It was hypothesized that divergent evolution of such duplicated genes may have played an important role in pigmentation diversity and complexity in teleost fishes, which therefore not only provide important insights into the evolution of the vertebrate pigmentary system but also allow us to study the significance of genome duplications for vertebrate biodiversity.     

Genetic parsimony: a factor in the evolution of complexity, order and emergence

Two conjectures, drawn from Gregory Chaitin's Algorithmic Information Theory, are examined with respect to the relationship between an algorithm and its product; in particular his finding that, where an algorithm is minimal, its length provides a measure of the complexity of the product. Algorithmic complexity is considered from the perspective of the relationship between genotype and phenotype, which Chaitin suggests is analogous to other algorithm-product systems. The first conjecture is that the genome is a minimal set of algorithms for the phenotype. Evidence is presented for a factor, here termed 'genetic parsimony', which is thought to have helped minimize the growth of genome size during evolution. Species that depend on rapid replication, such as prokaryotes which are generally r -selected are more likely to have small genomes, while the K -strategists accumulate introns and have large genomes. The second conjecture is that genome size could provide a measure of organism complexity. A surrogate index for coding DNA is in agreement with an established phenotypic index (number of cell types), in exhibiting an evolutionary trend of increasing organism complexity over time. Evidence for genetic parsimony indicates that simplicity in coding has been selected, and is responsible for phenotypic order. It is proposed that order evolved because order in the phenotype can be encoded more economically than disorder. Thus order arises due to selection for genetic parsimony, as does the evolution of other 'emergent' properties.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 295–308.     

When Less is More: Red Algae as Models for Studying Gene Loss and Genome Evolution in Eukaryotes

Genome evolution is usually viewed through the lens of growth in size and complexity over time, exemplified by plants and animals. In contrast, genome reduction is associated with a narrowing of ecological potential, such as in parasites and endosymbionts. But, can nuclear genome reduction also occur in, and potentially underpin a major radiation of free-living eukaryotes? An intriguing example of this phenomenon is provided by the red algae (Rhodophyta) that have lost many conserved pathways such as for flagellar motility, macroautophagy regulation, and phytochrome based light sensing. This anciently diverged, species-rich, and ecologically important algal lineage has undergone at least two rounds of large-scale genome reduction during its >1 billion-year evolutionary history. Here, using recent analyses of genome data, we review knowledge about the evolutionary trajectory of red algal nuclear and organelle gene inventories and plastid encoded autocatalytic introns. We compare and contrast Rhodophyta genome evolution to Viridiplantae (green algae and plants), both of which are members of the Archaeplastida, and highlight their divergent paths. We also discuss evidence for the speculative hypothesis that reduction in red algal plastid genome size through endosymbiotic gene transfer is counteracted by ongoing selection for compact nuclear genomes in red algae. Finally, we describe how the spliceosomal intron splicing apparatus provides an example of “evolution in action” in Rhodophyta and how the overall constraints on genome size in this lineage has left significant imprints on this key step in RNA maturation. Our review reveals the red algae to be an exciting, yet under-studied model that offers numerous novel insights as well as many unanswered questions that remain to be explored using modern genomic, genetic, and biochemical methods. The fact that a speciose lineage of free-living eukaryotes has spread throughout many aquatic habitats after having lost about 25% of its primordial gene inventory challenges us to elucidate the mechanisms underlying this remarkable feat.     

Three rounds (1R/2R/3R) of genome duplications and the evolution of the glycolytic pathway in vertebrates

Background  

Evolution of the deuterostome lineage was accompanied by an increase in systematic complexity especially with regard to highly specialized tissues and organs. Based on the observation of an increased number of paralogous genes in vertebrates compared with invertebrates, two entire genome duplications (2R) were proposed during the early evolution of vertebrates. Most glycolytic enzymes occur as several copies in vertebrate genomes, which are specifically expressed in certain tissues. Therefore, the glycolytic pathway is particularly suitable for testing theories of the involvement of gene/genome duplications in enzyme evolution.     

A computational approach to the functional screening of genomes

Comparative genomics usually involves managing the functional aspects of genomes, by simply comparing gene-by-gene functions. Following this approach, Mushegian and Koonin proposed a hypothetical minimal genome, Minimal Gene Set (MGS), aiming for a possible oldest ancestor genome. They obtained MGS by comparing the genomes of two simple bacteria and eliminating duplicated or functionally identical genes. The authors raised the fundamental question of whether a hypothetical organism possessing MGS is able to live or not. We attacked this viability problem specifying in silico the metabolic pathways of the MGS-based prokaryote. We then performed a dynamic simulation of cellular metabolic activities in order to check whether the MGS-prokaryote reaches some equilibrium state and produces the necessary biomass. We assumed these two conditions to be necessary for a living organism. Our simulations clearly show that the MGS does not express an organism that is able to live. We then iteratively proceeded with functional replacements in order to obtain a genome composition that gives rise to equilibrium. We ruled out 76 of the original 254 genes in the MGS, because they resulted in duplication from a functional point of view. We also added seven genes not present in the MGS. These genes encode for enzymes involved in critical nodes of the metabolic network. These modifications led to a genome composed of 187 elements expressing a virtually living organism, Virtual Cell (ViCe), that exhibits homeostatic capabilities and produces biomass. Moreover, the steady-state distribution of the concentrations of virtual metabolites that resulted was similar to that experimentally measured in bacteria. We conclude then that ViCe is able to “live in silico.”     

Evolutionary coupling between the deleteriousness of gene mutations and the amount of non-coding sequences

The phenotypic effects of random mutations depend on both the architecture of the genome and the gene-trait relationships. Both levels thus play a key role in the mutational variability of the phenotype, and hence in the long-term evolutionary success of the lineage. Here, by simulating the evolution of organisms with flexible genomes, we show that the need for an appropriate phenotypic variability induces a relationship between the deleteriousness of gene mutations and the quantity of non-coding sequences maintained in the genome. The more deleterious the gene mutations, the shorter the intergenic sequences. Indeed, in a shorter genome, fewer genes are affected by rearrangements (duplications, deletions, inversions, translocations) at each replication, which compensates for the higher impact of each gene mutation. This spontaneous adjustment of genome structure allows the organisms to retain the same average fitness loss per replication, despite the higher impact of single gene mutations. These results show how evolution can generate unexpected couplings between distinct organization levels.     

THE SMALLEST DINOFLAGELLATE GENOME IS YET TO BE FOUND: A COMMENT ON LAJEUNESSE ET AL. “SYMBIODINIUM (PYRRHOPHYTA) GENOME SIZES (DNA CONTENT) ARE SMALLEST AMONG DINOFLAGELLATES”1

LaJeunessse and colleagues (LaJeunesse et al. 2005) have recently documented small genome sizes of Symbiodinium and concluded that Symbiodinium is a dinoflagellate lineage with the smallest genome. The conclusion is inconsistent with recent discoveries of picoplanktonic dinoflagellates. The search for the smallest genome and the effort to understand the evolutionary history of dinoflagellate genome should be an area of research in the years to come, which can be greatly aided by an understanding on the current hypotheses regarding mechanisms of genome size evolution. Even the smallest dinoflagellate genome documented to date is too large to be sequenced with current technology, but sequencing of chromosomes or expressed genes of key representative species is feasible and can be very insightful for understanding genome composition and function in this important lineage of eukaryotes.     

The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The complete nucleotide sequences of the chloroplast (cp) and mitochondrial (mt) genomes of resurrection plant Boea hygrometrica (Bh, Gesneriaceae) have been determined with the lengths of 153,493 bp and 510,519 bp, respectively. The smaller chloroplast genome contains more genes (147) with a 72% coding sequence, and the larger mitochondrial genome have less genes (65) with a coding faction of 12%. Similar to other seed plants, the Bh cp genome has a typical quadripartite organization with a conserved gene in each region. The Bh mt genome has three recombinant sequence repeats of 222 bp, 843 bp, and 1474 bp in length, which divide the genome into a single master circle (MC) and four isomeric molecules. Compared to other angiosperms, one remarkable feature of the Bh mt genome is the frequent transfer of genetic material from the cp genome during recent Bh evolution. We also analyzed organellar genome evolution in general regarding genome features as well as compositional dynamics of sequence and gene structure/organization, providing clues for the understanding of the evolution of organellar genomes in plants. The cp-derived sequences including tRNAs found in angiosperm mt genomes support the conclusion that frequent gene transfer events may have begun early in the land plant lineage.     

Intron–Genome Size Relationship on a Large Evolutionary Scale

The intron–genome size relationship was studied across a wide evolutionary range (from slime mold and yeast to human and maize), as well as the relationship between genome size and the ratio of intervening/coding sequence size. The average intron size is scaled to genome size with a slope of about one-fourth for the log-transformed values; i.e., on the global scale its increase in evolution is lower than the increase in genome size by four orders of magnitude. There are exceptions to the general trend. In baker's yeast introns are extraordinarily long for its genome size. Tetrapods also have longer introns than expected for their genome sizes. In teleost fish the mean intron size does not differ significantly, notwithstanding the differences in genome size. In contrast to previous reports, avian introns were not found to be significantly shorter than introns of mammals, although avian genomes are smaller than genomes of mammals on average by about a factor of 2.5. The extra-/intragenic ratio of noncoding DNA can be higher in fungi than in animals, notwithstanding the smaller fungal genomes. In vertebrates and invertebrates taken separately, this ratio is increasing as the increase in genome size. Two hypotheses are proposed to explain the variation in the extra-/intragenic ratio of noncoding DNA in organisms with similar numbers of genes: transition (dynamic) and equilibrium (static). According to the transition model, this variation arises with the rapid shift of genome size because the bulk of extragenic DNA can be changed more rapidly than the finely interspersed intron sequences. The equilibrium model assumes that this variation is a result of selective adjustment of genome size with constraints imposed on the intron size due to its putative link to chromatin structure (and constraints of the splicing machinery). Received: 23 October 1997 / Accepted: 14 April 1999     

Mutational equilibrium model of genome size evolution

The paper describes a mutational equilibrium model of genome size evolution. This model is different from both adaptive and junk DNA models of genome size evolution in that it does not assume that genome size is maintained either by positive or stabilizing selection for the optimum genome size (as in adaptive theories) or by purifying selection against too much junk DNA (as in junk DNA theories). Instead the genome size is suggested to evolve until the loss of DNA through more frequent small deletions is equal to the rate of DNA gain through more frequent long insertions. The empirical basis for this theory is the finding of a strong correlation and of a clear power-function relationship between the rate of mutational DNA loss (per bp) through small deletions and genome size in animals. Genome size scales as a negative 1.3 power function of the deletion rate per nucleotide. Such a relationship is not predicted by either adaptive or junk DNA theories. However, if genome size is maintained at equilibrium by the balance of mutational forces, this empirilical relationship can be readily accommodated. Within this framework, this finding would imply that the rate of DNA gain through large insertions scales up a quarter-power function of genome size. On this view, as genome size grows, the rate of growth through large insertions is increasing as a quarter power function of genome size and the rate of DNA loss through small deletions increases linearly, until eventually, at the stable equilibrium genome size value, rates of growth and loss equal each other. The current data also suggest that the long-term variation is genome size in animals is brought about to a significant extent by changes in the intrinsic rates of DNA loss through small deletions. Both the origin of mutational biases and the adaptive consequences of such a mode of evolution of genome size are discussed.     

Lineage‐specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage‐specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage‐specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage‐specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.