Identifying species of moths (Lepidoptera) from Baihua Mountain,Beijing, China,using DNA barcodes

DNA barcoding has become a promising means for the identification of organisms of all life‐history stages. Currently, distance‐based and tree‐based methods are most widely used to define species boundaries and uncover cryptic species. However, there is no universal threshold of genetic distance values that can be used to distinguish taxonomic groups. Alternatively, DNA barcoding can deploy a “character‐based” method, whereby species are identified through the discrete nucleotide substitutions. Our research focuses on the delimitation of moth species using DNA‐barcoding methods. We analyzed 393 Lepidopteran specimens belonging to 80 morphologically recognized species with a standard cytochrome c oxidase subunit I (COI) sequencing approach, and deployed tree‐based, distance‐based, and diagnostic character‐based methods to identify the taxa. The tree‐based method divided the 393 specimens into 79 taxa (species), and the distance‐based method divided them into 84 taxa (species). Although the diagnostic character‐based method found only 39 so‐identifiable species in the 80 species, with a reduction in sample size the accuracy rate substantially improved. For example, in the Arctiidae subset, all 12 species had diagnostics characteristics. Compared with traditional morphological method, molecular taxonomy performed well. All three methods enable the rapid delimitation of species, although they have different characteristics and different strengths. The tree‐based and distance‐based methods can be used for accurate species identification and biodiversity studies in large data sets, while the character‐based method performs well in small data sets and can also be used as the foundation of species‐specific biochips.     

Deriving phylogenetic trees from the similarity analysis of metabolic pathways

MOTIVATION: Comparative analysis of metabolic pathways in different genomes can give insights into the understanding of evolutionary and organizational relationships among species. This type of analysis allows one to measure the evolution of complete processes (with different functional roles) rather than the individual elements of a conventional analysis. We present a new technique for the phylogenetic analysis of metabolic pathways based on the topology of the underlying graphs. A distance measure between graphs is defined using the similarity between nodes of the graphs and the structural relationship between them. This distance measure is applied to the enzyme-enzyme relational graphs derived from metabolic pathways. Using this approach, pathways and group of pathways of different organisms are compared to each other and the resulting distance matrix is used to obtain a phylogenetic tree. RESULTS: We apply the method to the Citric Acid Cycle and the Glycolysis pathways of different groups of organisms, as well as to the Carbohydrate metabolic networks. Phylogenetic trees obtained from the experiments were close to existing phylogenies and revealed interesting relationships among organisms.     

Molecular distance and divergence time in carnivores and primates

Numerous studies have used indices of genetic distance between species to reconstruct evolutionary relationships and to estimate divergence time. However, the empirical relationship between molecular-based indices of genetic divergence and divergence time based on the fossil record is poorly known. To date, the results of empirical studies conflict and are difficult to compare because they differ widely in their choice of taxa, genetic techniques, or methods for calibrating rates of molecular evolution. We use a single methodology to analyze the relationship of molecular distance and divergence time in 86 taxa (72 carnivores and 14 primates). These taxa have divergence times of 0.01-55 Myr and provide a graded series of phylogenetic divergences such that the shape of the curve relating genetic distance and divergence time is often well defined. The techniques used to obtain genetic distance estimates include one- and two-dimensional protein electrophoresis, DNA hybridization, and microcomplement fixation. Our results suggest that estimates of molecular distance and divergence time are highly correlated. However, rates of molecular evolution are not constant; rather, in general they decline with increasing divergence time in a linear fashion. The rate of decline may differ according to technique and taxa. Moreover, in some cases the variability in evolutionary rates changes with increasing divergence time such that the accuracy of nodes in a phylogenetic tree varies predictably with time.     

Endless forms most stupid,icky, and small: The preponderance of noncharismatic invertebrates as integral to a biologically sound view of life

Big, beautiful organisms are useful for biological education, increasing evolution literacy, and biodiversity conservation. But if educators gloss over the ubiquity of streamlined and miniaturized organisms, they unwittingly leave students and the public vulnerable to the idea that the primary evolutionary plot of every metazoan lineage is “progressive” and "favors" complexity. We show that simple, small, and intriguingly repulsive invertebrate animals provide a counterpoint to misconceptions about evolution. Our examples can be immediately deployed in biology courses and outreach. This context emphasizes that chordates are not the pinnacle of evolution. Rather, in the evolution of animals, miniaturization, trait loss, and lack of perfection are at least as frequent as their opposites. Teaching about invertebrate animals in a “tree thinking” context uproots evolution misconceptions (for students and the public alike), provides a mental scaffold for understanding all animals, and helps to cultivate future ambassadors and experts on these little‐known, weird, and fascinating taxa.     

