Understanding pathogenic Burkholderia glumae metabolic and signaling pathways within rice tissues through in vivo transcriptome analyses

Genome-guided comparative in planta transcriptome analyses for identifying cross-species common virulence factors in bacterial phytopathogens

Plant bacterial disease is a complex outcome achieved through a combination of virulence factors that are activated during infection. However, the common virulence factors across diverse plant pathogens are largely uncharacterized. Here, we established a pan-genome shared across the following plant pathogens: Burkholderia glumae, Ralstonia solanacearum, and Xanthomonas oryzae pv. oryzae. By overlaying in planta transcriptomes onto the pan-genome, we investigated the expression profiles of common genes during infection. We found over 70% of identical patterns for genes commonly expressed by the pathogens in different plant hosts or infection sites. Co-expression patterns revealed the activation of a signal transduction cascade to recognize and respond to external changes within hosts. Using mutagenesis, we uncovered a relationship between bacterial virulence and functions highly conserved and shared in the studied genomes of the bacterial phytopathogens, including flagellar biosynthesis protein, C4-dicarboxylate ABC transporter, 2-methylisocitrate lyase, and protocatechuate 3,4-dioxygenase (PCD). In particular, the disruption of PCD gene led to attenuated virulence in all pathogens and significantly affected phytotoxin production in B. glumae. This PCD gene was ubiquitously distributed in most plant pathogens with high homology. In conclusion, our results provide cross-species in planta models for identifying common virulence factors, which can be useful for the protection of crops against diverse pathogens.

Next Generation Sequencing and Comparative Genomic Analysis Reveal Extreme Plasticity of Two Burkholderia glumae Strains HN1 and HN2

Burkholderia glumae is an important rice pathogen, thus the genomic and evolutionary history may be helpful to control this notorious pathogen. Here, we present two complete genomes of the B. glumae strains HN1 and HN2, which were isolated from diseased rice seed in China. Average nucleotide identity (ANI) analysis shows greater than 99% similarity of the strains HN1 and HN2 with other published B. glumae genomes. Genomic annotation revealed that the genome of strain HN1 consists of five replicons (6,680,415 bp) with an overall G + C content of 68.06%, whereas the genome of strain HN2 comprises of three replicons (6,560,085 bp) with an overall G + C content of 68.34%. The genome of HN1 contains 5434 protein-coding genes, 351 pseudogenes, and 1 CRISPR, whereas the genome of HN2 encodes 5278 protein-coding genes, 357 pseudogenes, and 2 CRISPR. Both strains encode many pathogenic-associated genes (143 genes in HN1 vs. 141 genes in HN2). Moreover, comparative genomic analysis shows the extreme plasticity of B. glumae, which may contribute to its pathogenicity. In total, 259 single-copy genes were affected by positive selection. These genes may contribute to the adaption to different environments. Notably, six genes were characterized as virulence factors which may be an additional way to assist the pathogenicity of B. glumae.

Bacterial Panicle Blight and Burkholderia glumae: From Pathogen Biology to Disease Control

Bacterial panicle blight (BPB), caused by the bacterium Burkholderia glumae, has affected rice production worldwide. Despite its importance, neither the disease nor the causal agent are well understood. Moreover, methods to manage BPB are still lacking. Nevertheless, the emerging importance of this pathogen has stimulated research to identify the mechanisms of pathogenicity, to gain insight into plant disease resistance, and to develop strategies to manage the disease. In this review, we consolidate current information regarding the virulence factors that have been identified in B. glumae and present a model of the disease and the pathogen. We also provide an update on the current research status to develop methods to control the disease especially through biological control approaches and through the development of resistant cultivars.

Type VI secretion systems of plant-pathogenic Burkholderia glumae BGR1 play a functionally distinct role in interspecies interactions and virulence

