Agrobacterium rhizogenes-mediated transformation of soybean to study root biology

Functional Identification of miR2119 Targeting ADHs in Modulating Soybean Resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is a sedentary endoparasite nematode that results in severe economic losses in soybean crops. miRNAs play crucial roles in plant responses to nematode. However, the role of miR2119 responding to SCN stress in soybean. Here, we demonstrated that the transcript levels of polycistronic precursors containing miR2119 and miR398a were significantly reduced in soybean upon nematode infection. Promoter of the miR2119-398a precursor analysis was conducted containing a GUS reporter gene. GUS activity assays demonstrated a decrease in miR2119-398a promoter during SCN infection. Overexpression of polycistronic precursor miR2119-398a (OE-premiR2119-398a) and miR2119 precursor (OE-premiR2119) rendered soybean more susceptible to SCN. Conversely, silencing miR2119 (STTM2119) increased soybean resistance against SCN. Furthermore, RNA-seq analysis revealed that miR2119 is involved in many defense signaling pathways. GUS reporter gene assays demonstrated that miR2119 targets GmADH1.1a and GmADH1.1b. Functional analysis indicated that ADHs act as a major role in responding to H. glycines by modulating reactive oxygen species (ROS) levels. Together, the findings reveal a novel mechanism by which the polycistronic precursor miR2119-398a coordinately regulates in response to H. glycines. Additionally, miR2119 becomes an essential element contributing to H. glycines by modulating ADH activity and ROS homeostasis in soybean.

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Loss of Lateral suppressor gene is associated with evolution of root nodule symbiosis in Leguminosae

BACKGROUND

Root nodule symbiosis (RNS) is a fascinating evolutionary event. Given that limited genes conferring the evolution of RNS in Leguminosae have been functionally validated, the genetic basis of the evolution of RNS remains largely unknown. Identifying the genes involved in the evolution of RNS will help to reveal the mystery.

RESULTS

Here, we investigate the gene loss event during the evolution of RNS in Leguminosae through phylogenomic and synteny analyses in 48 species including 16 Leguminosae species. We reveal that loss of the Lateral suppressor gene, a member of the GRAS-domain protein family, is associated with the evolution of RNS in Leguminosae. Ectopic expression of the Lateral suppressor (Ls) gene from tomato and its homolog MONOCULM 1 (MOC1) and Os7 from rice in soybean and Medicago truncatula result in almost completely lost nodulation capability. Further investigation shows that Lateral suppressor protein, Ls, MOC1, and Os7 might function through an interaction with NODULATION SIGNALING PATHWAY 2 (NSP2) and CYCLOPS to repress the transcription of NODULE INCEPTION (NIN) to inhibit the nodulation in Leguminosae. Additionally, we find that the cathepsin H (CTSH), a conserved protein, could interact with Lateral suppressor protein, Ls, MOC1, and Os7 and affect the nodulation.

CONCLUSIONS

This study sheds light on uncovering the genetic basis of the evolution of RNS in Leguminosae and suggests that gene loss plays an essential role.

Phenoxyacetic acid enhances nodulation symbiosis during the rapid growth stage of soybean

Root exudates are known signaling agents that influence legume root nodulation, but the molecular mechanisms for nonflavonoid molecules remain largely unexplored. The number of soybean root nodules during the initial growth phase shows substantial discrepancies at distinct developmental junctures. Using a combination of metabolomics analyses on root exudates and nodulation experiments, we identify a pivotal role for certain root exudates during the rapid growth phase in promoting nodulation. Phenoxyacetic acid (POA) was found to activate the expression of GmGA2ox10 and thereby facilitate rhizobial infection and the formation of infection threads. Furthermore, POA exerts regulatory control on the miR172c-NNC1 module to foster nodule primordia development and consequently increase nodule numbers. These findings collectively highlight the important role of POA in enhancing nodulation during the accelerated growth phase of soybeans.

Genome-Wide Identification of the Soybean AlkB Homologue Gene Family and Functional Characterization of GmALKBH10Bs as RNA m6A Demethylases and Expression Patterns under Abiotic Stress

Soybean (Glycine max (L.) Merr) is one of the most important crops worldwide, but its yield is vulnerable to abiotic stresses. In Arabidopsis, the AlkB homologue (ALKBH) family genes plays a crucial role in plant development and stress response. However, the identification and functions of its homologous genes in soybean remain obscured. Here, we identified a total of 22 ALKBH genes in soybean and classified them into seven subfamilies according to phylogenetic analysis. Gene duplication events among the family members and gene structure, conserved domains, and motifs of all candidate genes were analyzed. By comparing the changes in the m6A levels on mRNA from hair roots between soybean seedlings harboring the empty vector and those harboring the GmALKBH10B protein, we demonstrated that all four GmALKBH10B proteins are bona fide m6A RNA demethylases in vivo. Subcellular localization and expression patterns of the GmALKBH10B revealed that they might be functionally redundant. Furthermore, an analysis of cis-elements coupled with gene expression data demonstrated that GmALKBH10B subfamily genes, including GmALKBH10B1, GmALKBH10B2, GmALKBH10B3, and GmALKBH10B4, are likely involved in the cis-elements' response to various environmental stimuli. In summary, our study is the first to report the genome-wide identification of GmALKBH family genes in soybean and to determine the function of GmALKBH10B proteins as m6A RNA demethylases, providing insights into GmALKBH10B genes in response to abiotic stresses.

Transcription factor GmAlfin09 regulates endoplasmic reticulum stress in soybean via peroxidase GmPRDX6

Soybean (Glycine max [L.] Merr.) is a valuable oil crop but is also highly susceptible to environmental stress. Thus, developing approaches to enhance soybean stress resistance is vital to soybean yield improvement. In previous studies, transcription factor Alfin has been shown to serve as an epigenetic regulator of plant growth and development. However, no studies on Alfin have yet been reported in soybean. In this study, the endoplasmic reticulum (ER) stress- and reactive oxygen species (ROS)-related GmAlfin09 was identified. Screening of genes co-expressed with GmAlfin09 unexpectedly led to the identification of soybean peroxidase 6 (GmPRDX6). Further analyses revealed that both GmAlfin09 and GmPRDX6 were responsive to ER stress, with GmPRDX6 localizing to the ER under stress. Promoter binding experiments confirmed the ability of GmAlfin09 to bind to the GmPRDX6 promoter directly. When GmAlfin09 and GmPRDX6 were overexpressed in soybean, enhanced ER stress resistance and decreased ROS levels were observed. Together, these findings suggest that GmAlfin09 promotes the upregulation of GmPRDX6, and GmPRDX6 subsequently localizes to the ER, reduces ROS levels, promotes ER homeostasis, and ensures the normal growth of soybean even under ER stress. This study highlights a vital target gene for future molecular breeding of stress-resistant soybean lines.

Identification and Expression Analysis of Sulfate Transporter Genes Family and Function Analysis of GmSULTR3;1a from Soybean

Sulfate transporters (SULTRs) are essential for the transport and absorption of sulfate in plants and serve as critical transport proteins within the sulfur metabolism pathway, significantly influencing plant growth, development, and stress adaptation. A bioinformatics analysis of SULTR genes in soybean was performed, resulting in the identification and classification of twenty-eight putative GmSULTRs into four distinct groups. In this study, the characteristics of the 28 GmSULTR genes, including those involved in collinearity, gene structure, protein motifs, cis-elements, tissue expression patterns, and the response to abiotic stress and plant hormone treatments, were systematically analyzed. This study focused on conducting a preliminary functional analysis of the GmSULTR3;1a gene, wherein a high expression level of GmSULTR3;1a in the roots, stems, and leaves was induced by a sulfur deficiency and GmSULTR3;1a improved the salt tolerance. A further functional characterization revealed that GmSULTR3;1a-overexpressing soybean hairy roots had higher SO42-, GSH, and methionine (Met) contents compared with the wild-type (WT) plant. These results demonstrate that the overexpression of GmSULTR3;1a may promote the sulfur assimilation metabolism and increase the content of sulfur-containing amino acids in plants.

Soybean symbiotic-nodule zonation and cell differentiation are defined by NIN2 signaling and GH3-dependent auxin homeostasis

Symbiotic nodules comprise two classes, indeterminate and determinate, defined by the presence/absence of apical meristem and developmental zonation. Why meristem and zonation are absent from determinate nodules remains unclear. Here, we define cell types in developing soybean nodules, highlighting the undifferentiated infection zones and differentiated nitrogen-fixation zones. Auxin governs infection zone maintenance. GRETCHEN HAGEN 3 (GH3) enzymes deactivate auxin by conjugation and promote cell differentiation. gh3 mutants increased undifferentiated cells and enlarged infection zones. The central symbiosis-transcription factor NIN2a activates GH3.1 to reduce auxin levels and facilitates cell differentiation. High auxin promotes NIN2a protein accumulation and enhances signaling, further deactivating auxin and depleting infection zones. Our findings shed light on the NIN2a-GH3-auxin module that drives soybean nodule cell differentiation. This study challenges our understanding of determinate nodule development and proposes that the regulation of nodule zonation offers valuable insights into broader mechanisms of cell differentiation across plant species.

