EuGene-PP: a next-generation automated annotation pipeline for prokaryotic genomes

Comparative genomics of Salmonella enterica serovars Paratyphi A, Typhi and Typhimurium reveals distinct profiles of their pangenome, mobile genetic elements, antimicrobial resistance and defense systems repertoire

Salmonella enterica (S. enterica) is a highly ubiquitous and diverse animal and human pathogen. Distinct S. enterica serovars may present varying host-specificity and cause different diseases. While the human-restricted serovars S. Typhi (STY) and S. Paratyphi A (SPA) cause in humans a systemic life-threatening enteric fever, the host-generalist serovar, S. Typhimurium (STM) causes in immunocompetent individuals a self-limited gastroenteritis. Here, we have performed whole-genome sequencing and hybrid assembly of new SPA and STY typhoidal strains and took a comparative genomics approach to examine their phylogeny, pangenome structure and accessory genome content in comparison to the reference non-typhoidal serovar, STM. Our results identified previously uncharacterized lineages of SPA and refined the presence and distribution of core pseudogenes in typhoidal serovars. Pangenome analysis showed that while these serovars have a relatively similar core-genome size, the accessory genome of STM is more than four times larger than those of typhoidal Salmonellae and that STY and SPA display a more closed pangenome than STM. Unexpectedly, we demonstrate that STY and SPA present distinct differences in their pangenome composition, with a noticeable lower number of prophages, conjugative elements and antimicrobial genes per genome in SPA vs. STY. These results suggest that although SPA and STY are closely related at the DNA level, share a similar lifestyle and cause a symptomatic-indistinguishable disease, their genomic evolution and accessory genomes are markedly different. Moreover, these results may provide genomic explanation to phenotypic and epidemiological differences in antimicrobial resistance profiles associated with these serovars globally.

MEGAnnotator2: a pipeline for the assembly and annotation of microbial genomes

The reconstruction of microbial genome sequences by bioinformatic pipelines and the consequent functional annotation of their genes' repertoire are fundamental activities aiming at unveiling their biological mechanisms, such as metabolism, virulence factors, and antimicrobial resistances. Here, we describe the development of the MEGAnnotator2 pipeline able to manage all next-generation sequencing methodologies producing short- and long-read DNA sequences. Starting from raw sequencing data, the updated pipeline can manage multiple analyses leading to the assembly of high-quality genome sequences and the functional classification of their genetic repertoire, providing the user with a useful report constituting features and statistics related to the microbial genome. The updated pipeline is fully automated from the installation to the delivery of the output, thus requiring minimal bioinformatics knowledge to be executed.

Intracellular Salmonella Paratyphi A is motile and differs in the expression of flagella-chemotaxis, SPI-1 and carbon utilization pathways in comparison to intracellular S. Typhimurium

Although Salmonella Typhimurium (STM) and Salmonella Paratyphi A (SPA) belong to the same phylogenetic species, share large portions of their genome and express many common virulence factors, they differ vastly in their host specificity, the immune response they elicit, and the clinical manifestations they cause. In this work, we compared their intracellular transcriptomic architecture and cellular phenotypes during human epithelial cell infection. While transcription induction of many metal transport systems, purines, biotin, PhoPQ and SPI-2 regulons was similar in both intracellular SPA and STM, we identified 234 differentially expressed genes that showed distinct expression patterns in intracellular SPA vs. STM. Surprisingly, clear expression differences were found in SPI-1, motility and chemotaxis, and carbon (mainly citrate, galactonate and ethanolamine) utilization pathways, indicating that these pathways are regulated differently during their intracellular phase. Concurring, on the cellular level, we show that while the majority of STM are non-motile and reside within Salmonella-Containing Vacuoles (SCV), a significant proportion of intracellular SPA cells are motile and compartmentalized in the cytosol. Moreover, we found that the elevated expression of SPI-1 and motility genes by intracellular SPA results in increased invasiveness of SPA, following exit from host cells. These findings demonstrate unexpected flagellum-dependent intracellular motility of a typhoidal Salmonella serovar and intriguing differences in intracellular localization between typhoidal and non-typhoidal salmonellae. We propose that these differences facilitate new cycles of host cell infection by SPA and may contribute to the ability of SPA to disseminate beyond the intestinal lamina propria of the human host during enteric fever.