In the environment, bacteria show close association, such as interspecies interaction, with other bacteria as well as host organisms. The type VI secretion system (T6SS) in gram-negative bacteria is involved in bacterial competition or virulence. The plant pathogen Burkholderia glumae BGR1, causing bacterial panicle blight in rice, has four T6SS gene clusters. The presence of at least one T6SS gene cluster in an organism indicates its distinct role, like in the bacterial and eukaryotic cell targeting system. In this study, deletion mutants targeting four tssD genes, which encode the main component of T6SS needle formation, were constructed to functionally dissect the four T6SSs in B. glumae BGR1. We found that both T6SS group_4 and group_5, belonging to the eukaryotic targeting system, act independently as bacterial virulence factors toward host plants. In contrast, T6SS group_1 is involved in bacterial competition by exerting antibacterial effects. The ΔtssD1 mutant lost the antibacterial effect of T6SS group_1. The ΔtssD1 mutant showed similar virulence as the wild-type BGR1 in rice because the ΔtssD1 mutant, like the wild-type BGR1, still has key virulence factors such as toxin production towards rice. However, metagenomic analysis showed different bacterial communities in rice infected with the ΔtssD1 mutant compared to wild-type BGR1. In particular, the T6SS group_1 controls endophytic plant-associated bacteria such as Luteibacter and Dyella in rice plants and may have an advantage in competing with endophytic plant-associated bacteria for settlement inside rice plants in the environment. Thus, B. glumae BGR1 causes disease using T6SSs with functionally distinct roles.

Characterization of Burkholderia glumae Putative Virulence Factor 11 (PVF11) via Yeast Two-Hybrid Interaction and Phenotypic Analysis

In this study, PVF11 was selected among 20 candidate pathogenesis-related genes in Burkholderia glumae based on its effect on virulence to rice. PVF11 was found to interact with several plant defense-related WRKY proteins as evidenced through yeast-two hybrid analysis (Y2H). Moreover, PVF11 showed interactions with abiotic and biotic stress response-related rice proteins, as shown by genome-wide Y2H screening employing PVF11 and a cDNA library from B. glumae-infected rice. To confirm the effect of PVF11 on B. glumae virulence, in planta assays were conducted at different stages of rice growth. As a result, a PVF11-defective mutant showed reduced virulence in rice seedlings and stems but not in rice panicles, indicating that PVF11 involvement in B. glumae virulence in rice is stage-dependent.

Plant-Growth Promotion and Biocontrol Properties of Three Streptomyces spp. Isolates to Control Bacterial Rice Pathogens

Bacterial Panicle Blight caused by Burkholderia glumae is a major disease of rice, which has dramatically affected rice production around the world in the last years. In this study we describe the assessment of three Streptomyces isolates as biocontrol agents for B. glumae. Additionally, the presence of other plant-growth promoting abilities and their possible beneficial effects upon their inoculation on rice plants was evaluated as an ecological analysis for their future inoculation in rice crops. Two isolates (A20 and 5.1) inhibited growth of virulent B. glumae strains, as well as a wide range of bacterial and fungal species, while a third strain (7.1) showed only antifungal activity. In vitro tests demonstrated the ability of these strains to produce siderophores, Indoleacetic acid (IAA), extracellular enzymes and solubilizing phosphate. Greenhouse experiments with two rice cultivars indicated that Streptomyces A20 is able to colonize rice plants and promote plant growth in both cultivars. Furthermore, an egfp tagged mutant was generated and colonization experiments were performed, indicating that Streptomyces A20 –GFP was strongly associated with root hairs, which may be related to the plant growth promotion observed in the gnotobiotic experiments. In order to characterize the antimicrobial compounds produced by strain A20 bacteria, mass spectrometry analyses were performed. This technique indicated that A20 produced several antimicrobial compounds with sizes below 3 kDa and three of these molecules were identified as Streptotricins D, E and F. These findings indicate the potential of Streptomyces A20 as a biocontrol inoculant to protect rice plants against bacterial diseases.

Sugarcane apoplast fluid modulates the global transcriptional profile of the diazotrophic bacteria Paraburkholderia tropica strain Ppe8

