PanExplorer: a web-based tool for exploratory analysis and visualization of bacterial pan-genomes.

MOTIVATION: As pan-genome approaches are largely employed for bacterial comparative genomics and evolution analyses, but still difficult to be carried out by non-bioinformatician biologists, there is a need for an innovative tool facilitating the exploration of bacterial pan-genomes. RESULTS: PanExplorer is a web application providing various genomic analyses and reports, giving intuitive views that enable a better understanding of bacterial pan-genomes. As an example, we produced the pan-genome for 121 Anaplasmataceae strains (including 30 Ehrlichia, 15 Anaplasma, 68 Wolbachia). AVAILABILITY AND IMPLEMENTATION: PanExplorer is written in Perl CGI and relies on several JavaScript libraries for visualization (hotmap.js, MauveViewer, CircosJS). It is freely available at http://panexplorer.southgreen.fr. The source code has been released in a GitHub repository https://github.com/SouthGreenPlatform/PanExplorer. A documentation section is available on PanExplorer website.

Complete Genome Sequencing of Three Clade-1 Xanthomonads Reveals Genetic Determinants for a Lateral Flagellin and the Biosynthesis of Coronatine-Like Molecules in Xanthomonas.

Evolutionarily, early-branching xanthomonads, also referred to as clade-1 xanthomonads, include major plant pathogens, most of which colonize monocotyledonous plants. Seven species have been validly described, among them the two sugarcane pathogens Xanthomonas albilineans and Xanthomonas sacchari, as well as Xanthomonas translucens, which infects small-grain cereals and diverse grasses but also asparagus and pistachio trees. Single-gene sequencing and genomic approaches have indicated that this clade likely contains more, yet-undescribed species. In this study, we sequenced representative strains of three novel species using long-read sequencing technology. Xanthomonas campestris pv. phormiicola strain CFBP 8444 causes bacterial streak on New Zealand flax, another monocotyledonous plant. Xanthomonas sp. strain CFBP 8443 has been isolated from common bean, and Xanthomonas sp. strain CFBP 8445 originated from banana. Complete assemblies of the chromosomes confirmed their unique phylogenetic position within clade 1 of Xanthomonas. Genome mining revealed novel genetic features, hitherto undescribed in other members of the Xanthomonas genus. In strain CFBP 8444, we identified genes related to the synthesis of coronatine-like compounds, a phytotoxin produced by several pseudomonads, which raises interesting questions about the evolution and pathogenicity of this pathogen. Furthermore, strain CFBP 8444 was found to contain a second, atypical flagellar gene cluster in addition to the canonical flagellar gene cluster. Overall, this research represents an important step toward better understanding the evolutionary history and biology of early-branching xanthomonads.

getSequenceInfo: a suite of tools allowing to get genome sequence information from public repositories.

BACKGROUND: Biological sequences are increasing rapidly and exponentially worldwide. Nucleotide sequence databases play an important role in providing meaningful genomic information on a variety of biological organisms. RESULTS: The getSequenceInfo software tool allows to access sequence information from various public repositories (GenBank, RefSeq, and the European Nucleotide Archive), and is compatible with different operating systems (Linux, MacOS, and Microsoft Windows) in a programmatic way (command line) or as a graphical user interface. getSequenceInfo or gSeqI v1.0 should help users to get some information on queried sequences that could be useful for specific studies (e.g. the country of origin/isolation or the release date of queried sequences). Queries can be made to retrieve sequence data based on a given kingdom and species, or from a given date. This program allows the separation between chromosomes and plasmids (or other genetic elements/components) by arranging each component in a given folder. Some basic statistics are also performed by the program (such as the calculation of GC content for queried assemblies). An empirically designed nucleotide ratio is calculated using nucleotide information in order to tentatively provide a "NucleScore" for studied genome assemblies. Besides the main gSeqI tool, other additional tools have been developed to perform various tasks related to sequence analysis. CONCLUSION: The aim of this study is to democratize the use of public repositories in programmatic ways, and to facilitate sequence data analysis in a pedagogical perspective. Output results are available in FASTA, FASTQ, Excel/TSV or HTML formats. The program is freely available at: https://github.com/karubiotools/getSequenceInfo . getSequenceInfo and supplementary tools are partly available through the recently released Galaxy KaruBioNet platform ( http://calamar.univ-ag.fr/c3i/galaxy_karubionet.html ).

