Persistent homology reveals strong phylogenetic signal in 3D protein structures.

Changes that occur in proteins over time provide a phylogenetic signal that can be used to decipher their evolutionary history and the relationships between organisms. Sequence comparison is the most common way to access this phylogenetic signal, while those based on 3D structure comparisons are still in their infancy. In this study, we propose an effective approach based on Persistent Homology Theory (PH) to extract the phylogenetic information contained in protein structures. PH provides efficient and robust algorithms for extracting and comparing geometric features from noisy datasets at different spatial resolutions. PH has a growing number of applications in the life sciences, including the study of proteins (e.g. classification, folding). However, it has never been used to study the phylogenetic signal they may contain. Here, using 518 protein families, representing 22,940 protein sequences and structures, from 10 major taxonomic groups, we show that distances calculated with PH from protein structures correlate strongly with phylogenetic distances calculated from protein sequences, at both small and large evolutionary scales. We test several methods for calculating PH distances and propose some refinements to improve their relevance for addressing evolutionary questions. This work opens up new perspectives in evolutionary biology by proposing an efficient way to access the phylogenetic signal contained in protein structures, as well as future developments of topological analysis in the life sciences.

Description of a new genus of the Pectobacteriaceae family isolated from water in coastal brackish wetlands of the French Camargue region, Prodigiosinella gen. nov., including the new species Prodigiosinella aquatilis sp. nov.

The Pectobacteriaceae family comprises plant pathogens able to provoke diverse diseases, including plant maceration due to the production of pectinases disrupting the plant cell wall. To better understand their diversity, a survey of pectinolytic bacteria was performed in brackish lakes of the French region La Camargue near the Mediterranean Sea. The genome of six atypical isolates was sequenced; their size is around 4.8 to 5.0 Mb, including a plasmid of 59 to 61 kb; their G+C values range from 49.1 to 49.3 mol%. Phylogenetic analyses indicated that the novel strains form a new clade of Pectobacteriaceae that branches at the basis of the group encompassing the genera Lonsdalea, Musicola, and Dickeya. Based on phenotypic, genomic and phylogenetic characteristics, we propose the creation of a new genus with the name Prodigiosinella gen. nov. Both the phenotypic and phylogenetic analyses separated the strains into two distinct subgroups, G1 and G2. The type strain LS101T (CFBP 8826T = LMG 32072T) and strain CE70 (CFBP 9054 = LMG 32867) are representative G1 and G2 members, respectively. Three genomic methods were used to analyze DNA-DNA relatedness: digital DNA-DNA hybridization (isDDH), average nucleotide identity (ANI), and genome alignment fraction (AF). They revealed a close relationship between genomes of the two groups, supporting their appurtenance to a same species for which we propose the name Prodigiosinella aquatilis sp. nov. Four strains previously designated as Serratia sp. (ATCC 39006), Brenneria "ulupoensis" (K61) or Erwinia sp. (MK01 and MK09) belong to the new genus Prodigiosinella.

Deciphering Evolutionary Trajectories of Lactate Dehydrogenases Provides New Insights into Allostery.

Lactate dehydrogenase (LDH, EC.1.1.127) is an important enzyme engaged in the anaerobic metabolism of cells, catalyzing the conversion of pyruvate to lactate and NADH to NAD+. LDH is a relevant enzyme to investigate structure-function relationships. The present work provides the missing link in our understanding of the evolution of LDHs. This allows to explain (i) the various evolutionary origins of LDHs in eukaryotic cells and their further diversification and (ii) subtle phenotypic modifications with respect to their regulation capacity. We identified a group of cyanobacterial LDHs displaying eukaryotic-like LDH sequence features. The biochemical and structural characterization of Cyanobacterium aponinum LDH, taken as representative, unexpectedly revealed that it displays homotropic and heterotropic activation, typical of an allosteric enzyme, whereas it harbors a long N-terminal extension, a structural feature considered responsible for the lack of allosteric capacity in eukaryotic LDHs. Its crystallographic structure was solved in 2 different configurations typical of the R-active and T-inactive states encountered in allosteric LDHs. Structural comparisons coupled with our evolutionary analyses helped to identify 2 amino acid positions that could have had a major role in the attenuation and extinction of the allosteric activation in eukaryotic LDHs rather than the presence of the N-terminal extension. We tested this hypothesis by site-directed mutagenesis. The resulting C. aponinum LDH mutants displayed reduced allosteric capacity mimicking those encountered in plants and human LDHs. This study provides a new evolutionary scenario of LDHs that unifies descriptions of regulatory properties with structural and mutational patterns of these important enzymes.