The stalk apoplast fluid of sugarcane contains different sugars, organic acids and amino acids that may supply the demand for carbohydrates by endophytic bacteria including diazotrophs P. tropica (syn. B. tropica) strain Ppe8, isolated from sugarcane, is part of the bacterial consortium recommended as inoculant to sugarcane. However, little information has been accumulated regarding this plant-bacterium interaction considering that it colonizes internal sugarcane tissues. Here, we made use of the RNA-Seq transcriptomic analysis to study the influence of sugarcane stalk apoplast fluid on Ppe8 gene expression. The bacterium was grown in JMV liquid medium (100 ml), divided equally and then supplemented with 50 ml of fresh JMV medium or 50 ml of apoplast fluid extracted from sugarcane variety RB867515. Total RNA was extracted 2 hours later, the rRNAs were depleted and mRNAs used to construct libraries to sequence the fragments using Ion Torrent technology. The mapping and statistical analysis were carried out with CLC Genomics Workbench software. The RNA-seq data was validated by RT-qPCR using the reference genes fliP1, paaF, and groL. The data analysis showed that 544 genes were repressed and 153 genes were induced in the presence of apoplast fluid. Genes that induce plant defense responses, genes related to chemotaxis and movements were repressed in the presence of apoplast fluid, indicating that strain Ppe8 recognizes the apoplast fluid as a plant component. The expression of genes involved in bacterial metabolism was regulated (up and down), suggesting that the metabolism of strain Ppe8 is modulated by the apoplast fluid. These results suggest that Ppe8 alters its gene expression pattern in the presence of apoplast fluid mainly in order to use compounds present in the fluid as well as to avoid the induction of plant defense mechanisms. This is a pioneer study showing the role played by the sugarcane apoplast fluid on the global modulation of genes in P. tropica strain Ppe8.

Whole-Transcriptome Sequence Analysis of Verbena bonariensis in Response to Drought Stress

Drought is an important abiotic factor that threatens the growth and development of plants. Verbena bonariensis is a widely used landscape plant with a very high ornamental value. We found that Verbena has drought tolerance in production practice, so in order to delve into its mechanism of drought resistance and screen out its drought-resistance genes, we used the RNA-Seq platform to perform a de novo transcriptome assembly to analyze Verbena transcription response to drought stress. By high-throughput sequencing with Illumina Hiseq Xten, a total of 44.59 Gb clean data was obtained from T01 (control group) and T02 (drought experiment group). After assembly, 111,313 unigenes were obtained, and 53,757 of them were annotated by compared databases. In this study, 4829 differentially expressed genes were obtained, of which 4165 were annotated. We performed GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyses, and explored a lot of differently expressed genes related to plant energy production, hormone synthesis, cell signal transduction, and metabolism to understand the stress response of Verbena in drought stress. In addition, we also found that a series of TFs related to drought-resistance of Verbena and provide excellent genetic resources for improving the drought tolerance of crops.

Trehalase plays a role in macrophage colonization and virulence of Burkholderia pseudomallei in insect and mammalian hosts

Trehalose is a disaccharide formed from two glucose molecules. This sugar molecule can be isolated from a range of organisms including bacteria, fungi, plants and invertebrates. Trehalose has a variety of functions including a role as an energy storage molecule, a structural component of glycolipids and plays a role in the virulence of some microorganisms. There are many metabolic pathways that control the biosynthesis and degradation of trehalose in different organisms. The enzyme trehalase forms part of a pathway that converts trehalose into glucose. In this study we set out to investigate whether trehalase plays a role in both stress adaptation and virulence of Burkholderia pseudomallei. We show that a trehalase deletion mutant (treA) had increased tolerance to thermal stress and produced less biofilm than the wild type B. pseudomallei K96243 strain. We also show that the ΔtreA mutant has reduced ability to survive in macrophages and that it is attenuated in both Galleria mellonella (wax moth larvae) and a mouse infection model. This is the first report that trehalase is important for bacterial virulence.

Rice Sheath Rot: An Emerging Ubiquitous Destructive Disease Complex