Salmonella enterica is a ubiquitous, facultative intracellular animal and human pathogen. Although non-typhoidal Salmonella (NTS) and typhoidal Salmonella serovars belong to the same phylogenetic species and share many virulence factors, the disease they cause in humans is very different. While the underlying mechanisms for these differences are not fully understood, one possible reason expected to contribute to their different pathogenicity is a distinct expression pattern of genes involved in host-pathogen interactions. Here, we compared the global gene expression and intracellular phenotypes, during human epithelial cell infection of S. Paratyphi A (SPA) and S. Typhimurium (STM), as prototypical serovars of typhoidal and NTS, respectively. Interestingly, we identified different expression patterns in key virulence and metabolic pathways, cytosolic motility and increased reinvasion of SPA, following exit from infected cells. We hypothesize that these differences contribute to the invasive and systemic disease developed following SPA infection in humans.

Cantaloupe melon genome reveals 3D chromatin features and structural relationship with the ancestral cucurbitaceae karyotype

Cucumis melo displays a large diversity of horticultural groups with cantaloupe melon the most cultivated type. Using a combination of single-molecule sequencing, 10X Genomics link-reads, high-density optical and genetic maps, and chromosome conformation capture (Hi-C), we assembled a chromosome scale C. melo var. cantalupensis Charentais mono genome. Integration of RNA-seq, MeDip-seq, ChIP-seq, and Hi-C data revealed a widespread compartmentalization of the melon genome, segregating constitutive heterochromatin and euchromatin. Genome-wide comparative and evolutionary analysis between melon botanical groups identified Charentais mono genome increasingly more divergent from Harukei-3 (reticulatus), Payzawat (inodorus), and HS (ssp. agrestis) genomes. To assess the paleohistory of the Cucurbitaceae, we reconstructed the ancestral Cucurbitaceae karyotype and compared it to sequenced cucurbit genomes. In contrast to other species that experienced massive chromosome shuffling, melon has retained the ancestral genome structure. We provide comprehensive genomic resources and new insights in the diversity of melon horticultural groups and evolution of cucurbits.

•

We provide a chromosome scale C. melo var. cantalupensis Charentais mono genome

•

Epigenomic analysis revealed a widespread compartmentalization of the melon genome

•

We reconstructed the ancestral Cucurbitaceae karyotype

•

Melon has retained the ancestral Cucurbitaceae genome structure

Pangenome of white lupin provides insights into the diversity of the species

White lupin is an old crop with renewed interest due to its seed high protein content and high nutritional value. Despite a long domestication history in the Mediterranean basin, modern breeding efforts have been fairly scarce. Recent sequencing of its genome has provided tools for further description of genetic resources but detailed characterization of genomic diversity is still missing. Here, we report the genome sequencing of 39 accessions that were used to establish a white lupin pangenome. We defined 32 068 core genes that are present in all individuals and 14 822 that are absent in some and may represent a gene pool for breeding for improved productivity, grain quality, and stress adaptation. We used this new pangenome resource to identify candidate genes for alkaloid synthesis, a key grain quality trait. The white lupin pangenome provides a novel genetic resource to better understand how domestication has shaped the genomic variability within this crop. Thus, this pangenome resource is an important step towards the effective and efficient genetic improvement of white lupin to help meet the rapidly growing demand for plant protein sources for human and animal consumption.

Prokaryotic Genome Annotation

In the last decade, the high-throughput and relatively low cost of short-read sequencing technologies have revolutionized prokaryotic genomics. This has led to an exponential increase in the number of bacterial and archaeal genome sequences available, as well as corresponding increase of genome assembly and annotation tools developed. Together, these hardware and software technologies have given scientists unprecedented options to study their chosen microbial systems without the need for large teams of bioinformaticists or supercomputing facilities. While these analysis tools largely fall into only a few categories, each may have different requirements, caveats and file formats, and some may be rarely updated or even abandoned. And so, despite the apparent ease in sequencing and analyzing a prokaryotic genome, it is no wonder that the budding genomicist may quickly find oneself overwhelmed. Here, we aim to provide the reader with an overview of genome annotation and its most important considerations, as well as an easy-to-follow protocol to get started with annotating a prokaryotic genome.

Genetic Variation in Host-Specific Competitiveness of the Symbiont Rhizobium leguminosarum Symbiovar viciae

Legumes of the Fabeae tribe form nitrogen-fixing root nodules resulting from symbiotic interaction with the soil bacteria Rhizobium leguminosarum symbiovar viciae (Rlv). These bacteria are all potential symbionts of the Fabeae hosts but display variable partner choice when co-inoculated in mixture. Because partner choice and symbiotic nitrogen fixation mostly behave as genetically independent traits, the efficiency of symbiosis is often suboptimal when Fabeae legumes are exposed to natural Rlv populations present in soil. A core collection of 32 Rlv bacteria was constituted based on the genomic comparison of a collection of 121 genome sequences, representative of known worldwide diversity of Rlv. A variable part of the nodD gene sequence was used as a DNA barcode to discriminate and quantify each of the 32 bacteria in mixture. This core collection was co-inoculated on a panel of nine genetically diverse Pisum sativum, Vicia faba, and Lens culinaris genotypes. We estimated the relative Early Partner Choice (EPC) of the bacteria with the Fabeae hosts by DNA metabarcoding on the nodulated root systems. Comparative genomic analyses within the bacterial core collection identified molecular markers associated with host-dependent symbiotic partner choice. The results revealed emergent properties of rhizobial populations. They pave the way to identify genes related to important symbiotic traits operating at this level.

