Integrating QTL mapping and transcriptomics identifies candidate genes underlying QTLs associated with soybean tolerance to low-phosphorus stress

Bioengineering and management for efficient and sustainable utilization of phosphorus in crops

Phosphorus (P) is an essential macronutrient for plant growth, but low P availability in soils is also a primary constraint to crop production. To meet the increasing demands for food, P fertilizer applications have been increased, causing the accumulation of surplus P in soils, which has led to the frequency and magnitude of associated risk effects on agroecosystems. Finding solutions for efficient and sustainable crop P utilization is, therefore, an urgent priority. This review summarizes recent progress in bioengineering approaches to improving crop P efficiency and highlights that modifying root architecture in P-deficient soils and reducing P accumulation in grains in soils with P surplus could offer a way forward for improving P use efficiency.

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.

Integrated Bulk Segregant Analysis, Fine Mapping, and Transcriptome Revealed QTLs and Candidate Genes Associated with Drought Adaptation in Wild Watermelon

Drought stress has detrimental effects on crop productivity worldwide. A strong root system is crucial for maintaining water and nutrients uptake under drought stress. Wild watermelons possess resilient roots with excellent drought adaptability. However, the genetic factors controlling this trait remain uninvestigated. In this study, we conducted a bulk segregant analysis (BSA) on an F2 population consisting of two watermelon genotypes, wild and domesticated, which differ in their lateral root development under drought conditions. We identified two quantitative trait loci (qNLR_Dr. Chr01 and qNLR_Dr. Chr02) associated with the lateral root response to drought. Furthermore, we determined that a small region (0.93 Mb in qNLR_Dr. Chr01) is closely linked to drought adaptation through quantitative trait loci (QTL) validation and fine mapping. Transcriptome analysis of the parent roots under drought stress revealed unique effects on numerous genes in the sensitive genotype but not in the tolerant genotype. By integrating BSA, fine mapping, and the transcriptome, we identified six genes, namely L-Ascorbate Oxidase (AO), Cellulose Synthase-Interactive Protein 1 (CSI1), Late Embryogenesis Abundant Protein (LEA), Zinc-Finger Homeodomain Protein 2 (ZHD2), Pericycle Factor Type-A 5 (PFA5), and bZIP transcription factor 53-like (bZIP53-like), that might be involved in the drought adaptation. Our findings provide valuable QTLs and genes for marker-assisted selection in improving water-use efficiency and drought tolerance in watermelon. They also lay the groundwork for the genetic manipulation of drought-adapting genes in watermelon and other Cucurbitacea species.

Integrated Analysis of Metabolome and Transcriptome Reveals Insights for Low Phosphorus Tolerance in Wheat Seedling

Low phosphorus (LP) stress leads to a significant reduction in wheat yield, primarily in the reduction of biomass, the number of tillers and spike grains, the delay in heading and flowering, and the inhibition of starch synthesis and grouting. However, the differences in regulatory pathway responses to low phosphorus stress among different wheat genotypes are still largely unknown. In this study, metabolome and transcriptome analyses of G28 (LP-tolerant) and L143 (LP-sensitive) wheat varieties after 72 h of normal phosphorus (CK) and LP stress were performed. A total of 181 and 163 differentially accumulated metabolites (DAMs) were detected for G28CK vs. G28LP and L143CK vs. L143LP, respectively. Notably, the expression of pilocarpine (C07474) in G28CK vs. G28LP was significantly downregulated 4.77-fold, while the expression of neochlorogenic acid (C17147) in L143CK vs. L143LP was significantly upregulated 2.34-fold. A total of 4023 differentially expressed genes (DEGs) were acquired between G28 and L143, of which 1120 DEGs were considered as the core DEGs of LP tolerance of wheat after LP treatment. The integration of metabolomics and transcriptomic data further revealed that the LP tolerance of wheat was closely related to 15 metabolites and 18 key genes in the sugar and amino acid metabolism pathway. The oxidative phosphorylation pathway was enriched to four ATPases, two cytochrome c reductase genes, and fumaric acid under LP treatment. Moreover, PHT1;1, TFs (ARFA, WRKY40, MYB4, MYB85), and IAA20 genes were related to the Pi starvation stress of wheat roots. Therefore, the differences in LP tolerance of different wheat varieties were related to energy metabolism, amino acid metabolism, phytohormones, and PHT proteins, and precisely regulated by the levels of various molecular pathways to adapt to Pi starvation stress. Taken together, this study may help to reveal the complex regulatory process of wheat adaptation to Pi starvation and provide new genetic clues for further study on improving plant Pi utilization efficiency.

Transcriptomic Analysis of the Response of the Toxic Dinoflagellate Prorocentrum lima to Phosphorous Limitation

Some dinoflagellates cause harmful algal blooms, releasing toxic secondary metabolites, to the detriment of marine ecosystems and human health. Phosphorus (P) is a limiting macronutrient for dinoflagellate growth in the ocean. Previous studies have been focused on the physiological response of dinoflagellates to ambient P changes. However, the whole-genome's molecular mechanisms are poorly understood. In this study, RNA-Seq was utilized to compare the global gene expression patterns of a marine diarrheic shellfish poisoning (DSP) toxin-producing dinoflagellate, Prorocentrum lima, grown in inorganic P-replete and P-deficient conditions. A total of 148 unigenes were significantly up-regulated, and 30 unigenes were down-regulated under 1/4 P-limited conditions, while 2708 unigenes were significantly up-regulated, and 284 unigenes were down-regulated under 1/16 P-limited conditions. KEGG enrichment analysis of the differentially expressed genes shows that genes related to ribosomal proteins, glycolysis, fatty acid biosynthesis, phagosome formation, and ubiquitin-mediated proteolysis are found to be up-regulated, while most of the genes related to photosynthesis are down-regulated. Further analysis shows that genes encoding P transporters, organic P utilization, and endocytosis are significantly up-regulated in the P-limited cells, indicating a strong ability of P. lima to utilize dissolved inorganic P as well as intracellular organic P. These transcriptomic data are further corroborated by biochemical and physiological analyses, which reveals that under P deficiency, cellular contents of starch, lipid, and toxin increase, while photosynthetic efficiency declines. Our results indicate that has P. lima evolved diverse strategies to acclimatize to low P environments. The accumulation of carbon sources and DSP toxins could provide protection for P. lima to cope with adverse environmental conditions.

