Genome-wide transcriptome analysis of the plant pathogen Xanthomonas identifies sRNAs with putative virulence functions

Genome-Wide Transcriptome Analysis of a Virulent sRNA, Trans217, in Xanthomonas oryzae pv. oryzae (Xoo), the Causative Agent of Rice Bacterial Blight

Small non-coding RNAs (sRNAs) act as post-transcriptional regulators to participate in many cellular processes. Among these, sRNA trans217 has been identified as a key virulent factor associated with pathogenicity in rice, triggering hypersensitive reactions in non-host tobacco and facilitating the secretion of the PthXo1 effector in Xanthomonas oryzae pv. oryzae (Xoo) strain PXO99A. Elucidating potential targets and downstream regulatory genes is crucial for understanding cellular networks governing pathogenicity and plant resistance. To explore the targets regulated by sRNA trans217, transcriptome sequencing was carried out to assess differential expression genes (DEGs) between the wild-type strain PXO99A and a mutant lacking the sRNA fragment under both virulence-inducing or normal growth conditions. DEG analysis revealed that sRNA trans217 was responsible for diverse functions, such as type III secretion system (T3SS), glutamate synthase activity, and oxidative stress response. Three genes were selected for further investigation due to their significant differential expression and biological relevance. Deletion of PXO_RS08490 attenuated the pathogenicity of Xoo in rice and reduced the tolerance level of PXO99A to hydrogen peroxide. These findings suggest a regulatory role of sRNA trans217 in modulating bacterial virulence through multiple gene targets, either directly or indirectly.

Genomic, Phylogenetic and Physiological Characterization of the PAH-Degrading Strain Gordonia polyisoprenivorans 135

The strain Gordonia polyisoprenivorans 135 is able to utilize a wide range of aromatic compounds. The aim of this work was to study the features of genetic organization and biotechnological potential of the strain G. polyisoprenivorans 135 as a degrader of aromatic compounds. The study of the genome of the strain 135 and the pangenome of the G. polyisoprenivorans species revealed that some genes, presumably involved in PAH catabolism, are atypical for Gordonia and belong to the pangenome of Actinobacteria. Analyzing the intergenic regions of strain 135 alongside the "panIGRome" of G. polyisoprenivorans showed that some intergenic regions in strain 135 also differ from those located between the same pairs of genes in related strains. The strain G. polyisoprenivorans 135 in our work utilized naphthalene (degradation degree 39.43%) and grew actively on salicylate. At present, this is the only known strain of G. polyisoprenivorans with experimentally confirmed ability to utilize these compounds.

Genome-Wide Computational Prediction and Analysis of Noncoding RNAs in Oleidesulfovibrio alaskensis G20

Noncoding RNAs (ncRNAs) play key roles in the regulation of important pathways, including cellular growth, stress management, signaling, and biofilm formation. Sulfate-reducing bacteria (SRB) contribute to huge economic losses causing microbial-induced corrosion through biofilms on metal surfaces. To effectively combat the challenges posed by SRB, it is essential to understand their molecular mechanisms of biofilm formation. This study aimed to identify ncRNAs in the genome of a model SRB, Oleidesulfovibrio alaskensis G20 (OA G20). Three in silico approaches revealed genome-wide distribution of 37 ncRNAs excluding tRNAs in the OA G20. These ncRNAs belonged to 18 different Rfam families. This study identified riboswitches, sRNAs, RNP, and SRP. The analysis revealed that these ncRNAs could play key roles in the regulation of several pathways of biosynthesis and transport involved in biofilm formation by OA G20. Three sRNAs, Pseudomonas P10, Hammerhead type II, and sX4, which were found in OA G20, are rare and their roles have not been determined in SRB. These results suggest that applying various computational methods could enrich the results and lead to the discovery of additional novel ncRNAs, which could lead to understanding the "rules of life of OA G20" during biofilm formation.

The Role of Iron in Phytopathogenic Microbe-Plant Interactions: Insights into Virulence and Host Immune Response

Iron is an essential element required for the growth and survival of nearly all forms of life. It serves as a catalytic component in multiple enzymatic reactions, such as photosynthesis, respiration, and DNA replication. However, the excessive accumulation of iron can result in cellular toxicity due to the production of reactive oxygen species (ROS) through the Fenton reaction. Therefore, to maintain iron homeostasis, organisms have developed a complex regulatory network at the molecular level. Besides catalyzing cellular redox reactions, iron also regulates virulence-associated functions in several microbial pathogens. Hosts and pathogens have evolved sophisticated strategies to compete against each other over iron resources. Although the role of iron in microbial pathogenesis in animals has been extensively studied, mechanistic insights into phytopathogenic microbe-plant associations remain poorly understood. Recent intensive research has provided intriguing insights into the role of iron in several plant-pathogen interactions. This review aims to describe the recent advances in understanding the role of iron in the lifestyle and virulence of phytopathogenic microbes, focusing on bacteria and host immune responses.

Small RNAs: Efficient and miraculous effectors that play key roles in plant-microbe interactions

Plants' response to pathogens is highly complex and involves changes at different levels, such as activation or repression of a vast array of genes. Recently, many studies have demonstrated that many RNAs, especially small RNAs (sRNAs), are involved in genetic expression and reprogramming affecting plant-pathogen interactions. The sRNAs, including short interfering RNAs and microRNAs, are noncoding RNA with 18-30 nucleotides, and are recognized as key genetic and epigenetic regulators. In this review, we summarize the new findings about defence-related sRNAs in the response to pathogens and our current understanding of their effects on plant-pathogen interactions. The main content of this review article includes the roles of sRNAs in plant-pathogen interactions, cross-kingdom sRNA trafficking between host and pathogen, and the application of RNA-based fungicides for plant disease control.

Mapping the Complex Transcriptional Landscape of the Phytopathogenic Bacterium Dickeya dadantii

Dickeya dadantii is a phytopathogenic bacterium that causes soft rot in a wide range of plant hosts worldwide and a model organism for studying virulence gene regulation. The present study provides a comprehensive and annotated transcriptomic map of D. dadantii obtained by a computational method combining five independent transcriptomic data sets: (i) paired-end RNA sequencing (RNA-seq) data for a precise reconstruction of the RNA landscape; (ii) DNA microarray data providing transcriptional responses to a broad variety of environmental conditions; (iii) long-read Nanopore native RNA-seq data for isoform-level transcriptome validation and determination of transcription termination sites; (iv) differential RNA sequencing (dRNA-seq) data for the precise mapping of transcription start sites; (v) in planta DNA microarray data for a comparison of gene expression profiles between in vitro experiments and the early stages of plant infection. Our results show that transcription units sometimes coincide with predicted operons but are generally longer, most of them comprising internal promoters and terminators that generate alternative transcripts of variable gene composition. We characterize the occurrence of transcriptional read-through at terminators, which might play a basal regulation role and explain the extent of transcription beyond the scale of operons. We finally highlight the presence of noncontiguous operons and excludons in the D. dadantii genome, novel genomic arrangements that might contribute to the basal coordination of transcription. The highlighted transcriptional organization may allow D. dadantii to finely adjust its gene expression program for a rapid adaptation to fast-changing environments. IMPORTANCE This is the first transcriptomic map of a Dickeya species. It may therefore significantly contribute to further progress in the field of phytopathogenicity. It is also one of the first reported applications of long-read Nanopore native RNA-seq in prokaryotes. Our findings yield insights into basal rules of coordination of transcription that might be valid for other bacteria and may raise interest in the field of microbiology in general. In particular, we demonstrate that gene expression is coordinated at the scale of transcription units rather than operons, which are larger functional genomic units capable of generating transcripts with variable gene composition for a fine-tuning of gene expression in response to environmental changes. In line with recent studies, our findings indicate that the canonical operon model is insufficient to explain the complexity of bacterial transcriptomes.

Small RNAs beyond Model Organisms: Have We Only Scratched the Surface?

