Nuclear and plastid markers reveal the persistence of genetic identity: A new perspective on the evolutionary history of Petunia exserta

Recently divergent species that can hybridize are ideal models for investigating the genetic exchanges that can occur while preserving the species boundaries. Petunia exserta is an endemic species from a very limited and specific area that grows exclusively in rocky shelters. These shaded spots are an inhospitable habitat for all other Petunia species, including the closely related and widely distributed species P. axillaris. Individuals with intermediate morphologic characteristics have been found near the rocky shelters and were believed to be putative hybrids between P. exserta and P. axillaris, suggesting a situation where Petunia exserta is losing its genetic identity. In the current study, we analyzed the plastid intergenic spacers trnS/trnG and trnH/psbA and six nuclear CAPS markers in a large sampling design of both species to understand the evolutionary process occurring in this biological system. Bayesian clustering methods, cpDNA haplotype networks, genetic diversity statistics, and coalescence-based analyses support a scenario where hybridization occurs while two genetic clusters corresponding to two species are maintained. Our results reinforce the importance of coupling differentially inherited markers with an extensive geographic sample to assess the evolutionary dynamics of recently diverged species that can hybridize. (C) 2013 Elsevier Inc. All rights reserved.

Molecular phylogenetic analyses identify Alpine differentiation and dysploid chromosome number changes as major forces for the evolution of the European endemic Phyteuma (Campanulaceae)

Phyteuma is a chromosomally and ecologically diverse vascular plant genus and constitutes an excellent system for studying both the role of chromosomal change for species diversification and the evolution of high-mountain biota. This kind of research is, however, hampered by the lack of a sound phylogenetic framework exacerbated by the notoriously low predictive power of traditional taxonomy with respect to phylogenetic relationships in Campanulaceae. Based on a comprehensive taxon sampling and analyses of nuclear and plastid sequence and AFLP fingerprint data, Phyteuma is confirmed as a monophyletic group sister to the monotypic Physoplexis, which is in line with their peculiar flower morphologies. Within Phyteuma two clades, largely corresponding to previously recognized sections, are consistently found. The traditional circumscription of taxonomic series is largely rejected. Whereas distinctness of the currently recognized species is mostly corroborated, some interspecific relationships remain ambiguous due to incongruences between nuclear and plastid data. Major forces for diversification and evolution of Phyteuma are descending dysploidy (i.e., a decrease in chromosome base number) as well as allopatric and ecological differentiation within the Alps, the genus' center of species diversity. (C) 2013 Elsevier Inc. All rights reserved.

‎101011 manuscripts: Perl and the Chronicle of Matthew of Edessa‎

The Chronicle of Matthew of Edessa is a twelfth-century Armenian history that survives in 43 manuscripts, held in ten libraries in eight countries. My task is to create a definitive text that is based on all of them. I will talk about the problems of medieval text editing, the ways in which Perl and phylogenetics have come to my rescue, and show a few pretty pictures of manuscripts

Complete genome of a Puumala virus strain from Central Europe

Puumala virus (PUUV) is one of the predominant hantavirus species in Europe causing mild to moderate cases of haemorrhagic fever with renal syndrome. Parts of Lower Saxony in north-western Germany are endemic for PUUV infections. In this study, the complete PUUV genome sequence of a bank vole-derived tissue sample from the 2007 outbreak was determined by a combined primer-walking and RNA ligation strategy. The S, M and L genome segments were 1,828, 3,680 and 6,550 nucleotides in length, respectively. Sliding-window analyses of the nucleotide sequences of all available complete PUUV genomes indicated a non-homogenous distribution of variability with hypervariable regions located at the 3′-ends of the S and M segments. The overall similarity of the coding genome regions to the other PUUV strains ranged between 80.1 and 84.7 % at the level of the nucleotide sequence and between 89.5 and 98.1 % for the deduced amino acid sequences. In comparison to the phylogenetic trees of the complete coding sequences, trees based on partial segments revealed a general drop in phylogenetic support and a lower resolution. The Astrup strain S and M segment sequences showed the highest similarity to sequences of strains from geographically close sites in the Osnabrück Hills region. In conclusion, a primer-walking-mediated strategy resulted in the determination of the first complete nucleotide sequence of a PUUV strain from Central Europe. Different levels of variability along the genome provide the opportunity to choose regions for analyses according to the particular research question, e.g., large-scale phylogenetics or within-host evolution.

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.

Phylogenetic analyses and toxigenic profiles of Fusarium equiseti and Fusarium acuminatum isolated from cereals from Southern Europe

