Molecular systematics of the Brassicaceae: evidence from coding plastidic matK and nuclear Chs sequences

Phylogenetic relationships were inferred using nucleotide sequence variation of the nuclear-encoded chalcone synthase gene (Chs) and the chloroplast gene matK for members of five tribes from the family Brassicaceae to analyze tribal and subtribal structures. Phylogenetic trees from individual data sets are mostly in congruence with the results from a combined matK-Chs analysis with a total of 2721 base pairs, but with greater resolution and higher statistical support for deeper branching patterns. The analysis indicates that tribes Lepidieae, Arabideae, and Sisymbrieae are not monophyletic. Among taxa under study four different lineages each were detected in tribes Arabideae and Lepidieae, interspersed with taxa from tribes Sisymbrieae, Hesperideae, and Brassiceae. It is concluded that tribe Brassiceae might be the only monophyletic group of the traditional tribes. From our data we estimated several divergence times for different lineages among cruciferous plants: 5.8 mya (million years ago) for the Arabidopsis-Cardaminopsis split, 20 mya for the Brassica-Arabidopsis split, and ～40 mya for the age of the deepest split between the most basal crucifer Aethionema and remaining cruciferous taxa.     

Estimating errors and confidence intervals for branch lengths in phylogenetic trees by a bootstrap approach

A method, based on the bootstrap procedure, is proposed for the estimation of branch-length errors and confidence intervals in a phylogenetic tree for which equal rates of substitution among lineages do not necessarily hold. The method can be used to test whether an estimated internodal distance is significantly greater than zero. In the application of the method, any estimator of genetic distances, as well as any tree reconstruction procedure (based on distance matrices), can be used. Also the method is not limited by the number of species involved in the phylogenetic tree. An example of the application of the method in the reconstruction of the phylogenetic tree for the four hominoid species—human, chimpanzee, gorilla, and orangutan—is shown. Correspondence to: J. Dopazo     

基于18S rDNA遗传距离与GC含量对游走类纤毛虫的分子系统学研究

为了揭示游走类纤毛虫的系统发生,对寄生于淡水鱼类的车轮虫科中的6种车轮虫进行了18S rDNA的测序并获得了9个序列。采用了最大似然法(ML)与贝叶斯法(BI)对GenBank中所有游走类纤毛虫的18S rDNA序列进行了系统树的构建,并首次将SPSS与18S rDNA遗传距离结合分析了游走类纤毛虫的系统发生。研究结果进一步证实了车轮虫属(Trichodina)的非单系发生与小车轮虫属 (Trichodinella) 的有效性。此外,研究结合18S rDNA 的GC含量与遗传距离分析提出了游走类纤毛虫科属及种间新的鉴定依据: 18S rDNA 的GC含量可用于游走类纤毛虫的科属区分,且与游走类纤毛虫的分化密切相关; 18S rDNA的遗传距离在游走类纤毛虫的不同阶元中具有一定的阈值范围,即通常种内遗传距离阈值范围为0.000-0.005,属种间阈值范围为0.005-0.150,当遗传距离大于0.150时,则达到了科间水平。     

Allopolyploid evolution in Geinae (Colurieae: Rosaceae) – building reticulate species trees from bifurcating gene trees

A previous phylogenetic study of paralogous nuclear low-copy granule-bound starch synthase (GBSSI) gene sequences from polyploid and diploid species in Geinae indicated that the clade has experienced two major allopolyploid events in its history. These were estimated to have occurred several million years ago. In this extended study we test if the reticulate phylogenetic hypothesis for Geinae can be maintained when additional sequences are added. The results are compatible with the hypothesis and strengthen it in minor aspects. We also attempt to identify extant members of one of the inferred ancestral lineages of the allopolyploids. On the basis of previous molecular phylogenies, one specific group has been proposed to be the descendants of this taxon. However, none of the additional paralogues belong to this ancestral lineage. A general method is proposed for converting a bifurcating gene tree, with multiple paralogous low-copy gene sequences from allopolyploid taxa, into a reticulate species tree.     

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.     

Using RAD seq for reconstructing phylogenies of highly diverged taxa: A test using the tribe Scandiceae (Apiaceae)