Construction of phylogenetic profiles based on the genetic distance of hundreds of genomes

Phylogenetic profiles have been widely applied in functional genomics research, especially in the prediction of protein-protein interactions (PPIs). A key issue in phylogenetic profiling is how to effectively select reference organisms from the available hundreds of genomes. In this study, we performed an assessment of reference organism selection based on the genetic distance between the target organism and 167 reference organisms. We found that inclusion of reference organisms from all distance levels had better performance in the prediction of PPIs than that at each distance level. The PPI prediction reached an optimal level when 70% of the reference organisms at all distance levels were selected; and this performance was similar to that in the optimal condition based on the taxonomy tree in our previous study. Because measurement of genetic distance is direct and simple compared to the topology of the taxonomy tree, we suggest selecting reference organisms based on genetic distance in the construction of phylogenetic profiles.     

Taxon ordering in phylogenetic trees: a workbench test

Background  

Phylogenetic trees are an important tool for representing evolutionary relationships among organisms. In a phylogram or chronogram, the ordering of taxa is not considered meaningful, since complete topological information is given by the branching order and length of the branches, which are represented in the root-to-node direction. We apply a novel method based on a (λ + μ)-Evolutionary Algorithm to give meaning to the order of taxa in a phylogeny. This method applies random swaps between two taxa connected to the same node, without changing the topology of the tree. The evaluation of a new tree is based on different distance matrices, representing non-phylogenetic information such as other types of genetic distance, geographic distance, or combinations of these. To test our method we use published trees of Vesicular stomatitis virus, West Nile virus and Rice yellow mottle virus.     

A new electrophoretic approach to biochemical systematics of bees

SDS electrophoresis of soluble proteins from 49 species of bees from 8 families produce a conservative character set. Differences between species appear to be actual differences in polypeptides and their concentrations rather than spurious differences. The data set based on the absorption values of bands proves to be more useful than does a data set based on the molecular weights of bands present or one based on the presence and absence of bands. Numerical analysis of a matrix of distance coefficients based on absorption values of bands shows much overlap between families of bees. Also, some taxa of questionable affinities when compared morphologically can be shown to be biochemically quite similar to other members of their traditional morphological families. Other taxa, which appear very similar morphologically, are highly differentiated biochemically. Examples are presented to compare and contrast the data from the SDS-separated character set to predictions based on the theories of regulatory and constant rate protein evolution. The results of the study largely agreed with predictions based on constant rate protein evolution. Exceptions may be understood in the framework of regulatory evolution. Consequences of discordance between morphological and biochemical phenetics are discussed as they affect systematic classifications.     

Identifying and reconstructing lateral transfers from distance matrices by combining the minimum contradiction method and neighbor-net

Identifying lateral gene transfers is an important problem in evolutionary biology. Under a simple model of evolution, the expected values of an evolutionary distance matrix describing a phylogenetic tree fulfill the so-called Kalmanson inequalities. The Minimum Contradiction method for identifying lateral gene transfers exploits the fact that lateral transfers may generate large deviations from the Kalmanson inequalities. Here a new approach is presented to deal with such cases that combines the Neighbor-Net algorithm for computing phylogenetic networks with the Minimum Contradiction method. A subset of taxa, prescribed using Neighbor-Net, is obtained by measuring how closely the Kalmanson inequalities are fulfilled by each taxon. A criterion is then used to identify the taxa, possibly involved in a lateral transfer between nonconsecutive taxa. We illustrate the utility of the new approach by applying it to a distance matrix for Archaea, Bacteria, and Eukaryota.     

Determining evolutionary distances from highly diverged nucleic acid sequences: Operator metrics

Summary Operator metrics are explicity designed to measure evolutionary distances from nucleic acid sequences when substitution rates differ greatly among the organisms being compared, or when substitutions have been extensive. Unlike lengths calculated by the distance matrix and parsimony methods, in which substitutions in one branch of a tree can alter the measured length of another branch, lengths determined by operator metrics are not affected by substitutions outside the branch.In the method, lengths (operator metrics) corresponding to each of the branches of an unrooted tree are calculated. The metric length of a branch reconstructs the number of (transversion) differences between sequences at a tip and a node (or between nodes) of a tree. The theory is general and is fundamentally independent of differences in substitution rates among the organisms being compared. Mathematically, the independence has been obtained becuase the metrics are eigen vectors of fundamental equations which describe the evolution of all unrooted trees.Even under conditions when both the distance matrix method or a simple parsimony length method are show to indicate lengths than are an order of magnitude too large or too small, the operator metrics are accurate. Examples, using data calculated with evolutionary rates and branchings designed to confuse the measurement of branch lengths and to camouflage the topology of the true tree, demonstrate the validity of operator metrics. The method is robust. Operator metric distances are easy to calculated, can be extended to any number of taxa, and provide a statistical estimate of their variances.The utility of the method is demonstrated by using it to analyze the origins and evolutionary of chloroplasts, mitochondria, and eubacteria.     

