Elusive relationships within order fabales: phylogenetic analyses using matK and rbcL sequence data

The order Fabales, including Leguminosae, Polygalaceae, Quillajaceae and Surianaceae, represents a novel hypothesis emerging from angiosperm molecular phylogenies. Despite good support for the order, molecular studies to date have suggested contradictory, poorly supported interfamilial relationships. Our reappraisal of relationships within Fabales addresses past taxon sampling deficiencies, and employs parsimony and Bayesian approaches using sequences from the plastid regions rbcL (166 spp.) and matK (78 spp.). Five alternative hypotheses for interfamilial relationships within Fabales were recovered. The Shimodaira-Hasegawa test found the likelihood of a resolved topology significantly higher than the one calculated for a polytomy, but did not favour any of the alternative hypotheses of relationship within Fabales. In the light of the morphological evidence available and the comparative behavior of rbcL and matK, the topology recovering Polygalaceae as sister to the rest of the order Fabales with Leguminosae more closely related to Quillajaceae + Surianaceae, is considered the most likely hypothesis of interfamilial relationships of the order. Dating of selected crown clades in the Fabales phylogeny using penalized likelihood suggests rapid radiation of the Leguminosae, Polygalaceae, and (Quillajaceae + Surianaceae) crown clades.

A phylogenetic reappraisal of the Peltophorum group (Caesalpinieae: Leguminosae) based on the chloroplast trnL-F, rbcL and rps16 sequence data

The monophyly of the Peltophorum group, one of nine informal groups recognized by Polhill in the Caesalpinieae, was tested using sequence data from the trnL-F, rbcL, and rps16 regions of the chloroplast genome. Exemplars were included from all 16 genera of the Peltophorum group, and from 15 genera representing seven of the other eight informal groups in the tribe. The data were analyzed separately and in combined analyses using parsimony and Bayesian methods. The analysis method had little effect on the topology of well-supported relationships. The molecular data recovered a generally well-supported phylogeny with many intergeneric relationships resolved. Results show that the Peltophorum group as currently delimited is polyphyletic, but that eight genera plus one undescribed genus form a core Peltophorum group, which is referred to here as the Peltophorum group sensu stricto. These genera are Bussea, Conzattia, Colvillea, Delonix, Heteroflorum (inedit.), Lemuropisum, Parkinsonia, Peltophorum, and Schizolobium. The remaining eight genera of the Peltophorum group s.l. are distributed across the Caesalpinieae. Morphological support for the redelimited Peltophorum group and the other recovered clades was assessed, and no unique synapomorphy was found for the Peltophorum group s.s. A proposal for the reclassification of the Peltophorum group s.l. is presented.

Phylogenetic systematics of Sternbergia (Amaryllidaceae) based on plastid and ITS sequence data

The phylogenetics of Sternbergia (Amaryllidaceae) were studied using DNA sequences of the plastid ndhF and matK genes and nuclear internal transcribed spacer (ITS) ribosomal region for 38, 37 and 32 ingroup and outgroup accessions, respectively. All members of Sternbergia were represented by at least one accession, except S. minoica and S. schubertii, with additional taxa from Narcissus and Pancratium serving as principal outgroups. Sternbergia was resolved and supported as sister to Narcissus and composed of two primary subclades: S. colchiciflora sister to S. vernalis, S. candida and S. clusiana, with this clade in turn sister to S. lutea and its allies in both Bayesian and bootstrap analyses. A clear relationship between the two vernal flowering members of the genus was recovered, supporting the hypothesis of a single origin of vernal flowering in Sternbergia. However, in the S. lutea complex, the DNA markers examined did not offer sufficient resolving power to separate taxa, providing some support for the idea that S. sicula and S. greuteriana are conspecific with S. lutea

Multi-target out-of-sequence data association: tracking using graphical models

When performing data fusion, one often measures where targets were and then wishes to deduce where targets currently are. There has been recent research on the processing of such out-of-sequence data. This research has culminated in the development of a number of algorithms for solving the associated tracking problem. This paper reviews these different approaches in a common Bayesian framework and proposes an architecture that orthogonalises the data association and out-of-sequence problems such that any combination of solutions to these two problems can be used together. The emphasis is not on advocating one approach over another on the basis of computational expense, but rather on understanding the relationships among the algorithms so that any approximations made are explicit. Results for a multi-sensor scenario involving out-of-sequence data association are used to illustrate the utility of this approach in a specific context.

Multi-target out-of-sequence data association

In data fusion systems, one often encounters measurements of past target locations and then wishes to deduce where the targets are currently located. Recent research on the processing of such out-of-sequence data has culminated in the development of a number of algorithms for solving the associated tracking problem. This paper reviews these different approaches in a common Bayesian framework and proposes an architecture that orthogonalises the data association and out-of-sequence problems such that any combination of solutions to these two problems can be used together. The emphasis is not on advocating one approach over another on the basis of computational expense, but rather on understanding the relationships between the algorithms so that any approximations made are explicit.

A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Population subdivision and molecular sequence variation: Theory and analysis of Drosophila ananassae data

