Convergent evolution of [D-Leucine1] microcystin-LR in taxonomically disparate cyanobacteria

Abstract Background Many important toxins and antibiotics are produced by non-ribosomal biosynthetic pathways. Microcystins are a chemically diverse family of potent peptide toxins and the end-products of a hybrid NRPS and PKS secondary metabolic pathway. They are produced by a variety of cyanobacteria and are responsible for the poisoning of humans as well as the deaths of wild and domestic animals around the world. The chemical diversity of the microcystin family is attributed to a number of genetic events that have resulted in the diversification of the pathway for microcystin assembly. Results Here, we show that independent evolutionary events affecting the substrate specificity of the microcystin biosynthetic pathway have resulted in convergence on a rare [D-Leu1] microcystin-LR chemical variant. We detected this rare microcystin variant from strains of the distantly related genera Microcystis, Nostoc, and Phormidium. Phylogenetic analysis performed using sequences of the catalytic domains within the mcy gene cluster demonstrated a clear recombination pattern in the adenylation domain phylogenetic tree. We found evidence for conversion of the gene encoding the McyA2 adenylation domain in strains of the genera Nostoc and Phormidium. However, point mutations affecting the substrate-binding sequence motifs of the McyA2 adenylation domain were associated with the change in substrate specificity in two strains of Microcystis. In addition to the main [D-Leu1] microcystin-LR variant, these two strains produced a new microcystin that was identified as [Met1] microcystin-LR. Conclusions Phylogenetic analysis demonstrated that both point mutations and gene conversion result in functional mcy gene clusters that produce the same rare [D-Leu1] variant of microcystin in strains of the genera Microcystis, Nostoc, and Phormidium. Engineering pathways to produce recombinant non-ribosomal peptides could provide new natural products or increase the activity of known compounds. Our results suggest that the replacement of entire adenylation domains could be a more successful strategy to obtain higher specificity in the modification of the non-ribosomal peptides than point mutations.

Distribution of nitrogen fixation and nitrogenase-like sequences amongst microbial genomes

Abstract Background The metabolic capacity for nitrogen fixation is known to be present in several prokaryotic species scattered across taxonomic groups. Experimental detection of nitrogen fixation in microbes requires species-specific conditions, making it difficult to obtain a comprehensive census of this trait. The recent and rapid increase in the availability of microbial genome sequences affords novel opportunities to re-examine the occurrence and distribution of nitrogen fixation genes. The current practice for computational prediction of nitrogen fixation is to use the presence of the nifH and/or nifD genes. Results Based on a careful comparison of the repertoire of nitrogen fixation genes in known diazotroph species we propose a new criterion for computational prediction of nitrogen fixation: the presence of a minimum set of six genes coding for structural and biosynthetic components, namely NifHDK and NifENB. Using this criterion, we conducted a comprehensive search in fully sequenced genomes and identified 149 diazotrophic species, including 82 known diazotrophs and 67 species not known to fix nitrogen. The taxonomic distribution of nitrogen fixation in Archaea was limited to the Euryarchaeota phylum; within the Bacteria domain we predict that nitrogen fixation occurs in 13 different phyla. Of these, seven phyla had not hitherto been known to contain species capable of nitrogen fixation. Our analyses also identified protein sequences that are similar to nitrogenase in organisms that do not meet the minimum-gene-set criteria. The existence of nitrogenase-like proteins lacking conserved co-factor ligands in both diazotrophs and non-diazotrophs suggests their potential for performing other, as yet unidentified, metabolic functions. Conclusions Our predictions expand the known phylogenetic diversity of nitrogen fixation, and suggest that this trait may be much more common in nature than it is currently thought. The diverse phylogenetic distribution of nitrogenase-like proteins indicates potential new roles for anciently duplicated and divergent members of this group of enzymes.

Cryptic species in the cosmopolitan Bugulan eritina complex (Bryozoa,Cheilostomata)

