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

Genome-Wide Association and RNA-Seq Analyses Reveal a Potential Candidate Gene Related to Oil Content in Soybean Seeds

Soybean is a crucial crop globally, serving as a significant source of unsaturated fatty acids and protein in the human diet. However, further enhancements are required for the related genes that regulate soybean oil synthesis. In this study, 155 soybean germplasms were cultivated under three different environmental conditions, followed by phenotypic identification and genome-wide association analysis using simplified sequencing data. Genome-wide association analysis was performed using SLAF-seq data. A total of 36 QTLs were significantly associated with oil content (-log10(p) > 3). Out of the 36 QTLs associated with oil content, 27 exhibited genetic overlap with previously reported QTLs related to oil traits. Further transcriptome sequencing was performed on extreme high-low oil soybean varieties. Combined with transcriptome expression data, 22 candidate genes were identified (|log2FC| ≥ 3). Further haplotype analysis of the potential candidate genes showed that three potential candidate genes had excellent haplotypes, including Glyma.03G186200, Glyma.09G099500, and Glyma.18G248900. The identified loci harboring beneficial alleles and candidate genes likely contribute significantly to the molecular network's underlying marker-assisted selection (MAS) and oil content.

Genome-Wide Identification of the Whirly Gene Family and Its Potential Function in Low Phosphate Stress in Soybean (Glycine max)

The Whirly (WHY) gene family, functioning as transcription factors, plays an essential role in the regulation of plant metabolic responses, which has been demonstrated across multiple species. However, the WHY gene family and its functions in soybean remains unclear. In this paper, we conducted genome-wide screening and identification to characterize the WHY gene family. Seven WHY members were identified and randomly distributed across six chromosomes. The phylogenetic evolutionary tree of WHY genes in soybean and other species was divided into five clades. An in-depth analysis revealed that segmental duplications significantly contributed to the expansion of GmWHYs, and the GmWHY gene members may have experienced evolutionary pressure for purifying selection in soybeans. The analysis of promoter Cis-elements in GmWHYs suggested their potential significance in addressing diverse stress conditions. The expression patterns of GmWHYs exhibited tissue-specific variations throughout the different stages of soybean development. Additionally, six GmWHY genes exhibited different responses to low phosphate stress. These findings will provide a theoretical basis and valuable reference for the future exploration of WHY gene function.

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.

The Soybean Expression Atlas v2: A comprehensive database of over 5000 RNA-seq samples

Soybean is a crucial crop worldwide, used as a source of food, feed, and industrial products due to its high protein and oil content. Previously, the rapid accumulation of soybean RNA-seq data in public databases and the computational challenges of processing raw RNA-seq data motivated us to develop the Soybean Expression Atlas, a gene expression database of over a thousand RNA-seq samples. Over the past few years, our database has allowed researchers to explore the expression profiles of important gene families, discover genes associated with agronomic traits, and understand the transcriptional dynamics of cellular processes. Here, we present the Soybean Expression Atlas v2, an updated version of our database with a fourfold increase in the number of samples, featuring transcript- and gene-level transcript abundance matrices for 5481 publicly available RNA-seq samples. New features in our database include the availability of transcript-level abundance estimates and equivalence classes to explore differential transcript usage, abundance estimates in bias-corrected counts to increase the accuracy of differential gene expression analyses, a new web interface with improved data visualization and user experience, and a reproducible and scalable pipeline available as an R package. The Soybean Expression Atlas v2 is available at https://soyatlas.venanciogroup.uenf.br/, and it will accelerate soybean research, empowering researchers with high-quality and easily accessible gene expression data.

Serendipita indica promotes P acquisition and growth in tea seedlings under P deficit conditions by increasing cytokinins and indoleacetic acid and phosphate transporter gene expression

The culturable endophytic fungus Serendipita indica has many beneficial effects on plants, but whether and how it affects physiological activities and phosphorus (P) acquisition of tea seedlings at low P levels is unclear. The objective of this study was to analyze the effects of inoculation with S. indica on growth, gas exchange, chlorophyll fluorescence, auxins, cytokinins, P levels, and expressions of two phosphate transporter (PT) genes in leaves of tea (Camellia sinensis L. cv. Fudingdabaicha) seedlings grown at 0.5 μM (P_0.5) and 50 μM (P₅₀) P levels. Sixteen weeks after the inoculation, S. indica colonized roots of tea seedlings, with root fungal colonization rates reaching 62.18% and 81.34% at P_0.5 and P₅₀ levels, respectively. Although plant growth behavior, leaf gas exchange, chlorophyll values, nitrogen balance index, and chlorophyll fluorescence parameters of tea seedlings were suppressed at P_0.5 versus P₅₀ levels, inoculation of S. indica mitigated the negative effects to some extent, along with more prominent promotion at P_0.5 levels. S. indica inoculation significantly increased leaf P and indoleacetic acid concentrations at P_0.5 and P₅₀ levels and leaf isopentenyladenine, dihydrozeatin, and transzeatin concentrations at P_0.5 levels, coupled with the reduction of indolebutyric acid at P₅₀ levels. Inoculation of S. indica up-regulated the relative expression of leaf CsPT1 at P_0.5 and P₅₀ levels and CsPT4 at P_0.5 levels. It is concluded that S. indica promoted P acquisition and growth in tea seedlings under P deficit conditions by increasing cytokinins and indoleacetic acid and CsPT1 and CsPT4 expression.

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.

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.

Comparative analysis of physiological variations and genetic architecture for cold stress response in soybean germplasm

Soybean (Glycine max L.) is susceptible to low temperatures. Increasing lines of evidence indicate that abiotic stress-responsive genes are involved in plant low-temperature stress response. However, the involvement of photosynthesis, antioxidants and metabolites genes in low temperature response is largely unexplored in Soybean. In the current study, a genetic panel of diverse soybean varieties was analyzed for photosynthesis, chlorophyll fluorescence and leaf injury parameters under cold stress and control conditions. This helps us to identify cold tolerant (V100) and cold sensitive (V45) varieties. The V100 variety outperformed for antioxidant enzymes activities and relative expression of photosynthesis (Glyma.08G204800.1, Glyma.12G232000.1), GmSOD (GmSOD01, GmSOD08), GmPOD (GmPOD29, GmPOD47), trehalose (GmTPS01, GmTPS13) and cold marker genes (DREB1E, DREB1D, SCOF1) than V45 under cold stress. Upon cold stress, the V100 variety showed reduced accumulation of H₂O₂ and MDA levels and subsequently showed lower leaf injury compared to V45. Together, our results uncovered new avenues for identifying cold tolerant soybean varieties from a large panel. Additionally, we identified the role of antioxidants, osmo-protectants and their posttranscriptional regulators miRNAs such as miR319, miR394, miR397, and miR398 in Soybean cold stress tolerance.

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.