Following the rivers: historical reconstruction of California voles Microtus californicus (Rodentia: Cricetidae) in the deserts of eastern California

The California vole, Microtus californicus, restricted to habitat patches where water is available nearly year-round, is a remnant of the mesic history of the southern Great Basin and Mojave deserts of eastern California. The history of voles in this region is a model for species-edge population dynamics through periods of climatic change. We sampled voles from the eastern deserts of California and examined variation in the mitochondrial cytb gene, three nuclear intron regions, and across 12 nuclear microsatellite markers. Samples are allocated to two mitochondrial clades: one associated with southern California and the other with central and northern California. The limited mtDNA structure largely recovers the geographical distribution, replicated by both nuclear introns and microsatellites. The most remote population, Microtus californicus scirpensis at Tecopa near Death Valley, was the most distinct. This population shares microsatellite alleles with both mtDNA clades, and both its northern clade nuclear introns and southern clade mtDNA sequences support a hybrid origin for this endangered population. The overall patterns support two major invasions into the desert through an ancient system of riparian corridors along streams and lake margins during the latter part of the Pleistocene followed by local in situ divergence subsequent to late Pleistocene and Holocene drying events. Changes in current water resource use could easily remove California voles from parts of the desert landscape

RevBayes: Bayesian Phylogenetic Inference Using Graphical Models and an Interactive Model-Specification Language.

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.]

Probabilistic Graphical Model Representation in Phylogenetics

Recent years have seen a rapid expansion of the model space explored in statistical phylogenetics, emphasizing the need for new approaches to statistical model representation and software development. Clear communication and representation of the chosen model is crucial for: (1) reproducibility of an analysis, (2) model development and (3) software design. Moreover, a unified, clear and understandable framework for model representation lowers the barrier for beginners and non-specialists to grasp complex phylogenetic models, including their assumptions and parameter/variable dependencies. Graphical modeling is a unifying framework that has gained in popularity in the statistical literature in recent years. The core idea is to break complex models into conditionally independent distributions. The strength lies in the comprehensibility, flexibility, and adaptability of this formalism, and the large body of computational work based on it. Graphical models are well-suited to teach statistical models, to facilitate communication among phylogeneticists and in the development of generic software for simulation and statistical inference. Here, we provide an introduction to graphical models for phylogeneticists and extend the standard graphical model representation to the realm of phylogenetics. We introduce a new graphical model component, tree plates, to capture the changing structure of the subgraph corresponding to a phylogenetic tree. We describe a range of phylogenetic models using the graphical model framework and introduce modules to simplify the representation of standard components in large and complex models. Phylogenetic model graphs can be readily used in simulation, maximum likelihood inference, and Bayesian inference using, for example, Metropolis-Hastings or Gibbs sampling of the posterior distribution

Inferring the demographic history of the North American firefly Photinus pyralis

The firefly Photinus pyralis inhabits a wide range of latitudinal and ecological niches, with populations living from temperate to tropical habitats. Despite its broad distribution, its demographic history is unknown. In this study, we modelled and inferred different demographic scenarios for North American populations of P. pyralis, which were collected from Texas to New Jersey. We used a combination of ABC techniques (for multi-population/colonization analyses) and likelihood inference (dadi, StairwayPlot2, PoMo) for single-population demographic inference, which proved useful with our RAD data. We uncovered that the most ancestral North American population lays in Texas, which further colonized the Central region of the US and more recently the North Eastern coast. Our study confidently rejects a demographic scenario where the North Eastern populations colonized more southern populations until reaching Texas. To estimate the age of divergence between of P. pyralis, which provides deeper insights into the history of the entire species, we assembled a multi-locus phylogenetic data covering the genus Photinus. We uncovered that the phylogenetic node leading to P. pyralis lies at the end of the Miocene. Importantly, modelling the demographic history of North American P. pyralis serves as a null model of nucleotide diversity patterns in a widespread native insect species, which will serve in future studies for the detection of adaptation events in this firefly species, as well as a comparison for future studies of other North American insect taxa

