Genome-wide identification and expression analysis of the WNK kinase gene family in soybean.

WNK kinases are a unique class of serine/threonine protein kinases that lack a conserved catalytic lysine residue in the kinase domain, hence the name WNK (with no K, i.e., lysine). WNK kinases are involved in various physiological processes in plants, such as circadian rhythm, flowering time, and stress responses. In this study, we identified 26 WNK genes in soybean and analyzed their phylogenetic relationships, gene structures, chromosomal distribution, cis-regulatory elements, expression patterns, and conserved protein motifs. The soybean WNK genes were unevenly distributed on 15 chromosomes and underwent 21 segmental duplication events during evolution. We detected 14 types of cis-regulatory elements in the promoters of the WNK genes, indicating their potential involvement in different signaling pathways. The transcriptome database revealed tissue-specific and salt stress-responsive expression of WNK genes in soybean, the second of which was confirmed by salt treatments and qRT-PCR analysis. We found that most WNK genes were significantly up-regulated by salt stress within 3 h in both roots and leaves, except for WNK5, which showed a distinct expression pattern. Our findings provide valuable insights into the molecular characteristics and evolutionary history of the soybean WNK gene family and lay a foundation for further analysis of WNK gene functions in soybean. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01440-5.

Genome-Wide Identification and Characterization of the Soybean Snf2 Gene Family and Expression Response to Rhizobia.

Sucrose nonfermenting 2 (Snf2) family proteins are the core component of chromatin remodeling complexes that can alter chromatin structure and nucleosome position by utilizing the energy of ATP, playing a vital role in transcription regulation, DNA replication, and DNA damage repair. Snf2 family proteins have been characterized in various species including plants, and they have been found to regulate development and stress responses in Arabidopsis. Soybean (Glycine max) is an important food and economic crop worldwide, unlike other non-leguminous crops, soybeans can form a symbiotic relationship with rhizobia for biological nitrogen fixation. However, little is known about Snf2 family proteins in soybean. In this study, we identified 66 Snf2 family genes in soybean that could be classified into six groups like Arabidopsis, unevenly distributed on 20 soybean chromosomes. Phylogenetic analysis with Arabidopsis revealed that these 66 Snf2 family genes could be divided into 18 subfamilies. Collinear analysis showed that segmental duplication was the main mechanism for expansion of Snf2 genes rather than tandem repeats. Further evolutionary analysis indicated that the duplicated gene pairs had undergone purifying selection. All Snf2 proteins contained seven domains, and each Snf2 protein had at least one SNF2_N domain and one Helicase_C domain. Promoter analysis revealed that most Snf2 genes had cis-elements associated with jasmonic acid, abscisic acid, and nodule specificity in their promoter regions. Microarray data and real-time quantitative PCR (qPCR) analysis revealed that the expression profiles of most Snf2 family genes were detected in both root and nodule tissues, and some of them were found to be significantly downregulated after rhizobial infection. In this study, we conducted a comprehensive analysis of the soybean Snf2 family genes and demonstrated their responsiveness to Rhizobia infection. This provides insight into the potential roles of Snf2 family genes in soybean symbiotic nodulation.

GmIAA27 Encodes an AUX/IAA Protein Involved in Dwarfing and Multi-Branching in Soybean.

Soybean plant height and branching affect plant architecture and yield potential in soybean. In this study, the mutant dmbn was obtained by treating the cultivar Zhongpin 661 with ethylmethane sulfonate. The dmbn mutant plants were shorter and more branched than the wild type. The genetic analysis showed that the mutant trait was controlled by a semi-dominant gene. The candidate gene was fine-mapped to a 91 kb interval on Chromosome 9 by combining BSA-seq and linkage analysis. Sequence analysis revealed that Glyma.09g193000 encoding an Aux/IAA protein (GmIAA27) was mutated from C to T in the second exon of the coding region, resulting to amino acid substitution of proline to leucine. Overexpression of the mutant type of this gene in Arabidopsis thaliana inhibited apical dominance and promoted lateral branch development. Expression analysis of GmIAA27 and auxin response genes revealed that some GH3 genes were induced. GmIAA27 relies on auxin to interact with TIR1, whereas Gmiaa27 cannot interact with TIR1 owing to the mutation in the degron motif. Identification of this unique gene that controls soybean plant height and branch development provides a basis for investigating the mechanisms regulating soybean plant architecture development.

The small GTPase Rab5a and its guanine nucleotide exchange factors are involved in post-Golgi trafficking of storage proteins in developing soybean cotyledon.

Storage protein is the most abundant nutritional component in soybean seed. Morphology-based evidence has verified that storage proteins are initially synthesized on the endoplasmic reticulum, and then follow the Golgi-mediated pathway to the protein storage vacuole. However, the molecular mechanisms of storage protein trafficking in soybean remain unknown. Here, we clone the soybean homologs of Rab5 and its guanine nucleotide exchange factor (GEF) VPS9. GEF activity combined with yeast two-hybrid assays demonstrated that GmVPS9a2 might specifically act as the GEF of the canonical Rab5, while GmVPS9b functions as a common activator for all Rab5s. Subcellular localization experiments showed that GmRab5a was dually localized to the trans-Golgi network and pre-vacuolar compartments in developing soybean cotyledon cells. Expression of a dominant negative variant of Rab5a, or RNAi of either Rab5a or GmVPS9s, significantly disrupted trafficking of mRFP-CT10, a cargo marker for storage protein sorting, to protein storage vacuoles in maturing soybean cotyledons. Together, our results systematically revealed the important role of GmRab5a and its GEFs in storage protein trafficking, and verified the transient expression system as an efficient approach for elucidating storage protein trafficking mechanisms in seed.

Gene-based SNP identification and validation in soybean using next-generation transcriptome sequencing.

Gene-based molecular markers are increasingly used in crop breeding programs for marker-assisted selection. However, identification of genetic variants associated with important agronomic traits has remained a difficult task in soybean. RNA-Seq provides an efficient way, other than assessing global expression variations of coding genes, to discover gene-based SNPs at the whole genome level. In this study, RNA isolated from four soybean accessions each with three replications was subjected to high-throughput sequencing and a range of 44.2-65.9 million paired-end reads were generated for each library. A total of 75,209 SNPs were identified among different genotypes after combination of replications, 89.1% of which were located in expressed regions and 27.0% resulted in amino acid changes. GO enrichment analysis revealed that most significant enriched genes with nonsynonymous SNPs were involved in ribonucleotide binding or catalytic activity. Of 22 SNPs subjected to PCR amplification and Sanger sequencing, all of them were validated. To test the utility of identified SNPs, these validated SNPs were also assessed by genotyping a relative large population with 393 wild and cultivated soybean accessions. These SNPs identified by RNA-Seq provide a useful resource for genetic and genomic studies of soybean. Moreover, the collection of nonsynonymous SNPs annotated with their predicted functional effects also provides a valuable asset for further discovery of genes, identification of gene variants, and development of functional markers.