Evolutionary History of the Plant Pathogenic Bacterium Xanthomonas axonopodis Nadia Mhedbi-Hajri 1,2,3 , Ahmed Hajri 1,2,3 , Tristan Boureau 1,2,3 , Armelle Darrasse 1,2,3 , Karine Durand 1,2,3 , Chrystelle Brin 1,2,3 , Marion Fischer-Le Saux 1,2,3 , Charles Manceau 1,2,3 , Ste ´ phane Poussier 1,2,3 , Olivier Pruvost 4 , Christophe Lemaire 1,2,3. , Marie-Agne ` s Jacques 1,2,3 * . 1 INRA, UMR1345 Institut de Recherche en Horticulture et Semences, Beaucouze ´, France, 2 Universite ´ d’Angers, UMR1345 Institut de Recherche en Horticulture et Semences, SFR4207 QUASAV, PRES L’UNAM, Angers, France, 3 AgroCampus-Ouest, UMR1345 Institut de Recherche en Horticulture et Semences, Angers, France, 4 CIRAD, UMR PVBMT, Saint Pierre, Cedex Re ´ union, France Abstract Deciphering mechanisms shaping bacterial diversity should help to build tools to predict the emergence of infectious diseases. Xanthomonads are plant pathogenic bacteria found worldwide. Xanthomonas axonopodis is a genetically heterogeneous species clustering, into six groups, strains that are collectively pathogenic on a large number of plants. However, each strain displays a narrow host range. We address the question of the nature of the evolutionary processes – geographical and ecological speciation – that shaped this diversity. We assembled a large collection of X. axonopodis strains that were isolated over a long period, over continents, and from various hosts. Based on the sequence analysis of seven housekeeping genes, we found that recombination occurred as frequently as point mutation in the evolutionary history of X. axonopodis. However, the impact of recombination was about three times greater than the impact of mutation on the diversity observed in the whole dataset. We then reconstructed the clonal genealogy of the strains using coalescent and genealogy approaches and we studied the diversification of the pathogen using a model of divergence with migration. The suggested scenario involves a first step of generalist diversification that spanned over the last 25 000 years. A second step of ecology-driven specialization occurred during the past two centuries. Eventually, secondary contacts between host- specialized strains probably occurred as a result of agricultural development and intensification, allowing genetic exchanges of virulence-associated genes. These transfers may have favored the emergence of novel pathotypes. Finally, we argue that the largest ecological entity within X. axonopodis is the pathovar. Citation: Mhedbi-Hajri N, Hajri A, Boureau T, Darrasse A, Durand K, et al. (2013) Evolutionary History of the Plant Pathogenic Bacterium Xanthomonas axonopodis. PLoS ONE 8(3): e58474. doi:10.1371/journal.pone.0058474 Editor: Keith A. Crandall, George Washington University, United States of America Received September 18, 2012; Accepted February 5, 2013; Published March 7, 2013 Copyright: ß 2013 Mhedbi-Hajri et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: N. Mhedbi-Hajri and A. Hajri were supported by grants from the Tunisian Government, from CG Maine et Loire and from INRA. This study was funded by the Region Pays de la Loire (Xanthost) and the Plant Health and Environment Department of INRA (SPE 2007-0077-01, AIP Bioressources 2008-077-01). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. * E-mail: marie-agnes.jacques@angers.inra.fr . These authors contributed equally to this work. Introduction Determining the relative part played by geography and ecological specialization in the reproductive isolation between populations of pathogens is still a great challenge [1], [2]. The extent of isolation between pathogen populations is conditioned by the divergence accumulated in allopatry or by the degree of adaptation to hosts. In complete allopatry, genetic reproductive barriers are less likely to be selected than in sympatry [3], and gene flow remains possible for populations that have diverged for a long time (several million years) [4]. Conversely, specialization on hosts greatly contributes to reproductive isolation and limits gene flow. Host specialization was promoted by plant domestication, which led to a reduction of host genetic variability, and by the development of agriculture, which has uniformized crops and increased host density [5]. More recently, the globalization of agriculture has strongly contributed to breaking natural barriers to dispersal, thereby reducing the geographic isolation and increasing opportunities for gene flow. Subsequent propagation of many pathogens allowed secondary contacts between populations, which have diverged in allopatry [6]. Would such a return to sympatry be sufficient to promote gene flow between these populations of pathogens? The study of the relative importance of gene flow and divergence times between populations occupying different hosts then becomes crucial for understanding evolutionary histories and emergences of pathogens [7], [8]. Unlike for sexual eukaryotes, the classical biological species definition does not apply for prokaryotes, which are asexually reproducing organisms. Bacteria, rather, form ecologically isolated units called ecotypes [9]. The pathovar is an infrasubspecific division that was created to group pathogenic bacteria that display the same symptomatology on the same host range [10]. Would the ecological cohesion of the ecotype concept be found in the pathovar concept? The extent of reproductive isolation among pathovars remains to be clarified, because ecological isolation may be associated with reproductive isolation. Bacteria are mainly clonal but can exchange genes or fragments of their genomes through horizontal gene transfer. Such transfers would be PLOS ONE | www.plosone.org 1 March 2013 | Volume 8 | Issue 3 | e58474