New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0

PhyML is a phylogeny software based on the maximum-likelihood principle. Early PhyML versions used a fast algorithm performing nearest neighbor interchanges to improve a reasonable starting tree topology. Since the original publication (Guindon S., Gascuel O. 2003. A simple, fast and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52:696-704), PhyML has been widely used (>2500 citations in ISI Web of Science) because of its simplicity and a fair compromise between accuracy and speed. In the meantime, research around PhyML has continued, and this article describes the new algorithms and methods implemented in the program. First, we introduce a new algorithm to search the tree space with user-defined intensity using subtree pruning and regrafting topological moves. The parsimony criterion is used here to filter out the least promising topology modifications with respect to the likelihood function. The analysis of a large collection of real nucleotide and amino acid data sets of various sizes demonstrates the good performance of this method. Second, we describe a new test to assess the support of the data for internal branches of a phylogeny. This approach extends the recently proposed approximate likelihood-ratio test and relies on a nonparametric, Shimodaira-Hasegawa-like procedure. A detailed analysis of real alignments sheds light on the links between this new approach and the more classical nonparametric bootstrap method. Overall, our tests show that the last version (3.0) of PhyML is fast, accurate, stable, and ready to use. A Web server and binary files are available from http://www.atgc-montpellier.fr/phym

E-Biothon : Une plate-forme pour accélérer les recherches en biologie, santé et environnement

National audienceE-Biothon : Une plate-forme pour accélérer les recherches en biologie, santé et environnemen

The Metabolic Consequences of Hepatic AMP-Kinase Phosphorylation in Rainbow Trout

AMP-activated protein kinase (AMPK), a phylogenetically conserved serine/threonine protein kinase, is proposed to function as a “fuel gauge” to monitor cellular energy status in response to nutritional environmental variations. However, in fish, few studies have addressed the metabolic consequences related to the activation of this kinase. This study demonstrates that the rainbow trout (Oncorhynchus mykiss) possesses paralogs of the three known AMPK subunits that co-diversified, that the AMPK protein is present in the liver and in isolated hepatocytes, and it does change in response to physiological (fasting-re-feeding cycle) and pharmacological (AICAR and metformin administration and incubations) manipulations. Moreover, the phosphorylation of AMPK results in the phosphorylation of acetyl-CoA carboxylase, a main downstream target of AMPK in mammals. Other findings include changes in hepatic glycogen levels and several molecular actors involved in hepatic glucose and lipid metabolism, including mRNA transcript levels for glucokinase, glucose-6-phosphatase and fatty acid synthase both in vivo and in vitro. The fact that most results presented in this study are consistent with the recognized role of AMPK as a master regulator of energy homeostasis in living organisms supports the idea that these functions are conserved in this piscine model

Novel Arenavirus Sequences in Hylomyscus sp. and Mus (Nannomys) setulosus from Côte d'Ivoire: Implications for Evolution of Arenaviruses in Africa

This study aimed to identify new arenaviruses and gather insights in the evolution of arenaviruses in Africa. During 2003 through 2005, 1,228 small mammals representing 14 different genera were trapped in 9 villages in south, east, and middle west of Côte d'Ivoire. Specimens were screened by pan-Old World arenavirus RT-PCRs targeting S and L RNA segments as well as immunofluorescence assay. Sequences of two novel tentative species of the family Arenaviridae, Menekre and Gbagroube virus, were detected in Hylomyscus sp. and Mus (Nannomys) setulosus, respectively. Arenavirus infection of Mus (Nannomys) setulosus was also demonstrated by serological testing. Lassa virus was not found, although 60% of the captured animals were Mastomys natalensis. Complete S RNA and partial L RNA sequences of the novel viruses were recovered from the rodent specimens and subjected to phylogenetic analysis. Gbagroube virus is a closely related sister taxon of Lassa virus, while Menekre virus clusters with the Ippy/Mobala/Mopeia virus complex. Reconstruction of possible virus–host co-phylogeny scenarios suggests that, within the African continent, signatures of co-evolution might have been obliterated by multiple host-switching events

Early Evolution of Ionotropic GABA Receptors and Selective Regimes Acting on the Mammalian-Specific Theta and Epsilon Subunits

BACKGROUND: The amino acid neurotransmitter GABA is abundant in the central nervous system (CNS) of both invertebrates and vertebrates. Receptors of this neurotransmitter play a key role in important processes such as learning and memory. Yet, little is known about the mode and tempo of evolution of the receptors of this neurotransmitter. Here, we investigate the phylogenetic relationships of GABA receptor subunits across the chordates and detail their mode of evolution among mammals. PRINCIPAL FINDINGS: Our analyses support two major monophyletic clades: one clade containing GABA(A) receptor alpha, gamma, and epsilon subunits, and another one containing GABA(A) receptor rho, beta, delta, theta, and pi subunits. The presence of GABA receptor subunits from each of the major clades in the Ciona intestinalis genome suggests that these ancestral duplication events occurred before the divergence of urochordates. However, while gene divergence proceeded at similar rates on most receptor subunits, we show that the mammalian-specific subunits theta and epsilon experienced an episode of positive selection and of relaxed constraints, respectively, after the duplication event. Sites putatively under positive selection are placed on a three-dimensional model obtained by homology-modeling. CONCLUSIONS: Our results suggest an early divergence of the GABA receptor subunits, before the split from urochordates. We show that functional changes occurred in the lineages leading to the mammalian-specific subunit theta, and we identify the amino acid sites putatively responsible for the functional divergence. We discuss potential consequences for the evolution of mammals and of their CNS