Non-monophyly and intricate morphological evolution within the avian family Cettiidae revealed by multilocus analysis of a taxonomically densely sampled dataset

<p>Abstract</p> <p>Background</p> <p>The avian family Cettiidae, including the genera <it>Cettia</it>, <it>Urosphena</it>, <it>Tesia</it>, <it>Abroscopus </it>and <it>Tickellia </it>and <it>Orthotomus cucullatus</it>, has recently been proposed based on analysis of a small number of loci and species. The close relationship of most of these taxa was unexpected, and called for a comprehensive study based on multiple loci and dense taxon sampling. In the present study, we infer the relationships of all except one of the species in this family using one mitochondrial and three nuclear loci. We use traditional gene tree methods (Bayesian inference, maximum likelihood bootstrapping, parsimony bootstrapping), as well as a recently developed Bayesian species tree approach (*BEAST) that accounts for lineage sorting processes that might produce discordance between gene trees. We also analyse mitochondrial DNA for a larger sample, comprising multiple individuals and a large number of subspecies of polytypic species.</p> <p>Results</p> <p>There are many topological incongruences among the single-locus trees, although none of these is strongly supported. The multi-locus tree inferred using concatenated sequences and the species tree agree well with each other, and are overall well resolved and well supported by the data. The main discrepancy between these trees concerns the most basal split. Both methods infer the genus <it>Cettia </it>to be highly non-monophyletic, as it is scattered across the entire family tree. Deep intraspecific divergences are revealed, and one or two species and one subspecies are inferred to be non-monophyletic (differences between methods).</p> <p>Conclusions</p> <p>The molecular phylogeny presented here is strongly inconsistent with the traditional, morphology-based classification. The remarkably high degree of non-monophyly in the genus <it>Cettia </it>is likely to be one of the most extraordinary examples of misconceived relationships in an avian genus. The phylogeny suggests instances of parallel evolution, as well as highly unequal rates of morphological divergence in different lineages. This complex morphological evolution apparently misled earlier taxonomists. These results underscore the well-known but still often neglected problem of basing classifications on overall morphological similarity. Based on the molecular data, a revised taxonomy is proposed. Although the traditional and species tree methods inferred much the same tree in the present study, the assumption by species tree methods that all species are monophyletic is a limitation in these methods, as some currently recognized species might have more complex histories.</p

Polymorphism‐aware estimation of species trees and evolutionary forces from genomic sequences with RevBayes

Funding: Funding information Austrian Science Fund, Grant/Award Number: P34524-B; Biotechnology and Biological Sciences Research Council, Grant/Award Number: BB/W000768/1; Deutsche Forschungsgemeinschaft, Grant/Award Number: HO 6201/1-1; Vienna Science and Technology Fund, Grant/Award Number: MA016-061.1. The availability of population genomic data through new sequencing technologies gives unprecedented opportunities for estimating important evolutionary forces such as genetic drift, selection and mutation biases across organisms. Yet, analytical methods that can handle polymorphisms jointly with sequence divergence across species are rare and not easily accessible to empiricists. 2. We implemented polymorphism-aware phylogenetic models (PoMos), an alternative approach for species tree estimation, in the Bayesian phylogenetic software RevBayes. PoMos naturally account for incomplete lineage sorting, which is known to cause difficulties for phylogenetic inference in species radiations, and scale well with genome-wide data. Simultaneously, PoMos can estimate mutation and selection biases. 3. We have applied our methods to resolve the complex phylogenetic relationships of a young radiation of Chorthippus grasshoppers, based on coding sequences. In addition to establishing a well-supported species tree, we found a mutation bias favouring AT alleles and selection bias promoting the fixation of GC alleles, the latter consistent with GC-biased gene conversion. The selection bias is two orders of magnitude lower than genetic drift, validating the critical role of nearly neutral evolutionary processes in species radiation. 4. PoMos offer a wide range of models to reconstruct phylogenies and can be easily combined with existing models in RevBayes—for example, relaxed clock and divergence time estimation—offering new insights into the evolutionary processes underlying molecular evolution and, ultimately, species diversification.Publisher PDFPeer reviewe

