GATA transcription factors are transcriptional regulatory proteins that contain a characteristic type-IV zinc finger DNA-binding domain and recognize the conserved GATA motif in the promoter sequence of target genes. Previous studies demonstrated that plant GATA factors possess critical functions in developmental control and responses to the environment. To date, the GATA factors in soybean (Glycine max) have yet to be characterized. Thus, this study identified 64 putative GATA factors from the entire soybean genomic sequence. The chromosomal distributions, gene structures, duplication patterns, phylogenetic tree, tissue expression patterns, and response to low nitrogen stress of the 64 GATA factors in soybean were analyzed to further investigate the functions of these factors. Results indicated that segmental duplication predominantly contributed to the expansion of the GATA factor gene family in soybean. These GATA proteins were phylogenetically clustered into four distinct subfamilies, wherein their gene structure and motif compositions were considerably conserved. A comparative phylogenetic analysis of the GATA factor zinc finger domain sequences in soybean, Arabidopsis (Arabidopsis thaliana), and rice (Oryza sativa) revealed four major classes. The GATA factors in soybean exhibited expression diversity among different tissues; some of these factors showed tissue-specific expression patterns. Numerous GATA factors displayed upregulation or downregulation in soybean leaf in response to low nitrogen stress, and two GATA factors GATA44 and GATA58 were likely to be involved in the regulation of nitrogen metabolism in soybean. Overexpression of GmGATA44 complemented the reduced chlorophyll phenotype of the Arabidopsis ortholog AtGATA21 mutant, implying that GmGATA44 played an important role in modulating chlorophyll biosynthesis. Overall, our study provides useful information for the further analysis of the biological functions of GATA factors in soybean and other crops.

Nod Factor Receptor5 (NFR5) is an atypical receptor-like kinase, having no activation loop in the protein kinase domain. It forms a heterodimer with NFR1 and is required for the early plant responses to Rhizobium infection. A Rho-like small GTPase from Lotus japonicus was identified as an NFR5-interacting protein. The amino acid sequence of this Rho-like GTPase is closest to the Arabidopsis (Arabidopsis thaliana) ROP6 and Medicago truncatula ROP6 and was designated as LjROP6. The interaction between Rop6 and NFR5 occurred both in vitro and in planta. No interaction between Rop6 and NFR1 was observed. Green fluorescent protein-tagged ROP6 was localized at the plasma membrane and cytoplasm. The interaction between ROP6 and NFR5 appeared to take place at the plasma membrane. The expression of the ROP6 gene could be detected in vascular tissues of Lotus roots. After inoculation with Mesorhizobium loti, elevated levels of ROP6 expression were found in the root hairs, root tips, vascular bundles of roots, nodule primordia, and young nodules. In transgenic hairy roots expressing ROP6 RNA interference constructs, Rhizobium entry into the root hairs did not appear to be affected, but infection thread growth through the root cortex were severely inhibited, resulting in the development of fewer nodules per plant. These data demonstrate a role of ROP6 as a positive regulator of infection thread formation and nodulation in L. japonicus.

Background The plant architecture has significant effects on grain yield of various crops, including soybean ( Glycine max ), but the knowledge on optimization of plant architecture in order to increase yield potential is still limited. Recently, CRISPR/Cas9 system has revolutionized genome editing, and has been widely utilized to edit the genomes of a diverse range of crop plants. Results In the present study, we employed the CRISPR/Cas9 system to mutate four genes encoding SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors of the SPL9 family in soybean. These four GmSPL9 genes are negatively regulated by GmmiR156b , a target for the improvement of soybean plant architecture and yields. The soybean Williams 82 was transformed with the binary CRISPR/Cas9 plasmid, assembled with four sgRNA expression cassettes driven by the Arabidopsis thaliana U3 or U6 promoter, targeting different sites of these four SPL9 genes via Agrobacterium tumefaciens -mediated transformation. A 1-bp deletion was detected in one target site of the GmSPL9a and one target site of the GmSPL9b , respectively, by DNA sequencing analysis of two T0-generation plants. T2-generation spl9a and spl9b homozygous single mutants exhibited no obvious phenotype changes; but the T2 double homozygous mutant spl9a / spl9b possessed shorter plastochron length. In T4 generation, higher-order mutant plants carrying various combinations of mutations showed increased node number on the main stem and branch number, consequently increased total node number per plants at different levels. In addition, the expression levels of the examined GmSPL9 genes were higher in the spl9b-1 single mutant than wild-type plants, which might suggest a feedback regulation on the expression of the investigated GmSPL9 genes in soybean. Conclusions Our results showed that CRISPR/Cas9-mediated targeted mutagenesis of four GmSPL9 genes in different combinations altered plant architecture in soybean. The findings demonstrated that GmSPL9a, GmSPL9b, GmSPL9c and GmSPL9 function as redundant transcription factors in regulating plant architecture in soybean. Electronic supplementary material The online version of this article (10.1186/s12870-019-1746-6) contains supplementary material, which is available to authorized users.

