The Myb73-GDPD2-GA2ox1 transcriptional regulatory module confers phosphate deficiency tolerance in soybean.

Phosphorus is indispensable in agricultural production. An increasing food supply requires more efficient use of phosphate due to limited phosphate resources. However, how crops regulate phosphate efficiency remains largely unknown. Here, we identified a major quantitative trait locus, qPE19, that controls 7 low-phosphate (LP)-related traits in soybean (Glycine max) through linkage mapping and genome-wide association studies. We identified the gene responsible for qPE19 as GLYCEROPHOSPHORYL DIESTER PHOSPHODIESTERASE2 (GmGDPD2), and haplotype 5 represents the optimal allele favoring LP tolerance. Overexpression of GmGDPD2 significantly affects hormone signaling and improves root architecture, phosphate efficiency and yield-related traits; conversely, CRISPR/Cas9-edited plants show decreases in these traits. GmMyb73 negatively regulates GmGDPD2 by directly binding to its promoter; thus, GmMyb73 negatively regulates LP tolerance. GmGDPD2 physically interacts with GA 2-oxidase 1 (GmGA2ox1) in the plasma membrane, and overexpressing GmGA2ox1 enhances LP-associated traits, similar to GmGDPD2 overexpression. Analysis of double mutants for GmGDPD2 and GmGA2ox1 demonstrated that GmGDPD2 regulates LP tolerance likely by influencing auxin and gibberellin dose-associated cell division in the root. These results reveal a regulatory module that plays a major role in regulating LP tolerance in soybeans and is expected to be utilized to develop phosphate-efficient varieties to enhance soybean production, particularly in phosphate-deficient soils.

Soybean type-B response regulator GmRR1 mediates phosphorus uptake and yield by modifying root architecture.

Phosphorus (P) plays a pivotal role in plant growth and development. Low P stress can greatly hamper plant growth. Here, we identified a QTL (named QPH-9-1), which is associated with P efficiency across multiple environments through linkage analysis and genome-wide association study. Furthermore, we successfully cloned the underlying soybean (Glycine max) gene GmRR1 (a soybean type-B Response Regulator 1) that encodes a type-B response regulator protein. Knockout of GmRR1 resulted in a substantial increase in plant height, biomass, P uptake efficiency, and yield-related traits due to the modification of root structure. In contrast, overexpression of GmRR1 in plants resulted in a decrease in these phenotypes. Further analysis revealed that knockout of GmRR1 substantially increased the levels of auxin and ethylene in roots, thereby promoting root hair formation and growth by promoting the formation of root hair primordium and lengthening the root apical meristem. Yeast two-hybrid, bimolecular fluorescence complementation, and dual-luciferase assays demonstrated an interaction between GmRR1 and Histidine-containing Phosphotransmitter protein 1. Expression analysis suggested that these proteins coparticipated in response to low P stress. Analysis of genomic sequences showed that GmRR1 underwent a selection during soybean domestication. Taken together, this study provides further insights into how plants respond to low P stress by modifying root architecture through phytohormone pathways.

Transcription factors GmERF1 and GmWRKY6 synergistically regulate low phosphorus tolerance in soybean.

Soybean (Glycine max) is a major grain and oil crop worldwide, but low phosphorus (LP) in soil severely limits the development of soybean production. Dissecting the regulatory mechanism of the phosphorus (P) response is crucial for improving the P use efficiency of soybean. Here, we identified a transcription factor, GmERF1 (ethylene response factor 1), that is mainly expressed in soybean root and localized in the nucleus. Its expression is induced by LP stress and differs substantially in extreme genotypes. The genomic sequences of 559 soybean accessions suggested that the allelic variation of GmERF1 has undergone artificial selection, and its haplotype is significantly related to LP tolerance. GmERF1 knockout or RNA interference resulted in significant increases in root and P uptake efficiency traits, while the overexpression of GmERF1 produced an LP-sensitive phenotype and affected the expression of 6 LP stress-related genes. In addition, GmERF1 directly interacted with GmWRKY6 to inhibit transcription of GmPT5 (phosphate transporter 5), GmPT7, and GmPT8, which affects plant P uptake and use efficiency under LP stress. Taken together, our results show that GmERF1 can affect root development by regulating hormone levels, thus promoting P absorption in soybean, and provide a better understanding of the role of GmERF1 in soybean P signal transduction. The favorable haplotypes from wild soybean will be conducive to the molecular breeding of high P use efficiency in soybean.

Genome-wide identification, evolution, and expression analysis of carbonic anhydrases genes in soybean (Glycine max).

Carbonic anhydrases (CAs), as zinc metalloenzymes, are ubiquitous in nature and play essential roles in diverse biological processes. Although CAs have been broadly explored and studied, comprehensive characteristics of CA gene family members in the soybean (Glycine max) are still lacking. A total of 35 CA genes (GmCAs) were identified; they distributed on sixteen chromosomes of the soybean genome and can be divided into three subfamilies (α-type, ß-type, and γ-type). Bioinformatics analysis showed that the specific GmCA gene subfamily or clade exhibited similar characteristics and that segmental duplications took the major role in generating new GmCAs. Furthermore, the synteny and evolutionary constraints analyses of CAs among soybean and distinct species provided more detailed evidence for GmCA gene family evolution. Cis-element analysis of promoter indicated that GmCAs may be responsive to abiotic stress and regulate photosynthesis. Moreover, the expression patterns of GmCAs varied in different tissues at diverse developmental stages in soybean. Additionally, we found that eight representative GmCAs may be involved in the response of soybean to low phosphorus stress. The systematic investigation of the GmCA gene family in this study will provide a valuable basis for further functional research on soybean CA genes.

