Inferring ancestral sequences in taxon-rich phylogenies

Statistical consistency in phylogenetics has traditionally referred to the accuracy of estimating phylogenetic parameters for a fixed number of species as we increase the number of characters. However, it is also useful to consider a dual type of statistical consistency where we increase the number of species, rather than characters. This raises some basic questions: what can we learn about the evolutionary process as we increase the number of species? In particular, does having more species allow us to infer the ancestral state of characters accurately? This question is particularly important when sequence evolution varies in a complex way from character to character, as methods applicable for i.i.d. models may no longer be valid. In this paper, we assemble a collection of results to analyse various approaches for inferring ancestral information with increasing accuracy as the number of taxa increases.     

Inferring species-level phylogenies and taxonomic distinctiveness using multilocus data in Sistrurus rattlesnakes

Phylogenetic relationships and taxonomic distinctiveness of closely related species and subspecies are most accurately inferred from data derived from multiple independent loci. Here, we apply several approaches for understanding species-level relationships using data from 18 nuclear DNA loci and 1 mitochondrial DNA locus within currently described species and subspecies of Sistrurus rattlesnakes. Collectively, these methods provide evidence that a currently described species, the massasauga rattlesnake (Sistrurus catenatus), consists of two well-supported clades, one composed of the two western subspecies (S. c. tergeminus and S. c. edwardsii) and the other the eastern subspecies (S. c. catenatus). Within pigmy rattlesnakes (S. miliarius), however, there is not strong support across methods for any particular grouping at the subspecific level. Monophyly based tests for taxonomic distinctiveness show evidence for distinctiveness of all subspecies but this support is strongest by far for the S. c. catenatus clade. Because support for the distinctiveness of S. c. catenatus is both strong and consistent across methods, and due to its morphological distinctiveness and allopatric distribution, we suggest that this subspecies be elevated to full species status, which has significant conservation implications. Finally, most divergence time estimates based upon a fossil-calibrated species tree are > 50% younger than those from a concatenated gene tree analysis and suggest that an active period of speciation within Sistrurus occurred within the late Pliocene/Pleistocene eras.     

Calculating structural complexity in phylogenies using ancestral ontologies

