Phylogenetic signal in AFLP data sets

AFLP markers provide a potential source of phylogenetic information for molecular systematic studies. However, there are properties of restriction fragment data that limit phylogenetic interpretation of AFLPs. These are (a) possible nonindependence of fragments, (b) problems of homology assignment of fragments, (c) asymmetry in the probability of losing and gaining fragments, and (d) problems in distinguishing heterozygote from homozygote bands. In the present study, AFLP data sets of Lactuca s.l. were examined for the presence of phylogenetic signal. An indication of this signal was provided by carrying out tree length distribution skewness (g1) tests, permutation tail probability (PTP) tests, and relative apparent synapomorphy analysis (RASA). A measure of the support for internal branches in the optimal parsimony tree (MPT) was made using bootstrap, jackknife, and decay analysis. Finally, the extent of congruence in MPTs for AFLP and internal transcribed spacer (ITS)-1 data sets for the same taxa was made using the partition homogeneity test (PHT) and the Templeton test. These analytical studies suggested the presence of phylogenetic signal in the AFLP data sets, although some incongruence was found between AFLP and ITS MPTs. An extensive literature survey undertaken indicated that authors report a general congruence of AFLP and ITS tree topologies across a wide range of taxonomic groups, suggesting that the present results and conclusions have a general bearing. In these earlier studies and those for Lactuca s.l., AFLP markers have been found to be informative at somewhat lower taxonomic levels than ITS sequences. Tentative estimates are suggested for the levels of ITS sequence divergence over which AFLP profiles are likely to be phylogenetically informative.     

壳斗科植物叶绿体基因组结构及变异分析

利用生物信息学方法比较壳斗科6个属14个物种的叶绿体基因组间差异,以近缘物种榛为外类群构建系统进化树,揭示壳斗科叶绿体基因组的结构特征及变异规律。结果显示,14种壳斗科植物的叶绿体基因组均为双链环状分子结构,大小在160 kB左右,差异较小,最大仅差1 366 bp;基因顺序基本一致,而基因数量有所差异,infA、petG、rpl22、ycf1、ycf15等多个基因在部分物种中发生丢失;主要有32个蛋白编码基因长度发生变异,其原因是内含子的丢失、内含子或者编码区的长度改变,华南锥基因长度变异较大;4个IR边界相对保守,但锥栗、Castanea pumila、华南锥3个物种由于边界扩张导致rps19基因部分序列进入到IR区;以榛为外类群构建的系统发育树,各进化支支持率较高,分辨率较好。研究结果表明,叶绿体基因组可以用于分析关系较近与进化较快物种的系统发生问题,为系统发育和进化研究提供依据。     

The phylogeny of the BEP clade in grasses revisited: evidence from the whole-genome sequences of chloroplasts

Despite the considerable efforts to reconstruct the phylogeny of grasses, the relationships among the subfamilies Bambusoideae, Pooideae and Ehrhartoideae in the BEP clade remain unresolved. Here we completely sequenced three chloroplast genomes of representative species from Bambusoideae and Ehrhartoideae and obtained 19 additional chloroplast genome sequences of other grasses from GenBank. Using sequences of 76 chloroplast protein-coding genes from the 22 grass species, we fully resolved the phylogeny of the BEP clade. Our results strongly supported the (B,P)E hypothesis, i.e., Bambusoideae and Pooideae are more closely related than Ehrhartoideae. This result was not biased by systematic or sampling errors and was impervious to phylogenetic methods or model specification. The divergence time estimate suggests that the initial diversification of the BEP clade into three subfamilies happened within a short time period (≈ 4 MY). The presence of these short internal branches may explain the inability of previous studies to achieve a confident resolution of the BEP clade. The combination of the sequences of the entire chloroplast genomes provided sufficient phylogenetic information to resolve the BEP phylogeny fully. These results provide a valuable evolutionary framework for comparative and functional genomic studies using the grass family as a model system.     

皱边喉毛花及其近缘物种的叶绿体基因组结构与进化分析

