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.

Phylogenetic inference's algorithms

Phylogenetic inference consist in the search of an evolutionary tree to explain the best way possible genealogical relationships of a set of species. Phylogenetic analysis has a large number of applications in areas such as biology, ecology, paleontology, etc. There are several criterias which has been defined in order to infer phylogenies, among which are the maximum parsimony and maximum likelihood. The first one tries to find the phylogenetic tree that minimizes the number of evolutionary steps needed to describe the evolutionary history among species, while the second tries to find the tree that has the highest probability of produce the observed data according to an evolutionary model. The search of a phylogenetic tree can be formulated as a multi-objective optimization problem, which aims to find trees which satisfy simultaneously (and as much as possible) both criteria of parsimony and likelihood. Due to the fact that these criteria are different there won't be a single optimal solution (a single tree), but a set of compromise solutions. The solutions of this set are called "Pareto Optimal". To find this solutions, evolutionary algorithms are being used with success nowadays.This algorithms are a family of techniques, which aren’t exact, inspired by the process of natural selection. They usually find great quality solutions in order to resolve convoluted optimization problems. The way this algorithms works is based on the handling of a set of trial solutions (trees in the phylogeny case) using operators, some of them exchanges information between solutions, simulating DNA crossing, and others apply aleatory modifications, simulating a mutation. The result of this algorithms is an approximation to the set of the “Pareto Optimal” which can be shown in a graph with in order that the expert in the problem (the biologist when we talk about inference) can choose the solution of the commitment which produces the higher interest. In the case of optimization multi-objective applied to phylogenetic inference, there is open source software tool, called MO-Phylogenetics, which is designed for the purpose of resolving inference problems with classic evolutionary algorithms and last generation algorithms. REFERENCES [1] C.A. Coello Coello, G.B. Lamont, D.A. van Veldhuizen. Evolutionary algorithms for solving multi-objective problems. Spring. Agosto 2007 [2] C. Zambrano-Vega, A.J. Nebro, J.F Aldana-Montes. MO-Phylogenetics: a phylogenetic inference software tool with multi-objective evolutionary metaheuristics. Methods in Ecology and Evolution. En prensa. Febrero 2016.

Phylogenetic Analysis Of The Microbial Community In Hypersaline Petroleum Produced Water From The Campos Basin.

In this work the archaea and eubacteria community of a hypersaline produced water from the Campos Basin that had been transported and discharged to an onshore storage facility was evaluated by 16S recombinant RNA (rRNA) gene sequence analysis. The produced water had a hypersaline salt content of 10 (w/v), had a carbon oxygen demand (COD) of 4,300 mg/l and contains phenol and other aromatic compounds. The high salt and COD content and the presence of toxic phenolic compounds present a problem for conventional discharge to open seawater. In previous studies, we demonstrated that the COD and phenolic content could be largely removed under aerobic conditions, without dilution, by either addition of phenol degrading Haloarchaea or the addition of nutrients alone. In this study our goal was to characterize the microbial community to gain further insight into the persistence of reservoir community members in the produced water and the potential for bioremediation of COD and toxic contaminants. Members of the archaea community were consistent with previously identified communities from mesothermic reservoirs. All identified archaea were located within the phylum Euryarchaeota, with 98 % being identified as methanogens while 2 % could not be affiliated with any known genus. Of the identified archaea, 37 % were identified as members of the strictly carbon-dioxide-reducing genus Methanoplanus and 59 % as members of the acetoclastic genus Methanosaeta. No Haloarchaea were detected, consistent with the need to add these organisms for COD and aromatic removal. Marinobacter and Halomonas dominated the eubacterial community. The presence of these genera is consistent with the ability to stimulate COD and aromatic removal with nutrient addition. In addition, anaerobic members of the phyla Thermotogae, Firmicutes, and unclassified eubacteria were identified and may represent reservoir organisms associated with the conversion hydrocarbons to methane.

Phylogenetic And Phylogeographic Mapping Of The Avian Coronavirus Spike Protein-encoding Gene In Wild And Synanthropic Birds.

