Search

Abstract Background Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches. Results We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome. Conclusions We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets.

Background: Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches. Results: We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome. Conclusions: We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets. [ABSTRACT FROM AUTHOR]

Abstract Background The plant growth promoting rhizobacterium Paenibacillus riograndensis SBR5 is a promising candidate to serve as crop inoculant. Despite its potential in providing environmental and economic benefits, the species P. riograndensis is poorly characterized. Here, we performed for the first time a detailed transcriptome analysis of P. riograndensis SBR5 using RNA-seq technology. Results RNA was isolated from P. riograndensis SBR5 cultivated under 15 different growth conditions and combined together in order to analyze an RNA pool representing a large set of expressed genes. The resultant total RNA was used to generate 2 different libraries, one enriched in 5′-ends of the primary transcripts and the other representing the whole transcriptome. Both libraries were sequenced and analyzed to identify the conserved sequences of ribosome biding sites and translation start motifs, and to elucidate operon structures present in the transcriptome of P. riograndensis. Sequence analysis of the library enriched in 5′-ends of the primary transcripts was used to identify 1082 transcription start sites (TSS) belonging to novel transcripts and allowed us to determine a promoter consensus sequence and regulatory sequences in 5′ untranslated regions including riboswitches. A putative thiamine pyrophosphate dependent riboswitch upstream of the thiamine biosynthesis gene thiC was characterized by translational fusion to a fluorescent reporter gene and shown to function in P. riograndensis SBR5. Conclusions Our RNA-seq analysis provides insight into the P. riograndensis SBR5 transcriptome at the systems level and will be a valuable basis for differential RNA-seq analysis of this bacterium.

Background: The plant growth promoting rhizobacterium Paenibacillus riograndensis SBR5 is a promising candidate to serve as crop inoculant. Despite its potential in providing environmental and economic benefits, the species P. riograndensis is poorly characterized. Here, we performed for the first time a detailed transcriptome analysis of P. riograndensis SBR5 using RNA-seq technology. Results: RNA was isolated from P. riograndensis SBR5 cultivated under 15 different growth conditions and combined together in order to analyze an RNA pool representing a large set of expressed genes. The resultant total RNA was used to generate 2 different libraries, one enriched in 5'-ends of the primary transcripts and the other representing the whole transcriptome. Both libraries were sequenced and analyzed to identify the conserved sequences of ribosome biding sites and translation start motifs, and to elucidate operon structures present in the transcriptome of P. riograndensis. Sequence analysis of the library enriched in 5'-ends of the primary transcripts was used to identify 1082 transcription start sites (TSS) belonging to novel transcripts and allowed us to determine a promoter consensus sequence and regulatory sequences in 5' untranslated regions including riboswitches. A putative thiamine pyrophosphate dependent riboswitch upstream of the thiamine biosynthesis gene thiC was characterized by translational fusion to a fluorescent reporter gene and shown to function in P. riograndensis SBR5. Conclusions: Our RNA-seq analysis provides insight into the P. riograndensis SBR5 transcriptome at the systems level and will be a valuable basis for differential RNA-seq analysis of this bacterium. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]