Genetic features and molecular epidemiology of Enterococcus faecium isolated in two university hospitals in Brazil

The global emergence of vancomycin-resistant Enterococcus faecium (VREfm) has been characterized by a clonal spread of strains belonging to clonal complex 17 (CC17). Genetic features and clonal relationships of 53 VREfm isolated from patients in 2 hospitals in Ribeirao Preto, Sao Paulo, Brazil, during 2005-2010 were determined as a contribution to the Brazilian evolutionary history of these nosocomial pathogens. All isolates were daptomycin susceptible, vancomycin-resistant, and had the vanA gene. The predominant virulence genes were acm and esp. Only 5 VREfm isolated in 2005-2006 had intact Tn1546, while 81% showed Tn1546 with deleted left extremity and insertion of IS1251 between the vanS and vanH genes. Multilocus sequence typing analysis permitted the identification of 9 different sequence types (STs), with 5 being new ones (656, 657, 658, 659, and 660). Predominant STs were ST412 and ST478, all belonging to CC17, except ST658. This is the first report of the ST78 in Brazil. (c) 2012 Elsevier Inc. All rights reserved.

Trypanosoma livingstonei: a new species from African bats supports the bat seeding hypothesis for the Trypanosoma cruzi clade

Abstract Background Bat trypanosomes have been implicated in the evolutionary history of the T. cruzi clade, which comprises species from a wide geographic and host range in South America, Africa and Europe, including bat-restricted species and the generalist agents of human American trypanosomosis T. cruzi and T. rangeli. Methods Trypanosomes from bats (Rhinolophus landeri and Hipposideros caffer) captured in Mozambique, southeast Africa, were isolated by hemoculture. Barcoding was carried out through the V7V8 region of Small Subunit (SSU) rRNA and Fluorescent Fragment Length barcoding (FFLB). Phylogenetic inferences were based on SSU rRNA, glyceraldehyde phosphate dehydrogenase (gGAPDH) and Spliced Leader (SL) genes. Morphological characterization included light, scanning and transmission electron microscopy. Results New trypanosomes from bats clustered together forming a clade basal to a larger assemblage called the T. cruzi clade. Barcoding, phylogenetic analyses and genetic distances based on SSU rRNA and gGAPDH supported these trypanosomes as a new species, which we named Trypanosoma livingstonei n. sp. The large and highly polymorphic SL gene repeats of this species showed a copy of the 5S ribosomal RNA into the intergenic region. Unique morphological (large and broad blood trypomastigotes compatible to species of the subgenus Megatrypanum and cultures showing highly pleomorphic epimastigotes and long and slender trypomastigotes) and ultrastructural (cytostome and reservosomes) features and growth behaviour (when co-cultivated with HeLa cells at 37°C differentiated into trypomastigotes resembling the blood forms and do not invaded the cells) complemented the description of this species. Conclusion Phylogenetic inferences supported the hypothesis that Trypanosoma livingstonei n. sp. diverged from a common ancestral bat trypanosome that evolved exclusively in Chiroptera or switched at independent opportunities to mammals of several orders forming the clade T. cruzi, hence, providing further support for the bat seeding hypothesis to explain the origin of T. cruzi and T. rangeli.

Desmodus rotundus and Artibeus spp. bats might present distinct rabies virus lineages

In Brazil, bats have been assigned an increasing importance in public health as they are important rabies reservoirs. Phylogenetic studies have shown that rabies virus (RABV) strains from frugivorous bats Artibeus spp. are closely associated to those from the vampire bat Desmodus rotundus, but little is known about the molecular diversity of RABV in Artibeus spp. The N and G genes of RABV isolated from Artibeus spp. and cattle infected by D. rotundus were sequenced, and phylogenetic trees were constructed. The N gene nucleotides tree showed three clusters: one for D. rotundus and two for Artibeus spp. Regarding putative N amino acid-trees, two clusters were formed, one for D. rotundus and another for Artibeus spp. RABV G gene phylogeny supported the distinction between D. rotundus and Artibeus spp. strains. These results show the intricate host relationship of RABV's evolutionary history, and are invaluable for the determination of RABV infection sources.

Trypanosoma livingstonei: a new species from African bats supports the bat seeding hypothesis for the Trypanosoma cruzi clade

