Transcriptome analysis of Acidovorax avenaesubsp. avenae cultivated in vivo and co-culture with Burkholderia seminalis

Response of bacterial pathogen to environmental bacteria and its host is critical for understanding of microbial adaption and pathogenesis. Here, we used RNA-Seq to comprehensively and quantitatively assess the transcriptional response of Acidovorax avenae subsp. avenae strain RS-1 cultivated in vitro, in vivo and in co-culture with rice rhizobacterium Burkholderia seminalis R456. Results revealed a slight response to other bacteria, but a strong response to host. In particular, a large number of virulence associated genes encoding Type I to VI secretion systems, 118 putative non-coding RNAs, and 7 genomic islands (GIs) were differentially expressed in vivo based on comparative genomic and transcriptomic analyses. Furthermore, the loss of virulence for knockout mutants of 11 differentially expressed T6SS genes emphasized the importance of these genes in bacterial pathogenicity. In addition, the reliability of expression data obtained by RNA-Seq was supported by quantitative real-time PCR of the 25 selected T6SS genes. Overall, this study highlighted the role of differentially expressed genes in elucidating bacterial pathogenesis based on combined analysis of RNA-Seq data and knockout of T6SS genes.