The angiosperm Apiaceae tribe Scandiceae includes four major clades—subtribes Daucinae, Ferulinae, Torilidinae, and Scandicinae—that originated ca. 20 Mya. Although all four subtribes are highly supported in molecular analyses, and morphological data indicate a sister relationship between Daucinae and Torilidinae, their branching order has not been resolved using standard Sanger multilocus data. Therefore, in this study, we test the utility of genomic RAD seq data in resolving deep phylogenetic relationships (up to 20 Mya) in Apiaceae subfamily Apioideae, with special emphasis on tribe Scandiceae using 12 representative species. We used two bioinformatic pipelines, pyRAD and RADIS (based on STACKS), to assemble RAD seq data and we tested the influence of various combinations of parameters on the robustness of the inferred tree topologies. Although different data processing approaches produced alignments with various amounts of missing data, they converged to two well‐supported topologies, irrespective of the phylogenetic method applied. Highly supported trees showed Scandicinae as sister to all other clades and indicated that Daucinae and Torilidinae are sister groups, thus confirming the relationship inferred from morphology. We conclude that the RAD seq method can be successfully used to resolve deep relationships formed 20 Mya within Apiaceae. We provide recommendations for parameter settings in RADIS and pyRAD for the analysis of taxa that have accumulated considerable genomic divergence.     

Mitochondrial DNA evolution in primates: Transition rate has been extremely low in the lemur

Summary Based on mitochondrial DNA (mt-DNA) sequence data from a wide range of primate species, branching order in the evolution of primates was inferred by the maximum likelihood method of Felsenstein without assuming rate constancy among lineages. Bootstrap probabilities for being the maximum likelihood tree topology among alternatives were estimated without performing a maximum likelihood estimation for each resampled data set. Variation in the evolutionary rate among lineages was examined for the maximum likelihood tree by a method developed by Kishino and Hasegawa. From these analyses it appears that the transition rate of mtDNA evolution in the lemur has been extremely low, only about 1/10 that in other primate lines, whereas the transversion rate does not differ significantly from that of other primates. Furthermore, the transition rate in catarrhines, except the gibbon, is higher than those in the tarsier and in platyrrhines, and the transition rate in the gibbon is lower than those in other catarrhines. Branching dates in primate evolution were estimated by a molecular clock analysis of mtDNA, taking into account the rate of variation among different lines, and the results were compared with those estimated from nuclear DNA. Under the most likely model, where the evolutionary rate of mtDNA has been unifrom within a great apes/human calde, human/chimpanzee clustering is preferred to the alternative branching orders among human, chimpanzee, and gorilla.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Southern hemisphere springtails: could any have survived glaciation of Antarctica?

Throughout the Southern Hemisphere many terrestrial taxa have circum-Antarctic distributions. This pattern is generally attributed to ongoing dispersal (by wind, water, or migrating birds) or relict Gondwanan distributions. Few of these terrestrial taxa have extant representatives in Antarctica, but such taxa would contribute to our understanding of the evolutionary origins of the continental Antarctic fauna. Either these taxa have survived the harsh climate cooling in Antarctica over the last 23 Myr (Gondwanan/vicariance origin) or they have dispersed there more recently (<2 MYA). In this context, we examined mtDNA (COI) sequence variation in Cryptopygus and related extant Antarctic and subantarctic terrestrial springtails (Collembola). Sequence divergence was estimated under a maximum likelihood model (general time reversible+I+Gamma) between individuals from subantarctic islands, Australia, New Zealand, Patagonia, Antarctic Peninsula, and continental Antarctica. Recent dispersal/colonization (<2 MYA) of Cryptopygus species was inferred between some subantarctic islands, and there was a close association between estimated times of divergences based on a molecular clock and proposed geological ages of islands. Most lineages generally grouped according to geographic proximity or by inferred dispersal/colonization pathways. In contrast, the deep divergences found for the four endemic Antarctic species indicate that they represent a continuous chain of descent dating from the break up of Gondwana to the present. We suggest that the diversification of these springtail species (21-11 MYA) in ice-free glacial refugia throughout the Trans-Antarctic Mountains was caused by the glaciation of the Antarctic continent during the middle to late Miocene.     

Statistical properties of bootstrap estimation of phylogenetic variability from nucleotide sequences: II. Four taxa without a molecular clock

Summary The statistical properties of sample estimation and bootstrap estimation of phylogenetic variability from a sample of nucleotide sequences were studied by considering model trees of three taxa with an outgroup. The cases of constant and varying rates of nucleotide substitution were compared. From sequences obtained by simulation, phylogenetic trees were constructed by using the maximum parsimony (MP) and neighbor joining (NJ) methods. The effectiveness and consistency of the MP method were studied in terms of proportions of informative sites. The results of simulation showed that bootstrap estimation of the confidence level for an inferred phylogeny can be used even under unequal rates of evolution if the rate differences are not large so that the MP method is not misleading. The condition under which the MP method becomes misleading (inconsistent) is more stringent for slowly evolving sequences than for rapidly evolving ones, and it also depends on the length of the internal branch. If the rate differences are large so that the MP method becomes consistently misleading, then bootstrap estimation will reinforce an erroneous conclusion on topology. Similar conclusions apply to the NJ method with uncorrected distances. The NJ method with corrected distances performs poorly when the sequence length is short but can avoid the inconsistency problem if the sequence length is long and if the distances can be estimated accurately.Offprint requests to: W.-H. Li     