Bradyrhizobium diazoefficiens USDA110 Nodulation of Aeschynomene afraspera Is Associated with Atypical Terminal Bacteroid Differentiation and Suboptimal Symbiotic Efficiency

Legume-rhizobium symbiosis is a major ecological process in the nitrogen cycle, responsible for the main input of fixed nitrogen into the biosphere. The efficiency of this symbiosis relies on the coevolution of the partners.

Legume plants can form root organs called nodules where they house intracellular symbiotic rhizobium bacteria. Within nodule cells, rhizobia differentiate into bacteroids, which fix nitrogen for the benefit of the plant. Depending on the combination of host plants and rhizobial strains, the output of rhizobium-legume interactions varies from nonfixing associations to symbioses that are highly beneficial for the plant. Bradyrhizobium diazoefficiens USDA110 was isolated as a soybean symbiont, but it can also establish a functional symbiotic interaction with Aeschynomene afraspera. In contrast to soybean, A. afraspera triggers terminal bacteroid differentiation, a process involving bacterial cell elongation, polyploidy, and increased membrane permeability, leading to a loss of bacterial viability while plants increase their symbiotic benefit. A combination of plant metabolomics, bacterial proteomics, and transcriptomics along with cytological analyses were used to study the physiology of USDA110 bacteroids in these two host plants. We show that USDA110 establishes a poorly efficient symbiosis with A. afraspera despite the full activation of the bacterial symbiotic program. We found molecular signatures of high levels of stress in A. afraspera bacteroids, whereas those of terminal bacteroid differentiation were only partially activated. Finally, we show that in A. afraspera, USDA110 bacteroids undergo atypical terminal differentiation hallmarked by the disconnection of the canonical features of this process. This study pinpoints how a rhizobium strain can adapt its physiology to a new host and cope with terminal differentiation when it did not coevolve with such a host.

IMPORTANCE Legume-rhizobium symbiosis is a major ecological process in the nitrogen cycle, responsible for the main input of fixed nitrogen into the biosphere. The efficiency of this symbiosis relies on the coevolution of the partners. Some, but not all, legume plants optimize their return on investment in the symbiosis by imposing on their microsymbionts a terminal differentiation program that increases their symbiotic efficiency but imposes a high level of stress and drastically reduces their viability. We combined multi-omics with physiological analyses to show that the symbiotic couple formed by Bradyrhizobium diazoefficiens USDA110 and Aeschynomene afraspera, in which the host and symbiont did not evolve together, is functional but displays a low symbiotic efficiency associated with a disconnection of terminal bacteroid differentiation features.

Host-specific competitiveness to form nodules in Rhizobium leguminosarum symbiovar viciae

Fabeae legumes such as pea and faba bean form symbiotic nodules with a large diversity of soil Rhizobium leguminosarum symbiovar viciae (Rlv) bacteria. However, bacteria competitive to form root nodules (CFN) are generally not the most efficient to fix dinitrogen, resulting in a decrease in legume crop yields. Here, we investigate differential selection by host plants on the diversity of Rlv. A large collection of Rlv was collected by nodule trapping with pea and faba bean from soils at five European sites. Representative genomes were sequenced. In parallel, diversity and abundance of Rlv were estimated directly in these soils using metabarcoding. The CFN of isolates was measured with both legume hosts. Pea/faba bean CFN were associated to Rlv genomic regions. Variations of bacterial pea and/or faba bean CFN explained the differential abundance of Rlv genotypes in pea and faba bean nodules. No evidence was found for genetic association between CFN and variations in the core genome, but variations in specific regions of the nod locus, as well as in other plasmid loci, were associated with differences in CFN. These findings shed light on the genetic control of CFN in Rlv and emphasise the importance of host plants in controlling Rhizobium diversity.

Zucchini Vein Clearing Disease Is Caused by Several Lineages Within Pseudomonas syringae Species Complex

