All lateral organ development in plants, such as nodulation in legumes, requires the temporal and spatial regulation of genes and gene networks. A total mRNA profiling approach using RNA-seq to target the specific soybean (Glycine max) root tissues responding to compatible rhizobia [i.e. the Zone Of Nodulation (ZON)] revealed a large number of novel, often transient, mRNA changes occurring during the early stages of nodulation. Focusing on the ZON enabled us to discard the majority of root tissues and their developmentally diverse gene transcripts, thereby highlighting the lowly and transiently expressed nodulation-specific genes. It also enabled us to concentrate on a precise moment in early nodule development at each sampling time. We focused on discovering genes regulated specifically by the Bradyrhizobium-produced Nod factor signal, by inoculating roots with either a competent wild-type or incompetent mutant (nodC-) strain of Bradyrhizobium japonicum. Collectively, 2915 genes were identified as being differentially expressed, including many known soybean nodulation genes. A number of unknown nodulation gene candidates and soybean orthologues of nodulation genes previously reported in other legume species were also identified. The differential expression of several candidates was confirmed and further characterized via inoculation time-course studies and qRT-PCR. The expression of many genes, including an endo-1,4-β-glucanase, a cytochrome P450 and a TIR-LRR-NBS receptor kinase, was transient, peaking quickly during the initiation of nodule ontogeny. Additional genes were found to be down-regulated. Significantly, a set of differentially regulated genes acting in the gibberellic acid (GA) biosynthesis pathway was discovered, suggesting a novel role of GAs in nodulation.