Joint QTL Mapping and Transcriptome Sequencing Analysis Reveal Candidate Seed-Shattering-Related Genes in Common Buckwheat

Common buckwheat (Fagopyrum esculentum M.) is an important traditional miscellaneous grain crop. However, seed-shattering is a significant problem in common buckwheat. To investigate the genetic architecture and genetic regulation of seed-shattering in common buckwheat, we constructed a genetic linkage map using the F₂ population of Gr (green-flower mutant and shattering resistance) and UD (white flower and susceptible to shattering), which included eight linkage groups with 174 loci, and detected seven QTLs of pedicel strength. RNA-seq analysis of pedicel in two parents revealed 214 differentially expressed genes DEGs that play roles in phenylpropanoid biosynthesis, vitamin B6 metabolism, and flavonoid biosynthesis. Weighted gene co-expression network analysis (WGCNA) was performed and screened out 19 core hub genes. Untargeted GC-MS analysis detected 138 different metabolites and conjoint analysis screened out 11 DEGs, which were significantly associated with differential metabolites. Furthermore, we identified 43 genes in the QTLs, of which six genes had high expression levels in the pedicel of common buckwheat. Finally, 21 candidate genes were screened out based on the above analysis and gene function. Our results provided additional knowledge for the identification and functions of causal candidate genes responsible for the variation in seed-shattering and would be an invaluable resource for the genetic dissection of common buckwheat resistance-shattering molecular breeding.

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.

Phosphorus Availability Affects the Photosynthesis and Antioxidant System of Contrasting Low-P-Tolerant Cotton Genotypes

Phosphorus (P) is an essential macronutrient, and an important component of plant metabolism. However, little is known about the effects of low P availability on P absorption, the photosynthetic electron transport chain, and the antioxidant system in cotton. This study used cotton genotypes (sensitive FJA and DLNTDH and tolerant BX014 and LuYuan343) with contrasting low-P tolerance in a hydroponic experiment under 15 µM, 50 µM, and 500 μM P concentrations. The results showed that low P availability reduced plant development and leaf area, shoot length, and dry weight in FJA and DLNADH, compared to BX014 and LuYuan343. The low P availability decreased the gas-exchange parameters such as the net photosynthetic rate, transpiration rate, and stomatal conductance, and increased the intercellular CO₂ concentration. Chlorophyll a fluorescence demonstrated that the leaves' absorption and trapped-energy flux were largely steady. In contrast, considerable gains in absorption and trapped-energy flux per reaction center resulted from decreases in the electron transport per reaction center under low-P conditions. In addition, low P availability reduced the activities of antioxidant enzymes and increased the content of malondialdehyde in the cotton genotypes, especially in FJA and DLNTDH. Moreover, low P availability reduced the activity of PEPC and generated a decline in the content of ATP and NADPH. Our research can provide a theoretical physiological basis for the growth and tolerance of cotton under low-P conditions.

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.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

ZmbZIP54 and ZmFDX5 cooperatively regulate maize seedling tolerance to lead by mediating ZmPRP1 transcription

Extensive lead (Pb) accumulation in plants exerts toxic effects on plant growth and development and enters the human food chain. Combining linkage mapping, transcriptome analysis, and association studies, we cloned the ZmbZIP54 transcription factor, which confers maize tolerance to Pb. Combined overexpression and knockdown confirmed that ZmbZIP54 mitigates Pb toxicity in maize by alleviating Pb absorption into the roots. Yeast one-hybrid and dual-luciferase assays revealed that ZmbZIP54 binds to the ZmPRP1 promoter and promotes its transcription. Yeast two-hybrid and bimolecular fluorescence complementation assays indicated that ZmFdx5 interacts with ZmbZIP54 in the nucleus. ZmFdx5 acts as a switch that controls the regulation of ZmPRP1 expression by ZmbZIP54 when maize encounters Pb stress. Furthermore, we revealed that variation in the 5'-UTR of ZmbZIP54 affects its expression level under Pb stress and contributes to the difference in Pb tolerance among maize lines. Finally, we proposed a model to summarize the role of ZmbZIP54 in Pb tolerance, which involves the cooperative effect of ZmbZIP54 and ZmFdx5 on the ZmPRP1 transcription in maize response to Pb. This study provides novel insights into the development of Pb-tolerant maize varieties and bioremediation of Pb-contaminated soils.

Genetic approaches for assessment of phosphorus use efficiency in groundnut (Arachis hypogaea L.)

Production of phosphorus efficient genotypes in groundnut can improve and also reduces environmental pollution. Identification of P-efficient groundnut genotypes is a need of the hour to sustain in P-deficient soils. The pot experiment showed significant differences between genotypes (G) and treatments (T) for all the traits and G × T interaction for majority of traits. The G × T × Y interaction effects were also significant for all the traits except leaf P% (LP%), leaf acid phosphatase (LAP) and root dry weight (RDW). In lysimeter experiment, the effect of G, T and G × T were significant for leaf dry weight (LDW), stem dry weight (SDW), total transpiration (TT) and transpiration efficiency (TE). For traits, LDW, SDW, TT, TE, ICGV 00351 and ICGS 76; for SDW, TT, ICGV 02266 are best performers under both P-sufficient and deficient conditions. Based on P-efficiency indices and surrogate traits of P-uptake, ICGV's 02266, 05155, 00308, 06040 and 06146 were considered as efficient P-responding genotypes. From GGE biplot, ICGV 06146 under P-deficient and TAG 24 under both P-sufficient and deficient conditions are portrayed as best performer. ICGV 06146 was identified as stable pod yielder and a promising genotype for P-deficient soils. The genotypes identified in this study can be used as a parent in developing mapping population to decipher the genetics and to devleop groundnut breeding lines suitable to P-deficient soils.

