CRISPR/Cas9 boosts wheat yield by reducing brassinosteroid signaling.

A modern green revolution is needed to ensure global food security. Recently, Song et al. reported a new strategy to create high-yielding, semi-dwarf wheat varieties with improved nitrogen-use efficiency by inhibiting brassinosteroid (BR) signaling through clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein9 (Cas9)-mediated knockout of the ZnF-B gene encoding a zinc-finger RING-type E3 ligase.

Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Toward clinical exomes in diagnostics and management of male infertility.

Infertility, affecting ∼10% of men, is predominantly caused by primary spermatogenic failure (SPGF). We screened likely pathogenic and pathogenic (LP/P) variants in 638 candidate genes for male infertility in 521 individuals presenting idiopathic SPGF and 323 normozoospermic men in the ESTAND cohort. Molecular diagnosis was reached for 64 men with SPGF (12%), with findings in 39 genes (6%). The yield did not differ significantly between the subgroups with azoospermia (20/185, 11%), oligozoospermia (18/181, 10%), and primary cryptorchidism with SPGF (26/155, 17%). Notably, 19 of 64 LP/P variants (30%) identified in 28 subjects represented recurrent findings in this study and/or with other male infertility cohorts. NR5A1 was the most frequently affected gene, with seven LP/P variants in six SPGF-affected men and two normozoospermic men. The link to SPGF was validated for recently proposed candidate genes ACTRT1, ASZ1, GLUD2, GREB1L, LEO1, RBM5, ROS1, and TGIF2LY. Heterozygous truncating variants in BNC1, reported in female infertility, emerged as plausible causes of severe oligozoospermia. Data suggested that several infertile men may present congenital conditions with less pronounced or pleiotropic phenotypes affecting the development and function of the reproductive system. Genes regulating the hypothalamic-pituitary-gonadal axis were affected in >30% of subjects with LP/P variants. Six individuals had more than one LP/P variant, including five with two findings from the gene panel. A 4-fold increased prevalence of cancer was observed in men with genetic infertility compared to the general male population (8% vs. 2%; p = 4.4 × 10-3). Expanding genetic testing in andrology will contribute to the multidisciplinary management of SPGF.

Positive effects of public breeding on US rice yields under future climate scenarios.

In this study, we model and predict rice yields by integrating molecular marker variation, varietal productivity, and climate, focusing on the Southern U.S. rice-growing region. This region spans the states of Arkansas, Louisiana, Texas, Mississippi, and Missouri and accounts for 85% of total U.S. rice production. By digitizing and combining four decades of county-level variety acreage data (1970 to 2015) with varietal information from genotyping-by-sequencing data, we estimate annual historical county-level allele frequencies. These allele frequencies are used together with county-level weather and yield data to develop ten machine learning models for yield prediction. A two-layer meta-learner ensemble model that combines all ten methods is externally evaluated against observations from historical Uniform Regional Rice Nursery trials (1980 to 2018) conducted in the same states. Finally, the ensemble model is used with forecasted weather from the Coupled Model Intercomparison Project across the 110 rice-growing counties to predict production in the coming decades for Composite Variety Groups assembled based on year of release, breeding program, and several breeding trends. Results indicate positive effects over time of public breeding on rice resilience to future climates, and potential reasons are discussed.

Genetic improvement of phosphate-limited photosynthesis for high yield in rice.

Low phosphate (Pi) availability decreases photosynthesis, with phosphate limitation of photosynthesis occurring particularly during grain filling of cereal crops; however, effective genetic solutions remain to be established. We previously discovered that rice phosphate transporter OsPHO1;2 controls seed (sink) development through Pi reallocation during grain filling. Here, we find that OsPHO1;2 regulates Pi homeostasis and thus photosynthesis in leaves (source). Loss-of-function of OsPHO1;2 decreased Pi levels in leaves, leading to decreased photosynthetic electron transport activity, CO2 assimilation rate, and early occurrence of phosphate-limited photosynthesis. Interestingly, ectopic expression of OsPHO1;2 greatly increased Pi availability, and thereby, increased photosynthetic rate in leaves during grain filling, contributing to increased yield. This was supported by the effect of foliar Pi application. Moreover, analysis of core rice germplasm resources revealed that higher OsPHO1;2 expression was associated with enhanced photosynthesis and yield potential compared to those with lower expression. These findings reveal that phosphate-limitation of photosynthesis can be relieved via a genetic approach, and the OsPHO1;2 gene can be employed to reinforce crop breeding strategies for achieving higher photosynthetic efficiency.

