Comparative genome analysis of rice-pathogenic Burkholderia provides insight into capacity to adapt to different environments and hosts

Background

In addition to human and animal diseases, bacteria of the genus Burkholderia can cause plant diseases. The representative species of rice-pathogenic Burkholderia are Burkholderia glumae, B. gladioli, and B. plantarii, which primarily cause grain rot, sheath rot, and seedling blight, respectively, resulting in severe reductions in rice production. Though Burkholderia rice pathogens cause problems in rice-growing countries, comprehensive studies of these rice-pathogenic species aiming to control Burkholderia-mediated diseases are only in the early stages.

Results

We first sequenced the complete genome of B. plantarii ATCC 43733^T. Second, we conducted comparative analysis of the newly sequenced B. plantarii ATCC 43733^T genome with eleven complete or draft genomes of B. glumae and B. gladioli strains. Furthermore, we compared the genome of three rice Burkholderia pathogens with those of other Burkholderia species such as those found in environmental habitats and those known as animal/human pathogens. These B. glumae, B. gladioli, and B. plantarii strains have unique genes involved in toxoflavin or tropolone toxin production and the clustered regularly interspaced short palindromic repeats (CRISPR)-mediated bacterial immune system. Although the genome of B. plantarii ATCC 43733^T has many common features with those of B. glumae and B. gladioli, this B. plantarii strain has several unique features, including quorum sensing and CRISPR/CRISPR-associated protein (Cas) systems.

Conclusions

The complete genome sequence of B. plantarii ATCC 43733^T and publicly available genomes of B. glumae BGR1 and B. gladioli BSR3 enabled comprehensive comparative genome analyses among three rice-pathogenic Burkholderia species responsible for tissue rotting and seedling blight. Our results suggest that B. glumae has evolved rapidly, or has undergone rapid genome rearrangements or deletions, in response to the hosts. It also, clarifies the unique features of rice pathogenic Burkholderia species relative to other animal and human Burkholderia species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1558-5) contains supplementary material, which is available to authorized users.     

Cucumis sativus L-type lectin receptor kinase (CsLecRK) gene family response to Phytophthora melonis, Phytophthora capsici and water immersion in disease resistant and susceptible cucumber cultivars

L-type lectin receptor kinase (LecRK) proteins are an important family involved in diverse biological processes such as pollen development, senescence, wounding, salinity and especially in innate immunity in model plants such as Arabidopsis and tobacco. Till date, LecRK proteins or genes of cucumber have not been reported. In this study, a total of 25 LecRK genes were identified in the cucumber genome, unequally distributed across its seven chromosomes. According to similarity comparison of their encoded proteins, the Cucumis sativus LecRK (CsLecRK) genes were classified into six major clades (from Clade I to CladeVI). Expression of CsLecRK genes were tested using QRT-PCR method and the results showed that 25 CsLecRK genes exhibited different responses to abiotic (water immersion) and biotic (Phytophthora melonis and Phytophthora capsici inoculation) stresses, as well as that between disease resistant cultivar (JSH) and disease susceptible cultivar (B80). Among the 25 CsLecRK genes, we found CsLecRK6.1 was especially induced by P. melonis and P. capsici in JSH plants. All these results suggested that CsLecRK genes may play important roles in biotic and abiotic stresses.     

Cloning and truncation modification of trehalose-6-phosphate synthase gene from Selaginella pulvinata

A homologous sequence was amplified from resurrection plant Selaginella pulvinta by RACE technique, proved to be the full-length cDNA of trehalose-6-phosphate synthase gene by homologous alignment and yeast complementation assay, and nominated as SpTPS1 gene. The open reading frame of this gene was truncated 225 bp at the 5′-end, resulting the N-terminal truncation modification of 75 amino acids for its encoding protein. The TPS1 deletion mutant strain YSH290 of the brewer's yeast transformed by the truncated gene SpTPS1Δ and its original full-length version restored growth on the medium with glucose as a sole carbon source and displayed growth curves with no significant difference, indicating their encoding proteins functioning as TPS enzyme. The TPS activity of the mutant strain transformed by the truncated gene SpTPS1Δ was about six fold higher than that transformed by its original version, reasoning that the extra N-terminal extension of the full-length amino acid sequence acts as an inhibitory domain to trehalose synthesis. However, the trehalose accumulation of the mutant strain transformed by the truncated gene SpTPS1Δ was only 8% higher than that transformed by its original version. This result is explained by the feedback balance of trehalose content coordinated by the comparative activities between trehalose synthase and trehalase. The truncated gene SpTPS1Δ is suggested to be used in transgenic operation, together with the inhibition of trehalase activity by the application of validamycin A or genetic deficiency of the endogenous trehalase gene, for the enhancement of trehalose accumulation and improvement of abiotic tolerance in transgenic plants.     

Genome-wide transcriptional changes of ramie (Boehmeria nivea L. Gaud) in response to root-lesion nematode infection

