Statistical Limitations on Molecular Evolution

Abstract

Complexity of functions evolving in an evolution process are expected to be limited by the time length of an evolution process among other factors. This paper outlines a general method of deriving function-complexity limitations based on mathematical statistics and independent from details of a biological or genetic mechanism of the evolution of the function. Limitations on the emergence of life are derived, these limitations indicate a possibility of a very fast evolution and are consistent with “RNA world” hypothesis. The discussed method is general and can be used to characterize evolution of more specific biological organism functions and relate functions to genetic structures. The derived general limitations indicate that a co-evolution of multiple functions and species could be a slow process, whereas an evolution of a specific function might proceed very fast, so that no trace of intermediate forms (species) is preserved in fossil records of phenotype or DNA structure; this is consistent with a picture of “punctuated equilibrium”.     

Breeding system variation, genetics and evolution in the Turneraceae

We review the genetics and evolution of breeding systems in the Turneraceae. Distyly occurs in seven of 10 genera and 81% of species. The remaining species are homostylous. Polyploid evolution has been significant in Turnera. Approximately 60% of species are polyploid ranging from diploid through decaploid. No relationship between breeding system and polyploidy is evident. The genetics of distyly involves a one-locus two-allele system (S and s). Evidence from crosses with homostylous species and mutants is consistent with the possibility that a "Primula-type" supergene underlies distyly but does not prove this to be the case. A polygalacturonase, and an alpha-dioxygenase specific to the transmitting tissue of short-styled plants both exhibit morph-limited expression in concert with predictions from an evolutionary model. The function of the proteins in distyly, if any, is unknown. We have begun constructing a fine-scale genetic map of Turnera. Two genetic markers lie within 0.2 cm of the distyly locus. This should provide a starting point for positional cloning of the distyly locus and reveal the genetic architecture and molecular basis of distyly.     

Evolution of structure,order and complexity in biological systems: Part III of a general theory

In this, Part III of a general theory, the large-scale features of evolution of structure, order, and complexity are considered as characteristic features of the biological state of matter. This starts with a rigorous formal definition of structure, classes of structural order, complexity, measures of complexity, and how these arise through evolution by a cumulative process of storing information in memory systems. Three such memory systems have evolved: the genetic memory, the immune memory, and the memories of the nervous system. The evolution, characteristic parameters and the limitations of these memory systems are explored. From these considerations emerge the large-scale features of the evolutionary pathways of biological structure, function, and complexity.     

Cancer Development and Progression: A Non-Adaptive Process Driven by Genetic Drift

The current mainstream in cancer research favours the idea that malignant tumour initiation is the result of a genetic mutation. Tumour development and progression is then explained as a sort of micro-evolutionary process, whereby an initial genetic alteration leads to abnormal proliferation of a single cell that leads to a population of clonally derived cells. It is widely claimed that tumour progression is driven by natural selection, based on the assumption that the initial tumour cells acquire some properties that endow such cells with a selective advantage over the normal cells from which the tumour cells are derived. The standard view assumes that the transformed bodily cell somehow acquires "responsiveness" to natural selection independently of the whole organism to which the cell belongs. Yet, it is never explained where such an acquired capacity to respond to natural selection by the individual bodily cell comes from. This situation poses many difficult questions that so far have been left unanswered. For example, there is no explanation why some cells belonging to an organised whole and as such having no independent capacity for survival, apparently become 'independent' entities, able to respond to selective pressures in an autonomous fashion and then to be evaluated by natural selection. Hereunder it is argued that such a qualitative change cannot be the consequence of specific genetic mutations. Moreover, it is shown that natural selection is unlikely to be acting within the organism during tumour development and progression and that tumour evolution is a random, non-adaptive process, driven by no fundamental biological principle. Thus, mutations in the so-called oncogenes and tumour suppressor genes observed in epithelial cancers (that constitute more than 90% of all cancers) are not the result of selection for better cellular growth or survival under restrictive conditions. Instead, here it is suggested that they are the consequence of genetic drift acting upon gene functions that become non-relevant, either for the individual or the species fitness, once the organism is past its reproductive prime and as such, they also become superfluous for cell survival in the short term. It is proposed that the origin of cancer is epigenetic and it is a consequence of the need for a continued turnover of the individuals that constitute a species.     