Zucchini (Cucurbita pepo) is worldwide affected by Pseudomonas syringae, inducing vein clearing, stunting, and necroses during plantlet development. A collection of 58 P. syringae strains isolated from diseased zucchini plantlets was characterized by multilocus sequence analysis (MLSA). A subset of 23 strains responsible for vein clearing of zucchini (VCZ) was evaluated for pathogenicity on zucchini, and their genomes were sequenced. The host range of six VCZ strains was evaluated on 11 cucurbit species. Most VCZ strains belong to clades 2a and 2b-a within phylogroup 2 of P. syringae species complex and are closely related to other strains previously isolated from cucurbits. Genome analyses revealed diversity among VCZ strains within each clade. One main cluster, once referred to by the invalid pathovar name (peponis), gathers VCZ strains presenting a narrow host range including zucchini and squashes. Other VCZ strains present a large host range including zucchini, squashes, cucumber, melons, and in some cases watermelon. The VCZ strain pathogenic features are strongly associated with type III effector repertoires. The presence of avrRpt2 and absence of hopZ5 are associated with a narrow host range, whereas the presence of hopZ5 and absence of avrRpt2 are most generally associated with a large host range. To better detect the different clusters identified with whole genome sequence and pathogenicity analyses, we used a specific-k-mers approach to refine the MLSA scheme. Using this novel MLSA scheme to type P. syringae isolates from diseased cucurbits would give insight into distribution of worldwide strains and origin of epidemics.

Xanthomonas citri pv. viticola Affecting Grapevine in Brazil: Emergence of a Successful Monomorphic Pathogen

The pathovar viticola of Xanthomonas citri causes bacterial canker of grapevine. This disease was first recorded in India in 1972, and later in Brazil in 1998, where its distribution is currently restricted to the northeastern region. A multilocus sequence analysis (MLSA) based on seven housekeeping genes and a multilocus variable number of tandem repeat analysis (MLVA) with eight loci were performed in order to assess the genetic relatedness among strains from India and Brazil. Strains isolated in India from three related pathovars affecting Vitaceae species and pathogenic strains isolated from Amaranthus sp. found in bacterial canker-infected vineyards in Brazil were also included. MLSA revealed lack of diversity in all seven genes and grouped grapevine and Amaranthus strains in a monophyletic group in X. citri. The VNTR (variable number of tandem repeat) typing scheme conducted on 107 strains detected 101 haplotypes. The total number of alleles per locus ranged from 5 to 12. A minimum spanning tree (MST) showed that Brazilian strains were clearly separated from Indian strains, which showed unique alleles at three loci. The two strains isolated from symptomatic Amaranthus sp. presented unique alleles at two loci. STRUCTURE analyses revealed three groups congruent with MST and a fourth group with strains from India and Brazil. Admixture among populations were observed in all groups. MST, STRUCTURE and e-BURST analyses showed that the strains collected in 1998 belong to two distinct groups, with predicted founder genotypes from two different vineyards in the same region. This suggest that one introduction of grape planting materials contaminated with genetically distinct strains took place, which was followed by pathogen adaptation. Genome sequencing of one Brazilian strain confirmed typical attributes of pathogenic xanthomonads and allowed the design of a complementary VNTR typing scheme dedicated to X. citri pv. viticola that will allow further epidemiological survey of this genetically monomorphic pathovar.

Identification of genetic relationships and subspecies signatures in Xylella fastidiosa

BACKGROUND

The phytopathogenic bacterium Xylella fastidiosa was thought to be restricted to the Americas where it infects and kills numerous hosts. Its detection worldwide has been blooming since 2013 in Europe and Asia. Genetically diverse, this species is divided into six subspecies but genetic traits governing this classification are poorly understood.

RESULTS

SkIf (Specific k-mers Identification) was designed and exploited for comparative genomics on a dataset of 46 X. fastidiosa genomes, including seven newly sequenced individuals. It was helpful to quickly check the synonymy between strains from different collections. SkIf identified specific SNPs within 16S rRNA sequences that can be employed for predicting the distribution of Xylella through data mining. Applied to inter- and intra-subspecies analyses, it identified specific k-mers in genes affiliated to differential gene ontologies. Chemotaxis-related genes more prevalently possess specific k-mers in genomes from subspecies fastidiosa, morus and sandyi taken as a whole group. In the subspecies pauca increased abundance of specific k-mers was found in genes associated with the bacterial cell wall/envelope/plasma membrane. Most often, the k-mer specificity occurred in core genes with non-synonymous SNPs in their sequences in genomes of the other subspecies, suggesting putative impact in the protein functions. The presence of two integrative and conjugative elements (ICEs) was identified, one chromosomic and an entire plasmid in a single strain of X. fastidiosa subsp. pauca. Finally, a revised taxonomy of X. fastidiosa into three major clades defined by the subspecies pauca (clade I), multiplex (clade II) and the combination of fastidiosa, morus and sandyi (clade III) was strongly supported by k-mers specifically associated with these subspecies.

CONCLUSIONS