Population subdivision complicates analysis of molecular variation. Even if neutrality is assumed, three evolutionary forces need to be considered: migration, mutation, and drift. Simplification can be achieved by assuming that the process of migration among and drift within subpopulations is occurring fast compared to Mutation and drift in the entire population. This allows a two-step approach in the analysis: (i) analysis of population subdivision and (ii) analysis of molecular variation in the migrant pool. We model population subdivision using an infinite island model, where we allow the migration/drift parameter Theta to vary among populations. Thus, central and peripheral populations can be differentiated. For inference of Theta, we use a coalescence approach, implemented via a Markov chain Monte Carlo (MCMC) integration method that allows estimation of allele frequencies in the migrant pool. The second step of this approach (analysis of molecular variation in the migrant pool) uses the estimated allele frequencies in the migrant pool for the study of molecular variation. We apply this method to a Drosophila ananassae sequence data set. We find little indication of isolation by distance, but large differences in the migration parameter among populations. The population as a whole seems to be expanding. A population from Bogor (Java, Indonesia) shows the highest variation and seems closest to the species center.

Extreme host plant conservatism during at least 20 million years of host plant pursuit by oak gallwasps

Diversification of insect herbivores is often associated with coevolution between plant toxins and insect countermeasures, resulting in a specificity that restricts host plant shifts. Gall inducers, however, bypass plant toxins and the factors influencing host plant associations in these specialized herbivores remain unclear. We reconstructed the evolution of host plant associations in Western Palaearctic oak gallwasps (Cynipidae: Cynipini), a species-rich lineage of specialist herbivores on oak (Quercus). (1) Bayesian analyses of sequence data for three genes revealed extreme host plant conservatism, with inferred shifts between major oak lineages (sections Cerris and Quercus) closely matching the minimum required to explain observed diversity. It thus appears that the coevolutionary demands of gall induction constrain host plant shifts, both in cases of mutualism (e.g., fig wasps, yucca moths) and parasitism (oak gallwasps). (2) Shifts between oak sections occurred independently in sexual and asexual generations of the gallwasp lifecycle, implying that these can evolve independently. (3) Western Palaearctic gallwasps associated with sections Cerris and Quercus diverged at least 20 million years ago (mya), prior to the arrival of oaks in the Western Palaearctic from Asia 5-7 mya. This implies an Asian origin for Western Palaearctic gallwasps, with independent westwards range expansion by multiple lineages.

Speciation and bursts of evolution

A longstanding debate in evolutionary biology concerns whether species diverge gradually through time or by rapid punctuational bursts at the time of speciation. The theory of punctuated equilibrium states that evolutionary change is characterised by short periods of rapid evolution followed by longer periods of stasis in which no change occurs. Despite years of work seeking evidence for punctuational change in the fossil record, the theory remains contentious. Further there is little consensus as to the size of the contribution of punctuational changes to overall evolutionary divergence. Here we review recent developments which show that punctuational evolution is common and widespread in gene sequence data.

DNA barcoding of a large genus, Aspalathus L. (Fabaceae)

The DNA barcode potential of three regions (the nuclear ribosomal ITS and the plastid psbA-trnH and trnT-trnL intergenic spacers) was investigated for the plant genus Aspalathus L. (Fabaceac: Crotalarieae). Aspalathus is a large genus (278 species) that revealed low levels of DNA variation in phylogenetic studies. In a 51-species dataset for the psbA-trnH and ITS regions, 45%, and 16% of sequences respectively were identical to the sequence of at least one other species, with two species undiscriminated even when the two regions were combined. In contrast, trnT-trnL, discriminated between all species in this dataset. In a larger ITS and trnT-trnL dataset. including a further 82 species. 7 species in five pairwise comparisons remained Undiscriminated when the two regions were combined. Four of the five pairs of species not discriminated by sequence data were readily distinguished using a combination of qualitative and quantitative morphological data. The difficulty of barcoding in this group is increased by the presence of intraspecific variation in all three regions studied. In the case of psbA-trnH, three intraspecific samples had a sequence identical to at least one other species. Overall, psbA-trnH. currently a candidate for plant barcoding, was the least discriminatory region in our study.

The Fox genes of Branchiostoma floridae

The Fox genes are united by encoding a fork head domain, a deoxyribonucleic acid (DNA)-binding domain of the winged-helix type that marks these genes as encoding transcription factors. Vertebrate Fox genes are classified into 23 subclasses named from FoxA to FoxS. We have surveyed the genome of the amphioxus Branchiostoma floridae, identifying 32 distinct Fox genes representing 21 of these 23 subclasses. The missing subclasses, FoxR and FoxS, are specific to vertebrates, and in addition, B. floridae has one further group, FoxAB, that is not found in vertebrates. Hence, we conclude B. floridae has maintained a high level of Fox gene diversity. Expressed sequence tag and complementary DNA sequence data support the expression of 23 genes. Several linkages between B. floridae Fox genes were noted, including some that have evolved relatively recently via tandem duplication in the amphioxus lineage and others that are more ancient.