为重建喉毛花属下系统发育关系,明晰属下皱边喉毛花及其近缘种之间的物种关系。本研究利用Illumina高通量测序平台对12 个叶绿体基因组进行双末端测序,获得大量高质量的Clean reads用于后续生物信息学分析。结果表明：（1）喉毛花属下物种的基因组差异较小,均在150 kb左右,基因总数为131 个,其中编码基因81 个。IR区核苷酸多态性比SC低,编码区比非编码区更保守。（2）进化分析结果显示,几乎所有的编码基因受到纯化选择的作用。（3）密码子偏好性分析表明有35 个密码子的RSCU值均大于1,说明使用这些密码子的频率较高,各项密码子偏好性衡量指标说明喉毛花属物种的密码子偏好性较弱。（4）系统发育分析表明CDS、密码子位置与基因间隔区数据集构建的系统发育树具有高度一致的拓扑结构,大部分分支的支持率高。这些结果表明皱边喉毛花及其近缘种的叶绿体基因组无明显差异,在系统发育树上无法按物种聚类,也为后续展开喉毛花属下群体遗传学研究提供科学依据。     

A phylogenetic model for understanding the effect of gene duplication on cancer progression

As biotechnology advances rapidly, a tremendous amount of cancer genetic data has become available, providing an unprecedented opportunity for understanding the genetic mechanisms of cancer. To understand the effects of duplications and deletions on cancer progression, two genomes (normal and tumor) were sequenced from each of five stomach cancer patients in different stages (I, II, III and IV). We developed a phylogenetic model for analyzing stomach cancer data. The model assumes that duplication and deletion occur in accordance with a continuous time Markov Chain along the branches of a phylogenetic tree attached with five extended branches leading to the tumor genomes. Moreover, coalescence times of the phylogenetic tree follow a coalescence process. The simulation study suggests that the maximum likelihood approach can accurately estimate parameters in the phylogenetic model. The phylogenetic model was applied to the stomach cancer data. We found that the expected number of changes (duplication and deletion) per gene for the tumor genomes is significantly higher than that for the normal genomes. The goodness-of-fit test suggests that the phylogenetic model with constant duplication and deletion rates can adequately fit the duplication data for the normal genomes. The analysis found nine duplicated genes that are significantly associated with stomach cancer.     

Examining the utility of categorical models and alleviating artifacts in phylogenetic reconstruction of the Squamata (Reptilia)

Reconstruction artifacts are a serious hindrance to the elucidation of phylogenetic relationships and a number of methods have been devised to alleviate them. Previous studies have demonstrated a striking disparity in the evolutionary rates of the mitochondrial (mt) genomes of squamate reptiles (lizards, worm lizards and snakes) and the reconstruction artifacts that may arise from this. Here, to examine basal squamate relationships, we have added the mt genome of the blind skink Dibamus novaeguineae to the mitogenomic dataset and applied different models for resolving the squamate tree. Categorical models were found to be less susceptible to artifacts than were the commonly used noncategorical phylogenetic models GTR and mtREV. The application of different treatments to the data showed that the removal of the fastest evolving sites in snakes improved phylogenetic signal in the dataset. Basal divergences remained, nevertheless, poorly resolved. The proportion of both fast-evolving and conserved sites in the squamate mt genomes relative to sites with intermediate rates of evolution suggests rapid early divergences among squamate taxa and at least partly explains the short internal relative to external branches in the squamate tree. Thus, mt and nuclear trees may never reach full agreement because of the short branches characterizing these divergences.     

Untangling long branches: identifying conflicting phylogenetic signals using spectral analysis, neighbor-net, and consensus networks

Long-branch attraction is a well-known source of systematic error that can mislead phylogenetic methods; it is frequently invoked post hoc, upon recovering a different tree from the one expected based on prior evidence. We demonstrate that methods that do not force the data onto a single tree, such as spectral analysis, Neighbor-Net, and consensus networks, can be used to detect conflicting signals within the data, including those caused by long-branch attraction. We illustrate this approach using a set of taxa from three unambiguously monophyletic families within the Pelecaniformes: the darters, the cormorants and shags, and the gannets and boobies. These three families are universally acknowledged as forming a monophyletic group, but the relationship between the families remains contentious. Using sequence data from three mitochondrial genes (12S, ATPase 6, and ATPase 8) we demonstrate that the relationship between these three families is difficult to resolve because they are separated by a short internal branch and there are conflicting signals due to long-branch attraction, which are confounded with nonhomogeneous sequence evolution across the different genes. Spectral analysis, Neighbor-Net, and consensus networks reveal conflicting signals regarding the placement of one of the darters, with support found for darter monophyly, but also support for a conflicting grouping with the outgroup, pelicans. Furthermore, parsimony and maximum-likelihood analyses produced different trees, with one of the two most parsimonious trees not supporting the monophyly of the darters. Monte Carlo simulations, however, were not sensitive enough to reveal long-branch attraction unless the branches are longer than those actually observed. These results indicate that spectral analysis, Neighbor-Net, and consensus networks offer a powerful approach to detecting and understanding the source of conflicting signals within phylogenetic data.     