Improvement of heat stress tolerance in soybean (Glycine max L), by using conventional and molecular tools

The soybean is a significant legume crop, providing several vital dietary components. Extreme heat stress negatively affects soybean yield and quality, especially at the germination stage. Continuous change in climatic conditions is threatening the global food supply and food security. Therefore, it is a critical need of time to develop heat-tolerant soybean genotypes. Different molecular techniques have been developed to improve heat stress tolerance in soybean, but until now complete genetic mechanism of soybean is not fully understood. Various molecular methods, like quantitative trait loci (QTL) mapping, genetic engineering, transcription factors (TFs), transcriptome, and clustered regularly interspaced short palindromic repeats (CRISPR), are employed to incorporate heat tolerance in soybean under the extreme conditions of heat stress. These molecular techniques have significantly improved heat stress tolerance in soybean. Besides this, we can also use specific classical breeding approaches and different hormones to reduce the harmful consequences of heat waves on soybean. In future, integrated use of these molecular tools would bring significant results in developing heat tolerance in soybean. In the current review, we have presented a detailed overview of the improvement of heat tolerance in soybean and highlighted future prospective. Further studies are required to investigate different genetic factors governing the heat stress response in soybean. This information would be helpful for future studies focusing on improving heat tolerance 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.

Mechanisms Underlying Soybean Response to Phosphorus Deficiency through Integration of Omics Analysis

Low phosphorus (P) availability limits soybean growth and yield. A set of potential strategies for plant responses to P deficiency have been elucidated in the past decades, especially in model plants such as Arabidopsis thaliana and rice (Oryza sativa). Recently, substantial efforts focus on the mechanisms underlying P deficiency improvement in legume crops, especially in soybeans (Glycine max). This review summarizes recent advances in the morphological, metabolic, and molecular responses of soybean to phosphate (Pi) starvation through the combined analysis of transcriptomics, proteomics, and metabolomics. Furthermore, we highlight the functions of the key factors controlling root growth and P homeostasis, base on which, a P signaling network in soybean was subsequently presumed. This review also discusses current barriers and depicts perspectives in engineering soybean cultivars with high P efficiency.

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.

Constraints and Prospects of Improving Cowpea Productivity to Ensure Food, Nutritional Security and Environmental Sustainability

Providing safe and secure food for an increasing number of people globally is challenging. Coping with such a human population by merely applying the conventional agricultural production system has not proved to be agro-ecologically friendly; nor is it sustainable. Cowpea (Vigna unguiculata (L) Walp) is a multi-purpose legume. It consists of high-quality protein for human consumption, and it is rich in protein for livestock fodder. It enriches the soil in that it recycles nutrients through the fixation of nitrogen in association with nodulating bacteria. However, the productivity of this multi-functional, indigenous legume that is of great value to African smallholder farmers and the rural populace, and also to urban consumers and entrepreneurs, is limited. Because cowpea is of strategic importance in Africa, there is a need to improve on its productivity. Such endeavors in Africa are wrought with challenges that include drought, salinity, the excessive demand among farmers for synthetic chemicals, the repercussions of climate change, declining soil nutrients, microbial infestations, pest issues, and so forth. Nevertheless, giant strides have already been made and there have already been improvements in adopting sustainable and smart biotechnological approaches that are favorably influencing the production costs of cowpea and its availability. As such, the prospects for a leap in cowpea productivity in Africa and in the enhancement of its genetic gain are good. Potential and viable means for overcoming some of the above-mentioned production constraints would be to focus on the key cowpea producer nations in Africa and to encourage them to embrace biotechnological techniques in an integrated approach to enhance for sustainable productivity. This review highlights the spectrum of constraints that limit the cowpea yield, but looks ahead of the constraints and seeks a way forward to improve cowpea productivity in Africa. More importantly, this review investigates applications and insights concerning mechanisms of action for implementing eco-friendly biotechnological techniques, such as the deployment of bio inoculants, applying climate-smart agricultural (CSA) practices, agricultural conservation techniques, and multi-omics smart technology in the spheres of genomics, transcriptomics, proteomics, and metabolomics, for improving cowpea yields and productivity to achieve sustainable agro-ecosystems, and ensuring their stability.

Transcriptome profiling reveals the spatial-temporal dynamics of gene expression essential for soybean seed development

BACKGROUND

Seeds are the economic basis of oilseed crops, especially soybeans, the most widely cultivated oilseed crop worldwide. Seed development is accompanied by a multitude of diverse cellular processes, and revealing the underlying regulatory activities is critical for seed improvement.

RESULTS

In this study, we profiled the transcriptomes of developing seeds at 20, 25, 30, and 40 days after flowering (DAF), as these stages represent critical time points of seed development from early to full development. We identified a set of highly abundant genes and highlighted the importance of these genes in supporting nutrient accumulation and transcriptional regulation for seed development. We identified 8925 differentially expressed genes (DEGs) that exhibited temporal expression patterns over the course and expression specificities in distinct tissues, including seeds and nonseed tissues (roots, stems, and leaves). Genes specific to nonseed tissues might have tissue-associated roles, with relatively low transcript abundance in developing seeds, suggesting their spatially supportive roles in seed development. Coexpression network analysis identified several underexplored genes in soybeans that bridge tissue-specific gene modules.

