Genome Context Influences Evolutionary Flexibility of Nearly Identical Type III Effectors in Two Phytopathogenic Pseudomonads

Beav: a bacterial genome and mobile element annotation pipeline

Comprehensive and accurate genome annotation is crucial for inferring the predicted functions of an organism. Numerous tools exist to annotate genes, gene clusters, mobile genetic elements, and other diverse features. However, these tools and pipelines can be difficult to install and run, be specialized for a particular element or feature, or lack annotations for larger elements that provide important genomic context. Integrating results across analyses is also important for understanding gene function. To address these challenges, we present the Beav annotation pipeline. Beav is a command-line tool that automates the annotation of bacterial genome sequences, mobile genetic elements, molecular systems and gene clusters, key regulatory features, and other elements. Beav uses existing tools in addition to custom models, scripts, and databases to annotate diverse elements, systems, and sequence features. Custom databases for plant-associated microbes are incorporated to improve annotation of key virulence and symbiosis genes in agriculturally important pathogens and mutualists. Beav includes an optional Agrobacterium-specific pipeline that identifies and classifies oncogenic plasmids and annotates plasmid-specific features. Following the completion of all analyses, annotations are consolidated to produce a single comprehensive output. Finally, Beav generates publication-quality genome and plasmid maps. Beav is on Bioconda and is available for download at https://github.com/weisberglab/beav.

IMPORTANCE

Annotation of genome features, such as the presence of genes and their predicted function, or larger loci encoding secretion systems or biosynthetic gene clusters, is necessary for understanding the functions encoded by an organism. Genomes can also host diverse mobile genetic elements, such as integrative and conjugative elements and/or phages, that are often not annotated by existing pipelines. These elements can horizontally mobilize genes encoding for virulence, antimicrobial resistance, or other adaptive functions and alter the phenotype of an organism. We developed a software pipeline, called Beav, that combines new and existing tools for the comprehensive annotation of these and other major features. Existing pipelines often misannotate loci important for virulence or mutualism in plant-associated bacteria. Beav includes custom databases and optional workflows for the improved annotation of plant-associated bacteria. Beav is designed to be easy to install and run, making comprehensive genome annotation broadly available to the research community.

Inter-species interactions between two bacterial flower commensals and a floral pathogen reduce disease incidence and alter pathogen activity

Flowers are colonized by a diverse community of microorganisms that can alter plant health and interact with floral pathogens. Erwinia amylovora is a flower-inhabiting bacterium and a pathogen that infects different plant species, including Malus × domestica (apple). Previously, we showed that the co-inoculation of two bacterial strains, members of the genera Pseudomonas and Pantoea, isolated from apple flowers, reduced disease incidence caused by this floral pathogen. Here, we decipher the ecological interactions between the two flower-associated bacteria and E. amylovora in field experimentation and in vitro co-cultures. The two flower commensal strains did not competitively exclude E. amylovora from the stigma habitat, as both bacteria and the pathogen co-existed on the stigma of apple flowers and in vitro. This suggests that plant protection might be mediated by other mechanisms than competitive niche exclusion. Using a synthetic stigma exudation medium, ternary co-culture of the bacterial strains led to a substantial alteration of gene expression in both the pathogen and the two microbiota members. Importantly, the gene expression profiles for the ternary co-culture were not just additive from binary co-cultures, suggesting that some functions only emerged in multipartite co-culture. Additionally, the ternary co-culture of the strains resulted in a stronger acidification of the growth milieu than mono- or binary co-cultures, pointing to another emergent property of co-inoculation. Our study emphasizes the critical role of emergent properties mediated by inter-species interactions within the plant holobiont and their potential impact on plant health and pathogen behavior.

IMPORTANCE

Fire blight, caused by Erwinia amylovora, is one of the most important plant diseases of pome fruits. Previous work largely suggested plant microbiota commensals suppressed disease by antagonizing pathogen growth. However, inter-species interactions of multiple flower commensals and their influence on pathogen activity and behavior have not been well studied. Here, we show that co-inoculating two bacterial strains that naturally colonize the apple flowers reduces disease incidence. We further demonstrate that the interactions between these two microbiota commensals and the floral pathogen led to the emergence of new gene expression patterns and a strong alteration of the external pH, factors that may modify the pathogen's behavior. Our findings emphasize the critical role of emergent properties mediated by inter-species interactions between plant microbiota and plant pathogens and their impact on plant health.