Biochemical and Structural Characterization of Chi-Class Glutathione Transferases: A Snapshot on the Glutathione Transferase Encoded by sll0067 Gene in the Cyanobacterium Synechocystis sp. Strain PCC 6803.

Glutathione transferases (GSTs) constitute a widespread superfamily of enzymes notably involved in detoxification processes and/or in specialized metabolism. In the cyanobacterium Synechocsytis sp. PCC 6803, SynGSTC1, a chi-class GST (GSTC), is thought to participate in the detoxification process of methylglyoxal, a toxic by-product of cellular metabolism. A comparative genomic analysis showed that GSTCs were present in all orders of cyanobacteria with the exception of the basal order Gloeobacterales. These enzymes were also detected in some marine and freshwater noncyanobacterial bacteria, probably as a result of horizontal gene transfer events. GSTCs were shorter of about 30 residues compared to most cytosolic GSTs and had a well-conserved SRAS motif in the active site (10SRAS13 in SynGSTC1). The crystal structure of SynGSTC1 in complex with glutathione adopted the canonical GST fold with a very open active site because the α4 and α5 helices were exceptionally short. A transferred multipolar electron-density analysis allowed a fine description of the solved structure. Unexpectedly, Ser10 did not have an electrostatic influence on glutathione as usually observed in serinyl-GSTs. The S10A variant was only slightly less efficient than the wild-type and molecular dynamics simulations suggested that S10 was a stabilizer of the protein backbone rather than an anchor site for glutathione.

Protein Conformational Space at the Edge of Allostery: Turning a Nonallosteric Malate Dehydrogenase into an "Allosterized" Enzyme Using Evolution-Guided Punctual Mutations.

We unveil the intimate relationship between protein dynamics and allostery by following the trajectories of model proteins in their conformational and sequence spaces. Starting from a nonallosteric hyperthermophilic malate dehydrogenase, we have tracked the role of protein dynamics in the evolution of the allosteric capacity. Based on a large phylogenetic analysis of the malate (MalDH) and lactate dehydrogenase (LDH) superfamily, we identified two amino acid positions that could have had a major role for the emergence of allostery in LDHs, which we targeted for investigation by site-directed mutagenesis. Wild-type MalDH and the single and double mutants were tested with respect to their substrate recognition profiles. The double mutant displayed a sigmoid-shaped profile typical of homotropic activation in LDH. By using molecular dynamics simulations, we showed that the mutations induce a drastic change in the protein sampling of its conformational landscape, making transiently T-like (inactive) conformers, typical of allosteric LDHs, accessible. Our data fit well with the seminal key concept linking protein dynamics and evolvability. We showed that the selection of a new phenotype can be achieved by a few key dynamics-enhancing mutations causing the enrichment of low-populated conformational substates.

A divide-and-conquer phylogenomic approach based on character supermatrices resolves early steps in the evolution of the Archaea.