MiR398b Targets Superoxide Dismutase Genes in Soybean in Defense Against Heterodera glycines via Modulating Reactive Oxygen Species Homeostasis

MicroRNAs play crucial roles in plant defense responses. However, the underlying mechanism by which miR398b contributes to soybean responses to soybean cyst nematode (Heterodera glycines) remains elusive. In this study, by using Agrobacterium rhizogenes-mediated transformation of soybean hairy roots, we observed that miR398b and target genes GmCCS and GmCSD1b played vital functions in soybean-H. glycines interaction. The study revealed that the abundance of miR398b was downregulated by H. glycines infection, and overexpression of miR398b enhanced the susceptibility of soybean to H. glycines. Conversely, silencing of miR398b improved soybean resistance to H. glycines. Detection assays revealed that miR398b rapidly senses stress-induced reactive oxygen species, leading to the repression of target genes GmCCS and GmCSD1b and regulating the accumulation of plant defense genes against nematode infection. Moreover, exogenous synthetic ds-miR398b enhanced soybean sensitivity to H. glycines by modulating H2O2 and O2- levels. Functional analysis demonstrated that overexpression of GmCCS and GmCSD1b in soybean enhanced resistance to H. glycines. RNA interference-mediated repression of GmCCS and GmCSD1b in soybean increased susceptibility to H. glycines. RNA sequencing revealed that a majority of differentially expressed genes in overexpressed GmCCS were associated with oxidative stress. Overall, the results indicate that miR398b targets superoxide dismutase genes, which negatively regulate soybean resistance to H. glycines via modulating reactive oxygen species levels and defense signals.

The role of GmHSP23.9 in regulating soybean nodulation under elevated CO2 condition

Legume-rhizobia symbiosis offers a unique approach to increase leguminous crop yields. Previous studies have indicated that the number of soybean nodules are increased under elevated CO2 concentration. However, the underlying mechanism behind this phenomenon remains elusive. In this study, transcriptome analysis was applied to identify candidate genes involved in regulating soybean nodulation mediated by elevated CO2 concentration. Among the different expression genes (DEGs), we identified a gene encoding small heat shock protein (sHSP) called GmHSP23.9, which mainly expressed in soybean roots and nodules, and its expression was significantly induced by rhizobium USDA110 infection at 14 days after inoculation (DAI) under elevated CO2 conditions. We further investigated the role of GmHSP23.9 by generating transgenic composite plants carrying GmHSP23.9 overexpression (GmHSP23.9-OE), RNA interference (GmHSP23.9-RNAi), and CRISPR-Cas9 (GmHSP23.9-KO), and these modifications resulted in notable changes in nodule number and the root hairs deformation and suggesting that GmHSP23.9 function as an important positive regulator in soybean. Moreover, we found that altering the expression of GmHSP23.9 influenced the expression of genes involved in the Nod factor signaling pathway and AON signaling pathway to modulate soybean nodulation. Interestingly, we found that knocking down of GmHSP23.9 prevented the increase in the nodule number of soybean in response to elevated CO2 concentration. This research has successfully identified a crucial regulator that influences soybean nodulation under elevated CO2 level and shedding new light on the role of sHSPs in legume nodulation.

A Highly Efficient Agrobacterium rhizogenes-Mediated Hairy Root Transformation Method of Idesia polycarpa and the Generation of Transgenic Plants

Idesia polycarpa is a promising woody oilseed species because of its high oil yield. However, its use is greatly limited due to the lack of varieties with good qualities; additionally, gene function has been less studied in this plant because an efficient transformation method has not been established yet. In this study, we established a rapid and efficient hairy root transformation method by infecting the whole seedling, the rootless seedling, and the leaf petiole with Agrobacterium rhizogenes using different infection methods. Among these transformation methods, a higher transformation efficiency was obtained using the whole seedling, which could reach up to 71.91%. Furthermore, we found that the seedling age significantly affected the transformation efficiency, either using whole or rootless seedlings. Additionally, we found that the transgenic roots could regenerate transgenic shoots. Taken together, our study lays the foundation for future study and for genetically modifying wood traits in the future.

GmMYB183, a R2R3-MYB Transcription Factor in Tamba Black Soybean (Glycine max. cv. Tamba), Conferred Aluminum Tolerance in Arabidopsis and Soybean

Aluminum (Al) toxicity is one of the environmental stress factors that affects crop growth, development, and productivity. MYB transcription factors play crucial roles in responding to biotic or abiotic stresses. However, the roles of MYB transcription factors in Al tolerance have not been clearly elucidated. Here, we found that GmMYB183, a gene encoding a R2R3 MYB transcription factor, is involved in Al tolerance. Subcellular localization studies revealed that GmMYB183 protein is located in the nucleus, cytoplasm and cell membrane. Overexpression of GmMYB183 in Arabidopsis and soybean hairy roots enhanced plant tolerance towards Al stress compared to the wild type, with higher citrate secretion and less Al accumulation. Furthermore, we showed that GmMYB183 binds the GmMATE75 gene promoter encoding for a plasma-membrane-localized citrate transporter. Through a dual-luciferase reporter system and yeast one hybrid, the GmMYB183 protein was shown to directly activate the transcription of GmMATE75. Furthermore, the expression of GmMATE75 may depend on phosphorylation of Ser36 residues in GmMYB183 and two MYB sites in P3 segment of the GmMATE75 promoter. In conclusion, GmMYB183 conferred Al tolerance by promoting the secretion of citrate, which provides a scientific basis for further elucidating the mechanism of plant Al resistance.

GmMATE100 Is Involved in the Import of Soyasaponins A and B into Vacuoles in Soybean Plants (Glycine max L.)

Triterpenoid saponins, synthesized via the mevalonic acid (MVA) pathway in the cytoplasm, provide protection against pathogens and pests in plants and health benefits for humans. However, the mechanisms by which triterpenoid saponins are transported between cellular compartments remain uncharacterized. Here, we characterize a tonoplast localized multidrug and toxic compound extrusion transporter, GmMATE100 (encoded by Glyma.18G143700), from soybean (Glycine max L.). GmMATE100 is co-expressed with soyasaponin biosynthetic genes, and its expression was induced by MeJA treatment, which also led to soyasaponin accumulation in soybean roots. GmMATE100 efficiently transports multiple type-B soyasaponins as well as type-A soyasaponins with low affinity from the cytosol to the vacuole in a yeast system. The GmMATE100 loss-of-function mutant showed a significant decrease in type-A and type-B soyasaponin contents in soybean roots. This study not only characterized the first soybean triterpenoid saponin transporter but also provided new knowledge for the rational engineering of soyasaponin content and composition in soybean plants to modulate their levels within crop environments.

Soybean MKK2 establishes intricate signalling pathways to regulate soybean response to cyst nematode infection

Mitogen-activated protein kinase (MPK) cascades play central signalling roles in plant immunity and stress response. The soybean orthologue of MPK kinase2 (GmMKK2) was recently identified as a potential signalling node whose expression is upregulated in the feeding site induced by soybean cyst nematode (SCN, Heterodera glycines). To investigate the role of GmMKK2 in soybean-SCN interactions, we overexpressed a catabolically inactive variant referred to as kinase-dead variant (KD-GmMKK2) using transgenic hairy roots. KD-GmMKK2 overexpression caused significant reduction in soybean susceptibility to SCN, while overexpression of the wild-type variant (WT-GmMKK2) exhibited no effect on susceptibility. Transcriptome analysis indicated that KD-GmMKK2 overexpressing plants are primed for SCN resistance via constitutive activation of defence signalling, particularly those related to chitin, respiratory burst, hydrogen peroxide and salicylic acid. Phosphoproteomic profiling of the WT-GmMKK2 and KD-GmMKK2 root samples upon SCN infection resulted in the identification of 391 potential targets of GmMKK2. These targets are involved in a broad range of biological processes, including defence signalling, vesicle fusion, chromatin remodelling and nuclear organization among others. Furthermore, GmMKK2 mediates phosphorylation of numerous transcriptional and translational regulators, pointing to the presence of signalling shortcuts besides the canonical MAPK cascades to initiate downstream signalling that eventually regulates gene expression and translation initiation. Finally, the functional requirement of specific phosphorylation sites for soybean response to SCN infection was validated by overexpressing phospho-mimic and phospho-dead variants of two differentially phosphorylated proteins SUN1 and IDD4. Together, our analyses identify GmMKK2 impacts on signalling modules that regulate soybean response to SCN infection.

