Swedish Museum of Natural History

Since its introduction in 2001, MrBayes has grown in popularity as a software package for Bayesian phylogenetic inference using Markov chain Monte Carlo (MCMC) methods. With this note, we announce the release of version 3.2, a major upgrade to the latest official release presented in 2003. The new version provides convergence diagnostics and allows multiple analyses to be run in parallel with convergence progress monitored on the fly. The introduction of new proposals and automatic optimization of tuning parameters has improved convergence for many problems. The new version also sports significantly faster likelihood calculations through streaming single-instruction-multiple-data extensions (SSE) and support of the BEAGLE library, allowing likelihood calculations to be delegated to graphics processing units (GPUs) on compatible hardware. Speedup factors range from around 2 with SSE code to more than 50 with BEAGLE for codon problems. Checkpointing across all models allows long runs to be completed even when an analysis is prematurely terminated. New models include relaxed clocks, dating, model averaging across time-reversible substitution models, and support for hard, negative, and partial (backbone) tree constraints. Inference of species trees from gene trees is supported by full incorporation of the Bayesian estimation of species trees (BEST) algorithms. Marginal model likelihoods for Bayes factor tests can be estimated accurately across the entire model space using the stepping stone method. The new version provides more output options than previously, including samples of ancestral states, site rates, site d(N)/d(S) rations, branch rates, and node dates. A wide range of statistics on tree parameters can also be output for visualization in FigTree and compatible software.

Abstract The main features of the phylogeny program TNT are discussed. Windows versions have a menu interface, while Macintosh and Linux versions are command‐driven. The program can analyze data sets with discrete (additive, non‐additive, step‐matrix) as well as continuous characters (evaluated with Farris optimization). Effective analysis of large data sets can be carried out in reasonable times, and a number of methods to help identifying wildcard taxa in the case of ambiguous data sets are implemented. A variety of methods for diagnosing trees and exploring character evolution is available in TNT, and publication‐quality tree‐diagrams can be saved as metafiles. Through the use of a number of native commands and a simple but powerful scripting language, TNT allows the user an enormous flexibility in phylogenetic analyses or simulations. © The Willi Hennig Society 2008.

Six DNA regions were evaluated as potential DNA barcodes for Fungi, the second largest kingdom of eukaryotic life, by a multinational, multilaboratory consortium. The region of the mitochondrial cytochrome c oxidase subunit 1 used as the animal barcode was excluded as a potential marker, because it is difficult to amplify in fungi, often includes large introns, and can be insufficiently variable. Three subunits from the nuclear ribosomal RNA cistron were compared together with regions of three representative protein-coding genes (largest subunit of RNA polymerase II, second largest subunit of RNA polymerase II, and minichromosome maintenance protein). Although the protein-coding gene regions often had a higher percent of correct identification compared with ribosomal markers, low PCR amplification and sequencing success eliminated them as candidates for a universal fungal barcode. Among the regions of the ribosomal cistron, the internal transcribed spacer (ITS) region has the highest probability of successful identification for the broadest range of fungi, with the most clearly defined barcode gap between inter- and intraspecific variation. The nuclear ribosomal large subunit, a popular phylogenetic marker in certain groups, had superior species resolution in some taxonomic groups, such as the early diverging lineages and the ascomycete yeasts, but was otherwise slightly inferior to the ITS. The nuclear ribosomal small subunit has poor species-level resolution in fungi. ITS will be formally proposed for adoption as the primary fungal barcode marker to the Consortium for the Barcode of Life, with the possibility that supplementary barcodes may be developed for particular narrowly circumscribed taxonomic groups.

Peer Reviewed

