Two Bradyrhizobium japonicum strains, TXVA and TXEA, were isolated for their desiccation tolerance and symbiotic performance with soybean as biofertilizers. Their genomes were sequenced and annotated using the Department of Energy Joint Genome Institute annotation pipeline. Sequencing yielded chromosomes of 9,193,770 and 9,339,455 bp for TXVA and TXEA, respectively.

Reference genes, also referred to as housekeeping genes (HKGs), play an important role in gene expression analysis by serving as an internal control. These HKGs are usually involved in basic cellular functions and their expression should remain at relatively constant levels. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) has been used to measure gene expression. Since the normalization of gene expression data depends on baseline expression of HKGs, it is important to identify and verify true HKGs for the qRT-PCR analysis. The goal of this study is to identify and confirm HKGs in Bradyrhizobium diazoefficiens, a nitrogen fixing bacterium which forms a symbiotic relationship with soybean. By revealing such HKGs, the normalization of gene expression would be more robust, reliable, and consistent. Here, we analyzed previous gene expression data for B. diazoefficiens under multiple environmental conditions. As a result, we identified seven constitutively expressed genes among 8453 genes across all conditions. Their fold-change values were within a range of −1.25-fold < x < 1.25-fold. We adopted GeNorm, NormFinder, and comparative ∆Ct methods to rank the seven candidate genes based on their expression stability. To validate these potential HKGs, we measured their expression in various experimental conditions, such as heat, pH, and heavy metal stress. The HKGs that were found in B. diazoefficiens were also applied in closely related species by identifying their homologs.