Simple stochastic birth and death models of genome evolution: was there enough time for us to evolve?

MOTIVATION: The distributions of many genome-associated quantities, including the membership of paralogous gene families can be approximated with power laws. We are interested in developing mathematical models of genome evolution that adequately account for the shape of these distributions and describe the evolutionary dynamics of their formation. RESULTS: We show that simple stochastic models of genome evolution lead to power-law asymptotics of protein domain family size distribution. These models, called Birth, Death and Innovation Models (BDIM), represent a special class of balanced birth-and-death processes, in which domain duplication and deletion rates are asymptotically equal up to the second order. The simplest, linear BDIM shows an excellent fit to the observed distributions of domain family size in diverse prokaryotic and eukaryotic genomes. However, the stochastic version of the linear BDIM explored here predicts that the actual size of large paralogous families is reached on an unrealistically long timescale. We show that introduction of non-linearity, which might be interpreted as interaction of a particular order between individual family members, allows the model to achieve genome evolution rates that are much better compatible with the current estimates of the rates of individual duplication/loss events.     

Are There Laws of Genome Evolution?

Research in quantitative evolutionary genomics and systems biology led to the discovery of several universal regularities connecting genomic and molecular phenomic variables. These universals include the log-normal distribution of the evolutionary rates of orthologous genes; the power law-like distributions of paralogous family size and node degree in various biological networks; the negative correlation between a gene's sequence evolution rate and expression level; and differential scaling of functional classes of genes with genome size. The universals of genome evolution can be accounted for by simple mathematical models similar to those used in statistical physics, such as the birth-death-innovation model. These models do not explicitly incorporate selection; therefore, the observed universal regularities do not appear to be shaped by selection but rather are emergent properties of gene ensembles. Although a complete physical theory of evolutionary biology is inconceivable, the universals of genome evolution might qualify as "laws of evolutionary genomics" in the same sense "law" is understood in modern physics.     

A theoretical framework for defining some concepts in evolution.

We present a theoretical framework for biological evolution with the intention of giving precise mathematical definitions of some concepts in evolutionary biology such as fitness, evolutionary pressure, specialization and natural selection. In this framework, such concepts are identified with well-known mathematical terms within the theory of dynamical systems. We also discuss some more general implications in evolution; for instance, the fact that our model naturally exhibits a frequency spectrum of the type 1/f for low frequencies of evolutionary events.     

Birth-and-death evolution of the internalin multigene family in Listeria

The birth-and-death model of multigene family evolution describes patterns of gene origination, diversification and loss within multigene families. Since it was first developed in the 1990s, the model has been found to characterize a large number of eukaryotic multigene families. In this paper, we report for the first time a bacterial multigene family that undergoes birth-and-death evolution. By analyzing the evolutionary relationships among internalins, a relatively large and diverse family of genes associated with key virulence functions in Listeria, we demonstrate the importance of birth-and-death evolution in the diversification of this important bacterial pathogen. We also detected two instances of lateral gene transfer within the internalins, but the estimated frequency would have been much higher had it not been analyzed within the context of birth-and-death evolutionary dynamics and a phenomenon that we term "paralog-sorting", which involves the unequal transmittal of gene duplicates during or subsequent to the speciation process. As such, in addition to providing the first demonstration of birth-and-death evolution within a bacterial multigene family, our results indicate that the extent of lateral transfer in bacterial multigene families should be re-examined in the light of birth-and-death evolution.     

Evolution and phylogeny of the heterogeneous cytosolic SSU rRNA genes in the genus Plasmodium