Small RNAs (sRNAs) are essential regulators in the adaptation of bacteria to environmental changes and act by binding targeted mRNAs through base complementarity. Approximately 550 distinct families of sRNAs have been identified since their initial characterization in the 1980s, accelerated by the emergence of RNA-sequencing. Small RNAs are found in a wide range of bacterial phyla, but they are more prominent in highly researched model organisms compared to the rest of the sequenced bacteria. Indeed, Escherichia coli and Salmonella enterica contain the highest number of sRNAs, with 98 and 118, respectively, with Enterobacteriaceae encoding 145 distinct sRNAs, while other bacteria families have only seven sRNAs on average. Although the past years brought major advances in research on sRNAs, we have perhaps only scratched the surface, even more so considering RNA annotations trail behind gene annotations. A distinctive trend can be observed for genes, whereby their number increases with genome size, but this is not observable for RNAs, although they would be expected to follow the same trend. In this perspective, we aimed at establishing a more accurate representation of the occurrence of sRNAs in bacteria, emphasizing the potential for novel sRNA discoveries.

Xanthomonas oryzae pv. oryzae Exoribonuclease R Is Required for Complete Virulence in Rice, Optimal Motility, and Growth Under Stress

Exoribonuclease R (RNase R) is a 3' hydrolytic exoribonuclease that can degrade structured RNA. Mutation in RNase R affects virulence of certain human pathogenic bacteria. The aim of this study was to determine whether RNase R is necessary for virulence of the phytopathogen that causes bacterial blight in rice, Xanthomonas oryzae pv. oryzae (Xoo). In silico analysis has indicated that RNase R is highly conserved among various xanthomonads. Amino acid sequence alignment of Xoo RNase R with RNase R from various taxa indicated that Xoo RNase R clustered with RNase R of order Xanthomonadales. To study its role in virulence, we generated a gene disruption mutant of Xoo RNase R. The Xoo rnr^- mutant is moderately virulence deficient, and the complementing strain (rnr^-/pHM1::rnr) rescued the virulence deficiency of the mutant. We investigated swimming and swarming motilities in both nutrient-deficient minimal media and nutrient-optimal media. We observed that RNase R mutation has adversely affected the swimming and swarming motilities of Xoo in optimal media. However, in nutrient-deficient media only swimming motility was noticeably affected. Growth curves in optimal media at suboptimal temperature (15°C cold stress) indicate that the Xoo rnr^- mutant grows more slowly than the Xoo wild type and complementing strain (rnr^-/pHM1::rnr). Given these findings, we report for the first time that RNase R function is necessary for complete virulence of Xoo in rice. It is also important for motility of Xoo in media and for growth of Xoo at suboptimal temperature.

Xanthomonas transcriptome inside cauliflower hydathodes reveals bacterial virulence strategies and physiological adaptations at early infection stages

Xanthomonas campestris pv. campestris (Xcc) is a seed-transmitted vascular pathogen causing black rot disease on cultivated and wild Brassicaceae. Xcc enters the plant tissues preferentially via hydathodes, which are organs localized at leaf margins. To decipher both physiological and virulence strategies deployed by Xcc during early stages of infection, the transcriptomic profile of Xcc was analysed 3 days after entry into cauliflower hydathodes. Despite the absence of visible plant tissue alterations and despite a biotrophic lifestyle, 18% of Xcc genes were differentially expressed, including a striking repression of chemotaxis and motility functions. The Xcc full repertoire of virulence factors had not yet been activated but the expression of the HrpG regulon composed of 95 genes, including genes coding for the type III secretion machinery important for suppression of plant immunity, was induced. The expression of genes involved in metabolic adaptations such as catabolism of plant compounds, transport functions, sulphur and phosphate metabolism was upregulated while limited stress responses were observed 3 days postinfection. We confirmed experimentally that high-affinity phosphate transport is needed for bacterial fitness inside hydathodes. This analysis provides information about the nutritional and stress status of bacteria during the early biotrophic infection stages and helps to decipher the adaptive strategy of Xcc to the hydathode environment.

Comparative Genomics of Typical and Atypical Aeromonas salmonicida Complete Genomes Revealed New Insights into Pathogenesis Evolution

Aeromonas salmonicida is a global distributed Gram-negative teleost pathogen, affecting mainly salmonids in fresh and marine environments. A. salmonicida strains are classified as typical or atypical depending on their origin of isolation and phenotype. Five subspecies have been described, where A. salmonicida subsp. salmonicida is the only typical subspecies, and the subsp. achromogenes, masoucida, smithia, and pectinolytica are considered atypical. Genomic differences between A. salmonicida subsp. salmonicida isolates and their relationship with the current classification have not been explored. Here, we sequenced and compared the complete closed genomes of four virulent strains to elucidate their molecular diversity and pathogenic evolution using the more accurate genomic information so far. Phenotypes, biochemical, and enzymatic profiles were determined. PacBio and MiSeq sequencing platforms were utilized for genome sequencing. Comparative genomics showed that atypical strains belong to the subsp. salmonicida, with 99.55% ± 0.25% identity with each other, and are closely related to typical strains. The typical strain A. salmonicida J223 is closely related to typical strains, with 99.17% identity with the A. salmonicida A449. Genomic differences between atypical and typical strains are strictly related to insertion sequences (ISs) activity. The absence and presence of genes encoding for virulence factors, transcriptional regulators, and non-coding RNAs are the most significant differences between typical and atypical strains that affect their phenotypes. Plasmidome plays an important role in A. salmonicida virulence and genome plasticity. Here, we determined that typical strains harbor a larger number of plasmids and virulence-related genes that contribute to its acute virulence. In contrast, atypical strains harbor a single, large plasmid and a smaller number of virulence genes, reflected by their less acute virulence and chronic infection. The relationship between phenotype and A. salmonicida subspecies’ taxonomy is not evident. Comparative genomic analysis based on completed genomes revealed that the subspecies classification is more of a reflection of the ecological niche occupied by bacteria than their divergences at the genomic level except for their accessory genome.

A centenary for bacterial spot of tomato and pepper

DISEASE SYMPTOMS

Symptoms include water-soaked areas surrounded by chlorosis turning into necrotic spots on all aerial parts of plants. On tomato fruits, small, water-soaked, or slightly raised pale-green spots with greenish-white halos are formed, ultimately becoming dark brown and slightly sunken with a scabby or wart-like surface.

HOST RANGE

Main and economically important hosts include different types of tomatoes and peppers. Alternative solanaceous and nonsolanaceous hosts include Datura spp., Hyoscyamus spp., Lycium spp., Nicotiana rustica, Physalis spp., Solanum spp., Amaranthus lividus, Emilia fosbergii, Euphorbia heterophylla, Nicandra physaloides, Physalis pubescens, Sida glomerata, and Solanum americanum.

TAXONOMIC STATUS OF THE PATHOGEN

Domain, Bacteria; phylum, Proteobacteria; class, Gammaproteobacteria; order, Xanthomonadales; family, Xanthomonadaceae; genus, Xanthomonas; species, X. euvesicatoria, X. hortorum, X. vesicatoria.

SYNONYMS (NONPREFERRED SCIENTIFIC NAMES)

Bacterium exitiosum, Bacterium vesicatorium, Phytomonas exitiosa, Phytomonas vesicatoria, Pseudomonas exitiosa, Pseudomonas gardneri, Pseudomonas vesicatoria, Xanthomonas axonopodis pv. vesicatoria, Xanthomonas campestris pv. vesicatoria, Xanthomonas cynarae pv. gardneri, Xanthomonas gardneri, Xanthomonas perforans.

MICROBIOLOGICAL PROPERTIES

Colonies are gram-negative, oxidase-negative, and catalase-positive and have oxidative metabolism. Pale-yellow domed circular colonies of 1-2 mm in diameter grow on general culture media.

DISTRIBUTION

The bacteria are widespread in Africa, Brazil, Canada and the USA, Australia, eastern Europe, and south-east Asia. Occurrence in western Europe is restricted.

PHYTOSANITARY CATEGORIZATION

A2 no. 157, EU Annex designation II/A2.

EPPO CODES

XANTEU, XANTGA, XANTPF, XANTVE.

Genome-wide screen and functional analysis in Xanthomonas reveal a large number of mRNA-derived sRNAs, including the novel RsmA-sequester RsmU