Which functional processes control the short-term effect of grazing on net primary production in grasslands?

For unicellular organisms, a lack of effects of local species richness on ecosystem function has been proposed due to their locally high species richness and their ubiquitous distribution. High dispersal ability and high individual numbers may enable unicellular taxa to occur everywhere. Using our own and published data sets on uni- and multicellular organisms, we conducted thorough statistical analyses to test whether (1) unicellular taxa show higher relative local species richness compared to multicellular taxa, (2) unicellular taxa show lower slopes of the species:area relationships and species:individuals relationships, and (3) the species composition of unicellular taxa is less influenced by geographic distance compared to multicellular taxa. We found higher local species richness compared to the global species pool for unicellular organisms than for metazoan taxa. The difference was significant if global species richness was conservatively estimated but not if extrapolated, and therefore higher richness estimates were used. Both microalgae and protozoans showed lower slopes between species richness and sample size (area or individuals) compared to macrozoobenthos, also indicating higher local species richness for unicellular taxa. The similarity of species composition of both benthic diatoms and ciliates decreased with increasing geographic distance. This indicated restricted dispersal ability of protists and the absence of ubiquity. However, a steeper slope between similarity and distance was found for polychaetes and corals, suggesting a stronger effect of distance on the dispersal of metazoans compared to unicellular taxa. In conclusion, we found partly different species richness patterns among uni- and multicellular eukaryotes, but no strict ubiquity of unicellular taxa. Therefore, the effect of local unicellular species richness on ecosystem function has to be reanalyzed. Macroecological patterns suggested for multicellular organisms may differ in unicellular communities.     

Prokaryote phylogeny based on ribosomal proteins and aminoacyl tRNA synthetases by using the compositional distance approach

1 Introduction The systematics of prokaryotes has been a chal-lenge in microbiology as there are too few morpho-logical characteristics that can be used for classifica-tion[1]. A major breakthrough took place in the 1970s when Carl Woese[2] and coworkers aligned the small subunit ribosomal RNA (SSU rRNA) sequences to infer phylogenetic relations. The recognition of Ar-chaea as a third domain of life in addition to Bacteria and Eukarya and the support to the endosymbiotic origin of chlor…     

Phylogenetic detection of recombination with a Bayesian prior on the distance between trees

Genomic regions participating in recombination events may support distinct topologies, and phylogenetic analyses should incorporate this heterogeneity. Existing phylogenetic methods for recombination detection are challenged by the enormous number of possible topologies, even for a moderate number of taxa. If, however, the detection analysis is conducted independently between each putative recombinant sequence and a set of reference parentals, potential recombinations between the recombinants are neglected. In this context, a recombination hotspot can be inferred in phylogenetic analyses if we observe several consecutive breakpoints. We developed a distance measure between unrooted topologies that closely resembles the number of recombinations. By introducing a prior distribution on these recombination distances, a Bayesian hierarchical model was devised to detect phylogenetic inconsistencies occurring due to recombinations. This model relaxes the assumption of known parental sequences, still common in HIV analysis, allowing the entire dataset to be analyzed at once. On simulated datasets with up to 16 taxa, our method correctly detected recombination breakpoints and the number of recombination events for each breakpoint. The procedure is robust to rate and transitionratiotransversion heterogeneities for simulations with and without recombination. This recombination distance is related to recombination hotspots. Applying this procedure to a genomic HIV-1 dataset, we found evidence for hotspots and de novo recombination.     

The phylogenetic status of arthropods, as inferred from 18S rRNA sequences