Revealing the factors that promote speciation

What biological attributes of organisms promote speciation, and ultimately, species diversity? This question has a long history of interest, with proposed diversity promoters including attributes such as sexual selection, ecological specialism and dispersability. However, such ideas are difficult to test because the time-scale of processes involved is too great for direct human observation and experimentation. An increasingly powerful solution is to investigate diversity patterns among extant groups to infer the nature of processes operating during the evolution of those groups. This approach relies on the use of robust, phylogenetically based null models to overcome some of the problems inherent in observational inference. We illustrate this area by (i) discussing recent advances in identifying correlates of diversity among higher taxa, and (ii) proposing new methods for analysing patterns in species-level phylogenies, drawing examples from a wide range of organisms.     

Topology improves phylogenetic motif functional site predictions

Prediction of protein functional sites from sequence-derived data remains an open bioinformatics problem. We have developed a phylogenetic motif (PM) functional site prediction approach that identifies functional sites from alignment fragments that parallel the evolutionary patterns of the family. In our approach, PMs are identified by comparing tree topologies of each alignment fragment to that of the complete phylogeny. Herein, we bypass the phylogenetic reconstruction step and identify PMs directly from distance matrix comparisons. In order to optimize the new algorithm, we consider three different distance matrices and 13 different matrix similarity scores. We assess the performance of the various approaches on a structurally nonredundant data set that includes three types of functional site definitions. Without exception, the predictive power of the original approach outperforms the distance matrix variants. While the distance matrix methods fail to improve upon the original approach, our results are important because they clearly demonstrate that the improved predictive power is based on the topological comparisons. Meaning that phylogenetic trees are a straightforward, yet powerful way to improve functional site prediction accuracy. While complementary studies have shown that topology improves predictions of protein-protein interactions, this report represents the first demonstration that trees improve functional site predictions as well.     

GroEL-induced topological dislocation of a substrate protein β-sheet core: a solution EPR spin–spin distance study

The Hsp60-type chaperonin GroEL assists in the folding of the enzyme human carbonic anhydrase II (HCA II) and protects it from aggregation. This study was aimed to monitor conformational rearrangement of the substrate protein during the initial GroEL capture (in the absence of ATP) of the thermally unfolded HCA II molten-globule. Single- and double-cysteine mutants were specifically spin-labeled at a topological breakpoint in the β-sheet rich core of HCA II, where the dominating antiparallel β-sheet is broken and β-strands 6 and 7 are parallel. Electron paramagnetic resonance (EPR) was used to monitor the GroEL-induced structural changes in this region of HCA II during thermal denaturation. Both qualitative analysis of the EPR spectra and refined inter-residue distance calculations based on magnetic dipolar interaction show that the spin-labeled positions F147C and K213C are in proximity in the native state of HCA II at 20 °C (as close as ～8 Å), and that this local structure is virtually intact in the thermally induced molten-globule state that binds to GroEL. In the absence of GroEL, the molten globule of HCA II irreversibly aggregates. In contrast, a substantial increase in spin–spin distance (up to >20 Å) was observed within minutes, upon interaction with GroEL (at 50 and 60 °C), which demonstrates a GroEL-induced conformational change in HCA II. The GroEL binding-induced disentanglement of the substrate protein core at the topological break-point is likely a key event for rearrangement of this potent aggregation initiation site, and hence, this conformational change averts HCA II misfolding.     

Functional Evolution of Proteins