Flowering time regulator qFT13-3 involved in soybean adaptation to high latitudes

Soybean is a short-day plant that typically flowers earlier when exposed to short-day conditions. However, the identification of genes associated with earlier flowering time but without a yield penalty is rare. In this study, we conducted genome-wide association studies (GWAS) using two re-sequencing datasets that included 113 wild soybeans (G. soja) and 1192 cultivated soybeans (G. max), respectively, and simultaneously identified a candidate flowering gene, qFT13-3, which encodes a protein homologous to the pseudo-response regulator (PRR) transcription factor. We identified four major haplotypes of qFT13-3 in the natural population, with haplotype H4 (qFT13-3H4) being lost during domestication, while qFT13-3H1 underwent natural and artificial selection, increasing in proportion from 4.5% in G. soja to 43.8% in landrace and to 81.9% in improve cultivars. Notably, most cultivars harbouring qFT13-3H1 were located in high-latitude regions. Knockout of qFT13-3 accelerated flowering and maturity time under long-day conditions, indicating that qFT13-3 functions as a flowering inhibitor. Our results also showed that qFT13-3 directly downregulates the expression of GmELF3b-2 which is a component of the circadian clock evening complex. Field trials revealed that the qft13-3 mutants shorten the maturity period by 11 days without a concomitant penalty on yield. Collectively, qFT13-3 can be utilized for the breeding of high-yield cultivars with a short maturity time suitable for high latitudes.

The MYB Transcription Factor GmMYB78 Negatively Regulates Phytophthora sojae Resistance in Soybean

Phytophthora root rot is a devastating disease of soybean caused by Phytophthora sojae. However, the resistance mechanism is not yet clear. Our previous studies have shown that GmAP2 enhances sensitivity to P. sojae in soybean, and GmMYB78 is downregulated in the transcriptome analysis of GmAP2-overexpressing transgenic hairy roots. Here, GmMYB78 was significantly induced by P. sojae in susceptible soybean, and the overexpressing of GmMYB78 enhanced sensitivity to the pathogen, while silencing GmMYB78 enhances resistance to P. sojae, indicating that GmMYB78 is a negative regulator of P. sojae. Moreover, the jasmonic acid (JA) content and JA synthesis gene GmAOS1 was highly upregulated in GmMYB78-silencing roots and highly downregulated in overexpressing ones, suggesting that GmMYB78 could respond to P. sojae through the JA signaling pathway. Furthermore, the expression of several pathogenesis-related genes was significantly lower in GmMYB78-overexpressing roots and higher in GmMYB78-silencing ones. Additionally, we screened and identified the upstream regulator GmbHLH122 and downstream target gene GmbZIP25 of GmMYB78. GmbHLH122 was highly induced by P. sojae and could inhibit GmMYB78 expression in resistant soybean, and GmMYB78 was highly expressed to activate downstream target gene GmbZIP25 transcription in susceptible soybean. In conclusion, our data reveal that GmMYB78 triggers soybean sensitivity to P. sojae by inhibiting the JA signaling pathway and the expression of pathogenesis-related genes or through the effects of the GmbHLH122-GmMYB78-GmbZIP25 cascade pathway.

The EIN3 transcription factor GmEIL1 improves soybean resistance to Phytophthora sojae

Phytophthora root and stem rot of soybean (Glycine max), caused by the oomycete Phytophthora sojae, is an extremely destructive disease worldwide. In this study, we identified GmEIL1, which encodes an ethylene-insensitive3 (EIN3) transcription factor. GmEIL1 was significantly induced following P. sojae infection of soybean plants. Compared to wild-type soybean plants, transgenic soybean plants overexpressing GmEIL1 showed enhanced resistance to P. sojae and GmEIL1-silenced RNA-interference lines showed more severe symptoms when infected with P. sojae. We screened for target genes of GmEIL1 and confirmed that GmEIL1 bound directly to the GmERF113 promoter and regulated GmERF113 expression. Moreover, GmEIL1 positively regulated the expression of the pathogenesis-related gene GmPR1. The GmEIL1-regulated defence response to P. sojae involved both ethylene biosynthesis and the ethylene signalling pathway. These findings suggest that the GmEIL1-GmERF113 module plays an important role in P. sojae resistance via the ethylene signalling pathway.

Tandemly duplicated MYB genes are functionally diverged in the regulation of anthocyanin biosynthesis in soybean

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean (Glycine max) genome is characterized by a large number of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that 4 MYB genes (GmMYBA5, GmMYBA2, GmMYBA1, and Glyma.09g235000)-presumably generated by tandem duplication specifically in the Phaseoleae lineage-exhibited a stronger purifying selection in soybean compared to common bean (Phaseolus vulgaris). To gain insights into the diverse functions of these tandemly duplicated MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, induced metabolites, and evolutionary history of these 4 MYB genes. Our data revealed that Glyma.09g235000 is a pseudogene, while the remaining 3 MYB genes exhibit strong transcriptional activation activity, promoting anthocyanin biosynthesis in different soybean tissues. GmMYBA5, GmMYBA2, and GmMYBA1 induced anthocyanin accumulation by upregulating the expression of anthocyanin pathway-related genes. Notably, GmMYBA5 showed a lower capacity for gene induction compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated that GmMYBA5 induced distinct anthocyanin accumulation in Nicotiana benthamiana leaves and soybean hairy roots compared to GmMYBA2 and GmMYBA1, suggesting their functional divergence leading to the accumulation of different metabolites accumulation following gene duplication. Together, our data provide evidence of functional divergence within the MYB gene cluster following tandem duplication, which sheds light on the potential evolutionary directions of gene duplications during legume evolution.

Identification of major QTLs for drought tolerance in soybean, together with a novel candidate gene, GmUAA6

Drought tolerance is a complex trait in soybean that is controlled by polygenetic quantitative trait loci (QTLs). In this study, wilting score, days-to-wilting, leaf relative water content, and leaf relative conductivity were used to identify QTLs associated with drought tolerance in recombinant inbred lines derived from a cross between a drought-sensitive variety, Lin, and a drought-tolerant variety, Meng. A total of 33 drought-tolerance QTLs were detected. Of these 17 were major QTLs. In addition, 15 were novel drought-tolerance QTLs. The most predominant QTL was on chromosome 11. This was detected in at least three environments. The overlapped mapping interval of the four measured traits was 0.2 cM in genetic distance (about 220 kb in physical length). Glyma.11g143500 (designated as GmUAA6), which encodes a UDP-N-acetylglucosamine transporter, was identified as the most likely candidate gene. The allele of GmUAA6 from Lin (GmUAA6Lin) was associated with improved soybean drought tolerance. Overexpression of GmUAA6Lin in Arabidopsis and soybean hairy roots enhanced drought tolerance. Furthermore, a 3-bp insertion/deletion (InDel) in the coding sequence of GmUAA6 explained up to 49.9% of the phenotypic variation in drought tolerance-related traits, suggesting that this InDel might be used in future marker-assisted selection of drought-tolerant lines in soybean breeding programs.

An engineered Cas12i nuclease that is an efficient genome editing tool in animals and plants

The type V-I CRISPR-Cas system is becoming increasingly more attractive for genome editing. However, natural nucleases of this system often exhibit low efficiency, limiting their application. Here, we used structure-guided rational design and protein engineering to optimize an uncharacterized Cas12i nuclease, Cas12i3. As a result, we developed Cas-SF01, a Cas12i3 variant that exhibits significantly improved gene editing activity in mammalian cells. Cas-SF01 shows comparable or superior editing performance compared to SpCas9 and other Cas12 nucleases. Compared to natural Cas12i3, Cas-SF01 has an expanded PAM range and effectively recognizes NTTN and noncanonical NATN and TTVN PAMs. In addition, we identified an amino acid substitution, D876R, that markedly reduced the off-target effect while maintaining high on-target activity, leading to the development of Cas-SF01HiFi (high-fidelity Cas-SF01). Finally, we show that Cas-SF01 has high gene editing activities in mice and plants. Our results suggest that Cas-SF01 can serve as a robust gene editing platform with high efficiency and specificity for genome editing applications in various organisms.

The Glycine soja cytochrome P450 gene GsCYP82C4 confers alkaline tolerance by promoting reactive oxygen species scavenging

