Rapid radiation in spiny lobsters (Palinurus spp) as revealed by classic and ABC methods using mtDNA and microsatellite data

Background: Molecular tools may help to uncover closely related and still diverging species from a wide variety of taxa and provide insight into the mechanisms, pace and geography of marine speciation. There is a certain controversy on the phylogeography and speciation modes of species-groups with an Eastern Atlantic-Western Indian Ocean distribution, with previous studies suggesting that older events (Miocene) and/or more recent (Pleistocene) oceanographic processes could have influenced the phylogeny of marine taxa. The spiny lobster genus Palinurus allows for testing among speciation hypotheses, since it has a particular distribution with two groups of three species each in the Northeastern Atlantic (P. elephas, P. mauritanicus and P. charlestoni) and Southeastern Atlantic and Southwestern Indian Oceans (P. gilchristi, P. delagoae and P. barbarae). In the present study, we obtain a more complete understanding of the phylogenetic relationships among these species through a combined dataset with both nuclear and mitochondrial markers, by testing alternative hypotheses on both the mutation rate and tree topology under the recently developed approximate Bayesian computation (ABC) methods. Results: Our analyses support a North-to-South speciation pattern in Palinurus with all the South-African species forming a monophyletic clade nested within the Northern Hemisphere species. Coalescent-based ABC methods allowed us to reject the previously proposed hypothesis of a Middle Miocene speciation event related with the closure of the Tethyan Seaway. Instead, divergence times obtained for Palinurus species using the combined mtDNA-microsatellite dataset and standard mutation rates for mtDNA agree with known glaciation-related processes occurring during the last 2 my. Conclusion: The Palinurus speciation pattern is a typical example of a series of rapid speciation events occurring within a group, with very short branches separating different species. Our results support the hypothesis that recent climate change-related oceanographic processes have influenced the phylogeny of marine taxa, with most Palinurus species originating during the last two million years. The present study highlights the value of new coalescent-based statistical methods such as ABC for testing different speciation hypotheses using molecular data.

The comparative toxicity to soil invertebrates of natural chemicals and their synthetic analogues

The introduction of Registration, Evaluation and Authorisation of Chemicals (REACH), requires companies to register and risk assess all substances produced or imported in volumes of >1 tonne per year. Extrapolation methods which use existing data for estimating the effects of chemicals are attractive to industry, and comparative data are therefore increasingly in demand. Data on natural toxic chemicals could be used for extrapolation methods Such as read-across. To test this hypothesis, the toxicity of natural chemicals and their synthetic analogues were compared using standardised toxicity tests. Two chemical pairs: the napthoquinones, juglone (natural) and 1,4-naphthoquinone (synthetic); and anthraquinones, emodin (natural) and quinizarin (synthetic) were chosen, and their comparative effects on the survival and reproduction of collembolans, earthworms, enchytraeids and predatory mites were assessed. Differences in sensitivity between the species were observed with the predatory mite (Hypoaspis aculeifer) showing the least sensitivity. Within the chemical pairs, toxicity to lethal and sub-lethal endpoints was very similar for the four invertebrate species. The exception was earthworm reproduction, which showed differential sensitivity to the chemicals in both naphthoquinone and anthraquinone pairs. Differences in toxicity identified in the present study may be related to degree of exposure and/or subtle differences in the mode of toxic action for the chemicals and species tested. It may be possible to predict differences by identifying functional groups which infer increased or decreased toxicity in one or other chemical. The development of such techniques would enable the use of read-across from natural to synthetic chemicals for a wider group of compounds. (C) 2009 Elsevier Ltd. All rights reserved.

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

Phylogenetic analysis of the pPT23A plasmid family of Pseudomonas syringae

The pPT23A plasmid family of Pseudomonas syringae contains members that contribute to the ecological and pathogenic fitness of their P. syringae hosts. In an effort to understand the evolution of these plasmids and their hosts, we undertook a comparative analysis of the phylogeny of plasmid genes and that of conserved chromosomal genes from P. syringae. In total, comparative sequence and phylogenetic analyses were done utilizing 47 pPT23A family plasmids (PFPs) from 16 pathovars belonging to six genomospecies. Our results showed that the plasmid replication gene (repA), the only gene currently known to be distributed among all the PFPs, had a phylogeny that was distinct from that of the P. syringae hosts of these plasmids and from those of other individual genes on PFPs. The phylogenies of two housekeeping chromosomal genes, those for DNA gyrase B subunit (gyrB) and primary sigma factor (rpoD), however, were strongly associated with genomospecies of P. syringae. Based on the results from this study, we conclude that the pPT23A plasmid family represents a dynamic genome that is mobile among P. syringae pathovars.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

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.