Unlike other eukaryotes, malaria parasites in the genus Plasmodium have structurally and functionally different paralogous copies of the cytosolic (cyto-) SSU rRNA (18S rRNA) gene that are expressed at different developmental stages. In P. falciparum, P. vivax, and P. berghei, A-type cyto-SSU rRNA is expressed in asexual stage, while S-type in sporozoite stage. A third type (O-type) has been described in P. vivax. It is expressed only in oocyst stage in the mosquito. Recently, it has been shown that the maintenance of heterogeneous cyto-SSU rRNAs in Plasmodium can be modeled as a birth-and-death process under strong purifying selection [Rooney, A.P., 2004. Mechanisms underlying the evolution and maintenance of functionally heterogeneous 18S rRNA genes in Apicomplexans. Mol. Biol. Evol. 21, 1704-1711]. In this study, we performed detailed phylogenetic analyses of Plasmodium cyto-SSU rRNAs with special emphasis on the evolution of multi-copy genes in simian Plasmodium species. We sequenced paralogous copies of the cyto-SSU rRNA genes from an African simian Plasmodium species, P. gonderi, and Asian simian Plasmodium species, P. fragile, P. coatneyi, P. inui, P. hylobati, P. fieldi, P. simiovale, and P. cynomolgi. Interestingly, all Asian simian Plasmodium species have a single S-type-like gene and several A-type-like genes. Alignment analysis demonstrated for the first time that an approximately 50-residue insertion in the V7 variable region near the stem 43 is shared exclusively by the S-type-like sequences of the Asian simian Plasmodium species and the S- and O-type sequences of P. vivax. We comprehensively analyzed all cyto-SSU rRNA sequences of the genus Plasmodium currently available in the database. Phylogenetic analyses of all publicly available cyto-SSU rRNA sequences for the genus Plasmodium clearly demonstrated that gene duplication events giving rise to A- and S-type-like sequences took place independently at least three times in the Plasmodium evolution, supporting the hypothesis that these genes evolve according to a birth-and-death model.     

Concerted and nonconcerted evolution of the Hsp70 gene superfamily in two sibling species of nematodes

We have identified the Hsp70 gene superfamily of the nematode Caenorhabditis briggsae and investigated the evolution of these genes in comparison with Hsp70 genes from C. elegans, Drosophila, and yeast. The Hsp70 genes are classified into three monophyletic groups according to their subcellular localization, namely, cytoplasm (CYT), endoplasmic reticulum (ER), and mitochondria (MT). The Hsp110 genes can be classified into the polyphyletic CYT group and the monophyletic ER group. The different Hsp70 and Hsp110 groups appeared to evolve following the model of divergent evolution. This model can also explain the evolution of the ER and MT genes. On the other hand, the CYT genes are divided into heat-inducible and constitutively expressed genes. The constitutively expressed genes have evolved more or less following the birth-and-death process, and the rates of gene birth and gene death are different between the two nematode species. By contrast, some heat-inducible genes show an intraspecies phylogenetic clustering. This suggests that they are subject to sequence homogenization resulting from gene conversion-like events. In addition, the heat-inducible genes show high levels of sequence conservation in both intra-species and inter-species comparisons, and in most cases, amino acid sequence similarity is higher than nucleotide sequence similarity. This indicates that purifying selection also plays an important role in maintaining high sequence similarity among paralogous Hsp70 genes. Therefore, we suggest that the CYT heat-inducible genes have been subjected to a combination of purifying selection, birth-and-death process, and gene conversion-like events.     

On the incidence of intron loss and gain in paralogous gene families

Understanding gene duplication and gene structure evolution are fundamental goals of molecular evolutionary biology. A previous study by Babenko et al. (2004. Prevalence of intron gain over intron loss in the evolution of paralogous gene families. Nucleic Acids Res. 32:3724-3733) employed Dollo parsimony to infer spliceosomal intron losses and gains in paralogous gene families and concluded that there was a general excess of gains over losses. This result contrasts with patterns in orthologous genes, in which most lineages show an excess of intron losses over gains, suggesting the possibility of fundamentally different modes of intron evolution between orthologous and paralogous genes. We further studied the data and found a low level of intron position conservation with outgroups, and this led to problems with using Dollo parsimony to analyze the data. Statistical reanalysis of the data suggests, instead, that intron losses have outnumbered intron gains in paralogous gene families.     

Darwinian evolution in the light of genomics

Comparative genomics and systems biology offer unprecedented opportunities for testing central tenets of evolutionary biology formulated by Darwin in the Origin of Species in 1859 and expanded in the Modern Synthesis 100 years later. Evolutionary-genomic studies show that natural selection is only one of the forces that shape genome evolution and is not quantitatively dominant, whereas non-adaptive processes are much more prominent than previously suspected. Major contributions of horizontal gene transfer and diverse selfish genetic elements to genome evolution undermine the Tree of Life concept. An adequate depiction of evolution requires the more complex concept of a network or ‘forest’ of life. There is no consistent tendency of evolution towards increased genomic complexity, and when complexity increases, this appears to be a non-adaptive consequence of evolution under weak purifying selection rather than an adaptation. Several universals of genome evolution were discovered including the invariant distributions of evolutionary rates among orthologous genes from diverse genomes and of paralogous gene family sizes, and the negative correlation between gene expression level and sequence evolution rate. Simple, non-adaptive models of evolution explain some of these universals, suggesting that a new synthesis of evolutionary biology might become feasible in a not so remote future.     