Multiscale embedded printing of engineered human tissue and organ equivalents.

Creating tissue and organ equivalents with intricate architectures and multiscale functional feature sizes is the first step toward the reconstruction of transplantable human tissues and organs. Existing embedded ink writing approaches are limited by achievable feature sizes ranging from hundreds of microns to tens of millimeters, which hinders their ability to accurately duplicate structures found in various human tissues and organs. In this study, a multiscale embedded printing (MSEP) strategy is developed, in which a stimuli-responsive yield-stress fluid is applied to facilitate the printing process. A dynamic layer height control method is developed to print the cornea with a smooth surface on the order of microns, which can effectively overcome the layered morphology in conventional extrusion-based three-dimensional bioprinting methods. Since the support bath is sensitive to temperature change, it can be easily removed after printing by tuning the ambient temperature, which facilitates the fabrication of human eyeballs with optic nerves and aortic heart valves with overhanging leaflets on the order of a few millimeters. The thermosensitivity of the support bath also enables the reconstruction of the full-scale human heart on the order of tens of centimeters by on-demand adding support bath materials during printing. The proposed MSEP demonstrates broader printable functional feature sizes ranging from microns to centimeters, providing a viable and reliable technical solution for tissue and organ printing in the future.

Dissecting pleiotropic functions of the wheat Green Revolution gene Rht-B1b in plant morphogenesis and yield formation.

The utilization of reduced plant height genes Rht-B1b and Rht-D1b, encoding homeologous DELLA proteins, led to the wheat Green Revolution (GR). However, the specific functions of GR genes in yield determination and the underlying regulatory mechanisms remained unknown. Here, we validated that Rht-B1b, as a representative of GR genes, affects plant architecture and yield component traits. Upregulation of Rht-B1b reduced plant height, leaf size and grain weight, but increased tiller number, tiller angle, spike number per unit area, and grain number per spike. Dynamic investigations showed that Rht-B1b increased spike number by improving tillering initiation rather than outgrowth, and enhanced grain number by promoting floret fertility. Rht-B1b reduced plant height by reducing cell size in the internodes, and reduced grain size or weight by decreasing cell number in the pericarp. Transcriptome analyses uncovered that Rht-B1b regulates many homologs of previously reported key genes for given traits and several putative integrators for different traits. These findings specify the pleiotropic functions of Rht-B1b in improving yield and provide new insights into the regulatory mechanisms underlying plant morphogenesis and yield formation.

A retrotransposon insertion in the Mao1 promoter results in erect pubescence and higher yield in soybean.

Adaptive changes in crops contribute to the diversity of agronomic traits, which directly or indirectly affect yield. The change of pubescence form from appressed to erect is a notable feature during soybean domestication. However, the biological significance and regulatory mechanism underlying this transformation remain largely unknown. Here, we identified a major-effect locus, PUBESCENCE FORM 1 (PF1), the upstream region of Mao1, that regulates pubescence form in soybean. The insertion of a Ty3/Gypsy retrotransposon in PF1 can recruit the transcription factor GAGA-binding protein to a GA-rich region, which up-regulates Mao1 expression, underpinning soybean pubescence evolution. Interestingly, the proportion of improved cultivars with erect pubescence increases gradually with increasing latitude, and erect-pubescence cultivars have a higher yield possibly through a higher photosynthetic rate and photosynthetic stability. These findings open an avenue for molecular breeding through either natural introgression or genome editing toward yield improvement and productivity.

Similar photosynthetic but different yield responses of C3 and C4 crops to elevated O3.