BACKGROUND: The recent rise in cultivation-independent genome sequencing has provided key material to explore uncharted branches of the Tree of Life. This has been particularly spectacular concerning the Archaea, projecting them at the center stage as prominently relevant to understand early stages in evolution and the emergence of fundamental metabolisms as well as the origin of eukaryotes. Yet, resolving deep divergences remains a challenging task due to well-known tree-reconstruction artefacts and biases in extracting robust ancient phylogenetic signal, notably when analyzing data sets including the three Domains of Life. Among the various strategies aimed at mitigating these problems, divide-and-conquer approaches remain poorly explored, and have been primarily based on reconciliation among single gene trees which however notoriously lack ancient phylogenetic signal. RESULTS: We analyzed sub-sets of full supermatrices covering the whole Tree of Life with specific taxonomic sampling to robustly resolve different parts of the archaeal phylogeny in light of their current diversity. Our results strongly support the existence and early emergence of two main clades, Cluster I and Cluster II, which we name Ouranosarchaea and Gaiarchaea, and we clarify the placement of important novel archaeal lineages within these two clades. However, the monophyly and branching of the fast evolving nanosized DPANN members remains unclear and worth of further study. CONCLUSIONS: We inferred a well resolved rooted phylogeny of the Archaea that includes all recently described phyla of high taxonomic rank. This phylogeny represents a valuable reference to study the evolutionary events associated to the early steps of the diversification of the archaeal domain. Beyond the specifics of archaeal phylogeny, our results demonstrate the power of divide-and-conquer approaches to resolve deep phylogenetic relationships, which should be applied to progressively resolve the entire Tree of Life.

Phylogenetics and biochemistry elucidate the evolutionary link between l-malate and l-lactate dehydrogenases and disclose an intermediate group of sequences with mix functional properties.

The NAD(P)-dependent malate dehydrogenases (MDH) (EC 1.1.1.37) and NAD-dependent lactate dehydrogenases (LDH) (EC. 1.1.1.27) form a large superfamily that has been characterized in organisms belonging to the three Domains of Life. MDH catalyzes the reversible conversion of the oxaloacetate into malate, while LDH operates at the late stage of glycolysis by converting pyruvate into lactate. Phylogenetic studies proposed that the LDH/MDH superfamily encompasses five main groups of enzymes. Here, starting from 16,052 reference proteomes, we reinvestigated the relationships between MDH and LDH. We showed that the LDH/MDH superfamily encompasses three main families: MDH1, MDH2, and a large family encompassing MDH3, LDH, and L-2-hydroxyisocaproate dehydrogenases (HicDH) sequences. An in-depth analysis of the phylogeny of the MDH3/LDH/HicDH family and of the nature of three important amino acids, located within the catalytic site and involved in binding and substrate discrimination, revealed a large group of sequences displaying unexpected combinations of amino acids at these three critical positions. This group branched in-between canonical MDH3 and LDH sequences. The functional characterization of several enzymes from this intermediate group disclosed a mix of functional properties, indicating that the MDH3/LDH/HicDH family is much more diverse than previously thought, and blurred the frontier between MDH3 and LDH enzymes. Present-days enzymes of the intermediate group are a valuable material to study the evolutionary steps that led to functional diversity and emergence of allosteric regulation within the LDH/MDH superfamily.

Plant genera Cannabis and Humulus share the same pair of well-differentiated sex chromosomes.

We recently described, in Cannabis sativa, the oldest sex chromosome system documented so far in plants (12-28 Myr old). Based on the estimated age, we predicted that it should be shared by its sister genus Humulus, which is known also to possess XY chromosomes. Here, we used transcriptome sequencing of an F1 family of H. lupulus to identify and study the sex chromosomes in this species using the probabilistic method SEX-DETector. We identified 265 sex-linked genes in H. lupulus, which preferentially mapped to the C. sativa X chromosome. Using phylogenies of sex-linked genes, we showed that a region of the sex chromosomes had already stopped recombining in an ancestor of both species. Furthermore, as in C. sativa, Y-linked gene expression reduction is correlated to the position on the X chromosome, and highly Y degenerated genes showed dosage compensation. We report, for the first time in Angiosperms, a sex chromosome system that is shared by two different genera. Thus, recombination suppression started at least 21-25 Myr ago, and then (either gradually or step-wise) spread to a large part of the sex chromosomes (c. 70%), leading to a degenerated Y chromosome.

