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.

Sympatric speciation in palms on an oceanic island.

The origin of species diversity has challenged biologists for over two centuries. Allopatric speciation, the divergence of species resulting from geographical isolation, is well documented. However, sympatric speciation, divergence without geographical isolation, is highly controversial. Claims of sympatric speciation must demonstrate species sympatry, sister relationships, reproductive isolation, and that an earlier allopatric phase is highly unlikely. Here we provide clear support for sympatric speciation in a case study of two species of palm (Arecaceae) on an oceanic island. A large dated phylogenetic tree shows that the two species of Howea, endemic to the remote Lord Howe Island, are sister taxa and diverged from each other well after the island was formed 6.9 million years ago. During fieldwork, we found a substantial disjunction in flowering time that is correlated with soil preference. In addition, a genome scan indicates that few genetic loci are more divergent between the two species than expected under neutrality, a finding consistent with models of sympatric speciation involving disruptive/divergent selection. This case study of sympatric speciation in plants provides an opportunity for refining theoretical models on the origin of species, and new impetus for exploring putative plant and animal examples on oceanic islands.

Spatial predictions of phylogenetic diversity in conservation decision making.

Considering genetic relatedness among species has long been argued as an important step toward measuring biological diversity more accurately, rather than relying solely on species richness. Some researchers have correlated measures of phylogenetic diversity and species richness across a series of sites and suggest that values of phylogenetic diversity do not differ enough from those of species richness to justify their inclusion in conservation planning. We compared predictions of species richness and 10 measures of phylogenetic diversity by creating distribution models for 168 individual species of a species-rich plant family, the Cape Proteaceae. When we used average amounts of land set aside for conservation to compare areas selected on the basis of species richness with areas selected on the basis of phylogenetic diversity, correlations between species richness and different measures of phylogenetic diversity varied considerably. Correlations between species richness and measures that were based on the length of phylogenetic tree branches and tree shape were weaker than those that were based on tree shape alone. Elevation explained up to 31% of the segregation of species rich versus phylogenetically rich areas. Given these results, the increased availability of molecular data, and the known ecological effect of phylogenetically rich communities, consideration of phylogenetic diversity in conservation decision making may be feasible and informative.

Incompletely resolved phylogenetic trees inflate estimates of phylogenetic conservatism.

The tendency for more closely related species to share similar traits and ecological strategies can be explained by their longer shared evolutionary histories and represents phylogenetic conservatism. How strongly species traits co-vary with phylogeny can significantly impact how we analyze cross-species data and can influence our interpretation of assembly rules in the rapidly expanding field of community phylogenetics. Phylogenetic conservatism is typically quantified by analyzing the distribution of species values on the phylogenetic tree that connects them. Many phylogenetic approaches, however, assume a completely sampled phylogeny: while we have good estimates of deeper phylogenetic relationships for many species-rich groups, such as birds and flowering plants, we often lack information on more recent interspecific relationships (i.e., within a genus). A common solution has been to represent these relationships as polytomies on trees using taxonomy as a guide. Here we show that such trees can dramatically inflate estimates of phylogenetic conservatism quantified using S. P. Blomberg et al.'s K statistic. Using simulations, we show that even randomly generated traits can appear to be phylogenetically conserved on poorly resolved trees. We provide a simple rarefaction-based solution that can reliably retrieve unbiased estimates of K, and we illustrate our method using data on first flowering times from Thoreau's woods (Concord, Massachusetts, USA).

SlimCodeML: An Optimized Version of CodeML for the Branch-Site Model

CodeML (part of the PAML package) im- plements a maximum likelihood-based approach to de- tect positive selection on a specific branch of a given phylogenetic tree. While CodeML is widely used, it is very compute-intensive. We present SlimCodeML, an optimized version of CodeML for the branch-site model. Our performance analysis shows that SlimCodeML substantially outperforms CodeML (up to 9.38 times faster), especially for large-scale genomic analyses.

Plant functional and phylogenetic turnover correlate with climate and land use in the Western Swiss Alps