Phylogenetic congruence and discordance among one morphological and three molecular data sets from Pontederiaceae

A morphological data set and three sources of data from the chloroplast genome (two genes and a restriction site survey) were used to reconstruct the phylogenetic history of the pickerelweed family Pontederiaceae. The chloroplast data converged towards a single tree, presumably the true chloroplast phylogeny of the family. Unrooted trees estimated from each of the three chloroplast data sets were identical or extremely similar in shape to each other and mostly robustly supported. There was no evidence of significant heterogeneity among the data sets, and the few topological differences seen among unrooted trees from each chloroplast data set are probably artifacts of sampling error on short branches. Despite well-documented differences in rates of evolution for different characters in individual data sets, equally weighted parsimony permits accurate reconstructions of chloroplast relationships in Pontederiaceae. A separate morphology-based data set yielded trees that were very different from the chloroplast trees. Although there was substantial support from the morphological evidence for several major clades supported by chloroplast trees, most of the conflicting phylogenetic structure on the morphology trees was not robust. Nonetheless, several statistical tests of incongruence indicate significant heterogeneity between molecules and morphology. The source of this apparent incongruence appears to be a low ratio of phylogenetic signal to noise in the morphological data.     

The root of the mammalian tree inferred from whole mitochondrial genomes

Morphological and molecular data are currently contradictory over the position of monotremes with respect to marsupial and placental mammals. As part of a re-evaluation of both forms of data we examine complete mitochondrial genomes in more detail. There is a particularly large discrepancy in the frequencies of thymine and cytosine (T-C) between mitochondrial genomes that appears to affect some deep divergences in the mammalian tree. We report that recoding nucleotides to RY-characters, and partitioning maximum-likelihood analyses among subsets of data reduces such biases, and improves the fit of models to the data, respectively. RY-coding also increases the signal on the internal branches relative to external, and thus increases the phylogenetic signal. In contrast to previous analyses of mitochondrial data, our analyses favor Theria (marsupials plus placentals) over Marsupionta (monotremes plus marsupials). However, a short therian stem lineage is inferred, which is at variance with the traditionally deep placement of monotremes on morphological data.     

Origin and phylogeny of chloroplasts revealed by a simple correlation analysis of complete genomes

The complete sequenced genomes of chloroplast have provided much information on the origin and evolution of this organelle. In this paper we attempt to use these sequences to test a novel approach for phylogenetic analysis of complete genomes based on correlation analysis of compositional vectors. All protein sequences from 21 complete chloroplast genomes are analyzed in comparison with selected archaea, eubacteria, and eukaryotes. The distance-based analysis shows that the chloroplast genomes are most closely related to cyanobacteria, consistent with the endosymbiotic origin of chloroplasts. The chloroplast genomes are separated to two major clades corresponding to chlorophytes (green plants) s.l. and rhodophytes (red algae) s.l. The interrelationships among the chloroplasts are largely in agreement with the current understanding on chloroplast evolution. For instance, the analysis places the chloroplasts of two chromophytes (Guillardia and Odontella) within the rhodophyte lineage, supporting secondary endosymbiosis as the source of these chloroplasts. The relationships among the green algae and land plants in our tree also agree with results from traditional phylogenetic analyses. Thus, this study establishes the value of our simple correlation analysis in elucidating the evolutionary relationships among genomes. It is hoped that this approach will provide insights on comparative genome analysis.     

沙冬青属植物叶绿体基因组对比和系统发育分析