Previous analyses of the mitochondrial gene cytochrome c oxidase subunit 1 (COI) and γ-proteobacterial endosymbiont diversity have suggested that the marine bryozoan Bugula neritina is a complex of three cryptic species, namely Types S, D and N. Types D and N were previously reported to have restricted distributions along California (western USA) and Delaware and Connecticut (eastern USA), respectively, whereas Type S is considered widespread in tropical, subtropical and temperate regions due to anthropogenic transport. Here, Bayesian species delimitation analysis of a data set composed of two mitochondrial (COI and large ribosomal RNA subunit [16S]) and two nuclear genes (dynein light chain roadblock type-2 protein [DYN] and voltage-dependent anion-selective channel protein [VDAC]) demonstrated that Types S, D and N correspond to three biological species. This finding was significantly supported, in spite of the combinations of priors applied for ancestral population size and root age. Furthermore, COI sequences were used to assess the introduction patterns of the cosmopolitan Type S species. Two COI haplotypes of Type S (S1a and S1d) were found occurring at a global scale. Mantel tests showed correlation between these haplotypes and local sea surface temperature tolerance. Accordingly, the distributions of Type S haplotypes may reflect intraspecific temperature tolerance variation, in addition to the role of introduction vectors. Finally, we show that the Type N may also have been introduced widely, as this species was found for the first time in Central California and north-eastern Australia.

Variational second order flow estimation for PIV sequences

[EN] We present in this paper a variational approach to accurately estimate simultaneously the velocity field and its derivatives directly from PIV image sequences. Our method differs from other techniques that have been presented in the literature in the fact that the energy minimization used to estimate the particles motion depends on a second order Taylor development of the flow. In this way, we are not only able to compute the motion vector field, but we also obtain an accurate estimation of their derivatives. Hence, we avoid the use of numerical schemes to compute the derivatives from the estimated flow that usually yield to numerical amplification of the inherent uncertainty on the estimated flow. The performance of our approach is illustrated with the estimation of the motion vector field and the vorticity on both synthetic and real PIV datasets.

An evolutive approach for smile recognition in video sequences

Facial expression recognition is one of the most challenging research areas in the image recognition ¯eld and has been actively studied since the 70's. For instance, smile recognition has been studied due to the fact that it is considered an important facial expression in human communication, it is therefore likely useful for human–machine interaction. Moreover, if a smile can be detected and also its intensity estimated, it will raise the possibility of new applications in the future

Species identification, distributioin and reproductive interactions of common smoot-hound and blackspotted smooth-hound shark (Genus mustelus) in the Adriatic Sea

Longstanding taxonomic ambiguity and uncertainty exist in the identification of the common (M. mustelus) and blackspotted (M. punctulatus) smooth-hound in the Adriatic Sea. The lack of a clear and accurate method of morphological identification, leading to frequent misidentification, prevents the collation of species-specific landings and survey data for these fishes and hampers the delineation of the distribution ranges and stock boundaries of the species. In this context, adequate species-specific conservation and management strategies can not be applied without risks of population declining and local extinction. In this thesis work I investigated the molecular ecology of the two smooth-hound sharks which are abundant in the demersal trawl surveys carried out in the NC Adriatic Sea to monitor and assess the fishery resources. Ecological and evolutionary relationships were assessed by two molecular tests: a DNA barcoding analysis to improve species identification (and consequently the knowledge of their spatial ecology and taxonomy) and a hybridization assay based on the nuclear codominant marker ITS2 to evaluate reproductive interactions (hybridization or gene introgression). The smooth-hound sharks (N=208) were collected during the MEDITS 2008 and 2010 campaigns along the Italian and Croatian coasts of the Adriatic Sea, in the Sicilian Channel and in the Algerian fisheries. Since the identification based on morphological characters is not strongly reliable, I performed a molecular identification of the specimens producing for each one the cytochrome oxidase subunit 1 (COI) gene sequence (ca. 640 bp long) and compared them with reference sequences from different databases (GenBank and BOLD). From these molecular ID data I inferred the distribution of the two target species in the NC Adriatic Sea. In almost the totality of the MEDITS hauls I found no evidence of species sympatry. The data collected during the MEDITS survey showed an almost different distribution of M. mustelus (confined along the Italian coasts) and M. punctulatus (confined along the Croatian coasts); just one sample (Gulf of Venice, where probably the ranges of the species overlap) was found to have catches of both the species. Despite these data results suggested no interaction occurred between my two target species at least during the summertime (the period in which MEDITS survey is carried out), I still wanted to know if there were inter-species reproductive interactions so I developed a simple molecular genetic method to detect hybridization. This method is based on DNA sequence polymorphism among species in the nuclear ribosomal Internal Transcribed Spacer 2 locus (ITS2). Its application to the 208 specimens collected raised important questions regarding the ecology of this two species in the Adriatic Sea. In fact results showed signs of hybridization and/or gene introgression in two sharks collected during the trawl survey of 2008 and one collected during the 2010 one along the Italian and Croatian coasts. In the case that it will be confirmed the hybrid nature of these individuals, a spatiotemporal overlapping of the mating behaviour and ecology must occur. At the spatial level, the northern part of the Adriatic Sea (an area where the two species occur with high frequency of immature individuals) could likely play the role of a common nursery area for both species.