CONCLUSIONS

Our study provides a global view of gene activities and biological processes critical for seed formation in soybeans and prioritizes a set of genes for further study. The results of this study help to elucidate the mechanism controlling seed development and storage reserves.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting soybean genotypes subjected to phosphate starvation

BACKGROUND

Phosphorus (P) is essential for plant growth and development, and low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Long noncoding RNAs (lncRNAs) have recently been reported to be key regulators in the responses of plants to stress conditions, but the mechanism through which LP stress mediates the biogenesis of lncRNAs in soybean remains unclear.

RESULTS

In this study, to explore the response mechanisms of lncRNAs to LP stress, we used the roots of two representative soybean genotypes that present opposite responses to P deficiency, namely, a P-sensitive genotype (Bogao) and a P-tolerant genotype (NN94156), for the construction of RNA sequencing (RNA-seq) libraries. In total, 4,166 novel lncRNAs, including 525 differentially expressed (DE) lncRNAs, were identified from the two genotypes at different P levels. GO and KEGG analyses indicated that numerous DE lncRNAs might be involved in diverse biological processes related to phosphate, such as lipid metabolic processes, catalytic activity, cell membrane formation, signal transduction, and nitrogen fixation. Moreover, lncRNA-mRNA-miRNA and lncRNA-mRNA networks were constructed, and the results identified several promising lncRNAs that might be highly valuable for further analysis of the mechanism underlying the response of soybean to LP stress.

CONCLUSIONS

These results revealed that LP stress can significantly alter the genome-wide profiles of lncRNAs, particularly those of the P-sensitive genotype Bogao. Our findings increase the understanding of and provide new insights into the function of lncRNAs in the responses of soybean to P stress.

Genomic Architecture of Phenotypic Plasticity in Response to Water Stress in Tetraploid Wheat

Phenotypic plasticity is one of the main mechanisms of adaptation to abiotic stresses via changes in critical developmental stages. Altering flowering phenology is a key evolutionary strategy of plant adaptation to abiotic stresses, to achieve the maximum possible reproduction. The current study is the first to apply the linear regression residuals as drought plasticity scores while considering the variation in flowering phenology and traits under non-stress conditions. We characterized the genomic architecture of 17 complex traits and their drought plasticity scores for quantitative trait loci (QTL) mapping, using a mapping population derived from a cross between durum wheat (Triticum turgidum ssp. durum) and wild emmer wheat (T. turgidum ssp. dicoccoides). We identified 79 QTLs affected observed traits and their plasticity scores, of which 33 reflected plasticity in response to water stress and exhibited epistatic interactions and/or pleiotropy between the observed and plasticity traits. Vrn-B3 (TaTF1) residing within an interval of a major drought-escape QTL was proposed as a candidate gene. The favorable alleles for most of the plasticity QTLs were contributed by wild emmer wheat, demonstrating its high potential for wheat improvement. Our study presents a new approach for the quantification of plant adaptation to various stresses and provides new insights into the genetic basis of wheat complex traits under water-deficit stress.

Combining QTL Mapping and Transcriptomics to Decipher the Genetic Architecture of Phenolic Compounds Metabolism in the Conifer White Spruce

Conifer forests worldwide are becoming increasingly vulnerable to the effects of climate change. Although the production of phenolic compounds (PCs) has been shown to be modulated by biotic and abiotic stresses, the genetic basis underlying the variation in their constitutive production level remains poorly documented in conifers. We used QTL mapping and RNA-Seq to explore the complex polygenic network underlying the constitutive production of PCs in a white spruce (Picea glauca) full-sib family for 2 years. QTL detection was performed for nine PCs and differentially expressed genes (DEGs) were identified between individuals with high and low PC contents for five PCs exhibiting stable QTLs across time. A total of 17 QTLs were detected for eight metabolites, including one major QTL explaining up to 91.3% of the neolignan-2 variance. The RNA-Seq analysis highlighted 50 DEGs associated with phenylpropanoid biosynthesis, several key transcription factors, and a subset of 137 genes showing opposite expression patterns in individuals with high levels of the flavonoids gallocatechin and taxifolin glucoside. A total of 19 DEGs co-localized with QTLs. Our findings represent a significant step toward resolving the genomic architecture of PC production in spruce and facilitate the functional characterization of genes and transcriptional networks responsible for differences in constitutive production of PCs in conifers.

Genome-wide identification of low phosphorus responsive microRNAs in two soybean genotypes by high-throughput sequencing

MicroRNAs (miRNAs) have been reported to be correlated with various stress responses in soybean, but only a few miRNAs have been demonstrated to respond to low phosphorus (LP) stress. To unravel the response mechanisms of miRNAs to low-P stress, the roots of two representative soybean genotypes with different P efficiency, Nannong94-156 (a LP-tolerant genotype) and Bogao (a LP-sensitive genotype), were used for the construction of RNA sequencing (RNA-seq) libraries under low/normal-P treatment by high-throughput sequencing. In total, 603 existing miRNAs and 1699 novel miRNAs belonging to 248 and 1582 families in all samples were identified, respectively. Among these miRNAs, 777 miRNAs were differentially expressed (DE) across different P levels and genotypes. Furthermore, putative targets of DE miRNAs were predicted, and these miRNAs mainly targeted ERF (ethylene responsive factor), auxin response factors (ARF), zinc finger protein, MYB, and NAC domain transcription factors. Gene ontology (GO) analysis showed that targets of DE miRNAs were significantly enriched in binding, metabolic processes, biological regulation, response to stress, and phosphorus metabolic processes. In addition, the expression profiles of chosen P-responsive miRNAs and target genes were validated by quantitative real-time PCR (qRT-PCR). Our study focused on genome-wide miRNA identification in two representative soybean genotypes under low-P stress. Overall, the DE miRNAs across different P levels and genotypes and their putative target genes will provide useful information for further study of miRNAs mediating low-P response and facilitate improvements in soybean breeding.

