Search

Comparative genomics combined with phylogenetic reconstructions are powerful approaches to study the evolution of genes and genomes. However, the current rapid expansion of the volume of genomic information makes it increasingly difficult to interrogate, integrate and synthesize comparative genome data while taking into account the maximum breadth of information available. GenomicusPlants (http://www.genomicus.biologie.ens.fr/genomicus-plants) is an extension of the Genomicus webserver that addresses this issue by allowing users to explore flowering plant genomes in an intuitive way, across the broadest evolutionary scales. Extant genomes of 26 flowering plants can be analyzed, as well as 23 ancestral reconstructed genomes. Ancestral gene order provides a long-term chronological view of gene order evolution, greatly facilitating comparative genomics and evolutionary studies. Four main interfaces (‘views’) are available where: (i) PhyloView combines phylogenetic trees with comparisons of genomic loci across any number of genomes; (ii) AlignView projects loci of interest against all other genomes to visualize its topological conservation; (iii) MatrixView compares two genomes in a classical dotplot representation; and (iv) Karyoview visualizes chromosome karyotypes ‘painted’ with colours of another genome of interest. All four views are interconnected and benefit from many customizable features.

Background The persistent and growing gap between the availability of sequenced genomes and the ability to assign functions to sequenced genes led us to explore ways to maximize the information content of automated annotation for studies of anopheline mosquitos. Specifically, we use genome content analysis of a large number of previously sequenced anopheline mosquitos to follow the loss and gain of protein families over the evolutionary history of this group. The importance of this endeavor lies in the potential for comparative genomic studies between Anopheles and closely related non-vector species to reveal ancestral genome content dynamics involved in vector competence. In addition, comparisons within Anopheles could identify genome content changes responsible for variation in the vectorial capacity of this family of important parasite vectors. Results The competence and capacity of P. falciparum vectors do not appear to be phylogenetically constrained within the Anophelinae. Instead, using ancestral reconstruction methods, we suggest that a previously unexamined component of vector biology, anopheline nucleotide metabolism, may contribute to the unique status of anophelines as P. falciparum vectors. While the fitness effects of nucleotide co-option by P. falciparum parasites on their anopheline hosts are not yet known, our results suggest that anopheline genome content may be responding to selection pressure from P. falciparum. Whether this response is defensive, in an attempt to redress improper nucleotide balance resulting from P. falciparum infection, or perhaps symbiotic, resulting from an as-yet-unknown mutualism between anophelines and P. falciparum, is an open question that deserves further study. Conclusions Clearly, there is a wealth of functional information to be gained from detailed manual genome annotation, yet the rapid increase in the number of available sequences means that most researchers will not have the time or resources to manually annotate all the sequence data they generate. We believe that efforts to maximize the amount of information obtained from automated annotation can help address the functional annotation deficit that most evolutionary biologists now face, and here demonstrate the value of such an approach. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3590-0) contains supplementary material, which is available to authorized users.

Background Although resurrecting ancestral proteins is a powerful tool for understanding the molecular-functional evolution of gene families, nearly all studies have examined proteins functioning in relatively stable biological processes. The extent to which more dynamic systems obey the same ‘rules’ governing stable processes is unclear. Here we present the first detailed investigation of the functional evolution of the RIG-like receptors (RLRs), a family of innate immune receptors that detect viral RNA in the cytoplasm. Results Using kinetic binding assays and molecular dynamics simulations of ancestral proteins, we demonstrate how a small number of adaptive protein-coding changes repeatedly shifted the RNA preference of RLRs throughout animal evolution by reorganizing the shape and electrostatic distribution across the RNA binding pocket, altering the hydrogen bond network between the RLR and its RNA target. In contrast to observations of proteins involved in metabolism and development, we find that RLR-RNA preference ‘flip flopped’ between two functional states, and shifts in RNA preference were not always coupled to gene duplications or speciation events. We demonstrate at least one reversion of RLR-RNA preference from a derived to an ancestral function through a novel structural mechanism, indicating multiple structural implementations of similar functions. Conclusions Our results suggest a model in which frequent shifts in selection pressures imposed by an evolutionary arms race preclude the long-term functional optimization observed in stable biological systems. As a result, the evolutionary dynamics of immune receptors may be less constrained by structural epistasis and historical contingency. Electronic supplementary material The online version of this article (doi:10.1186/s12862-016-0818-6) contains supplementary material, which is available to authorized users.