Fusarium equiseti and Fusarium acuminatum are toxigenic species that contaminate cereal crops from diverse climatic regions. They are common in Spanish cereals. The information available on their phylogenetics and toxigenic profiles is, however, insufficient to assist risk evaluation. In this work, phylogenetic analyses were performed using partial sequences of the translation elongation factor gene (EF-1a) of F. equiseti and F. acuminatum strains isolated from barley and wheat from Spain and other countries. The Northern and Southern European F. equiseti strains largely separated into two phylogenetically distinct clusters. This suggests the existence of two distinct populations within this species, explaining its presence in these regions of markedly different climate. Production of type A and B trichothecenes by the Spanish strains, examined in wheat cultures using a multitoxin analytical method, indicated that F. equiseti could produce deoxynivalenol and nivalenol and other trichothecenes, at concentrations that might represent a significant risk of toxin contamination for Southern European cereals. F. acuminatum showed low intraspecific genetic variability and 58% of the strains could produce deoxynivalenol at low level. Neither species was found to produce T-2 or HT-2 toxins. The present results provide important phylogenetic and toxigenic information essential for the accurate prediction of toxigenic risk.

Resolving taxonomic uncertainty in vulnerable elasmobranchs : are the Madeira skate (Raja maderensis) and the thornback ray (Raja clavata) distinct species?

We are grateful to all those who helped with sample collection. This includes the MEDITS survey programme (IEO Mallorca) for Mediterranean samples. Portugal mainland samples were collected under the EU Data Collection Framework (DCF, PNAB). Azores specimens from the Department of Oceanography and Fisheries (DOP) of the University of the Azores (UAc) were collected under the project DEMERSAIS “Monitorização das espécies demersais dos Açores” financed by the Azorean government, and the project DEECON “Unravelling population connectivity for sustainable fisheries in the Deep Sea” project approved by the European Science Foundation (ESF) under the EUROCORES programme (proposal No 06-EuroDEEP-FP-008 & SFRH-EuroDEEP/0002/2007). This study was funded by the UK Natural Environment Research Council (NERC) Oceans 2025 Strategic Research Programme Theme 6 (Science for Sustainable Marine Resources). REB was supported by the Fisheries Society of the British Isles (FSBI), BSP was funded by the Fundação para a Ciência e a Tecnologia, SFRH/BPD/72351/2010 and JL was supported by The Alasdair Downes Marine Conservation Fund.

Discovery of a new source of rifamycin antibiotics in marine sponge actinobacteria by phylogenetic prediction

Phylogenetic analysis of the ketosynthase (KS) gene sequences of marine sponge-derived Salinispora strains of actinobacteria indicated that the polyketide synthase (PKS) gene sequence most closely related to that of Salinispora was the rifamycin B synthase of Amycolatopsis mediterranei. This result was not expected from taxonomic species tree phylogenetics using 16S rRNA sequences. From the PKS sequence data generated from our sponge-derived Salinispora strains, we predicted that such strains might synthesize rifamycin-like compounds. Liquid chromatography-tandem mass spectrometry (LC/MS/MS) analysis was applied to one sponge-derived Salinispora strain to test the hypothesis of rifamycin synthesis. The analysis reported here demonstrates that this Salinispora isolate does produce compounds of the rifamycin class, including rifamycin B and rifamycin SV. A rifamycin-specific KS primer set was designed, and that primer set increased the number of rifamycin-positive strains detected by PCR screening relative to the number detectable using a conserved KS-specific set. Thus, the Salinispora group of actinobacteria represents a potential new source of rifamycins outside the genus Amycolatopsis and the first recorded source of rifamycins from marine bacteria.

Diversity of polyketide synthase genes from bacteria associated with the marine sponge Pseudoceratina clavata: culture-dependent and culture-independent approaches

Diverse ketosynthase (KS) genes were retrieved from the microbial community associated with the Great Barrier Reef sponge Pseudoceratina clavata. Bacterial isolation and metagenomic approaches were employed. Phylogenetic analysis of 16S rRNA of culturable sponge-associated bacterial communities comprised eight groups over four phyla. Ten KS domains were amplified from four genera of isolates and phylogenetics demonstrated that these KS domains were located in three clusters (actinobacterial, cyanobacterial and trans-AT type). Metagenomic DNA of the sponge microbial community was extracted to explore community KS genes by two approaches: direct amplification of KS domains and construction of fosmid libraries for KS domain screening. Five KS domains were retrieved from polymerase chain reaction (PCR) amplification using sponge metagenome DNA as template and five fosmid clones containing KS domains found using multiplex PCR screening. Analysis of selected polyketide synthase (PKS) from one fosmid showed that the PKS consists of two modules. Open reading frames located up- and downstream of the PKS displayed similarity with membrane synthesis-related proteins such as cardiolipin synthase. Metagenome approaches did not detect KS domains found in sponge isolates. All KS domains from both metagenome approaches formed a single cluster with KS domains originating from metagenomes derived from other sponge species from other geographical regions.