Around one century ago, a rice disease characterized mainly by rotting of sheaths was reported in Taiwan. The causal agent was identified as Acrocylindrium oryzae, later known as Sarocladium oryzae. Since then it has become clear that various other organisms can cause similar disease symptoms, including Fusarium sp. and fluorescent pseudomonads. These organisms have in common that they produce a range of phytotoxins that induce necrosis in plants. The same agents also cause grain discoloration, chaffiness, and sterility and are all seed-transmitted. Rice sheath rot disease symptoms are found in all rice-growing areas of the world. The disease is now getting momentum and is considered as an important emerging rice production threat. The disease can lead to variable yield losses, which can be as high as 85%. This review aims at improving our understanding of the disease etiology of rice sheath rot and mainly deals with the three most reported rice sheath rot pathogens: S. oryzae, the Fusarium fujikuroi complex, and Pseudomonas fuscovaginae. Causal agents, pathogenicity determinants, interactions among the various pathogens, epidemiology, geographical distribution, and control options will be discussed.

Genome-wide RNA sequencing analysis of quorum sensing-controlled regulons in the plant-associated Burkholderia glumae PG1 strain

Burkholderia glumae PG1 is a soil-associated motile plant-pathogenic bacterium possessing a cell density-dependent regulation system called quorum sensing (QS). Its genome contains three genes, here designated bgaI1 to bgaI3, encoding distinct autoinducer-1 (AI-1) synthases, which are capable of synthesizing QS signaling molecules. Here, we report on the construction of B. glumae PG1 ΔbgaI1, ΔbgaI2, and ΔbgaI3 mutants, their phenotypic characterization, and genome-wide transcriptome analysis using RNA sequencing (RNA-seq) technology. Knockout of each of these bgaI genes resulted in strongly decreased motility, reduced extracellular lipase activity, a reduced ability to cause plant tissue maceration, and decreased pathogenicity. RNA-seq analysis of all three B. glumae PG1 AI-1 synthase mutants performed in the transition from exponential to stationary growth phase revealed differential expression of a significant number of predicted genes. In comparison with the levels of gene expression by wild-type strain B. glumae PG1, 481 genes were differentially expressed in the ΔbgaI1 mutant, 213 were differentially expressed in the ΔbgaI2 mutant, and 367 were differentially expressed in the ΔbgaI3 mutant. Interestingly, only a minor set of 78 genes was coregulated in all three mutants. The majority of the QS-regulated genes were linked to metabolic activities, and the most pronounced regulation was observed for genes involved in rhamnolipid and Flp pilus biosynthesis and the type VI secretion system and genes linked to a clustered regularly interspaced short palindromic repeat (CRISPR)-cas gene cluster.

Metatranscriptomics reveals metabolic adaptation and induction of virulence factors by Haemophilus parasuis during lung infection

Haemophilus parasuis is a common inhabitant of the upper respiratory tract of pigs, and the causative agent of Glässer’s disease. This disease is characterized by polyserositis and arthritis, produced by the severe inflammation caused by the systemic spread of the bacterium. After an initial colonization of the upper respiratory tract, H. parasuis enters the lung during the early stages of pig infection. In order to study gene expression at this location, we sequenced the ex vivo and in vivo H. parasuis Nagasaki transcriptome in the lung using a metatranscriptomic approach. Comparison of gene expression under these conditions with that found in conventional plate culture showed generally reduced expression of genes associated with anabolic and catabolic pathways, coupled with up-regulation of membrane-related genes involved in carbon acquisition, iron binding and pathogenesis. Some of the up-regulated membrane genes, including ABC transporters, virulence-associated autotransporters (vtaAs) and several hypothetical proteins, were only present in virulent H. parasuis strains, highlighting their significance as markers of disease potential. Finally, the analysis also revealed the presence of numerous antisense transcripts with possible roles in gene regulation. In summary, this data sheds some light on the scarcely studied in vivo transcriptome of H. parasuis, revealing nutritional virulence as an adaptive strategy for host survival, besides induction of classical virulence factors.