Inverse drug screens: a rapid and inexpensive method for implicating molecular targets

Identification of gene products that function in some specific process of interest is a common goal in developmental biology. Although use of drug compounds to probe biological systems has a very long history in teratology and toxicology, systematic hierarchical drug screening has not been capitalized upon by the developmental biology community. This "chemical genetics" approach can greatly benefit the study of embryonic and regenerative systems, and we have formalized a strategy for using known pharmacological compounds to implicate specific molecular candidates in any chosen biological phenomenon. Taking advantage of a hierarchical structure that can be imposed on drug reagents in a number of fields such as ion transport, neurotransmitter function, metabolism, and cytoskeleton, any assay can be carried out as a binary search algorithm. This inverse drug screen methodology is much more efficient than exhaustive testing of large numbers of drugs, and reveals the identity of a manageable number of specific molecular candidates that can then be validated and targeted using more expensive and specific molecular reagents. Here, we describe the process of this loss-of-function screen and illustrate its use in uncovering novel bioelectrical and serotonergic mechanisms in embryonic patterning. This technique is an inexpensive and rapid complement to existing molecular screening strategies. Moreover, it is applicable to maternal proteins, and model species in which traditional genetic screens are not feasible, significantly extending the opportunities to identify key endogenous players in biological processes.     

Chromosomal phylogeny and evolution of gibbons (Hylobatidae)

Although human and gibbons are classified in the same primate superfamily (Hominoidae), their karyotypes differ by extensive chromosome reshuffling. To date, there is still limited understanding of the events that shaped extant gibbon karyotypes. Further, the phylogeny and evolution of the twelve or more extant gibbon species (lesser apes, Hylobatidae) is poorly understood, and conflicting phylogenies have been published. We present a comprehensive analysis of gibbon chromosome rearrangements and a phylogenetic reconstruction of the four recognized subgenera based on molecular cytogenetics data. We have used two different approaches to interpret our data: (1) a cladistic reconstruction based on the identification of ancestral versus derived chromosome forms observed in extant gibbon species; (2) an approach in which adjacent homologous segments that have been changed by translocations and intra-chromosomal rearrangements are treated as discrete characters in a parsimony analysis (PAUP). The orangutan serves as an "outgroup", since it has a karyotype that is supposed to be most similar to the ancestral form of all humans and apes. Both approaches place the subgenus Bunopithecus as the most basal group of the Hylobatidae, followed by Hylobates, with Symphalangus and Nomascus as the last to diverge. Since most chromosome rearrangements observed in gibbons are either ancestral to all four subgenera or specific for individual species and only a few common derived rearrangements at subsequent branching points have been recorded, all extant gibbons may have diverged within relatively short evolutionary time. In general, chromosomal rearrangements produce changes that should be considered as unique landmarks at the divergence nodes. Thus, molecular cytogenetics could be an important tool to elucidate phylogenies in other species in which speciation may have occurred over very short evolutionary time with not enough genetic (DNA sequence) and other biological divergence to be picked up.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Gene-based sequence diversity analysis of field pea (Pisum)