Recent studies have demonstrated the crucial role of Cytochrome P450 enzymes (CYPs) in the production of secondary metabolites, phytohormones and antioxidants in plants. However, their functional characterization specifically under alkaline stress remains elusive. CYP82C4 was the key gene screened from a family of wild soybean CYPs in our previous studies. The aim of this present study was to clone the Glycine soja GsCYP82C4 gene and characterize its functions in Arabidopsis and Glycine max. The results showed that the GsCYP82C4 gene displayed a high expression in different plant tissues at mature stages compared to young stages. Further, higher temporal expression of the GsCYP82C4 gene was noted at 6, 12 and 24 h time points after alkali treatment in leaves compared to roots. In addition, overexpression of GsCYP82C4 improved alkaline stress tolerance in Arabidopsis via increased root lengths and fresh biomass and strengthened the antioxidant defense system via a reduction in superoxide radicals in transgenic lines compared to wild type (WT) and atcyp82c4 mutants. Further, the expression levels of stress-related marker genes were up-regulated in GsCYP82C4 OX lines under alkali stress. The functional analysis of GsCYP82C4 overexpression in soybean displayed better hairy root growth, increased fresh weight, higher antioxidant enzyme activities and reduced lipid peroxidation rates in OX lines compared to the soybean WT (K599) line. In total, our study displayed positive roles of GsCYP82C4 overexpression in both Arabidopsis and Glycine max to alleviate alkaline stress via altering expression abundance of stress responsive genes, stronger roots, higher antioxidant enzyme activities as well as reduced rates of lipid peroxidation and superoxide radicals.

Identification of GmPT proteins and investigation of their expressions in response to abiotic stress in soybean

MAIN CONCLUSION

Investigation the expression patterns of GmPT genes in response to various abiotic stresses and overexpression of GmPT11 in soybean hairy roots and Arabidopsis exhibited hypersensitivity to salt stress. Soybean is considered to be one of the significant oil crops globally, as it offers a diverse range of essential nutrients that contribute to human health. Salt stress seriously affects the yield of soybean through negative impacts on the growth, nodulation, reproduction, and other agronomy traits. The phosphate transporters 1(PHT1) subfamily, which is a part of the PHTs family in plants, is primarily found in the cell membrane and responsible for the uptake and transport of phosphorus. However, the role of GmPT (GmPT1-GmPT14) genes in response to salt stress has not been comprehensively studied. Here, we conducted a systematic analysis to ascertain the distribution and genomic duplications of GmPT genes, as well as their expression patterns in response to various abiotic stresses. Promoter analysis of GmPT genes revealed that six stress-related cis-elements were enriched in these genes. The overexpression of GmPT11 in soybean hairy roots and Arabidopsis exhibited hypersensitivity to salt stress, while no significant change was observed under low phosphate treatment, suggesting a crucial role in the response to salt stress. These findings provide novel insights into enhancing plant tolerance to salt stress.

Alfalfa MsbHLH115 confers tolerance to cadmium stress through activating the iron deficiency response in Arabidopsis thaliana

Cadmium (Cd) pollution severely affects plant growth and development, posing risks to human health throughout the food chain. Improved iron (Fe) nutrients could mitigate Cd toxicity in plants, but the regulatory network involving Cd and Fe interplay remains unresolved. Here, a transcription factor gene of alfalfa, MsbHLH115 was verified to respond to iron deficiency and Cd stress. Overexpression of MsbHLH115 enhanced tolerance to Cd stress, showing better growth and less ROS accumulation in Arabidopsis thaliana. Overexpression of MsbHLH115 significantly enhanced Fe and Zn accumulation and did not affect Cd, Mn, and Cu concentration in Arabidopsis. Further investigations revealed that MsbHLH115 up-regulated iron homeostasis regulation genes, ROS-related genes, and metal chelation and detoxification genes, contributing to attenuating Cd toxicity. Y1H, EMSA, and LUC assays confirmed the physical interaction between MsbHLH115 and E-box, which is present in the promoter regions of most of the above-mentioned iron homeostasis regulatory genes. The transient expression experiment showed that MsbHLH115 interacted with MsbHLH121pro. The results suggest that MsbHLH115 may directly regulate the iron-deficiency response system and indirectly regulate the metal detoxification response mechanism, thereby enhancing plant Cd tolerance. In summary, enhancing iron accumulation through transcription factor regulation holds promise for improving plant tolerance to Cd toxicity, and MsbHLH115 is a potential candidate for addressing Cd toxicity issues.

A novel natural variation in the promoter of GmCHX1 regulates conditional gene expression to improve salt tolerance in soybean

Identification and characterization of soybean germplasm and gene(s)/allele(s) for salt tolerance is an effective way to develop improved varieties for saline soils. Previous studies identified GmCHX1 (Glyma03g32900) as a major salt tolerance gene in soybean, and two main functional variations were found in the promoter region (148/150 bp insertion) and the third exon with a retrotransposon insertion (3.78 kb). In the current study, we identified four salt-tolerant soybean lines, including PI 483460B (Glycine soja), carrying the previously identified salt-sensitive variations at GmCHX1, suggesting new gene(s) or new functional allele(s) of GmCHX1 in these soybean lines. Subsequently, we conducted quantitative trait locus (QTL) mapping in a recombinant-inbred line population (Williams 82 (salt-sensitive) × PI 483460B) to identify the new salt tolerance loci/alleles. A new locus, qSalt_Gm18, was mapped on chromosome 18 associated with leaf scorch score. Another major QTL, qSalt_Gm03, was identified to be associated with chlorophyll content ratio and leaf scorch score in the same chromosomal region of GmCHX1 on chromosome 3. Novel variations in a STRE (stress response element) cis-element in the promoter region of GmCHX1 were found to regulate the salt-inducible expression of the gene in these four newly identified salt-tolerant lines including PI 483460B. This new allele of GmCHX1 with salt-inducible expression pattern provides an energy cost efficient (conditional gene expression) strategy to protect soybean yield in saline soils without yield penalty under non-stress conditions. Our results suggest that there might be no other major salt tolerance locus similar to GmCHX1 in soybean germplasm, and further improvement of salt tolerance in soybean may rely on gene-editing techniques instead of looking for natural variations.

The mechanism of low blue light-induced leaf senescence mediated by GmCRY1s in soybean

Leaf senescence is a crucial trait that has a significant impact on crop quality and yield. Previous studies have demonstrated that light is a key factor in modulating the senescence process. However, the precise mechanism by which plants sense light and control senescence remains largely unknown, particularly in crop species. In this study, we reveal that the reduction in blue light under shading conditions can efficiently induce leaf senescence in soybean. The blue light receptors GmCRY1s rather than GmCRY2s, primarily regulate leaf senescence in response to blue light signals. Our results show that GmCRY1s interact with DELLA proteins under light-activated conditions, stabilizing them and consequently suppressing the transcription of GmWRKY100 to delay senescence. Conversely, LBL reduces the interaction between GmCRY1s and the DELLA proteins, leading to their degradation and premature senescence of leaves. Our findings suggest a GmCRY1s-GmDELLAs-GmWRKY100 regulatory cascade that is involved in mediating LBL-induced leaf senescence in soybean, providing insight into the mechanism of how light signals regulate leaf senescence. Additionally, we generate GmWRKY100 knockout soybeans that show delayed leaf senescence and improved yield under natural field conditions, indicating potential applications in enhancing soybean production by manipulating the leaf senescence trait.

Agrobacterium rhizogenes-mediated marker-free transformation and gene editing system revealed that AeCBL3 mediates the formation of calcium oxalate crystal in kiwifruit

The transformation and gene editing of the woody species kiwifruit are difficult and time-consuming. The fast and marker-free genetic modification system for kiwifruit has not been developed yet. Here, we establish a rapid and efficient marker-free transformation and gene editing system mediated by Agrobacterium rhizogenes for kiwifruit. Moreover, a removing-root-tip method was developed to significantly increase the regeneration efficiency of transgenic hairy roots. Through A. rhizogenes-mediated CRISPR/Cas9 gene editing, the editing efficiencies of CEN4 and AeCBL3 achieved 55 and 50%, respectively. And several homozygous knockout lines for both genes were obtained. Our method has been successfully applied in the transformation of two different species of kiwifruit (Actinidia chinensis 'Hongyang' and A.eriantha 'White'). Next, we used the method to study the formation of calcium oxalate (CaOx) crystals in kiwifruit. To date, little is known about how CaOx crystal is formed in plants. Our results indicated that AeCBL3 overexpression enhanced CaOx crystal formation, but its knockout via CRISPR/Cas9 significantly impaired crystal formation in kiwifruit. Together, we developed a fast maker-free transformation and highly efficient CRISPR-Cas9 gene editing system for kiwifruit. Moreover, our work revealed a novel gene mediating CaOx crystal formation and provided a clue to elaborate the underlying mechanisms.

Construction of a high-density genetic map for yardlong bean and identification of ANT1 as a regulator of anthocyanin biosynthesis