Genotyping by sequencing can reveal the complex mosaic genomes in gene pools resulting from reticulate evolution: a case study in diploid and polyploid citrus.

BACKGROUND AND AIMS: Reticulate evolution, coupled with reproductive features limiting further interspecific recombinations, results in admixed mosaics of large genomic fragments from the ancestral taxa. Whole-genome sequencing (WGS) data are powerful tools to decipher such complex genomes but still too costly to be used for large populations. The aim of this work was to develop an approach to infer phylogenomic structures in diploid, triploid and tetraploid individuals from sequencing data in reduced genome complexity libraries. The approach was applied to the cultivated Citrus gene pool resulting from reticulate evolution involving four ancestral taxa, C. maxima, C. medica, C. micrantha and C. reticulata. METHODS: A genotyping by sequencing library was established with the restriction enzyme ApeKI applying one base (A) selection. Diagnostic single nucleotide polymorphisms (DSNPs) for the four ancestral taxa were mined in 29 representative varieties. A generic pipeline based on a maximum likelihood analysis of the number of read data was established to infer ancestral contributions along the genome of diploid, triploid and tetraploid individuals. The pipeline was applied to 48 diploid, four triploid and one tetraploid citrus accessions. KEY RESULTS: Among 43 598 mined SNPs, we identified a set of 15 946 DSNPs covering the whole genome with a distribution similar to that of gene sequences. The set efficiently inferred the phylogenomic karyotype of the 53 analysed accessions, providing patterns for common accessions very close to that previously established using WGS data. The complex phylogenomic karyotypes of 21 cultivated citrus, including bergamot, triploid and tetraploid limes, were revealed for the first time. CONCLUSIONS: The pipeline, available online, efficiently inferred the phylogenomic structures of diploid, triploid and tetraploid citrus. It will be useful for any species whose reproductive behaviour resulted in an interspecific mosaic of large genomic fragments. It can also be used for the first generations of interspecific breeding schemes.

Development of the Automated Primer Design Workflow Uniqprimer and Diagnostic Primers for the Broad-Host-Range Plant Pathogen Dickeya dianthicola.

Uniqprimer, a software pipeline developed in Python, was deployed as a user-friendly internet tool in Rice Galaxy for comparative genome analyses to design primer sets for PCRassays capable of detecting target bacterial taxa. The pipeline was trialed with Dickeya dianthicola, a destructive broad-host-range bacterial pathogen found in most potato-growing regions. Dickeya is a highly variable genus, and some primers available to detect this genus and species exhibit common diagnostic failures. Upon uploading a selection of target and nontarget genomes, six primer sets were rapidly identified with Uniqprimer, of which two were specific and sensitive when tested with D. dianthicola. The remaining four amplified a minority of the nontarget strains tested. The two promising candidate primer sets were trialed with DNA isolated from 116 field samples from across the United States that were previously submitted for testing. D. dianthicola was detected in 41 samples, demonstrating the applicability of our detection primers and suggesting widespread occurrence of D. dianthicola in North America.

daTALbase: A Database for Genomic and Transcriptomic Data Related to TAL Effectors.

Transcription activator-like effectors (TALEs) are proteins found in the genus Xanthomonas of phytopathogenic bacteria. These proteins enter the nucleus of cells in the host plant and can induce the expression of susceptibility genes (S genes), triggering disease. TALEs bind the promoter region of S genes following a specific code, which allows the prediction of binding sites based on TALEs amino acid sequences. New candidate S genes can then be discovered by finding the intersection between genes induced in the presence of TALEs and genes containing predicted effector binding elements. By contrasting differential expression data and binding site predictions across different datasets, patterns of TALE diversification or convergence may be unveiled, but this requires the seamless integration of different genomic and transcriptomic data. With this in mind, we present daTALbase, a curated relational database that integrates TALE-related data including bacterial TALE sequences, plant promoter sequences, predicted TALE binding sites, transcriptomic data of host plants in response to TALE-harboring bacteria, and other associated data. The database can be explored to uncover new candidate S genes as well as to study variation in TALE repertories and their corresponding targets. The first version of the database here presented includes data for Oryza sp.-Xanthomonas pv. oryzae interactions. Future versions of the database will incorporate information for other pathosystems involving TALEs.