Sequence diversity of 39 dispersed gene loci was analyzed in 48 diverse individuals representative of the genus Pisum. The different genes show large variation in diversity parameters, suggesting widely differing levels of selection and a high overall diversity level for the species. The data set yields a genetic diversity tree whose deep branches, involving wild samples, are preserved in a tree derived from a polymorphic retrotransposon insertions in an identical sample set. Thus, gene regions and intergenic "junk DNA" share a consistent picture for the genomic diversity of Pisum, despite low linkage disequilibrium in wild and landrace germplasm, which might be expected to allow independent evolution of these very different DNA classes. Additional lines of evidence indicate that recombination has shuffled gene haplotypes efficiently within Pisum, despite its high level of inbreeding and widespread geographic distribution. Trees derived from individual gene loci show marked differences from each other, and genetic distance values between sample pairs show high standard deviations. Sequence mosaic analysis of aligned sequences identifies nine loci showing evidence for intragenic recombination. Lastly, phylogenetic network analysis confirms the non-treelike structure of Pisum diversity and indicates the major germplasm classes involved. Overall, these data emphasize the artificiality of simple tree structures for representing genomic sequence variation within Pisum and emphasize the need for fine structure haplotype analysis to accurately define the genetic structure of the species.     

Genetic variation: molecular mechanisms and impact on microbial evolution

On the basis of established knowledge of microbial genetics one can distinguish three major natural strategies in the spontaneous generation of genetic variations in bacteria. These strategies are: (1) small local changes in the nucleotide sequence of the genome, (2) intragenomic reshuffling of segments of genomic sequences and (3) the acquisition of DNA sequences from another organism. The three general strategies differ in the quality of their contribution to microbial evolution. Besides a number of non-genetic factors, various specific gene products are involved in the generation of genetic variation and in the modulation of the frequency of genetic variation. The underlying genes are called evolution genes. They act for the benefit of the biological evolution of populations as opposed to the action of housekeeping genes and accessory genes which are for the benefit of individuals. Examples of evolution genes acting as variation generators are found in the transposition of mobile genetic elements and in so-called site-specific recombination systems. DNA repair systems and restriction-modification systems are examples of modulators of the frequency of genetic variation. The involvement of bacterial viruses and of plasmids in DNA reshuffling and in horizontal gene transfer is a hint for their evolutionary functions. Evolution genes are thought to undergo biological evolution themselves, but natural selection for their functions is indirect, at the level of populations, and is called second-order selection. In spite of an involvement of gene products in the generation of genetic variations, evolution genes do not programmatically direct evolution towards a specific goal. Rather, a steady interplay between natural selection and mixed populations of genetic variants gives microbial evolution its direction.     

Inferring phylogenetic networks by the maximum parsimony criterion: a case study

Horizontal gene transfer (HGT) may result in genes whose evolutionary histories disagree with each other, as well as with the species tree. In this case, reconciling the species and gene trees results in a network of relationships, known as the "phylogenetic network" of the set of species. A phylogenetic network that incorporates HGT consists of an underlying species tree that captures vertical inheritance and a set of edges which model the "horizontal" transfer of genetic material. In a series of papers, Nakhleh and colleagues have recently formulated a maximum parsimony (MP) criterion for phylogenetic networks, provided an array of computationally efficient algorithms and heuristics for computing it, and demonstrated its plausibility on simulated data. In this article, we study the performance and robustness of this criterion on biological data. Our findings indicate that MP is very promising when its application is extended to the domain of phylogenetic network reconstruction and HGT detection. In all cases we investigated, the MP criterion detected the correct number of HGT events required to map the evolutionary history of a gene data set onto the species phylogeny. Furthermore, our results indicate that the criterion is robust with respect to both incomplete taxon sampling and the use of different site substitution matrices. Finally, our results show that the MP criterion is very promising in detecting HGT in chimeric genes, whose evolutionary histories are a mix of vertical and horizontal evolution. Besides the performance analysis of MP, our findings offer new insights into the evolution of 4 biological data sets and new possible explanations of HGT scenarios in their evolutionary history.     

Evolutionary plasticity and functional versatility of CRISPR systems