Abstract Global climate change impacts can already be tracked in many physical and biological systems; in particular, terrestrial ecosystems provide a consistent picture of observed changes. One of the preferred indicators is phenology, the science of natural recurring events, as their recorded dates provide a high‐temporal resolution of ongoing changes. Thus, numerous analyses have demonstrated an earlier onset of spring events for mid and higher latitudes and a lengthening of the growing season. However, published single‐site or single‐species studies are particularly open to suspicion of being biased towards predominantly reporting climate change‐induced impacts. No comprehensive study or meta‐analysis has so far examined the possible lack of evidence for changes or shifts at sites where no temperature change is observed. We used an enormous systematic phenological network data set of more than 125 000 observational series of 542 plant and 19 animal species in 21 European countries (1971–2000). Our results showed that 78% of all leafing, flowering and fruiting records advanced (30% significantly) and only 3% were significantly delayed, whereas the signal of leaf colouring/fall is ambiguous. We conclude that previously published results of phenological changes were not biased by reporting or publication predisposition: the average advance of spring/summer was 2.5 days decade −1 in Europe. Our analysis of 254 mean national time series undoubtedly demonstrates that species' phenology is responsive to temperature of the preceding months (mean advance of spring/summer by 2.5 days°C −1 , delay of leaf colouring and fall by 1.0 day°C −1 ). The pattern of observed change in spring efficiently matches measured national warming across 19 European countries (correlation coefficient r =−0.69, P <0.001).

Amino acid sequence data are available for ribulose biphosphate carboxylase, plastocyanin, cytochrome c, and ferredoxin for a number of angiosperm families. Cladistic analysis of the data, including evaluation of all equally or almost equally parsimonious cladograms, shows that much homoplasy (parallelisms and reversals) is present and that few or no well supported monophyletic groups of families can be demonstrated. In one analysis of nine angiosperm families and 40 variable amino acid positions from three proteins, the most parsimonious cladograms were 151 steps long and contained 63 parallelisms and reversals (consistency index = 0.583). In another analysis of six families and 53 variable amino acid positions from four proteins, the most parsimonious cladogram was 161 steps long and contained 50 parallelisms and reversals (consistency index = 0.689). Single changes in both data matrices could yield most parsimonious cladograms with quite different topologies and without common monophyletic groups. Presently, amino acid sequence data are not comprehensive enough for phylogenetic reconstruction among angiosperms. More informative positions are needed, either from sequencing longer parts of the proteins or from sequencing more proteins from the same taxa.

Carbonate communities: The activity of anaerobic methane oxidizing microbes facilitates precipitation of vast quantities of authigenic carbonate at methane seeps. Here we demonstrate the significant role of carbonate rocks in promoting diversity by providing unique habitat and food resources for macrofaunal assemblages at seeps on the Costa Rica margin (400-1850 m). The attendant fauna is surprisingly similar to that in rocky intertidal shores, with numerous grazing gastropods (limpets and snails) as dominant taxa. However, the community feeds upon seep-associated microbes. Macrofaunal density, composition, and diversity on carbonates vary as a function of seepage activity, biogenic habitat and location. The macrofaunal community of carbonates at non-seeping (inactive) sites is strongly related to the hydrography (depth, temperature, O2) of overlying water, whereas the fauna at sites of active seepage is not. Densities are highest on active rocks from tubeworm bushes and mussel beds, particularly at the Mound 12 location (1000 m). Species diversity is higher on rocks exposed to active seepage, with multiple species of gastropods and polychaetes dominant, while crustaceans, cnidarians, and ophiuroids were better represented on rocks at inactive sites. Macro-infauna (larger than 0.3 mm) from tube cores taken in nearby seep sediments at comparable depths exhibited densities similar to those on carbonate rocks, but had lower diversity and different taxonomic composition. Seep sediments had higher densities of ampharetid, dorvilleid, hesionid, cirratulid and lacydoniid polychaetes, whereas carbonates had more gastropods, as well as syllid, chrysopetalid and polynoid polychaetes. Stable isotope signatures and metrics: The stable isotope signatures of carbonates were heterogeneous, as were the food sources and nutrition used by the animals. Carbonate δ13Cinorg values (mean = -26.98‰) ranged from -53.3‰ to +10.0‰, and were significantly heavier than carbonate δ13Corg (mean = -33.83‰), which ranged from -74.4‰ to -20.6‰. Invertebrates on carbonates had average δ13C (per rock) = -31.0‰ (range -18.5‰ to -46.5‰) and δ15N = 5.7‰ (range -4.5‰ to +13.4‰). Average δ13C values did not differ between active and inactive sites; carbonate fauna from both settings depend on chemosynthesis-based nutrition. Community metrics reflecting trophic diversity (SEAc, total Hull Area, ranges of δ13C and δ15N) and species packing (mean distance to centroid, nearest neighbor distance) also did not vary as a function of seepage activity or site. However, distinct isotopic signatures were observed among related, co-occurring species of gastropods and polychaetes, reflecting intense microbial resource partitioning. Overall, the substrate and nutritional heterogeneity introduced by authigenic seep carbonates act to promote diverse, uniquely adapted assemblages, even after seepage ceases. The macrofauna in these ecosystems remain largely overlooked in most surveys, but are major contributors to biodiversity of chemosynthetic ecosystems and the deep sea in general.

