Michael A. Charleston

Insights into the evolutionary history of an emerging livestock pathogen: porcine circovirus 2.

ABSTRACT Porcine circovirus 2 (PCV2) is the primary etiological agent of postweaning multisystemic wasting syndrome (PMWS), one of the most economically important emerging swine diseases worldwide. Virulent PCV2 was first identified following nearly simultaneous outbreaks of PMWS in North America and Europe in the 1990s and has since achieved global distribution. However, the processes responsible for the emergence and spread of PCV2 remain poorly understood. Here, phylogenetic and cophylogenetic inferences were utilized to address key questions on the time scale, processes, and geographic diffusion of emerging PCV2. The results of these analyses suggest that the two genotypes of PCV2 (PCV2a and PCV2b) are likely to have emerged from a common ancestor approximately 100 years ago and have been on independent evolutionary trajectories since that time, despite cocirculating in the same host species and geographic regions. The patterns of geographic movement of PCV2 that we recovered appear to mimic those of the global pig trade and suggest that the movement of asymptomatic animals is likely to have facilitated the rapid spread of virulent PCV2 around the globe. We further estimated the rate of nucleotide substitution for PCV2 to be on the order of 1.2 × 10−3 substitutions/site/year, the highest yet recorded for a single-stranded DNA virus. This high rate of evolution may allow PCV2 to maintain evolutionary dynamics closer to those of single-stranded RNA viruses than to those of double-stranded DNA viruses, further facilitating the rapid emergence of PCV2 worldwide.

Modelling nutritional interactions: from individuals to communities

Nutrient acquisition is a major context for ecological interactions among species but ecologists and nutritionists have developed theory in isolation from each other. Developments in agent-based modelling, state-space modelling of nutrition and multi-scale modelling of landscape ecology provide the components for a new synthesis in nutritional ecology linking the nutritional biology of individual organisms to population- and community-level processes across multiple scales within an evolutionary context. We review the core elements for such a synthesis and set out the principles for a generic modelling framework that could be used to test specific ecological hypotheses.

Differential variability analysis of gene expression and its application to human diseases

Motivation: Current microarray analyses focus on identifying sets of genes that are differentially expressed (DE) or differentially coexpressed (DC) in different biological states (e.g. diseased versus non-diseased). We observed that in many human diseases, some genes have a significantincrease or decrease in expression variability (variance). Asthese observed changes in expression variability may be caused by alteration of the underlying expression dynamics, such differential variability (DV) patterns are also biologically interesting. Results: Here we propose a novel analysis for changes in gene expression variability between groups of amples, which we call differential variability analysis. We introduce the concept of differential variability (DV), and present a simple procedure for identifying DV genes from microarray data. Our procedure is evaluated with simulated and real microarray datasets. The effect of data preprocessing methods on identification of DV gene is investigated. The biological significance of DV analysis is demonstrated with four human disease datasets. The relationships among DV, DE and DC genes are investigated. The results suggest that changes in expression variability are associated with changes in coexpression pattern, which imply that DV is not merely stochastic noise, but informative signal. Availability: The R source code for differential variability analysis is available from the contact authors upon request. Contact: joshua@it.usyd.edu.au; mcharleston@it.usyd.edu.au

On the Computational Complexity of the Reticulate Cophylogeny Reconstruction Problem

The cophylogeny reconstruction problem is that of finding minimal cost explanations of differences between evolutionary histories of ecologically linked groups of biological organisms. We present a proof that shows that the general problem of reconciling evolutionary histories is NP-complete and provide a sharp boundary where this intractability begins. We also show that a related problem, that of finding Pareto optimal solutions, is NP-hard. As a byproduct of our results, we give a framework by which meta-heuristics can be applied to find good solutions to this problem.

Traversing the tangle: algorithms and applications for cophylogenetic studies

Cophylogenetic analysis supposes that two or more phylogenetic trees for linked groups have been constructed, and explores the relationships the trees have with each other. These types of analyses are most commonly used to assess relationships between hosts and their parasites, however the methodology can also be applied to diverse types of problems such as an examination of the phylogenies of genes with respect to those of organisms or those of geographic areas and the organisms that reside there. The working hypothesis is that the trees are correct, though sometimes attempts are made to take into account their uncertainty. Cophylogeny is computationally hard: that is, there are no known fast methods to compute relationships among such trees for any but the simplest of models. A review of methodology that has been developed to examine cophylogenetic relationships is presented and a brief discussion of some medically relevant examples is given.

