Evaluating alternative hypotheses for the origin of eusociality in corbiculate bees

We use a likelihood-based statistical test to evaluate the extent to which the available molecular data sets can be used to falsify alternative phylogenetic hypotheses describing the inter-relationship among corbiculate bee tribes. Based on the results of this test, we explore three alternative models of behavioural character state evolution and evaluate the support each model has for single-origin versus dual-origin hypotheses for 'highly' eusocial behaviour. We show that only one of four data sets could statistically reject any of the 15 possible outgroup-rooted phylogenetic hypotheses. However, a cytochrome b data set rejected all but three alternative topologies. Using this information, a simple model of behavioural character state evolution, in which transitions between solitary/communal, 'primitively' eusocial, and 'highly' eusocial are unconstrained, supports single-origin hypotheses for 'highly' eusocial behaviour, in spite of phylogenetic uncertainty. By contrast, an ordered model, in which 'highly' eusocial is constrained to be an evolutionarily terminal state, supports a dual-origins hypothesis. Our results show that the molecular phylogenetic evidence favouring a dual-origins hypothesis for 'highly' eusocial behaviour is, at present, conditional on information from one gene (cyt b) and on specific, though likely realistic, assumptions regarding the nature of eusocial evolution.     

Improving evolutionary models of protein interaction networks

MOTIVATION: Theoretical models of biological networks are valuable tools in evolutionary inference. Theoretical models based on gene duplication and divergence provide biologically plausible evolutionary mechanics. Similarities found between empirical networks and their theoretically generated counterpart are considered evidence of the role modeled mechanics play in biological evolution. However, the method by which these models are parameterized can lead to questions about the validity of the inferences. Selecting parameter values in order to produce a particular topological value obfuscates the possibility that the model may produce a similar topology for a large range of parameter values. Alternately, a model may produce a large range of topologies, allowing (incorrect) parameter values to produce a valid topology from an otherwise flawed model. In order to lend biological credence to the modeled evolutionary mechanics, parameter values should be derived from the empirical data. Furthermore, recent work indicates that the timing and fate of gene duplications are critical to proper derivation of these parameters. RESULTS: We present a methodology for deriving evolutionary rates from empirical data that is used to parameterize duplication and divergence models of protein interaction network evolution. Our method avoids shortcomings of previous methods, which failed to consider the effect of subsequent duplications. From our parameter values, we find that concurrent and existing existing duplication and divergence models are insufficient for modeling protein interaction network evolution. We introduce a model enhancement based on heritable interaction sites on the surface of a protein and find that it more closely reflects the high clustering found in the empirical network.     

The systematic position of Aspidytidae, the diversification of Dytiscoidea (Coleoptera, Adephaga) and the phylogenetic signal of third codon positions

Characters of the newly discovered larvae of the South African Cliff Water Beetle Aspidytes niobe were examined and integrated into a data matrix including all families of Dytiscoidea as well as Haliplidae. Fifty-three morphological characters of adults and larvae were analysed separately and combined with molecular data from six nuclear and mitochondrial genes. The phylogeny of the group is reconstructed for the study of the evolution of swimming behaviour and larval feeding habits, as well as the shift in diversification rates leading to the two most speciose lineages. The parsimony analysis of all equally weighted morphological and molecular characters combined resulted in a single well supported tree with the topology (Noteridae (Hygrobiidae ((Aspidytidae, Amphizoidae) Dytiscidae))), in agreement with the molecular data alone, but in contradiction to the morphological data, which favoured a topology in which Hygrobiidae is sister to Dytiscidae. The exclusion of third codon positions of the three protein coding genes resulted in a topology identical to that obtained with the morphological data alone, but the use of Bayesian probabilities or the amino acid sequence resulted in the same topology as that of the tree obtained with parsimony using all equally weighted characters. We concluded that interactions of third codon positions with the other data are complex, and their removal is not justified. There was a significant increase in the diversification rate at the base of the richest families (Noteridae and Dytiscidae), which could be associated with the development of simultaneous stroke and higher swimming performance, although data on the swimming behaviour of some basal groups of Noteridae are incomplete. The presence of larval mandibular sucking channels may have contributed to the diversification of Dytiscidae and the species-rich noterid genera Hydrocanthus and Canthydrus .     

