Reference genes for RT-qPCR studies in Corynebacterium pseudotuberculosis identified through analysis of RNA-seq data

Reference genes presenting stable expression profiles over a wide variety of conditions are required in relative expression studies of specific bacterial genes by quantitative reverse transcription PCR (RT-qPCR). High-throughput sequencing of bacterial transcriptomes using the RNA-seq methodology now provides a wealth of data that may be searched for identification of the most stably expressed genes of a given bacterium. Herein, we searched a RNA-seq dataset from various experiments with the pathogenic bacterium Corynebacterium pseudotuberculosis, grown under different stress conditions, in order to select appropriate candidate reference genes for this species. Nineteen genes involved in maintenance of basic cellular functions, so-called housekeeping genes, were chosen for study and their expression profiles in C. pseudotuberculosis were evaluated throughout all growth conditions. Eight of these genes (atpA, dnaG, efp, fusA, gyrA, gyrB, rpoB, and rpoC), mostly participating in DNA replication and transcription, matched the defined criteria to be included as candidate reference genes. Transcriptional levels of these genes were quantified by RT-qPCR assays after growth of C. pseudotuberculosis under two additional conditions. Expression stability analysis by NormFinder indicated the combination of genes encoding DNA gyrase subunit A (gyrA) and elongation factor P (fusA) as the most suitable for normalization of RT-qPCR studies in C. pseudotuberculosis.