A Comprehensive Evolutionary Scenario of Cell Division and Associated Processes in the Firmicutes.

The cell cycle is a fundamental process that has been extensively studied in bacteria. However, many of its components and their interactions with machineries involved in other cellular processes are poorly understood. Furthermore, most knowledge relies on the study of a few models, but the real diversity of the cell division apparatus and its evolution are largely unknown. Here, we present a massive in-silico analysis of cell division and associated processes in around 1,000 genomes of the Firmicutes, a major bacterial phylum encompassing models (i.e. Bacillus subtilis, Streptococcus pneumoniae, and Staphylococcus aureus), as well as many important pathogens. We analyzed over 160 proteins by using an original approach combining phylogenetic reconciliation, phylogenetic profiles, and gene cluster survey. Our results reveal the presence of substantial differences among clades and pinpoints a number of evolutionary hotspots. In particular, the emergence of Bacilli coincides with an expansion of the gene repertoires involved in cell wall synthesis and remodeling. We also highlight major genomic rearrangements at the emergence of Streptococcaceae. We establish a functional network in Firmicutes that allows identifying new functional links inside one same process such as between FtsW (peptidoglycan polymerase) and a previously undescribed Penicilin-Binding Protein or between different processes, such as replication and cell wall synthesis. Finally, we identify new candidates involved in sporulation and cell wall synthesis. Our results provide a previously undescribed view on the diversity of the bacterial cell cycle, testable hypotheses for further experimental studies, and a methodological framework for the analysis of any other biological system.

The Molecular Determinants of Thermoadaptation: Methanococcales as a Case Study.

Previous reports have shown that environmental temperature impacts proteome evolution in Bacteria and Archaea. However, it is unknown whether thermoadaptation mainly occurs via the sequential accumulation of substitutions, massive horizontal gene transfers, or both. Measuring the real contribution of amino acid substitution to thermoadaptation is challenging, because of confounding environmental and genetic factors (e.g., pH, salinity, genomic G + C content) that also affect proteome evolution. Here, using Methanococcales, a major archaeal lineage, as a study model, we show that optimal growth temperature is the major factor affecting variations in amino acid frequencies of proteomes. By combining phylogenomic and ancestral sequence reconstruction approaches, we disclose a sequential substitutional scheme in which lysine plays a central role by fine tuning the pool of arginine, serine, threonine, glutamine, and asparagine, whose frequencies are strongly correlated with optimal growth temperature. Finally, we show that colonization to new thermal niches is not associated with high amounts of horizontal gene transfers. Altogether, although the acquisition of a few key proteins through horizontal gene transfer may have favored thermoadaptation in Methanococcales, our findings support sequential amino acid substitutions as the main factor driving thermoadaptation.

Structural and functional analysis of the Francisella lysine decarboxylase as a key actor in oxidative stress resistance.

Francisella tularensis is one of the most virulent pathogenic bacteria causing the acute human respiratory disease tularemia. While the mechanisms underlying F. tularensis pathogenesis are largely unknown, previous studies have shown that a F. novicida transposon mutant with insertions in a gene coding for a putative lysine decarboxylase was attenuated in mouse spleen, suggesting a possible role of its protein product as a virulence factor. Therefore, we set out to structurally and functionally characterize the F. novicida lysine decarboxylase, which we termed LdcF. Here, we investigate the genetic environment of ldcF as well as its evolutionary relationships with other basic AAT-fold amino acid decarboxylase superfamily members, known as key actors in bacterial adaptative stress response and polyamine biosynthesis. We determine the crystal structure of LdcF and compare it with the most thoroughly studied lysine decarboxylase, E. coli LdcI. We analyze the influence of ldcF deletion on bacterial growth under different stress conditions in dedicated growth media, as well as in infected macrophages, and demonstrate its involvement in oxidative stress resistance. Finally, our mass spectrometry-based quantitative proteomic analysis enables identification of 80 proteins with expression levels significantly affected by ldcF deletion, including several DNA repair proteins potentially involved in the diminished capacity of the F. novicida mutant to deal with oxidative stress. Taken together, we uncover an important role of LdcF in F. novicida survival in host cells through participation in oxidative stress response, thereby singling out this previously uncharacterized protein as a potential drug target.