The evolution and population dynamics of avian coronaviruses (AvCoVs) remain underexplored. In the present study, in-depth phylogenetic and Bayesian phylogeographic studies were conducted to investigate the evolutionary dynamics of AvCoVs detected in wild and synanthropic birds. A total of 500 samples, including tracheal and cloacal swabs collected from 312 wild birds belonging to 42 species, were analysed using molecular assays. A total of 65 samples (13%) from 22 bird species were positive for AvCoV. Molecular evolution analyses revealed that the sequences from samples collected in Brazil did not cluster with any of the AvCoV S1 gene sequences deposited in the GenBank database. Bayesian framework analysis estimated an AvCoV strain from Sweden (1999) as the most recent common ancestor of the AvCoVs detected in this study. Furthermore, the analysis inferred an increase in the AvCoV dynamic demographic population in different wild and synanthropic bird species, suggesting that birds may be potential new hosts responsible for spreading this virus.

The Role Of Edaphic Environment And Climate In Structuring Phylogenetic Pattern In Seasonally Dry Tropical Plant Communities.

Seasonally dry tropical plant formations (SDTF) are likely to exhibit phylogenetic clustering owing to niche conservatism driven by a strong environmental filter (water stress), but heterogeneous edaphic environments and life histories may result in heterogeneity in degree of phylogenetic clustering. We investigated phylogenetic patterns across ecological gradients related to water availability (edaphic environment and climate) in the Caatinga, a SDTF in Brazil. Caatinga is characterized by semiarid climate and three distinct edaphic environments - sedimentary, crystalline, and inselberg -representing a decreasing gradient in soil water availability. We used two measures of phylogenetic diversity: Net Relatedness Index based on the entire phylogeny among species present in a site, reflecting long-term diversification; and Nearest Taxon Index based on the tips of the phylogeny, reflecting more recent diversification. We also evaluated woody species in contrast to herbaceous species. The main climatic variable influencing phylogenetic pattern was precipitation in the driest quarter, particularly for herbaceous species, suggesting that environmental filtering related to minimal periods of precipitation is an important driver of Caatinga biodiversity, as one might expect for a SDTF. Woody species tended to show phylogenetic clustering whereas herbaceous species tended towards phylogenetic overdispersion. We also found phylogenetic clustering in two edaphic environments (sedimentary and crystalline) in contrast to phylogenetic overdispersion in the third (inselberg). We conclude that while niche conservatism is evident in phylogenetic clustering in the Caatinga, this is not a universal pattern likely due to heterogeneity in the degree of realized environmental filtering across edaphic environments. Thus, SDTF, in spite of a strong shared environmental filter, are potentially heterogeneous in phylogenetic structuring. Our results support the need for scientifically informed conservation strategies in the Caatinga and other SDTF regions that have not previously been prioritized for conservation in order to take into account this heterogeneity.

Systematics of Lamontichthys Miranda-Ribeiro (Siluriformes: Loricariidae), with the description of two new species

The taxonomic revision of the genus Lamontichthys Miranda-Ribeiro, based on the examination of 164 specimens of different river drainages throughout the Amazon basin, revealed the presence of six species of which two are new. Lamontichthys filamentosus occurs in the upper and middle portions of the rio Amazonas basin; L. llanero in the río Orinoco basin; L. maracaibero in the lago Maracaibo basin; and L. stibaros in the upper río Amazonas basin. Lamontichthys avacanoeiro, new species, occurs in the upper rio Tocantins basin; and L. parakana, new species, in the lower rio Tocantins basin. The new species represent a considerable extension in the so far known distribution of the genus. A parsimony analysis, including 87 osteological and external morphological characters from Lamontichthys and related taxa (total of 16), resulted in three most parsimonious trees with 194 steps (CI = 0.73 and RI = 0.78). The hypothesis of monophyly of Lamontichthys is corroborated and supported by six derived characters. Within Lamontichthys two monophyletic assemblages are recognized, one includes L. avacanoeiro and L. stibaros, the other includes L. maracaibero and the clade formed by L. filamentosus and L. llanero. The relationships of Lamontichthys parakana, a species that was not included in the phylogenetic analysis is discussed. The monophyly and relationships of the monotypic genus Pterosturisoma microps are also discussed.

A phylogenetic analysis of Brycon and Henochilus (Characiformes, Characidae, Bryconinae) based on the mitochondrial gene 16S rRNA