Analysis of eruptive and seismic sequences to improve the short-and long-term eruption forecasting

Forecasting the time, location, nature, and scale of volcanic eruptions is one of the most urgent aspects of modern applied volcanology. The reliability of probabilistic forecasting procedures is strongly related to the reliability of the input information provided, implying objective criteria for interpreting the historical and monitoring data. For this reason both, detailed analysis of past data and more basic research into the processes of volcanism, are fundamental tasks of a continuous information-gain process; in this way the precursor events of eruptions can be better interpreted in terms of their physical meanings with correlated uncertainties. This should lead to better predictions of the nature of eruptive events. In this work we have studied different problems associated with the long- and short-term eruption forecasting assessment. First, we discuss different approaches for the analysis of the eruptive history of a volcano, most of them generally applied for long-term eruption forecasting purposes; furthermore, we present a model based on the characteristics of a Brownian passage-time process to describe recurrent eruptive activity, and apply it for long-term, time-dependent, eruption forecasting (Chapter 1). Conversely, in an effort to define further monitoring parameters as input data for short-term eruption forecasting in probabilistic models (as for example, the Bayesian Event Tree for eruption forecasting -BET_EF-), we analyze some characteristics of typical seismic activity recorded in active volcanoes; in particular, we use some methodologies that may be applied to analyze long-period (LP) events (Chapter 2) and volcano-tectonic (VT) seismic swarms (Chapter 3); our analysis in general are oriented toward the tracking of phenomena that can provide information about magmatic processes. Finally, we discuss some possible ways to integrate the results presented in Chapters 1 (for long-term EF), 2 and 3 (for short-term EF) in the BET_EF model (Chapter 4).

Transposable elements: structure and dynamic of the autonomous retrotransposon R2 in Arthropoda with non-canonical reproduction