Although bacterial small noncoding RNAs (sRNAs) are known to play a critical role in various cellular processes, including pathogenesis, the identity and action of such sRNAs are still poorly understood in many organisms. Here we have performed a genome-wide screen and functional analysis of the sRNAs in Xanthomonas campestris pv. campestris (Xcc), an important phytopathogen. The 50-500-nt RNA fragments isolated from the wild-type strain grown in a virulence gene-inducing condition were sequenced and a total of 612 sRNA candidates (SRCs) were identified. The majority (82%) of the SRCs were derived from mRNA, rather than specific sRNA genes. A representative panel of 121 SRCs were analysed by northern blotting; 117 SRCs were detected, supporting the contention that the overwhelming majority of the 612 SRCs identified are indeed sRNAs. Phenotypic analysis of strains overexpressing different candidates showed that a particular sRNA, RsmU, acts as a negative regulator of virulence, the hypersensitive response, and cell motility in Xcc. In vitro electrophoretic mobility shift assay and in vivo coimmunoprecipitation analyses indicated that RsmU interacted with the global posttranscriptional regulator RsmA, although sequence analysis displayed that RsmU is not a member of the sRNAs families known to antagonize RsmA. Northern blotting analyses demonstrated that RsmU has two isoforms that are processed from the 3'-untranslated region of the mRNA of XC1332 predicted to encode ComEA, a periplasmic protein required for DNA uptake in bacteria. This work uncovers an unexpected major sRNA biogenesis strategy in bacteria and a hidden layer of sRNA-mediated virulence regulation in Xcc.

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

Genome-Scale Metabolic Model of Xanthomonas phaseoli pv. manihotis: An Approach to Elucidate Pathogenicity at the Metabolic Level