The genus Brycon, the largest subunit of the Bryconinae, has 42 valid species distributed from southern Mexico to the La Plata River in Argentina. Henochilus is a monotypic genus, comprising a single species (H. wheatlandii) found in the upper Rio Doce basin. In the present study, partial sequences of the mitochondrial gene 16S were obtained for fifteen species of Brycon and for Henochilus wheatlandii. The results showed that the genus Brycon is paraphyletic, since Henochilus is the sister-group of B. ferox and B. insignis. The most basal species analyzed were the trans-Andean species B. henni, B. petrosus, and B. chagrensis.

A new species of sand-dwelling catfish, with a phylogenetic diagnosis of Pygidianops Myers (Siluriformes: Trichomycteridae: Glanapteryginae)

A new species of sand-dwelling catfish genus Pygidianops, P. amphioxus, is described from the Negro and lower Amazon basins. The new species differs from its three congeners in the elongate eel-like body, the short barbels, and the small caudal fin, continuous with the body, among other traits of internal anatomy. The absence of anal fin further distinguishes P. amphioxus from all other Pygidianops species except P. magoi and the presence of eyes from all except P. cuao. The new Pygidianops seems to be the sister species to P. magoi, the two species sharing a unique mesethmoid with a dorsally-bent tip lacking cornua, and a produced articular process in the palatine for the articulation with the neurocranium. Pygidianops amphioxus is a permanent and highly-specialized inhabitant of psammic environments. Additional characters are proposed as synapomorphies of Pygidianops, including a hypertrophied symphyseal joint and associated ligament in the lower jaw; an elongate, laterally-directed, process on the dorsal surface of the premaxilla; and a rotated lower jaw, where the surface normally facing laterally in other glanapterygines is instead directed ventrally. These and other characters are incorporated into a revised phylogenetic diagnosis of Pygidianops.

A phylogenetic study of the subtribe Dicrepidiina (Elateridae, Elaterinae, Ampedini)

It is presented a cladistic analysis of the Dicrepidiina aiming to test the monophyletism of the subtribe and to establish the relationships among the genera. The subtribe is composed by 36 genera and all of them, except Asebis, Lamononia, Neopsephus, Semiotopsis and Spilomorphus were included in the analysis. Fifty two species, especially the type-species of each genus were studied: Achrestus flavocinctus (Candèze, 1859), A. venustus Champion, 1895, Adiaphorus gracilis Schwarz, 1901, A. ponticerianus Candèze, 1859, Anoplischiopsis bivittatus Champion, 1895, Anoplischius bicarinatus Candèze, 1859, A. conicus Candèze, 1900, A. haematopus Candèze, 1859, A. pyronotus Candèze, 1859, Atractosomus flavescens (Germar, 1839), Blauta cribraria (Germar, 1844), Calopsephus apicalis (Schwarz, 1903), Catalamprus angustus (Fleutiaux, 1902), Crepidius flabellifer (Erichson, 1847), C. resectus Candèze, 1859, Cyathodera auripilosus Costa, 1968, C. lanugicollis (Candèze, 1859), C. longicornis Blanchard, 1843, Dayakus angularis Candèze, 1893, Dicrepidius ramicornis (Palisot de Beauvois, 1805), Dipropus brasilianus (Germar, 1824), D. factuellus Candèze, 1859, D. laticollis (Eschscholtz, 1829), D. pinguis (Candèze, 1859), D. schwarzi (Becker, 1961), Elius birmanicus Candèze, 1893, E. dilatatus Candèze, 1878, Heterocrepidius gilvellus Candèze, 1859, H. ventralis Guérin-Méneville, 1838, Lampropsephus cyaneus (Candèze, 1878), Loboederus appendiculatus (Perty, 1830), Olophoeus gibbus Candèze, 1859, Ovipalpus pubescens Solier, 1851, Pantolamprus ligneus Candèze, 1896, P. mirabilis Candèze, 1896, P. perpulcher Westwood, 1842, Paraloboderus glaber Golbach, 1990, Proloboderus crassipes Fleutiaux, 1912, Propsephus beniensis (Candèze, 1859), P. cavifrons (Erichson, 1843), Pseudolophoeus guineensis (Candèze, 1881), Rhinopsephus apicalis (Schwarz, 1903), Sephilus formosanus Schwarz, 1912, S. frontalis Candèze, 1878, Singhalenus gibbus Candèze, 1892, S. taprobanicus Candèze, 1859, Sphenomerus antennalis Candèze, 1859, S. brunneus Candèze, 1865, Spilus atractomorphus Candèze, 1859, S. nitidus Candèze, 1859, Stenocrepidius simonii Fleutiaux, 1891 and Trielasmus varians Blanchard, 1846. Chalcolepidius zonatus (Hemirhipini, Agrypninae), Ctenicera silvatica (Prosternini, Prosterninae), and species of the other subtribes of Ampedini (Elaterinae): Ampedus sanguineus (Ampedina), Melanotus spernendus (Melanotina) and Anchastus digittatus and Physorhinus xanthocephalus (Physorhinina) were used as outgroups. The results of the phylogenetic analysis demonstrated that Dicrepidiina, as formerly defined, does not form a monophyletic group. One genus, represented by Ovipalpus pubescens, was removed from the subtribe. The subtribe is characterized by presence of lamella under 2nd and 3rd tarsomeres of all legs. Also, it was revealed that the genera Achrestus, Anoplischius, Dipropus and Propsephus are not monophyletic. Due to the scarcity of information, all the studied species are redescribed and illustrated.