Eukaryotic ribosomal DNA constitutes a multi gene family organized in a cluster called nucleolar organizer region (NOR); this region is composed usually by hundreds to thousands of tandemly repeated units. Ribosomal genes, being repeated sequences, evolve following the typical pattern of concerted evolution. The autonomous retroelement R2 inserts in the ribosomal gene 28S, leading to defective 28S rDNA genes. R2 element, being a retrotransposon, performs its activity in the genome multiplying its copy number through a “copy and paste” mechanism called target primed reverse transcription. It consists in the retrotranscription of the element’s mRNA into DNA, then the DNA is integrated in the target site. Since the retrotranscription can be interrupted, but the integration will be carried out anyway, truncated copies of the element will also be present in the genome. The study of these truncated variants is a tool to examine the activity of the element. R2 phylogeny appears, in general, not consistent with that of its hosts, except some cases (e.g. Drosophila spp. and Reticulitermes spp.); moreover R2 is absent in some species (Fugu rubripes, human, mouse, etc.), while other species have more R2 lineages in their genome (the turtle Mauremys reevesii, the Japanese beetle Popilia japonica, etc). R2 elements here presented are isolated in 4 species of notostracan branchiopods and in two species of stick insects, whose reproductive strategies range from strict gonochorism to unisexuality. From sequencing data emerges that in Triops cancriformis (Spanish gonochoric population), in Lepidurus arcticus (two putatively unisexual populations from Iceland) and in Bacillus rossius (gonochoric population from Capalbio) the R2 elements are complete and encode functional proteins, reflecting the general features of this family of transposable elements. On the other hand, R2 from Italian and Austrian populations of T. cancriformis (respectively unisexual and hermaphroditic), Lepidurus lubbocki (two elements within the same Italian population, gonochoric but with unfunctional males) and Bacillus grandii grandii (gonochoric population from Ponte Manghisi) have sequences that encode incomplete or non-functional proteins in which it is possible to recognize only part of the characteristic domains. In Lepidurus couesii (Italian gonochoric populations) different elements were found as in L. lubbocki, and the sequencing is still in progress. Two hypothesis are given to explain the inconsistency of R2/host phylogeny: vertical inheritance of the element followed by extinction/diversification or horizontal transmission. My data support previous study that state the vertical transmission as the most likely explanation; nevertheless horizontal transfer events can’t be excluded. I also studied the element’s activity in Spanish populations of T. cancriformis, in L. lubbocki, in L. arcticus and in gonochoric and parthenogenetic populations of B. rossius. In gonochoric populations of T. cancriformis and B. rossius I found that each individual has its own private set of truncated variants. The situation is the opposite for the remaining hermaphroditic/parthenogenetic species and populations, all individuals sharing – in the so far analyzed samples - the majority of variants. This situation is very interesting, because it isn’t concordant with the Muller’s ratchet theory that hypothesizes the parthenogenetic populations being either devoided of transposable elements or TEs overloaded. My data suggest a possible epigenetic mechanism that can block the retrotransposon activity, and in this way deleterious mutations don’t accumulate.

A molecular phylogeny of bivalve mollusks: ancient radiations and divergences as revealed by mitochondrial genes

The main scope of my PhD is the reconstruction of the large-scale bivalve phylogeny on the basis of four mitochondrial genes, with samples taken from all major groups of the class. To my knowledge, it is the first attempt of such a breadth in Bivalvia. I decided to focus on both ribosomal and protein coding DNA sequences (two ribosomal encoding genes -12s and 16s -, and two protein coding ones - cytochrome c oxidase I and cytochrome b), since either bibliography and my preliminary results confirmed the importance of combined gene signals in improving evolutionary pathways of the group. Moreover, I wanted to propose a methodological pipeline that proved to be useful to obtain robust results in bivalves phylogeny. Actually, best-performing taxon sampling and alignment strategies were tested, and several data partitioning and molecular evolution models were analyzed, thus demonstrating the importance of molding and implementing non-trivial evolutionary models. In the line of a more rigorous approach to data analysis, I also proposed a new method to assess taxon sampling, by developing Clarke and Warwick statistics: taxon sampling is a major concern in phylogenetic studies, and incomplete, biased, or improper taxon assemblies can lead to misleading results in reconstructing evolutionary trees. Theoretical methods are already available to optimize taxon choice in phylogenetic analyses, but most involve some knowledge about genetic relationships of the group of interest, or even a well-established phylogeny itself; these data are not always available in general phylogenetic applications. The method I proposed measures the "phylogenetic representativeness" of a given sample or set of samples and it is based entirely on the pre-existing available taxonomy of the ingroup, which is commonly known to investigators. Moreover, it also accounts for instability and discordance in taxonomies. A Python-based script suite, called PhyRe, has been developed to implement all analyses.

Systematics, phylogeny and biogeography of Cousinia (Asteraceae)