SkIf is a robust and rapid software, freely available, that can be dedicated to the comparison of sequence datasets and is applicable to any field of research. Applied to X. fastidiosa, an emerging pathogen in Europe, it provided an important resource to mine for identifying genetic markers of subspecies to optimize the strategies attempted to limit the pathogen dissemination in novel areas.

EuGene: An Automated Integrative Gene Finder for Eukaryotes and Prokaryotes

EuGene is an integrative gene finder applicable to both prokaryotic and eukaryotic genomes. EuGene annotated its first genome in 1999. Starting from genomic DNA sequences representing a complete genome, EuGene is able to predict the major transcript units in the genome from a variety of sources of information: statistical information, similarities with known transcripts and proteins, but also any GFF3 structured information supporting the presence or absence of specific types of elements. EuGene has been used to find genes in the plants Arabidopsis thaliana, Medicago truncatula, and Theobroma cacao; tomato, sunflower, and Rosa genomes; and in the nematode Meloidogyne incognita genome, among many others. The large fraction of plant in this list probably influenced EuGene development, especially in its capacities to withstand a genome with a large number of repeated regions and transposable elements.Depending on the sources of information used for prediction, EuGene can be considered as purely ab initio, purely similarity based, or hybrid. With the general availability of NGS-transcribed sequence data in genome projects, EuGene adopts a default hybrid behavior that strongly relies on similarity information. Initially targeted at eukaryotic genomes, EuGene has also been extended to offer integrative gene prediction for bacteria, allowing for richer and robust predictions than either purely statistical or homology-based prokaryotic gene finders.This text has been written as a practical guide that will give you the capacity to train and execute EuGene on your favorite eukaryotic genome. As the prokaryotic case is simpler and has already been described, only the main differences with the eukaryotic version were reported.

ANNOgesic: a Swiss army knife for the RNA-seq based annotation of bacterial/archaeal genomes

To understand the gene regulation of an organism of interest, a comprehensive genome annotation is essential. While some features, such as coding sequences, can be computationally predicted with high accuracy based purely on the genomic sequence, others, such as promoter elements or noncoding RNAs, are harder to detect. RNA sequencing (RNA-seq) has proven to be an efficient method to identify these genomic features and to improve genome annotations. However, processing and integrating RNA-seq data in order to generate high-resolution annotations is challenging, time consuming, and requires numerous steps. We have constructed a powerful and modular tool called ANNOgesic that provides the required analyses and simplifies RNA-seq-based bacterial and archaeal genome annotation. It can integrate data from conventional RNA-seq and differential RNA-seq and predicts and annotates numerous features, including small noncoding RNAs, with high precision. The software is available under an open source license (ISCL) at https://pypi.org/project/ANNOgesic/.

Horizontal gene transfer plays a major role in the pathological convergence of Xanthomonas lineages on common bean

BACKGROUND

Host specialization is a hallmark of numerous plant pathogens including bacteria, fungi, oomycetes and viruses. Yet, the molecular and evolutionary bases of host specificity are poorly understood. In some cases, pathological convergence is observed for individuals belonging to distant phylogenetic clades. This is the case for Xanthomonas strains responsible for common bacterial blight of bean, spread across four genetic lineages. All the strains from these four lineages converged for pathogenicity on common bean, implying possible gene convergences and/or sharing of a common arsenal of genes conferring the ability to infect common bean.

RESULTS

To search for genes involved in common bean specificity, we used a combination of whole-genome analyses without a priori, including a genome scan based on k-mer search. Analysis of 72 genomes from a collection of Xanthomonas pathovars unveiled 115 genes bearing DNA sequences specific to strains responsible for common bacterial blight, including 20 genes located on a plasmid. Of these 115 genes, 88 were involved in successive events of horizontal gene transfers among the four genetic lineages, and 44 contained nonsynonymous polymorphisms unique to the causal agents of common bacterial blight.

CONCLUSIONS

Our study revealed that host specificity of common bacterial blight agents is associated with a combination of horizontal transfers of genes, and highlights the role of plasmids in these horizontal transfers.

Xanthomonas prunicola sp. nov., a novel pathogen that affects nectarine (Prunus persica var. nectarina) trees