Heaven it's my wife! Male canaries conceal extra-pair courtships but increase aggressions when their mate watches

Many animals live in a communication network, an environment where individuals can obtain information about competitors or potential mates by observing interactions between conspecifics. In such an environment, interactants might benefit by changing their signalling behaviour in the presence of an audience. This audience effect seems widespread among species, has been observed during various types of interaction (e.g. intra-sexual vs. inter-sexual interaction) and varies according to the social context (e.g. gender, hierarchical or mating status of the audience). However, the way individuals might adapt their signalling behaviour to a combination of these factors remains poorly understood. To address this question, we studied how the presence of an audience affects the behaviour of male domestic canaries Serinus canaria during two types of interactions: (i) an extra-pair interaction and (ii) a male-male competition for food. Males were observed under three conditions: (a) in the absence of audience, (b) in the presence of their mate or (c) of a familiar female. Our results show that male domestic canaries minutely adapt their courting and agonistic behaviours to a combination of: (i) the type of interaction (extra-pair interaction/male-male competition), (ii) the social context (mate, familiar female or nobody in audience) and (iii) the behaviours of both the audience and the interactant. These results highlight the ability of animals to subtly adapt their behaviour to the social environment. This also raises questions about the cognitive foundations and evolution of these processes especially considering that canaries are known neither for having high cognitive abilities nor for being a typical example for the social intelligence hypothesis.     

Modularity in the evolution of yeast protein interaction network

Protein interaction networks are known to exhibit remarkable structures: scale-free and small-world and modular structures. To explain the evolutionary processes of protein interaction networks possessing scale-free and small-world structures, preferential attachment and duplication-divergence models have been proposed as mathematical models. Protein interaction networks are also known to exhibit another remarkable structural characteristic, modular structure. How the protein interaction networks became to exhibit modularity in their evolution? Here, we propose a hypothesis of modularity in the evolution of yeast protein interaction network based on molecular evolutionary evidence. We assigned yeast proteins into six evolutionary ages by constructing a phylogenetic profile. We found that all the almost half of hub proteins are evolutionarily new. Examining the evolutionary processes of protein complexes, functional modules and topological modules, we also found that member proteins of these modules tend to appear in one or two evolutionary ages. Moreover, proteins in protein complexes and topological modules show significantly low evolutionary rates than those not in these modules. Our results suggest a hypothesis of modularity in the evolution of yeast protein interaction network as systems evolution.     

Questioning the Ubiquity of Neofunctionalization

Gene duplication provides much of the raw material from which functional diversity evolves. Two evolutionary mechanisms have been proposed that generate functional diversity: neofunctionalization, the de novo acquisition of function by one duplicate, and subfunctionalization, the partitioning of ancestral functions between gene duplicates. With protein interactions as a surrogate for protein functions, evidence of prodigious neofunctionalization and subfunctionalization has been identified in analyses of empirical protein interactions and evolutionary models of protein interactions. However, we have identified three phenomena that have contributed to neofunctionalization being erroneously identified as a significant factor in protein interaction network evolution. First, self-interacting proteins are underreported in interaction data due to biological artifacts and design limitations in the two most common high-throughput protein interaction assays. Second, evolutionary inferences have been drawn from paralog analysis without consideration for concurrent and subsequent duplication events. Third, the theoretical model of prodigious neofunctionalization is unable to reproduce empirical network clustering and relies on untenable parameter requirements. In light of these findings, we believe that protein interaction evolution is more persuasively characterized by subfunctionalization and self-interactions.     

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.     

Evolution of genes and taxa: a primer