Description and phylogenetic analysis of the Calycopidina (Lepidoptera, Lycaenidae, Theclinae, Eumaeini): a subtribe of detritivores

Description and phylogenetic analysis of the Calycopidina (Lepidoptera, Lycaenidae, Theclinae, Eumaeini): a subtribe of detritivores. The purpose of this paper is to establish a phylogenetic basis for a new Eumaeini subtribe that includes those lycaenid genera in which detritivory has been recorded. Morphological characters were coded for 82 species of the previously proposed "Lamprospilus Section" of the Eumaeini (19 of these had coding identical to another species), and a phylogenetic analysis was performed using the 63 distinct ingroup terminal taxa and six outgroups belonging to four genera. Taxonomic results include the description in the Eumaeini of Calycopidina Duarte & Robbins new subtribe (type genus Calycopis Scudder, 1876), which contains Lamprospilus Geyer, Badecla Duarte & Robbins new genus (type species Thecla badaca Hewitson), Arzecla Duarte & Robbins new genus (type species Thecla arza Hewitson), Arumecla Robbins & Duarte, Camissecla Robbins & Duarte, Electrostrymon Clench, Rubroserrata K. Johnson & Kroenlein revalidated status, Ziegleria K. Johnson, Kisutam K. Johnson & Kroenlein revalidated status, and Calycopis. Previous "infratribe" names Angulopina K. Johnson & Kroenlein, 1993, and Calycopina K. Johnson & Kroenlein, 1993, are nomenclaturally unavailable and polyphyletic as proposed. New combinations include Badecla badaca (Hewitson), Badecla picentia (Hewitson), Badecla quadramacula (Austin & K. Johnson), Badecla lanckena (Schaus), Badecla argentinensis (K. Johnson & Kroenlein), Badecla clarissa (Draudt), Arzecla arza (Hewitson), Arzecla tarpa (Godman & Salvin), Arzecla canacha (Hewitson), Arzecla calatia (Hewitson), Arzecla tucumanensis (K. Johnson & Kroenlein), Arzecla sethon (Godman & Salvin), Arzecla nubilum (H. H. Druce), Arzecla paralus (Godman & Salvin), Arzecla taminella (Schaus), Arzecla albolineata (Lathy), Electrostrymon denarius (Butler & H.Druce), Electrostrymon guzanta (Schaus), Electrostrymon perisus (H. H. Druce), Rubroserrata mathewi (Hewitson), Rubroserrata ecbatana (Hewitson), Kisutam micandriana (K. Johnson), and Kisutam syllis (Godman & Salvin). The structure of the male genitalia lateral window, labides, and brush organs are described and discussed, as are the female genitalia signa of the corpus bursae and 8th abdominal tergum. Widespread wing pattern sexual dimorphism in the Calycopidina is noted and illustrated, and the presence of alternating dark and light bands on the ventral wings of both sexes is discussed. The evidence for detritivory in Lamprospilus, Badecla, Arzecla, Arumecla, Camissecla, Electrostrymon, Ziegleria, Kisutam, and Calycopis is summarized using the new classification.