该研究以沙冬青和矮沙冬青叶绿体基因组为研究对象,比较分析其基因组结构和系统发育关系。结果显示: (1)沙冬青和矮沙冬青的叶绿体基因组具有典型的四段式结构,全长为153 935 bp和154 140 bp,其中大单拷贝区(LSC)分别为83 891 bp和84 126 bp,小单拷贝区(SSC)分别为18 022 bp和18 014 bp,以及长度各自为26 011 bp和26 000 bp的成对反向重复区(IRs);沙冬青和矮沙冬青均注释130个基因,包括85个蛋白编码基因(PCGs),37个转运RNA(tRNA)以及8个核糖体RNA(rRNA)。(2)沙冬青和矮沙冬青的叶绿体基因组中分别检测出26和15个回文重复序列,39和50个串联重复序列,23和34个散在重复序列。同时都鉴定出96个SSRs位点,包括74和73个单核苷酸重复,5和6个二核苷酸重复,以及各有17个复合SSR位点;边界分析显示两者IR区相似,但仍有一定差异。(3)通过近邻结合法(NJ)对沙冬青和矮沙冬青在内的17种蝶形花亚科以及2种云实亚科植物的叶绿体基因组构建的系统发育树显示,沙冬青和矮沙冬青以较高的支持率聚为一个独立分枝。该研究结果为沙冬青属的种间鉴别、SSR分子标记开发、保育工作、种群动态以及进一步研究坡塔里族的演化过程奠定了理论基础。     

丝兰属6种植物叶绿体基因组特征分析

为了理清丝兰属（Yucca）叶绿体基因组特征和序列变异情况,进行丝兰属植物叶绿体比较基因组学分析,并构建基于叶绿体基因组的系统发育树。利用高通量测序技术获得无刺龙舌兰（Y. treculeana）叶绿体基因组序列,结合丝兰属现已发表的叶绿体基因组,使用生物信息学方法对6种丝兰属植物叶绿体全基因组进行基本结构、重复序列、边界收缩与扩张以及序列变异分析等在内的比较基因组学研究,并进行系统发育分析。结果表明：6种丝兰属植物叶绿体基因组大小、基因的类型及数目相近,种间基因组结构比较保守;从丝兰属植物叶绿体基因组中检测到多条重复序列,其中SSR位点多是由单核苷酸、双核苷酸和四核苷酸组成,且偏好使用A、T碱基;根据核酸多态性指数π≥0.008,在6种丝兰属植物叶绿体基因组中筛选出了psbK-psbl-trnS-GCU、rpl20-rps12和ccsA-ndhD 3个高变异区域;基于叶绿体全基因组和LSC+SSC区序列构建的系统发育关系基本一致,确定了6种丝兰属植物间的系统发育关系,其中无刺龙舌兰与克雷塔罗丝兰（Y. queretaroensis）的亲缘关系最近。本研究测序获得了无刺龙舌兰叶绿体基因组,揭示了6种丝兰属植物叶绿体基因组特征和序列变异情况,明确了各物种间的亲缘关系,研究结果可为后续丝兰属植物分子标记开发及系统发育研究提供参考。     

A phylogenetic analysis of Doronicum (Asteraceae, Senecioneae) based on morphological, nuclear ribosomal (ITS), and chloroplast (trnL-F) evidence.

A phylogenetic analysis of the Old World genus Doronicum (26 species, 4 subspecies) based on sequence data of the internal transcribed spacer (ITS) region of the nuclear ribosomal DNA, the chloroplast spacer trnL-F, and morphology is presented. Congruence among the three data sets was explored by the computing of several indices, all of which suggest homogeneity between only the two molecular matrices. We argue that the morphological data set contains poor phylogenetic signal and advocate simultaneous analysis of the three data sets (total evidence approach) so that morphological characters are tested for homology by congruence with molecular data. The resulting phylogenetic hypothesis allows several well-supported conclusions including the placement of a Corsican endemic (D. corsicum), sister to the remainder of the genus, and the inference that an early southern European or Mediterranean diversification took place in the genus. Shifts in morphological characters (e.g., homocarpy to heterocarpy) are confirmed to have evolved several times. Results from comparative studies of sequence data of the chloroplast gene ndhF support inclusion of Doronicum in tribe Senecioneae.     