Genome-Wide Association Analysis for Phosphorus Use Efficiency Traits in Mungbean (Vigna radiata L. Wilczek) Using Genotyping by Sequencing Approach

Mungbean (Vigna radiata L. Wilczek) is an annual grain legume crop affected by low availability of phosphorus. Phosphorus deficiency mainly affects the growth and development of plants along with changes in root morphology and increase in root-to-shoot ratio. Deciphering the genetic basis of phosphorus use efficiency (PUE) traits can benefit our understanding of mungbean tolerance to low-phosphorus condition. To address this issue, 144 diverse mungbean genotypes were evaluated for 12 PUE traits under hydroponics with optimum- and low-phosphorus levels. The broad sense heritability of traits ranged from 0.63 to 0.92 and 0.58 to 0.92 under optimum- and low-phosphorus conditions, respectively. This study, reports for the first time such a large number of genome wide Single nucleotide polymorphisms (SNPs) (76,160) in mungbean. Further, genome wide association study was conducted using 55,634 SNPs obtained by genotyping-by-sequencing method. The results indicated that total 136 SNPs shared by both GLM and MLM models were associated with tested PUE traits under different phosphorus regimes. We have identified SNPs with highest p value (-log10(p)) for some traits like, TLA and RDW with p value (-log10(p)) of more than 6.0 at LP/OP and OP condition. We have identified nine SNPs (three for TLA and six for RDW trait) which was found to be present in chromosomes 8, 4, and 7. One SNP present in Vradi07g06230 gene contains zinc finger CCCH domain. In total, 71 protein coding genes were identified, of which 13 genes were found to be putative candidate genes controlling PUE by regulating nutrient uptake and root architectural development pathways in mungbean. Moreover, we identified three potential candidate genes VRADI11G08340, VRADI01G05520, and VRADI04G10750 with missense SNPs in coding sequence region, which results in significant variation in protein structure at tertiary level. The identified SNPs and candidate genes provide the essential information for genetic studies and marker-assisted breeding program for improving low-phosphorus tolerance in mungbean.

High-density QTL mapping of leaf-related traits and chlorophyll content in three soybean RIL populations

BACKGROUND

Leaf size and shape, which affect light capture, and chlorophyll content are important factors affecting photosynthetic efficiency. Genetic variation of these components significantly affects yield potential and seed quality. Identification of the genetic basis for these traits and the relationship between them is of great practical significance for achieving ideal plant architecture and high photosynthetic efficiency for improved yield.

RESULTS

Here, we undertook a large-scale linkage mapping study using three mapping populations to determine the genetic interplay between soybean leaf-related traits and chlorophyll content across two environments. Correlation analysis revealed a significant negative correlation between leaf size and shape, while both traits were positively correlated with chlorophyll content. This phenotypic relationship was verified across the three mapping populations as determined by principal component analysis, suggesting that these traits are under the control of complex and interrelated genetic components. The QTLs for leaf-related traits and chlorophyll are partly shared, which further supports the close genetic relationship between the two traits. The largest-effect major loci, q20, was stably identified across all population and environments and harbored the narrow leaflet gene Gm-JAG1 (Ln/ln), which is a key regulator of leaflet shape in soybean.

CONCLUSION

Our results uncover several major QTLs (q4-1, q4-2, q11, q13, q18 and q20) and its candidate genes specific or common to leaf-related traits and chlorophyll, and also show a complex epistatic interaction between the two traits. The SNP markers closely linked to these valuable QTLs could be used for molecular design breeding with improved plant architecture, photosynthetic capacity and even yield.

Genome-Wide Analysis Reveals Dynamic Epigenomic Differences in Soybean Response to Low-Phosphorus Stress

Low-phosphorus (low-P) stress has a significant limiting effect on crop yield and quality. Although the molecular mechanisms of the transcriptional level responsible for the low-P stress response have been studied in detail, the underlying epigenetic mechanisms in gene regulation remain largely unknown. In this study, we evaluated the changes in DNA methylation, gene expression and small interfering RNAs (siRNAs) abundance genome-wide in response to low-P stress in two representative soybean genotypes with different P-efficiencies. The DNA methylation levels were slightly higher under low-P stress in both genotypes. Integrative methylation and transcription analysis suggested a complex regulatory relationship between DNA methylation and gene expression that may be associated with the type, region, and extent of methylation. Association analysis of low-P-induced differential methylation and gene expression showed that transcriptional alterations of a small part of genes were associated with methylation changes. Dynamic methylation alterations in transposable element (TE) regions in the CHH methylation context correspond with changes in the amount of siRNA under low-P conditions, indicating an important role of siRNAs in modulating TE activity by guiding CHH methylation in TE regions. Together, these results could help to elucidate the epigenetic regulation mechanisms governing the responses of plants to abiotic stresses.

Up-regulating GmETO1 improves phosphorus uptake and use efficiency by promoting root growth in soybean