Understanding the relative importance of historical and environmental processes in the structure and composition of communities is one of the longest quests in ecological research. Increasingly, researchers are relying on the functional and phylogenetic β-diversity of natural communities to provide concise explanations on the mechanistic basis of community assembly and the drivers of trait variation among species. The present study investigated how plant functional and phylogenetic β-diversity change along key environmental and spatial gradients in the Western Swiss Alps. Methods Using the quadratic diversity measure based on six functional traits: specific leaf area (SLA), leaf dry matter content (LDMC), plant height (H), leaf carbon content (C), leaf nitrogen content (N), and leaf carbon to nitrogen content (C/N) alongside a species-resolved phylogenetic tree, we relate variations in climate, spatial geographic, land use and soil gradients to plant functional and phylogenetic turnover in mountain communities of the Western Swiss Alps. Important findings Our study highlights two main points. First, climate and land use factors play an important role in mountain plant community turnover. Second, the overlap between plant functional and phylogenetic turnover along these gradients correlates with the low phylogenetic signal in traits, suggesting that in mountain landscapes, trait lability is likely an important factor in driving plant community assembly. Overall, we demonstrate the importance of climate and land use factors in plant functional and phylogenetic community turnover, and provide valuable complementary insights into understanding patterns of β-diversity along several ecological gradients.

Prediction of direct and indirect genetic gains and genotypic correlations in rubber tree progenies

The objective of this work was to estimate the genetic parameters, genotypic and phenotypic correlations, and direct and indirect genetic gains among and within rubber tree (Hevea brasiliensis) progenies. The experiment was set up at the Municipality of Jaú, SP, Brazil. A randomized complete block design was used, with 22 treatments (progenies), 6 replicates, and 10 plants per plot at a spacing of 3x3 m. Three‑year‑old progenies were assessed for girth, rubber yield, and bark thickness by direct and indirect gains and genotypic correlations. The number of latex vessel rings showed the best correlations, correlating positively and significantly with girth and bark thickness. Selection gains among progenies were greater than within progeny for all the variables analyzed. Total gains obtained were high, especially for girth increase and rubber yield, which were 93.38 and 105.95%, respectively. Young progeny selection can maximize the expected genetic gains, reducing the rubber tree selection cycle.

In silico phylogenetic and virulence gene profile analyses of avian pathogenic Escherichia coli genome sequences

Avian pathogenic Escherichia coli (APEC) infections are responsible for significant losses in the poultry industry worldwide. A zoonotic risk has been attributed to APEC strains because they present similarities to extraintestinal pathogenic E. coli (ExPEC) associated with illness in humans, mainly urinary tract infections and neonatal meningitis. Here, we present in silico analyses with pathogenic E. coli genome sequences, including recently available APEC genomes. The phylogenetic tree, based on multi-locus sequence typing (MLST) of seven housekeeping genes, revealed high diversity in the allelic composition. Nevertheless, despite this diversity, the phylogenetic tree was able to cluster the different pathotypes together. An in silico virulence gene profile was also determined for each of these strains, through the presence or absence of 83 well-known virulence genes/traits described in pathogenic E. coli strains. The MLST phylogeny and the virulence gene profiles demonstrated a certain genetic similarity between Brazilian APEC strains, APEC isolated in the United States, UPEC (uropathogenic E. coli) and diarrheagenic strains isolated from humans. This correlation corroborates and reinforces the zoonotic potential hypothesis proposed to APEC.

Phylogenetic Analysis of Avidins and Expression of Novel Avidins from Lactrodectus hesperus and Hoeflea phototrophica

Avidins (Avds) are homotetrameric or homodimeric glycoproteins with typically less than 130 amino acid residues per monomer. They form a highly stable, non-covalent complex with biotin (vitamin H) with Kd = 10-15 M (for chicken Avd). The best-studied Avds are the chicken Avd from Gallus gallus and streptavidin from Streptomyces avidinii, although other Avd studies have also included Avds from various origins, e.g., from frogs, fishes, mushrooms and from many different bacteria. Several engineered Avds have been reported as well, e.g., dual-chain Avds (dcAvds) and single-chain Avds (scAvds), circular permutants with up to four simultaneously modifiable ligand-binding sites. These engineered Avds along with the many native Avds have potential to be used in various nanobiotechnological applications. In this study, we made a structure-based alignment representing all currently available sequences of Avds and studied the evolutionary relationship of Avds using phylogenetic analysis. First, we created an initial multiple sequence alignment of Avds using 42 closely related sequences, guided by the known Avd crystal structures. Next, we searched for non-redundant Avd sequences from various online databases, including National Centre for Biotechnology Information and the Universal Protein Resource; the identified sequences were added to the initial alignment to expand it to a final alignment of 242 Avd sequences. The MEGA software package was used to create distance matrices and a phylogenetic tree. Bootstrap reproducibility of the tree was poor at multiple nodes and may reflect on several possible issues with the data: the sequence length compared is relatively short and, whereas some positions are highly conserved and functional, others can vary without impinging on the structure or the function, so there are few informative sites; it may be that periods of rapid duplication have led to paralogs and that the differences among them are within the error limit of the data; and there may be other yet unknown reasons. Principle component analysis applied to alternative distance data did segregate the major groups, and success is likely due to the multivariate consideration of all the information. Furthermore, based on our extensive alignment and phylogenetic analysis, we expressed two novel Avds, lacavidin from Lactrodectus Hesperus, a western black widow spider, and hoefavidin from Hoeflea phototrophica, an aerobic marine bacterium, the ultimate aim being to determine their X-ray structures. These Avds were selected because of their unique sequences: lacavidin has an N-terminal Avd-like domain but a long C-terminal overhang, whereas hoefavidin was thought to be a dimeric Avd. Both these Avds could be used as novel scaffolds in biotechnological applications.