The rapidly growing fields of molecular evolution and systematics have much to offer to molecular biology, but like any field have their own repertoire of terms and concepts. Homology, for example, is a central theme in evolutionary biology whose definition is complex and often controversial. Homology extends to multigene families, where the distinction between orthology and paralogy is key. Nucleotide sequence alignment is also a homology issue, and is a key stage in any evolutionary analysis of sequence data. Models based on our understanding of the processes of nucleotide substitution are used both in the estimation of the number of evolutionary changes between aligned sequences and in phylogeny reconstruction from sequence data. The three common methods of phylogeny reconstruction – parsimony, distance and maximum likelihood – differ in their use of these models. All three face similar problems in finding optimal – and reliable – solutions among the vast number of possible trees. Moreover, even optimal trees for a given gene may not reflect the relationships of the organisms from which the gene was sampled. Knowledge of how genes evolve and at what rate is critical for understanding gene function across species or within gene families. The Neutral Theory of Molecular Evolution serves as the null model of molecular evolution and plays a central role in data analysis. Three areas in which the Neutral Theory plays a vital role are: interpreting ratios of nonsynonymous to synonymous nucleotide substitutions, assessing the reliability of molecular clocks, and providing a foundation for molecular population genetics.     

Detecting phylogenetic signal in mutualistic interaction networks using a Markov process model

Ecological interaction networks, such as those describing the mutualistic interactions between plants and their pollinators or between plants and their frugivores, exhibit non‐random structural properties that cannot be explained by simple models of network formation. One factor affecting the formation and eventual structure of such a network is its evolutionary history. We argue that this, in many cases, is closely linked to the evolutionary histories of the species involved in the interactions. Indeed, empirical studies of interaction networks along with the phylogenies of the interacting species have demonstrated significant associations between phylogeny and network structure. To date, however, no generative model explaining the way in which the evolution of individual species affects the evolution of interaction networks has been proposed. We present a model describing the evolution of pairwise interactions as a branching Markov process, drawing on phylogenetic models of molecular evolution. Using knowledge of the phylogenies of the interacting species, our model yielded a significantly better fit to 21% of a set of plant–pollinator and plant–frugivore mutualistic networks. This highlights the importance, in a substantial minority of cases, of inheritance of interaction patterns without excluding the potential role of ecological novelties in forming the current network architecture. We suggest that our model can be used as a null model for controlling evolutionary signals when evaluating the role of other factors in shaping the emergence of ecological networks.     

Nuttalliella namaqua: a living fossil and closest relative to the ancestral tick lineage: implications for the evolution of blood-feeding in ticks

Ticks are monophyletic and composed of the hard (Ixodidae) and soft (Argasidae) tick families, as well as the Nuttalliellidae, a family with a single species, Nuttalliella namaqua. Significant biological differences in lifestyle strategies for hard and soft ticks suggest that various blood-feeding adaptations occurred after their divergence. The phylogenetic relationships between the tick families have not yet been resolved due to the lack of molecular data for N. namaqua. This tick possesses a pseudo-scutum and apical gnathostoma as observed for ixodids, has a leathery cuticle similar to argasids and has been considered the evolutionary missing link between the two families. Little knowledge exists with regard to its feeding biology or host preferences. Data on its biology and systematic relationship to the other tick families could therefore be crucial in understanding the evolution of blood-feeding behaviour in ticks. Live specimens were collected and blood meal analysis showed the presence of DNA for girdled lizards from the Cordylid family. Feeding of ticks on lizards showed that engorgement occurred rapidly, similar to argasids, but that blood meal concentration occurs via malpighian excretion of water. Phylogenetic analysis of the 18S nuclear and 16S mitochondrial genes indicate that N. namaqua grouped basal to the main tick families. The data supports the monophyly of all tick families and suggests the evolution of argasid-like blood-feeding behaviour in the ancestral tick lineage. Based on the data and considerations from literature we propose an origin for ticks in the Karoo basin of Gondwanaland during the late Permian. The nuttalliellid family almost became extinct during the End Permian event, leaving N. namaqua as the closest living relative to the ancestral tick lineage and the evolutionary missing link between the tick families.     