Estimating the duration of speciation from phylogenies

Speciation is not instantaneous but takes time. The protracted birth-death diversification model incorporates this fact and predicts the often observed slowdown of lineage accumulation toward the present. The mathematical complexity of the protracted speciation model has barred estimation of its parameters until recently a method to compute the likelihood of phylogenetic branching times under this model was outlined (Lambert et al. ). Here, we implement this method and study using simulated phylogenies of extant species how well we can estimate the model parameters (rate of initiation of speciation, rate of extinction of incipient and good species, and rate of completion of speciation) as well as the duration of speciation, which is a combination of the aforementioned parameters. We illustrate our approach by applying it to a primate phylogeny. The simulations show that phylogenies often do not contain enough information to provide unbiased estimates of the speciation-initiation rate and the extinction rate, but the duration of speciation can be estimated without much bias. The estimate of the duration of speciation for the primate clade is consistent with literature estimates. We conclude that phylogenies combined with the protracted speciation model provide a promising way to estimate the duration of speciation.</p

PM10- und PM2.5- Emissionspotentiale von Substraten der Tagebaue im Lausitzer Revier

Im Lausitzer Revier werden aktuell 4 Braukohlen-Tagebaue betrieben, die als Quellen für Feinstaub-Emissionen gelten und somit zur lokalen Luftbelastungen beitragen. Berechnungen von möglichen Zusatzbelastungen durch den Tagebaubetrieb ergaben jedoch große Differenzen zu Messungen der Behörden vor Ort. Die Ursachen hierfür liegen in der starren Handhabung von Emissionsfaktoren, die vor allem die durch Winderosion hervorgerufene flächenhafte Emission von PM10 und PM2.5 stark überschätzen. Die Substrate der Hauptarbeitsebenen aller Tagebaue wurden untersucht, um die Materialeigenschaften als auch die Oberflächeneigenschaften, die die Emissionen beeinflussen, zu charakterisieren. Im ersten Schritt wurde mittels Horizontal-Querstromsichtung das Emissionspotential aller Substrate im luftgetrockneten Zustand ermittelt. Hierbei wird bei einer Windgeschwindigkeit von 3 ms-1 das Probenmaterial am Anfang des Windkanals von oben zugeführt und durch die Schwerkraft und die horizontale Strömung nach Größe und aerodynamischen Eigenschaften über die 7 m lange Messstrecke sortiert. Am Ende des Windkanals erfolgte die Messung der Partikelgrößenverteilung der Staubfraktion. Einzelne Proben wurden behutsam rückbefeuchtet und ebenfalls auf diese Weise untersucht. Für Untersuchungen zum Einfluss der Winderosion auf die PM-Emissionen wurde die Messstrecke in voller Länge mit den Substraten befüllt und mit Windgeschwindigkeiten von 6, 8 und 10 ms-1 abgeblasen. Die abgetragene Sedimentmenge als auch die PM- Emissionen wurden am Ende der Messstrecke erfasst. Die Emissionspotentiale der Substrate nahmen in folgender Reihung ab: homogene Kohle &gt; homogene Feinsande &gt; heterogene Feinsande &gt; heterogene Grobsande &gt; heterogene (faserige) Kohle und lagen in den Bereichen 475 µgg-1 bis 22 µgg-1. Die Befeuchtung der sandigen Substrate auf ca. 2 M% erbrachte eine Reduzierung der PM-Emissionen um 95%, die der Kohle um 45%. Für die durch Winderosion ausgelösten PM-Emissionen ergab sich eine andere Reihenfolge der sandigen Substrate: Kohle &gt; heterogene Feinsande &gt; heterogene Grobsande &gt; homogene Feinsande. Hier wurden vor allem durch den Impakt saltierender Sandkörner Staubpartikel freigesetzt. Für jede der Windgeschwindigkeiten ergab sich über die Zeit eine maximale Abtrags- und PM-Emissionsrate. Wurde diese erreicht, blieb die Oberfläche stabil und es erfolgten keine weiteren PM-Emissionen