The deleterious effects of ozone (O3) pollution on crop physiology, yield, and productivity are widely acknowledged. It has also been assumed that C4 crops with a carbon concentrating mechanism and greater water use efficiency are less sensitive to O3 pollution than C3 crops. This assumption has not been widely tested. Therefore, we compiled 46 journal articles and unpublished datasets that reported leaf photosynthetic and biochemical traits, plant biomass, and yield in five C3 crops (chickpea, rice, snap bean, soybean, and wheat) and four C4 crops (sorghum, maize, Miscanthus × giganteus, and switchgrass) grown under ambient and elevated O3 concentration ([O3]) in the field at free-air O3 concentration enrichment (O3-FACE) facilities over the past 20 y. When normalized by O3 exposure, C3 and C4 crops showed a similar response of leaf photosynthesis, but the reduction in chlorophyll content, fluorescence, and yield was greater in C3 crops compared with C4 crops. Additionally, inbred and hybrid lines of rice and maize showed different sensitivities to O3 exposure. This study quantitatively demonstrates that C4 crops respond less to elevated [O3] than C3 crops. This understanding could help maintain cropland productivity in an increasingly polluted atmosphere.

GmEID1 modulates light signaling through the Evening Complex to control flowering time and yield in soybean.

Soybean (Glycine max) morphogenesis and flowering time are accurately regulated by photoperiod, which determine the yield potential and limit soybean cultivars to a narrow latitudinal range. The E3 and E4 genes, which encode phytochrome A photoreceptors in soybean, promote the expression of the legume-specific flowering repressor E1 to delay floral transition under long-day (LD) conditions. However, the underlying molecular mechanism remains unclear. Here, we show that the diurnal expression pattern of GmEID1 is opposite to that of E1 and targeted mutations in the GmEID1 gene delay soybean flowering regardless of daylength. GmEID1 interacts with J, a key component of circadian Evening Complex (EC), to inhibit E1 transcription. Photoactivated E3/E4 interacts with GmEID1 to inhibit GmEID1-J interaction, promoting J degradation resulting in a negative correlation between daylength and the level of J protein. Notably, targeted mutations in GmEID1 improved soybean adaptability by enhancing yield per plant up to 55.3% compared to WT in field trials performed in a broad latitudinal span of more than 24°. Together, this study reveals a unique mechanism in which E3/E4-GmEID1-EC module controls flowering time and provides an effective strategy to improve soybean adaptability and production for molecular breeding.

Eco-physiology of maize crops under combined stresses.

The yield of maize (Zea mays L.) crops depends on their ability to intercept sunlight throughout the growing cycle, transform this energy into biomass and allocate it to the kernels. Abiotic stresses affect these eco-physiological determinants, reducing crop grain yield below the potential of each environment. Here we analyse the impact of combined abiotic stresses, such as water restriction and nitrogen deficiency or water restriction and elevated temperatures. Crop yield depends on the product of kernel yield per plant and the number of plants per unit soil area, but increasing plant population density imposes a crowding stress that reduces yield per plant, even within the range that maximises crop yield per unit soil area. Therefore, we also analyse the impact of abiotic stresses under different plant densities. We show that the magnitude of the detrimental effects of two combined stresses on field-grown plants can be lower, similar or higher than the sum of the individual stresses. These patterns depend on the timing and intensity of each one of the combined stresses and on the effects of one of the stresses on the status of the resource whose limitation causes the other. The analysis of the eco-physiological determinants of crop yield is useful to guide and prioritise the rapidly progressing studies aimed at understanding the molecular mechanisms underlying plant responses to combined stresses.

Green Revolution dwarfing Rht genes negatively affected wheat floral traits related to cross-pollination efficiency.

Hybrid breeding is a promising strategy to quickly improve wheat yield and stability. Due to the usefulness of the Rht 'Green Revolution' dwarfing alleles, it is important to gain a better understanding of their impact on traits related to hybrid development. Traits associated with cross-pollination efficiency were studied using Near Isogenic Lines carrying the different sets of alleles in Rht genes: Rht1 (semi-dwarf), Rht2 (semi-dwarf), Rht1 + 2 (dwarf), Rht3 (extreme dwarf), Rht2 + 3 (extreme dwarf), and rht (tall) during four growing seasons. Results showed that the extreme dwarfing alleles Rht2 + 3, Rht3, and Rht1 + 2 presented the greatest effects in all the traits analyzed. Plant height showed reductions up to 64% (Rht2 + 3) compared to rht. Decreases up to 20.2% in anther length and 33% in filament length (Rht2 + 3) were observed. Anthers extrusion decreased from 40% (rht) to 20% (Rht1 and Rht2), 11% (Rht3), 8.3% (Rht1 + 2), and 6.5% (Rht2 + 3). Positive correlations were detected between plant height and anther extrusion, anther, and anther filament lengths, suggesting the negative effect of dwarfing alleles. Moreover, the magnitude of these negative impacts depends on the combination of the alleles: Rht2 + 3 > Rht3/Rht1 + 2 > Rht2/Rht1 > rht (tall). Reductions were consistent across genotypes and environments with interactions due to magnitude effects. Our results indicate that Rht alleles are involved in multiple traits of interest for hybrid wheat production and the need to select alternative sources for reduced height/lodging resistance for hybrid breeding programs.