Diversité et structure du zooplancton : importance pour la conservation de la biodiversité aquatique en zone urbaine

L'écologie urbaine est un nouveau champ de recherche qui cherche à comprendre les structures et les patrons des communautés et des écosystèmes situés dans des paysages urbains. Les petits plans d’eau sont connus comme des écosystèmes aquatiques qui peuvent contenir une biodiversité considérable pour plusieurs groupes taxonomiques (oiseaux, amphibiens, macroinvertébrés), ce qui en fait des écosystèmes intéressants pour les études de conservation. Cependant, la biodiversité du zooplancton, un élément central des réseaux trophiques aquatiques, n’est pas entièrement connue pour les plans d’eaux urbains et devrait être mieux décrite et comprise. Cette étude a évalué les patrons de biodiversité des communautés zooplanctoniques dans des plans d’eau urbains sur l’Ile de Montréal et leurs sources de variation. Des suggestions pour l’évaluation et la conservation de la biodiversité sont aussi discutées. La biodiversité zooplanctonique des plans d’eaux urbains s’est avérée être assez élevée, avec les cladocères et les rotifères montrant les contributions à la diversité gamma et bêta les plus élevées. Sur l’ensemble des plans d’eau, il y avait une corrélation négative entre les contributions à la bêta diversité des cladocères et des rotifères. Au niveau de chaque plan d'eau, la zone littorale colonisée par des macrophytes s'est avérée être un habitat important pour la biodiversité zooplactonique, contribuant considérablement à la richesse en taxons, souvent avec une différente composition en espèces. Les communautés zooplanctoniques répondaient aux facteurs ascendants et descendants, mais aussi aux pratiques d’entretien, car le fait de vider les plans d’eau en hiver affecte la composition des communautés zooplanctoniques. Les communautés de cladocères dans ces plans d’eau possédaient des quantités variables de diversité phylogénétique, ce qui permet de les classer afin de prioriser les sites à préserver par rapport à la diversité phylogénétique. Le choix des sites à préserver afin de maximiser la diversité phylogénétique devrait être correctement établi, afin d’eviter de faire des choix sous-optimaux. Cependant, pour des taxons tels que les cladocères, pour lesquels les relations phylogénétiques demeurent difficiles à établir, placer une confiance absolue dans un seul arbre est une procédure dangereuse. L’incorporation de l’incertitude phylogénétique a démontré que, lorsqu’elle est prise en compte, plusieurs différences potentielles entre la diversité phylogenétique ne sont plus supportées. Les patrons de composition des communautés différaient entre les plans d’eau, les mois et les zones d’échantillonnage. Etant donné les intéractions sont significatives entres ces facters; ceci indique que tous ces facteurs devraient êtres considérés. L’urbanisation ne semblait pas sélectionner pour un type unique de composition des groupes alimentaires, étant donné que les communautés pouvaient changer entres des assemblages de types alimentaires différents. Les variables environnementales, surtout la couverture du plan d’eau en macrophytes, étaient des facteurs importants pour la biodiversité zooplanctonique, affectant la richesse spécifique de divers groupes taxonomiques et alimentaires. Ces variables affectaient aussi la composition des communautés, mais dans une moindre mesure, étant des variables explicatives modestes, ce qui indiquerait le besoin de considérer d’autres processus.

Plasmodium vivax lineages: geographical distribution, tandem repeat polymorphism, and phylogenetic relationship