Partial 18S rRNA sequences of five chelicerate arthropods plus a crustacean, myriapod, insect, chordate, echinoderm, annelid, and platyhelminth were compared. The sequence data were used to infer phylogeny by using a maximum-parsimony method, an evolutionary-distance method, and the evolutionary-parsimony method. The phylogenetic inferences generated by maximum-parsimony and distance methods support both monophyly of the Arthropoda and monophyly of the Chelicerata within the Arthropoda. These results are congruent with phylogenies based on rigorous cladistic analyses of morphological characters. Results support the inclusion of the Arthropoda within a spiralian or protostome coelomate clade that is the sister group of a deuterostome clade, refuting the hypothesis that the arthropods represent the "primitive" sister group of a protostome coelomate clade. Bootstrap analyses and consideration of all trees within 1% of the length of the most parsimonious tree suggest that relationships between the nonchelicerate arthropods and relationships within the chelicerate clade cannot be reliably inferred with the partial 18S rRNA sequence data. With the evolutionary-parsimony method, support for monophyly of the Arthropoda is found in the majority of the combinations analyzed if the coelomates are used as "outgroups." Monophyly of the Chelicerata is supported in most combinations assessed. Our analyses also indicate that the evolutionary-parsimony method, like distance and parsimony, may be biased by taxa with long branches. We suggest that a previous study's inference of the Arthropoda as paraphyletic may be the result of (a) having two few arthropod taxa available for analysis and (b) including long-branched taxa.      

Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals

There have been increasing efforts to relate drug efficacy and disease predisposition with genetic polymorphisms. We present statistical tests for association of haplotype frequencies with discrete and continuous traits in samples of unrelated individuals. Haplotype frequencies are estimated through the expectation-maximization algorithm, and each individual in the sample is expanded into all possible haplotype configurations with corresponding probabilities, conditional on their genotype. A regression-based approach is then used to relate inferred haplotype probabilities to the response. The relationship of this technique to commonly used approaches developed for case-control data is discussed. We confirm the proper size of the test under H(0) and find an increase in power under the alternative by comparing test results using inferred haplotypes with single-marker tests using simulated data. More importantly, analysis of real data comprised of a dense map of single nucleotide polymorphisms spaced along a 12-cM chromosomal region allows us to confirm the utility of the haplotype approach as well as the validity and usefulness of the proposed statistical technique. The method appears to be successful in relating data from multiple, correlated markers to response.     

Phylogeny of Trichomonads Inferred from Small-Subunit rRNA Sequences

ABSTRACT. Small subunit (16S-like) ribosomal RNA sequences were obtained from representatives of all four families constituting the order Trichomonadida. Comparative sequence analysis revealed that the Trichomonadida are a monophyletic lineage and a deep branch of the eukaryotic tree. Relative to other early divergent eukaryotic assemblages the branching pattern within the Trichomonadida is very shallow. This pattern suggests the Trichomonadida radiated recently, perhaps in conjunction with their animal hosts. From a morphological perspective the Devescovinidae and Calonymphidae are considered more derived than the Monocercomonadidae and Trichomonadidae. Molecular trees inferred by distance, parsimony and likelihood techniques consistently show the Devescovinidae and Calonymphidae are the earliest diverging lineages within the Trichomonadida, however bootstrap values do not strongly support a particular branching order. In an analysis of all known 16S-like ribosomal RNA sequences, the Trichomonadida share most recent common ancestry with unidentified protists from the hindgut of the termite Reticulitermes flavipes. The position of two putative free-living trichomonads in the tree is indicative of derivation from symbionts rather than direct descent from some free-living ancestral trichomonad.     

Evolutionary Relationships Among the Eukaryotic Crown Taxa Taking into Account Site-to-Site Rate Variation in 18S rRNA

In this study we constructed a bootstrapped distance tree of 500 small subunit ribosomal RNA sequences from organisms belonging to the so-called crown of eukaryote evolution. Taking into account the substitution rate of the individual nucleotides of the rRNA sequence alignment, our results suggest that (1) animals, true fungi, and choanoflagellates share a common origin: The branch joining these taxa is highly supported by bootstrap analysis (bootstrap support [BS] > 90%), (2) stramenopiles and alveolates are sister groups (BS = 75%), (3) within the alveolates, dinoflagellates and apicomplexans share a common ancestor BS > 95%), while in turn they both share a common origin with the ciliates (BS > 80%), and (4) within the stramenopiles, heterokont algae, hyphochytriomycetes, and oomycetes form a monophyletic grouping well supported by bootstrap analysis (BS > 85%), preceded by the well-supported successive divergence of labyrinthulomycetes and bicosoecids. On the other hand, many evolutionary relationships between crown taxa are still obscure on the basis of 18S rRNA. The branching order between the animal-fungal-choanoflagellates clade and the chlorobionts, the alveolates and stramenopiles, red algae, and several smaller groups of organisms remains largely unresolved. When among-site rate variation is not considered, the inferred tree topologies are inferior to those where the substitution rate spectrum for the 18S rRNA is taken into account. This is primarily indicated by the erroneous branching of fast-evolving sequences. Moreover, when different substitution rates among sites are not considered, the animals no longer appear as a monophyletic grouping in most distance trees. Received: 11 June 1997 / Accepted: 21 July 1997     