Hepatitis B virus preS2Δ38-55 variants: A newly identified risk factor for hepatocellular carcinoma.

BACKGROUND & AIMS: Although HBV is a major cause of death in Africa, its genetic variability has been poorly documented. This study aimed to address whether HBV genotype and surface gene variants are associated with HBV-related liver disease in The Gambia. METHODS: We conducted a case-control study nested in the Prevention of Liver Fibrosis and Cancer in Africa programme. Consecutive treatment-naive patients with chronic HBV infection and detectable viral load were recruited: 211 controls with no significant liver disease and 91 cases (56 cirrhosis and 35 HCC cases). HBV genotypes and surface gene variants were determined by Sanger sequencing or next-generation sequencing (NGS) in serum DNA. Aflatoxin B1 (AFB1)-specific codon 249 TP53 mutation was determined by NGS in circulating cell-free plasma DNA. RESULTS: In phylogenetic analysis, 85% of individuals carried HBV genotype E, 14% genotype A, and 1% A/E recombinant viruses. Surface gene variants were more frequently observed in cases (43% and 57% in cirrhosis and HCC cases, respectively) than controls (25%; p <0.001), with preS2 deletions between nucleotides 38-55 (preS2Δ38-55) being the main genetic variant detected. In multivariable analysis, HBeAg seropositivity, low HBsAg levels, and HDV seropositivity were significantly associated with cirrhosis and HCC, whilst older age, higher viral load, genotype A, preS2Δ38-55, and AFB1 exposure were only associated with HCC. There was a multiplicative joint effect of preS2Δ38-55 variants with HBeAg seropositivity (odds ratio [OR] 43.1 [10.4-177.7]), high viral load >2,000 IU/ml (OR 22.7 [8.0-64.9]), HBsAg levels <10,000 IU/ml (OR 19.0 [5.5-65.3]), and AFB1 exposure (OR 29.3 [3.7-230.4]) on HCC risk. CONCLUSIONS: This study identified a hotspot for HBV preS2 deletions as a strong independent factor for HCC in The Gambia, with HBV genotypes and AFB1 exposure contributing to the high liver cancer risk. LAY SUMMARY: Although HBV-related liver disease is highly prevalent in sub-Saharan Africa, the associated virological characteristics are poorly studied. Using clinical data from African patients chronically infected with HBV, an assessment of the virological variability (genotypes and mutations) and exposure to AFB1, a toxin often contaminating food, was carried out. Our results show that HBV genotypes, the presence of a highly prevalent mutant form of HBV, and AFB1 exposure contribute to the high liver cancer risk in this population.

Dickeya poaceiphila sp. nov., a plant-pathogenic bacterium isolated from sugar cane (Saccharum officinarum).

The genus Dickeya is an important group of plant pathogens that currently comprises 10 recognized species. Although most Dickeya isolates originated from infected cultivated plants, they are also isolated from water. The genomic sequence of the Australian strain NCPPB 569T clearly established its separation from the previously characterized Dickeya species. The average nucleotide identity and digital DNA-DNA hybridization values obtained by comparing strain NCPPB 569T with strains of characterized Dickeya species were lower than 87 and 32 %, respectively, supporting the delineation of a new species. The name Dickeya poaceiphila sp. nov. is proposed for this taxon with the type strain NCPPB 569T (=CFBP 8731T). Two other strains isolated in Australia, CFBP 1537 and CFBP 2040, also belong to this species. Phenotypic and genomic comparisons enabled the identification of traits distinguishing D. poaceiphila isolates from strains of other Dickeya species.