BACKGROUND: Bat trypanosomes have been implicated in the evolutionary history of the T. cruzi clade, which comprises species from a wide geographic and host range in South America, Africa and Europe, including bat-restricted species and the generalist agents of human American trypanosomosis T. cruzi and T. rangeli. METHODS: Trypanosomes from bats (Rhinolophus landeri and Hipposideros caffer) captured in Mozambique, southeast Africa, were isolated by hemoculture. Barcoding was carried out through the V7V8 region of Small Subunit (SSU) rRNA and Fluorescent Fragment Length barcoding (FFLB). Phylogenetic inferences were based on SSU rRNA, glyceraldehyde phosphate dehydrogenase (gGAPDH) and Spliced Leader (SL) genes. Morphological characterization included light, scanning and transmission electron microscopy. RESULTS: New trypanosomes from bats clustered together forming a clade basal to a larger assemblage called the T. cruzi clade. Barcoding, phylogenetic analyses and genetic distances based on SSU rRNA and gGAPDH supported these trypanosomes as a new species, which we named Trypanosoma livingstonei n. sp. The large and highly polymorphic SL gene repeats of this species showed a copy of the 5S ribosomal RNA into the intergenic region. Unique morphological (large and broad blood trypomastigotes compatible to species of the subgenus Megatrypanum and cultures showing highly pleomorphic epimastigotes and long and slender trypomastigotes) and ultrastructural (cytostome and reservosomes) features and growth behaviour (when co-cultivated with HeLa cells at 37°C differentiated into trypomastigotes resembling the blood forms and do not invaded the cells) complemented the description of this species. CONCLUSION: Phylogenetic inferences supported the hypothesis that Trypanosoma livingstonei n. sp. diverged from a common ancestral bat trypanosome that evolved exclusively in Chiroptera or switched at independent opportunities to mammals of several orders forming the clade T. cruzi, hence, providing further support for the bat seeding hypothesis to explain the origin of T. cruzi and T. rangeli.

Dengue virus type 3 adaptive changes during epidemics in São Jose de Rio Preto, Brazil, 2006–2007

Global dengue virus spread in tropical and sub-tropical regions has become a major international public health concern. It is evident that DENV genetic diversity plays a significant role in the immunopathology of the disease and that the identification of polymorphisms associated with adaptive responses is important for vaccine development. The investigation of naturally occurring genomic variants may play an important role in the comprehension of different adaptive strategies used by these mutants to evade the human immune system. In order to elucidate this role we sequenced the complete polyprotein-coding region of thirty-three DENV-3 isolates to characterize variants circulating under high endemicity in the city of São José de Rio Preto, Brazil, during the onset of the 2006-07 epidemic. By inferring the evolutionary history on a local-scale and estimating rates of synonymous (dS) and nonsynonimous (dN) substitutions, we have documented at least two different introductions of DENV-3 into the city and detected 10 polymorphic codon sites under significant positive selection (dN/dS > 1) and 8 under significant purifying selection (dN/dS < 1). We found several polymorphic amino acid coding sites in the envelope (15), NS1 (17), NS2A (11), and NS5 (24) genes, which suggests that these genes may be experiencing relatively recent adaptive changes. Furthermore, some polymorphisms correlated with changes in the immunogenicity of several epitopes. Our study highlights the existence of significant and informative DENV variability at the spatio-temporal scale of an urban outbreak.

El uso de la información genética en la gestión de la biodiversidad marina de islas oceánicas : el cangrejo moro "Grapsus adscensionis" en la Macaronesia

Máster Oficial en Gestión Costera

Image analysis in the morphological and morphometric study of teeth

The subject of this doctoral dissertation concerns the definition of a new methodology for the morphological and morphometric study of fossilized human teeth, and therefore strives to provide a contribution to the reconstruction of human evolutionary history that proposes to extend to the different species of hominid fossils. Standardized investigative methodologies are lacking both regarding the orientation of teeth subject to study and in the analysis that can be carried out on these teeth once they are oriented. The opportunity to standardize a primary analysis methodology is furnished by the study of certain early Neanderthal and preneanderthal molars recovered in two caves in southern Italy [Grotta Taddeo (Taddeo Cave) and Grotta del Poggio (Poggio Cave), near Marina di Camerata, Campania]. To these we can add other molars of Neanderthal and modern man of the upper Paleolithic era, specifically scanned in the paleoanthropology laboratory of the University of Arkansas (Fayetteville, Arkansas, USA), in order to increase the paleoanthropological sample data and thereby make the final results of the analyses more significant. The new analysis methodology is rendered as follows: 1. Standardization of an orientation system for primary molars (superior and inferior), starting from a scan of a sample of 30 molars belonging to modern man (15 M1 inferior and 15 M1 superior), the definition of landmarks, the comparison of various systems and the choice of a system of orientation for each of the two dental typologies. 2. The definition of an analysis procedure that considers only the first 4 millimeters of the dental crown starting from the collar: 5 sections parallel to the plane according to which the tooth has been oriented are carried out, spaced 1 millimeter between them. The intention is to determine a method that allows for the differentiation of fossilized species even in the presence of worn teeth. 3. Results and Conclusions. The new approach to the study of teeth provides a considerable quantity of information that can better be evaluated by increasing the fossil sample data. It has been demonstrated to be a valid tool in evolutionary classification that has allowed (us) to differentiate the Neanderthal sample from that of modern man. In a particular sense the molars of Grotta Taddeo, which up until this point it has not been possible to determine with exactness their species of origin, through the present research they are classified as Neanderthal.