The principal biological function of bacterial and archaeal CRISPR systems is RNA-guided adaptive immunity against viruses and other mobile genetic elements (MGEs). These systems show remarkable evolutionary plasticity and functional versatility at multiple levels, including both the defense mechanisms that lead to direct, specific elimination of the target DNA or RNA and those that cause programmed cell death (PCD) or induction of dormancy. This flexibility is also evident in the recruitment of CRISPR systems for nondefense functions. Defective CRISPR systems or individual CRISPR components have been recruited by transposons for RNA-guided transposition, by plasmids for interplasmid competition, and by viruses for antidefense and interviral conflicts. Additionally, multiple highly derived CRISPR variants of yet unknown functions have been discovered. A major route of innovation in CRISPR evolution is the repurposing of diverged repeat variants encoded outside CRISPR arrays for various structural and regulatory functions. The evolutionary plasticity and functional versatility of CRISPR systems are striking manifestations of the ubiquitous interplay between defense and “normal” cellular functions.

The CRISPR systems show remarkable functional versatility beyond their principal function as an adaptive immune mechanism. This Essay discusses how derived CRISPR systems have been recruited by transposons on multiple occasions and mediate RNA-guided transposition; derived CRISPR RNAs are frequently recruited for regulatory functions.     

Inference of miRNA targets using evolutionary conservation and pathway analysis

Background

MicroRNAs have emerged as important regulatory genes in a variety of cellular processes and, in recent years, hundreds of such genes have been discovered in animals. In contrast, functional annotations are available only for a very small fraction of these miRNAs, and even in these cases only partially.

Results

We developed a general Bayesian method for the inference of miRNA target sites, in which, for each miRNA, we explicitly model the evolution of orthologous target sites in a set of related species. Using this method we predict target sites for all known miRNAs in flies, worms, fish, and mammals. By comparing our predictions in fly with a reference set of experimentally tested miRNA-mRNA interactions we show that our general method performs at least as well as the most accurate methods available to date, including ones specifically tailored for target prediction in fly. An important novel feature of our model is that it explicitly infers the phylogenetic distribution of functional target sites, independently for each miRNA. This allows us to infer species-specific and clade-specific miRNA targeting. We also show that, in long human 3' UTRs, miRNA target sites occur preferentially near the start and near the end of the 3' UTR. To characterize miRNA function beyond the predicted lists of targets we further present a method to infer significant associations between the sets of targets predicted for individual miRNAs and specific biochemical pathways, in particular those of the KEGG pathway database. We show that this approach retrieves several known functional miRNA-mRNA associations, and predicts novel functions for known miRNAs in cell growth and in development.

Conclusion

We have presented a Bayesian target prediction algorithm without any tunable parameters, that can be applied to sequences from any clade of species. The algorithm automatically infers the phylogenetic distribution of functional sites for each miRNA, and assigns a posterior probability to each putative target site. The results presented here indicate that our general method achieves very good performance in predicting miRNA target sites, providing at the same time insights into the evolution of target sites for individual miRNAs. Moreover, by combining our predictions with pathway analysis, we propose functions of specific miRNAs in nervous system development, inter-cellular communication and cell growth. The complete target site predictions as well as the miRNA/pathway associations are accessible on the ElMMo web server.     

The utility and limitations of genetic data for stock identification and management of North Pacific rockfish (Sebastes spp.)