Three isolates obtained from symptomatic nectarine trees (Prunus persica var. nectarina) cultivated in Murcia, Spain, which showed yellow and mucoid colonies similar to Xanthomonas arboricola pv. pruni, were negative after serological and real-time PCR analyses for this pathogen. For that reason, these isolates were characterized following a polyphasic approach that included both phenotypic and genomic methods. By sequence analysis of the 16S rRNA gene, these novel strains were identified as members of the genus Xanthomonas, and by multilocus sequence analysis (MLSA) they were clustered together in a distinct group that showed similarity values below 95 % with the rest of the species of this genus. Whole-genome comparisons of the average nucleotide identity (ANI) of genomes of the strains showed less than 91 % average nucleotide identity with all other species of the genus Xanthomonas. Additionally, phenotypic characterization based on API 20 NE, API 50 CH and BIOLOG tests differentiated the strains from the species of the genus Xanthomonas described previously. Moreover, the three strains were confirmed to be pathogenic on peach (Prunus persica), causing necrotic lesions on leaves. On the basis of these results, the novel strains represent a novel species of the genus Xanthomonas, for which the name Xanthomonas prunicola is proposed. The type strain is CFBP 8353 (=CECT 9404=IVIA 3287.1).

Ancestral acquisitions, gene flow and multiple evolutionary trajectories of the type three secretion system and effectors in Xanthomonas plant pathogens

Deciphering the evolutionary history and transmission patterns of virulence determinants is necessary to understand the emergence of novel pathogens. The main virulence determinant of most pathogenic proteobacteria is the type three secretion system (T3SS). The Xanthomonas genus includes bacteria responsible for numerous epidemics in agroecosystems worldwide and represents a major threat to plant health. The main virulence factor of Xanthomonas is the Hrp2 family T3SS; however, this system is not conserved in all strains and it has not been previously determined whether the distribution of T3SS in this bacterial genus has resulted from losses or independent acquisitions. Based on comparative genomics of 82 genome sequences representing the diversity of the genus, we have inferred three ancestral acquisitions of the Hrp2 cluster during Xanthomonas evolution followed by subsequent losses in some commensal strains and re‐acquisition in some species. While mutation was the main force driving polymorphism at the gene level, interspecies homologous recombination of large fragments expanding through several genes shaped Hrp2 cluster polymorphism. Horizontal gene transfer of the entire Hrp2 cluster also occurred. A reduced core effectome composed of xopF1, xopM, avrBs2 and xopR was identified that may allow commensal strains overcoming plant basal immunity. In contrast, stepwise accumulation of numerous type 3 effector genes was shown in successful pathogens responsible for epidemics. Our data suggest that capacity to intimately interact with plants through T3SS would be an ancestral trait of xanthomonads. Since its acquisition, T3SS has experienced a highly dynamic evolutionary history characterized by intense gene flux between species that may reflect its role in host adaptation.

Xanthomonas adaptation to common bean is associated with horizontal transfers of genes encoding TAL effectors

Background

Common bacterial blight is a devastating bacterial disease of common bean (Phaseolus vulgaris) caused by Xanthomonas citri pv. fuscans and Xanthomonas phaseoli pv. phaseoli. These phylogenetically distant strains are able to cause similar symptoms on common bean, suggesting that they have acquired common genetic determinants of adaptation to common bean. Transcription Activator-Like (TAL) effectors are bacterial type III effectors that are able to induce the expression of host genes to promote infection or resistance. Their capacity to bind to a specific host DNA sequence suggests that they are potential candidates for host adaption.

Results

To study the diversity of tal genes from Xanthomonas strains responsible for common bacterial blight of bean, whole genome sequences of 17 strains representing the diversity of X. citri pv. fuscans and X. phaseoli pv. phaseoli were obtained by single molecule real time sequencing. Analysis of these genomes revealed the existence of four tal genes named tal23A, tal20F, tal18G and tal18H, respectively. While tal20F and tal18G were chromosomic, tal23A and tal18H were carried on plasmids and shared between phylogenetically distant strains, therefore suggesting recent horizontal transfers of these genes between X. citri pv. fuscans and X. phaseoli pv. phaseoli strains. Strikingly, tal23A was present in all strains studied, suggesting that it played an important role in adaptation to common bean. In silico predictions of TAL effectors targets in the common bean genome suggested that TAL effectors shared by X. citri pv. fuscans and X. phaseoli pv. phaseoli strains target the promoters of genes of similar functions. This could be a trace of convergent evolution among TAL effectors from different phylogenetic groups, and comforts the hypothesis that TAL effectors have been implied in the adaptation to common bean.

Conclusions

Altogether, our results favour a model where plasmidic TAL effectors are able to contribute to host adaptation by being horizontally transferred between distant lineages.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4087-6) contains supplementary material, which is available to authorized users.

Transcriptomic analysis of nickel exposure in Sphingobium sp. ba1 cells using RNA-seq