Fungal identification using a Bayesian classifier and the Warcup training set of internal transcribed spacer sequences

Fungi are key organisms in many ecological processes and communities. Rapid and low cost surveys of the fungal members of a community can be undertaken by isolating and sequencing a taxonomically informative genomic region, such as the ITS (internal transcribed spacer), from DNA extracted from a metagenomic sample, and then classifying these sequences to determine which organisms are present. This paper announces the availability of the Warcup ITS training set and shows how it can be used with the Ribosomal Database Project (RDP) Bayesian Classifier to rapidly and accurately identify fungi using ITS sequences. The classifications can be down to species level and use conventional literature-based mycological nomenclature and taxonomic assignments.

Nutritional ecology beyond the individual: a conceptual framework for integrating nutrition and social interactions.

Over recent years, modelling approaches from nutritional ecology (known as Nutritional Geometry) have been increasingly used to describe how animals and some other organisms select foods and eat them in appropriate amounts in order to maintain a balanced nutritional state maximising fitness. These nutritional strategies profoundly affect the physiology, behaviour and performance of individuals, which in turn impact their social interactions within groups and societies. Here, we present a conceptual framework to study the role of nutrition as a major ecological factor influencing the development and maintenance of social life. We first illustrate some of the mechanisms by which nutritional differences among individuals mediate social interactions in a broad range of species and ecological contexts. We then explain how studying individual- and collective-level nutrition in a common conceptual framework derived from Nutritional Geometry can bring new fundamental insights into the mechanisms and evolution of social interactions, using a combination of simulation models and manipulative experiments.

Spectrum: spectral analysis of phylogenetic data.

MOTIVATION Spectrum is a new Macintosh and Microsoft Windows program designed to read in phylogenetic four-state or binary data in NEXUS format, and output the corresponding bipartition spectra. It can be used to find the tree whose expected spectrum is closest to the observed spectrum (the closest tree, Hendy, Discr, Math., 96, 51-58,. 1991).

The Effects of Sequence Length, Tree Topology, and Number of Taxa on the Performance of Phylogenetic Methods

Simulations were used to study the performance of several character-based and distance-based phylogenetic methods in obtaining the correct tree from pseudo-randomly generated input data. The study included all the topologies of unrooted binary trees with from 4 to 10 pendant vertices (taxa) inclusive. The length of the character sequences used ranged from 10 to 10(5) characters exponentially. The methods studied include Closest Tree, Compatibility, Lis method, Maximum Parsimony, Neighbor-joining, Neighborliness, and UPGMA. We also provide a modification to Lis method (SimpLi) which is consistent with additive data. We give estimations of the sequence lengths required for given confidence in the output of these methods under the assumptions of molecular evolution used in this study. A notation for characterizing all tree topologies is described. We show that when the number of taxa, the maximum path length, and the minimum edge length are held constant, there it little but significant dependence of the performance of the methods on the tree topology. We show that those methods that are consistent with the model used perform similarly, whereas the inconsistent methods, UPGMA and Lis method, perform very poorly.

Phylogenetic Codivergence Supports Coevolution of Mimetic Heliconius Butterflies

The unpalatable and warning-patterned butterflies Heliconius erato and Heliconius melpomene provide the best studied example of mutualistic Müllerian mimicry, thought–but rarely demonstrated–to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of H. erato and H. melpomene, and this was initially hailed as one of the most striking known cases of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of H. erato and H. melpomene, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of H. erato and H. melpomene. Cophylogenetic historical reconstructions support repeated codivergence of mimetic populations, from the base of the sampled radiations. Pairwise distance correlation tests, based on our coalescent analyses plus recently published AFLP and wing colour pattern gene data, also suggest that the phylogenies of H. erato and H. melpomene show significant topological congruence. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of H. erato and H. melpomene occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. Our results suggest that differences in within-species genetic divergence are the result of a greater overall effective population size for H. erato relative to H. melpomene and do not imply incongruence in the timing of their phylogenetic radiations. Repeated codivergence between Müllerian co-mimics, predicted to exert mutual selection pressures, strongly suggests coevolution. Our results therefore support a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change.