Metodi molecolari nello studio della variabilità genetica di Tuber borchii

Tuber borchii (Ascomycota, order Pezizales) is highly valued truffle sold in local markets in Italy. Despite its economic importance, knowledge on its distribution and population variation is scarce. The objective of this work was to investigate the evolutionary forces shaping the genetic structure of this fungus using coalescent and phylogenetic methods to reconstruct the evolutionary history of populations in Italy. To assess population structure, 61 specimens were collected from 11 different Provinces of Italy. Sampling was stratified across hosts and habitats to maximize coverage in native oak and pine stands and both mychorrizae and fruiting bodies were collected. Samples were identified considering anatomo-morphological characters. DNA was extracted and both multilocus (AFLP) and single-locus (18 loci from rDNA, nDNA, and mtDNA) approaches were used to look for polymorphisms. Screening AFLP profiles, both Jaccard and Dice coefficients of similarity were utilized to transform binary matrix into a distance matrix and then to desume Neighbour-Joining trees. Though these are only preliminary examinations, phylogenetic trees were totally concordant with those deriving from single locus analyses. Phylogenetic analyses of the nuclear loci were performed using maximum likelihood with PAUP and a combined phylogenetic inference, using Bayesian estimation with all nuclear gene regions, was carried out. To reconstruct the evolutionary history, we estimated recurrent migration, migration across the history of the sample, and estimated the mutation and approximate age of mutations in each tree using SNAP Workbench. The combined phylogenetic tree using Bayesian estimation suggests that there are two main haplotypes that are difficult to be differentiated on the basis of morphology, of ecological parameters and symbiontic tree. Between these two lineages, that occur in sympatry within T. borchii populations, there is no evidence of recurrent migration. However, migration over the history of the sample was asymmetrical suggesting that isolation was a result of interrupted gene flow followed by range expansion. Low levels of divergence between the haplotypes indicate that there are likely to be two cryptic species within the T. borchii population sampled. Our results suggest that isolation between populations of T. borchii could have led to reproductive isolation between two lineages. This isolation is likely due to sympatric speciation caused by a multiple colonization from different refugia or a recent isolation. In attempting to determinate whether these haplotypes represent separate species or a partition of the same species we applied Biological and Mechanistic species Concepts. Notwithstanding, further analyses are necessary to evaluate if selection favoured premating or post-mating isolation.

Computational methods for the analysis of protein structure and function