Soybean is a high inorganic phosphate (Pi) demanding crop; its production is strongly suppressed when Pi is deficient in soil. However, the regulatory mechanism of Pi deficiency tolerance in soybean is still largely unclear. Here, our findings highlighted the pivotal role of the ethylene-associated pathway in soybean tolerance to Pi deficiency by comparatively studying transcriptome changes between a representative Pi-deficiency-tolerant soybean genotype NN94156 and a sensitive genotype Bogao under different Pi supplies. By further integrating high-confident linkage and association mapping, we identified that Ethylene-Overproduction Protein 1 (GmETO1), an essential ethylene-biosynthesis regulator, underlies the major quantitative trait locus (QTL) q14-2 controlling Pi uptake. GmETO1 was also the representative member of ETO1 family members that was strongly induced by Pi deficiency. Overexpressing GmETO1 significantly enhanced Pi deficiency tolerance by increasing proliferation and elongation of hairy roots, Pi uptake and use efficiency, and conversely, silencing of GmETO1 led to opposite findings. We further demonstrated that Pi-deficiency inducible genes critical for root morphological and physiological traits including GmACP1/2, Pht1;4, Expansin-A7 and Root Primordium Defective 1 functioned downstream of GmETO1. Our study provides comprehensive insight into the complex regulatory mechanism of Pi deficiency tolerance in soybean and a potential way to genetically improve soybean low-Pi tolerance.

Scion genotypes exert long distance control over rootstock transcriptome responses to low phosphate in grafted grapevine

BACKGROUND

Grafting is widely used in horticulture and rootstocks are known to modify scion growth and adaptation to soil conditions. However, the role of scion genotype in regulating rootstock development and functioning has remained largely unexplored. In this study, reciprocal grafts of two grapevine genotypes were produced as well as the corresponding homo-graft controls. These plants were subjected to a low phosphate (LP) treatment and transcriptome profiling by RNA sequencing was done on root samples collected 27 h after the onset of the LP treatment.

RESULTS

A set of transcripts responsive to the LP treatment in all scion/rootstock combinations was identified. Gene expression patterns associated with genetic variation in response to LP were identified by comparing the response of the two homo-grafts. In addition, the scion was shown to modify root transcriptome responses to LP in a rootstock dependent manner. A weighted gene co-expression network analysis identified modules of correlated genes; the analysis of the association of these modules with the phosphate treatment, and the scion and rootstock genotype identified potential hub genes.

CONCLUSIONS

This study provides insights into the response of grafted grapevine to phosphate supply and identifies potential shoot-to-root signals that could vary between different grapevine genotypes.

Dynamic gene expression changes in response to micronutrient, macronutrient, and multiple stress exposures in soybean

Preserving crop yield is critical for US soybean production and the global economy. Crop species have been selected for increased yield for thousands of years with individual lines selected for improved performance in unique environments, constraints not experienced by model species such as Arabidopsis. This selection likely resulted in novel stress adaptations, unique to crop species. Given that iron deficiency is a perennial problem in the soybean growing regions of the USA and phosphate deficiency looms as a limitation to global agricultural production, nutrient stress studies in crop species are critically important. In this study, we directly compared whole-genome expression responses of leaves and roots to iron (Fe) and phosphate (Pi) deficiency, representing a micronutrient and macronutrient, respectively. Conducting experiments side by side, we observed soybean responds to both nutrient deficiencies within 24 h. While soybean responds largely to -Fe deficiency, it responds strongly to Pi resupply. Though the timing of the responses was different, both nutrient stress signals used the same molecular pathways. Our study is the first to demonstrate the speed and diversity of the soybean stress response to multiple nutrient deficiencies. We also designed the study to examine gene expression changes in response to multiple stress events. We identified 865 and 3375 genes that either altered their direction of expression after a second stress exposure or were only differentially expressed after a second stress event. Understanding the molecular underpinnings of these responses in crop species could have major implications for improving stress tolerance and preserving yield.

Transcriptome-wide identification of novel circular RNAs in soybean in response to low-phosphorus stress

Low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Circular RNAs (circRNAs) are novel noncoding RNAs that play a crucial role in plant responses to abiotic stress. However, how LP stress mediates the biogenesis of circRNAs in soybean remains unclear. Here, to explore the response mechanisms of circRNAs to LP stress, the roots of two representative soybean genotypes with different P-use efficiency, Bogao (a LP-sensitive genotype) and Nannong 94156 (a LP-tolerant genotype), were used for the construction of RNA sequencing (RNA-seq) libraries and circRNA identification. In total, 371 novel circRNA candidates, including 120 significantly differentially expressed (DE) circRNAs, were identified across different P levels and genotypes. More DE circRNAs were significantly regulated by LP stress in Bogao than in NN94156, suggesting that the tolerant genotype was less affected by LP stress than the sensitive genotype was; in other words, NN94156 may have a better ability to maintain P homeostasis under LP stress. Moreover, a positive correlation was observed between the expression patterns of P stress-induced circRNAs and their circRNA-host genes. Gene Ontology (GO) enrichment analysis of these circRNA-host genes and microRNA (miRNA)-targeted genes indicated that these DE circRNAs were involved mainly in defense responses, ADP binding, nucleoside binding, organic substance catabolic processes, oxidoreductase activity, and signal transduction. Together, our results revealed that LP stress can significantly alter the genome-wide profiles of circRNAs and indicated that the regulation of circRNAs was both genotype and environment specific in response to LP stress. LP-induced circRNAs might provide a rich resource for LP-responsive circRNA candidates for future studies.

Divergent metabolic adjustments in nodules are indispensable for efficient N2 fixation of soybean under phosphate stress