Nickel acts as cofactor for a number of enzymes of many bacteria species. Its homeostasis is ensured by proteins working as ion efflux or accumulation systems. These mechanisms are also generally adopted to counteract life-threatening high extra-cellular Ni²⁺ concentrations. Little is known regarding nickel tolerance in the genus Sphingobium. We studied the response of the novel Sphingobium sp. ba1 strain, able to adapt to high Ni²⁺ concentrations. Differential gene expression in cells cultured in 10 mM Ni²⁺, investigated by RNA-seq analysis, identified 118 differentially expressed genes. Among the 90 up-regulated genes, a cluster including genes coding for nickel and other metal ion efflux systems (similar to either cnrCBA, nccCBA or cznABC) and for a NreB-like permease was found. Comparative analyses among thirty genomes of Sphingobium species show that this cluster is conserved only in two cases, while in the other genomes it is partially present or even absent. The differential expression of genes encoding proteins which could also work as Ni²⁺-accumulators (HupE/UreJ-like protein, NreA and components of TonB-associated transport and copper-homeostasis systems) was also detected. The identification of Sphingobium sp. ba1 strain adaptive mechanisms to nickel ions, can foster its possible use for biodegradation of poly-aromatic compounds in metal-rich environments.

Multilocus Variable Number of Tandem Repeat Analysis Reveals Multiple Introductions in Spain of Xanthomonas arboricola pv. pruni, the Causal Agent of Bacterial Spot Disease of Stone Fruits and Almond

Xanthomonas arboricola pv. pruni is the causal agent of the bacterial spot disease of stone fruits, almond and some ornamental Prunus species. In Spain it was first detected in 2002 and since then, several outbreaks have occurred in different regions affecting mainly Japanese plum, peach and almond, both in commercial orchards and nurseries. As the origin of the introduction(s) was unknown, we have assessed the genetic diversity of 239 X. arboricola pv. pruni strains collected from 11 Spanish provinces from 2002 to 2013 and 25 reference strains from international collections. We have developed an optimized multilocus variable number of tandem repeat analysis (MLVA) scheme targeting 18 microsatellites and five minisatellites. A high discriminatory power was achieved since almost 50% of the Spanish strains were distinguishable, confirming the usefulness of this genotyping technique at small spatio-temporal scales. Spanish strains grouped in 18 genetic clusters (conservatively delineated so that each cluster contained haplotype networks linked by up to quadruple-locus variations). Furthermore, pairwise comparisons among populations from different provinces showed a strong genetic differentiation. Our results suggest multiple introductions of this pathogen in Spain and redistribution through contaminated nursery propagative plant material.

Expression of the Sinorhizobium meliloti small RNA gene mmgR is controlled by the nitrogen source

Small non-coding regulatory RNAs (sRNAs) are key players in post-transcriptional regulation of gene expression. Hundreds of sRNAs have been identified in Sinorhizobium meliloti, but their biological function remains unknown for most of them. In this study, we characterized the expression pattern of the gene encoding the 77-nt sRNA MmgR in S. meliloti strain 2011. A chromosomal transcriptional reporter fusion (PmmgR-gfp) showed that the mmgR promoter is active along different stages of the interaction with alfalfa roots. In pure cultures, PmmgR-gfp activity paralleled the sRNA abundance indicating that mmgR expression is primarily controlled at the level of transcriptional initiation. PmmgR-gfp activity was higher during growth in rhizobial defined medium (RDM) than in TY medium. Furthermore, PmmgR-gfp was induced at 60 min after shifting growing cells from TY to RDM medium, i.e. shorter than the cell doubling time. In defined RDM medium containing NO3 (-), both PmmgR-gfp and MmgR level were repressed by the addition of tryptone or single amino acids, suggesting that mmgR expression depends on the cellular nitrogen (N) status. In silico analysis failed to detect conserved motifs upstream the promoter RNA polymerase binding site, but revealed a strongly conserved motif centered at -28 that may be linked to the observed regulatory pattern by the N source.

New Coffee Plant-Infecting Xylella fastidiosa Variants Derived via Homologous Recombination

Xylella fastidiosa is a xylem-limited phytopathogenic bacterium endemic to the Americas that has recently emerged in Asia and Europe. Although this bacterium is classified as a quarantine organism in the European Union, importation of plant material from contaminated areas and latent infection in asymptomatic plants have engendered its inevitable introduction. In 2012, four coffee plants (Coffea arabica and Coffea canephora) with leaf scorch symptoms growing in a confined greenhouse were detected and intercepted in France. After identification of the causal agent, this outbreak was eradicated. Three X. fastidiosa strains were isolated from these plants, confirming a preliminary identification based on immunology. The strains were characterized by multiplex PCR and by multilocus sequence analysis/typing (MLSA-MLST) based on seven housekeeping genes. One strain, CFBP 8073, isolated from C. canephora imported from Mexico, was assigned to X. fastidiosa subsp. fastidiosa/X. fastidiosa subsp. sandyi. This strain harbors a novel sequence type (ST) with novel alleles at two loci. The two other strains, CFBP 8072 and CFBP 8074, isolated from Coffea arabica imported from Ecuador, were allocated to X. fastidiosa subsp. pauca. These two strains shared a novel ST with novel alleles at two loci. These MLST profiles showed evidence of recombination events. We provide genome sequences for CFBP 8072 and CFBP 8073 strains. Comparative genomic analyses of these two genome sequences with publicly available X. fastidiosa genomes, including the Italian strain CoDiRO, confirmed these phylogenetic positions and provided candidate alleles for coffee plant adaptation. This study demonstrates the global diversity of X. fastidiosa and highlights the diversity of strains isolated from coffee plants.