Recent and simultaneous origins of eusociality in halictid bees

Eusocial organisms are characterized by cooperative brood care, generation overlap and reproductive division of labour. Traits associated with eusociality are most developed in ants, termites, paper wasps and corbiculate bees; the fossil record indicates that each of these advanced eusocial taxa evolved in the Late Cretaceous or earlier (greater than 65 Myr ago). Halictid bees also include a large and diverse number of eusocial members, but, in contrast to advanced eusocial taxa, they are characterized by substantial intra- and inter-specific variation in social behaviour, which may be indicative of more recent eusocial evolution. To test this hypothesis, we used over 2400 bp of DNA sequence data gathered from three protein-coding nuclear genes (opsin, wingless and EF-1a) to infer the phylogeny of eusocial halictid lineages and their relatives. Results from relaxed molecular clock dating techniques that utilize a combination of molecular and fossil data indicate that the three independent origins of eusociality in halictid bees occurred within a narrow time frame between approximately 20 and 22 Myr ago. This relatively recent evolution helps to explain the pronounced levels of social variation observed within these bees. The three origins of eusociality appear to be temporally correlated with a period of global warming, suggesting that climate may have had an important role in the evolution and maintenance of eusociality in these bees.     

Domain combinations in archaeal, eubacterial and eukaryotic proteomes.

There is a limited repertoire of domain families that are duplicated and combined in different ways to form the set of proteins in a genome. Proteins are gene products, and at the level of genes, duplication, recombination, fusion and fission are the processes that produce new genes. We attempt to gain an overview of these processes by studying the evolutionary units in proteins, domains, in the protein sequences of 40 genomes. The domain and superfamily definitions in the Structural Classification of Proteins Database are used, so that we can view all pairs of adjacent domains in genome sequences in terms of their superfamily combinations. We find 783 out of the 859 superfamilies in SCOP in these genomes, and the 783 families occur in 1307 pairwise combinations. Most families are observed in combination with one or two other families, while a few families are very versatile in their combinatorial behaviour; 209 families do not make combinations with other families. This type of pattern can be described as a scale-free network. We also study the N to C-terminal orientation of domain pairs and domain repeats. The phylogenetic distribution of domain combinations is surveyed, to establish the extent of common and kingdom-specific combinations. Of the kingdom-specific combinations, significantly more combinations consist of families present in all three kingdoms than of families present in one or two kingdoms. Hence, we are led to conclude that recombination between common families, as compared to the invention of new families and recombination among these, has also been a major contribution to the evolution of kingdom-specific and species-specific functions in organisms in all three kingdoms. Finally, we compare the set of the domain combinations in the genomes to those in the RCSB Protein Data Bank, and discuss the implications for structural genomics.     

Protein folds, functions and evolution.

The evolution of proteins and their functions is reviewed from a structural perspective in the light of the current database. Protein domain families segregate unequally between the three major classes, the 32 different architectures and almost 700 folds observed to date. We find that the number of new topologies is still increasing, although 25 new structures are now determined for each new topology. The corresponding analysis and classification of function is only just beginning, fuelled by the genome data. The structural data revealed unexpected conservations and divergence of function both within and between families. The next five years will see the compilation of a definitive dictionary of protein families and their related functions, based on structural data which reveals relationships hidden at the sequence level. Such information will provide the foundation to build a better understanding of the molecular basis of biological complexity and hopefully to facilitate rational molecular design.     

Role of gene duplication in evolution