Modeling the evolution of interconnected processes: It is the song and the singers

Recently, Doolittle and Inkpen formulated a thought provoking theory, asserting that evolution by natural selection was responsible for the sideways evolution of two radically different kinds of selective units (also called Domains). The former entities, termed singers, correspond to the usual objects studied by evolutionary biologists (gene, genomes, individuals, species, etc.), whereas the later, termed songs, correspond to re‐produced biological and ecosystemic functions, processes, information, and memes. Singers perform songs through selected patterns of interactions, meaning that a wealth of critical phenomena might receive novel evolutionary explanations. However, this theory did not provide an empirical approach to study evolution in such a broadened context. Here, we show that analyzing songs and singers, using patterns of interaction networks as a common ontology for both, offers a novel, actionable, inclusive and mathematical way to analyze not only the re‐production but also the evolution and fitness of biological and ecosystemic interconnected processes.     

Mechanisms underlying the evolution and maintenance of functionally heterogeneous 18S rRNA genes in Apicomplexans

In many species of the protist phylum Apicomplexa, ribosomal RNA (rRNA) gene copies are structurally and functionally heterogeneous, owing to distinct requirements for rRNA-expression patterns at different developmental stages. The genomic mechanisms underlying the maintenance of this system over long-term evolutionary history are unclear. Therefore, the aim of this study was to investigate what processes underlie the long-term evolution of apicomplexan 18S genes in representative species. The results show that these genes evolve according to a birth-and-death model under strong purifying selection, thereby explaining how divergent 18S genes are generated over time while continuing to maintain their ability to produce fully functional rRNAs. In addition, it was found that Cryptosporidium parvum undergoes a rapid form of birth-and-death evolution that may facilitate host-specific adaptation, including that of type I and II strains found in humans. This represents the first case in which an rRNA gene family has been found to evolve under the birth-and-death model.     

The cell as nexus: connections between the history, philosophy and science of cell biology

Although the cell is commonly addressed as the unit of life, historians and philosophers have devoted relatively little attention to this concept in comparison to other fundamental concepts of biology such as the gene or species. As a partial remedy to this neglect, we introduce the cell as a major point of connection between various disciplinary approaches, epistemic strategies, technological vectors and overarching biological processes such as metabolism, growth, reproduction and evolution. We suggest that the role of the cell as a nexus forms the basis for a new philosophical and historical appreciation of cell biology. This perspective focuses less on the cell as a well-defined, stable object and places more emphasis on its role as a mediator of fundamental biological processes.     

分析水稻APETALA2类基因的系统发育和表达模式

多基因家族经历着生与死的进化并因此有助于生物的创新性进化。APETALA2类基因属于AP2转录因子家族中euAP2分支。该类基因具有特征性基元,并存在于蕨类植物、裸子植物和被子植物中。系统发育分析表明这类基因经历了生与死的进化。原位杂交和RT-PCR结果表明水稻的5个APETALA2类基因在表达模式上既是冗余又是不同的。根据基因系统发育树和表达模式讨论了这些基因的功能。     

Origins and divergence times of mammalian class II MHC gene clusters

The class I and II major histocompatibility complex (MHC) genes are apparently subject to evolution by a birth-and-death process. The rate of gene turnover is much slower in the latter genes than in the former. In placental mammals, the class II region can be subdivided into different orthologous subregions or gene clusters (DR, DQ, DO, and DN), but the origins and evolutionary relationships of these gene clusters are not well established. Here we report the results of our study of the times of origin and evolutionary relationships of these gene clusters in mammals. Our analysis suggests that both class II alpha-chain and beta-chain gene clusters are shared by placental mammals and marsupials, but the gene clusters from nonmammalian species are paralogous to mammalian gene clusters. We estimated the times of divergence between gene clusters in placental mammals using the linearized tree and distance regression methods. Our results indicate that most gene clusters originated 170-200 million years (MY) ago, but that DO beta-chain genes diverged from the other beta-chain gene clusters approximately 210-260 MY ago. The phylogenetic trees for the alpha- and beta-chain genes were not congruent, suggesting that the evolutionary history of the class II gene clusters is more complex than previously thought.     

Evolutionary dynamics of genes controlling floral development