The TERB1-TERB2-MAJIN complex of mouse meiotic telomeres dates back to the common ancestor of metazoans.

BACKGROUND: Meiosis is essential for sexual reproduction and generates genetically diverse haploid gametes from a diploid germ cell. Reduction of ploidy depends on active chromosome movements during early meiotic prophase I. Chromosome movements require telomere attachment to the nuclear envelope. This attachment is mediated by telomere adaptor proteins. Telomere adaptor proteins have to date been identified in fission yeast and mice. In the mouse, they form a complex composed of the meiotic proteins TERB1, TERB2, and MAJIN. No sequence similarity was observed between these three mouse proteins and the adaptor proteins of fission yeast, raising the question of the evolutionary history and significance of this specific protein complex. RESULT: Here, we show the TERB1, TERB2, and MAJIN proteins are found throughout the Metazoa and even in early-branching non-bilateral phyla such as Cnidaria, Placozoa and Porifera. Metazoan TERB1, TERB2, and MAJIN showed comparable domain architecture across all clades. Furthermore, the protein domains involved in the formation of the complex as well as those involved for the interaction with the telomere shelterin protein and the LINC complexes revealed high sequence similarity. Finally, gene expression in the cnidarian Hydra vulgaris provided evidence that the TERB1-TERB2-MAJIN complex is selectively expressed in the germ line. CONCLUSION: Our results indicate that the TERB1-TERB2-MAJIN complex has an ancient origin in metazoans, suggesting conservation of meiotic functions.

Draft genome sequences for three unisolated Alnus-infective Frankia Sp+ strains, AgTrS, AiOr and AvVan, the first sequenced Frankia strains able to sporulate in-planta.

Actinobacteria from genus Frankia are able to form symbiotic associations with actinorhizal plants including alders. Among them, Sp+ strains are characterized by their ability to differentiate numerous sporangia inside host plant cells (unlike "Sp-" strains unable of in-planta sporulation). Here, we report the first genome sequences of three unisolated Sp+ strains: AgTrS, AiOr and AvVan obtained from Alnus glutinosa, A. incana and A. alnobetula (previously known as viridis), respectively (with genome completeness estimated at more than 98%). They represent new Frankia species based on Average Nucleotide Identity (ANI) calculations, and the smallest Alnus-infective Frankia genomes so far sequenced (~5 Mbp), with 5,178, 6,192 and 5,751 candidate protein-encoding genes for AgTrS, AiOr and AvVan, respectively.

Structure, Function, and Evolution of the Pseudomonas aeruginosa Lysine Decarboxylase LdcA.

The only enzyme responsible for cadaverine production in the major multidrug-resistant human pathogen Pseudomonas aeruginosa is the lysine decarboxylase LdcA. This enzyme modulates the general polyamine homeostasis, promotes growth, and reduces bacterial persistence during carbenicillin treatment. Here we present a 3.7-Å resolution cryoelectron microscopy structure of LdcA. We introduce an original approach correlating phylogenetic signal with structural information and reveal possible recombination among LdcA and arginine decarboxylase subfamilies within structural domain boundaries. We show that LdcA is involved in full virulence in an insect pathogenesis model. Furthermore, unlike its enterobacterial counterparts, LdcA is regulated neither by the stringent response alarmone ppGpp nor by the AAA+ ATPase RavA. Instead, the P. aeruginosa ravA gene seems to play a defensive role. Altogether, our study identifies LdcA as an important player in P. aeruginosa physiology and virulence and as a potential drug target.

Taxonomic assignment of uncultured prokaryotes with long range PCR targeting the spectinomycin operon.