Die Systematik, Phylogenie und Biogeographie der Gattung Cousinia (Asteraceae, Cardueae) als größter Gattung der Tribus Cardueae mit mehr als 600 Arten wurde untersucht. Diese Dissertation umfasst drei Hauptteile: Im ersten Teil wurde die Phylogenie und Evolution des Arctium-Cousinia-Komplexes Untersucht. Dieser Gattungskomplex enthält Arctium, Cousinia, Hypacanthium und Schmalhausenia und zeigt die höchste Diversität in der Irano-Turanischen Region und in den Gebirgen Zentralasiens. Es wurden ITS und rpS4-trnT-trnL-Sequenzen für insgesamt 138 Arten generiert, darunter von 129 (von ca. 600) Arten von Cousinia. Wie in früheren Analysen bereits gefunden, ist Cousinia nicht monophyletisch. Stattdessen sind Cousinia subg. Cynaroides und subg. Hypacanthodes mit insgesamt ca. 30 Arten enger mit Arctium, Hypacanthium und Schmalhausenia (Arctioid Clade) als mit subg. Cousinia (Cousinioid Clade) verwandt. Die Arctioid und Cousiniod clades werden auch durch Pollenmorphologie und Chromosomenzahl unterstützt, wie bereits früher bekannt war. In dem Arctioid Clade entsprechen morphologische Gattungsgrenzen, basierend auf Blattform, Blattbedornung und Morphologie der Involukralblätter, nicht den in der molekularen analyse gefundenen clades. Es kann keine taxonomische Lösung für diesen Konflikt gefunden werden, und die gennanten Merkmale wurden als homoplastisch betrachtet. Obwohl die phylogenetische Auflösung in dem Cousinioid Clade schlecht ist, enthalten die ITS und rpS4-trnT-trnL-Sequenzen phylogenetische Information. So gruppierten z.B. die sechs annuellen Arten in zwei Gruppen. Schlechte phylogenetische Auflösung resultiert wahrscheinlich aus dem Mangel an Merkmalen und der großen Artenzahl in dieser artenreichen und vergleichsweise jungen (ca. 8,7 mya) Linie. Artbildung in dem Cousinioid Clade scheint hauptsächlich allopatrisch zu sein. Der zweite Teil der Dissertation untersucht die Rolle der Hybridisierung in der Evolution von Cousinia s.s. Die in der Vergangenheit publizierteten 28 Hybrid-Kombinationen und 11 Zwischenformen wurden kritisch geprüft, und zwei Hybridindividuen wurden morphologisch und molekular untersucht. Die vermutlichen oder nachgewiesenen Eltern der Hybriden und Zwischenformen wurden auf die aus einer Bayesischen Analyse der ITS-Sequenzen von 216 Arten von Cousinia und verwandten Gattungen resultierenden Phylogenie aufgetragen. Weder Hybriden zwischen dem Cousinioid Clade und anderen Haupt-Claden des Arctium-Cousinia-Komplexes noch zwischen annuellen und perennirenden Arten von Cousinia s.s. wurden beobachtet. Die Ergebnisse zeigen eindeutig, dass Hybridisierung in Cousinia möglich is, und dass ca. 10,7% der Arten an interspezifischer Hybridisierung beteiligt sind. Obwohl Hybridisierung in Cousinia s.s. stattfindet und zu den Schwierigkeiten bei der Rekonstruktion ihrer phylogenetischen Geschichte beitragen könnte, war ihre Rolle für die Entwicklung und Diversität der Gruppe offenbar gering. Im dritten Teil wird eine taxonomische Revision der C. sect. Cynaroideae präsentiert. Cousinia sect. Cynaroideae, die größte Sektion der Gattung mit 110 veröffentlichten Arten, zeichnet sich durch eine Chromosomenzahl von 2n = 24 und durch ± herablaufende Blätter und Hüllblätter mit Anhängseln aus. Sie kommt im Iran, Irak, dem Kaukasus, der Türkei, Turkmenistan, Afghanistan, Pakistan, dem Libanon und Anti-Libanon vor und hat ihre Hauptzentren der Artdiversität im westlichen und nordwestlichen Iran, im Irak und in der südöstlichen Türkei. Die Revision dieser Gruppe, hauptsächlich basierend auf der Untersuchung von ca. 2250 Herbarbögen, führte zu einer Verringerung der Artenzahl auf 31 Arten mit acht Unterarten. Alle Arten werden typifiziert und ausgeschlüsselt, und Beschreibungen, Abbildungen und Verbreitungskarten werden für jede Art angegeben.