Abstract We present a revised classification for extant ferns, with emphasis on ordinal and familial ranks, and a synopsis of included genera. Our classification reflects recently published phylogenetic hypotheses based on both morphological and molecular data. Within our new classification, we recognize four monophyletic classes, 11 monophyletic orders, and 37 families, 32 of which are strongly supported as monophyletic. One new family, Cibotiaceae Korall, is described. The phylogenetic affinities of a few genera in the order Polypodiales are unclear and their familial placements are therefore tentative. Alphabetical lists of accepted genera (including common synonyms), families, orders, and taxa of higher rank are provided.

Abstract- Because they are designed to produced just one tree, neighbor-joining programs can obscure ambiguities in data. Ambiguities can be uncovered by resampling, but existing neighbor-joining programs may give misleading bootstrap frequencies because they do not suppress zero-length branches and/or are sensitive to the order of terminals in the data. A new procedure, parsimony jackknifing, overcomes these problems while running hundreds of times faster than existing programs for neighbor-joining bootstrapping. For analysis of large matrices, parsimony jackknifing is hundreds of thousands of times faster than extensive branch-swapping, yet is better able to screen out poorly-supported groups.

This paper reports the results from a second characterisation of the 91500 zircon, including data from electron probe microanalysis, laser ablation inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS), secondary ion mass spectrometry (SIMS) and laser fluorination analyses. The focus of this initiative was to establish the suitability of this large single zircon crystal for calibrating in situ analyses of the rare earth elements and oxygen isotopes, as well as to provide working values for key geochemical systems. In addition to extensive testing of the chemical and structural homogeneity of this sample, the occurrence of banding in 91500 in both backscattered electron and cathodoluminescence images is described in detail. Blind intercomparison data reported by both LA‐ICP‐MS and SIMS laboratories indicate that only small systematic differences exist between the data sets provided by these two techniques. Furthermore, the use of NIST SRM 610 glass as the calibrant for SIMS analyses was found to introduce little or no systematic error into the results for zircon. Based on both laser fluorination and SIMS data, zircon 91500 seems to be very well suited for calibrating in situ oxygen isotopic analyses.

Granitic plutonism is the principal agent of crustal differentiation, but linking granite emplacement to crust formation requires knowledge of the magmatic evolution, which is notoriously difficult to reconstruct from bulk rock compositions. We unlocked the plutonic archive through hafnium (Hf) and oxygen (O) isotope analysis of zoned zircon crystals from the classic hornblende-bearing (I-type) granites of eastern Australia. This granite type forms by the reworking of sedimentary materials by mantle-like magmas instead of by remelting ancient metamorphosed igneous rocks as widely believed. I-type magmatism thus drives the coupled growth and differentiation of continental crust.

TEST 02 - Elsevier's Scopus, the largest abstract and citation database of peer-reviewed literature. Search and access research from the science, technology, medicine, social sciences and arts and humanities fields.