The taxonomic assignment of uncultured prokaryotes to known taxa is a major challenge in microbial systematics. This relies usually on the phylogenetic analysis of the ribosomal small subunit RNA or a few housekeeping genes. Recent works have disclosed ribosomal proteins as valuable markers for systematics and, due to the boom in complete genome sequencing, their use has become widespread. Yet, in the case of uncultured strains, for which complete genome sequences cannot be easily obtained, sequencing many markers is complicated and time consuming. Taking the advantage of the organization of ribosomal protein coding genes in large gene clusters, we amplified a 32 kb conserved region encompassing the spectinomycin (spc) operon using long range PCR from isolated and from uncultured nodular endophytic Frankia strains. The phylogenetic analysis of the 27 ribosomal protein genes contained in this region provided a robust phylogenetic tree consistent with phylogenies based on larger set of markers, indicating that this subset of ribosomal proteins contains enough phylogenetic signal to address systematic issues. This work shows that using long range PCR could break down the barrier preventing the use of ribosomal proteins as phylogenetic markers when complete genome sequences cannot be easily obtained.

Francisella tularensis: FupA mutation contributes to fluoroquinolone resistance by increasing vesicle secretion and biofilm formation.

Francisella tularensis is the causative agent in tularemia for which the high prevalence of treatment failure and relapse is a major concern. Directed-evolution experiments revealed that acquisition of fluoroquinolone (FQ) resistance was linked to factors in addition to mutations in DNA gyrase. Here, using F. tularensis live vaccine strain (LVS) as a model, we demonstrated that FupA/B (Fer-Utilization Protein) expression is linked to FQ susceptibility, and that the virulent strain F. tularensis subsp. tularensis SCHU S4 deleted for the homologous FupA protein exhibited even higher FQ resistance. In addition to an increased FQ minimal inhibitory concentration, LVSΔfupA/B displayed tolerance toward bactericidal compounds including ciprofloxacin and gentamicin. Interestingly, the FupA/B deletion was found to promote increased secretion of outer membrane vesicles (OMVs). Mass spectrometry-based quantitative proteomic characterization of vesicles from LVS and LVS∆fupA/B identified 801 proteins, including a subset of 23 proteins exhibiting differential abundance between both strains which may therefore contribute to the reduced antibiotic susceptibility of the FupA/B-deleted strain. We also demonstrated that OMVs are key structural elements of LVSΔfupA/B biofilms providing protection against FQ. These results provide a new basis for understanding and tackling antibiotic resistance and/or persistence of Francisella and other pathogenic members of the Thiotrichales class.

The phytopathogenic nature of Dickeya aquatica 174/2 and the dynamic early evolution of Dickeya pathogenicity.

Dickeya is a genus of phytopathogenic enterobacterales causing soft rot in a variety of plants (e.g. potato, chicory, maize). Among the species affiliated to this genus, Dickeya aquatica, described in 2014, remained particularly mysterious because it had no known host. Furthermore, while D. aquatica was proposed to represent a deep-branching species among Dickeya genus, its precise phylogenetic position remained elusive. Here, we report the complete genome sequence of the D. aquatica type strain 174/2. We demonstrate the affinity of D. aquatica strain 174/2 for acidic fruits such as tomato and cucumber and show that exposure of this bacterium to acidic pH induces twitching motility. An in-depth phylogenomic analysis of all available Dickeya proteomes pinpoints D. aquatica as the second deepest branching lineage within this genus and reclassifies two lineages that likely correspond to new genomospecies (gs.): Dickeya gs. poaceaephila (Dickeya sp NCPPB 569) and Dickeya gs. undicola (Dickeya sp 2B12), together with a new putative genus, tentatively named Prodigiosinella. Finally, from comparative analyses of Dickeya proteomes, we infer the complex evolutionary history of this genus, paving the way to study the adaptive patterns and processes of Dickeya to different environmental niches and hosts. In particular, we hypothesize that the lack of xylanases and xylose degradation pathways in D. aquatica could reflect adaptation to aquatic charophyte hosts which, in contrast to land plants, do not contain xyloglucans.