Ramie fiber extracted from stem bark is one of the most important natural fibers. The root-lesion nematode (RLN) Pratylenchus coffeae is a major ramie pest and causes large fiber yield losses in China annually. The response mechanism of ramie to RLN infection is poorly understood. In this study, we identified genes that are potentially involved in the RLN-resistance in ramie using Illumina pair-end sequencing in two RLN-infected plants (Inf1 and Inf2) and two control plants (CO1 and CO2). Approximately 56.3, 51.7, 43.4, and 45.0 million sequencing reads were generated from the libraries of CO1, CO2, Inf1, and Inf2, respectively. De novo assembly for these 196 million reads yielded 50,486 unigenes with an average length of 853.3 bp. A total of 24,820 (49.2%) genes were annotated for their function. Comparison of gene expression levels between CO and Inf ramie revealed 777 differentially expressed genes (DEGs). The expression levels of 12 DEGs were further confirmed by real-time quantitative PCR (qRT-PCR). Pathway enrichment analysis showed that three pathways (phenylalanine metabolism, carotenoid biosynthesis, and phenylpropanoid biosynthesis) were strongly influenced by RLN infection. A series of candidate genes and pathways that may contribute to the defense response against RLN in ramie will be helpful for further improving resistance to RLN infection.     

An antibacterial pyrazole derivative from Burkholderia glumae, a bacterial pathogen of rice

Burkholderia glumae, a bacterial pathogen on rice, produced compounds in liquid culture that, in agar diffusion assays, gave strong inhibitory action against Erwinia amylovora, the bacterium responsible for fire blight disease of apple and pear trees. Products were isolated from culture medium by cation exchange and then purified by bioassay-guided chromatographic methods. Two major products were obtained, one of which was not active when fully purified. Each product showed a single ninhydrin-staining spot on TLC and a single HPLC peak. The non-active product was deduced from NMR, MS, and chemical data, to be the tripeptide l-alanyl-l-homoserinyl-l-aspartate. The NMR data for the active product demonstrated that it contained the same tripeptide, but functionalised at the β-carboxyl of the C-terminal aspartate, by a moiety that provided an additional 98 mass units to the parent tripeptide. Various data led to the interpretation that this moiety was a highly unusual oxygenated pyrazole structure, and thus the bioactive product was deduced to be 3-[l-alanyl-l-homoserinyl-l-aspartyl-β-carboxy]-4-hydroxy-5-oxopyrazole. This compound was found to inhibit the growth of a number of different bacterial species.     

Molecular evolution of the HD-ZIP I gene family in legume genomes

Homeodomain leucine zipper I (HD-ZIP I) genes were used to increase the plasticity of plants by mediating external signals and regulating growth in response to environmental conditions. The way genomic histories drove the evolution of the HD-ZIP I family in legume species was described; HD-ZIP I genes were searched in Lotus japonicus, Medicago truncatula, Cajanus cajan and Phaseolus vulgaris, and then divided into five clades through phylogenetic analysis. Microsynteny analysis was made based on genomic segments containing the HD-ZIP I genes. Some pairs turned out to conform with syntenic genome regions, while others corresponded to those that were inverted, expanded, or contracted after the divergence of legumes. Besides, we dated their duplications by Ks analysis and demonstrated that all the blocks were formed after the monocot–dicot split; we observed Ka/Ks ratios representing strong purifying selections in the four legume species which might have been followed by gene loss and rearrangement.     

Comparative analysis of the virulence of invertebrate and mammalian pathogenic bacteria in the oral insect infection model Galleria mellonella

Infection of Galleria mellonella by feeding a mixture of Bacillus thuringiensis spores or vegetative bacteria in association with the toxin Cry1C results in high levels of larval mortality. Under these conditions the toxin or bacteria have minimal effects on the larva when inoculated separately. In order to evaluate whether G. mellonella can function as an oral infection model for human and entomo-bacterial pathogens, we tested strains of Bacillus cereus, Bacillus anthracis, Enterococcus faecalis, Listeria monocytogenes, Pseudomonas aeruginosa and a Drosophila targeting Pseudomonas entomophila strain. Six B. cereus strains (5 diarrheal, 1 environmental isolate) were first screened in 2nd instar G. mellonella larvae by free ingestion and four of them were analyzed by force-feeding 5th instar larvae. The virulence of these B. cereus strains did not differ from the B. thuringiensis virulent reference strain 407Cry⁻ with the exception of strain D19 (NVH391/98) that showed a lower virulence. Following force-feeding, 5th instar G. mellonella larvae survived infection with B. anthracis, L. monocytogenes, E. faecalis and P. aeruginosa strains in contrast to the P. entomophila strain which led to high mortality even without Cry1C toxin co-ingestion. Thus, specific virulence factors adapted to the insect intestine might exist in B. thuringiensis/B. cereus and P. entomophila. This suggests a co-evolution between host and pathogens and supports the close links between B. thuringiensis and B. cereus and more distant links to their relative B. anthracis.     