The problem of reconstructing ancestral genomes in a given phylogenetic tree arises in many different comparative genomics fields. Here, we focus on reconstructing the gene order of ancestral genomes, a problem that has been largely studied in the past 20 years, especially with the increasing availability of whole genome DNA sequences. There are two main approaches to this problem: event-based methods, that try to find the ancestral genomes that minimize the number of rearrangement events in the tree; and homology-based, that look for conserved structures, such as adjacent genes in the extant genomes, to build the ancestral genomes. We propose algorithms that use the concept of intermediate genomes, arising in optimal pairwise rearrangement scenarios. We show that intermediate genomes have combinatorial properties that make them easy to reconstruct, and develop fast algorithms with better reconstructed ancestral genomes than current event-based methods. The proposed framework is also designed to accept extra information, such as results from homology-based approaches, giving rise to combined algorithms with better results than the original methods.

To reconstruct the evolution history of DNA sequences, novel models of increasing complexity regarding the number of free parameters taken into account in the sequence evolution, as well as faster and more accurate algorithms, and statistical and computational methods, are needed. More particularly, as the principal forces that have led to major structural changes are genome rearrangements (such as translocations, fusions, and so on), understanding their underlying mechanisms, among other things via the ancestral genome reconstruction, are essential. In this problem, since finding the ancestral genomes that minimize the number of rearrangements in a phylogenetic tree is known to be NP-hard for three or more genomes, heuristics are commonly chosen to obtain approximations of the exact solution. The aim of this work is to show that another path is possible. Various algorithms and software already deal with the difficult nature of the problem of reconstruction of the ancestral genome, but they do not function with precision, in particular when indels or single nucleotide polymorphisms fall into repeated sequences. In this article, and despite the theoretical NP-hardness of the ancestral reconstruction problem, we show that an exact solution can be found in practice in various cases, encompassing organelles and some bacteria. A practical example proves that an accurate reconstruction, which also allows to highlight homoplasic events, can be obtained. This is illustrated by the reconstruction of ancestral genomes of two bacterial pathogens, belonging in Mycobacterium and Brucella genera. By putting together automatically reconstructed ancestral regions with handmade ones for problematic cases, we show that an accurate reconstruction of ancestors of the Brucella genus and of the Mycobacterium tuberculosis complex is possible. By doing so, we are able to investigate the evolutionary history of each pathogen by computing their common ancestors. They can be investigated extensively, by studying the gene content evolution over time, the resistance acquisition, and the impacts of mobile elements on genome plasticity.

Background Invasive species can have devastating effects on native ecosystems and therefore impose a significant threat to human welfare. The introduction rate of invasive species has accelerated dramatically in recent times due to human activity (anthropogenic effects), with a steadily growing pool of widespread tramp species. We present an in-depth analysis of four pantropical species of Xyleborus ambrosia beetles (Xyleborus volvulus, Xyleborus perforans, Xyleborus ferrugineus, and Xyleborus affinis) with similar ecology (fungus cultivation in dead wood), reproductive biology (permanent inbreeding) and genetic system (haplodiploidy). The unique combination of reproductive traits and broad host plant usage pre-adapts these beetles for colonizing of new areas. Results We found that all four species were broadly distributed long before human-assisted dispersal became common, and that the impact of anthropogenic effects varied among the species. For X. volvulus, X. perforans, and X. affinis there was evidence of ancient establishment in numerous regions, but also of abundant recent introductions into previously colonized areas. For X. ferrugineus, we found clear biogeographical structuring of old clades, but little evidence for recent successful introductions. Conclusions Our results indicate that current human-aided transoceanic dispersal has strongly affected the genetic makeup of three of the species in this study. However, current biogeographical patterns of all four species are equally, if not more strongly, influenced by ancient establishment on different continents. Electronic supplementary material The online version of this article (doi:10.1186/s12862-016-0610-7) contains supplementary material, which is available to authorized users.