Comparative expression analysis of transcription factor genes in the endostyle of invertebrate chordates

The endostyle of invertebrate chordates is a pharyngeal organ that is thought to be homologous with the follicular thyroid of vertebrates. Although thyroid-like features such as iodine-concentrating and peroxidase activities are located in the dorsolateral part of both ascidian and amphioxus endostyles, the structural organization and numbers of functional units are different. To estimate phylogenetic relationships of each functional zone with special reference to the evolution of the thyroid, we have investigated, in ascidian and amphioxus, the expression patterns of thyroid-related transcription factors such as TTF-2/MoxE4 and Pax2/5/8, as well as the forkhead transcription factors FoxQ1 and FoxA. Comparative gene expression analyses depicted an overall similarity between ascidians and amphioxus endostyles, while differences in expression patterns of these genes might be specifically related to the addition or elimination of a pair of glandular zones. Expressions of Ci-FoxE and BbFoxE4 suggest that the ancestral FoxE class might have been recruited for the formation of thyroid-like region in a possible common ancestor of chordates. Furthermore, coexpression of FoxE4, Pax2/5/8, and TPO in the dorsolateral part of both ascidian and amphioxus endostyles suggests that genetic basis of the thyroid function was already in place before the vertebrate lineage. (c) 2005 Wiley-Liss, Inc.

Bayesian estimation of ancestral character states on phylogenies

Biologists frequently attempt to infer the character states at ancestral nodes of a phylogeny from the distribution of traits observed in contemporary organisms. Because phylogenies are normally inferences from data, it is desirable to account for the uncertainty in estimates of the tree and its branch lengths when making inferences about ancestral states or other comparative parameters. Here we present a general Bayesian approach for testing comparative hypotheses across statistically justified samples of phylogenies, focusing on the specific issue of reconstructing ancestral states. The method uses Markov chain Monte Carlo techniques for sampling phylogenetic trees and for investigating the parameters of a statistical model of trait evolution. We describe how to combine information about the uncertainty of the phylogeny with uncertainty in the estimate of the ancestral state. Our approach does not constrain the sample of trees only to those that contain the ancestral node or nodes of interest, and we show how to reconstruct ancestral states of uncertain nodes using a most-recent-common-ancestor approach. We illustrate the methods with data on ribonuclease evolution in the Artiodactyla. Software implementing the methods ( BayesMultiState) is available from the authors.

Comparative characterisation by atomic force microscopy and ellipsometry of soft and solid thin films

Ellipsometry and atomic force microscopy (AFM) were used to study the film thickness and the surface roughness of both 'soft' and solid thin films. 'Soft' polymer thin films of polystyrene and poly(styrene-ethylene/butylene-styrene) block copolymer were prepared by spin-coating onto planar silicon wafers. Ellipsometric parameters were fitted by the Cauchy approach using a two-layer model with planar boundaries between the layers. The smooth surfaces of the prepared polymer films were confirmed by AFM. There is good agreement between AFM and ellipsometry in the 80-130 nm thickness range. Semiconductor surfaces (Si) obtained by anisotropic chemical etching were investigated as an example of a randomly rough surface. To define roughness parameters by ellipsometry, the top rough layers were treated as thin films according to the Bruggeman effective medium approximation (BEMA). Surface roughness values measured by AFM and ellipsometry show the same tendency of increasing roughness with increased etching time, although AFM results depend on the used window size. The combined use of both methods appears to offer the most comprehensive route to quantitative surface roughness characterisation of solid films. Copyright (c) 2007 John Wiley & Sons, Ltd.

Corynebacterium suicordis sp. nov., from pigs