The Genome of Borrelia recurrentis, the Agent of Deadly Louse-Borne Relapsing Fever, Is a Degraded Subset of Tick-Borne Borrelia duttonii

In an effort to understand how a tick-borne pathogen adapts to the body louse, we sequenced and compared the genomes of the recurrent fever agents Borrelia recurrentis and B. duttonii. The 1,242,163–1,574,910-bp fragmented genomes of B. recurrentis and B. duttonii contain a unique 23-kb linear plasmid. This linear plasmid exhibits a large polyT track within the promoter region of an intact variable large protein gene and a telomere resolvase that is unique to Borrelia. The genome content is characterized by several repeat families, including antigenic lipoproteins. B. recurrentis exhibited a 20.4% genome size reduction and appeared to be a strain of B. duttonii, with a decaying genome, possibly due to the accumulation of genomic errors induced by the loss of recA and mutS. Accompanying this were increases in the number of impaired genes and a reduction in coding capacity, including surface-exposed lipoproteins and putative virulence factors. Analysis of the reconstructed ancestral sequence compared to B. duttonii and B. recurrentis was consistent with the accelerated evolution observed in B. recurrentis. Vector specialization of louse-borne pathogens responsible for major epidemics was associated with rapid genome reduction. The correlation between gene loss and increased virulence of B. recurrentis parallels that of Rickettsia prowazekii, with both species being genomic subsets of less-virulent strains

Comparative Analysis of Acinetobacters: Three Genomes for Three Lifestyles

Acinetobacter baumannii is the source of numerous nosocomial infections in humans and therefore deserves close attention as multidrug or even pandrug resistant strains are increasingly being identified worldwide. Here we report the comparison of two newly sequenced genomes of A. baumannii. The human isolate A. baumannii AYE is multidrug resistant whereas strain SDF, which was isolated from body lice, is antibiotic susceptible. As reference for comparison in this analysis, the genome of the soil-living bacterium A. baylyi strain ADP1 was used. The most interesting dissimilarities we observed were that i) whereas strain AYE and A. baylyi genomes harbored very few Insertion Sequence elements which could promote expression of downstream genes, strain SDF sequence contains several hundred of them that have played a crucial role in its genome reduction (gene disruptions and simple DNA loss); ii) strain SDF has low catabolic capacities compared to strain AYE. Interestingly, the latter has even higher catabolic capacities than A. baylyi which has already been reported as a very nutritionally versatile organism. This metabolic performance could explain the persistence of A. baumannii nosocomial strains in environments where nutrients are scarce; iii) several processes known to play a key role during host infection (biofilm formation, iron uptake, quorum sensing, virulence factors) were either different or absent, the best example of which is iron uptake. Indeed, strain AYE and A. baylyi use siderophore-based systems to scavenge iron from the environment whereas strain SDF uses an alternate system similar to the Haem Acquisition System (HAS). Taken together, all these observations suggest that the genome contents of the 3 Acinetobacters compared are partly shaped by life in distinct ecological niches: human (and more largely hospital environment), louse, soil

Search for dark matter produced in association with bottom or top quarks in √s = 13 TeV pp collisions with the ATLAS detector

A search for weakly interacting massive particle dark matter produced in association with bottom or top quarks is presented. Final states containing third-generation quarks and miss- ing transverse momentum are considered. The analysis uses 36.1 fb−1 of proton–proton collision data recorded by the ATLAS experiment at √s = 13 TeV in 2015 and 2016. No significant excess of events above the estimated backgrounds is observed. The results are in- terpreted in the framework of simplified models of spin-0 dark-matter mediators. For colour- neutral spin-0 mediators produced in association with top quarks and decaying into a pair of dark-matter particles, mediator masses below 50 GeV are excluded assuming a dark-matter candidate mass of 1 GeV and unitary couplings. For scalar and pseudoscalar mediators produced in association with bottom quarks, the search sets limits on the production cross- section of 300 times the predicted rate for mediators with masses between 10 and 50 GeV and assuming a dark-matter mass of 1 GeV and unitary coupling. Constraints on colour- charged scalar simplified models are also presented. Assuming a dark-matter particle mass of 35 GeV, mediator particles with mass below 1.1 TeV are excluded for couplings yielding a dark-matter relic density consistent with measurements