The impact of multifactorial stress combination on reproductive tissues and grain yield of a crop plant.

Global warming, climate change, and industrial pollution are altering our environment subjecting plants, microbiomes, and ecosystems to an increasing number and complexity of abiotic stress conditions, concurrently or sequentially. These conditions, termed, "multifactorial stress combination" (MFSC), can cause a significant decline in plant growth and survival. However, the impacts of MFSC on reproductive tissues and yield of major crop plants are largely unknown. We subjected soybean (Glycine max) plants to a MFSC of up to five different stresses (water deficit, salinity, low phosphate, acidity, and cadmium), in an increasing level of complexity, and conducted integrative transcriptomic-phenotypic analysis of their reproductive and vegetative tissues. We reveal that MFSC has a negative cumulative effect on soybean yield, that each set of MFSC condition elicits a unique transcriptomic response (that is different between flowers and leaves), and that selected genes expressed in leaves or flowers of soybean are linked to the effects of MFSC on different vegetative, physiological, and/or reproductive parameters. Our study identified networks and pathways associated with reactive oxygen species, ascorbic acid and aldarate, and iron/copper signaling/metabolism as promising targets for future biotechnological efforts to augment the resilience of reproductive tissues of major crop plants to MFSC. In addition, we provide unique phenotypic and transcriptomic datasets for dissecting the mechanistic effects of MFSC on the vegetative, physiological, and reproductive processes of a crop plant.

Yield reduction caused by elevated temperatures and high nitrogen fertilization is mitigated by SP6A overexpression in potato (Solanum tuberosum L.).

Potatoes (Solanum tuberosum L.) are a fundamental staple for millions of people worldwide. They provide essential amino acids, vitamins, and starch - a vital component of the human diet, providing energy and serving as a source of fiber. Unfortunately, global warming is posing a severe threat to this crop, leading to significant yield losses, and thereby endangering global food security. Industrial agriculture traditionally relies on excessive nitrogen (N) fertilization to boost yields. However, it remains uncertain whether this is effective in combating heat-related yield losses of potato. Therefore, our study aimed to investigate the combinatory effects of heat stress and N fertilization on potato tuber formation. We demonstrate that N levels and heat significantly impact tuber development. The combination of high N and heat delays tuberization, while N deficiency initiates early tuberization, likely through starvation-induced signals, independent of SELF-PRUNING 6A (SP6A), a critical regulator of tuberization. We also found that high N levels in combination with heat reduce tuber yield rather than improve it. However, our study revealed that SP6A overexpression can promote tuberization under these inhibiting conditions. By utilizing the excess of N for accumulating tuber biomass, SP6A overexpressing plants exhibit a shift in biomass distribution towards the tubers. This results in an increased yield compared to wild-type plants. Our results highlight the role of SP6A overexpression as a viable strategy for ensuring stable potato yields in the face of global warming. As such, our findings provide insights into the complex factors impacting potato crop productivity.

Yield-related quantitative trait loci identification and lint percentage hereditary dissection under salt stress in upland cotton.