Phylogenies are usually dated by calibrating interior nodes against the fossil record. This relies on indirect methods that, in the worst case, misrepresent the fossil information. Here, we contrast such node dating with an approach that includes fossils along with the extant taxa in a Bayesian total-evidence analysis. As a test case, we focus on the early radiation of the Hymenoptera, mostly documented by poorly preserved impression fossils that are difficult to place phylogenetically. Specifically, we compare node dating using nine calibration points derived from the fossil record with total-evidence dating based on 343 morphological characters scored for 45 fossil (4--20 complete) and 68 extant taxa. In both cases we use molecular data from seven markers (∼5 kb) for the extant taxa. Because it is difficult to model speciation, extinction, sampling, and fossil preservation realistically, we develop a simple uniform prior for clock trees with fossils, and we use relaxed clock models to accommodate rate variation across the tree. Despite considerable uncertainty in the placement of most fossils, we find that they contribute significantly to the estimation of divergence times in the total-evidence analysis. In particular, the posterior distributions on divergence times are less sensitive to prior assumptions and tend to be more precise than in node dating. The total-evidence analysis also shows that four of the seven Hymenoptera calibration points used in node dating are likely to be based on erroneous or doubtful assumptions about the fossil placement. With respect to the early radiation of Hymenoptera, our results suggest that the crown group dates back to the Carboniferous, ∼309 Ma (95% interval: 291--347 Ma), and diversified into major extant lineages much earlier than previously thought, well before the Triassic. [Bayesian inference; fossil dating; morphological evolution; relaxed clock; statistical phylogenetics.].

Programs for Bayesian inference of phylogeny currently implement a unique and fixed suite of models. Consequently, users of these software packages are simultaneously forced to use a number of programs for a given study, while also lacking the freedom to explore models that have not been implemented by the developers of those programs. We developed a new open-source software package, RevBayes, to address these problems. RevBayes is entirely based on probabilistic graphical models, a powerful generic framework for specifying and analyzing statistical models. Phylogenetic-graphical models can be specified interactively in RevBayes, piece by piece, using a new succinct and intuitive language called Rev. Rev is similar to the R language and the BUGS model-specification language, and should be easy to learn for most users. The strength of RevBayes is the simplicity with which one can design, specify, and implement new and complex models. Fortunately, this tremendous flexibility does not come at the cost of slower computation; as we demonstrate, RevBayes outperforms competing software for several standard analyses. Compared with other programs, RevBayes has fewer black-box elements. Users need to explicitly specify each part of the model and analysis. Although this explicitness may initially be unfamiliar, we are convinced that this transparency will improve understanding of phylogenetic models in our field. Moreover, it will motivate the search for improvements to existing methods by brazenly exposing the model choices that we make to critical scrutiny. RevBayes is freely available at http://www.RevBayes.com [Bayesian inference; Graphical models; MCMC; statistical phylogenetics.].

Phylogenetic inference is fundamental to our understanding of most aspects of the origin and evolution of life, and in recent years, there has been a concentration of interest in statistical approaches such as Bayesian inference and maximum likelihood estimation. Yet, for large data sets and realistic or interesting models of evolution, these approaches remain computationally demanding. High-throughput sequencing can yield data for thousands of taxa, but scaling to such problems using serial computing often necessitates the use of nonstatistical or approximate approaches. The recent emergence of graphics processing units (GPUs) provides an opportunity to leverage their excellent floating-point computational performance to accelerate statistical phylogenetic inference. A specialized library for phylogenetic calculation would allow existing software packages to make more effective use of available computer hardware, including GPUs. Adoption of a common library would also make it easier for other emerging computing architectures, such as field programmable gate arrays, to be used in the future. We present BEAGLE, an application programming interface (API) and library for high-performance statistical phylogenetic inference. The API provides a uniform interface for performing phylogenetic likelihood calculations on a variety of compute hardware platforms. The library includes a set of efficient implementations and can currently exploit hardware including GPUs using NVIDIA CUDA, central processing units (CPUs) with Streaming SIMD Extensions and related processor supplementary instruction sets, and multicore CPUs via OpenMP. To demonstrate the advantages of a common API, we have incorporated the library into several popular phylogenetic software packages. The BEAGLE library is free open source software licensed under the Lesser GPL and available from http://beagle-lib.googlecode.com. An example client program is available as public domain software.