Molecular systematics of the Philippine malaria vector Anopheles flavirostris

Allozyme and molecular sequence data from the malaria vector Anopheles flavirostris (Ludlow) (Diptera: Culicidae) were analysed from 34 sites throughout the Philippines, including the type locality, to test the hypothesis that this taxon is a single panmictic species. A finer-scaled allozyme study, of mainly Luzon samples, revealed no fixed genetic differences in sympatric sites and only low levels of variation. We obtained data from partial sequences for the internal transcribed spacer 2 (ITS2) (483 bp), the third domain (D3) (330 bp) of the 28S ribosomal DNA subunit and cytochrome c oxidase subunit I (COI) of mitochondrial DNA (261 bp). No sequence variation was observed for ITS2, only a one base pair difference was observed between Philippine and Indonesian D3 sequences and An. flavirostris sequences were unique, confirming their diagnostic value for this taxon. Sixteen COI haplotypes were identified, giving 25 parsimony informative sites. Neighbour-Joining, Maximum Parsimony, Maximum Likelihood and Bayesian phylogenetic analysis of COI sequences for An. flavirostris and outgroup taxa revealed strong branch support for the monophyly of An. flavirostris, thus confirming that Philippine populations of this taxon comprise a single separate species within the Minimus Subgroup of the Funestus Group. Variation in the behaviour of An. flavirostris is likely to be intraspecific rather than interspecific in origin. © 2006 The Royal Entomological Society.

Searching for convergence in phylogenetic Markov chain Monte Carlo

Markov chain Monte Carlo (MCMC) is a methodology that is gaining widespread use in the phylogenetics community and is central to phylogenetic software packages such as MrBayes. An important issue for users of MCMC methods is how to select appropriate values for adjustable parameters such as the length of the Markov chain or chains, the sampling density, the proposal mechanism, and, if Metropolis-coupled MCMC is being used, the number of heated chains and their temperatures. Although some parameter settings have been examined in detail in the literature, others are frequently chosen with more regard to computational time or personal experience with other data sets. Such choices may lead to inadequate sampling of tree space or an inefficient use of computational resources. We performed a detailed study of convergence and mixing for 70 randomly selected, putatively orthologous protein sets with different sizes and taxonomic compositions. Replicated runs from multiple random starting points permit a more rigorous assessment of convergence, and we developed two novel statistics, delta and epsilon, for this purpose. Although likelihood values invariably stabilized quickly, adequate sampling of the posterior distribution of tree topologies took considerably longer. Our results suggest that multimodality is common for data sets with 30 or more taxa and that this results in slow convergence and mixing. However, we also found that the pragmatic approach of combining data from several short, replicated runs into a metachain to estimate bipartition posterior probabilities provided good approximations, and that such estimates were no worse in approximating a reference posterior distribution than those obtained using a single long run of the same length as the metachain. Precision appears to be best when heated Markov chains have low temperatures, whereas chains with high temperatures appear to sample trees with high posterior probabilities only rarely. [Bayesian phylogenetic inference; heating parameter; Markov chain Monte Carlo; replicated chains.]

Pattern-based phylogenetic distance estimation and tree reconstruction

We have developed an alignment-free method that calculates phylogenetic distances using a maximum-likelihood approach for a model of sequence change on patterns that are discovered in unaligned sequences. To evaluate the phylogenetic accuracy of our method, and to conduct a comprehensive comparison of existing alignment-free methods (freely available as Python package decaf+py at http://www.bioinformatics.org.au), we have created a data set of reference trees covering a wide range of phylogenetic distances. Amino acid sequences were evolved along the trees and input to the tested methods; from their calculated distances we infered trees whose topologies we compared to the reference trees. We find our pattern-based method statistically superior to all other tested alignment-free methods. We also demonstrate the general advantage of alignment-free methods over an approach based on automated alignments when sequences violate the assumption of collinearity. Similarly, we compare methods on empirical data from an existing alignment benchmark set that we used to derive reference distances and trees. Our pattern-based approach yields distances that show a linear relationship to reference distances over a substantially longer range than other alignment-free methods. The pattern-based approach outperforms alignment-free methods and its phylogenetic accuracy is statistically indistinguishable from alignment-based distances.