Salinity is frequently mentioned as a major constraint in worldwide agricultural production. Lint percentage (LP) is a crucial yield-component in cotton lint production. While the genetic factors affect cotton yield in saline soils are still unclear. Here, we employed a recombinant inbred line population in upland cotton (Gossypium hirsutum L.) and investigated the effects of salt stress on five yield and yield component traits, including seed cotton yield per plant, lint yield per plant, boll number per plant, boll weight, and LP. Between three datasets of salt stress (E1), normal growth (E2), and the difference values dataset of salt stress and normal conditions (D-value), 87, 82, and 55 quantitative trait loci (QTL) were detectable, respectively. In total, five QTL (qLY-Chr6-2, qBNP-Chr4-1, qBNP-Chr12-1, qBNP-Chr15-5, qLP-Chr19-2) detected in both in E1 and D-value were salt related QTL, and three stable QTL (qLP-Chr5-3, qLP-Chr13-1, qBW-Chr5-5) were detected both in E1 and E2 across 3 years. Silencing of nine genes within a stable QTL (qLP-Chr5-3) highly expressed in fiber developmental stages increased LP and decreased fiber length (FL), indicating that multiple minor-effect genes clustered on Chromosome 5 regulate LP and FL. Additionally, the difference in LP caused by Gh_A05G3226 is mainly in transcription level rather than in the sequence difference. Moreover, silencing of salt related gene (GhDAAT) within qBNP-Chr4-1 decreased salt tolerance in cotton. Our findings shed light on the regulatory mechanisms underlining cotton salt tolerance and fiber initiation.

IbNF-YA1 is a key factor in the storage root development of sweet potato.

As a major worldwide root crop, the mechanism underlying storage root yield formation has always been a hot topic in sweet potato [Ipomoea batatas (L.) Lam.]. Previously, we conducted the transcriptome database of differentially expressed genes between the cultivated sweet potato cultivar "Xushu18," its diploid wild relative Ipomoea triloba without storage root, and their interspecific somatic hybrid XT1 with medium-sized storage root. We selected one of these candidate genes, IbNF-YA1, for subsequent analysis. IbNF-YA1 encodes a nuclear transcription factor Y subunit alpha (NF-YA) gene, which is significantly induced by the natural auxin indole-3-acetic acid (IAA). The storage root yield of the IbNF-YA1 overexpression (OE) plant decreased by 29.15-40.22% compared with the wild type, while that of the RNAi plant increased by 10.16-21.58%. Additionally, IAA content increased significantly in OE plants. Conversely, the content of IAA decreased significantly in RNAi plants. Furthermore, real-time quantitative reverse transcription-PCR (qRT-PCR) analysis demonstrated that the expressions of the key genes IbYUCCA2, IbYUCCA4, and IbYUCCA8 in the IAA biosynthetic pathway were significantly changed in transgenic plants. The results indicated that IbNF-YA1 could directly target IbYUCCA4 and activate IbYUCCA4 transcription. The IAA content of IbYUCCA4 OE plants increased by 71.77-98.31%. Correspondingly, the storage root yield of the IbYUCCA4 OE plant decreased by 77.91-80.52%. These findings indicate that downregulating the IbNF-YA1 gene could improve the storage root yield in sweet potato.

Light recovery after maize harvesting promotes soybean flowering in a maize-soybean relay strip intercropping system.

Moving from sole cropping to intercropping is a transformative change in agriculture, contributing to yield. Soybeans adapt to light conditions in intercropping by adjusting the onset of reproduction and the inflorescence architecture to optimize reproductive success. Maize-soybean strip intercropping (MS), maize-soybean relay strip intercropping (IS), and sole soybean (SS) systems are typical soybean planting systems with significant differences in light environments during growth periods. To elucidate the effect of changes in the light environment on soybean flowering processes and provide a theoretical basis for selecting suitable varieties in various planting systems to improve yields, field experiments combining planting systems (IS, MS, and SS) and soybean varieties (GQ8, GX7, ND25, and NN996) were conducted in 2021 and 2022. Results showed that growth recovery in the IS resulted in a balance in the expression of TERMINAL FLOWER 1 (TFL1) and FLOWERING LOCUS T (FT) in the meristematic tissues of soybeans, which promoted the formation of new branches or flowers. IS prolonged the flowering time (2-7 days) and increased the number of forming flowers compared with SS (93.0 and 169%) and MS (67.3 and 103.3%) at the later soybean flowering stage. The higher carbon and nitrogen content in the middle and bottom canopies of soybean contributed to decreased flower abscission by 26.7 and 30.2%, respectively, compared with SS. Canopy light environment recovery promoted branch and flower formation and transformation of flowers into pods with lower flower-pod abscission, which contributed to elevating soybean yields in late-maturing and multibranching varieties (ND25) in IS.

