The Bradyrhizobium japonicum hsfA gene exhibits a unique developmental expression pattern in cowpea nodules

Integrated single-nucleus and spatial transcriptomics captures transitional states in soybean nodule maturation

Legumes form symbiosis with rhizobium leading to the development of nitrogen-fixing nodules. By integrating single-nucleus and spatial transcriptomics, we established a cell atlas of soybean nodules and roots. In central infected zones of nodules, we found that uninfected cells specialize into functionally distinct subgroups during nodule development, and revealed a transitional subtype of infected cells with enriched nodulation-related genes. Overall, our results provide a single-cell perspective for understanding rhizobium-legume symbiosis.

GmSPX8, a nodule-localized regulator confers nodule development and nitrogen fixation under phosphorus starvation in soybean

BACKGROUND

Biological nitrogen fixation (BNF) is an important nitrogen source for legume plants, and highly efficient nitrogen fixation requires sufficient phosphorus (P). However, the mechanism of maintaining nitrogen fixation of the legume nodules under low P concentration remains largely unknown.

RESULTS

A nodule-localized SPX protein, GmSPX8, was discovered by transcriptome and functional analysis of its role in N2 fixation was characterized in soybean nodules. GmSPX8 was preferentially expressed in nodules and its expression was gradually increased during nodule development. And also the expression pattern was investigated using reporter gene β-glucuronidase (GUS) driven by the promoter of GmSPX8. GmSPX8 was greatly induced and the GUS activity was increased by 12.2% under P deficiency. Overexpression of GmSPX8 in transgenic plants resulted in increased nodule number, nodule fresh weight and nitrogenase activity by 15.0%, 16.0%, 42.5%, subsequently leading to increased N and P content by 17.0% and 19.0%, while suppression of GmSPX8 showed significantly impaired nodule development and nitrogen fixation efficiency under low P stress. These data indicated that GmSPX8 conferred nodule development and nitrogen fixation under low P condition. By yeast two-hybrid screening, GmPTF1 was identified as a potential interacting protein of GmSPX8, which was further confirmed by BiFC, Y2H and pull down assay. Transcript accumulation of GmPTF1 and its downstream genes such as GmEXLB1 and EXPB2 were increased in GmSPX8 overexpressed transgenic nodules, and in the presence of GmSPX8, the transcriptional activity of GmPTF1 in yeast cells and tobacco leaves was greatly enhanced.

CONCLUSIONS

In summary, these findings contribute novel insights towards the role of GmSPX8 in nodule development and nitrogen fixation partly through interacting with GmPTF1 in soybean under low P condition.

A Nodule-Localized Small Heat Shock Protein GmHSP17.1 Confers Nodule Development and Nitrogen Fixation in Soybean

Small heat shock proteins (sHSPs) are ubiquitous proteins present in all organisms. The sHSPs are not only upregulated under heat shock as well as other stresses but also are expressed in unstressed cells, indicating quite diverse functions of sHSPs. However, there is little known about the role of sHSPs in nodulation and nitrogen fixation in soybean. In this study, we cloned a candidate protein of sHSP, GmHSP17.1, from proteome of nodule and analyzed its function in soybean nodulation. We found that GmHSP17.1 was a cytosolic protein and preferentially expressed during nodule development. An overexpression of GmHSP17.1 in composite transgenic plants showed increases in nodule number, fresh weight, nodule size, area of infection cells, and nitrogenase activity, and subsequently promoted the content of nitrogen and growth of soybean plants. While GmHSP17.1 RNA interference (RNAi) lines showed significantly impaired nodule development and nitrogen fixation efficiency. Through liquid chromatography-tandem mass spectrometry (LC-MS/MS), GmRIP1 was identified as the first potential target of GmHSP17.1, and was shown to be specifically expressed in soybean nodules. The interaction between GmHSP17.1 and GmRIP1 was further confirmed by yeast-two hybrid (Y2H), bimolecular fluorescence complementation (BiFC) in vivo and pull-down assay in vitro. Furthermore, peroxidase activity was markedly increased in GmHSP17.1 overexpressed nodules and decreased in RNAi lines. As a result, the reactive oxygen species (ROS) content greatly decreased in GmHSP17.1 overexpression lines and increased in suppression lines. Taken together, we conclude that GmHSP17.1 plays an important role in soybean nodulation through interacting with GmRIP1. Our results provide foundation for studying the mechanism of nitrogen fixation and for the genetics improvement of legume plants.

GmPAP12 Is Required for Nodule Development and Nitrogen Fixation Under Phosphorus Starvation in Soybean