Because its long, tender pods supply essential proteins, vitamins, and fibers to humans, yardlong bean (Vigna unguiculata ssp. sesquipedalis) is a commonly consumed vegetable, especially in Southeast Asia. To provide insights into the genetic bases of key agricultural traits in yardlong bean, we here created a high-density bin-map with 2084 bin markers using 514 227 SNPs from a recombinant-inbred line (RIL) population. Quantitative trait loci (QTL) mapping was carried out to identify loci associated with anthocyanin content (ANT), vitamin E content (VE), total soluble protein content (TSP), pod length (PL), hundred-seed weight (HSW), seed length and width (SL and SW, respectively), and seed coat color (SCC). In total, 20 related QTLs were isolated, explaining 7.58-56.03% of the phenotypic variation. Of these, five major QTLs (qANT5, qTSP11, qVE7, qPL3, and qSCC9) were detected in 2020, 2021, and the combined environment, explaining 11.96-56.03% of the phenotypic variation. VuANT1 was identified as a causal gene for the QTL qANT5, which regulated anthocyanin content; VuANT1 was highly expressed in immature purple pods but barely detectable in white pods. VuANT1 overexpression in tobacco leaves and yardlong bean hairy roots resulted in purple coloration as a result of anthocyanin accumulation. These findings suggested that VuANT1 was a key regulator of anthocyanin accumulation in yardlong bean. Our results lay a firm foundation for target agricultural trait improvement and clarification of the genetic mechanisms underlying agricultural traits in yardlong bean.

Tracking organelle activities through efficient and stable root genetic transformation system in woody plants

Due to the protracted transgenic timeline and low efficiency in stable genetic transformation of woody plants, there has been limited exploration of real-time organelle imaging within stable transgenic woody plant cells. Here, we established an efficient in vivo genetic transformation system for woody plants using an Agrobacterium rhizogenes-mediated approach. This system was successfully validated in multiple perennial woody species. Using citrus as a model, we introduced organelle-targeted fluorescent reporters via genetic transformation and investigated their subcellular localization and dynamics using advanced imaging techniques, such as confocal microscopy and live-cell imaging. Moreover, we subjected transgenic MT-GFP-labeled mitochondria in root cells to stress conditions simulating agricultural adversities faced by fruit crops. The stress-induced experiments revealed notable alterations in mitochondrial morphology. Our study contributes novel insights into membrane trafficking processes, protein localization dynamics, and cellular physiology in woody plants, while also providing stable and efficient genetic transformation methods for perennial woody species.

Technology Invention and Mechanism Analysis of Rapid Rooting of Taxus × media Rehder Branches Induced by Agrobacterium rhizogenes

Taxus, a vital source of the anticancer drug paclitaxel, grapples with a pronounced supply-demand gap. Current efforts to alleviate the paclitaxel shortage involve expanding Taxus cultivation through cutting propagation. However, traditional cutting propagation of Taxus is difficult to root and time-consuming. Obtaining the roots with high paclitaxel content will cause tree death and resource destruction, which is not conducive to the development of the Taxus industry. To address this, establishing rapid and efficient stem rooting systems emerges as a key solution for Taxus propagation, facilitating direct and continuous root utilization. In this study, Agrobacterium rhizogenes were induced in the 1-3-year-old branches of Taxus × media Rehder, which has the highest paclitaxel content. The research delves into the rooting efficiency induced by different A. rhizogenes strains, with MSU440 and C58 exhibiting superior effects. Transcriptome and metabolome analyses revealed A. rhizogenes' impact on hormone signal transduction, amino acid metabolism, zeatin synthesis, and secondary metabolite synthesis pathways in roots. LC-MS-targeted quantitative detection showed no significant difference in paclitaxel and baccatin III content between naturally formed and induced roots. These findings underpin the theoretical framework for T. media rapid propagation, contributing to the sustainable advancement of the Taxus industry.

H2O2-dependent oxidation of the transcription factor GmNTL1 promotes salt tolerance in soybean

Reactive oxygen species (ROS) play an essential role in plant growth and responses to environmental stresses. Plant cells sense and transduce ROS signaling directly via hydrogen peroxide (H2O2)-mediated posttranslational modifications (PTMs) on protein cysteine residues. Here, we show that the H2O2-mediated cysteine oxidation of NAC WITH TRANS-MEMBRANE MOTIF1-LIKE 1 (GmNTL1) in soybean (Glycine max) during salt stress promotes its release from the endoplasmic reticulum (ER) membrane and translocation to the nucleus. We further show that an oxidative posttranslational modification on GmNTL1 residue Cys-247 steers downstream amplification of ROS production by binding to and activating the promoters of RESPIRATORY BURST OXIDASE HOMOLOG B (GmRbohB) genes, thereby creating a feed-forward loop to fine-tune GmNTL1 activity. In addition, oxidation of GmNTL1 Cys-247 directly promotes the expression of CATION H+ EXCHANGER 1 (GmCHX1)/SALT TOLERANCE-ASSOCIATED GENE ON CHROMOSOME 3 (GmSALT3) and Na+/H+ Antiporter 1 (GmNHX1). Accordingly, transgenic overexpression of GmNTL1 in soybean increases the H2O2 levels and K+/Na+ ratio in the cell, promotes salt tolerance, and increases yield under salt stress, while an RNA interference-mediated knockdown of GmNTL1 elicits the opposite effects. Our results reveal that the salt-induced oxidation of GmNTL1 promotes its relocation and transcriptional activity through an H2O2-mediated posttranslational modification on cysteine that improves resilience of soybean against salt stress.

An ex vitro hairy root system from petioles of detached soybean leaves for in planta screening of target genes and CRISPR strategies associated with nematode bioassays

MAIN CONCLUSION

The ex vitro hairy root system from petioles of detached soybean leaves allows the functional validation of genes using classical transgenesis and CRISPR strategies (e.g., sgRNA validation, gene activation) associated with nematode bioassays. Agrobacterium rhizogenes-mediated root transformation has been widely used in soybean for the functional validation of target genes in classical transgenesis and single-guide RNA (sgRNA) in CRISPR-based technologies. Initial data showed that in vitro hairy root induction from soybean cotyledons and hypocotyls were not the most suitable strategies for simultaneous performing genetic studies and nematode bioassays. Therefore, an ex vitro hairy root system was developed for in planta screening of target molecules during soybean parasitism by root-knot nematodes (RKNs). Applying this method, hairy roots were successfully induced by A. rhizogenes from petioles of detached soybean leaves. The soybean GmPR10 and GmGST genes were then constitutively overexpressed in both soybean hairy roots and tobacco plants, showing a reduction in the number of Meloidogyne incognita-induced galls of up to 41% and 39%, respectively. In addition, this system was evaluated for upregulation of the endogenous GmExpA and GmExpLB genes by CRISPR/dCas9, showing high levels of gene activation and reductions in gall number of up to 58.7% and 67.4%, respectively. Furthermore, morphological and histological analyses of the galls were successfully performed. These collective data validate the ex vitro hairy root system for screening target genes, using classical overexpression and CRISPR approaches, directly in soybean in a simple manner and associated with nematode bioassays. This system can also be used in other root pathosystems for analyses of gene function and studies of parasite interactions with plants, as well as for other purposes such as studies of root biology and promoter characterization.

NopC/T/L Signal Crosstalk Gene GmPHT1-4

Symbiotic nodulation between leguminous plants and rhizobia is a critical biological interaction. The type III secretion system (T3SS) employed by rhizobia manipulates the host's nodulation signaling, analogous to mechanisms used by certain bacterial pathogens for effector protein delivery into host cells. This investigation explores the interactive signaling among type III effectors HH103ΩNopC, HH103ΩNopT, and HH103ΩNopL from SinoRhizobium fredii HH103. Experimental results revealed that these effectors positively regulate nodule formation. Transcriptomic analysis pinpointed GmPHT1-4 as the key gene facilitating this effector-mediated signaling. Overexpression of GmPHT1-4 enhances nodulation, indicating a dual function in nodulation and phosphorus homeostasis. This research elucidates the intricate regulatory network governing Rhizobium-soybean (Glycine max (L.) Merr) interactions and the complex interplay between type III effectors.

Agrobacterium rhizogenes: paving the road to research and breeding for woody plants

Woody plants play a vital role in global ecosystems and serve as valuable resources for various industries and human needs. While many woody plant genomes have been fully sequenced, gene function research and biotechnological breeding advances have lagged behind. As a result, only a limited number of genes have been elucidated, making it difficult to use newer tools such as CRISPR-Cas9 for biotechnological breeding purposes. The use of Agrobacterium rhizogenes as a transformative tool in plant biotechnology has received considerable attention in recent years, particularly in the research field on woody plants. Over the past three decades, numerous woody plants have been effectively transformed using A. rhizogenes-mediated techniques. Some of these transformed plants have successfully regenerated. Recent research on A. rhizogenes-mediated transformation of woody plants has demonstrated its potential for various applications, including gene function analysis, gene expression profiling, gene interaction studies, and gene regulation analysis. The introduction of the Ri plasmid has resulted in the emergence of several Ri phenotypes, such as compact plant types, which can be exploited for Ri breeding purposes. This review paper presents recent advances in A. rhizogenes-mediated basic research and Ri breeding in woody plants. This study highlights various aspects of A. rhizogenes-mediated transformation, its multiple applications in gene function analysis, and the potential of Ri lines as valuable breeding materials.