During evolution, genomes are modified by large scale structural events, such as rearrangements, deletions or insertions of large blocks of DNA. Of particular interest, in order to better understand how this type of genomic evolution happens, is the reconstruction of ancestral genomes, given a phylogenetic tree with extant genomes at its leaves. One way of solving this problem is to assume a rearrangement model, such as Double Cut and Join (DCJ), and find a set of ancestral genomes that minimizes the number of events on the input tree. Since this problem is NP-hard for most rearrangement models, exact solutions are practical only for small instances, and heuristics have to be used for larger datasets. This type of approach can be called event-based. Another common approach is based on finding conserved structures between the input genomes, such as adjacencies between genes, possibly also assigning weights that indicate a measure of confidence or probability that this particular structure is present on each ancestral genome, and then finding a set of non conflicting adjacencies that optimize some given function, usually trying to maximize total weight and minimizing character changes in the tree. We call this type of methods homology-based. In previous work, we proposed an ancestral reconstruction method that combines homology- and event-based ideas, using the concept of intermediate genomes, that arise in DCJ rearrangement scenarios. This method showed better rate of correctly reconstructed adjacencies than other methods, while also being faster, since the use of intermediate genomes greatly reduces the search space. Here, we generalize the intermediate genome concept to genomes with unequal gene content, extending our method to account for gene insertions and deletions of any length. In many of the simulated datasets, our proposed method had better results than MLGO and MGRA, two state-of-the-art algorithms for ancestral reconstruction with unequal gene content, while running much faster, making it more scalable to larger datasets. Studing ancestral reconstruction problems under a new light, using the concept of intermediate genomes, allows the design of very fast algorithms by greatly reducing the solution search space, while also giving very good results. The algorithms introduced in this paper were implemented in an open-source software called RINGO (ancestral Reconstruction with INtermediate GenOmes), available at https://github.com/pedrofeijao/RINGO .

Background Geckos of the genus Pachydactylus and their close relatives comprise the most species-rich clade of lizards in sub-Saharan Africa. Many explanations have been offered to explain species richness patterns of clades. In the Pachydactylus group, one possible explanation is a history of diversification via geographic isolation. If geographic isolation has played a key role in facilitating diversification, then we expect species in more species-rich subclades to have smaller ranges than species in less diverse subclades. We also expect traits promoting geographic isolation to be correlated with small geographic ranges. In order to test these expectations, we performed phylogenetic analyses and tested for correlations among body size, habitat choice, range sizes, and diversification rates in the Pachydactylus group. Results Both body size and habitat use are inferred to have shifted multiple times across the phylogeny of the Pachydactylus group, with large size and generalist habitat use being ancestral for the group. Geographic range size is correlated with both of these traits. Small-bodied species have more restricted ranges than large-bodied species, and rock-dwelling species have more restricted ranges than either terrestrial or generalist species. Rock-dwelling and small body size are also associated with higher rates of diversification, and subclades retaining ancestral conditions for these traits are less species rich than subclades in which shifts to small body size and rocky habitat use have occurred. The phylogeny also illustrates inadequacies of the current taxonomy of the group. Conclusions The results are consistent with a model in which lineages more likely to become geographically isolated diversify to a greater extent, although some patterns also resemble those expected of an adaptive radiation in which ecological divergence acts as a driver of speciation. Therefore, the Pachydactylus group may represent an intermediate between clades in which radiation is adaptive versus those in which it is non-adaptive.

by other methods of ancestral reconstruction. Critics also argue that it is not possible to determine which of several alternative topologies, differing by a small number of NR's, most closely reflects cladogeny since the statistical significance of the NR score differences is unknown. A careful investigation of the method of these critics reveals that their method is nonrigorous, nonreplicable and produces erroneous results. The need for rigorous, automated systems of data analysis in studies of molecular evolution is made apparent, particularly in light of the voluminous nature of the data. We here describe a series of parsimony tests conducted with globin amino acid sequence data, focusing on the position of Tarsius within Primates, in order to illustrate the importance of interpreting NR score differences within a context which includes the number of lineages represented in a data set, the density of the lineages in particular regions of a dendrogram and the exact nature of the changes attributed to specific NR savings. Analyzed within the context of these changing parameters, small NR score differences between alternative topologies can be meaningful and do produce results with significant confidence. The phylogenetic relationship of Tarsius to Anthropoidea, as well as other important issues in primate systematics, are expected to be greatly clarified with the gathering and rigorous analysis of nucleotide sequence data which can resolve the difficulties presented by convergent evolution and paralogous gene lines. [Tarsier phylogeny; maximum parsimony methodology; molecular evolution; primate evolution; molecular anthropology; hemoglobin phylogeny.]