Deep learning-based association analysis of root image data and cucumber yield.

The root system is important for the absorption of water and nutrients by plants. Cultivating and selecting a root system architecture (RSA) with good adaptability and ultrahigh productivity have become the primary goals of agricultural improvement. Exploring the correlation between the RSA and crop yield is important for cultivating crop varieties with high-stress resistance and productivity. In this study, 277 cucumber varieties were collected for root system image analysis and yield using germination plates and greenhouse cultivation. Deep learning tools were used to train ResNet50 and U-Net models for image classification and segmentation of seedlings and to perform quality inspection and productivity prediction of cucumber seedling root system images. The results showed that U-Net can automatically extract cucumber root systems with high quality (F1_score ≥ 0.95), and the trained ResNet50 can predict cucumber yield grade through seedling root system image, with the highest F1_score reaching 0.86 using 10-day-old seedlings. The root angle had the strongest correlation with yield, and the shallow- and steep-angle frequencies had significant positive and negative correlations with yield, respectively. RSA and nutrient absorption jointly affected the production capacity of cucumber plants. The germination plate planting method and automated root system segmentation model used in this study are convenient for high-throughput phenotypic (HTP) research on root systems. Moreover, using seedling root system images to predict yield grade provides a new method for rapidly breeding high-yield RSA in crops such as cucumbers.

Trehalose-6-phosphate synthase 8 increases photosynthesis and seed yield in Brassica napus.

Trehalose-6-phosphate (T6P) functions as a vital proxy for assessing carbohydrate status in plants. While class II T6P synthases (TPS) do not exhibit TPS activity, they are believed to play pivotal regulatory roles in trehalose metabolism. However, their precise functions in carbon metabolism and crop yield have remained largely unknown. Here, BnaC02.TPS8, a class II TPS gene, is shown to be specifically expressed in mature leaves and the developing pod walls of Brassica napus. Overexpression of BnaC02.TPS8 increased photosynthesis and the accumulation of sugars, starch, and biomass compared to wild type. Metabolomic analysis of BnaC02.TPS8 overexpressing lines and CRISPR/Cas9 mutants indicated that BnaC02.TPS8 enhanced the partitioning of photoassimilate into starch and sucrose, as opposed to glycolytic intermediates and organic acids, which might be associated with TPS activity. Furthermore, the overexpression of BnaC02.TPS8 not only increased seed yield but also enhanced seed oil accumulation and improved the oil fatty acid composition in B. napus under both high nitrogen (N) and low N conditions in the field. These results highlight the role of class II TPS in impacting photosynthesis and seed yield of B. napus, and BnaC02.TPS8 emerges as a promising target for improving B. napus seed yield.

The vacuolar sulfate transporter PsSULTR4 is a key determinant of seed yield and protein composition in pea.

Pea is a grain legume crop with a high potential to accelerate the food transition due to its high seed protein content and relatively well-balanced amino acid composition. The critical role of external sulfur (S) supply in determining seed yield and seed quality in pea makes it essential to understand the impact of whole plant S management on the trade-off between these two traits. Here, we investigated the physiological relevance of vacuolar sulfate remobilization by targeting PsSULTR4, the only pea sulfate transporter showing substantial similarity to the vacuolar sulfate exporter AtSULTR4;1. Five mutations in PsSULTR4 were identified by TILLING (Targeting Induced Local Lesions IN Genomes), two of which, a loss of function (W78*) and a missense (E568K), significantly decreased seed yield under S deprivation. We demonstrate that PsSULTR4 triggers S distribution from source tissues, especially lower leaves, to reproductive organs to maintain seed yield under S deficiency. Under sufficient S supply, sultr4 seeds display lower levels of the S-rich storage protein PA1 at maturity. They also overaccumulate sulfate in the endosperm at the onset of seed filling. These findings uncover a role of PsSULTR4 in the remobilization of vacuolar sulfate during embryo development, allowing the efficient synthesis of S-rich proteins. Our study uncovers that PsSULTR4 functions (i) in source tissues to remobilize stored vacuolar sulfate for seed production under low S availability and (ii) in developing seeds well supplied with S to fine-tune sulfate remobilization from the endosperm as a critical control point for storage activities in the embryo.