Phylogenetic placement of the enigmatic parasite, Polypodium hydriforme, within the Phylum Cnidaria

Background: Polypodium hydriforme is a parasite with an unusual life cycle and peculiar morphology, both of which have made its systematic position uncertain. Polypodium has traditionally been considered a cnidarian because it possesses nematocysts, the stinging structures characteristic of this phylum. However, recent molecular phylogenetic studies using 18S rDNA sequence data have challenged this interpretation, and have shown that Polypodium is a close relative to myxozoans and together they share a closer affinity to bilaterians than cnidarians. Due to the variable rates of 18S rDNA sequences, these results have been suggested to be an artifact of long-branch attraction ( LBA). A recent study, using multiple protein coding markers, shows that the myxozoan Buddenbrockia, is nested within cnidarians. Polypodium was not included in this study. To further investigate the phylogenetic placement of Polypodium, we have performed phylogenetic analyses of metazoans with 18S and partial 28S rDNA sequences in a large dataset that includes Polypodium and a comprehensive sampling of cnidarian taxa. Results: Analyses of a combined dataset of 18S and partial 28S sequences, and partial 28S alone, support the placement of Polypodium within Cnidaria. Removal of the long-branched myxozoans from the 18S dataset also results in Polypodium being nested within Cnidaria. These results suggest that previous reports showing that Polypodium and Myxozoa form a sister group to Bilateria were an artifact of long-branch attraction. Conclusion: By including 28S rDNA sequences and a comprehensive sampling of cnidarian taxa, we demonstrate that previously conflicting hypotheses concerning the phylogenetic placement of Polypodium can be reconciled. Specifically, the data presented provide evidence that Polypodium is indeed a cnidarian and is either the sister taxon to Hydrozoa, or part of the hydrozoan clade, Leptothecata. The former hypothesis is consistent with the traditional view that Polypodium should be placed in its own cnidarian class, Polypodiozoa.

WOOD DENSITY VARIATION AND TREE RING DEMARCATION IN Gmelina arborea TREES USING X-RAY DENSITOMETRY

Due to its relationship with other properties, wood density is the main wood quality parameter. Modern, accurate methods - such as X-ray densitometry - are applied to determine the spatial distribution of density in wood sections and to evaluate wood quality. The objectives of this study were to determinate the influence of growing conditions on wood density variation and tree ring demarcation of gmelina trees from fast growing plantations in Costa Rica. The wood density was determined by X-ray densitometry method. Wood samples were cut from gmelina trees and were exposed to low X-rays. The radiographic films were developed and scanned using a 256 gray scale with 1000 dpi resolution and the wood density was determined by CRAD and CERD software. The results showed tree-ring boundaries were distinctly delimited in trees growing in site with rainfall lower than 25 10 mm/year. It was demonstrated that tree age, climatic conditions and management of plantation affects wood density and its variability. The specific effect of variables on wood density was quantified by for multiple regression method. It was determined that tree year explained 25.8% of the total variation of density and 19.9% were caused by climatic condition where the tree growing. Wood density was less affected by the intensity of forest management with 5.9% of total variation.

Variation in the wood anatomical structure of Gmelina arborea (Verbenaceae) trees at different ecological conditions in Costa Rica

The tree Gmelina arborea has been widely introduced in Costa Rica for commercial purposes. This new conditions for melina cause variations on anatomy in secondary xylem of the trees growing in plantations. The objective of the present research was to determine the variation in the anatomy of xylem caused by the ecological conduction variation. Dimensions of fiber, axial parenchyma percentage of cross sections, parameters of vessels and the ray were measured. The results showed that some anatomical characteristics remained stable despite variations of ecological conditions, especially radial parenchyma and anatomical features which were less affected by the altitude. On the other hand, the vessels, axial parenchyma and fiber were less stable because they were affected significantly by the longitude, latitude, altitude and precipitation. Latitude significantly affected vessel percentage, length and diameter of the fiber and lumen. Longitude affected vessel percentage and fiber diameter. Altitude had a significant correlation with the amount of cells at my height. Annual average precipitation affected vessel percentage and diameter, not only of the fiber, but also of the lumen. These results suggest that the new growth conditions of G. arborea trees in Costa Rica have produced an anatomic adaptation.