Root nodule symbiosis (RNS) is one of the most productive and economical systems for nitrogen fixation, and previous studies have shown that several nodule-specific C2H2-zinc finger proteins (ZFPs) play important roles in symbiosis establishment and nodule function. However, C2H2-ZFPs are the most widespread ZFPs in eukaryotes, and a great variation of structure and function exist among the family members. It remains largely unclear whether or not special types of C2H2-ZF genes participate in symbiosis, especially in soybean. In the present study, we performed a genome-wide survey of soybean C2H2-ZF genes, and 321 soybean C2H2-ZF genes were identified and classified into 11 clearly distinguishable subsets (Gm-t1-SF, Gm-t2-SF, Gm-1i-Q-SF, Gm-1i-M-SF, Gm-1i-Z-SF, Gm-1i-D-SF, Gm-2i-Q-SF, Gm-2i-M-SF, Gm-2i-Mix-SF, Gm-3i-SF, and Gm-4i-SF) based on the arrangements, numbers, and types of C2H2-ZF domains. Phylogenetic and gene ontology analyses were carried out to assess the conserved sequence and GO function among these subsets, and the results showed that the classification of soybean C2H2-ZFPs was reasonable. The expression profile of soybean C2H2-ZFPs in multiple tissues showed that nearly half of soybean C2H2-ZFPs within different subsets had expressions in nodules, including a clustering branch consisting of 11 Gm-1i-Q-SF genes specifically expressed in symbiotic-relative tissues. RNA-Seq was used to identify symbiosis-related soybean C2H2-ZFPs, and the expression pattern of the soybean C2H2-ZFPs in roots and nodules at different development stages showed that soybean C2H2-ZFPs mainly played roles in nodule development or nodule function rather than nodulation signal transduction, and nearly half of these genes had high expressions and/or different expression patterns during soybean nodule development, especially for the six clustering branches of genes consisting of different subsets of C2H2-ZFPs. Furthermore, the selected symbiosis-related soybean C2H2-ZFPs might function in legume-rhizobium symbiosis through regulating or interacting with other key proteins. Taken together, our findings provided useful information for the study on classification and conservative function of C2H2-ZFPs, and offered solid evidence for investigation of rhizobium symbiosis-related C2H2-ZFPs in soybean or other legumes.

The symbiosis receptor kinase (SymRK) is required for morphological changes of legume root hairs triggered by rhizobial infection. How protein turnover of SymRK is regulated and how the nodulation factor signals are transduced downstream of SymRK are not known. In this report, a SymRK-interacting E3 ubiquitin ligase (SIE3) was shown to bind and ubiquitinate SymRK. The SIE3-SymRK interaction and the ubiquitination of SymRK were shown to occur in vitro and in planta. SIE3 represents a new class of plant-specific E3 ligases that contain a unique pattern of the conserved CTLH (for C-terminal to LisH), CRA (for CT11-RanBPM), and RING (for Really Interesting New Gene) domains. Expression of SIE3 was detected in all tested tissues of Lotus japonicus plants, and its transcript level in roots was enhanced by rhizobial infection. The SIE3 protein was localized to multiple subcellular locations including the nuclei and plasma membrane, where the SIE3-SymRK interaction took place. Overexpression of SIE3 promoted nodulation in transgenic hairy roots, whereas downregulation of SIE3 transcripts by RNA interference inhibited infection thread development and nodule organogenesis. These results suggest that SIE3 represents a new class of E3 ubiquitin ligase, acts as a regulator of SymRK, and is involved in rhizobial infection and nodulation in L. japonicus.