Nineteen strains of Gram-positive, non-motile, non-spore-forming, catalase-positive, rod-shaped bacteria isolated from pigs were characterized by using biochemical, molecular chemical and molecular genetic methods. Two distinct groups of organisms were discerned, based on their colonial morphology, CAMP (Christie-Atkins-Munch-Petersen) reaction and numerical profile by using the API Coryne system. The first group (113 strains) gave a doubtful discrimination between Corynebacterium striatum and Corynebacterium amycolatum, whilst the second group (six strains) were identified tentatively as Corynebacterium urealyticum. Comparative 16S rRNA gene sequencing studies demonstrated that all of the isolates belonged phylogenetically to the genus Corynebacterium. The first group of organisms was highly similar to Corynebacterium testudinoris with respect to 16S rRNA gene sequences and physiological characteristics, whereas the remaining six isolates formed a hitherto unknown subline within the genus, associated with a small subcluster of species that included Corynebacterium auriscanis and its close relatives. The unknown Corynebacterium sp. was distinguished readily from these and other species of the genus by biochemical tests. Based on both phenotypic and phylogenetic evidence, it is proposed that the new isolates from pigs should be classified as a novel species, Corynebacterium suicordis sp. nov. The type strain is P81/02(T) (=CECT 5724(T) =CCUG 46963(T)).

Uruburuella suis gen. nov., sp nov., isolated from clinical specimens of pigs

Five strains of an unusual Gram-negative, catalase-positive, oxidase-positive, coccobacillus-shaped bacterium isolated from the lungs and heart of pigs with pneumonia and pericarditis were characterized by phenotypic and molecular genetic methods. On the basis of cellular morphology and biochemical criteria, the isolates were tentatively assigned to the family Neisseriaceae, although they did not appear to correspond to any recognized genus or species. Comparative 16S rRNA gene sequencing showed that the five unidentified strains were phylogenetically highly related to each other and represent a hitherto unknown subline within the family Neisseriaceae. On the basis of both phenotypic and phylogenetic evidence, it is proposed that the unknown isolates from pigs be classified as a novel genus and species within the family Neisseriaceae, for which the name Uruburuella suis gen. nov., sp. nov. is proposed. The type strain of U. suis is 1258/02(T) (=CCUG 47806(T) =CECT 5685(T)).

Clostridium bolteae sp. nov., isolated from human sources

Seven obligately anaerobic, gram-positive, rod-shaped, spore-forming organisms isolated from human sources were characterized using phenotypic and molecular taxonomic methods. Comparative 16S rRNA gene sequencing showed that the strains were genetically highly related to each other (displaying >99% sequence similarity) and represent a previously unknown sub-line within the Clostridium coccoides rRNA group of organisms. Strains of the unidentified bacterium used carbohydrate as fermentable substrates, producing acetic acid and lactic acid as the major products of glucose metabolism. The closest described species to the novel bacterium corresponded to Clostridium clostridioforme, although a 16S rRNA sequence divergence of 3% demonstrated they represent different species. Genomic DNA-DNA pairing studies confirmed the separateness of the unknown species and Clostridium clostridioforme. Based on phenotypic and phylogenetic evidence, it is therefore proposed that the unknown bacterium, be classified as Clostridium bolteae sp. nov. The type strain of Clostridium bolteae is WAL 16351(T) (= ATCC(T) = BAA-613(T), CCUG(T) = 46953(T)).

Comparative effects of fatty acids on endothelial inflammatory gene expression

Background: Endothelial dysfunction may be related to adverse effects of some dietary fatty acids (FAs). Although in vitro studies have failed to show consistent findings, this may reflect the diverse experimental protocols employed and the limited range of FAs and end points studied. Aims: To investigate the effect of dietary FA type (saturated, monounsaturated, n-6 and n-3 polyunsaturated fatty acids), concentration, incubation time and cell stimulation state, on a broad spectrum of endothelial inflammatory gene expression. Methods: Using human umbilical vein endothelial cells, with and without stimulation (+/- 10 ng/ml TNF alpha), the effects of arachidonic (AA), docosahexaenoic (DHA), eicosapentaenoic (EPA), linoleic (LA), oleic (OA) and palmitic acids (PA) (10, 25 and 100 mu M), on the expression of genes encoding a number of inflammatory proteins and transcription factors were assessed by quantitative real time RT-PCR. Results: Individual FAs differentially affect endothelial inflammatory gene expression in a gene-specific manner. EPA, LA and OA significantly up-regulated MCP-1 gene expression compared to AA (p = 0.001, 0.013, 0.008, respectively) and DHA (p < 0.0005, = 0.004, 0.002, respectively). Furthermore, cell stimulation state and FA incubation time significantly influenced reported FA effects on gene expression. Conclusions: The comparative effects of saturated, monounsaturated, n-6 and n-3 polyunsaturated FAs on endothelial gene expression depend on the specific FA investigated, its length of incubation, cell stimulation state and the gene investigated. These findings may explain existing disparity in the literature. This work was funded by the EC, Framework Programme 6 via the LIPGENE project (FOOD-CT-2003-505944).