A recently evolved BAHD acetyltransferase, responsible for bitter soyasaponin A production, is indispensable for soybean seed germination

Soyasaponins are major small molecules that accumulate in soybean (Glycine max) seeds. Among them, type-A soyasaponins, fully acetylated at the terminal sugar of their C22 sugar chain, are responsible for the bitter taste of soybean-derived foods. However, the molecular basis for the acetylation of type-A soyasaponins remains unclear. Here, we identify and characterize GmSSAcT1, encoding a BADH-type soyasaponin acetyltransferase that catalyzes three or four consecutive acetylations on type-A soyasaponins in vitro and in planta. Phylogenetic analysis and biochemical assays suggest that GmSSAcT1 likely evolved from acyltransferases present in leguminous plants involved in isoflavonoid acylation. Loss-of-function mutants of GmSSAcT1 exhibited impaired seed germination, which attribute to the excessive accumulation of null-acetylated type-A soyasaponins. We conclude that GmSSAcT1 not only functions as a detoxification gene for high accumulation of type-A soyasaponins in soybean seeds but is also a promising target for breeding new soybean varieties with lower bitter soyasaponin content.

MicroRNA 4407 modulates nodulation in soybean by repressing a root-specific ISOPENTENYLTRANSFERASE (GmIPT3)

MicroRNAs (miRNAs) are important regulators of plant biological processes, including soybean nodulation. One miRNA, miR4407, was identified in soybean roots and nodules. However, the function of miR4407 in soybean is still unknown. MiR4407, unique to soybean, positively regulates lateral root emergence and root structures and represses a root-specific ISOPENTENYLTRANSFERASE (GmIPT3). By altering the expression of miR4407 and GmIPT3, we investigated the role of miR4407 in lateral root and nodule development. Both miR4407 and GmIPT3 are expressed in the inner root cortex and nodule primordia. Upon rhizobial inoculation, miR4407 was downregulated while GmIPT3 was upregulated. Overexpressing miR4407 reduced the number of nodules in transgenic soybean hairy roots while overexpressing the wild-type GmIPT3 or a miR4407-resistant GmIPT3 mutant (mGmIPT3) significantly increased the nodule number. The mechanism of miR4407 and GmIPT3 functions was also linked to autoregulation of nodulation (AON), where miR4407 overexpression repressed miR172c and activated its target, GmNNC1, turning on AON. Exogenous CK mimicked the effects of GmIPT3 overexpression on miR172c, supporting the notion that GmIPT3 regulates nodulation by enhancing root-derived CK. Overall, our data revealed a new miRNA-mediated regulatory mechanism of nodulation in soybean. MiR4407 showed a dual role in lateral root and nodule development.

Highly efficient Agrobacterium rhizogenes-mediated hairy root transformation in citrus seeds and its application in gene functional analysis

Highly efficient genetic transformation technology is beneficial for plant gene functional research and molecular improvement breeding. However, the most commonly used Agrobacterium tumefaciens-mediated genetic transformation technology is time-consuming and recalcitrant for some woody plants such as citrus, hampering the high-throughput functional analysis of citrus genes. Thus, we dedicated to develop a rapid, simple, and highly efficient hairy root transformation system induced by Agrobacterium rhizogenes to analyze citrus gene function. In this report, a rapid, universal, and highly efficient hairy root transformation system in citrus seeds was described. Only 15 days were required for the entire workflow and the system was applicable for various citrus genotypes, with a maximum transformation frequency of 96.1%. After optimization, the transformation frequency of Citrus sinensis, which shows the lowest transformation frequency of 52.3% among four citrus genotypes initially, was increased to 71.4% successfully. To test the applicability of the hairy roots transformation system for gene functional analysis of citrus genes, we evaluated the subcellular localization, gene overexpression and gene editing in transformed hairy roots. Compared with the traditional transient transformation system performed in tobacco leaves, the transgenic citrus hairy roots displayed a more clear and specific subcellular fluorescence localization. Transcript levels of genes were significantly increased in overexpressing transgenic citrus hairy roots as compared with wild-type (WT). Additionally, hairy root transformation system in citrus seeds was successful in obtaining transformants with knocked out targets, indicating that the Agrobacterium rhizogenes-mediated transformation enables the CRISPR/Cas9-mediated gene editing. In summary, we established a highly efficient genetic transformation technology with non-tissue-culture in citrus that can be used for functional analysis such as protein subcellular localization, gene overexpression and gene editing. Since the material used for genetic transformation are roots protruding out of citrus seeds, the process of planting seedlings prior to transformation of conventional tissue culture or non-tissue-culture was eliminated, and the experimental time was greatly reduced. We anticipate that this genetic transformation technology will be a valuable tool for routine research of citrus genes in the future.

Genome-Wide Identification of the CDPK Gene Family and Their Involvement in Taproot Cracking in Radish

Taproot cracking, a severe and common physiological disorder, markedly reduces radish yield and commercial value. Calcium-dependent protein kinase (CDPK) plays a pivotal role in various plant developmental processes; however, its function in radish taproot cracking remains largely unknown. Here, 37 RsCDPK gene members were identified from the long-read radish genome "QZ-16". Phylogenetic analysis revealed that the CDPK members in radish, tomato, and Arabidopsis were clustered into four groups. Additionally, synteny analysis identified 13 segmental duplication events in the RsCDPK genes. Analysis of paraffin-embedded sections showed that the density and arrangement of fleshy taproot cortex cells are important factors that affect radish cracking. Transcriptome sequencing of the fleshy taproot cortex revealed 5755 differentially expressed genes (DEGs) (3252 upregulated and 2503 downregulated) between non-cracking radish "HongYun" and cracking radish "505". These DEGs were significantly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interaction KEGG pathways. Furthermore, when comparing the 37 RsCDPK gene family members and RNA-seq DEGs, we identified six RsCDPK genes related to taproot cracking in radish. Soybean hairy root transformation experiments showed that RsCDPK21 significantly and positively regulates root length development. These findings provide valuable insights into the relationship between radish taproot cracking and RsCDPK gene function.

The Bax inhibitor GmBI-1α interacts with a Nod factor receptor and plays a dual role in the legume-rhizobia symbiosis

The gene networks surrounding Nod factor receptors that govern the symbiotic process between legumes and rhizobia remain largely unexplored. Here, we identify 13 novel GmNFR1α-associated proteins by yeast two-hybrid screening, and describe a potential interacting protein, GmBI-1α. GmBI-1α had the highest positive correlation with GmNFR1α in a co-expression network analysis, and its expression at the mRNA level in roots was enhanced by rhizobial infection. Moreover, GmBI-1α-GmNFR1α interaction was shown to occur in vitro and in vivo. The GmBI-1α protein was localized to multiple subcellular locations, including the endoplasmic reticulum and plasma membrane. Overexpression of GmBI-1α increased the nodule number in transgenic hairy roots or transgenic soybean, whereas down-regulation of GmBI-1α transcripts by RNA interference reduced the nodule number. In addition, the nodules in GmBI-1α-overexpressing plants became smaller in size and infected area with reduced nitrogenase activity. In GmBI-1α-overexpressing transgenic soybean, the elevated GmBI-1α also promoted plant growth and suppressed the expression of defense signaling-related genes. Infection thread analysis of GmBI-1α-overexpressing plants showed that GmBI-1α promoted rhizobial infection. Collectively, our findings support a GmNFR1α-associated protein in the Nod factor signaling pathway and shed new light on the regulatory mechanism of GmNFR1α in rhizobial symbiosis.

Highly efficient Agrobacterium rhizogenes-mediated transformation for functional analysis in woodland strawberry

BACKGROUND

The diploid woodland strawberry (Fragaria vesca) is an excellent model plant for investigating economically significant traits and several genetic resources within the Rosaceae family. Agrobacterium rhizogenes-mediated hairy root transformation is an alternative for exploring gene functions, especially the genes specifically expressed in roots. However, the hairy root transformation has not been established in strawberry.

RESULTS