The vast majority of known proteins have not yet been experimentally characterized and little is known about their function. The design and implementation of computational tools can provide insight into the function of proteins based on their sequence, their structure, their evolutionary history and their association with other proteins. Knowledge of the three-dimensional (3D) structure of a protein can lead to a deep understanding of its mode of action and interaction, but currently the structures of <1% of sequences have been experimentally solved. For this reason, it became urgent to develop new methods that are able to computationally extract relevant information from protein sequence and structure. The starting point of my work has been the study of the properties of contacts between protein residues, since they constrain protein folding and characterize different protein structures. Prediction of residue contacts in proteins is an interesting problem whose solution may be useful in protein folding recognition and de novo design. The prediction of these contacts requires the study of the protein inter-residue distances related to the specific type of amino acid pair that are encoded in the so-called contact map. An interesting new way of analyzing those structures came out when network studies were introduced, with pivotal papers demonstrating that protein contact networks also exhibit small-world behavior. In order to highlight constraints for the prediction of protein contact maps and for applications in the field of protein structure prediction and/or reconstruction from experimentally determined contact maps, I studied to which extent the characteristic path length and clustering coefficient of the protein contacts network are values that reveal characteristic features of protein contact maps. Provided that residue contacts are known for a protein sequence, the major features of its 3D structure could be deduced by combining this knowledge with correctly predicted motifs of secondary structure. In the second part of my work I focused on a particular protein structural motif, the coiled-coil, known to mediate a variety of fundamental biological interactions. Coiled-coils are found in a variety of structural forms and in a wide range of proteins including, for example, small units such as leucine zippers that drive the dimerization of many transcription factors or more complex structures such as the family of viral proteins responsible for virus-host membrane fusion. The coiled-coil structural motif is estimated to account for 5-10% of the protein sequences in the various genomes. Given their biological importance, in my work I introduced a Hidden Markov Model (HMM) that exploits the evolutionary information derived from multiple sequence alignments, to predict coiled-coil regions and to discriminate coiled-coil sequences. The results indicate that the new HMM outperforms all the existing programs and can be adopted for the coiled-coil prediction and for large-scale genome annotation. Genome annotation is a key issue in modern computational biology, being the starting point towards the understanding of the complex processes involved in biological networks. The rapid growth in the number of protein sequences and structures available poses new fundamental problems that still deserve an interpretation. Nevertheless, these data are at the basis of the design of new strategies for tackling problems such as the prediction of protein structure and function. Experimental determination of the functions of all these proteins would be a hugely time-consuming and costly task and, in most instances, has not been carried out. As an example, currently, approximately only 20% of annotated proteins in the Homo sapiens genome have been experimentally characterized. A commonly adopted procedure for annotating protein sequences relies on the "inheritance through homology" based on the notion that similar sequences share similar functions and structures. This procedure consists in the assignment of sequences to a specific group of functionally related sequences which had been grouped through clustering techniques. The clustering procedure is based on suitable similarity rules, since predicting protein structure and function from sequence largely depends on the value of sequence identity. However, additional levels of complexity are due to multi-domain proteins, to proteins that share common domains but that do not necessarily share the same function, to the finding that different combinations of shared domains can lead to different biological roles. In the last part of this study I developed and validate a system that contributes to sequence annotation by taking advantage of a validated transfer through inheritance procedure of the molecular functions and of the structural templates. After a cross-genome comparison with the BLAST program, clusters were built on the basis of two stringent constraints on sequence identity and coverage of the alignment. The adopted measure explicity answers to the problem of multi-domain proteins annotation and allows a fine grain division of the whole set of proteomes used, that ensures cluster homogeneity in terms of sequence length. A high level of coverage of structure templates on the length of protein sequences within clusters ensures that multi-domain proteins when present can be templates for sequences of similar length. This annotation procedure includes the possibility of reliably transferring statistically validated functions and structures to sequences considering information available in the present data bases of molecular functions and structures.

The pre-alpine evolution of the basement of the Pelagonian Zone and the Vardar Zone, Greece

The Pelagonian Zone and the Vardar Zone in Greece represent the western part of the Hellenide hinterland (Internal Hellenides). While the Pelagonian Zone comprises predominantly crystalline basement and sedimentary cover rocks, the Vardar Zone has long been regarded as an ophiolite-decorated suture zone separating the Pelagonian Zone from the Serbo-Macedonian Massif to the east. Felsic basement rocks from both areas, with the main focus put on the Pelagonian Zone, were dated in order to identify the major crust-forming episodes and to improve the understanding of the evolutionary history of the region. The interpretation of the single-zircon geochronology results was aided by geochemical investigations. The majority of the basement rocks from the Pelagonian Zone yielded Permo-Carboniferous intrusion ages around 300 Ma, underlining the importance of this crust-forming event for the Internal Hellenides of Greece. Geochemically these basement rocks are classified as subduction-related granitoids, which formed in an active continental margin setting. An important result was the identification of a Precambrian crustal unit within the crystalline basement of the Pelagonian Zone. Orthogneisses from the NW Pelagonian Zone yielded Neoproterozoic ages of c. 700 Ma and are so far the oldest known rocks in Greece. These basement rocks, which are also similar to active margin granitoids, were interpreted as remnants of a terrane, the Florina Terrane, which can be correlated to a Pan-African or Cadomian arc. Since the gneisses contain inherited zircons of Middle to Late Proterozoic ages, the original location of the Florina Terrane was probably at the northwestern margin of Gondwana. In the Vardar Zone an important phase of Upper Jurassic felsic magmatism is documented by igneous formation ages ranging from 155 to 164 Ma. The chemical and isotopic composition of these rocks is also in accord with their formation in a volcanic-arc setting at an active continental margin. Older continental material incorporated in the Vardar Zone is documented by 319-Ma-old gneisses and by inherited zircons of mainly Middle Palaeozoic ages. The prevalence of subduction-related igneous rocks indicates that arc formation and accretion orogeny were the most important processes during the evolution of this part of the Internal Hellenides. The geochronological results demonstrate that most of the Pelagonian Zone and the Vardar Zone crystalline basement formed during distinct pre-Alpine episodes at c. 700, 300 and 160 Ma with a predominance of the Permo-Carboniferous magmatic phase.

