A global coexpression network of soybean genes gives insight into the evolution of nodulation in non‐legumes and legumes

A coexpression network is a powerful tool for revealing genes’ relationship with many biological processes. Mass transcriptomic and genomic data from different plant species provide the foundation for understanding the evolution of nodulation across the Viridiplantae at a systematic level. We used weighted coexpression network analysis (WGCNA) to mine a nodule‐related module (NRM) in Glycine max. Comparative genomic analysis of 78 green plant species revealed that NRM genes are recruited from different evolutionary nodes along with gene duplication events. A set of core coexpressed genes within legumes may play vital roles in regulating nodule environments essential for nitrogen fixation, including oxygen concentrations, sulfur transport, and iron homeostasis (such as GmCHY). The regulation of these genes occurred mainly at the transcription level, although some of them, such as sulfate transporters, may also undergo positive selection at protein level. We revealed that ancient orthologs and duplication events before the origin of legumes were preadapted for symbiosis. Conserved coregulated genes found within legumes paved the way for nodule formation and nitrogen fixation. These findings provide significant insights into the evolution of nodulation and indicate promising candidates for identifying other key components of legume nodulation and nitrogen fixation.