Novel high-resolution characterization of ancient DNA reveals C > U-type base modification events as the sole cause of post mortem miscoding lesions

Ancient DNA (aDNA) research has long depended on the power of PCR to amplify trace amounts of surviving genetic material from preserved specimens. While PCR permits specific loci to be targeted and amplified, in many ways it can be intrinsically unsuited to damaged and degraded aDNA templates. PCR amplification of aDNA can produce highly-skewed distributions with significant contributions from miscoding lesion damage and non-authentic sequence artefacts. As traditional PCR-based approaches have been unable to fully resolve the molecular nature of aDNA damage over many years, we have developed a novel single primer extension (SPEX)-based approach to generate more accurate sequence information. SPEX targets selected template strands at defined loci and can generate a quantifiable redundancy of coverage; providing new insights into the molecular nature of aDNA damage and fragmentation. SPEX sequence data reveals inherent limitations in both traditional and metagenomic PCR-based approaches to aDNA, which can make current damage analyses and correct genotyping of ancient specimens problematic. In contrast to previous aDNA studies, SPEX provides strong quantitative evidence that C U-type base modifications are the sole cause of authentic endogenous damage-derived miscoding lesions. This new approach could allow ancient specimens to be genotyped with unprecedented accuracy.

Diversification of a large genus in a continental biodiversity hotspot: temporal and spatial origin of Muraltia (Polygalaceae) in the Cape of South Africa

Phylogenetic relationships in the largely South African genus Muraltia (Polygalaceae) are assessed based on DNA sequence data (nuclear ribosomal ITS, plastid atpB-rbcL spacer, trnL intron, and trnL-F spacer) for 73 of the 117 currently recognized species in the genus. The previously recognised subgenus Muraltia is monophyletic, but the South African endemic genus Nylandtia is embedded in Muraltia subgenus Psiloclada. Subgenus Muraltia is found to be sister to subgenus Psiloclada. Estimates show the beginning of diversification of the two subgenera in the early Miocene (Psiloclada, 19.3+/-3.4 Ma; Muraltia, 21.0+/-3.5 Ma) pre-dating the establishment of the Benguela current (intermittent in the middle to late Oligocene and markedly intensifying in the late Miocene), and summer-dry climate in the Cape region. However, the later increase in species numbers is contemporaneous with these climatic phenomena. Results of dispersal-vicariance analyses indicate that major clades in Muraltia diversified from the southwestern and northwestern Cape, where most of the species are found today.

Intercontinental dispersal prior to human translocation revealed in a cryptogenic invasive tree

In this study, complementary species-level and intraspecific phylogenies were used to better circumscribe the original native range and history of translocation of the invasive tree Parkinsonia aculeata. Species-level phylogenies were reconstructed using three chloroplast gene regions, and amplified fragment length polymorphism (AFLP) markers were used to reconstruct the intraspecific phylogeny. Together, these phylogenies revealed the timescale of transcontinental lineage divergence and the likely source of recent introductions of the invasive. The sequence data showed that divergence between North American and Argentinean P. aculeata occurred at least 5.7 million years ago, refuting previous hypotheses of recent dispersal between North and South America. AFLP phylogenies revealed the most likely sources of naturalized populations. The AFLP data also identified putatively introgressed plants, underlining the importance of wide sampling of AFLPs and of comparison with uniparentally inherited marker data when investigating hybridizing groups. Although P. aculeata has generally been considered North American, these data show that the original native range of P. aculeata included South America; recent introductions to Africa and Australia are most likely to have occurred from South American populations.

Phylogeny of Pelargonium (Geraniaceae) based on DNA sequences from three genomes

Phylogenetic hypotheses for the largely South African genus Pelargonium L'Hér. (Geraniaceae) were derived based on DNA sequence data from nuclear, chloroplast and mitochondrial encoded regions. The datasets were unequally represented and comprised cpDNA trnL-F sequences for 152 taxa, nrDNA ITS sequences for 55 taxa, and mtDNA nad1 b/c exons for 51 taxa. Phylogenetic hypotheses derived from the separate three datasets were overall congruent. A single hypothesis synthesising the information in the three datasets was constructed following a total evidence approach and implementing dataset specific stepmatrices in order to correct for substitution biases. Pelargonium was found to consist of five main clades, some with contrasting evolutionary patterns with respect to biogeographic distributions, dispersal capacity, pollination biology and karyological diversification. The five main clades are structured in two (subgeneric) clades that correlate with chromosome size. One of these clades includes a "winter rainfall clade" containing more than 70% of all currently described Pelargonium species, and all restricted to the South African Cape winter rainfall region. Apart from (woody) shrubs and small herbaceous rosette subshrubs, this clade comprises a large "xerophytic" clade including geophytes, stem and leaf succulents, harbouring in total almost half of the genus. This clade is considered to be the result of in situ proliferation, possibly in response to late-Miocene and Pliocene aridification events. Nested within it is a radiation comprising c. 80 species from the geophytic Pelargonium section Hoarea, all characterised by the possession of (a series of) tunicate tubers.