Mobile Genetic Element Flexibility as an Underlying Principle to Bacterial Evolution

Mobile genetic elements are key to the evolution of bacteria and traits that affect host and ecosystem health. Here, we use a framework of a hierarchical and modular system that scales from genes to populations to synthesize recent findings on mobile genetic elements (MGEs) of bacteria. Doing so highlights the role that emergent properties of flexibility, robustness, and genetic capacitance of MGEs have on the evolution of bacteria. Some of their traits can be stored, shared, and diversified across different MGEs, taxa of bacteria, and time. Collectively, these properties contribute to maintaining functionality against perturbations while allowing changes to accumulate in order to diversify and give rise to new traits. These properties of MGEs have long challenged our abilities to study them. Implementation of new technologies and strategies allows for MGEs to be analyzed in new and powerful ways.

Comparative sequence analysis of pPATH pathogenicity plasmids in Pantoea agglomerans gall-forming bacteria

Acquisition of the pathogenicity plasmid pPATH that encodes a type III secretion system (T3SS) and effectors (T3Es) has likely led to the transition of a non-pathogenic bacterium into the tumorigenic pathogen Pantoea agglomerans. P. agglomerans pv. gypsophilae (Pag) forms galls on gypsophila (Gypsophila paniculata) and triggers immunity on sugar beet (Beta vulgaris), while P. agglomerans pv. betae (Pab) causes galls on both gypsophila and sugar beet. Draft sequences of the Pag and Pab genomes were previously generated using the MiSeq Illumina technology and used to determine partial T3E inventories of Pab and Pag. Here, we fully assembled the Pab and Pag genomes following sequencing with PacBio technology and carried out a comparative sequence analysis of the Pab and Pag pathogenicity plasmids pPATHpag and pPATHpab. Assembly of Pab and Pag genomes revealed a ~4 Mbp chromosome with a 55% GC content, and three and four plasmids in Pab and Pag, respectively. pPATHpag and pPATHpab share 97% identity within a 74% coverage, and a similar GC content (51%); they are ~156 kb and ~131 kb in size and consist of 198 and 155 coding sequences (CDSs), respectively. In both plasmids, we confirmed the presence of highly similar gene clusters encoding a T3SS, as well as auxin and cytokinins biosynthetic enzymes. Three putative novel T3Es were identified in Pab and one in Pag. Among T3SS-associated proteins encoded by Pag and Pab, we identified two novel chaperons of the ShcV and CesT families that are present in both pathovars with high similarity. We also identified insertion sequences (ISs) and transposons (Tns) that may have contributed to the evolution of the two pathovars. These include seven shared IS elements, and three ISs and two transposons unique to Pab. Finally, comparative sequence analysis revealed plasmid regions and CDSs that are present only in pPATHpab or in pPATHpag. The high similarity and common features of the pPATH plasmids support the hypothesis that the two strains recently evolved into host-specific pathogens.

Single-cell profiling of Arabidopsis leaves to Pseudomonas syringae infection

Plant response to pathogen infection varies within a leaf, yet this heterogeneity is not well resolved. We expose Arabidopsis to Pseudomonas syringae or mock treatment and profile >11,000 individual cells using single-cell RNA sequencing. Integrative analysis of cell populations from both treatments identifies distinct pathogen-responsive cell clusters exhibiting transcriptional responses ranging from immunity to susceptibility. Pseudotime analyses through pathogen infection reveals a continuum of disease progression from an immune to a susceptible state. Confocal imaging of promoter-reporter lines for transcripts enriched in immune cell clusters shows expression surrounding substomatal cavities colonized or in close proximity to bacterial colonies, suggesting that cells within immune clusters represent sites of early pathogen invasion. Susceptibility clusters exhibit more general localization and are highly induced at later stages of infection. Overall, our work shows cellular heterogeneity within an infected leaf and provides insight into plant differential response to infection at a single-cell level.