Background: Multi-drug resistance and severe/ complicated cases are the emerging phenotypes of vivax malaria, which may deteriorate current anti-malarial control measures. The emergence of these phenotypes could be associated with either of the two Plasmodium vivax lineages. The two lineages had been categorized as Old World and New World, based on geographical sub-division and genetic and phenotypical markers. This study revisited the lineage hypothesis of P. vivax by typing the distribution of lineages among global isolates and evaluated their genetic relatedness using a panel of new mini-satellite markers. Methods: 18S SSU rRNA S-type gene was amplified from 420 Plasmodium vivax field isolates collected from different geographical regions of India, Thailand and Colombia as well as four strains each of P. vivax originating from Nicaragua, Panama, Thailand (Pak Chang), and Vietnam (ONG). A mini-satellite marker panel was then developed to understand the population genetic parameters and tested on a sample subset of both lineages. Results: 18S SSU rRNA S-type gene typing revealed the distribution of both lineages (Old World and New World) in all geographical regions. However, distribution of Plasmodium vivax lineages was highly variable in every geographical region. The lack of geographical sub-division between lineages suggests that both lineages are globally distributed. Ten mini-satellites were scanned from the P. vivax genome sequence; these tandem repeats were located in eight of the chromosomes. Mini-satellites revealed substantial allelic diversity (7-21, AE = 14.6 +/- 2.0) and heterozygosity (He = 0.697-0.924, AE = 0.857 +/- 0.033) per locus. Mini-satellite comparison between the two lineages revealed high but similar pattern of genetic diversity, allele frequency, and high degree of allele sharing. A Neighbour-Joining phylogenetic tree derived from genetic distance data obtained from ten mini-satellites also placed both lineages together in every cluster. Conclusions: The global lineage distribution, lack of genetic distance, similar pattern of genetic diversity, and allele sharing strongly suggested that both lineages are a single species and thus new emerging phenotypes associated with vivax malaria could not be clearly classified as belonging to a particular lineage on basis of their geographical origin.

Genes of cathepsin L-like proteases in Trypanosoma rangeli isolates: Markers for diagnosis, genotyping and phylogenetic relationships

We have sequenced genes encoding cathepsin L-like (CatL-like) cysteine proteases from isolates of Trypanosoma rangeli from humans, wild mammals and Rhodnius species of Central and South America. Phylogenetic trees of sequences encoding mature CatL-like enzymes of T rangeli and homologous genes from other trypanosomes, Leishmania spp. and bodonids positioned sequences of T rangeli (rangelipain) closest to T cruzi (cruzipain). Phylogenetic tree of kinetoplastids based on sequences of CatL-like was totally congruent with those derived from SSU rRNA and gGAPDH genes. Analysis of sequences from the CatL-like catalytic domains of 17 isolates representative of the overall phylogenetic diversity and geographical range of T rangeli supported all the lineages (A-D) previously defined using ribosomal and spliced leader genes. Comparison of the proteolytic activities of T rangeli isolates revealed heterogeneous banding profiles of cysteine proteases in gelatin gels, with differences even among isolates of the same lineage. CatL-like sequences proved to be excellent targets for diagnosis and genotyping of T rangeli by PCR. Data from CatL-like encoding genes agreed with results from previous studies of kDNA markers, and ribosomal and spliced leader genes, thereby corroborating clonal evolution, independent transmission cycles and the divergence of T rangeli lineages associated with sympatric species of Rhodnius. (c) 2009 Elsevier B.V. All rights reserved.

Sequence, evolution and ligand binding properties of mammalian Duffy antigen/receptor for chemokines

The Duffy antigen/receptor for chemokine, DARC, acts as a widely expressed promiscuous chemokine receptor and as the erythrocyte receptor for Plasmodium vivax. To gain insight into the evolution and structure/function relations of DARC, we analyzed the binding of anti-human Fy monoclonal antibodies (mAbs) and human chemokines to red blood cells (RBCs) from 11 nonhuman primates and two nonprimate mammals, and we elucidated the structures of the DARC genes from gorilla, gibbon, baboon, marmoset, tamarin, night monkey and cattle. CXCL-8 and CCL-5 chemokine binding analysis indicated that the promiscuous binding profile characteristic of DARC is conserved across species. Among three mAbs that detected the Fy6 epitope by flow cytometric analysis of human and chimpanzee RBCs, only one reacted with night monkey and squirrel monkey. Only chimpanzee RBCs bound a significant amount of the anti-Fy3 mAb. Fy3 was also poorly detected on RBCs from gorilla, baboon and rhesus monkey, but not from new world monkeys. Alignment of DARC homologous sequences allowed us to construct a phylogenetic tree in which all branchings were in accordance with current knowledge of primate phylogeny. Although DARC was expected to be under strong internal and external selection pressure, in order to maintain chemokine binding and avoid Plasmodium vivax binding, respectively, our present study did not provide arguments in favor of a selection pressure on the extracellular domains involved in ligand specificity. The amino acid variability of DARC-like polypeptides was found to be well correlated with the hydrophylicity indexes, with the highest divergence on the amino-terminal extracellular domain. Analysis of the deduced amino acid sequences highlighted the conservation of some amino acid residues, which should prove to be critical for the structural and functional properties of DARC.