Nodulation process in legume plants is essential for biological nitrogen fixation during which process a large amount of phosphorus (P) is required. Under P deficiency, nodule formation is greatly affected, and induction of purple acid phosphatases (PAPs) is an adaptive strategy for nodules to acquire more P. However, regulation roles of PAPs in nodules remain largely understood. In this study, by transcriptome sequencing technology, five PAP genes were found to be differentially expressed, which led to the greatly increased acid phosphatase (APase) and phytase activities in soybean mature nodules under P starvation conditions; and among the five PAP genes, GmPAP12 had the highest transcript level, and RT-PCR indicated expression of GmPAP12 was gradually increasing during nodule development. GUS activity driven by GmPAP12 promoter was also significantly induced in low phosphorus conditions. Further functional analysis showed that under low phosphorus stress, overexpression of GmPAP12 resulted in higher nodule number, fresh weight, and nitrogenase activity as well as the APase activity than those of control plant nodules, whereas the growth performance and APase activity of nodules on hairy roots were greatly lower when GmPAP12 was suppressed, indicating that GmPAP12 may promote P utilization in soybean nodules under low P stress, which thus played an important role in nodulation and biological nitrogen fixation. Moreover, P1BS elements were found in the promoter of GmPAP12, and yeast one-hybrid experiment further proved the binding of P1BS by transcription factor GmPHR1 in the promoter of GmPAP12. At last, overexpression and suppression of GmPHR1 in nodules indeed caused highly increased and decreased expression of GmPAP12, respectively, indicating that GmPAP12 is regulated by GmPHR1 in soybean nodules. Taken together, these data suggested that GmPAP12 was a novel soybean PAP involved in the P utilization and metabolism in soybean root nodules and played an important role in the growth and development of root nodules and biological nitrogen fixation.

Genetic divergence of bradyrhizobium strains nodulating soybeans as revealed by multilocus sequence analysis of genes inside and outside the symbiosis island

The genus Bradyrhizobium has been considered to be a taxonomically difficult group. In this study, phylogenetics and evolutionary genetics analyses were used to investigate divergence levels among Bradyrhizobium strains nodulating soybeans in China. Eleven genospecies were identified by sequence analysis of three phylogenetic and taxonomic markers (SMc00019, thrA, and truA). This was also supported by analyses of eight genes outside the symbiosis island ("off-island" genes; SMc00019, thrA, truA, fabB, glyA, phyR, exoN, and hsfA). However, seven genes inside the symbiosis island ("island" genes; nifA, nifH, nodC, nodV, fixA, trpD, and rhcC2) showed contrasting lower levels of nucleotide diversity and recombination rates than did off-island genes. Island genes had significantly incongruent gene phylogenies compared to the species tree. Four phylogenetic clusters were observed in island genes, and the epidemic cluster IV (harbored by Bradyrhizobium japonicum, Bradyrhizobium diazoefficiens, Bradyrhizobium huanghuaihaiense, Bradyrhizobium liaoningense, Bradyrhizobium daqingense, Bradyrhizobium sp. I, Bradyrhizobium sp. III, and Bradyrhizobium sp. IV) was not found in Bradyrhizobium yuanmingense, Bradyrhizobium sp. II, or Bradyrhizobium elkanii. The gene flow level of island genes among genospecies is discussed in the context of the divergence level of off-island genes.

Functional implication of β-carotene hydroxylases in soybean nodulation

Legume-Rhizobium spp. symbiosis requires signaling between the symbiotic partners and differential expression of plant genes during nodule development. Previously, we cloned a gene encoding a putative β-carotene hydroxylase (GmBCH1) from soybean (Glycine max) whose expression increased during nodulation with Bradyrhizobium japonicum. In this work, we extended our study to three GmBCHs to examine their possible role(s) in nodule development, as they were additionally identified as nodule specific, along with the completion of the soybean genome. In situ hybridization revealed the expression of three GmBCHs (GmBCH1, GmBCH2, and GmBCH3) in the infected cells of root nodules, and their enzymatic activities were confirmed by functional assays in Escherichia coli. Localization of GmBCHs by transfecting Arabidopsis (Arabidopsis thaliana) protoplasts with green fluorescent protein fusions and by electron microscopic immunogold detection in soybean nodules indicated that GmBCH2 and GmBCH3 were present in plastids, while GmBCH1 appeared to be cytosolic. RNA interference of the GmBCHs severely impaired nitrogen fixation as well as nodule development. Surprisingly, we failed to detect zeaxanthin, a product of GmBCH, or any other carotenoids in nodules. Therefore, we examined the possibility that most of the carotenoids in nodules are converted or cleaved to other compounds. We detected the expression of some carotenoid cleavage dioxygenases (GmCCDs) in wild-type nodules and also a reduced amount of zeaxanthin in GmCCD8-expressing E. coli, suggesting cleavage of the carotenoid. In view of these findings, we propose that carotenoids such as zeaxanthin synthesized in root nodules are cleaved by GmCCDs, and we discuss the possible roles of the carotenoid cleavage products in nodulation.