The main objective of the present study was to characterize the symbiotic N2 fixation (SNF) capacity and to elucidate the underlying mechanisms for low-Pi acclimation in soybean plants grown in association with two Bradyrhizobium diazoefficiens strains which differ in SNF capacity (USDA110 vs. CB1809). In comparison with the USDA110-soybean, the CB1809-soybean association revealed a greater SNF capacity in response to Pi starvation, as evidenced by relative higher plant growth and higher expression levels of the nifHDK genes. This enhanced Pi acclimation was partially related to the efficient utilization to the overall carbon (C) budget of symbiosis in the CB1809-induced nodules compared with that of the USDA110-induced nodules under low-Pi provision. In contrast, the USDA110-induced nodules favored other metabolic acclimation mechanisms that expend substantial C cost, and consequently cause negative implications on nodule C expenditure during low-Pi conditions. Fatty acids, phytosterols and secondary metabolites are characterized among the metabolic pathways involved in nodule acclimation under Pi starvation. While USDA110-soybean association performed better under Pi sufficiency, it is very likely that the CB1809-soybean association is better acclimatized to cope with Pi deficiency owing to the more effective functional plasticity and lower C cost associated with these nodular metabolic arrangements.

Artificial selection on GmOLEO1 contributes to the increase in seed oil during soybean domestication

Increasing seed oil content is one of the most important breeding goals for soybean due to a high global demand for edible vegetable oil. However, genetic improvement of seed oil content has been difficult in soybean because of the complexity of oil metabolism. Determining the major variants and molecular mechanisms conferring oil accumulation is critical for substantial oil enhancement in soybean and other oilseed crops. In this study, we evaluated the seed oil contents of 219 diverse soybean accessions across six different environments and dissected the underlying mechanism using a high-resolution genome-wide association study (GWAS). An environmentally stable quantitative trait locus (QTL), GqOil20, significantly associated with oil content was identified, accounting for 23.70% of the total phenotypic variance of seed oil across multiple environments. Haplotype and expression analyses indicate that an oleosin protein-encoding gene (GmOLEO1), colocated with a leading single nucleotide polymorphism (SNP) from the GWAS, was significantly correlated with seed oil content. GmOLEO1 is predominantly expressed during seed maturation, and GmOLEO1 is localized to accumulated oil bodies (OBs) in maturing seeds. Overexpression of GmOLEO1 significantly enriched smaller OBs and increased seed oil content by 10.6% compared with those of control seeds. A time-course transcriptomics analysis between transgenic and control soybeans indicated that GmOLEO1 positively enhanced oil accumulation by affecting triacylglycerol metabolism. Our results also showed that strong artificial selection had occurred in the promoter region of GmOLEO1, which resulted in its high expression in cultivated soybean relative to wild soybean, leading to increased seed oil accumulation. The GmOLEO1 locus may serve as a direct target for both genetic engineering and selection for soybean oil improvement.

The EIL transcription factor family in soybean: Genome-wide identification, expression profiling and genetic diversity analysis

The ETHYLENE INSENSITIVE3-LIKE (EIL) transcription factor family plays a critical role in the ethylene signaling pathway, which regulates a broad spectrum of plant growth and developmental processes, as well as defenses to myriad stresses. Although genome-wide analysis of this family has been carried out for several plant species, no comprehensive analysis of the EIL gene family in soybean has been reported so far. Furthermore, there are few studies on the functions of EIL genes in soybean. In this study, we identified 12 soybean (Gm) EIL genes, which we divided into three groups based on their phylogenetic relationships. We then detected their duplication status and found that most of the GmEIL genes have duplicated copies derived from two whole-genome duplication events. These duplicated genes underwent strong negative selection during evolution. We further analyzed the transcript profiles of GmEIL genes using the transcriptome data and found that their spatio-temporal and stress expression patterns varied considerably. For example, GmEIL1-GmEIL5 were found to be strongly expressed in almost every sample, while GmEIL8-GmEIL12 exhibited low expression, or were not expressed at all. Additionally, these genes showed different responses to dehydration, salinity and phosphate starvation. Finally, we surveyed genetic variations of these genes in 302 resequenced wild soybeans, landraces and improved soybean cultivars. Our data showed that most GmEIL genes are well conserved, and are not modified in domesticated or improved cultivars. Together, these findings provide a potentially valuable resource for characterizing the GmEIL gene family and lay the basis for further elucidation of their molecular mechanisms.

Combining QTL mapping with transcriptome and metabolome profiling reveals a possible role for ABA signaling in resistance against the cabbage whitefly in cabbage

Whiteflies are among the world's most significant agricultural pests and chemical insecticides are extensively used to reduce crop damage to acceptable levels. However, nearly all insecticides pose a threat to the environment and alternative control methods, such as breeding of crop varieties that are inherently insect-resistant, are needed. Previously, a strong source of plant-age dependent resistance to the cabbage whitefly (Aleyrodes proletella) has been identified in the modern white cabbage (Brassica oleracea var. capitata) variety Rivera. However, nothing is known about the molecular mechanisms or the genes involved in this resistance. In the present study, a multidisciplinary approach combining transcriptome and metabolome profiling with genetic mapping was used to identify the molecular players of whitefly resistance in cabbage. Transcriptome profiles of young (susceptible) and older (resistant) Rivera plants were analyzed using RNA sequencing. While many genes involved in general processes were differentially expressed between both ages, several defense-related processes were overrepresented in the transcriptome profile of older plants. Hormone measurements revealed that jasmonic acid (JA) levels decreased upon whitefly infestation at both plant ages. Interestingly, abscisic acid (ABA) levels showed contrasting effects in response to whitefly infestation: ABA levels were reduced in young plants but induced in older plants upon whitefly feeding. Auxin levels were significantly lower in older plants compared with young plants, independent of whitefly presence, while glucosinolate levels were higher. Additionally, whitefly performance was monitored in an F2 population derived from a cross between Rivera and the susceptible white cabbage variety Christmas Drumhead. Significant QTL intervals were mapped on chromosome 2 and 9 for oviposition rate and whitefly adult survival, respectively. Several genes that were higher expressed in older plants and located in the identified QTL intervals were orthologous to Arabidopsis genes that have been related to ABA signaling, suggesting a role for ABA in the regulation of resistance towards whiteflies. Our results show that combining different omics approaches is a useful strategy to identify candidate genes underlying insect resistance.