Comparative Genomics of Pathogenic and Nonpathogenic Strains of Xanthomonas arboricola Unveil Molecular and Evolutionary Events Linked to Pathoadaptation

The bacterial species Xanthomonas arboricola contains plant pathogenic and nonpathogenic strains. It includes the pathogen X. arboricola pv. juglandis, causing the bacterial blight of Juglans regia. The emergence of a new bacterial disease of J. regia in France called vertical oozing canker (VOC) was previously described and the causal agent was identified as a distinct genetic lineage within the pathovar juglandis. Symptoms on walnut leaves and fruits are similar to those of a bacterial blight but VOC includes also cankers on trunk and branches. In this work, we used comparative genomics and physiological tests to detect differences between four X. arboricola strains isolated from walnut tree: strain CFBP 2528 causing walnut blight (WB), strain CFBP 7179 causing VOC and two nonpathogenic strains, CFBP 7634 and CFBP 7651, isolated from healthy walnut buds. Whole genome sequence comparisons revealed that pathogenic strains possess a larger and wider range of mobile genetic elements than nonpathogenic strains. One pathogenic strain, CFBP 7179, possessed a specific integrative and conjugative element (ICE) of 95 kb encoding genes involved in copper resistance, transport and regulation. The type three effector repertoire was larger in pathogenic strains than in nonpathogenic strains. Moreover, CFBP 7634 strain lacked the type three secretion system encoding genes. The flagellar system appeared incomplete and nonfunctional in the pathogenic strain CFBP 2528. Differential sets of chemoreceptor and different repertoires of genes coding adhesins were identified between pathogenic and nonpathogenic strains. Besides these differences, some strain-specific differences were also observed. Altogether, this study provides valuable insights to highlight the mechanisms involved in ecology, environment perception, plant adhesion and interaction, leading to the emergence of new strains in a dynamic environment.

Whole-Genome Sequencing of 10 Pseudomonas syringae Strains Representing Different Host Range Spectra

Pseudomonas syringae is a ubiquitous bacterium that readily persists in environmental habitats as a saprophyte and also is responsible for numerous diseases of crops. Here, we report the whole-genome sequences of 10 strains isolated from both woody and herbaceous plants that will contribute to the elucidation of the determinants of their host ranges.

Comparative Genomics of Pathogenic and Nonpathogenic Strains of Xanthomonas arboricola Unveil Molecular and Evolutionary Events Linked to Pathoadaptation

The bacterial species Xanthomonas arboricola contains plant pathogenic and nonpathogenic strains. It includes the pathogen X. arboricola pv. juglandis, causing the bacterial blight of Juglans regia. The emergence of a new bacterial disease of J. regia in France called vertical oozing canker (VOC) was previously described and the causal agent was identified as a distinct genetic lineage within the pathovar juglandis. Symptoms on walnut leaves and fruits are similar to those of a bacterial blight but VOC includes also cankers on trunk and branches. In this work, we used comparative genomics and physiological tests to detect differences between four X. arboricola strains isolated from walnut tree: strain CFBP 2528 causing walnut blight (WB), strain CFBP 7179 causing VOC and two nonpathogenic strains, CFBP 7634 and CFBP 7651, isolated from healthy walnut buds. Whole genome sequence comparisons revealed that pathogenic strains possess a larger and wider range of mobile genetic elements than nonpathogenic strains. One pathogenic strain, CFBP 7179, possessed a specific integrative and conjugative element (ICE) of 95 kb encoding genes involved in copper resistance, transport and regulation. The type three effector repertoire was larger in pathogenic strains than in nonpathogenic strains. Moreover, CFBP 7634 strain lacked the type three secretion system encoding genes. The flagellar system appeared incomplete and nonfunctional in the pathogenic strain CFBP 2528. Differential sets of chemoreceptor and different repertoires of genes coding adhesins were identified between pathogenic and nonpathogenic strains. Besides these differences, some strain-specific differences were also observed. Altogether, this study provides valuable insights to highlight the mechanisms involved in ecology, environment perception, plant adhesion and interaction, leading to the emergence of new strains in a dynamic environment.