RNA interference-mediated repression of S6 kinase 1 impairs root nodule development in soybean

Symbiotic nodule formation on legume roots is characterized with a series of developmental reprograming in root tissues, including extensive proliferation of cortical cells. We examined a possible involvement of the target of rapamycin (TOR) pathway, a central regulator of cell growth and proliferation in animals and yeasts, during soybean nodule development. Our results show that transcription of both GmTOR and its key downstream effector, GmS6K1, are activated during nodulation, which is paralleled with higher kinase activities of these gene products as well. RNAi-mediated knockdown of GmS6K1 impaired the nodule development with severely reduced nodule weight and numbers. In addition, expression of a few nodulins including leghemoglobin was also decreased, and consequently nitrogen fixation was found to be reduced by half. Proteomic analysis of the GmS6K1-RNAi nodules identified glutamine synthetase (GS), an essential enzyme for nitrogen assimilation in nodules, as one of the proteins that are significantly down regulated. These results appear to provide solid evidence for a functional link between GmS6K1 and nodule development.

An atypical soybean leucine-rich repeat receptor-like kinase, GmLRK1, may be involved in the regulation of cell elongation

A leucine-rich repeat receptor-like kinase (LRR-RLK) encoded by one of the genes highly expressed in a specific stage of soybean seed development, referred to as GmLRK1, was identified and characterized. Examination of its kinase domain indicated that GmLRK1 may be a catalytically inactive atypical receptor kinase. An autophosphorylation assay confirmed that GmLRK1 is incapable of autophosphorylation in vitro. However, the phosphorylation of GmRLK1 could be induced after incubation with plant protein extracts, suggesting that some plant proteins may interact with GmLRK1 and phosphorylate the protein in vivo. Analyses of the expression profiles of GmLRK1 and its Arabidopsis ortholog At2g36570 revealed that they may be involved in regulation of more fundamental metabolic and/or developmental pathways, rather than a specialized developmental program such as seed development. Our results further indicate that the GmLRK1 and At2g36570 may play a role in the regulation of certain cellular processes that lead to cell elongation and expansion.

Overexpression of GmAKR1, a stress-induced aldo/keto reductase from soybean, retards nodule development

Development of symbiotic root nodules in legumes involves the induction and repression of numerous genes in conjunction with changes in the level of phytohormones. We have isolated several genes that exhibit differential expression patterns during the development of soybean nodules. One of such genes, which were repressed in mature nodules, was identified as a putative aldo/keto reductase and thus named Glycine max aldo/keto reductase 1 (GmAKR1). GmAKR1 appears to be a close relative of a yeast aldo/keto reductase YakC whose in vivo substrate has not been identified yet. The expression of GmAKR1 in soybean showed a root-specific expression pattern and inducibility by a synthetic auxin analogue 2,4-D, which appeared to be corroborated by presence of the root-specific element and the stress-response element in the promoter region. In addition, constitutive overexpression of GmAKR1 in transgenic soybean hairy roots inhibited nodule development, which suggests that it plays a negative role in the regulation of nodule development. One of the Arabidopsis orthologues of GmAKR1 is the ARF-GAP domain 2 protein, which is a potential negative regulator of vesicle trafficking; therefore GmAKR1 may have a similar function in the roots and nodules of legume plants.

Induction of thioredoxin is required for nodule development to reduce reactive oxygen species levels in soybean roots

Nodules are formed on legume roots as a result of signaling between symbiotic partners and in response to the activities of numerous genes. We cloned fragments of differentially expressed genes in spot-inoculated soybean (Glycine max) roots. Many of the induced clones were similar to known genes related to oxidative stress, such as thioredoxin and beta-carotene hydroxylase. The deduced amino acid sequences of full-length soybean cDNAs for thioredoxin and beta-carotene hydroxylase were similar to those in other species. In situ RNA hybridization revealed that the thioredoxin gene is expressed on the pericycle of 2-d-old nodules and in the infected cells of mature nodules, suggesting that thioredoxin is involved in nodule development. The thioredoxin promoter was found to contain a sequence resembling an antioxidant responsive element. When a thioredoxin mutant of yeast was transformed with the soybean thioredoxin gene it became hydrogen peroxide tolerant. These observations prompted us to measure reactive oxygen species levels. These were decreased by 3- to 5-fold in 7-d-old and 27-d-old nodules, coincident with increases in the expression of thioredoxin and beta-carotene hydroxylase genes. Hydrogen peroxide-producing regions identified with cerium chloride were found in uninoculated roots and 2-d-old nodules, but not in 7-d-old and 27-d-old nodules. RNA interference-mediated repression of the thioredoxin gene severely impaired nodule development. These data indicate that antioxidants such as thioredoxin are essential to lower reactive oxygen species levels during nodule development.