Xanthomonas phaseoli pv. manihotis (Xpm) is the causal agent of cassava bacterial blight, the most important bacterial disease in this crop. There is a paucity of knowledge about the metabolism of Xanthomonas and its relevance in the pathogenic process, with the exception of the elucidation of the xanthan biosynthesis route. Here we report the reconstruction of the genome-scale model of Xpm metabolism and the insights it provides into plant-pathogen interactions. The model, iXpm1556, displayed 1,556 reactions, 1,527 compounds, and 890 genes. Metabolic maps of central amino acid and carbohydrate metabolism, as well as xanthan biosynthesis of Xpm, were reconstructed using Escher (https://escher.github.io/) to guide the curation process and for further analyses. The model was constrained using the RNA-seq data of a mutant of Xpm for quorum sensing (QS), and these data were used to construct context-specific models (CSMs) of the metabolism of the two strains (wild type and QS mutant). The CSMs and flux balance analysis were used to get insights into pathogenicity, xanthan biosynthesis, and QS mechanisms. Between the CSMs, 653 reactions were shared; unique reactions belong to purine, pyrimidine, and amino acid metabolism. Alternative objective functions were used to demonstrate a trade-off between xanthan biosynthesis and growth and the re-allocation of resources in the process of biosynthesis. Important features altered by QS included carbohydrate metabolism, NAD(P)+ balance, and fatty acid elongation. In this work, we modeled the xanthan biosynthesis and the QS process and their impact on the metabolism of the bacterium. This model will be useful for researchers studying host-pathogen interactions and will provide insights into the mechanisms of infection used by this and other Xanthomonas species.

Xanthomonas diversity, virulence and plant-pathogen interactions

Xanthomonas spp. encompass a wide range of plant pathogens that use numerous virulence factors for pathogenicity and fitness in plant hosts. In this Review, we examine recent insights into host-pathogen co-evolution, diversity in Xanthomonas populations and host specificity of Xanthomonas spp. that have substantially improved our fundamental understanding of pathogen biology. We emphasize the virulence factors in xanthomonads, such as type III secreted effectors including transcription activator-like effectors, type II secretion systems, diversity resulting in host specificity, evolution of emerging strains, activation of susceptibility genes and strategies of host evasion. We summarize the genomic diversity in several Xanthomonas spp. and implications for disease outbreaks, management strategies and breeding for disease resistance.

Structure of a Core Promoter in Bifidobacterium longum NCC2705

Bacterial promoters consist of core sequence motifs termed -35 and -10 boxes. The consensus motifs are TTGACA and TATAAT, respectively, which were identified from leading investigations on Escherichia coli However, the consensus sequences are not likely to fit genetically divergent bacteria. The sigma factor of the genus Bifidobacterium has a characteristic polar domain in the N terminus, suggesting the possibility of specific promoter recognition. We reevaluated the structure of Bifidobacterium longum NCC2705 promoters and compared them to other bacteria. Transcriptional start sites (TSSs) of the B. longum NCC2705 strain were identified using transcriptome sequencing (RNA-Seq) analysis to extract promoter regions. Conserved motifs of a bifidobacterial promoter were determined using regions upstream of TSSs and a hidden Markov model. As a result, consensus motifs of the -35 and -10 boxes were TTGTGC and TACAAT, respectively. To assess each base of both motifs, we constructed 37 plasmids based on pKO403-TPCTcon, including the hup promoter connected with a chloramphenicol acetyltransferase as a reporter gene. This reporter assay showed two optimal motifs of the -35 and -10 boxes, namely, TTGNNN and TANNNT, respectively. We further analyzed spacer lengths between the -35 and -10 boxes via a bioinformatics approach. The spacer lengths predominant in bacteria have been generally reported to be approximately 17 bp. In contrast, the predominant spacer lengths in the genus Bifidobacterium and related species were 11 bp, in addition to 17 bp. A reporter assay to assess the spacer lengths indicated that the 11-bp spacer length produced unusually high activity.IMPORTANCE The structures of sigma factors vary among bacterial strains, indicating that recognition rules may also vary. Therefore, we investigated the promoter structure of Bifidobacterium longum NCC2705 using a bioinformatics approach and wet analyses. The most frequent and optimal motifs were similar to other bacterial consensus motifs. The optimal spacer length between the two boxes was reported to be 17 bp. It is widely applied to a bioinformatics approach for other bacteria. Unexpectedly, conserved spacer lengths were 11 bp as well as 17 bp in the genus Bifidobacterium Moreover, the sigma factor of the genus Bifidobacterium has a characteristic domain in the N terminus which may contribute to the additional functions. Hence, it would be valuable to reevaluate the promoter in other organisms.

Deep genome annotation of the opportunistic human pathogen Streptococcus pneumoniae D39

A precise understanding of the genomic organization into transcriptional units and their regulation is essential for our comprehension of opportunistic human pathogens and how they cause disease. Using single-molecule real-time (PacBio) sequencing we unambiguously determined the genome sequence of Streptococcus pneumoniae strain D39 and revealed several inversions previously undetected by short-read sequencing. Significantly, a chromosomal inversion results in antigenic variation of PhtD, an important surface-exposed virulence factor. We generated a new genome annotation using automated tools, followed by manual curation, reflecting the current knowledge in the field. By combining sequence-driven terminator prediction, deep paired-end transcriptome sequencing and enrichment of primary transcripts by Cappable-Seq, we mapped 1015 transcriptional start sites and 748 termination sites. We show that the pneumococcal transcriptional landscape is complex and includes many secondary, antisense and internal promoters. Using this new genomic map, we identified several new small RNAs (sRNAs), RNA switches (including sixteen previously misidentified as sRNAs), and antisense RNAs. In total, we annotated 89 new protein-encoding genes, 34 sRNAs and 165 pseudogenes, bringing the S. pneumoniae D39 repertoire to 2146 genetic elements. We report operon structures and observed that 9% of operons are leaderless. The genome data are accessible in an online resource called PneumoBrowse (https://veeninglab.com/pneumobrowse) providing one of the most complete inventories of a bacterial genome to date. PneumoBrowse will accelerate pneumococcal research and the development of new prevention and treatment strategies.

Modular Cloning of the Type III Secretion Gene Cluster from the Plant-Pathogenic Bacterium Xanthomonas euvesicatoria

Type III secretion (T3S) systems are essential pathogenicity factors of most Gram-negative bacteria and translocate effector proteins into plant or animal cells. T3S systems can, therefore, be used as tools for protein delivery into eukaryotic cells, for instance after transfer of the T3S gene cluster into nonpathogenic recipient strains. Here, we report the modular cloning of the T3S gene cluster from the plant-pathogenic bacterium Xanthomonas euvesicatoria. The resulting multigene construct encoded a functional T3S system and delivered effector proteins into plant cells. The modular design of the T3S gene cluster allowed the efficient replacement and rearrangement of single genes or operons and the insertion of reporter genes for functional studies. In the present study, we used the modular T3S system to analyze the assembly of a fluorescent fusion of the predicted cytoplasmic ring protein HrcQ. Our studies demonstrate the use of the modular T3S gene cluster for functional analyses and mutant approaches in X. euvesicatoria. A potential application of the modular T3S system as protein delivery tool is discussed.

Bacterial Small RNAs in the Genus Herbaspirillum spp

The genus Herbaspirillum includes several strains isolated from different grasses. The identification of non-coding RNAs (ncRNAs) in the genus Herbaspirillum is an important stage studying the interaction of these molecules and the way they modulate physiological responses of different mechanisms, through RNA⁻RNA interaction or RNA⁻protein interaction. This interaction with their target occurs through the perfect pairing of short sequences (cis-encoded ncRNAs) or by the partial pairing of short sequences (trans-encoded ncRNAs). However, the companion Hfq can stabilize interactions in the trans-acting class. In addition, there are Riboswitches, located at the 5' end of mRNA and less often at the 3' end, which respond to environmental signals, high temperatures, or small binder molecules. Recently, CRISPR (clustered regularly interspaced palindromic repeats), in prokaryotes, have been described that consist of serial repeats of base sequences (spacer DNA) resulting from a previous exposure to exogenous plasmids or bacteriophages. We identified 285 ncRNAs in Herbaspirillum seropedicae (H. seropedicae) SmR1, expressed in different experimental conditions of RNA-seq material, classified as cis-encoded ncRNAs or trans-encoded ncRNAs and detected RNA riboswitch domains and CRISPR sequences. The results provide a better understanding of the participation of this type of RNA in the regulation of the metabolism of bacteria of the genus Herbaspirillum spp.

Two virulent sRNAs identified by genomic sequencing target the type III secretion system in rice bacterial blight pathogen

BACKGROUND

Small non-coding RNA (sRNA) short sequences regulate various biological processes in all organisms, including bacteria that are animal or plant pathogens. Virulent or pathogenicity-associated sRNAs have been increasingly elucidated in animal pathogens but little is known about similar category of sRNAs in plant-pathogenic bacteria. This is particularly true regarding rice bacterial blight pathogen Xanthomonas oryzae pathovar oryzae (Xoo) as studies on the virulent role of Xoo sRNAs is very limited at present.

RESULTS

The number and genomic distribution of sRNAs in Xoo were determined by bioinformatics analysis based on high throughput sequencing (sRNA-Seq) of the bacterial cultures from virulence-inducing and standard growth media, respectively. A total of 601 sRNAs were identified in the Xoo genome and ten virulent sRNA candidates were screened out based on significant differences of their expression levels between the culture conditions. In addition, trans3287 and trans3288 were also selected as candidates due to high expression levels in both media. The differential expression of 12 sRNAs evidenced by the sRNA-Seq data was confirmed by a convincing quantitative method. Based on genetic analysis of Xoo ΔsRNA mutants generated by deletion of the 12 single sRNAs, trans217 and trans3287 were characterized as virulent sRNAs. They are essential not only for the formation of bacterial blight in a susceptible rice variety Nipponbare but also for the induction of hypersensitive response (HR) in nonhost plant tobacco. Xoo Δtrans217 and Δtrans3287 mutants fail to induce bacterial blight in Nipponbare and also fail to induce the HR in tobacco, whereas, genetic complementation restores both mutants to the wild type in the virulent performance and HR induction. Similar effects of gene knockout and complementation were found in the expression of hrpG and hrpX genes, which encode regulatory proteins of the type III secretion system. Consistently, secretion of a type III effector, PthXo1, is blocked in Δtrans217 or Δtrans3287 bacterial cultures but retrieved by genetic complementation to both mutants.

CONCLUSIONS

The genetic analysis characterizes trans217 and trans3287 as pathogenicity-associated sRNAs essential for the bacterial virulence on the susceptible rice variety and for the HR elicitation in the nonhost plant. The molecular evidence suggests that both virulent sRNAs regulate the bacterial virulence by targeting the type III secretion system.

Genome-scale analysis of Methicillin-resistant Staphylococcus aureus USA300 reveals a tradeoff between pathogenesis and drug resistance

Staphylococcus aureus infection is a rising public health care threat. S. aureus is believed to have elaborate regulatory networks that orchestrate its virulence. Despite its importance, the systematic understanding of the transcriptional landscape of S. aureus is limited. Here, we describe the primary transcriptome landscape of an epidemic USA300 isolate of community-acquired methicillin-resistant S. aureus. We experimentally determined 1,861 transcription start sites with their principal promoter elements, including well-conserved -35 and -10 elements and weakly conserved -16 element and 5' untranslated regions containing AG-rich Shine-Dalgarno sequence. In addition, we identified 225 genes whose transcription was initiated from multiple transcription start sites, suggesting potential regulatory functions at transcription level. Along with the transcription unit architecture derived by integrating the primary transcriptome analysis with operon prediction, the measurement of differential gene expression revealed the regulatory framework of the virulence regulator Agr, the SarA-family transcriptional regulators, and β-lactam resistance regulators. Interestingly, we observed a complex interplay between virulence regulation, β-lactam resistance, and metabolism, suggesting a possible tradeoff between pathogenesis and drug resistance in the USA300 strain. Our results provide platform resource for the location of transcription initiation and an in-depth understanding of transcriptional regulation of pathogenesis, virulence, and antibiotic resistance in S. aureus.

Transcriptome analysis reveals the host selection fitness mechanisms of the Rhizoctonia solani AG1IA pathogen

Rhizoctonia solani AG1IA is a major generalist pathogen that causes sheath blight. Its genome, which was the first to be sequenced from the Rhizoctonia genus, may serve as a model for studying pathogenic mechanisms. To explore the pathogen-host fitness mechanism of sheath-blight fungus, a comprehensive comparative transcriptome ecotype analysis of R. solani AG1IA isolated from rice, soybean and corn during infection was performed. Special characteristics in gene expression, gene ontology terms and expression of pathogenesis-associated genes, including genes encoding secreted proteins, candidate effectors, hydrolases, and proteins involved in secondary metabolite production and the MAPK pathway, were revealed. Furthermore, as an important means of pathogenic modulation, diverse alternative splicing of key pathogenic genes in Rhizoctonia solani AG1IA during infections of the abovementioned hosts was uncovered for the first time. These important findings of key factors in the pathogenicity of R. solani AG1IA ecotypes during infection of various hosts explain host preference and provide novel insights into the pathogenic mechanisms and host-pathogen selection. Furthermore, they provide information on the fitness of Rhizoctonia, a severe pathogen with a wide host range.

Identification of new protein-coding genes with a potential role in the virulence of the plant pathogen Xanthomonas euvesicatoria

Background

Bacteria of the genus Xanthomonas are economically important plant pathogens. Pathogenicity of Xanthomonas spp. depends on the type III-secretion system and additional virulence determinants. The number of sequenced Xanthomonas genomes increases rapidly, however, accurate annotation of these genomes is difficult, because it relies on gene prediction programs. In this study, we used a mass-spectrometry (MS)-based approach to identify the proteome of Xanthomonas euvesicatoria (Xe) strain 85–10 also known as X. campestris pv. vesicatoria, a well-studied member of plant-pathogenic Xanthomonadaceae.

Results

Using different culture conditions, MS-datasets were searched against a six-frame-translated genome database of Xe. In total, we identified 2588 proteins covering 55% of the Xe genome, including 764 hitherto hypothetical proteins. Our proteogenomic approach identified 30 new protein-coding genes and allowed correction of the N-termini of 50 protein-coding genes. For five novel and two N-terminally corrected genes the corresponding proteins were confirmed by immunoblot. Furthermore, our data indicate that two putative type VI-secretion systems encoded in Xe play no role in bacterial virulence which was experimentally confirmed.

Conclusions

The discovery and re-annotation of numerous genes in the genome of Xe shows that also a well-annotated genome can be improved. Additionally, our proteogenomic analyses validates “hypothetical” proteins and will improve annotation of Xanthomonadaceae genomes, providing a solid basis for further studies.

Electronic supplementary material

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

Genome-Wide Detection of Small Regulatory RNAs in Deep-Sea Bacterium Shewanella piezotolerans WP3

Shewanella are one of the most abundant Proteobacteria in the deep-sea and are renowned for their versatile electron accepting capacities. The molecular mechanisms involved in their adaptation to diverse and extreme environments are not well understood. Small non-coding RNAs (sRNAs) are known for modulating the gene expression at transcriptional and posttranscriptional levels, subsequently playing a key role in microbial adaptation. To understand the potential roles of sRNAs in the adaptation of Shewanella toward deep-sea environments, here an in silico approach was utilized to detect the sRNAs in the genome of Shewanella piezotolerans WP3, a piezotolerant and psychrotolerant deep-sea iron reducing bacterium. After scanning 3673 sets of 5' and 3' UTRs of orthologous genes, 209 sRNA candidates were identified with high confidence in S. piezotolerans WP3. About 92% (193 out of 209) of these putative sRNAs belong to the class trans-encoded RNAs, suggesting that trans-regulatory RNAs are the dominant class of sRNAs in S. piezotolerans WP3. The remaining 16 cis-regulatory RNAs were validated through quantitative polymerase chain reaction. Five cis-sRNAs were further shown to act as cold regulated sRNAs. Our study provided additional evidence at the transcriptional level to decipher the microbial adaptation mechanisms to extreme environmental conditions.

Host-Pathogen interactions modulated by small RNAs

Biological processes such as defense mechanisms and microbial offense strategies are regulated through RNA induced interference in eukaryotes. Genetic mutations are modulated through biogenesis of small RNAs which directly impacts upon host development. Plant defense mechanisms are regulated and supported by a diversified group of small RNAs which are involved in streamlining several RNA interference pathways leading toward the initiation of pathogen gene silencing mechanisms. In the similar context, pathogens also utilize the support of small RNAs to launch their offensive attacks. Also there are strong evidences about the active involvement of these RNAs in symbiotic associations. Interestingly, small RNAs are not limited to the individuals in whom they are produced; they also show cross kingdom influences through variable interactions with other species thus leading toward the inter-organismic gene silencing. The phenomenon is understandable in the microbes which utilize these mechanisms to overcome host defense line. Understanding the mechanism of triggering host defense strategies can be a valuable step toward the generation of disease resistant host plants. We think that the cross kingdom trafficking of small RNA is an interesting insight that is needed to be explored for its vitality.

If it transcribes, we can sequence it: mining the complexities of host-pathogen-environment interactions using RNA-seq

Host-pathogen interactions are exceedingly complex because they involve multiple host tissues, often occur in the context of normal microflora, and can span diverse microenvironments. Although decades of gene expression studies have provided detailed insights into infection processes, technical challenges have restricted experiments to single pathogenic species or host tissues. RNA-sequencing (RNA-seq) has revolutionized the study of gene expression because in addition to quantifying transcriptional output, it allows detection and characterization of all transcripts in a genome. Here, we review how refined approaches to RNA-seq are used to map the transcriptional networks that control host-pathogen interactions. These enhanced techniques include dRNA-seq and term-seq for the fine-scale mapping of transcriptional start and termination sites, and dual RNA-seq for simultaneous sequencing of host and bacterial pathogen transcriptomes. Dual RNA-seq experiments are currently limited to in vitro infection systems that do not fully reflect the complexities of the in vivo environment, thus a challenge is to develop in vivo model systems and experimental approaches that address the biological heterogeneity of host environments, followed by the integration of RNA-seq with other genome-scale datasets to identify the transcriptional networks that mediate host-pathogen interactions.

Perspectives on the Transition From Bacterial Phytopathogen Genomics Studies to Applications Enhancing Disease Management: From Promise to Practice

The advent of genomics has advanced science into a new era, providing a plethora of "toys" for researchers in many related and disparate fields. Genomics has also spawned many new fields, including proteomics and metabolomics, furthering our ability to gain a more comprehensive view of individual organisms and of interacting organisms. Genomic information of both bacterial pathogens and their hosts has provided the critical starting point in understanding the molecular bases of how pathogens disrupt host cells to cause disease. In addition, knowledge of the complete genome sequence of the pathogen provides a potentially broad slate of targets for the development of novel virulence inhibitors that are desperately needed for disease management. Regarding plant bacterial pathogens and disease management, the potential for utilizing genomics resources in the development of durable resistance is enhanced because of developing technologies that enable targeted modification of the host. Here, we summarize the role of genomics studies in furthering efforts to manage bacterial plant diseases and highlight novel genomics-enabled strategies heading down this path.

Active Shiga-Like Toxin Produced by Some Aeromonas spp., Isolated in Mexico City

RNA silencing is a conserved mechanism that utilizes small RNAs (sRNAs) to direct the regulation of gene expression at the transcriptional or post-transcriptional level. Plants utilizing RNA silencing machinery to defend pathogen infection was first identified in plant–virus interaction and later was observed in distinct plant–pathogen interactions. RNA silencing is not only responsible for suppressing RNA accumulation and movement of virus and viroid, but also facilitates plant immune responses against bacterial, oomycete, and fungal infection. Interestingly, even the same plant sRNA can perform different roles when encounters with different pathogens. On the other side, pathogens counteract by generating sRNAs that directly regulate pathogen gene expression to increase virulence or target host genes to facilitate pathogen infection. Here, we summarize the current knowledge of the characterization and biogenesis of host- and pathogen-derived sRNAs, as well as the different RNA silencing machineries that plants utilize to defend against different pathogens. The functions of these sRNAs in defense and counter-defense and their mechanisms for regulation during different plant–pathogen interactions are also discussed.

Using Ecology, Physiology, and Genomics to Understand Host Specificity in Xanthomonas

How pathogens coevolve with and adapt to their hosts are critical to understanding how host jumps and/or acquisition of novel traits can lead to new disease emergences. The Xanthomonas genus includes Gram-negative plant-pathogenic bacteria that collectively infect a broad range of crops and wild plant species. However, individual Xanthomonas strains usually cause disease on only a few plant species and are highly adapted to their hosts, making them pertinent models to study host specificity. This review summarizes our current understanding of the molecular basis of host specificity in the Xanthomonas genus, with a particular focus on the ecology, physiology, and pathogenicity of the bacterium. Despite our limited understanding of the basis of host specificity, type III effectors, microbe-associated molecular patterns, lipopolysaccharides, transcriptional regulators, and chemotactic sensors emerge as key determinants for shaping host specificity.

Comparative transcriptomics across the prokaryotic tree of life

Whole-transcriptome sequencing studies from recent years revealed an unexpected complexity in transcriptomes of bacteria and archaea, including abundant non-coding RNAs, cis-antisense transcription and regulatory untranslated regions (UTRs). Understanding the functional relevance of the plethora of non-coding RNAs in a given organism is challenging, especially since some of these RNAs were attributed to ‘transcriptional noise’. To allow the search for conserved transcriptomic elements we produced comparative transcriptome maps for multiple species across the microbial tree of life. These transcriptome maps are detailed in annotations, comparable by gene families, and BLAST-searchable by user provided sequences. Our transcriptome collection includes 18 model organisms spanning 10 phyla/subphyla of bacteria and archaea that were sequenced using standardized RNA-seq methods. The utility of the comparative approach, as implemented in our web server, is demonstrated by highlighting genes with exceptionally long 5′UTRs across species, which correspond to many known riboswitches and further suggest novel putative regulatory elements. Our study provides a standardized reference transcriptome to major clinically and environmentally important microbial phyla. The viewer is available at http://exploration.weizmann.ac.il/TCOL, setting a framework for comparative studies of the microbial non-coding genome.

Genome-wide transcription start site mapping of Bradyrhizobium japonicum grown free-living or in symbiosis - a rich resource to identify new transcripts, proteins and to study gene regulation

BACKGROUND

Differential RNA-sequencing (dRNA-seq) is indispensable for determination of primary transcriptomes. However, using dRNA-seq data to map transcriptional start sites (TSSs) and promoters genome-wide is a bioinformatics challenge. We performed dRNA-seq of Bradyrhizobium japonicum USDA 110, the nitrogen-fixing symbiont of soybean, and developed algorithms to map TSSs and promoters.

RESULTS

A specialized machine learning procedure for TSS recognition allowed us to map 15,923 TSSs: 14,360 in free-living bacteria, 4329 in symbiosis with soybean and 2766 in both conditions. Further, we provide proteomic evidence for 4090 proteins, among them 107 proteins corresponding to new genes and 178 proteins with N-termini different from the existing annotation (72 and 109 of them with TSS support, respectively). Guided by proteomics evidence, previously identified TSSs and TSSs experimentally validated here, we assign a score threshold to flag 14 % of the mapped TSSs as a class of lower confidence. However, this class of lower confidence contains valid TSSs of low-abundant transcripts. Moreover, we developed a de novo algorithm to identify promoter motifs upstream of mapped TSSs, which is publicly available, and found motifs mainly used in symbiosis (similar to RpoN-dependent promoters) or under both conditions (similar to RpoD-dependent promoters). Mapped TSSs and putative promoters, proteomic evidence and updated gene annotation were combined into an annotation file.

CONCLUSIONS

The genome-wide TSS and promoter maps along with the extended genome annotation of B. japonicum represent a valuable resource for future systems biology studies and for detailed analyses of individual non-coding transcripts and ORFs. Our data will also provide new insights into bacterial gene regulation during the agriculturally important symbiosis between rhizobia and legumes.

Genome-wide transcriptomic analysis reveals correlation between higher WRKY61 expression and reduced symptom severity in Turnip crinkle virus infected Arabidopsis thaliana

Turnip crinkle virus (TCV) is a carmovirus that infects many Arabidopsis ecotypes. Most studies mainly focused on discovery of resistance genes against TCV infection, and there is no Next Generation Sequencing based comparative genome wide transcriptome analysis reported. In this study, RNA-seq based transcriptome analysis revealed that 238 (155 up-regulated and 83 down-regulated) significant differentially expressed genes with at least 15-fold change were determined. Fifteen genes (including upregulated, unchanged and downregulated) were selected for RNA-seq data validation using quantitative real-time PCR, which showed consistencies between these two sets of data. GO enrichment analysis showed that numerous terms such as stress, immunity, defence and chemical stimulus were affected in TCV-infected plants. One putative plant defence related gene named WRKY61 was selected for further investigation. It showed that WRKY61 overexpression plants displayed reduced symptoms and less virus accumulation, as compared to wild type (WT) and WRKY61 deficient lines, suggesting that higher WRKY61 expression level reduced TCV viral accumulation. In conclusion, our transcriptome analysis showed that global gene expression was detected in TCV-infected Arabidopsis thaliana. WRKY61 gene was shown to be negatively correlated with TCV infection and viral symptoms, which may be connected to plant immunity pathways.

Genome wide transcription start sites analysis of Xanthomonas campestris pv. campestris B100 with insights into the gum gene cluster directing the biosynthesis of the exopolysaccharide xanthan

Xanthomonas campestris pv. campestris (Xcc) is the major producer of the exopolysaccharide xanthan, the commercially most important natural polysaccharide of microbial origin. The current work provides deeper insights into the yet uncharacterized transcriptomic features of the xanthan producing strain Xcc-B100. Towards this goal, RNA sequencing of a library based on the selective enrichment of the 5' ends of native transcripts was performed. This approach resulted in the genome wide identification of 3067 transcription start sites (TSSs) that were further classified based on their genomic positions. Among them, 1545 mapped upstream of an actively transcribed CDS and 1363 were classified as novel TSSs representing antisense, internal, and TSSs belonging to previously unidentified genomic features. Analyzing the transcriptional strength of primary and antisense TSSs revealed that in some instances antisense transcription seemed to be initiated at a higher level than its sense counterpart. Mapping the exact positions of TSSs aided in the identification of promoter consensus motifs, ribosomal binding sites, and enhanced the genome annotation of 159 in silico predicted translational start (TLS) sites. The global view on length distribution of the 5' untranslated regions (5'-UTRs) deduced from the data pointed to the occurrence of leaderless transcripts and transcripts with unusually long 5'-UTRs, in addition to identifying seven putative riboswitch elements for Xcc-B100. Concerning the biosynthesis of xanthan, we focused on the transcriptional organization of the gum gene cluster. Under the conditions tested, we present evidence for a complex transcription pattern of the gum genes with multiple TSSs and an obvious considerable role of antisense transcription. The gene gumB, encoding an outer membrane xanthan exporter, is presented here as an example for genes that possessed a strong antisense TSS.

Novel Essential Role of Ethanol Oxidation Genes at Low Temperature Revealed by Transcriptome Analysis in the Antarctic Bacterium Pseudomonas extremaustralis

Temperature is one of the most important factors for bacterial growth and development. Cold environments are widely distributed on earth, and psychrotolerant and psychrophilic microorganisms have developed different adaptation strategies to cope with the stress derived from low temperatures. Pseudomonas extremaustralis is an Antarctic bacterium able to grow under low temperatures and to produce high amounts of polyhydroxyalkanoates (PHAs). In this work, we analyzed the genome-wide transcriptome by RNA deep-sequencing technology of early exponential cultures of P. extremaustralis growing in LB (Luria Broth) supplemented with sodium octanoate to favor PHA accumulation at 8°C and 30°C. We found that genes involved in primary metabolism, including tricarboxylic acid cycle (TCA) related genes, as well as cytochromes and amino acid metabolism coding genes, were repressed at low temperature. Among up-regulated genes, those coding for transcriptional regulatory and signal transduction proteins were over-represented at cold conditions. Remarkably, we found that genes involved in ethanol oxidation, exaA, exaB and exaC, encoding a pyrroloquinoline quinone (PQQ)-dependent ethanol dehydrogenase, the cytochrome c550 and an aldehyde dehydrogenase respectively, were up-regulated. Along with RNA-seq experiments, analysis of mutant strains for pqqB (PQQ biosynthesis protein B) and exaA were carried out. We found that the exaA and pqqB genes are essential for growth under low temperature in LB supplemented with sodium octanoate. Additionally, p-rosaniline assay measurements showed the presence of alcohol dehydrogenase activity at both 8°C and 30°C, while the activity was abolished in a pqqB mutant strain. These results together with the detection of ethanol by gas chromatography in P. extremaustralis cultures grown at 8°C support the conclusion that this pathway is important under cold conditions. The obtained results have led to the identification of novel components involved in cold adaptation mechanisms in this bacterium, suggesting for the first time a role of the ethanol oxidation pathway for bacterial growth at low temperatures.

Bacterial spot of tomato and pepper: diverse Xanthomonas species with a wide variety of virulence factors posing a worldwide challenge

TAXONOMIC STATUS

Bacteria; Phylum Proteobacteria; Class Gammaproteobacteria; Order Xanthomonadales; Family Xanthomonadaceae; Genus Xanthomonas; Species Xanthomonas euvesicatoria, Xanthomonas vesicatoria, Xanthomonas perforans and Xanthomonas gardneri.

MICROBIOLOGICAL PROPERTIES

Gram-negative, rod-shaped bacterium, aerobic, motile, single polar flagellum.

HOST RANGE

Causes bacterial spot disease on plants belonging to the Solanaceae family, primarily tomato (Solanum lycopersicum), pepper (Capsicum annuum) and chilli peppers (Capsicum frutescens).

DISEASE SYMPTOMS

Necrotic lesions on all above-ground plant parts.

DISTRIBUTION

Worldwide distribution of X. euvesicatoria and X. vesicatoria on tomato and pepper; X. perforans and X. gardneri increasingly being isolated from the USA, Canada, South America, Africa and Europe. A wide diversity within the bacterial spot disease complex, with an ability to cause disease at different temperatures, makes this pathogen group a worldwide threat to tomato and pepper production. Recent advances in genome analyses have revealed the evolution of the pathogen with a plethora of novel virulence factors. Current management strategies rely on the use of various chemical control strategies and sanitary measures to minimize pathogen spread through contaminated seed. Chemical control strategies have been a challenge because of resistance by the pathogen. Breeding programmes have been successful in developing commercial lines with hypersensitive and quantitative resistance. However, durability of resistance has been elusive. Recently, a transgenic approach has resulted in the development of tomato genotypes with significant levels of resistance and improved yield that hold promise. In this article, we discuss the current taxonomic status, distribution of the four species, knowledge of virulence factors, detection methods and strategies for disease control with possible directions for future research.

Genomics and transcriptomics of Xanthomonas campestris species challenge the concept of core type III effectome

BACKGROUND

The bacterial species Xanthomonas campestris infects a wide range of Brassicaceae. Specific pathovars of this species cause black rot (pv. campestris), bacterial blight of stock (pv. incanae) or bacterial leaf spot (pv. raphani).

RESULTS

In this study, we extended the genomic coverage of the species by sequencing and annotating the genomes of strains from pathovar incanae (CFBP 1606R and CFBP 2527R), pathovar raphani (CFBP 5828R) and a pathovar formerly named barbareae (CFBP 5825R). While comparative analyses identified a large core ORFeome at the species level, the core type III effectome was limited to only three putative type III effectors (XopP, XopF1 and XopAL1). In Xanthomonas, these effector proteins are injected inside the plant cells by the type III secretion system and contribute collectively to virulence. A deep and strand-specific RNA sequencing strategy was adopted in order to experimentally refine genome annotation for strain CFBP 5828R. This approach also allowed the experimental definition of novel ORFs and non-coding RNA transcripts. Using a constitutively active allele of hrpG, a master regulator of the type III secretion system, a HrpG-dependent regulon of 141 genes co-regulated with the type III secretion system was identified. Importantly, all these genes but seven are positively regulated by HrpG and 56 of those encode components of the Hrp type III secretion system and putative effector proteins.

CONCLUSIONS

This dataset is an important resource to mine for novel type III effector proteins as well as for bacterial genes which could contribute to pathogenicity of X. campestris.

dRNA-Seq Reveals Genomewide TSSs and Noncoding RNAs of Plant Beneficial Rhizobacterium Bacillus amyloliquefaciens FZB42

Bacillus amyloliquefaciens subsp. plantarum FZB42 is a representative of Gram-positive plant-growth-promoting rhizobacteria (PGPR) that inhabit plant root environments. In order to better understand the molecular mechanisms of bacteria-plant symbiosis, we have systematically analyzed the primary transcriptome of strain FZB42 grown under rhizosphere-mimicking conditions using differential RNA sequencing (dRNA-seq). Our analysis revealed 4,877 transcription start sites for protein-coding genes, identified genes differentially expressed under different growth conditions, and corrected many previously mis-annotated genes. We also identified a large number of riboswitches and cis-encoded antisense RNAs, as well as trans-encoded small noncoding RNAs that may play important roles in the gene regulation of Bacillus. Overall, our analyses provided a landscape of Bacillus primary transcriptome and improved the knowledge of rhizobacteria-host interactions.

Leaderless Transcripts and Small Proteins Are Common Features of the Mycobacterial Translational Landscape

RNA-seq technologies have provided significant insight into the transcription networks of mycobacteria. However, such studies provide no definitive information on the translational landscape. Here, we use a combination of high-throughput transcriptome and proteome-profiling approaches to more rigorously understand protein expression in two mycobacterial species. RNA-seq and ribosome profiling in Mycobacterium smegmatis, and transcription start site (TSS) mapping and N-terminal peptide mass spectrometry in Mycobacterium tuberculosis, provide complementary, empirical datasets to examine the congruence of transcription and translation in the Mycobacterium genus. We find that nearly one-quarter of mycobacterial transcripts are leaderless, lacking a 5’ untranslated region (UTR) and Shine-Dalgarno ribosome-binding site. Our data indicate that leaderless translation is a major feature of mycobacterial genomes and is comparably robust to leadered initiation. Using translational reporters to systematically probe the cis-sequence requirements of leaderless translation initiation in mycobacteria, we find that an ATG or GTG at the mRNA 5’ end is both necessary and sufficient. This criterion, together with our ribosome occupancy data, suggests that mycobacteria encode hundreds of small, unannotated proteins at the 5’ ends of transcripts. The conservation of small proteins in both mycobacterial species tested suggests that some play important roles in mycobacterial physiology. Our translational-reporter system further indicates that mycobacterial leadered translation initiation requires a Shine Dalgarno site in the 5’ UTR and that ATG, GTG, TTG, and ATT codons can robustly initiate translation. Our combined approaches provide the first comprehensive view of mycobacterial gene structures and their non-canonical mechanisms of protein expression.

The current paradigm for bacterial translation is based on an mRNA that includes an untranslated leader sequence containing the ribosome-binding site upstream of the initiation codon. We applied genome-scale approaches to map the protein-coding regions in the genomes of Mycobacterium smegmatis and Mycobacterium tuberculosis. We found that nearly one-quarter of mycobacterial transcripts are leaderless in mycobacterial species, thus indicating that ribosomes must recognize these mRNAs by a novel mechanism and suggesting that there are alternative modes of bacterial translation beyond the Escherichia coli paradigm. Our translational profiling showed that many mycobacterial proteins are mis-annotated, and also found many new genes encoding small proteins that had been previously overlooked, which are likely to play novel roles in diverse cellular processes. We also developed a new reporter system that provides mechanistic insights into translation initiation through deep sequencing. Our data show that leaderless translation is a robust process that is conserved in mycobacteria, that leaderless translation only requires that the mRNA begin with a start codon, and predict that mycobacteria encode hundreds of small proteins. This work will help us understand gene structure, genome organization and protein expression in bacteria, and how the translational machinery differs in different organisms.

The complexity of bacterial transcriptomes

For eukaryotes there seems to be no doubt that differences on the trancriptomic level substantially contribute to the process of species diversification, whereas for bacteria this is thought to be less important. Recent years saw a significant increase in full transcriptome studies for bacteria, which provided deep insight into the architecture of bacterial transcriptomes. Most notably, it became evident that, in contrast to previous scientific consensus, bacterial transcriptomes are quite complex. There exist a large number of cis-antisense RNAs, non-coding RNAs, overlapping transcripts and RNA elements that regulate transcription, such as riboswitches. Furthermore, processing and degradation of RNA has gained interest, because it has a significant impact on the composition of the transcriptome. In this review, we summarize recent findings and put them into a broader context with respect to the complexity of bacterial transcriptomes and its putative biological meanings.

Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors

Background

The fungus Cronartium ribicola (Cri) is an economically and ecologically important forest pathogen that causes white pine blister rust (WPBR) disease on five-needle pines. To cause stem cankers and kill white pine trees the fungus elaborates a life cycle with five stages of spore development on five-needle pines and the alternate host Ribes plants. To increase our understanding of molecular WP-BR interactions, here we report genome-wide transcriptional profile analysis of C. ribicola using RNA-seq.

Results

cDNA libraries were constructed from aeciospore, urediniospore, and western white pine (Pinus monticola) tissues post Cri infection. Over 200 million RNA-seq 100-bp paired-end (PE) reads from rust fungal spores were de novo assembled and a reference transcriptome was generated with 17,880 transcripts that were expressed from 13,629 unigenes. A total of 734 unique proteins were predicted as a part of the Cri secretome from complete open reading frames (ORFs), and 41 % of them were Cronartium-specific. This study further identified a repertoire of candidate effectors and other pathogenicity determinants. Differentially expressed genes (DEGs) were identified to gain an understanding of molecular events important during the WPBR fungus life cycle by comparing Cri transcriptomes at different infection stages. Large-scale changes of in planta gene expression profiles were observed, revealing that multiple fungal biosynthetic pathways were enhanced during mycelium growth inside infected pine stem tissues. Conversely, many fungal genes that were up-regulated at the urediniospore stage appeared to be signalling components and transporters. The secreted fungal protein genes that were up-regulated in pine needle tissues during early infection were primarily associated with cell wall modifications, possibly to mask the rust pathogen from plant defenses.

Conclusion

This comprehensive transcriptome profiling substantially improves our current understanding of molecular WP-BR interactions. The repertoire of candidate effectors and other putative pathogenicity determinants identified here are valuable for future functional analysis of Cri virulence and pathogenicity.

Electronic supplementary material

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

Use of siRNA molecular beacons to detect and attenuate mycobacterial infection in macrophages

Tuberculosis is one of the leading infectious diseases plaguing mankind and is mediated by the facultative pathogen, Mycobacterium tuberculosis (MTB). Once the pathogen enters the body, it subverts the host immune defenses and thrives for extended periods of time within the host macrophages in the lung granulomas, a condition called latent tuberculosis (LTB). Persons with LTB are prone to reactivation of the disease when the body’s immunity is compromised. Currently there are no reliable and effective diagnosis and treatment options for LTB, which necessitates new research in this area. The mycobacterial proteins and genes mediating the adaptive responses inside the macrophage is largely yet to be determined. Recently, it has been shown that the mce operon genes are critical for host cell invasion by the mycobacterium and for establishing a persistent infection in both in vitro and in mouse models of tuberculosis. The YrbE and Mce proteins which are encoded by the MTB mce operons display high degrees of homology to the permeases and the surface binding protein of the ABC transports, respectively. Similarities in structure and cell surface location impute a role in cell invasion at cholesterol rich regions and immunomodulation. The mce4 operon is also thought to encode a cholesterol transport system that enables the mycobacterium to derive both energy and carbon from the host membrane lipids and possibly generating virulence mediating metabolites, thus enabling the bacteria in its long term survival within the granuloma. Various deletion mutation studies involving individual or whole mce operon genes have shown to be conferring varying degrees of attenuation of infectivity or at times hypervirulence to the host MTB, with the deletion of mce4A operon gene conferring the greatest degree of attenuation of virulence. Antisense technology using synthetic siRNAs has been used in knocking down genes in bacteria and over the years this has evolved into a powerful tool for elucidating the roles of various genes mediating infectivity and survival in mycobacteria. Molecular beacons are a newer class of antisense RNA tagged with a fluorophore/quencher pair and their use for in vivo detection and knockdown of mRNA is rapidly gaining popularity.

Small RNAs--the secret agents in the plant-pathogen interactions

Eukaryotic regulatory small RNAs (sRNAs) that induce RNA interference (RNAi) are involved in a plethora of biological processes, including host immunity and pathogen virulence. In plants, diverse classes of sRNAs contribute to the regulation of host innate immunity. These immune-regulatory sRNAs operate through distinct RNAi pathways that trigger transcriptional or post-transcriptional gene silencing. Similarly, many pathogen-derived sRNAs also regulate pathogen virulence. Remarkably, the influence of regulatory sRNAs is not limited to the individual organism in which they are generated. It can sometimes extend to interacting species from even different kingdoms. There they trigger gene silencing in the interacting organism, a phenomenon called cross-kingdom RNAi. This is exhibited in advanced pathogens and parasites that produce sRNAs to suppress host immunity. Conversely, in host-induced gene silencing (HIGS), diverse plants are engineered to trigger RNAi against pathogens and pests to confer host resistance. Cross-kingdom RNAi opens up a vastly unexplored area of research on mobile sRNAs in the battlefield between hosts and pathogens.

Small RNA Deep-Sequencing Analyses Reveal a New Regulator of Virulence in Agrobacterium fabrum C58

Novel ways of regulating Ti plasmid functions were investigated by studying small RNAs (sRNAs) that are known to act as posttranscriptional regulators in plant pathogenic bacteria. sRNA-seq analyses of Agrobacterium fabrum C58 allowed us to identify 1,108 small transcripts expressed in several growth conditions that could be sRNAs. A quarter of them were confirmed by bioinformatics or by biological experiments. Antisense RNAs represent 24% of the candidates and they are over-represented on the pTi (with 62% of pTi sRNAs), suggesting differences in the regulatory mechanisms between the essential and accessory replicons. Moreover, a large number of these pTi antisense RNAs are transcribed opposite to those genes involved in virulence. Others are 5'- and 3'-untranslated region RNAs and trans-encoded RNAs. We have validated, by rapid amplification of cDNA ends polymerase chain reaction, the transcription of 14 trans-encoded RNAs, among which RNA1111 is expressed from the pTiC58. Its deletion decreased the aggressiveness of A. fabrum C58 on tomatoes, tobaccos, and kalanchoe, suggesting that this sRNA activates virulence. The identification of its putative target mRNAs (6b gene, virC2, virD3, and traA) suggests that this sRNA may coordinate two of the major pTi functions, the infection of plants and its dissemination among bacteria.

Transcriptomic profiling of Yersinia pseudotuberculosis reveals reprogramming of the Crp regulon by temperature and uncovers Crp as a master regulator of small RNAs

One hallmark of pathogenic yersiniae is their ability to rapidly adjust their life-style and pathogenesis upon host entry. In order to capture the range, magnitude and complexity of the underlying gene control mechanisms we used comparative RNA-seq-based transcriptomic profiling of the enteric pathogen Y. pseudotuberculosis under environmental and infection-relevant conditions. We identified 1151 individual transcription start sites, multiple riboswitch-like RNA elements, and a global set of antisense RNAs and previously unrecognized trans-acting RNAs. Taking advantage of these data, we revealed a temperature-induced and growth phase-dependent reprogramming of a large set of catabolic/energy production genes and uncovered the existence of a thermo-regulated ‘acetate switch’, which appear to prime the bacteria for growth in the digestive tract. To elucidate the regulatory architecture linking nutritional status to virulence we also refined the CRP regulon. We identified a massive remodelling of the CRP-controlled network in response to temperature and discovered CRP as a transcriptional master regulator of numerous conserved and newly identified non-coding RNAs which participate in this process. This finding highlights a novel level of complexity of the regulatory network in which the concerted action of transcriptional regulators and multiple non-coding RNAs under control of CRP adjusts the control of Yersinia fitness and virulence to the requirements of their environmental and virulent life-styles.

Many bacterial pathogens cycle between environmental sources and mammalian hosts. Adaptation to the different natural habitats and host niches is achieved through complex regulatory networks which adjust synthesis of the large repertoire of crucial virulence factors and fitness determinants. To uncover underlying control circuits, we determined the first in-depth single-nucleotide resolution transcriptome of Yersinia. This revealed important novel genetic information, such as global locations of transcriptional start sites, non-coding RNAs, potential riboswitches and provided a set of virulence-relevant expression profiles, which constitute a valuable tool for the research community. The analysis further uncovered a temperature-induced global reprogramming of central metabolic functions, likely to support intestinal colonization of the pathogen. This is accompanied by a major reorganization of the CRP regulon, which involves a multitude of regulatory RNAs. The primary consequence is a fine-tuned, coordinated control of metabolism and virulence through a plethora of environmentally controlled regulatory RNAs allowing rapid adaptation and high flexibility during life-style changes.

Transcriptome analysis of thermophilic methylotrophic Bacillus methanolicus MGA3 using RNA-sequencing provides detailed insights into its previously uncharted transcriptional landscape

Background

Bacillus methanolicus MGA3 is a thermophilic, facultative ribulose monophosphate (RuMP) cycle methylotroph. Together with its ability to produce high yields of amino acids, the relevance of this microorganism as a promising candidate for biotechnological applications is evident. The B. methanolicus MGA3 genome consists of a 3,337,035 nucleotides (nt) circular chromosome, the 19,174 nt plasmid pBM19 and the 68,999 nt plasmid pBM69. 3,218 protein-coding regions were annotated on the chromosome, 22 on pBM19 and 82 on pBM69. In the present study, the RNA-seq approach was used to comprehensively investigate the transcriptome of B. methanolicus MGA3 in order to improve the genome annotation, identify novel transcripts, analyze conserved sequence motifs involved in gene expression and reveal operon structures. For this aim, two different cDNA library preparation methods were applied: one which allows characterization of the whole transcriptome and another which includes enrichment of primary transcript 5′-ends.

Results

Analysis of the primary transcriptome data enabled the detection of 2,167 putative transcription start sites (TSSs) which were categorized into 1,642 TSSs located in the upstream region (5′-UTR) of known protein-coding genes and 525 TSSs of novel antisense, intragenic, or intergenic transcripts. Firstly, 14 wrongly annotated translation start sites (TLSs) were corrected based on primary transcriptome data. Further investigation of the identified 5′-UTRs resulted in the detailed characterization of their length distribution and the detection of 75 hitherto unknown cis-regulatory RNA elements. Moreover, the exact TSSs positions were utilized to define conserved sequence motifs for translation start sites, ribosome binding sites and promoters in B. methanolicus MGA3. Based on the whole transcriptome data set, novel transcripts, operon structures and mRNA abundances were determined. The analysis of the operon structures revealed that almost half of the genes are transcribed monocistronically (940), whereas 1,164 genes are organized in 381 operons. Several of the genes related to methylotrophy had highly abundant transcripts.

Conclusion

The extensive insights into the transcriptional landscape of B. methanolicus MGA3, gained in this study, represent a valuable foundation for further comparative quantitative transcriptome analyses and possibly also for the development of molecular biology tools which at present are very limited for this organism.

Electronic supplementary material

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