The root nodule symbiosis (RNS) between legume plants and rhizobia is the most efficient and productive source of nitrogen fixation, and has critical importance in agriculture and mesology. Soybean (Glycine max), one of the most important legume crops in the world, establishes a nitrogen-fixing symbiosis with different types of rhizobia, and the efficiency of symbiotic nitrogen fixation in soybean greatly depends on the symbiotic host-specificity. Although, it has been reported that rhizobia use surface polysaccharides, secretion proteins of the type-three secretion systems and nod factors to modulate host range, the host control of nodulation specificity remains poorly understood. In this report, the soybean roots of two symbiotic systems (Bradyrhizobium japonicum strain 113-2-soybean and Sinorhizobium fredii USDA205-soybean)with notable different nodulation phenotypes and the control were studied at five different post-inoculation time points (0.5, 7–24 h, 5, 16, and 21 day) by RNA-seq (Quantification). The results of qPCR analysis of 11 randomly-selected genes agreed with transcriptional profile data for 136 out of 165 (82.42%) data points and quality assessment showed that the sequencing library is of quality and reliable. Three comparisons (control vs. 113-2, control vs. USDA205 and USDA205 vs. 113-2) were made and the differentially expressed genes (DEGs) between them were analyzed. The number of DEGs at 16 days post-inoculation (dpi) was the highest in the three comparisons, and most of the DEGs in USDA205 vs. 113-2 were found at 16 dpi and 21 dpi. 44 go function terms in USDA205 vs. 113-2 were analyzed to evaluate the potential functions of the DEGs, and 10 important KEGG pathway enrichment terms were analyzed in the three comparisons. Some important genes induced in response to different strains (113-2 and USDA205) were identified and analyzed, and these genes primarily encoded soybean resistance proteins, NF-related proteins, nodulins and immunity defense proteins, as well as proteins involving flavonoids/flavone/flavonol biosynthesis and plant-pathogen interaction. Besides, 189 candidate genes are largely expressed in roots and\or nodules. The DEGs uncovered in this study provides molecular candidates for better understanding the mechanisms of symbiotic host-specificity and explaining the different symbiotic effects between soybean roots inoculated with different strains (113-2 and USDA205).

In the Rhizobium-legume symbiosis, calcium/calmodulin-dependent protein kinase (CCaMK) is a key regulator for both rhizobial infection and nodule organogenesis. Deregulation of CCaMK by either a point mutation in the autophosphorylation site or the deletion of the carboxyl-terminal regulatory domain results in spontaneous nodule formation without rhizobia. However, the underlying biochemical mechanisms are poorly understood. Here, using the kinase domain of CCaMK as a bait in yeast two-hybrid screening, we identify a novel protein, CIP73 (for CCaMK-interacting protein of approximately 73 kD), that interacts with CCaMK. CIP73 contains a Scythe_N ubiquitin-like domain and belongs to the large ubiquitin superfamily.Deletion and mutagenesis analysis demonstrate that CIP73 could only interact with CCaMK when the calmodulin-binding domain and three EF-hand motifs are removed from the kinase domain. The amino-terminal 80 amino acid residues (80-160) of CCaMK are required for interacting with CIP73 in yeast cells. On the other hand, protein pull-down assay and bimolecular fluorescence complementation assay in Nicotiana benthamiana show that the full-length CCaMK could interact with CIP73 in vitro and in planta. Importantly, CCaMK phosphorylates the amino terminus of CIP73 in a Ca 2+ /calmodulin-dependent manner in vitro. CIP73 transcripts are preferentially expressed in roots, and very low expression is detected in leaves, stems, and nodules. The expression in roots is significantly decreased after inoculation of Mesorhizobium loti. RNA interference knockdown of CIP73 expression by hairy root transformation in Lotus japonicus led to decreased nodule formation, suggesting that CIP73 performed an essential role in nodulation.

Summary MYB transcription factors (TFs) have been reported to regulate the biosynthesis of secondary metabolites, as well as to mediate plant adaption to abiotic stresses, including drought. However, the roles of MYB TFs in regulating plant architecture and yield potential remain poorly understood. Here, we studied the roles of the dehydration‐inducible GmMYB14 gene in regulating plant architecture, high‐density yield and drought tolerance through the brassinosteroid (BR) pathway in soybean. GmMYB14 was shown to localize to nucleus and has a transactivation activity. Stable GmMYB14‐overexpressing (GmMYB14‐OX) transgenic soybean plants displayed a semi‐dwarfism and compact plant architecture associated with decreased cell size, resulting in a decrease in plant height, internode length, leaf area, leaf petiole length and leaf petiole angle, and improved yield in high density under field conditions. Results of the transcriptome sequencing suggested the involvement of BRs in regulating GmMYB14‐OX plant architecture. Indeed, GmMYB14‐OX plants showed reduced endogenous BR contents, while exogenous application of brassinolide could partly rescue the phenotype of GmMYB14‐OX plants. Furthermore, GmMYB14 was shown to directly bind to the promoter of GmBEN1 and up‐regulate its expression, leading to reduced BR content in GmMYB14‐OX plants. GmMYB14‐OX plants also displayed improved drought tolerance under field conditions. GmBEN1 expression was also up‐regulated in the leaves of GmMYB14‐OX plants under polyethylene glycol treatment, indicating that the GmBEN1‐mediated reduction in BR level under stress also contributed to drought/osmotic stress tolerance of the transgenic plants. Our findings provided a strategy for stably increasing high‐density yield and drought tolerance in soybean using a single TF‐encoding gene.