Molecular phylogeography of the Woodland Ringlet (Erebia medusa [Denis and Schiffermüller] 1775) in Europe

Phylogeography is a recent field of biological research that links phylogenetics to biogeography through deciphering the imprint that evolutionary history has left on the genetic structure of extant populations. During the cold phases of the successive ice ages, which drastically shaped species’ distributions since the Pliocene, populations of numerous species were isolated in refugia where many of them evolved into different genetic lineages. My dissertation deals with the phylogeography of the Woodland Ringlet (Erebia medusa [Denis and Schiffermüller] 1775) in Central and Eastern Europe. This Palaearctic butterfly species is currently distributed from central France and south eastern Belgium over large parts of Central Europe and southern Siberia to the Pacific. It is absent from those parts of Europe with mediterranean, oceanic and boreal climates. It was supposed to be a Siberian faunal element with a rather homogeneous population structure in Central Europe due to its postglacial expansion out of a single eastern refugium. An already existing evolutionary scenario for the Woodland Ringlet in Central and Eastern Europe is based on nuclear data (allozymes). To know if this is corroborated by organelle evolutionary history, I sequenced two mitochondrial markers (part of the cytochrome oxydase subunit one and the control region) for populations sampled over the same area. Phylogeography largely relies on the construction of networks of uniparentally inherited haplotypes that are compared to geographic haplotype distribution thanks to recent developed methods such as nested clade phylogeographic analysis (NCPA). Several ring-shaped ambiguities (loops) emerged from both haplotype networks in E. medusa. They can be attributed to recombination and homoplasy. Such loops usually avert the straightforward extraction of the phylogeographic signal contained in a gene tree. I developed several new approaches to extract phylogeographic information in the presence of loops, considering either homoplasy or recombination. This allowed me to deduce a consistent evolutionary history for the species from the mitochondrial data and also adds plausibility for the occurrence of recombination in E. medusa mitochondria. Despite the fact that the control region is assumed to have a lack of resolving power in other species, I found a considerable genetic variation of this marker in E. medusa which makes it a useful tool for phylogeographic studies. In combination with the allozyme data, the mitochondrial genome supports the following phylogeographic scenario for E. medusa in Europe: (i) a first vicariance, due to the onset of the Würm glaciation, led to the formation of several major lineages, and is mirrored in the NCPA by restricted gene flow, (ii) later on further vicariances led to the formation of two sub-lineages in the Western lineage and two sub-lineages in the Eastern lineage during the Last Glacial Maximum or Older Dryas; additionally the NCPA supports a restriction of gene flow with isolation by distance, (iii) finally, vicariance resulted in two secondary sub-lineages in the area of Germany and, maybe, to two other secondary sub-lineages in the Czech Republic. The last postglacial warming was accompanied by strong range expansions in most of the genetic lineages. The scenario expected for a presumably Siberian faunal element such as E. medusa is a continuous loss of genetic diversity during postglacial westward expansion. Hence, the pattern found in this thesis contradicts a typical Siberian origin of E. medusa. In contrast, it corroboratess the importance of multiple extra-Mediterranean refugia for European fauna as it was recently assumed for other continental species.

Die Überarbeitung der Fernsehrichtlinie unter Berücksichtigung der EuGH-Entscheidungen zur bisherigen Fassung