A clear picture of animal relationships is a prerequisite to understand how the morphological and ecological diversity of animals evolved over time. Among others, the placement of the acoelomorph flatworms, Acoela and Nemertodermatida, has fundamental implications for the origin and evolution of various animal organ systems. Their position, however, has been inconsistent in phylogenetic studies using one or several genes. Furthermore, Acoela has been among the least stable taxa in recent animal phylogenomic analyses, which simultaneously examine many genes from many species, while Nemertodermatida has not been sampled in any phylogenomic study. New sequence data are presented here from organisms targeted for their instability or lack of representation in prior analyses, and are analysed in combination with other publicly available data. We also designed new automated explicit methods for identifying and selecting common genes across different species, and developed highly optimized supercomputing tools to reconstruct relationships from gene sequences. The results of the work corroborate several recently established findings about animal relationships and provide new support for the placement of other groups. These new data and methods strongly uphold previous suggestions that Acoelomorpha is sister clade to all other bilaterian animals, find diminishing evidence for the placement of the enigmatic Xenoturbella within Deuterostomia, and place Cycliophora with Entoprocta and Ectoprocta. The work highlights the implications that these arrangements have for metazoan evolution and permits a clearer picture of ancestral morphologies and life histories in the deep past.

Patterns of diversification and timing of evolution within Neoaves, which includes almost 95% of all bird species, are virtually unknown. On the other hand, molecular data consistently indicate a Cretaceous origin of many neoavian lineages and the fossil record seems to support an Early Tertiary diversification. Here, we present the first well-resolved molecular phylogeny for Neoaves, together with divergence time estimates calibrated with a large number of stratigraphically and phylogenetically well-documented fossils. Our study defines several well-supported clades within Neoaves. The calibration results suggest that Neoaves, after an initial split from Galloanseres in Mid-Cretaceous, diversified around or soon after the K/T boundary. Our results thus do not contradict palaeontological data and show that there is no solid molecular evidence for an extensive pre-Tertiary radiation of Neoaves.

International audience

Although extensive data exist on avian influenza in wild birds in North America, limited information is available from elsewhere, including Europe. Here, molecular diagnostic tools were employed for high-throughput surveillance of migratory birds, as an alternative to classical labor-intensive methods of virus isolation in eggs. This study included 36,809 samples from 323 bird species belonging to 18 orders, of which only 25 species of three orders were positive for influenza A virus. Information on species, locations, and timing is provided for all samples tested. Seven previously unknown host species for avian influenza virus were identified: barnacle goose, bean goose, brent goose, pink-footed goose, bewick's swan, common gull, and guillemot. Dabbling ducks were more frequently infected than other ducks and Anseriformes; this distinction was probably related to bird behavior rather than population sizes. Waders did not appear to play a role in the epidemiology of avian influenza in Europe, in contrast to the Americas. The high virus prevalence in ducks in Europe in spring as compared with North America could explain the differences in virus-host ecology between these continents. Most influenza A virus subtypes were detected in ducks, but H13 and H16 subtypes were detected primarily in gulls. Viruses of subtype H6 were more promiscuous in host range than other subtypes. Temporal and spatial variation in influenza virus prevalence in wild birds was observed, with influenza A virus prevalence varying by sampling location; this is probably related to migration patterns from northeast to southwest and a higher prevalence farther north along the flyways. We discuss the ecology and epidemiology of avian influenza A virus in wild birds in relation to host ecology and compare our results with published studies. These data are useful for designing new surveillance programs and are particularly relevant due to increased interest in avian influenza in wild birds.

In a scanning electron microscope (SEM) an electron beam sets up an omni-directional source of scattered electrons within a specimen. Diffraction of these electrons will occur simultaneously on all lattice planes in the sample and the backscattered electrons (BSE), which escape from the specimen, will form a diffraction pattern that can be imaged on a phosphor screen. This is the basis of electron backscatter diffraction (EBSD). Similar diffraction effects cause individual grains of different orientations to give different total BSE. SEM images that exploit this effect will show orientation contrast (OC). EBSD and OC imaging are SEM-based crystallographic tools.