Allelic variation for broad-spectrum resistance and susceptibility to bacterial pathogens identified in a rice MAGIC population.

Quantitative trait loci (QTL) that confer broad-spectrum resistance (BSR), or resistance that is effective against multiple and diverse plant pathogens, have been elusive targets of crop breeding programmes. Multiparent advanced generation intercross (MAGIC) populations, with their diverse genetic composition and high levels of recombination, are potential resources for the identification of QTL for BSR. In this study, a rice MAGIC population was used to map QTL conferring BSR to two major rice diseases, bacterial leaf streak (BLS) and bacterial blight (BB), caused by Xanthomonas oryzae pathovars (pv.) oryzicola (Xoc) and oryzae (Xoo), respectively. Controlling these diseases is particularly important in sub-Saharan Africa, where no sources of BSR are currently available in deployed varieties. The MAGIC founders and lines were genotyped by sequencing and phenotyped in the greenhouse and field by inoculation with multiple strains of Xoc and Xoo. A combination of genomewide association studies (GWAS) and interval mapping analyses revealed 11 BSR QTL, effective against both diseases, and three pathovar-specific QTL. The most promising BSR QTL (qXO-2-1, qXO-4-1 and qXO-11-2) conferred resistance to more than nine Xoc and Xoo strains. GWAS detected 369 significant SNP markers with distinguishable phenotypic effects, allowing the identification of alleles conferring disease resistance and susceptibility. The BSR and susceptibility QTL will improve our understanding of the mechanisms of both resistance and susceptibility in the long term and will be immediately useful resources for rice breeding programmes.

Identification of a Hypervirulent Pathotype of Rice yellow mottle virus: A Threat to Genetic Resistance Deployment in West-Central Africa.

Rice yellow mottle virus (RYMV) causes high losses to rice production in Africa. Several sources of varietal high resistance are available but the emergence of virulent pathotypes that are able to overcome one or two resistance alleles can sometimes occur. Both resistance spectra and viral adaptability have to be taken into account to develop sustainable rice breeding strategies against RYMV. In this study, we extended previous resistance spectrum analyses by testing the rymv1-4 and rymv1-5 alleles that are carried by the rice accessions Tog5438 and Tog5674, respectively, against isolates that are representative of RYMV genetic and pathogenic diversity. Our study revealed a hypervirulent pathotype, named thereafter pathotype T', that is able to overcome all known sources of high resistance. This pathotype, which is spatially localized in West-Central Africa, appears to be more abundant than previously suspected. To better understand the adaptive processes of pathotype T', molecular determinants of resistance breakdown were identified via Sanger sequencing and validated through directed mutagenesis of an infectious clone. These analyses confirmed the key role of convergent nonsynonymous substitutions in the central part of the viral genome-linked protein to overcome RYMV1-mediated resistance. In addition, deep-sequencing analyses revealed that resistance breakdown does not always coincide with fixed mutations. Actually, virulence mutations that are present in a small proportion of the virus population can be sufficient for resistance breakdown. Considering the spatial distribution of RYMV strains in Africa and their ability to overcome the RYMV resistance genes and alleles, we established a resistance-breaking risk map to optimize strategies for the deployment of sustainable and resistant rice lines in Africa.

Transcriptomic and histological responses of African rice (Oryza glaberrima) to Meloidogyne graminicola provide new insights into root-knot nematode resistance in monocots.

Background and Aims: The root-knot nematode Meloidogyne graminicola is responsible for production losses in rice ( Oryza sativa ) in Asia and Latin America. The accession TOG5681 of African rice, O. glaberrima , presents improved resistance to several biotic and abiotic factors, including nematodes. The aim of this study was to assess the cytological and molecular mechanisms underlying nematode resistance in this accession. Methods: Penetration and development in M. graminicola in TOG5681 and the susceptible O. sativa genotype 'Nipponbare' were compared by microscopic observation of infected roots and histological analysis of galls. In parallel, host molecular responses to M. graminicola were assessed by root transcriptome profiling at 2, 4 and 8 d post-infection (dpi). Specific treatments with hormone inhibitors were conducted in TOG5681 to assess the impact of the jasmonic acid and salicylic acid pathways on nematode penetration and reproduction. Key Results: Penetration and development of M. graminicola juveniles were reduced in the resistant TOG5681 in comparison with the susceptible accession, with degeneration of giant cells observed in the resistant genotype from 15 dpi onwards. Transcriptome changes were observed as early as 2 dpi, with genes predicted to be involved in defence responses, phenylpropanoid and hormone pathways strongly induced in TOG5681, in contrast to 'Nipponbare'. No specific hormonal pathway could be identified as the major determinant of resistance in the rice-nematode incompatible interaction. Candidate genes proposed as involved in resistance to M. graminicola in TOG5681 were identified based on their expression pattern and quantitative trait locus (QTL) position, including chalcone synthase, isoflavone reductase, phenylalanine ammonia lyase, WRKY62 transcription factor, thionin, stripe rust resistance protein, thaumatins and ATPase3. Conclusions: This study provides a novel set of candidate genes for O. glaberrima resistance to nematodes and highlights the rice- M. graminicola pathosystem as a model to study plant-nematode incompatible interactions.

