RNA-Seq analysis of the T3SA regulon in Shigella flexneri reveals two new chromosomal genes upregulated in the on-state

Shigella spp. are enterobacteria that invade human colonic mucosal cells using their Type Three Secretion Apparatus (T3SA). Shigella spp. possess a large plasmid that encodes most of its virulence factors and has been the focus of seminal work that defined the T3SA regulon. Thus, a global assessment of the transcriptional response regulated by the T3SA has been lacking. Herein we used RNA-Seq to identify genes that are differentially expressed when the T3SA is active (on-state) versus inactive (off-state). The quality of the RNA-Seq dataset was validated by its correlation with a prior microarray study. Using novel insights about the expression of non-coding regions, bioinformatic tools and experimentations, we demonstrated the existence of six operons and evidence that ipaH2.5 is a pseudogene. In addition, 86 chromosomal genes were downregulated in the on-state including several non-coding transcripts corresponding to short antisense RNA embedded in the 16S and 23S RNA genes, and 40 coding transcripts, whose cognate proteins were highly connected at the genetic and bio-chemical levels. Finally, we identified two novel chromosomal genes dubbed gem1 and gem3, which were upregulated in the on-state similarly to genes belonging to the T3SA regulon. The latter findings were validated on biological triplicates by droplet digital PCR. To our knowledge gem1 and gem3 are the first chromosomal members of the T3SA regulon that have no homologs on the plasmid. Our approach provides a path to optimizing RNA-Seq studies in case of bacterial models that had previously been the subject of medium to large scale studies.