Early Trichinella spiralis and Trichinella nativa infections induce similar gene expression profiles in rat jejunal mucosa

Trichinella spiralis causes a significantly higher parasite burden in rat muscle than Trichinella nativa. To assess whether the difference in infectivity is due to the early intestinal response, we analyzed gene expression changes in the rat jejunum during Trichinella infection with a whole-genome microarray. The rats were euthanized on day five of infection, and their jejunal mucosa was sampled for microarray analysis. In addition, intestinal histology and hematology were examined. Against our expectations, the gene expression changes were similar in both T.nativa- and T. spiralis-infected groups. The two groups were hence pooled, and in the combined Trichinella-infected group, 551 genes were overexpressed and 427 underexpressed when compared to controls (false discovery rate ?0.001 and fold change at least 2 in either direction). Pathway analysis identified seven pathways significantly associated with Trichinella infection (p < 0.05). The microarray data suggested nonspecific damage and an inflammatory response in the jejunal mucosa. Histological findings, including hyperemia, hemorrhage and a marked infiltration of inflammatory cells, supported the microarray data. Trichinella infection caused complex gene expression changes that indicate a host response to tissue damage in the mucosa of the jejunum, but the changes were not notably dependent on the studied species of Trichinella.     

VP2118 has major roles in Vibrio parahaemolyticus response to oxidative stress

Background

Reactive oxygen species (ROS), including superoxide anion radical, induce chronic risk of oxidative damage to many cellular macromolecules resulting in damage to cells. Superoxide dismutases (SODs) catalyze the dismutation of superoxide to oxygen and hydrogen peroxide and are a primary defense against ROS. Vibrio parahaemolyticus, a marine bacterium that causes acute gastroenteritis following consumption of raw or undercooked seafood, can survive ROS generated by intestinal inflammatory cells. However, there is little information concerning SODs in V. parahaemolyticus. This study aims to clarify the role of V. parahaemolyticus SODs against ROS.

Methods

V. parahaemolyticus SOD gene promoter activities were measured by a GFP reporter assay. Mutants of V. parahaemolyticus SOD genes were constructed and their SOD activity and resistance to oxidative stresses were measured.

Results

Bioinformatic analysis showed that V. parahaemolyticus SODs were distinguished by their metal cofactors, FeSOD (VP2118), MnSOD (VP2860), and CuZnSOD (VPA1514). VP2118 gene promoter activity was significantly higher than the other SOD genes. In a VP2118 gene deletion mutant, SOD activity was significantly decreased and could be recovered by VP2118 gene complementation. The absence of VP2118 resulted in significantly lowered resistance to ROS generated by hydrogen peroxide, hypoxanthine–xanthine oxidase, or Paraquat. Furthermore, both the N- and C-terminal SOD domains of VP2118 were necessary for ROS resistance.

Conclusion

VP2118 is the primary V. parahaemolyticus SOD and is vital for anti-oxidative stress responses.

General significance

The V. parahaemolyticus FeSOD VP2118 may enhance ROS resistance and could promote its survival in the intestinal tract to facilitate host tissue infection.     

AFN-1252 is a potent inhibitor of enoyl-ACP reductase from Burkholderia pseudomallei—Crystal structure,mode of action,and biological activity

Melioidosis is a tropical bacterial infection caused by Burkholderia pseudomallei (B. pseudomallei; Bpm), a Gram-negative bacterium. Current therapeutic options are largely limited to trimethoprim-sulfamethoxazole and β-lactam drugs, and the treatment duration is about 4 months. Moreover, resistance has been reported to these drugs. Hence, there is a pressing need to develop new antibiotics for Melioidosis. Inhibition of enoyl-ACP reducatase (FabI), a key enzyme in the fatty acid biosynthesis pathway has shown significant promise for antibacterial drug development. FabI has been identified as the major enoyl-ACP reductase present in B. pseudomallei. In this study, we evaluated AFN-1252, a Staphylococcus aureus FabI inhibitor currently in clinical development, for its potential to bind to BpmFabI enzyme and inhibit B. pseudomallei bacterial growth. AFN-1252 stabilized BpmFabI and inhibited the enzyme activity with an IC₅₀ of 9.6 nM. It showed good antibacterial activity against B. pseudomallei R15 strain, isolated from a melioidosis patient (MIC of 2.35 mg/L). X-ray structure of BpmFabI with AFN-1252 was determined at a resolution of 2.3 Å. Complex of BpmFabI with AFN-1252 formed a symmetrical tetrameric structure with one molecule of AFN-1252 bound to each monomeric subunit. The kinetic and thermal melting studies supported the finding that AFN-1252 can bind to BpmFabI independent of cofactor. The structural and mechanistic insights from these studies might help the rational design and development of new FabI inhibitors.