Gegenstand der vorliegenden Arbeit ist die Überarbeitung der Richtlinie 89/552/EWG des Rates zur Koordinierung bestimmter Rechts- und Verwaltungsvorschriften der Mitgliedstaaten über die Ausübung der Fernsehtätigkeit, welche aus praktikablen Gründen meist als „(EG-)Fernsehrichtlinie“ bezeichnet wird. Sie bildet den Eckpfeiler der audiovisuellen Politik der EU. Seit Erlass der Fernsehrichtlinie im Jahre 1989 bewirkt der technologische Fortschritt jedoch zunehmend enorme Veränderungen nicht nur im Bereich des klassischen Fernsehens, sondern auch und vor allem im Bereich der neuen Medien. Ausgangspunkt hierfür ist die Verbesserung der Digitaltechnologie, die ihrerseits wiederum technische Konvergenzprozesse begünstigt. Diese Entwicklungen führen nicht nur zu einer Vervielfachung von Übertragungskapazitäten und –techniken, sondern ermöglichen neben neuen Formen audiovisueller Angebote auch die Entstehung neuer Dienste. Unsere Medienlandschaft steht vor „epochalen Umbrüchen“. Im Hinblick auf diese Vorgänge wird seit geraumer Zeit eine Überarbeitung der EG-Fernsehrichtlinie angestrebt, um dem technologischen Fortschritt auch „regulatorisch“ gerecht werden zu können. Diesem Überarbeitungsprozess möchte sich die vorliegende Arbeit widmen, indem sie die Fernsehrichtlinie in einem ersten Teil sowohl inhaltlich wie auch hinsichtlich ihrer Entstehungsgeschichte und der zu ihr ergangenen EuGH-Entscheidungen erläutert. Anschließend werden alle Überarbeitungsvorgänge der Fernsehrichtlinie seit 1997 dargestellt, um sodann die aktuellen Reformansätze analysieren und bewerten zu können. Aus zeitlichen Gründen (der neue Richtlinienvorschlag der Kommission vom 13. Dezember 2005 wurde ca. 2 Wochen vor dem Abgabetermin der Arbeit verabschiedet) sind die Ausführungen zum Entwurf der neuen „Richtlinie über audiovisuelle Mediendienste“ allerdings relativ knapp gehalten.

Environmental influence on the phenotype: Morphological variation in human dentition

This work is about the role that environment plays in the production of evolutionary significant variations. It starts with an historical introduction about the concept of variation and the role of environment in its production. Then, I show how a lack of attention to these topics may lead to serious mistakes in data interpretation. A statistical re-analysis of published data on the effects of malnutrition on dental eruption, shows that what has been interpreted as an increase in the mean value, is actually linked to increase of variability. In Chapter 3 I present the topic of development as a link between variability and environmental influence, giving a review of the possible mechanisms by which development influences evolutionary dynamics. Chapter 4 is the core chapter of the thesis; I investigated the role of environment in the development of dental morphology. I used dental hypoplasia as a marker of stress, characterizing two groups. Comparing the morphology of upper molars in the two groups, three major results came out: (i) there is a significant effect of environmental stressors on the overall morphology of upper molars; (ii) the developmental response increases morphological variability of the stressed population; (iii) increase of variability is directional: stressed individuals have increased cusps dimensions and number. I also hypothesized the molecular mechanisms that could be responsible of the observed effects. In Chapter 5, I present future perspectives for developing this research. The direction of dental development response is the same direction of the trend in mammalian dental evolution. Since malnutrition triggers the developmental response, and this particular kind of stressor must have been very common in our class evolutionary history, I propose the possibility that environmental stress actively influenced mammals evolution. Moreover, I discuss the possibility of reconsidering the role of natural selection in the evolution of dental morphology.

Species identification of archived fish bones collected from the Mediterranean Sea 100 years ago using molecular techniques

With the discovery that DNA can be successfully recovered from museum collections, a new source of genetic information has been provided to extend our comprehension of the evolutionary history of species. However, historical specimens are often mislabeled or report incorrect information of origin, thus accurate identification of specimens is essential. Due to the highly damaged nature of ancient DNA many pitfalls exist and particular precautions need to be considered in order to perform genetic analysis. In this study we analyze 208 historical remains of pelagic fishes collected in the beginning of the 20th century. Through the adaptation of existing protocols, usually applied to human remains, we manage to successfully retrieve valuable genetic material from almost all of the examined samples using a guanidine and silica column-based approach. The combined use of two mitochondrial markers cytochrome-oxidase-1(mtDNA COI) and Control Region (mtDNA CR), and the nuclear marker first internal transcriber space (ITS1) allowed us to identify the majority of the examined specimens using traditional PCR and Sanger sequencing techniques. The creation of primers capable of amplifying heavily degraded DNA have great potential for future uses, both in ancient and in modern investigation. The methodologies developed in this study can in fact be applied for other ancient fish specimens as well as cooked or canned samples.

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.