Here, we described an efficient and rapid hairy root transgenic system for strawberry using A. rhizogenes. Strain of A. rhizogenes MSU440 or C58C1 was the most suitable for hairy root transformation. The transformation efficiency was highest when tissues contained hypocotyls as explants. The optimal procedure involves A. rhizogenes at an optical density (OD600) of 0.7 for 10 min and co-cultivation duration for four days, achieving a transgenic efficiency of up to 71.43%. An auxin responsive promoter DR5ver2 carrying an enhanced green fluorescent protein (eGFP) marker was transformed by A. rhizogenes MSU440, thereby generating transgenic hairy roots capable of high eGFP expression in root tip and meristem of strawberry where auxin accumulated. Finally, this system was applied for functional analysis using jGCaMP7c, which could sense calcium signals. A significant upsurge in eGFP expression in the transgenic hairy roots was displayed after adding calcium chloride. The results suggested that this approach was feasible for studying specific promoters and could be a tool to analyze gene functions in the roots of strawberries.

CONCLUSION

We established a rapid and efficient hairy root transformation in strawberry by optimizing parameters, which was adequate for promoter analysis and functional characterization of candidate genes in strawberry and other rosaceous plants.

The soybean immune receptor GmBIR1 regulates host transcriptome, spliceome, and immunity during cyst nematode infection

BAK1-INTERACTING RECEPTOR LIKE KINASE1 (BIR1) is a negative regulator of various aspects of disease resistance and immune responses. Here, we investigated the functional role of soybean (Glycine max) BIR1 (GmBIR1) during soybean interaction with soybean cyst nematode (SCN, Heterodera glycines) and the molecular mechanism through which GmBIR1 regulates plant immunity. Overexpression of wild-type variant of GmBIR1 (WT-GmBIR1) using transgenic soybean hairy roots significantly increased soybean susceptibility to SCN, whereas overexpression of kinase-dead variant (KD-GmBIR1) significantly increased plant resistance. Transcriptome analysis revealed that genes oppositely regulated in WT-GmBIR1 and KD-GmBIR1 upon SCN infection were enriched primarily in defense and immunity-related functions. Quantitative phosphoproteomic analysis identified 208 proteins as putative substrates of the GmBIR1 signaling pathway, 114 of which were differentially phosphorylated upon SCN infection. In addition, the phosphoproteomic data pointed to a role of the GmBIR1 signaling pathway in regulating alternative pre-mRNA splicing. Genome-wide analysis of splicing events provided compelling evidence supporting a role of the GmBIR1 signaling pathway in establishing alternative splicing during SCN infection. Our results provide novel mechanistic insights into the function of the GmBIR1 signaling pathway in regulating soybean transcriptome and spliceome via differential phosphorylation of splicing factors and regulation of splicing events of pre-mRNA decay- and spliceosome-related genes.

GsPKS24, a calcineurin B-like protein-interacting protein kinase gene from Glycine soja, positively regulates tolerance to pH stress and ABA signal transduction

SOS2-like protein kinases (PKS/CIPK) family genes are known to be involved in various abiotic stresses in plants. Even though, its functions have been well characterized under salt and drought stresses. The roles of PKS genes associated with alkaline stress response are not fully established yet. In this study, we identified 56 PKS family genes which could be mainly classified into three groups in wild soybean (Glycine soja). PKS family genes transcript profiles revealed different expression patterns under alkali stress. Furthermore, we confirmed the regulatory roles of GsPKS24 in response to NaHCO3, pH and ABA treatments. Overexpression of GsPKS24 enhanced plant tolerance to pH stress in Arabidopsis and soybean hairy roots but conferred suppressed pH tolerance in Arabidopsis atpks mutant. Additionally, Overexpression of GsPKS24 decreased the ABA sensitivity compared to Arabidopsis atpks mutant which displayed more sensitivity towards ABA. Moreover, upregulated expression of stress responsive and ABA signal-related genes were detected in GsPKS24 overexpression lines. In conclusion, we identified the wild soybean PKS family genes, and explored the roles of GsPKS24 in positive response to pH stress tolerance, and in alleviation of ABA sensitivity.

GmNAC039 and GmNAC018 activate the expression of cysteine protease genes to promote soybean nodule senescence

Root nodules are major sources of nitrogen for soybean (Glycine max (L.) Merr.) growth, development, production, and seed quality. Symbiotic nitrogen fixation is time-limited, as the root nodule senesces during the reproductive stage of plant development, specifically during seed development. Nodule senescence is characterized by the induction of senescence-related genes, such as papain-like cysteine proteases (CYPs), which ultimately leads to the degradation of both bacteroids and plant cells. However, how nodule senescence-related genes are activated in soybean is unknown. Here, we identified 2 paralogous NAC transcription factors, GmNAC039 and GmNAC018, as master regulators of nodule senescence. Overexpression of either gene induced soybean nodule senescence with increased cell death as detected using a TUNEL assay, whereas their knockout delayed senescence and increased nitrogenase activity. Transcriptome analysis and nCUT&Tag-qPCR assays revealed that GmNAC039 directly binds to the core motif CAC(A)A and activates the expression of 4 GmCYP genes (GmCYP35, GmCYP37, GmCYP39, and GmCYP45). Similar to GmNAC039 and GmNAC018, overexpression or knockout of GmCYP genes in nodules resulted in precocious or delayed senescence, respectively. These data provide essential insights into the regulatory mechanisms of nodule senescence, in which GmNAC039 and GmNAC018 directly activate the expression of GmCYP genes to promote nodule senescence.

Optimization of in vitro and ex vitro Agrobacterium rhizogenes-mediated hairy root transformation of soybean for visual screening of transformants using RUBY

In vitro and ex vitro Agrobacterium rhizogenes-mediated hairy root transformation (HRT) assays are key components of the plant biotechnology and functional genomics toolkit. In this report, both in vitro and ex vitro HRT were optimized in soybean using the RUBY reporter. Different parameters including A. rhizogenes strain, optical density of the bacterial cell culture (OD₆₀₀), co-cultivation media, soybean genotype, explant age, and acetosyringone addition and concentration were evaluated. Overall, the in vitro assay was more efficient than the ex vitro assay in terms of the percentage of induction of hairy roots and transformed roots (expressing RUBY). Nonetheless, the ex vitro technique was deemed faster and a less complicated approach. The highest transformation of RUBY was observed on 7-d-old cotyledons of cv. Bert inoculated for 30 minutes with the R1000 resuspended in ¼ B5 medium to OD₆₀₀ (0.3) and 150 µM of acetosyringone. The parameters of this assay also led to the highest percentage of RUBY through two-step ex vitro hairy root transformation. Finally, using machine learning-based modeling, optimal protocols for both assays were further defined. This study establishes efficient and reliable hairy root transformation protocols applicable for functional studies in soybean.

Adaptation of high-efficiency CRISPR/Cas9-based multiplex genome editing system in white lupin by using endogenous promoters

White lupin (Lupinus albus L.) is an important crop with high phosphorus (P) use efficiency; however, technologies for its functional genomic and molecular analyses are limited. Cluster regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) (CRISPR/Cas9) system has been applied to gene editing and function genomics in many crops, but its application in white lupin has not been well documented. Here, we adapted the CRISPR/Cas9-based multiplex genome editing system by using the native U3/U6 and ubiquitin (UBQ) promoters to drive sgRNAs and Cas9. Two target sites (T1 and T2) within the Lalb_Chr05g0223881 gene, encoding a putative trehalase, were selected to validate its efficacy in white lupin based on the Agrobacterium rhizogenes-mediated transformation. We found that the T0 hairy roots were efficiently mutated at T1 and T2 with a frequency of 6.25%-35% and 50%-92.31%, respectively. The mutation types include nucleotide insertion, deletion, substitution, and complicated variant. Simultaneous mutations of the two targets were also observed with a range of 6.25%-35%. The combination of LaU6.6 promoter for sgRNA and LaUBQ12 promoter for Cas9 generated the highest frequency of homozygous/biallelic mutations at 38.46%. In addition, the target-sgRNA sequence also contributes to the editing efficiency of the CRISPR/Cas9 system in white lupin. In conclusion, our results expand the toolbox of the CRISPR/Cas9 system and benefit the basic research in white lupin.

Natural variation of GmRj2/Rfg1 determines symbiont differentiation in soybean

Symbiotic nitrogen fixation (SNF) provides much of the N utilized by leguminous plants throughout growth and development. Legumes may simultaneously establish symbiosis with different taxa of microbial symbionts. Yet, the mechanisms used to steer associations toward symbionts that are most propitious across variations in soil types remain mysterious. Here, we demonstrate that GmRj2/Rfg1 is responsible for regulating symbiosis with multiple taxa of soybean symbionts. In our experiments, the GmRj2/Rfg1SC haplotype favored association with Bradyrhizobia, which is mostly distributed in acid soils, whereas the GmRj2/Rfg1HH haplotype and knockout mutants of GmRj2/Rfg1SC associated equally with Bradyrhizobia and Sinorhizobium. Association between GmRj2/Rfg1 and NopP, furthermore, appeared to be involved in symbiont selection. Furthermore, geographic distribution analysis of 1,821 soybean accessions showed that GmRj2/Rfg1SC haplotypes were enriched in acidic soils where Bradyrhizobia were the dominant symbionts, whereas GmRj2/Rfg1HH haplotypes were most prevalent in alkaline soils dominated by Sinorhizobium, and neutral soils harbored no apparent predilections toward either haplotype. Taken together, our results suggest that GmRj2/Rfg1 regulates symbiosis with different symbionts and is a strong determinant of soybean adaptability across soil regions. As a consequence, the manipulation of the GmRj2/Rfg1 genotype or application of suitable symbionts according to the haplotype at the GmRj2/Rfg1 locus might be suitable strategies to explore for increasing soybean yield through the management of SNF.