Rockfishes (Sebastes spp.) represent a speciose and ecologically important group of marine fishes found in both the Pacific and Atlantic oceans, with approximately 105 species found world-wide (Hyde and Vetter 2007). They also comprise the majority of species found in the Pacific groundfish fishery. Thorough species assessments in terms of harvest management have been done for only 11 species, and of the 11 species, seven have been declared overfished. Having accurate genetic information is critical to the continuing effort at stock assessments, but sampling is often difficult in marine fishes. Genetic techniques are a powerful tool in the effort to better characterize the ecology of these species. These techniques can be used to investigate multiple biological traits, including species identity, intra- and interspecific genetic variation, migration patterns, and effective population size. There are important caveats and limitations when applying specific genetic methods, especially in marine species that lack discrete spawning aggregates. Nevertheless, it is clear from a review of recent literature that genetic tools have already provided very specific insight regarding rockfish population dynamics. The results are diverse and difficult to synthesize; however, existing studies show five primary patterns to population groupings in rockfishes: no obvious pattern of structure, structure consistent with isolation by distance, structure evident but inconsistent with isolation by distance, structure that correlates to oceanographic features, and potential genetic introgression. Clearly the study of rockfish population genetics is poised for rapid expansion that will unquestionably aid management of the rockfish fisheries and general understanding of rockfish evolutionary systematics. A principle challenge at this point is to derive generalized inferences from such a diverse array of study results across the vast North Pacific range of Sebastes. This review summarizes existing genetic studies in Sebastes spp. in the North Pacific to assist in identifying knowledge gaps for this ecologically important and diverse group.     

Molecular clock and gene function

Molecular phylogenies based on the molecular clock require the comparison of orthologous genes. Orthologous and paralogous genes usually have very different evolutionary fates. In general, orthologs keep the same functions in species, whereas, particularly over a long time span, paralogs diverge functionally and may become pseudogenes or get lost. In eukaryotic genomes, because of the degree of redundancy of genetic information, homologous genes are grouped in gene families, the evolution of which may differ greatly between the various organisms. This implies that each gene in a species does not always have an ortholog in another species and thus, due to multiple duplication events following a speciation, many orthologous clades of paralogs are generated. We are often dealing with a one-to-many or many-to-many relationship between genes. In this paper, we analyze the evolution of two gene families, the p53 gene family and the porin gene family. The evolution of the p53 family shows a one-to-many gene relationship going from invertebrates to vertebrates. In invertebrates only a single gene has been found, while in vertebrates three members of the family, namely p53, p63, and p73, are present. The evolution of porin (VDAC) genes (VDAC1, VDAC2, and VDAC3) is an example of a many-to-many gene relationship going from yeast to mammals. However, the porin gene redundancy found in invertebrates and possibly in some fishes may indicate a tendency to duplicate the genetic material, rather than a real need for function innovation.     

Clustering and visualizing similarity networks of membrane proteins

We proposed a fast and unsupervised clustering method, minimum span clustering (MSC), for analyzing the sequence–structure–function relationship of biological networks, and demonstrated its validity in clustering the sequence/structure similarity networks (SSN) of 682 membrane protein (MP) chains. The MSC clustering of MPs based on their sequence information was found to be consistent with their tertiary structures and functions. For the largest seven clusters predicted by MSC, the consistency in chain function within the same cluster is found to be 100%. From analyzing the edge distribution of SSN for MPs, we found a characteristic threshold distance for the boundary between clusters, over which SSN of MPs could be properly clustered by an unsupervised sparsification of the network distance matrix. The clustering results of MPs from both MSC and the unsupervised sparsification methods are consistent with each other, and have high intracluster similarity and low intercluster similarity in sequence, structure, and function. Our study showed a strong sequence–structure–function relationship of MPs. We discussed evidence of convergent evolution of MPs and suggested applications in finding structural similarities and predicting biological functions of MP chains based on their sequence information. Proteins 2015; 83:1450–1461. © 2015 Wiley Periodicals, Inc.     

Evolution of longevity in primates