Phylogenomische Analysen bei Metazoen - Zur Stellung der Xenoturbellida und Syndermata

Im Rahmen der vorliegenden Dissertation wurde die phylogenetischen Stellungen der Xenoturbellida (Deuterostomia) und der Syndermata (Protostomia) mit phylogenomischen Techniken untersucht. Auf methodischer Ebene konnte gezeigt werden, dass ribosomale Proteine aufgrund ihres mittleren bis hohen Konservierungsgrades, ihrer Häufigkeit in kleineren EST-Projekten, damit verbunden ihrer Häufigkeit in Datenbanken und ihres phylogenetischen Informationsgehalts nützliche Werkzeuge für phylogenetische Fragestellungen sind. Es konnte durch phylogenetische Rekonstruktionen und Hypothesentests auf Basis eines 11.912 Aminosäuren langen Datensatzes gezeigt werden, dass die Xenoturbellida innerhalb der Deuterostomia eine Schwestergruppenbeziehung zu den Ambulacraria eingehen. Diese Arbeit zeigt im Vergleich aller bisher durchgeführten Arbeiten die beste statistische Unterstützung für diese Topologie. Weiterhin konnte untermauert werden, dass die Urochordata vermutlich anstelle der Cephalochordata die Schwestergruppe der Vertebrata sind. Der Vergleich der publizierten Xenoturbella EST-Datensätze mit dem eigenen Datensatz ließ den Rückschluß zu, dass ESTs offenbar klar weniger anfällig gegen Kontaminationen mit Erbmaterial (DNA+RNA) anderer Spezies sind als PCR-Amplifikate genomischer oder mitochondrialer Gene. Allerdings bestimmt anscheinend der physiologische Zustand der Tiere die Repräsentation von Transkriptklassen wie Stressproteine und mitochondriale Transkripte. Die bakteriellen Transkripte in einem der EST-Datensätze stammen vermutlich von Chlamydien, die möglicherweise symbiontisch in Xenoturbella bocki leben. Im Bereich der Protostomia wurden drei EST-Projekte für Vertreter der Syndermata durchgeführt. Basierend auf drei verschiedenen Proteinalignment-Datensätzen von ca. 11.000 Aminosäuren Länge konnte gezeigt werden, dass die Syndermata innerhalb der Spiralia einzugruppieren sind und dass sie mit den Gnathostomulida das monophyletische Supertaxon Gnathifera bilden. Die genaue phylogenetische Position der Syndermata innerhalb der Spiralia konnte hingegen noch nicht eindeutig geklärt werden, ebenso wie kein kongruenter Beweis für die Existenz des Supertaxons Platyzoa gefunden werden konnte. Im Rahmen der Untersuchung der internen Phylogenie der Syndermata konnten drei der fünf konkurrierenden Hypothesen aufgrund der Paraphylie der Eurotatoria ausgeschlossen werden. Da keine Daten der Seisonidea in den Analysen implementiert waren, bleibt die Frage der internen Phylogenie der Syndermata letztlich offen. Klar ist jedoch, dass die Eurotatoria nicht wie bislang angenommen monophyletisch sind, da die räderorgantragenden Bdelloidea keinesfalls den morphologisch diesbezüglich ähnlichen Monogononta ähnlich sind, sondern den räderorganlosen Acanthocephala näher stehen. Die Abbildung der molekularen Phylogenie auf die morphologischen Verhältnisse zeigt, dass das Räderorgan (partiell oder komplett) offenbar kurz nach der Aufspaltung der Syndermata in Monogononta und Acanthocephala + Bdelloidea in der Acanthocephala + Bdelloidea-Linie reduziert wurde. Die Entstehung des einziehbaren hinteren Körperteils (Rostrum bei Bdelloidea bzw. Proboscis bei Acanthocephala) in der Acanthocephala + Bdelloidea-Linie könnte das Schlüsselereignis zur Entstehung des Endoparasitismus der Acanthocephala gewesen sein.