Establishment of an efficient transformation system and its application in regulatory mechanism analysis of biological macromolecules in tea plants

Tea (Camellia sinensis), one of the most important beverage crops originated from China and is now cultivated worldwide, provides numerous secondary metabolites that account for its health benefits and rich flavor. However, the lack of an efficient and reliable genetic transformation system has seriously hindered the gene function investigation and precise breeding of C. sinensis. In this study, we established a highly efficient, labor-saving, and cost-effective Agrobacterium rhizogenes-mediated hairy roots genetic transformation system for C. sinensis, which can be used for gene overexpression and genome editing. The established transformation system was simple to operate, bypassing tissue culture and antibiotic screening, and only took two months to complete. We used this system to conduct function analysis of transcription factor CsMYB73 and found that CsMYB73 negatively regulates L-theanine synthesis in tea plant. Additionally, callus formation was successfully induced using transgenic roots, and the transgenic callus exhibited normal chlorophyll production, enabling the study of the corresponding biological functions. Furthermore, this genetic transformation system was effective for multiple C. sinensis varieties and other woody plant species. By overcoming technical obstacles such as low efficiency, long experimental periods, and high costs, this genetic transformation will be a valuable tool for routine gene investigation and precise breeding in tea plants.

Genetic Dissection of Extreme Seed-Flooding Tolerance in a Wild Soybean PI342618B by Linkage Mapping and Candidate Gene Analysis

Seed-flooding stress is one of major abiotic constraints that adversely affects soybean production worldwide. Identifying tolerant germplasms and revealing the genetic basis of seed-flooding tolerance are imperative goals for soybean breeding. In the present study, high-density linkage maps of two inter-specific recombinant inbred line (RIL) populations, named NJIRNP and NJIR4P, were utilized to identify major quantitative trait loci (QTLs) for seed-flooding tolerance using three parameters viz., germination rate (GR), normal seedling rate (NSR), and electrical conductivity (EC). A total of 25 and 18 QTLs were detected by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 12 common QTLs were identified through both methods. All favorable alleles for the tolerance are notably from the wild soybean parent. Moreover, four digenic epistatic QTL pairs were identified, and three of them showed no main effects. In addition, the pigmented soybean genotypes exhibited high seed-flooding tolerance compared with yellow seed coat genotypes in both populations. Moreover, out of five identified QTLs, one major region containing multiple QTLs associated with all three traits was identified on Chromosome 8, and most of the QTLs within this hotspot were major loci (R² > 10) and detectable in both populations and multiple environments. Based on the gene expression and functional annotation information, 10 candidate genes from QTL "hotspot 8-2" were screened for further analysis. Furthermore, the results of qRT-PCR and sequence analysis revealed that only one gene, GmDREB2 (Glyma.08G137600), was significantly induced under flooding stress and displayed a TTC tribasic insertion mutation of the nucleotide sequence in the tolerant wild parent (PI342618B). GmDREB2 encodes an ERF transcription factor, and the subcellular localization analysis using green fluorescent protein (GFP) revealed that GmDREB2 protein was localized in the nucleus and plasma membrane. Furthermore, overexpression of GmDREB2 significantly promoted the growth of soybean hairy roots, which might indicate its critical role in seed-flooding stress. Thus, GmDREB2 was considered as the most possible candidate gene for seed-flooding tolerance.

Genome-wide identification and analysis of glyceraldehyde-3-phosphate dehydrogenase family reveals the role of GmGAPDH14 to improve salt tolerance in soybean (Glycine max L.)

Introduction

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is an essential key enzyme in the glycolytic pathway and plays an important role in stress responses. Although GAPDH family genes have been found in different plant species, the determination of their gene family analysis and their functional roles in soybean are still unknown.

Methods

In this study, gene sequence and expression data were obtained using online tools, and systematic evolution, expression profile analysis, and qRT-PCR analysis were conducted.

Results and Discussion

Here a total of 16 GmGAPDH genes were identified on nine chromosomes, which were classified into three clusters. Additionally, all GmGAPDH genes harbor two highly conserved domains, including Gp_dh_N (PF00044) and Gp_dh_C (PF02800). The qRTPCR analysis also showed that most GmGAPDH genes significantly responded to multiple abiotic stresses, including NaHCO3, polyethylene glycol, cold, and salt. Among them, GmGAPDH14 was extraordinarily induced by salt stress. The GmGAPDH14 gene was cloned and overexpressed through soybean hair roots. The overexpressed transgenic soybean plants of the GmGAPDH14 gene have also shown better growth than that of control plants. Moreover, the overexpressed transgenic plants of GmGAPDH14 gene had higher activities of superoxide dismutase but lower malonaldehyde (MDA) content than those of control plants under salt stress. Meanwhile, a total of four haplotypes were found for the GmGAPDH14 gene, and haplotypes 2, 3, and 4 were beneficial for the tolerance of soybean to salt stress. These results suggest that the GmGAPDH14 gene might be involved in the process of soybean tolerance to salt stress. The results of this study will be valuable in understanding the role of GAPDH genes in the abiotic stress response of soybean.

Critical role for uricase and xanthine dehydrogenase in soybean nitrogen fixation and nodule development

De novo purine biosynthesis is required for the incorporation of fixed nitrogen in ureide exporting nodules, as formed on soybean [Glycine max (L.) Merr.] roots. However, in many cases, the enzymes involved in this pathway have been deduced strictly from genome annotations with little direct genetic evidence, such as mutant studies, to confirm their biochemical function or importance to nodule development. While efforts to develop large mutant collections of soybean are underway, research on this plant is still hampered by the inability to obtain mutations in any specific gene of interest. Using a forward genetic approach, as well as CRISPR/Cas9 gene editing via Agrobacterium rhizogenes-mediated hairy root transformation, we identified and characterized the role of GmUOX (Uricase) and GmXDH (Xanthine Dehydrogenase) in nitrogen fixation and nodule development in soybean. The gmuox knockout soybean mutants displayed nitrogen deficiency chlorosis and early nodule senescence, as exemplified by the reduced nitrogenase (acetylene reduction) activity in nodules, the internal greenish-white internal appearance of nodules, and diminished leghemoglobin production. In addition, gmuox1 nodules showed collapsed infected cells with degraded cytoplasm, aggregated bacteroids with no discernable symbiosome membranes, and increased formation of poly-β-hydroxybutyrate granules. Similarly, knockout gmxdh mutant nodules, generated with the CRISPR/Cas9 system, also exhibited early nodule senescence. These genetic studies confirm the critical role of the de novo purine metabolisms pathway not only in the incorporation of fixed nitrogen but also in the successful development of a functional, nitrogen-fixing nodule. Furthermore, these studies demonstrate the great utility of the CRISPR/Cas9 system for studying root-associated gene traits when coupled with hairy root transformation.

Genome-wide identification of genes encoding SWI/SNF components in soybean and the functional characterization of GmLFR1 in drought-stressed plants

ATP-dependent SWI/SNF chromatin remodeling complexes (CRCs) are evolutionarily conserved multi-component machines that regulate transcription, replication, and genome stability in eukaryotes. SWI/SNF components play pivotal roles in development and various stress responses in plants. However, the compositions and biological functions of SWI/SNF complex subunits remain poorly understood in soybean. In this study, we used bioinformatics to identify 39 genes encoding SWI/SNF subunit distributed on the 19 chromosomes of soybean. The promoter regions of the genes were enriched with several cis-regulatory elements that are responsive to various hormones and stresses. Digital expression profiling and qRT-PCR revealed that most of the SWI/SNF subunit genes were expressed in multiple tissues of soybean and were sensitive to drought stress. Phenotypical, physiological, and molecular genetic analyses revealed that GmLFR1 (Leaf and Flower-Related1) plays a negative role in drought tolerance in soybean and Arabidopsis thaliana. Together, our findings characterize putative components of soybean SWI/SNF complex and indicate possible roles for GmLFR1 in plants under drought stress. This study offers a foundation for comprehensive analyses of soybean SWI/SNF subunit and provides mechanistic insight into the epigenetic regulation of drought tolerance in soybean.