Complexity is an important aspect of evolutionary biology, but there are many reasonable concepts of complexity, and its objective measurement is an elusive matter. Here we develop a simple measure of complexity based on counts of elements, incorporating the hierarchical information as represented in anatomical ontologies. Neomorphic and transformational characters are used to identify novelties and individuated morphological regions, respectively. By linking the characters to terms in an anatomical ontology a node‐driven approach is implemented, where a node ontology and a complexity score are inferred from the optimization of individual characters on each ancestral or terminal node. From the atomized vector of character scorings, the anatomical ontology is used to integrate the hierarchical structure of morphology in terminals and ancestors. These node ontologies are used to calculate a measure of complexity that can be traced on phylogenetic trees and is harmonious with usual phylogenetic operations. This strategy is compared with a terminal‐driven approach, in which the complexity scores are calculated only for terminals, and optimized as a continuous character on the internal nodes. These ideas are applied to a real dataset of 166 araneomorph spider species scored for 393 characters, using Spider Ontology (SPD, https://bioportal.bioontology.org/ontologies/SPD ); complexity scores and transitions are calculated for each node and branch, respectively. This result in a distribution of transitions skewed towards simplification; the transitions in complexity have no apparent correlation with character branch lengths. The node‐driven and terminal‐driven estimations are generally correlated in the complexity scores, but have higher divergence in the transition values. The structure of the ontology is used to provide complexity scores for organ systems and body parts of the focal groups.     

Bayesian estimation of ancestral character states on phylogenies

Biologists frequently attempt to infer the character states at ancestral nodes of a phylogeny from the distribution of traits observed in contemporary organisms. Because phylogenies are normally inferences from data, it is desirable to account for the uncertainty in estimates of the tree and its branch lengths when making inferences about ancestral states or other comparative parameters. Here we present a general Bayesian approach for testing comparative hypotheses across statistically justified samples of phylogenies, focusing on the specific issue of reconstructing ancestral states. The method uses Markov chain Monte Carlo techniques for sampling phylogenetic trees and for investigating the parameters of a statistical model of trait evolution. We describe how to combine information about the uncertainty of the phylogeny with uncertainty in the estimate of the ancestral state. Our approach does not constrain the sample of trees only to those that contain the ancestral node or nodes of interest, and we show how to reconstruct ancestral states of uncertain nodes using a most-recent-common-ancestor approach. We illustrate the methods with data on ribonuclease evolution in the Artiodactyla. Software implementing the methods (BayesMultiState) is available from the authors.     

Molecular phylogeny of coleoid cephalopods (Mollusca: Cephalopoda) using a multigene approach; the effect of data partitioning on resolving phylogenies in a Bayesian framework

The resolution of higher level phylogeny of the coleoid cephalopods (octopuses, squids, and cuttlefishes) has been hindered by homoplasy among morphological characters in conjunction with a very poor fossil record. Initial molecular studies, based primarily on small fragments of single mitochondrial genes, have produced little resolution of the deep relationships amongst coleoid cephalopod families. The present study investigated this issue using 3415 base pairs (bp) from three nuclear genes (octopine dehydrogenase, pax-6, and rhodopsin) and three mitochondrial genes (12S rDNA, 16S rDNA, and cytochrome oxidase I) from a total of 35 species (including representatives of each of the higher level taxa). Bayesian analyses were conducted on mitochondrial and nuclear genes separately and also all six genes together. Separate analyses were conducted with the data partitioned by gene, codon/rDNA, gene+codon/rDNA or not partitioned at all. In the majority of analyses partitioning the data by gene+codon was the appropriate model with partitioning by codon the second most selected model. In some instances the topology varied according to the model used. Relatively high posterior probabilities and high levels of congruence were present between the topologies resulting from the analysis of all Octopodiform (octopuses and vampire "squid") taxa for all six genes, and independently for the datasets of mitochondrial and nuclear genes. In contrast, the highest levels of resolution within the Decapodiformes (squids and cuttlefishes) resulted from analysis of nuclear genes alone. Different higher level Decapodiform topologies were obtained through the analysis of only the 1st+2nd codon positions of nuclear genes and of all three codon positions. It is notable that there is strong evidence of saturation among the 3rd codon positions within the Decapodiformes and this may contribute spurious signal. The results suggest that the Decapodiformes may have radiated earlier and/or had faster rates of evolution than the Octopodiformes. The following taxonomic conclusions are drawn from our analyses: (1) the order Octopoda and suborders Cirrata, Incirrata, and Oegopsida are monophyletic groups; (2) the family Spirulidae (Ram's horn squids) are the sister taxon to the family Sepiidae (cuttlefishes); (3) the family Octopodidae, as currently defined, is paraphyletic; (4) the superfamily Argonautoidea are basal within the suborder Incirrata; and (5) the benthic octopus genera Benthoctopus and Enteroctopus are sister taxa.     

On the impossibility of reconstructing ancestral data and phylogenies.

We prove that it is impossible to reconstruct ancestral data at the root of "deep" phylogenetic trees with high mutation rates. Moreover, we prove that it is impossible to reconstruct the topology of "deep" trees with high mutation rates from a number of characters smaller than a low-degree polynomial in the number of leaves. Our impossibility results hold for all reconstruction methods. The proofs apply tools from information theory and percolation theory.     

Molecular phylogeny of Elephantidae. Extreme divergence of the extant forest African elephant.

A phylogenetic study of the Elephantidae (Proboscidea, Mammalia) is based on the cytochrome b mitochondrial gene: 31 terminals, that is, all known sequences, one non-elephantid proboscidean, the extinct American mastodon, and four outgroups. The data set includes 11 new sequences with the first published sequence of the forest African elephant, L. a. cyclotis. The analyses of extant taxa only and of both extant and extinct taxa show that L. a. cyclotis is highly divergent from L. a. africana. It is as divergent from L. a. africana as Loxodonta is divergent from Elephas. Southern L. a. africana form a clade. The continental subspecies E. m. indicus is paraphyletic with individuals from India and Thailand closer to E. m. maximus (Sri-Lanka). Members of Mammuthus primigenius are more closely related to Loxodonta although they do not form a clade; two specimens of M. primigenius are closer to L. a. africana making the genus Loxodonta paraphyletic. The latter conclusion may be partly due to unequal length of the various polymorphic mammoth sequences.     

被子植物系统发育深层关系研究: 进展与挑战

被子植物系统发育学是研究被子植物及其各类群间亲缘关系与进化历史的学科。从20世纪90年代起, 核苷酸和氨基酸序列等分子数据开始被广泛运用于被子植物系统发育研究, 经过20多年的发展, 从使用单个或联合少数几个细胞器基因, 到近期应用整个叶绿体基因组来重建被子植物的系统发育关系, 目、科水平上的被子植物系统发育框架已被广泛接受。在这个框架中, 基部类群、主要的5个分支(即真双子叶植物、单子叶植物、木兰类、金粟兰目和金鱼藻目)、每个分支所包含的目以及几个大分支包括的核心类群等都具有高度支持。与此同时, 细胞器基因还存在一些固有的问题, 例如单亲遗传、系统发育信息量有限等, 因此近年来双亲遗传的核基因在被子植物系统发育研究中的重要性逐渐得到关注, 并在不同分类阶元的研究中都取得了一定进展。但是, 被子植物系统发育中仍然存在一些难以确定的关系, 例如被子植物5个分支之间的关系、真双子叶植物内部某些类群的位置等。本文简述了20多年来被子植物系统发育深层关系的主要研究进展, 讨论了被子植物系统发育学常用的细胞器基因和核基因的选用, 已经确定和尚未确定系统发育位置的主要类群, 以及研究中尚存在的问题和可能的解决方法。     

Molecular phylogenies in angiosperm evolution

We have cloned and sequenced cDNAs for the glyceraldehyde-3-phosphate dehydrogenase of glycolysis, gapC, from a bryophyte, a gymnosperm, and three angiosperms. Phylogenetic analyses are presented for these data in the context of other gapC sequences and in parallel with published nucleotide sequences for the chloroplast encoded gene for the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (rbcL). Relative-rate tests were performed for these genes in order to assess variation in substitution rate for coding regions, along individual plant lineages studied. The results of both gene analyses suggest that the deepest dichotomy within the angiosperms separates not magnoliids from remaining angiosperms, but monocotyledons from dicotyledons, in sharp contrast to prediction from the Euanthial theory for angiosperm evolution. Furthermore, these chloroplast and nuclear sequence data taken together suggest that the separation of monocotyledonous and dicotyledonous lineages took place in late Carboniferous times [approximately 300 Myr before the present (Mybp)]. This date would exceed but be compatible with the late-Triassic (approximately 220 Mybp) occurrence of fossil reproductive structures of the primitive angiosperm Sanmiguelia lewisii.      

The phylogeny of Cetartiodactyla: the importance of dense taxon sampling, missing data, and the remarkable promise of cytochrome b to provide reliable species-level phylogenies

We perform Bayesian phylogenetic analyses on cytochrome b sequences from 264 of the 290 extant cetartiodactyl mammals (whales plus even-toed ungulates) and two recently extinct species, the 'Mouse Goat' and the 'Irish Elk'. Previous primary analyses have included only a small portion of the species diversity within Cetartiodactyla, while a complete supertree analysis lacks resolution and branch lengths limiting its utility for comparative studies. The benefits of using a single-gene approach include rapid phylogenetic estimates for a large number of species. However, single-gene phylogenies often differ dramatically from studies involving multiple datasets suggesting that they often are unreliable. However, based on recovery of benchmark clades-clades supported in prior studies based on multiple independent datasets-and recovery of undisputed traditional taxonomic groups, Cytb performs extraordinarily well in resolving cetartiodactyl phylogeny when taxon sampling is dense. Missing data, however, (taxa with partial sequences) can compromise phylogenetic accuracy, suggesting a tradeoff between the benefits of adding taxa and introducing question marks. In the full data, a few species with a short sequences appear misplaced, however, sequence length alone seems a poor predictor of this phenomenon as other taxa with equally short sequences were not conspicuously misplaced. Although we recommend awaiting a better supported phylogeny based on more character data to reconsider classification and taxonomy within Cetartiodactyla, the new phylogenetic hypotheses provided here represent the currently best available tool for comparative species-level studies within this group. Cytb has been sequenced for a large percentage of mammals and appears to be a reliable phylogenetic marker as long as taxon sampling is dense. Therefore, an opportunity exists now to reconstruct detailed phylogenies of most of the major mammalian clades to rapidly provide much needed tools for species-level comparative studies.     

Synergistic effects of combining morphological and molecular data in resolving the phylogeny of butterflies and skippers

Phylogenetic relationships among major clades of butterflies and skippers have long been controversial, with no general consensus even today. Such lack of resolution is a substantial impediment to using the otherwise well studied butterflies as a model group in biology. Here we report the results of a combined analysis of DNA sequences from three genes and a morphological data matrix for 57 taxa (3258 characters, 1290 parsimony informative) representing all major lineages from the three putative butterfly super-families (Hedyloidea, Hesperioidea and Papilionoidea), plus out-groups representing other ditrysian Lepidoptera families. Recently, the utility of morphological data as a source of phylogenetic evidence has been debated. We present the first well supported phylogenetic hypothesis for the butterflies and skippers based on a total-evidence analysis of both traditional morphological characters and new molecular characters from three gene regions (COI, EF-1alpha and wingless). All four data partitions show substantial hidden support for the deeper nodes, which emerges only in a combined analysis in which the addition of morphological data plays a crucial role. With the exception of Nymphalidae, the traditionally recognized families are found to be strongly supported monophyletic clades with the following relationships: (Hesperiidae+(Papilionidae+(Pieridae+(Nymphalidae+(Lycaenidae+Riodinidae))))). Nymphalidae is recovered as a monophyletic clade but this clade does not have strong support. Lycaenidae and Riodinidae are sister groups with strong support and we suggest that the latter be given family rank. The position of Pieridae as the sister taxon to nymphalids, lycaenids and riodinids is supported by morphology and the EF-1alpha data but conflicted by the COI and wingless data. Hedylidae are more likely to be related to butterflies and skippers than geometrid moths and appear to be the sister group to Papilionoidea+Hesperioidea.     

Molecular evidence for species-level distinctions in clouded leopards

Among the 37 living species of Felidae, the clouded leopard (Neofelis nebulosa) is generally classified as a monotypic genus basal to the Panthera lineage of great cats. This secretive, mid-sized (16-23 kg) carnivore, now severely endangered, is traditionally subdivided into four southeast Asian subspecies (Figure 1A). We used molecular genetic methods to re-evaluate subspecies partitions and to quantify patterns of population genetic variation among 109 clouded leopards of known geographic origin (Figure 1A, Tables S1 ans S2 in the Supplemental Data available online). We found strong phylogeographic monophyly and large genetic distances between N. n. nebulosa (mainland) and N. n. diardi (Borneo; n = 3 individuals) with mtDNA (771 bp), nuclear DNA (3100 bp), and 51 microsatellite loci. Thirty-six fixed mitochondrial and nuclear nucleotide differences and 20 microsatellite loci with nonoverlapping allele-size ranges distinguished N. n. nebulosa from N. n. diardi. Along with fixed subspecies-specific chromosomal differences, this degree of differentiation is equivalent to, or greater than, comparable measures among five recognized Panthera species (lion, tiger, leopard, jaguar, and snow leopard). These distinctions increase the urgency of clouded leopard conservation efforts, and if affirmed by morphological analysis and wider sampling of N. n. diardi in Borneo and Sumatra, would support reclassification of N. n. diardi as a new species (Neofelis diardi).     

Dinoflagellate phylogeny revisited: reconciling morphological and molecular based phylogenies

Ultrastructural and molecular phylogenetic data suggest that dinoflagellates diverged as a lineage possibly as early as the Precambrian. However, the fossil record is problematic before the Mesozoic. From the mid Triassic, though, the fossil record of dinoflagellates is a rich source of information on Mesozoic-Cenozoic dinoflagellates, especially the gonyaulacoids and peridinioids. From the sequence of appearance of species and tabulation types and the impression of early morphological experimentation and later stabilization, the early Mesozoic radiation of dinoflagellates appears to be a real evolutionary event: indeed, dinoflagellate morphology as we know it today may originate in that event. This would explain why it is so difficult to interpret earlier fossils as dinoflagellates. However, that the dinoflagellate lineage existed in some form in the pre-Mesozoic is supported by biogeochemical data, early results of which indicate that certain early Paleozoic acanthomorph acritarchs may belong to the lineage.

A surprising degree of consistency is observed between ultrastructural (including tabulational), coarse biochemical and molecular sequence data. For example, sequence data provided by small subunit (SSU) rRNA support the hypothesis of progressive loss of histones within the dinoflagellates. Gymnodinioids have long been considered to be polyphyletic but are thought of generally as forerunners to the strongly thecate groups such as gonyaulacoids and peridinioids. In molecular trees they appear in both early-derived and late-derived positions, but mostly the latter. SSU data clearly support the gonyaulacoid/peridinioid ordinal separation, as does the fossil record. Prorocentroids are now thought to be the among the most derived dinoflagellates (and presumably the morphologically similar dinophysoids), but SSU sequences have so far failed to resolve the relationships of most gymnodinioids, peridinioids and prorocentroids (the so-called GPP complex) to one another. However, they do suggest the origin of prorocentroids from peridinioids rather than gonyaulacoids and that gymnodinioids probably had several origins.     

Molecular phylogeny of extant horseshoe crabs (Xiphosura, Limulidae) indicates Paleogene diversification of Asian species

Horseshoe crabs are marine invertebrates well known for their exceptionally low rates of diversification during their entire evolutionary history. Despite the low species diversity in the group, the phylogenetic relationships among the extant species, especially among the three Asian species are still unresolved. Here we apply a new set of molecular genetic data in combination with a wide geographic sampling of the intra-specific diversity to reinvestigate the evolutionary history among the four living limulid xiphosurans. Our analysis of the intraspecific diversity reveals low levels of connectivity among Carcinoscorpius rotundicauda lineages, which can be explained by the estuarine-bound ecology of this species. Moreover, a clear genetic break across the Thai-Malay Peninsula suggests the presence of cryptic species in C. rotundicauda. The limulid phylogeny finds strong support for a monophyletic genus Tachypleus and a diversification of the three Asian species during the Paleogene period, with speciation events well separated in time by several million years. The tree topology suggests that the three Asian species originated in central South East Asia from a marine stem group that inhabited the shallow coastal waters between the Andaman Sea, Vietnam, and Borneo. In this region C. rotundicauda probably separated from the Tachypleus stem group by invading estuarine habitats, while Tachypleus tridentatus most likely migrated northeast along the Southern coast of China and towards Japan.     

Dogs,cats, and kin: A molecular species-level phylogeny of Carnivora

Phylogenies underpin comparative biology as high-utility tools to test evolutionary and biogeographic hypotheses, inform on conservation strategies, and reveal the age and evolutionary histories of traits and lineages. As tools, most powerful are those phylogenies that contain all, or nearly all, of the taxa of a given group. Despite their obvious utility, such phylogenies, other than summary ‘supertrees’, are currently lacking for most mammalian orders, including the order Carnivora. Carnivora consists of about 270 extant species including most of the world’s large terrestrial predators (e.g., the big cats, wolves, bears), as well as many of man’s favorite wild (panda, cheetah, tiger) and domesticated animals (dog, cat). Distributed globally, carnivores are highly diverse ecologically, having occupied all major habitat types on the planet and being diverse in traits such as sociality, communication, body/brain size, and foraging ecology. Thus, numerous studies continue to address comparative questions within the order, highlighting the need for a detailed species-level phylogeny. Here we present a phylogeny of Carnivora that increases taxon sampling density from 28% in the most detailed primary-data study to date, to 82% containing 243 taxa (222 extant species, 17 subspecies). In addition to extant species, we sampled four extinct species: American cheetah, saber-toothed cat, cave bear and the giant short-faced bear. Bayesian analysis of cytochrome b sequences data-mined from GenBank results in a phylogenetic hypothesis that is largely congruent with prior studies based on fewer taxa but more characters. We find support for the monophyly of Carnivora, its major division into Caniformia and Feliformia, and for all but one family within the order. The only exception is the placement of the kinkajou outside Procyonidae, however, prior studies have already cast doubt on its family placement. In contrast, at the subfamily and genus level, our results indicate numerous problems with current classification. Our results also propose new, controversial hypotheses, such as the possible placement of the red panda (Ailuridae) sister to canids (Canidae). Our results confirm previous findings suggesting that the dog was domesticated from the Eurasian wolf (Canis lupus lupus) and are congruent with the Near East domestication of the cat. In sum, this study presents the most detailed species-level phylogeny of Carnivora to date and a much needed tool for comparative studies of carnivoran species. To demonstrate one such use, we perform a phylogenetic analysis of evolutionary distinctiveness (EDGE), which can be used to help establish conservation priorities. According with those criteria, and under one of the many possible sets of parameters, the highest priority Carnivora species for conservation of evolutionary diversity include: monk seals, giant and red panda, giant otter, otter civet, Owston’s palm civet, sea otter, Liberian mongoose, spectacled bear, walrus, binturong, and the fossa.     

Mitogenomics and phylogenomics reveal priapulid worms as extant models of the ancestral Ecdysozoan

Research into arthropod evolution is hampered by the derived nature and rapid evolution of the best-studied out-group: the nematodes. We consider priapulids as an alternative out-group. Priapulids are a small phylum of bottom-dwelling marine worms; their tubular body with spiny proboscis or introvert has changed little over 520 million years and recognizable priapulids are common among exceptionally preserved Cambrian fossils. Using the complete mitochondrial genome and 42 nuclear genes from Priapulus caudatus, we show that priapulids are slowly evolving ecdysozoans; almost all these priapulid genes have evolved more slowly than nematode orthologs and the priapulid mitochondrial gene order may be unchanged since the Cambrian. Considering their primitive bodyplan and embryology and the great conservation of both nuclear and mitochondrial genomes, priapulids may deserve the popular epithet of "living fossil." Their study is likely to yield significant new insights into the early evolution of the Ecdysozoa and the origins of the arthropods and their kin as well as aiding inference of the morphology of ancestral Ecdysozoa and Bilateria and their genomes.     

The effect of habitat type on speciation rates and range movements in aquatic beetles: inferences from species-level phylogenies

Most aquatic beetles in the family Dytiscidae are tightly associated either with running (lotic) or stagnant (lentic) water bodies. The range size of lotic species is known to be, on average, much smaller than that of lentic species, presumably as a result of differences in dispersal strategies in each habitat type. We explored possible effects of these differences on clade evolution and speciation rates by comparing species-level phylogenies based on cytochrome oxidase I (COI) and 16S rRNA mitochondrial genes for two genera, the lentic Ilybius and the lotic Deronectes. The expectation that species turnover is higher in lotic lineages due to their lower dispersal propensity compared to lentic species was not strongly supported. Deronectes displays a higher frequency of recent splits than Ilybius, consistent with the hypothesis, but the difference was not significant compared to expected patterns under a constant speciation rate null model. Similarly, when the degree of sympatry was plotted against relative node age, more allopatric splits were evident in the lentic Deronectes, suggesting a slower rate of range movement since speciation, but the differences were not significant. We discuss two explanations for our failure to detect differences between the two clades. First, current methods for analysing species-level phylogenies may be sensitive to taxonomic and sampling artefacts. Second, lentic and lotic clades may indeed display similar levels of species turnover despite occupying very different habitats at different spatial scales. More work is needed to investigate the effects of population level processes and spatial scale on macroevolutionary dynamics.     

Mitochondrial phylogeny of extant Hawaiian tree snails (Achatinellinae)

Hawaiian tree snails in the endemic subfamily Achatinellinae display a staggering variety of shell colors and banding patterns. Despite numerous attempts to classify this morphological variation, a conclusive phylogeny has not been proposed. To improve conservation efforts, we sought to better understand the species identities and phylogenetic relationships among the extant species of Achatinella and Partulina using partial mitochondrial 16S ribosomal DNA sequences. The reconstructed phylogeny showed a high degree of support for more recent branches, but gave little support to deeper nodes. The most confident branches challenge previous systematic arrangements of these snails, grouping species that previously had been placed into different subgenera. High levels of sequence divergence within some species may reflect the long-term isolation of subpopulations. Rapid rates of sequence divergence may have saturated base substitutions and contributed to the lack of resolution of higher-order relationships. We did not find support for the monophyly of the Achatinella species, nor thus for a single colonization of Oahu from Maui Nui.     

Molecular phylogenies and evolution of crt genes in algae

The carotenoids constitute the most widespread class of pigments in nature. Most previous work has concentrated on the identification and characterization of their chemical physical properties and bioavailability. In recent years, significant amounts of research have been conducted in an attempt to analyze the genes and the molecular regulation of the genes involved in the biosynthesis of carotenoids. However, it is important not to lose sight of the early evolution of carotenoid biosynthesis. One of the major obstacles in understanding the evolution of the respective enzymes and their patterns of selection is a lack of a well-supported phylogenic analysis. In the present research, a major long-term objective was to provide a clearer picture of the evolutionary history of genes, together with an evaluation of the patterns of selection in algae. These phylogenies will be important in studies characterizing the evolution of algae. The gene sequences of the enzymes involved in the major steps of the carotenoid biosynthetic pathway in algae (cyanobacteria, rhofophyta, chlorophyta) have been analyzed. Phylogenetic relationships among protein-coding DNA sequences were reconstructed by neighbor-joining (NJ) analysis for the respective carotenoid biosynthetic pathway genes (crt) in algae. The analysis also contains an estimation of the rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)), synonymous nucleotide substitution per synonymous site (d(S)), and the ratio of nonsynonmous (d(N)/d(S)) for the test of selection patterns. The phylogenetic trees show that the taxa of some genera have a closer evolutionary relationship with other genera in some gene sequences, which suggests a common ancient origin and that lateral gene transfer has occurred among unrelated genera. The d(N) values of crt genes in the early pathway are relatively low, while those of the following steps are slightly higher, while the d(N) values of crt genes in chlorophyta are higher than those in cyanobacteria. Most of the d(N)/d(S) values exceed 1. The phylogenetic analysis revealed that lateral gene transfer may have taken place across algal genomes and the d(N) values suggest that most of the early crt genes are well conserved compared to the later crt genes. Furthermore, d(N) values also revealed that the crt genes of chlorophyta are more evolutionary than cyanobacteria. The amino acids' changes are mostly adaptive evolution under the influence of positive diversity selection.     

MLGO: phylogeny reconstruction and ancestral inference from gene-order data

Background

The rapid accumulation of whole-genome data has renewed interest in the study of using gene-order data for phylogenetic analyses and ancestral reconstruction. Current software and web servers typically do not support duplication and loss events along with rearrangements.

Results

MLGO (Maximum Likelihood for Gene-Order Analysis) is a web tool for the reconstruction of phylogeny and/or ancestral genomes from gene-order data. MLGO is based on likelihood computation and shows advantages over existing methods in terms of accuracy, scalability and flexibility.

Conclusions

To the best of our knowledge, it is the first web tool for analysis of large-scale genomic changes including not only rearrangements but also gene insertions, deletions and duplications. The web tool is available from http://www.geneorder.org/server.php.