It is now known that many multigene and supergene families exist in eukaryote genomes: multigene families with uniform copy members like genes for ribosomal RNA, those with variable members like immunoglobulin genes, and supergene families such as those for various growth factor and hormone receptors. Many such examples indicate that gene duplication and subsequent differentiation are extremely important for organismal evolution. In particular, gene duplication could well have been the primary mechanism for the evolution of complexity in higher organisms. Population genetic models for the origin of gene families with diverse functions are presented, in which natural selection favors those genomes with more useful mutants in duplicated genes. Since any gene has a certain probability of degenerating by mutation, success versus failure in acquiring a new gene by duplication may be expressed as the ratio of probabilities of spreading of useful versus detrimental mutations in redundant gene copies. Also examined are the effects of gene duplication on evolution by compensatory advantageous mutations. Results of the analyses show that both natural selection and random drift are important for the origin of gene families. In addition, interaction between molecular mechanisms such as unequal crossing-over and gene conversion, and selection or drift is found to have a large effect on evolution by gene duplication.     

Influence of metabolic network structure and function on enzyme evolution

Background  

Most studies of molecular evolution are focused on individual genes and proteins. However, understanding the design principles and evolutionary properties of molecular networks requires a system-wide perspective. In the present work we connect molecular evolution on the gene level with system properties of a cellular metabolic network. In contrast to protein interaction networks, where several previous studies investigated the molecular evolution of proteins, metabolic networks have a relatively well-defined global function. The ability to consider fluxes in a metabolic network allows us to relate the functional role of each enzyme in a network to its rate of evolution.     

Multidimensional vector space representation for convergent evolution and molecular phylogeny

With growing amounts of genome data and constant improvement of models of molecular evolution, phylogenetic reconstruction became more reliable. However, our knowledge of the real process of molecular evolution is still limited. When enough large-sized data sets are analyzed, any subtle biases in statistical models can support incorrect topologies significantly because of the high signal-to-noise ratio. We propose a procedure to locate sequences in a multidimensional vector space (MVS), in which the geometry of the space is uniquely determined in such a way that the vectors of sequence evolution are orthogonal among different branches. In this paper, the MVS approach is developed to detect and remove biases in models of molecular evolution caused by unrecognized convergent evolution among lineages or unexpected patterns of substitutions. Biases in the estimated pairwise distances are identified as deviations (outliers) of sequence spatial vectors from the expected orthogonality. Modifications to the estimated distances are made by minimizing an index to quantify the deviations. In this way, it becomes possible to reconstruct the phylogenetic tree, taking account of possible biases in the model of molecular evolution. The efficacy of the modification procedure was verified by simulating evolution on various topologies with rate heterogeneity and convergent change. The phylogeny of placental mammals in previous analyses of large data sets has varied according to the genes being analyzed. Systematic deviations caused by convergent evolution were detected by our procedure in all representative data sets and were found to strongly affect the tree structure. However, the bias correction yielded a consistent topology among data sets. The existence of strong biases was validated by examining the sites of convergent evolution between the hedgehog and other species in mitochondrial data set. This convergent evolution explains why it has been difficult to determine the phylogenetic placement of the hedgehog in previous studies.     

Conserved flanking microsatellite sequences (ReFS) differentiate between Lepidoptera species,and provide insight into microsatellite evolution

High rates of mutation and homoplasy mean that microsatellites generally are not considered to be useful molecular markers for inferring systematic relationships between species. However, an earlier pilot study suggested that conserved flanking microsatellite sequences, also known as repetitive flanking sequences (ReFS), may form a basis for a dominant marker that can differentiate between species of Lepidoptera. We present data that demonstrate that ReFS are quick and easy to use, and generate highly repeatable banding patterns from a range of Lepidoptera species. Sequence data from a subset of ReFS‐amplified bands revealed microsatellite families with flanking sequences that are more conserved within than among species: this is probably attributable to recombination‐mediated events, transposition of mobile elements or a combination of the two. Our data support the use of ReFS as dominant interspecific molecular markers, and add to the growing literature on the evolution of microsatellites in Lepidoptera.