Genomic, molecular evolution, and expression analysis of NOX genes in soybean (Glycine max)

Reactive oxygen species (ROS) are versatile signaling molecules in sensing stresses and play critical roles in signaling and development. Plasma membrane NADPH oxidases (NOXs) are key producers of ROS, and play important roles in the regulation of plant-pathogen interactions. Here, we performed a comprehensive analysis of the NOX gene family in the soybean genome (Glycine max) and 17 NOX (GmNOX) genes were identified. Structural analysis revealed that the GmNOX proteins in soybean were as conserved as those in other plants. 8 duplicated gene pairs were formed by a Glycine-specific whole-genome duplication (WGD) event approximately 13 million years ago (Mya). The Ka/Ks ratios of GmNOX genes ranged from 0.04 to 0.28, suggesting that the GmNOX family had undergone purifying selection in soybean. Gene expression patterns showed different expression of these duplicate genes, suggesting that the GmNOXs were retained by substantial subfunctionalization during the soybean evolutionary processes. Subsequently, the expression of GmNOXs in response to drought and phytohormones were characterized via qPCR. Importantly, four GmNOXs showed strong expression in nodules, pointing to their probable involvement in nodulation. Thus, our results shed light on the evolutionary history of this family in soybean and contribute to the functional characterization of GmNOX genes in soybean.