Advances in the understanding of floral developmental genetics in model species such as Arabidopsis continue to provide an important foundation for comparative studies in other flowering plants. In particular, floral organ identity genes are the focus of many projects that are addressing both ancient and recent evolutionary questions. Expanded analyses of the evolution of these gene lineages have highlighted the dynamic nature of the gene birth-and-death process, and may have significant implications for the evolution of genetic pathways. Crucial functional studies of floral organ identity genes in diverse taxa are allowing the first real insight into the conservation of gene function, while findings on the genetic control of organ elaboration offer to open up new avenues for investigation. Taken together, these trends show that the field of floral developmental evolution continues to make significant progress towards elucidating the processes that have shaped the evolution of flower development and morphology.     

花、基因、禾本科

进化发育生物学的一个重要任务就是揭示形态多样性的分子基础, 该领域的研究包含形态、形态发育相关基因和形态所属类群等三个要素。花/花序是进化发育生物学研究的首要对象, 系统发育重建和个体发育剖析的结合将促进认知花的形态进化。发育相关基因的进化表现为等位基因遗传或表观遗传的突变, 基因家族生与死的进化, 不同基因组拥有独特的基因。运用形态学或序列分析方法很大程度揭示了禾本科植物花进化过程中的基因进化。试从学科问题、思路方法以及具体例子介绍植物进化发育生物学。     

花、基因、禾本科

进化发育生物学的一个重要任务就是揭示形态多样性的分子基础,该领域的研究包含形态、形态发育相关基因和形态所属类群等三个要素。花/花序是进化发育生物学研究的首要对象,系统发育重建和个体发育剖析的结合将促进认知花的形态进化。发育相关基因的进化表现为等位基因遗传或表观遗传的突变,基因家族生与死的进化,不同基因组拥有独特的基因。运用形态学或序列分析方法很大程度揭示了禾本科植物花进化过程中的基因进化。试从学科问题、思路方法以及具体例子介绍植物进化发育生物学。     

On a Birth-and-Death Process Induced Distribution

A type of distribution induced by the linear birth-and-death processes is useful in modeling certain biological phenomena. In addition to extending some existing results, new findings are presented concerning the mathematical properties of the distribution.     

Systematic investigation of Amphioxus (Branchiostoma floridae) microRNAs

MicroRNAs (miRNAs) participate in various biological processes via controlling gene activity. Amphioxus is the best available stand-in as the proximate invertebrate ancestor of the vertebrates. Here, we systematically investigated the miRNAs in amphioxus. First, we identified 245 candidate amphioxus miRNAs, in which 183 miRNAs were firstly reported. Second, we gave evidences to support a birth-and-death process of miRNA genes in some families and gave implications for the functional diversification of miRNA during evolution. Third, we identified 47 development-specific expression miRNAs. We found that only 19 miRNAs were expressed in all developmental stages, 16 miRNAs were neurula-specific and 13 miRNAs were larva-specific. In addition, these potential miRNA-targeting genes were mainly classified into development, muscle formation, cell adhesion, and gene regulation categories. Finally, we found 79 immune related genes targeted by 136 miRNAs in amphioxus. In conclusion, our results take an insight into both the function and evolution of the amphioxus miRNAs.     

Character identification in evolutionary biology: The role of the organism

Summary In this article we argue that an organismic perspective in character identification can alleviate a structural deficiency of mathematical models in biology relative to the ones in the physical sciences. The problem with many biological theories is that they do not contain the conditions of their validity or a method of identifying objects that are appropriate instances of the models. Here functionally important biological characters are introduced as conceptual abstractions derived within the context of an ontologically prior object, such as a cell or an organism. To illustrate this approach, we present an analytical method of character decomposition based on the notion of the quasi-independence of traits. Two cases are analyzed: context dependent units of inheritance and a model of character identification in adaptive evolution. We demonstrate that in each case the biological process as represented by a mathematical theory entails the conditions for the individualization of characters. Our approach also requires a conceptual re-orientation in the way we build biological models. Rather than defining a set of biological characters a priori, functionally relevant characters are identified in the context of a higher level biological process.     

A unifying view of 21st century systems biology

The idea that multi-scale dynamic complex systems formed by interacting macromolecules and metabolites, cells, organs and organisms underlie some of the most fundamental aspects of life was proposed by a few visionaries half a century ago. We are witnessing a powerful resurgence of this idea made possible by the availability of nearly complete genome sequences, ever improving gene annotations and interactome network maps, the development of sophisticated informatic and imaging tools, and importantly, the use of engineering and physics concepts such as control and graph theory. Alongside four other fundamental “great ideas” as suggested by Sir Paul Nurse, namely, the gene, the cell, the role of chemistry in biological processes, and evolution by natural selection, systems-level understanding of “What is Life” may materialize as one of the major ideas of biology.