Analysis of HCV quasispecies dynamic under selective pressure of combined therapy

Background: The quasispecies composition of Hepatitis C virus (HCV) could have important implications with regard to viral persistence and response to interferon-based therapy. The complete NS5A was analyzed to evaluate whether the composition of NS5A quasispecies of HCV 1a/1b is related to responsiveness to combined interferon pegylated (PEG-IFN) and ribavirin therapy.Methods: Viral RNA was isolated from serum samples collected before, during and after treatment from virological sustained responder (SVR), non-responder (NR) and the end-of-treatment responder patients (ETR). NS5A region was amplified, cloned and sequenced. Six hundred and ninety full-length NS5A sequences were analyzed.Results: This study provides evidence that lower nucleotide diversity of the NS5A region pre-therapy is associated with viral clearance. Analysis of samples of NRs and the ETRs time points showed that genetic diversity of populations tend to decrease over time. Post-therapy population of ETRs presented higher genetic distance from baseline probably due to the bottleneck phenomenon observed for those patients in the end of treatment. The viral effective population of those patients also showed a strong decrease after therapy. Otherwise, NRs demonstrated a continuous variation or stability of effective populations and genetic diversity over time that did not seem to be related to therapy. Phylogenetic relationships concerning complete NS5A sequences obtained from patients did not demonstrate clustering associated with specific response patterns. However, distinctive clustering of pre/post-therapy sequences was observed. In addition, the evolution of quasispecies over time was subjected to purifying or relaxed purifying selection. Codons 157 (P03), 182 and 440 (P42), 62 and 404 (P44) were found to be under positive selective pressure but it failed to be related to the therapy.Conclusion: These results confirm the hypothesis that a relationship exists between NS5A heterogeneity and response to therapy in patients infected with chronic hepatitis C. © 2013 Jardim et al.; licensee BioMed Central Ltd.

Circulation of Different Lineages of Dengue Virus 2, Genotype American/Asian in Brazil: Dynamics and Molecular and Phylogenetic Characterization

The American/Asian genotype of Dengue virus type 2 (DENV-2) was introduced into the Americas in the 80′s. Although there is no data showing when this genotype was first introduced into Brazil, it was first detected in Brazil in 1990. After which the virus spread throughout the country and major epidemics occurred in 1998, 2007/08 and 2010. In this study we sequenced 12 DENV-2 genomes obtained from serum samples of patients with dengue fever residing in São José do Rio Preto, São Paulo (SJRP/SP), Brazil, in 2008. The whole open reading frame or envelope sequences were used to perform phylogenetic, phylogeographic and evolutionary analyses. Isolates from SJRP/SP were grouped within one lineage (BR3) close to isolates from Rio de Janeiro, Brazil. Isolates from SJRP were probably introduced there at least in 2007, prior to its detection in the 2008 outbreak. DENV-2 circulation in Brazil is characterized by the introduction, displacement and circulation of three well-defined lineages in different times, most probably from the Caribbean. Thirty-seven unique amino acid substitutions were observed among the lineages, including seven amino acid differences in domains I to III of the envelope protein. Moreover, we dated here, for the first time, the introduction of American/Asian genotype into Brazil (lineage BR1) to 1988/89, followed by the introduction of lineages BR2 (1998-2000) and BR3 (2003-05). Our results show a delay between the introduction and detection of DENV-2 lineages in Brazil, reinforcing the importance and need for surveillance programs to detect and trace the evolution of these viruses. Additionally, Brazilian DENV-2 differed in genetic diversity, date of introduction and geographic origin and distribution in Brazil, and these are important factors for the evolution, dynamics and control of dengue. © 2013 Drumond et al.