Acid phosphatase gene GmHAD1 linked to low phosphorus tolerance in soybean, through fine mapping

KEY MESSAGE

Map-based cloning identified GmHAD1, a gene which encodes a HAD-like acid phosphatase, associated with soybean tolerance to low phosphorus stress. Phosphorus (P) deficiency in soils is a major limiting factor for crop growth worldwide. Plants may adapt to low phosphorus (LP) conditions via changes to root morphology, including the number, length, orientation, and branching of the principal root classes. To elucidate the genetic mechanisms for LP tolerance in soybean, quantitative trait loci (QTL) related to root morphology responses to LP were identified via hydroponic experiments. In total, we identified 14 major loci associated with these traits in a RIL population. The log-likelihood scores ranged from 2.81 to 7.43, explaining 4.23-13.98% of phenotypic variance. A major locus on chromosome 08, named qP8-2, was co-localized with an important P efficiency QTL (qPE8), containing phosphatase genes GmACP1 and GmACP2. Another major locus on chromosome 10 named qP10-2 explained 4.80-13.98% of the total phenotypic variance in root morphology. The qP10-2 contains GmHAD1, a gene which encodes an acid phosphatase. In the transgenic soybean hairy roots, GmHAD1 overexpression increased P efficiency by 8.4-16.5% relative to the control. Transgenic Arabidopsis plants had higher biomass than wild-type plants across both short- and long-term P reduction. These results suggest that GmHAD1, an acid phosphatase gene, improved the utilization of organic phosphate by soybean and Arabidopsis plants.

Early Transcriptomic Response to Phosphate Deprivation in Soybean Leaves as Revealed by RNA-Sequencing

Low phosphate (Pi) availability is an important limiting factor affecting soybean production. However, the underlying molecular mechanisms responsible for low Pi stress response and tolerance remain largely unknown, especially for the early signaling events under low Pi stress. Here, a genome-wide transcriptomic analysis in soybean leaves treated with a short-term Pi-deprivation (24 h) was performed through high-throughput RNA sequencing (RNA-seq) technology. A total of 533 loci were found to be differentially expressed in response to Pi deprivation, including 36 mis-annotated loci and 32 novel loci. Among the differentially expressed genes (DEGs), 303 were induced and 230 were repressed by Pi deprivation. To validate the reliability of the RNA-seq data, 18 DEGs were randomly selected and analyzed by quantitative RT-PCR (reverse transcription polymerase chain reaction), which exhibited similar fold changes with RNA-seq. Enrichment analyses showed that 29 GO (Gene Ontology) terms and 8 KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were significantly enriched in the up-regulated DEGs and 25 GO terms and 16 KEGG pathways were significantly enriched in the down-regulated DEGs. Some DEGs potentially involved in Pi sensing and signaling were up-regulated by short-term Pi deprivation, including five SPX-containing genes. Some DEGs possibly associated with water and nutrient uptake, hormonal and calcium signaling, protein phosphorylation and dephosphorylation and cell wall modification were affected at the early stage of Pi deprivation. The cis-elements of PHO (phosphatase) element, PHO-like element and P responsive element were present more frequently in promoter regions of up-regulated DEGs compared to that of randomly-selected genes in the soybean genome. Our transcriptomic data showed an intricate network containing transporters, transcription factors, kinases and phosphatases, hormone and calcium signaling components is involved in plant responses to early Pi deprivation.

Physiological and Proteomics Analyses Reveal Low-Phosphorus Stress Affected the Regulation of Photosynthesis in Soybean

Previous studies have revealed a significant genetic relationship between phosphorus (P)-efficiency and photosynthesis-related traits in soybean. In this study, we used proteome profiling in combination with expression analysis, biochemical investigations, and leaf ultrastructural analysis to identify the underlying physiological and molecular responses. The expression analysis and ultrastructural analysis showed that the photosynthesis key genes were decreased at transcript levels and the leaf mesophyll and chloroplast were severely damaged after low-P stress. Approximately 55 protein spots showed changes under low-P condition by mass spectrometry, of which 17 were involved in various photosynthetic processes. Further analysis revealed the depression of photosynthesis caused by low-P stress mainly involves the regulation of leaf structure, adenosine triphosphate (ATP) synthesis, absorption and transportation of CO₂, photosynthetic electron transport, production of assimilatory power, and levels of enzymes related to the Calvin cycle. In summary, our findings indicated that the existence of a stringent relationship between P supply and the genomic control of photosynthesis in soybean. As an important strategy to protect soybean photosynthesis, P could maintain the stability of cell structure, up-regulate the enzymes’ activities, recover the process of photosystem II (PSII), and induce the expression of low-P responsive genes and proteins.

Unraveling lipid metabolism in maize with time-resolved multi-omics data

Maize is the cereal crop with the highest production worldwide, and its oil is a key energy resource. Improving the quantity and quality of maize oil requires a better understanding of lipid metabolism. To predict the function of maize genes involved in lipid biosynthesis, we assembled transcriptomic and lipidomic data sets from leaves of B73 and the high-oil line By804 in two distinct time-series experiments. The integrative analysis based on high-dimensional regularized regression yielded lipid-transcript associations indirectly validated by Gene Ontology and promoter motif enrichment analyses. The co-localization of lipid-transcript associations using the genetic mapping of lipid traits in leaves and seedlings of a B73 × By804 recombinant inbred line population uncovered 323 genes involved in the metabolism of phospholipids, galactolipids, sulfolipids and glycerolipids. The resulting association network further supported the involvement of 50 gene candidates in modulating levels of representatives from multiple acyl-lipid classes. Therefore, the proposed approach provides high-confidence candidates for experimental testing in maize and model plant species.