Background Because prokaryotic genomes experience a rapid flux of genes, selection may act at a higher level than an individual genome. We explore a quantitative model of the distributed genome whereby groups of genomes evolve by acquiring genes from a fixed reservoir which we denote as supergenome. Previous attempts to understand the nature of the supergenome treated genomes as random, independent collections of genes and assumed that the supergenome consists of a small number of homogeneous sub-reservoirs. Here we explore the consequences of relaxing both assumptions. Results We surveyed several methods for estimating the size and composition of the supergenome. The methods assumed that genomes were either random, independent samples of the supergenome or that they evolved from a common ancestor along a known tree via stochastic sampling from the reservoir. The reservoir was assumed to be either a collection of homogeneous sub-reservoirs or alternatively composed of genes with Gamma distributed gain probabilities. Empirical gene frequencies were used to either compute the likelihood of the data directly or first to reconstruct the history of gene gains and then compute the likelihood of the reconstructed numbers of gains. Conclusions Supergenome size estimates using the empirical gene frequencies directly are not robust with respect to the choice of the model. By contrast, using the gene frequencies and the phylogenetic tree to reconstruct multiple gene gains produces reliable estimates of the supergenome size and indicates that a homogeneous supergenome is more consistent with the data than a supergenome with Gamma distributed gain probabilities.

Background The problem of reconstructing ancestral genomes in a given phylogenetic tree arises in many different comparative genomics fields. Here, we focus on reconstructing the gene order of ancestral genomes, a problem that has been largely studied in the past 20 years, especially with the increasing availability of whole genome DNA sequences. There are two main approaches to this problem: event-based methods, that try to find the ancestral genomes that minimize the number of rearrangement events in the tree; and homology-based, that look for conserved structures, such as adjacent genes in the extant genomes, to build the ancestral genomes. Results We propose algorithms that use the concept of intermediate genomes, arising in optimal pairwise rearrangement scenarios. We show that intermediate genomes have combinatorial properties that make them easy to reconstruct, and develop fast algorithms with better reconstructed ancestral genomes than current event-based methods. The proposed framework is also designed to accept extra information, such as results from homology-based approaches, giving rise to combined algorithms with better results than the original methods.

Introduction Many mammals have evolved highly adapted hearing associated with ecological specialisation. Of these, bats possess the widest frequency range of vocalisations and associated hearing sensitivities, with frequencies of above 200 kHz in some lineages that use laryngeal echolocation. High frequency hearing in bats appears to have evolved via structural modifications of the inner ear, however, studying these minute features presents considerable challenges and hitherto few such attempts have been made. To understand these adaptations more fully, as well as gain insights into the evolutionary origins of ultrasonic hearing and echolocation in bats, we undertook micro-computed tomography (μCT) scans of the cochleae of representative bat species from 16 families, encompassing their broad range of ecological diversity. To characterise cochlear gross morphology, we measured the relative basilar membrane length and number of turns, and compared these values between echolocating and non-echolocating bats, as well as other mammals. Results We found that hearing and echolocation call frequencies in bats correlated with both measures of cochlear morphology. In particular, relative basilar membrane length was typically longer in echolocating species, and also correlated positively with the number of cochlear turns. Ancestral reconstructions of these parameters suggested that the common ancestor of all extant bats was probably capable of ultrasonic hearing; however, we also found evidence of a significant decrease in the rate of morphological evolution of the basilar membrane in multiple ancestral branches within the Yangochiroptera suborder. Within the echolocating Yinpterochiroptera, there was some evidence of an increase in the rate of basilar membrane evolution in some tips of the tree, possibly associated with reported shifts in call frequency associated with recent speciation events. Conclusions The two main groups of echolocating bat were found to display highly variable inner ear morphologies. Ancestral reconstructions and rate shift analyses of ear morphology point to a complex evolutionary history, with the former supporting ultrasonic hearing in the common bat ancestor but the latter suggesting that morphological changes associated with echolocation might have occurred later. These findings are consistent with theories that sophisticated laryngeal echolocation, as seen in modern lineages, evolved following the divergence of the two main suborders.