The pitfalls of molecular phylogeny based on four species,as illustrated by the Cetacea/Artiodactyla relationships

We study the reliability of phylogeny based on four taxa, when the internal, ancestral, branch is short. Such a quartet approach has been broadly used for inferring phylogenetic patterns. The question of branching pattern between the suborders Ruminantia and Suiformes (order Artiodactyla) and the order Cetacea is chosen as an example. All the combinations of four taxa were generated by taking on and only one species per group under study (three ingroups and one outgroup). Using real sequences, the analysis of these combinations demonstrates that the quartet approach is seriously misleading. Using both maximum parsimony and distance methods, it is possible to find a quartet of species which provided a high bootstrap proportion for each of the three possible unrooted trees. With the same set of sequences, we used all the available species simultaneously to construct a molecular phylogeny. This approach proved much more reliable than the quartet approach. When the number of informative sites is rather low, the branching patterns are not supported through bootstrap analysis, preventing us from false inference due to the lack of information. The reliable resolution of the phylogenetic relationships among Ruminantia, Suiformes, and Cetacea will therefore require a large number of nucleotides, such as the complete mitochondrial genomes of at least 30 species.     

Relative age inference in paleontology

Fossil taxa often occur in a regular vertical order in strata, yet, this regularity does not necessarily imply that the taxa succeeded one another in the same order in time. An argument for the time significance of such a regular pattern is strengthened if (1) the taxa involved were preserved in a wide range of depositional environments, and (2) the Eldredge-Gould model of evolution is correct: by contrast, (3) the areal extent of the pattern and (4) the inferred phylogeny of the taxa involved are not relevant. Paleontologists commonly do conclude that a homotaxial pattern, i.e. a definite orderly succession of fossil taxa, results from the taxa having occurred in essentially the same order in time. The conclusion may be justified as follows: If it is false, there should have been some local areas studied where the taxa (or first and last occurrences of taxa) occur in reverse order, or where two or more allegedly sequential taxa occur together in the same strata. But this is not the case: so the taxa probably succeeded one another in essentially the same order in time. The argument is broadly applicable, yet compelling only if (1) certain constraints involving the prior probability of the conclusion apply, and (2) transgression/regression can be ruled out as a cause of the pattern observed. Subsidiary defenses of the time significance of a homotaxial pattern involve other-group, radiometric, marker-bed, magnetic, seismic, or stable-isotope data. The use of fossils for time correlation of strata does not involve circular reasoning.     

Molecular evolution of a tandemly repeated trnF(GAA) gene in the chloroplast genomes of Microseris (Asteraceae) and the use of structural mutations in phylogenetic analyses

We sequenced the first ca. 900 bp of the 5'-trnL(UAA)-trnV(UAC)/ndhJ region of the chloroplast DNA of different Microseris accessions in order to resolve homoplasious length variation detected in the trnL(UAA)-trnF(GAA) region. We found two to four tandemly repeated trnF genes in the species of Microseris (Asteraceae, Lactuceae) and two in their sister genus Uropappus. Sequences indicated nonhomologous transitions between two, three, and four trnF genes in different Microseris taxa. Independent origins of similar trnF copy numbers were inferred from a chloroplast phylogeny of Microseris. The taxa involved grow on separate continents, supporting parallel origins of similar length variants. The changes in trnF copy numbers were best explained by interchromosomal recombination with unequal crossing over. The 5' copies of the repeats showed the highest sequence conservation, suggesting that these copies are likely to be functional trnF genes, whereas the other ones probably represent pseudogenes. Our results show that length polymorphisms accumulate once a duplicated sequence has become incorporated. Due to parallel gains of similar trnF copy numbers, homoplasious length variation was introduced into the data matrix. The data demonstrate that length polymorphisms cannot be used as indicators for phylogenetic distance unless they can be analyzed at the sequence level.     

Prokaryote phylogeny based on ribosomal proteins and aminoacyl tRNA synthetases by using the compositional distance approach

In order to show that the newly developed K-string composition distance method, based on counting oligopeptide frequencies, for inferring phylogenetic relations of prokaryotes works equally well without requiring the whole proteome data, we used all ribosomal proteins and the set of aminoacyl tRNA synthetases for each species. The latter group has been known to yield inconsistent trees if used individually. Our trees are obtained without making any sequence alignment. Altogether 16 Archaea, 105 Bacteria and 2 Eucarya are represented on the tree. Most of the lower branchings agree well with the latest, 2003, Outline of the second edition of the Bergey’s Manual of Systematic Bacteriology and the trees also suggest some relationships among higher taxa.