GmEIL4 enhances soybean (Glycine max) phosphorus efficiency by improving root system development.

Phosphorus (P) deficiency seriously affects plant growth and development and ultimately limits the quality and yield of crops. Here, a new P efficiency-related major quantitative trait locus gene, GmEIL4 (encoding an ethylene-insensitive 3-like 1 protein), was cloned at qP2, which was identified by linkage analysis and genome-wide association study across four environments. Overexpressing GmEIL4 significantly improved the P uptake efficiency by increasing the number, length and surface area of lateral roots of hairy roots in transgenic soybeans, while interfering with GmEIL4 resulted in poor root phenotypic characteristics compared with the control plants under low P conditions. Interestingly, we found that GmEIL4 interacted with EIN3-binding F box protein 1 (GmEBF1), which may regulate the root response to low P stress. We conclude that the expression of GmEIL4 was induced by low-P stress and that overexpressing GmEIL4 improved P accumulation by regulating root elongation and architecture. Analysis of allele variation of GmEIL4 in 894 soybean accessions suggested that GmEIL4 is undergoing artificial selection during soybean evolution, which will benefit soybean production. Together, this study further elucidates how plants respond to low P stress by modifying root structure and provides insight into the great potential of GmEIL4 in crop P-efficient breeding.

Identification of soybean phosphorous efficiency QTLs and genes using chlorophyll fluorescence parameters through GWAS and RNA-seq.

MAIN CONCLUSION: Soybean phosphorous efficiency QTLs were identified and candidate genes were predicted using chlorophyll fluorescence parameters through GWAS and RNA-seq. Phosphorus (P) is an essential nutrient element for crop growth and development, lack of P uptake seriously affects yield in various crops. Photosynthesis is the basis of crop production, while it is very sensitive to P deficiency. It is of great importance to study the genetic relationship between photosynthesis and P efficiency to provide genetic insight for soybean improvement. In this study, a genome-wide association study (GWAS) was performed using 292,035 SNPs and the ratios of four main chlorophyll fluorescence parameters of 219 diverse soybean accessions under P deficiency and normal P across three experiments. In total, 52 SNPs in 12 genomic regions were detected in association with the four main chlorophyll fluorescence parameters under sufficient or deficient P levels. Combined it with RNA-seq analysis, we predicted three candidate genes for the significant genomic regions. For example, the expression level of the candidate gene (Glyma.18g092900) in P deficiency tolerant accession was three times higher than that of P deficiency sensitive one under phosphorous deficiency condition. This study provides insight into genetic links between photosynthetic and phosphorous efficiency and further functional analysis will provide valuable information for understanding the underlying genetic mechanism to facilitate marker-assisted breeding in soybean.

The purple acid phosphatase GmPAP17 predominantly enhances phosphorus use efficiency in soybean.

Purple acid phosphatase (PAP) is an important plant acid phosphatase, which can secrete to the rhizosphere to decompose organophosphorus, promote phosphorus use efficiency, plant growth and development. However, little is known about the functions of intracellular PAP in plants, especially for soybean. Our previous study integrating QTL mapping and transcriptome analysis identified an promising low phosphorus (LP)-induced gene GmPAP17. Here, we determined that GmPAP17 was mainly expressed in roots and had a strong response to LP stress. Furthermore, and the relative expression in the root of LP tolerant genotypes NN94-156 was significantly greater than that of LP sensitive genotype Bogao after LP stress treatment. The overexpression of GmPAP17 significantly enhanced both acid phosphatase activity and growth performance of hairy roots under LP stress condition, it was vice versa for RNAi interference of GmPAP17, indicating that GmPAP17 plays an important role in P use efficiency. Moreover, yeast two-hybrid and bimolecular fluorescence complementation analysis showed that GmRAP2.2 was involved in the regulation network of GmPAP17. Taken together, our results suggest that GmPAP17 is a novel plant PAP that functions in the adaptation of soybean to LP stress, possibly through its involvement in P recycling in plants.

GmWRKY46, a WRKY transcription factor, negatively regulates phosphorus tolerance primarily through modifying root morphology in soybean.

Phosphorus (P) deficiency affects soybean growth and development, resulting in significant reduction of yields. However, the regulatory mechanism of P deficiency tolerance in soybean is still largely unclear. WRKY transcription factors are a family of regulators involved in a variety of abiotic stresses in plants while rarely reported in P deficiency. Here, we demonstrated that a soybean GmWRKY46 gene, belonging to group III of WRKY TF family, was involved in the regulation of P deficiency tolerance in soybean. The expression of GmWRKY46 in low P sensitive soybean varieties was significantly higher than that in tolerant soybean varieties. It was primarily expressed in roots and strongly induced by P deprivation. GmWRKY46 was localized in the nucleus. Compared with the control expressing the empty vector, overexpression of GmWRKY46 in soybean hairy roots exhibited more sensitive phenotypes to low P stress, while the RNA interfered GmWRKY46 significantly enhanced P deficiency tolerance by increasing the proliferation, elongation and P absorption efficiency of hairy roots. Expression patterns of a number of P-responsive genes (GmPht1;1, GmPht1;4, GmPTF1, GmACP1, GmPAP21 and GmExpansin-A7) were altered in both overexpression and gene silenced plants. The results provided a novel insight into how soybean responds to low P stress and new gene that may be used to improve soybean low P tolerance through gene editing approach.