Relationships within Cornales and circumscription of Cornaceae-matK and rbcL sequence data and effects of outgroups and long branches

Phylogenetic relationships in Cornales were assessed using sequences rbcL and matK. Various combinations of outgroups were assessed for their suitability and the effects of long branches and outgroups on tree topology were examined using RASA 2.4 prior to conducting phylogenetic analyses. RASA identified several potentially problematic taxa having long branches in individual data sets that may have obscured phylogenetic signal, but when data sets were combined RASA no longer detected long branch problems. t(RASA) provides a more conservative measurement for phylogenetic signal than the PTP and skewness tests. The separate matK and rbcL sequence data sets were measured as not containing phylogenetic signal by RASA, but PTP and skewness tests suggested the reverse [corrected]. Nonetheless, the matK and rbcL sequence data sets suggested relationships within Cornales largely congruent with those suggested by the combined matK-rbcL sequence data set that contains significant phylogenetic signal as measured by t(RASA), PTP, and skewness tests. Our analyses also showed that a taxon having a long branch on the tree may not be identified as a "long-branched" taxon by RASA. The long branches identified by RASA had little effect on the arrangement of other taxa in the tree, but the placements of the long-branched taxa themselves were often problematic. Removing the long-branched taxa from analyses generally increased bootstrap support, often substantially. Use of non-optimal outgroups (as identified by RASA) decreased phylogenetic resolution in parsimony analyses and suggested different relationships in maximum likelihood analyses, although usually weakly supported clades (less than 50% support) were impacted. Our results do not recommend using t(RASA) as a sole criterion to discard data or taxa in phylogenetic analyses, but t(RASA) and the taxon variance ratio obtained from RASA may be useful as a guide for improved phylogenetic analyses. Results of parsimony and ML analyses of the sequence data using optimal outgroups suggested by RASA revealed four major clades within Cornales: (1) Curtisia-Grubbia, (2) Cornus-Alangium, (3) Nyssa-Camptotheca-Davidia-Mastixia-Diplopanax, and (4) Hydrangeaceae-Loasaceae, with clades (2) and (3) forming a monophyletic group sister to clade (4) and clade (1) sister to the remainder of Cornales. However, there was not strong bootstrap support for relationships among the major clades. The placement of Hydrostachys could not be reliably determined, although most analyses place the genus within Hydrangeaceae; ML analyses, for example, placed the genus as the sister of Hydrangeeae. Our results supported a Cornales including the systematically problematic Hydrostachys, a Cornaceae consisting of Cornus and Alangium, a Nyssaceae consisting of Nyssa and Camptotheca, a monogeneric Davidiaceae, a Mastixiaceae consisting of Mastixia and Diplopanax, and an expanded Grubbiaceae consisting of Grubbia and Curtisia, and two larger families, Hydrangeaceae and Loasaceae.     

Detecting and overcoming systematic errors in genome-scale phylogenies

Genome-scale data sets result in an enhanced resolution of the phylogenetic inference by reducing stochastic errors. However, there is also an increase of systematic errors due to model violations, which can lead to erroneous phylogenies. Here, we explore the impact of systematic errors on the resolution of the eukaryotic phylogeny using a data set of 143 nuclear-encoded proteins from 37 species. The initial observation was that, despite the impressive amount of data, some branches had no significant statistical support. To demonstrate that this lack of resolution is due to a mutual annihilation of phylogenetic and nonphylogenetic signals, we created a series of data sets with slightly different taxon sampling. As expected, these data sets yielded strongly supported but mutually exclusive trees, thus confirming the presence of conflicting phylogenetic and nonphylogenetic signals in the original data set. To decide on the correct tree, we applied several methods expected to reduce the impact of some kinds of systematic error. Briefly, we show that (i) removing fast-evolving positions, (ii) recoding amino acids into functional categories, and (iii) using a site-heterogeneous mixture model (CAT) are three effective means of increasing the ratio of phylogenetic to nonphylogenetic signal. Finally, our results allow us to formulate guidelines for detecting and overcoming phylogenetic artefacts in genome-scale phylogenetic analyses.     

Molecular phylogeny of chloroplast DNA of Nicotiana species