Evolution of maximum lifespan potential (MLP) and its possible relation to the evolution of extra-brain functions are investigated in the primates. MLP is related to a species' characteristic aging rate and is considered a basic biological property of an organism. MLP may be of more importance than realized in the past in determining the evolutionary success of a species. It is related to the postnatal development rate, the length of time general vigor is maintained, the length of reproductive period, generation time and the time available for learning and teaching behavior. Three other parameters are considered to be importantly related to MLP: MLP calorie consumption (MCC), encephalization quotient (EQ) and extra number of cortical neurons (N_c). MCC is calculated as the product of MLP and specific metabolic rate (SMR) and is considered to represent total “life-capacity” of an organism. It is of potential value in studying the biological mechanisms involved in the evolution of MLP. Brain function is estimated by the Jerison EQ and N_c parameters. These parameters estimate the “extra” brain capacity involved in functions beyond normal body requirements. The rate of change in MLP and N_c per unit time occurring during an ancestral-descendant sequence is used to estimate the biological complexity of the genetic processes which have evolved in governing the rate of expression of the general aging process and increasing brain function. The average rate of change of MLP during the emergence of the primates was analyzed by the difference between MLP in closely related living primate species and the evolutionary time of appearance of a common ancestor. MLP, SMR, MCC, EQ and N_c were estimated from fossil cranial capacity and body weight measurements. The rate of change in these values was calculated according to the time of appearance of the fossil species. MLP and N_c were found to increase together and reached their highest rate of increase approximately 200,000 years ago along the hominid ancestral-descendant sequence leading to modern man. The high rate of increase of these parameters suggests that few genetic changes were responsible. The general increase in MLP during the evolution of the primate species indicates that a corresponding general decrease in mutation rate may have occurred. The high levels of MLP, MCC, EQ and N_c represented in the living primates, as compared to other mammals, are considered to represent a major characteristic determining their evolutionary success.     

Transposable elements, gene creation and genome rearrangement in flowering plants

Plant genome structure is largely derived from the differing specificities, abundances and activities of transposable elements. Recent studies indicate that both the amplification and the removal of transposons are rapid processes in plants, accounting for the general lack of intergenic homology between species that last shared a common ancestor more than 10 million years ago. Two newly discovered transposon varieties, Helitrons and Pack-MULEs, acquire and fuse fragments of plant genes, creating the raw material for the evolution of new genes and new genetic functions. Many of these recently assembled, chimeric gene-candidates are expressed, suggesting that some might escape epigenetic silencing and mutational decay, but a proven case of gene creation by any transposable element activity in plants remains to be demonstrated.     

The fittest versus the flattest: experimental confirmation of the quasispecies effect with subviral pathogens

The "survival of the fittest" is the paradigm of Darwinian evolution in which the best-adapted replicators are favored by natural selection. However, at high mutation rates, the fittest organisms are not necessarily the fastest replicators but rather are those that show the greatest robustness against deleterious mutational effects, even at the cost of a low replication rate. This scenario, dubbed the "survival of the flattest", has so far only been shown to operate in digital organisms. We show that "survival of the flattest" can also occur in biological entities by analyzing the outcome of competition between two viroid species coinfecting the same plant. Under optimal growth conditions, a viroid species characterized by fast population growth and genetic homogeneity outcompeted a viroid species with slow population growth and a high degree of variation. In contrast, the slow-growth species was able to outcompete the fast species when the mutation rate was increased. These experimental results were supported by an in silico model of competing viroid quasispecies.     

Genomic exploration of the hemiascomycetous yeasts: 21. Comparative functional classification of genes

We explored the biological diversity of hemiascomycetous yeasts using a set of 22000 newly identified genes in 13 species through BLASTX searches. Genes without clear homologue in Saccharomyces cerevisiae appeared to be conserved in several species, suggesting that they were recently lost by S. cerevisiae. They often identified well-known species-specific traits. Cases of gene acquisition through horizontal transfer appeared to occur very rarely if at all. All identified genes were ascribed to functional classes. Functional classes were differently represented among species. Species classification by functional clustering roughly paralleled rDNA phylogeny. Unequal distribution of rapidly evolving, ascomycete-specific, genes among species and functions was shown to contribute strongly to this clustering. A few cases of gene family amplification were documented, but no general correlation could be observed between functional differentiation of yeast species and variations of gene family sizes. Yeast biological diversity seems thus to result from limited species-specific gene losses or duplications, and for a large part from rapid evolution of genes and regulatory factors dedicated to specific functions.