SNiPlay3: a web-based application for exploration and large scale analyses of genomic variations.

SNiPlay is a web-based tool for detection, management and analysis of genetic variants including both single nucleotide polymorphisms (SNPs) and InDels. Version 3 now extends functionalities in order to easily manage and exploit SNPs derived from next generation sequencing technologies, such as GBS (genotyping by sequencing), WGRS (whole gre-sequencing) and RNA-Seq technologies. Based on the standard VCF (variant call format) format, the application offers an intuitive interface for filtering and comparing polymorphisms using user-defined sets of individuals and then establishing a reliable genotyping data matrix for further analyses. Namely, in addition to the various scaled-up analyses allowed by the application (genomic annotation of SNP, diversity analysis, haplotype reconstruction and network, linkage disequilibrium), SNiPlay3 proposes new modules for GWAS (genome-wide association studies), population stratification, distance tree analysis and visualization of SNP density. Additionally, we developed a suite of Galaxy wrappers for each step of the SNiPlay3 process, so that the complete pipeline can also be deployed on a Galaxy instance using the Galaxy ToolShed procedure and then be computed as a Galaxy workflow. SNiPlay is accessible at http://sniplay.southgreen.fr.

The coffee genome hub: a resource for coffee genomes.

The whole genome sequence of Coffea canephora, the perennial diploid species known as Robusta, has been recently released. In the context of the C. canephora genome sequencing project and to support post-genomics efforts, we developed the Coffee Genome Hub (http://coffee-genome.org/), an integrative genome information system that allows centralized access to genomics and genetics data and analysis tools to facilitate translational and applied research in coffee. We provide the complete genome sequence of C. canephora along with gene structure, gene product information, metabolism, gene families, transcriptomics, syntenic blocks, genetic markers and genetic maps. The hub relies on generic software (e.g. GMOD tools) for easy querying, visualizing and downloading research data. It includes a Genome Browser enhanced by a Community Annotation System, enabling the improvement of automatic gene annotation through an annotation editor. In addition, the hub aims at developing interoperability among other existing South Green tools managing coffee data (phylogenomics resources, SNPs) and/or supporting data analyses with the Galaxy workflow manager.

Genetic diversity, linkage disequilibrium and power of a large grapevine (Vitis vinifera L) diversity panel newly designed for association studies.

BACKGROUND: As for many crops, new high-quality grapevine varieties requiring less pesticide and adapted to climate change are needed. In perennial species, breeding is a long process which can be speeded up by gaining knowledge about quantitative trait loci linked to agronomic traits variation. However, due to the long juvenile period of these species, establishing numerous highly recombinant populations for high resolution mapping is both costly and time-consuming. Genome wide association studies in germplasm panels is an alternative method of choice, since it allows identifying the main quantitative trait loci with high resolution by exploiting past recombination events between cultivars. Such studies require adequate panel design to represent most of the available genetic and phenotypic diversity. Assessing linkage disequilibrium extent and panel power is also needed to determine the marker density required for association studies. RESULTS: Starting from the largest grapevine collection worldwide maintained in Vassal (France), we designed a diversity panel of 279 cultivars with limited relatedness, reflecting the low structuration in three genetic pools resulting from different uses (table vs wine) and geographical origin (East vs West), and including the major founders of modern cultivars. With 20 simple sequence repeat markers and five quantitative traits, we showed that our panel adequately captured most of the genetic and phenotypic diversity existing within the entire Vassal collection. To assess linkage disequilibrium extent and panel power, we genotyped single nucleotide polymorphisms: 372 over four genomic regions and 129 distributed over the whole genome. Linkage disequilibrium, measured by correlation corrected for kinship, reached 0.2 for a physical distance between 9 and 458 Kb depending on genetic pool and genomic region, with varying size of linkage disequilibrium blocks. This panel achieved reasonable power to detect associations between traits with high broad-sense heritability (> 0.7) and causal loci with intermediate allelic frequency and strong effect (explaining > 10 % of total variance). CONCLUSIONS: Our association panel constitutes a new, highly valuable resource for genetic association studies in grapevine, and deserves dissemination to diverse field and greenhouse trials to gain more insight into the genetic control of many agronomic traits and their interaction with the environment.

Genome rearrangements derived from homoeologous recombination following allopolyploidy speciation in coffee.

Allopolyploidization is widespread and has played a major role in flowering plant diversification. Genomic changes are common consequences of allopolyploidization, but their mechanisms of occurrence and dynamics over time are still poorly understood. Coffea arabica, a recently formed allotetraploid, was chosen as a model to investigate genetic changes in allopolyploid using an approach that exploits next-generation sequencing technologies. Genes affected by putative homoeolog loss were inferred by comparing the numbers of single-nucleotide polymorphisms detected using RNA-seq in individual accessions of C. arabica, and between accessions of its two diploid progenitor species for common sequence positions. Their physical locations were investigated and clusters of genes exhibiting homoeolog loss were identified. To validate these results, genome sequencing data were generated from one accession of C. arabica and further analyzed. Genomic rearrangements involving homoeologous exchanges appear to occur in C. arabica and to be a major source of genetic diversity. At least 5% of the C. arabica genes were inferred to have undergone homoeolog loss. The detection of a large number of homoeologous exchange events (HEEs) shared by all accessions of C. arabica strongly reinforces the assumption of a single allopolyploidization event. Furthermore, HEEs were specific to one or a few accessions, suggesting that HEE accumulates gradually. Our results provide evidence for the important role of HEE in allopolyploid genome evolution.

The Greater Phenotypic Homeostasis of the Allopolyploid Coffea arabica Improved the Transcriptional Homeostasis Over that of Both Diploid Parents.

Polyploidy impacts the diversity of plant species, giving rise to novel phenotypes and leading to ecological diversification. In order to observe adaptive and evolutionary capacities of polyploids, we compared the growth, primary metabolism and transcriptomic expression level in the leaves of the newly formed allotetraploid Coffea arabica species compared with its two diploid parental species (Coffea eugenioides and Coffea canephora), exposed to four thermal regimes (TRs; 18-14, 23-19, 28-24 and 33-29°C). The growth rate of the allopolyploid C. arabica was similar to that of C. canephora under the hottest TR and that of C. eugenioides under the coldest TR. For metabolite contents measured at the hottest TR, the allopolyploid showed similar behavior to C. canephora, the parent which tolerates higher growth temperatures in the natural environment. However, at the coldest TR, the allopolyploid displayed higher sucrose, raffinose and ABA contents than those of its two parents and similar linolenic acid leaf composition and Chl content to those of C. eugenioides. At the gene expression level, few differences between the allopolyploid and its parents were observed for studied genes linked to photosynthesis, respiration and the circadian clock, whereas genes linked to redox activity showed a greater capacity of the allopolyploid for homeostasis. Finally, we found that the overall transcriptional response to TRs of the allopolyploid was more homeostatic compared with its parents. This better transcriptional homeostasis of the allopolyploid C. arabica afforded a greater phenotypic homeostasis when faced with environments that are unsuited to the diploid parental species.

BAC-end sequences analysis provides first insights into coffee (Coffea canephora P.) genome composition and evolution.

Coffee is one of the world's most important agricultural commodities. Coffee belongs to the Rubiaceae family in the euasterid I clade of dicotyledonous plants, to which the Solanaceae family also belongs. Two bacterial artificial chromosome (BAC) libraries of a homozygous doubled haploid plant of Coffea canephora were constructed using two enzymes, HindIII and BstYI. A total of 134,827 high quality BAC-end sequences (BESs) were generated from the 73,728 clones of the two libraries, and 131,412 BESs were conserved for further analysis after elimination of chloroplast and mitochondrial sequences. This corresponded to almost 13 % of the estimated size of the C. canephora genome. 6.7 % of BESs contained simple sequence repeats, the most abundant (47.8 %) being mononucleotide motifs. These sequences allow the development of numerous useful marker sites. Potential transposable elements (TEs) represented 11.9 % of the full length BESs. A difference was observed between the BstYI and HindIII libraries (14.9 vs. 8.8 %). Analysis of BESs against known coding sequences of TEs indicated that 11.9 % of the genome corresponded to known repeat sequences, like for other flowering plants. The number of genes in the coffee genome was estimated at 41,973 which is probably overestimated. Comparative genome mapping revealed that microsynteny was higher between coffee and grapevine than between coffee and tomato or Arabidopsis. BESs constitute valuable resources for the first genome wide survey of coffee and provide new insights into the composition and evolution of the coffee genome.

Contribution of subgenomes to the transcriptome and their intertwined regulation in the allopolyploid Coffea arabica grown at contrasted temperatures.

Polyploidy has occurred throughout the evolutionary history of plants and led to diversification and plant ecological adaptation. Functional plasticity of duplicate genes is believed to play a major role in the environmental adaptation of polyploids. In this context, we characterized genome-wide homoeologous gene expression in Coffea arabica, a recent allopolyploid combining two subgenomes that derive from two closely related diploid species, and investigated its variation in response to changing environment. The transcriptome of leaves of C. arabica cultivated at different growing temperatures suitable for one or the other parental species was examined using RNA-sequencing. The relative contribution of homoeologs to gene expression was estimated for 9959 and 10,628 genes in warm and cold conditions, respectively. Whatever the growing conditions, 65% of the genes showed equivalent levels of homoeologous gene expression. In 92% of the genes, relative homoeologous gene expression varied < 10% between growing temperatures. The subgenome contributions to the transcriptome appeared to be only marginally altered by the different conditions (involving intertwined regulations of homeologs) suggesting that C. arabica's ability to tolerate a broader range of growing temperatures than its diploid parents does not result from differential use of homoeologs.

Meiotic Behaviors of Allotetraploid Citrus Drive the Interspecific Recombination Landscape, the Genetic Structures, and Traits Inheritance in Tetrazyg Progenies Aiming to Select New Rootstocks.

Sexual breeding at the tetraploid level is a promising strategy for rootstock breeding in citrus. Due to the interspecific origin of most of the conventional diploid citrus rootstocks that produced the tetraploid germplasm, the optimization of this strategy requires better knowledge of the meiotic behavior of the tetraploid parents. This work used Genotyping By Sequencing (GBS) data from 103 tetraploid hybrids to study the meiotic behavior and generate a high-density recombination landscape for their tetraploid intergenic Swingle citrumelo and interspecific Volkamer lemon progenitors. A genetic association study was performed with root architecture traits. For citrumelo, high preferential chromosome pairing was revealed and led to an intermediate inheritance with a disomic tendency. Meiosis in Volkamer lemon was more complex than that of citrumelo, with mixed segregation patterns from disomy to tetrasomy. The preferential pairing resulted in low interspecific recombination levels and high interspecific heterozygosity transmission by the diploid gametes. This meiotic behavior affected the efficiency of Quantitative Trait Loci (QTL) detection. Nevertheless, it enabled a high transmission of disease and pest resistance candidate genes from P. trifoliata that are heterozygous in the citrumelo progenitor. The tetrazyg strategy, using doubled diploids of interspecific origin as parents, appears to be efficient in transferring the dominant traits selected at the parental level to the tetraploid progenies.

Clustered Regularly Interspaced Short Palindromic Repeats in Xanthomonas citri-Witnesses to a Global Expansion of a Bacterial Pathogen over Time.

Xanthomonas citri pv. citri, a Gram-negative bacterium, is the causal agent of citrus canker, a significant threat to citrus production. Understanding of global expansion of the pathogen and monitoring introduction into new regions are of interest for integrated disease management at the local and global level. Genetic diversity can be assessed using genomic approaches or information from partial gene sequences, satellite markers or clustered regularly interspaced short palindromic repeats (CRISPR). Here, we compared CRISPR loci from 355 strains of X. citri pv. citri, including a sample from ancient DNA, and generated the genealogy of the spoligotypes, i.e., the absence/presence patterns of CRISPR spacers. We identified 26 novel spoligotypes and constructed their likely evolutionary trajectory based on the whole-genome information. Moreover, we analyzed ~30 additional pathovars of X. citri and found that the oldest part of the CRISPR array was present in the ancestor of several pathovars of X. citri. This work presents a framework for further analyses of CRISPR loci and allows drawing conclusions about the global spread of the citrus canker pathogen, as exemplified by two introductions in West Africa.

Naegleria genus pangenome reveals new structural and functional insights into the versatility of these free-living amoebae.

Introduction: Free-living amoebae of the Naegleria genus belong to the major protist clade Heterolobosea and are ubiquitously distributed in soil and freshwater habitats. Of the 47 Naegleria species described, N. fowleri is the only one being pathogenic to humans, causing a rare but fulminant primary amoebic meningoencephalitis. Some Naegleria genome sequences are publicly available, but the genetic basis for Naegleria diversity and ability to thrive in diverse environments (including human brain) remains unclear. Methods: Herein, we constructed a high-quality Naegleria genus pangenome to obtain a comprehensive catalog of genes encoded by these amoebae. For this, we first sequenced, assembled, and annotated six new Naegleria genomes. Results and Discussion: Genome architecture analyses revealed that Naegleria may use genome plasticity features such as ploidy/aneuploidy to modulate their behavior in different environments. When comparing 14 near-to-complete genome sequences, our results estimated the theoretical Naegleria pangenome as a closed genome, with 13,943 genes, including 3,563 core and 10,380 accessory genes. The functional annotations revealed that a large fraction of Naegleria genes show significant sequence similarity with those already described in other kingdoms, namely Animalia and Plantae. Comparative analyses highlighted a remarkable genomic heterogeneity, even for closely related strains and demonstrate that Naegleria harbors extensive genome variability, reflected in different metabolic repertoires. If Naegleria core genome was enriched in conserved genes essential for metabolic, regulatory and survival processes, the accessory genome revealed the presence of genes involved in stress response, macromolecule modifications, cell signaling and immune response. Commonly reported N. fowleri virulence-associated genes were present in both core and accessory genomes, suggesting that N. fowleri's ability to infect human brain could be related to its unique species-specific genes (mostly of unknown function) and/or to differential gene expression. The construction of Naegleria first pangenome allowed us to move away from a single reference genome (that does not necessarily represent each species as a whole) and to identify essential and dispensable genes in Naegleria evolution, diversity and biology, paving the way for further genomic and post-genomic studies.

Emergence of a Novel Lineage and Wide Spread of a blaCTX-M-15/IncHI2/ST1 Plasmid among Nosocomial Enterobacter in Guadeloupe.

Between April 2018 and August 2019, a total of 135 strains of Enterobacter cloacae complex (ECC) were randomly collected at the University Hospital Center of Guadeloupe to investigate the structure and diversity of the local bacterial population. These nosocomial isolates were initially identified genetically by the hsp60 typing method, which revealed the clinical relevance of E. xiangfangensis (n = 69). Overall, 57/94 of the third cephalosporin-resistant strains were characterized as extended-spectrum-ß-lactamase (ESBL) producers, and their whole-genome was sequenced using Illumina technology to determine the clonal relatedness and diffusion of resistance genes. We found limited genetic diversity among sequence types (STs). ST114 (n = 13), ST1503 (n = 9), ST53 (n = 5) and ST113 (n = 4), which belong to three different Enterobacter species, were the most prevalent among the 57 ESBL producers. The blaCTXM-15 gene was the most prevalent ESBL determinant (56/57) and was in most cases associated with IncHI2/ST1 plasmid replicon carriage (36/57). To fully characterize this predominant blaCTXM-15/IncHI2/ST1 plasmid, four isolates from different lineages were also sequenced using Oxford Nanopore sequencing technology to generate long-reads. Hybrid sequence analyses confirmed the circulation of a well-conserved plasmid among ECC members. In addition, the novel ST1503 and its associated species (ECC taxon 4) were analyzed, in view of its high prevalence in nosocomial infections. These genetic observations confirmed the overall incidence of nosocomial ESBL Enterobacteriaceae infections acquired in this hospital during the study period, which was clearly higher in Guadeloupe (1.59/1000 hospitalization days) than in mainland France (0.52/1,000 hospitalization days). This project revealed issues and future challenges for the management and surveillance of nosocomial and multidrug-resistant Enterobacter in the Caribbean.