The functional evolution of proteins advances through gene duplication followed by functional drift, whereas molecular evolution occurs through random mutational events. Over time, protein active-site structures or functional epitopes remain highly conserved, which enables relationships to be inferred between distant orthologs or paralogs. In this study, we present the first functional clustering and evolutionary analysis of the RCSB Protein Data Bank (RCSB PDB) based on similarities between active-site structures. All of the ligand-bound proteins within the RCSB PDB were scored using our Comparison of Protein Active-site Structures (CPASS) software and database ( http://cpass.unl.edu/ ). Principal component analysis was then used to identify 4431 representative structures to construct a phylogenetic tree based on the CPASS comparative scores ( http://itol.embl.de/shared/jcatazaro ). The resulting phylogenetic tree identified a sequential, step-wise evolution of protein active-sites and provides novel insights into the emergence of protein function or changes in substrate specificity based on subtle changes in geometry and amino acid composition.     

Diskontinuierlich aktive Gene und Evolution Eine Diskussion am Beispiel der Archaeognatha (Insecta)

Discontinuous activity of genes and evolution - a discussion based on examples of the Archaeognatha (Insecta) On the basis of examples, mainly from the group of Archaeognatha (Insecta) it seems probable that, in multicellular organisms, there is not only a discontinuous activity of genes and groups of genes during the individual development, but also during the phylogenesis. The possibility that silent genes can be reactivated is probably an essential reason for the existence of parallelisms and homologous tendencies. Both conceptions have a complementary importance for an establishment of taxa ana the reconstruction of phylogenesis and are defined according to evolutionary genetics. The examples are discussed and the necessity of similar comparative studies in other supraspecific taxa is emphasized.     

Resolving tylenchid evolutionary relationships through multiple gene analysis derived from EST data

Sequence-based phylogenetic analyses typically are based on a small number of character sets and report gene trees which may not reflect the true species tree. We employed an EST mining strategy to suppress such incongruencies, and recovered the most robust phylogeny for five species of plant-parasitic nematode (Meloidogyne arenaria, M. chitwoodi, M. hapla, M. incognita, and M. javanica), three closely related tylenchid taxa (Heterodera glycines, Globodera pallida, and G. rostochiensis) and a distant taxon, Caenorhabditis elegans. Our multiple-gene approach is based on sampling more than 80,000 publicly available tylenchid EST sequences to identify phylum-wide orthologues. Bayesian inference, minimum evolution, maximum likelihood and protein distance methods were employed for phylogenetic reconstruction and hypothesis tests were constructed to elucidate differential selective pressures across the phylogeny for each gene. Our results place M. incognita and M. javanica as sister taxa, with M. arenaria as the next closely related nematode. Significant differences in selective pressure were revealed for some genes under some hypotheses, though all but one gene are exclusively under purifying selection, indicating conservation across the orthologous groups. This EST-based multi-gene analysis is a first step towards accomplishing genome-wide coverage for tylenchid evolutionary analyses.     

Phylogenetic comparison of metabolic capacities of organisms at genome level

Horizontal gene transfer (HGT) has been shown to widely spread in organisms by comparative genomic studies. However, its effect on the phylogenetic relationship of organisms, especially at a system level of different cellular functions, is still not well understood. In this work, we have constructed phylogenetic trees based on the enzyme, reaction, and gene contents of metabolic networks reconstructed from annotated genome information of 82 sequenced organisms. Results from different phylogenetic distance definitions and based on three different functional subsystems (i.e., metabolism, cellular processes, information storage and processing) were compared. Results based on the three different functional subsystems give different pictures on the phylogenetic relationship of organisms, reflecting the different extents of HGT in the different functional systems. In general, horizontal transfer is prevailing in genes for metabolism, but less in genes for information processing. Nevertheless, the major results of metabolic network-based phylogenetic trees are in good agreement with the tree based on 16S rRNA and genome trees, confirming the three domain classification and the close relationship between eukaryotes and archaea at the level of metabolic networks. These results strongly support the hypothesis that although HGT is widely distributed, it is nevertheless constrained by certain pre-existing metabolic organization principle(s) during the evolution. Further research is needed to identify the organization principle and constraints of metabolic network on HGT which have large impacts on understanding the evolution of life and in purposefully manipulating cellular metabolism.     

Towards building the tree of life: a simulation study for all angiosperm genera

Comprehensive phylogenetic trees are essential tools to better understand evolutionary processes. For many groups of organisms or projects aiming to build the Tree of Life, comprehensive phylogenetic analysis implies sampling hundreds to thousands of taxa. For the tree of all life this task rises to a highly conservative 13 million. Here, we assessed the performances of methods to reconstruct large trees using Monte Carlo simulations with parameters inferred from four large angiosperm DNA matrices, containing between 141 and 567 taxa. For each data set, parameters of the HKY85+G model were estimated and used to simulate 20 new matrices for sequence lengths from 100 to 10,000 base pairs. Maximum parsimony and neighbor joining were used to analyze each simulated matrix. In our simulations, accuracy was measured by counting the number of nodes in the model tree that were correctly inferred. The accuracy of the two methods increased very quickly with the addition of characters before reaching a plateau around 1000 nucleotides for any sizes of trees simulated. An increase in the number of taxa from 141 to 567 did not significantly decrease the accuracy of the methods used, despite the increase in the complexity of tree space. Moreover, the distribution of branch lengths rather than the rate of evolution was found to be the most important factor for accurately inferring these large trees. Finally, a tree containing 13,000 taxa was created to represent a hypothetical tree of all angiosperm genera and the efficiency of phylogenetic reconstructions was tested with simulated matrices containing an increasing number of nucleotides up to a maximum of 30,000. Even with such a large tree, our simulations suggested that simple heuristic searches were able to infer up to 80% of the nodes correctly.     

生命之树及其应用

生命之树的概念由达尔文在1859年提出, 用以反映分类群的亲缘关系和进化历史。近30年来, 随着建树性状种类的多样化、数据量的快速增长以及建树方法的不断发展和完善, 生命之树的规模越来越大, 可信度也越来越高。分子生物学、生态学、基因组学、生物信息学及计算机科学等的快速发展, 使得生命之树成为开展学科间交叉研究的桥梁, 其用途日益广泛。本文综述了生命之树研究的历史和现状, 介绍了生命之树在以下几个方面的应用: (1)通过构建不同尺度的生命之树, 理解生物类群间的系统发育关系; (2)通过时间估算和地理分布区重建, 推测现存生物的起源和地理分布格局及其成因; (3)基于时间树, 结合生态、环境因子及关键创新性状, 探讨生物的多样化进程和成因; (4)揭示生物多样性的来源和格局, 预测生物多样性动态变化, 并提出相应的保护策略。最后, 本文评估了生命之树在目前海量数据情况下遇到的序列比对困难、基因树冲突、“流浪类群”干扰等建树难题, 并指出了构建“超大树”的发展趋势。     

Modularity and integration in the hominoid scapula

In this paper, several hypotheses of morphological integration within the hominoid (ape) scapula are tested. In particular, whether the scapula represents a set of developmental tissues sharing tight correlations between constituent parts (i.e., highly integrated) or is more modularly organized (i.e., covariation is greater within regions than between) is tested. Whether the patterns of integration in the scapula have changed over phylogenetic time or in response to selective forces is also examined. Results from two different analyses (matrix correlations and edge deviance) indicate traits comprising the blade and acromion, and to a weaker degree the glenoid, correlate highly with each other. The coracoid exhibits more independence from other parts of the scapula, perhaps reflecting its distinct evolutionary developmental history. Overall, similarity in species-specific patterns of correlation was high between all taxa. Correlation matrix similarity was significantly correlated with functional similarity and morphological distance, but not with phylogenetic distance. These results are congruent with other studies of integration that suggest correlation patterns remain stable over evolutionary time. There are changes associated with phylogeny, but the tight link between functional similarity and phylogenetic distance at this level of comparison presents possible challenges to interpretation. Overall similarities in the pattern of integration in all taxa might be better interpreted as relative strengthening or weakening of trait correlations rather than broadscale changes in the pattern of relationship between developmental regions. Larger sample sizes with greater taxonomic/functional breadth, and finer scale analyses of patterns of correlation are needed to test these hypotheses further.     

Long branch attraction effects and the status of "basal eukaryotes": phylogeny and structural analysis of the ribosomal RNA gene cluster of the free-living diplomonad Trepomonas agilis

The three taxa emerging at the base of the eukaryotic ribosomal RNA phylogenetic tree (Diplomonadida, Microspora, and Parabasalia) include a wide array of parasitic species. and some free-living organisms that appear to be derived from a parasitic ancestry. The basal position of these taxa, which lack mitochondria, has recently been questioned. I sequenced most of the ribosomal RNA gene cluster of the free-living diplomonad Trepomonas agilis and a secondary structure model was reconstructed for the SSU rRNA. I conducted a RASA matrix analysis to identify, independently from tree reconstruction, putative long branch attraction effects in the data matrix. The results show that each of the basal clades and the euglenozoan clade act, indeed, as long branches and may have been engaged in a process of accelerated rate of evolution. A nucleotide signature analysis was conducted in the conserved regions for positions defining the three great domains of life (Eubacteria, Archea, and Eukaryota). For the three basal taxa, this analysis showed the presence of a significant number of different non-eukaryotic nucleotides. A precise study of the nature and location of these nucleotides led to conclusions supporting the results of the RASA analysis. Altogether, these findings suggest that the basal placement of these taxa in the SSU ribosomal RNA phylogenetic tree is artifactual, and flawed by long branch attraction effects.