Amplified fragment length polymorphism (AFLP) analysis of chloroplast DNA was used to study the relationships within the genus Nicotiana. Resulting phylogenetic tree, reconstructed using the UPGMA method, generally agreed with the existing taxonomic classification based on morphological and cytogenetic data, as well as sequence comparison of the internal transcribed spacer of the nuclear ribosomal DNA.     

Evolution of ITS1 rDNA in the Digenea (Platyhelminthes: Trematoda): 3′ End Sequence Conservation and Its Phylogenetic Utility

A comparison of ribosomal internal transcribed spacer 1 (ITS1) elements of digenetic trematodes (Platyhelminthes) including unidentified digeneans isolated from Cyathura carinata (Crustacea: Isopoda) revealed DNA sequence similarities at more than half of the spacer at its 3′ end. Primary sequence similarity was shown to be associated with secondary structure conservation, which suggested that similarity is due to identity by descent and not chance. Using an analysis of apomorphies, the sequence data were shown to produce a distinct phylogenetic signal. This was confirmed by the consistency of results of different tree reconstruction methods such as distance approaches, maximum parsimony, and maximum likelihood. Morphological evidence additionally supported the phylogenetic tree based on ITS1 data and the inferred phylogenetic position of the unidentified digeneans of C. carinata met the expectations from known trematode life-cycle patterns. Although ribosomal ITS1 elements are generally believed to be too variable for phylogenetic analysis above the species or genus level, the overall consistency of the results of this study strongly suggests that this is not the case in digenetic trematodes. Here, 3′ end ITS1 sequence data seem to provide a valuable tool for elucidating phylogenetic relationships of a broad range of phylogenetically distinct taxa. Received: 20 October 1997 / Accepted: 24 March 1998     

Phylogenomics reveal a robust fungal tree of life

Our understanding of the tree of life (TOL) is still fragmentary. Until recently, molecular phylogeneticists have built trees based on ribosomal RNA sequences and selected protein sequences, which, however, usually suffered from lack of support for the deeper branches and inconsistencies probably due to limited subsampling of the entire genome. Now, phylogenetic hypotheses can be based on the analysis of full genomes. We used available complete genome data as well as the eukaryote orthologous group (KOG) proteins to reconstruct with confidence basal branches of the fungal TOL. Phylogenetic analysis of a core of 531 KOGs shared among 21 fungal genomes, three animal genomes and one plant genome showed a single tree with high support resulting from four different methods of phylogenetic reconstruction. The single tree that we inferred from our dataset showed excellent nodal support for each branch, suggesting that it reflects the true phylogenetic relationships of the species involved.     

Is the Felsenstein zone a fly trap?

Although long-branch attraction, the incorrect grouping of long lineages in a phylogeny because of systematic error, has been identified as a potential source of error in phylogenetic analysis for almost two decades, no empirical examples of the phenomenon exist. Here, I outline several criteria for identifying long-branch attraction and apply these criteria to 18S ribosomal DNA (rDNA) sequence data for 13 insects. Parsimony and minimum evolution with p distances group the two longest branches together (those leading to Strepsiptera and Diptera). Simulation studies show that the long branches are long enough to attract. When a tree is assumed in which Strepsiptera and Diptera are separated and many data sets are simulated for that tree (using the parameter estimates for that tree for the original data), parsimony analysis of the simulated data consistently groups Strepsiptera and Diptera. Analyses of the 18S rDNA sequences using methods that are less sensitive to the problem of long-branch attraction estimate trees in which the long branches are separate.     

The Consistent Phylogenetic Signal in Genome Trees Revealed by Reducing the Impact of Noise

Phylogenetic trees based on gene repertoires are remarkably similar to the current consensus of life history. Yet it has been argued that shared gene content is unreliable for phylogenetic reconstruction because of convergence in gene content due to horizontal gene transfer and parallel gene loss. Here we test this argument, by filtering out as noise those orthologous groups that have an inconsistent phylogenetic distribution, using two independent methods. The resulting phylogenies do indeed contain small but significant improvements. More importantly, we find that the majority of orthologous groups contain some phylogenetic signal and that the resulting phylogeny is the only detectable signal present in the gene distribution across genomes. Horizontal gene transfer or parallel gene loss does not cause systematic biases in the gene content tree.