An RNA-Seq transcriptome analysis of orthophosphate-deficient white lupin reveals novel insights into phosphorus acclimation in plants

Transcriptomics Provide Insights into Early Responses to Sucrose Signaling in Lupinus albus, a Model Plant for Adaptations to Phosphorus and Iron Deficiency

Phosphorus (P) and iron (Fe) deficiency are major limiting factors for plant productivity worldwide. White lupin (Lupinus albus L.) has become a model plant for understanding plant adaptations to P and Fe deficiency, because of its ability to form cluster roots, bottle-brush-like root structures play an important role in the uptake of P and Fe from soil. However, little is known about the signaling pathways involved in sensing and responding to P and Fe deficiency. Sucrose, sent in increased concentrations from the shoot to the root, has been identified as a long-distance signal of both P and Fe deficiency. To unravel the responses to sucrose as a signal, we performed Oxford Nanopore cDNA sequencing of white lupin roots treated with sucrose for 10, 15, or 20 min compared to untreated controls. We identified a set of 17 genes, including 2 bHLH transcription factors, that were up-regulated at all three time points of sucrose treatment. GO (gene ontology) analysis revealed enrichment of auxin and gibberellin responses as early as 10 min after sucrose addition, as well as the emerging of ethylene responses at 20 min of sucrose treatment, indicating a sequential involvement of these hormones in plant responses to sucrose.

Exploring the dynamic adaptive responses of Epimedium pubescens to phosphorus deficiency by Integrated transcriptome and miRNA analysis

Phosphorus, a crucial macronutrient essential for plant growth and development. Due to widespread phosphorus deficiency in soils, phosphorus deficiency stress has become one of the major abiotic stresses that plants encounter. Despite the evolution of adaptive mechanisms in plants to address phosphorus deficiency, the specific strategies employed by species such as Epimedium pubescens remain elusive. Therefore, this study observed the changes in the growth, physiological reponses, and active components accumulation in E. pubescensunder phosphorus deficiency treatment, and integrated transcriptome and miRNA analysis, so as to offer comprehensive insights into the adaptive mechanisms employed by E. pubescens in response to phosphorus deficiency across various stages of phosphorus treatment. Remarkably, our findings indicate that phosphorus deficiency induces root growth stimulation in E. pubescens, while concurrently inhibiting the growth of leaves, which are of medicinal value. Surprisingly, this stressful condition results in an augmented accumulation of active components in the leaves. During the early stages (30 days), leaves respond by upregulating genes associated with carbon metabolism, flavonoid biosynthesis, and hormone signaling. This adaptive response facilitates energy production, ROS scavenging, and morphological adjustments to cope with short-term phosphorus deficiency and sustain its growth. As time progresses (90 days), the expression of genes related to phosphorus cycling and recycling in leaves is upregulated, and transcriptional and post-transcriptional regulation (miRNA regulation and protein modification) is enhanced. Simultaneously, plant growth is further suppressed, and it gradually begins to discard and decompose leaves to resist the challenges of long-term phosphorus deficiency stress and sustain survival. In conclusion, our study deeply and comprehensively reveals adaptive strategies utilized by E. pubescens in response to phosphorus deficiency, demonstrating its resilience and thriving potential under stressful conditions. Furthermore, it provides valuable information on potential target genes for the cultivation of E. pubescens genotypes tolerant to low phosphorus.

Revealing the Mechanistic Basis of Regulation of Phosphorus Uptake in Soybean (Glycine max) Roots by Molybdenum: An Integrated Omics Approach

While molybdenum (Mo) application can improve phosphorus (P) availability to plants by changing P speciation in the rhizosphere, the mechanistic basis of this process remains unclear. This work investigated the impact of various combinations of Mo and P treatments on root morphology, P and Mo uptake, and root transcriptome and metabolome. Mo application significantly increased soybean biomass and the number of lateral roots at both low (5 μmol) or normal (500 μmol) P levels and significantly improved P concentration and accumulation in Normal P treatment. Compared with the Normal P treatment, Low P significantly increased the number of roots, root surface area, and root acid phosphatase secretion. A total of 6811 Mo-responsive differentially expressed genes and 135 differential metabolites were identified at two P levels. At Low P, transcriptional changes significantly increased root synthesis and secretion of succinic acid, methylmalonic acid, and other organic acids as well as acid phosphatase, thereby increasing the conversion of soil aluminum-bound P and organic P into available P. At Normal P, Mo application increased P uptake mainly by increasing the number of lateral roots. Thus, Mo helps crops adapt to different P levels by regulating root anatomy and transcriptional and metabolic profiles of their roots.

Carbon Fluxes in Potato (Solanum tuberosum) Remain Stable in Cell Cultures Exposed to Nutritional Phosphate Deficiency

Nutritional phosphate deficiency is a major limitation to plant growth. Here, we monitored fluxes in pathways supporting respiratory metabolism in potato (Solanum tuberosum) cell cultures growing in control or limiting phosphate conditions. Sugar uptake was quantified using [U-14C]sucrose as precursor. Carbohydrate degradation through glycolysis and respiratory pathways was estimated using the catabolism of [U-14C]sucrose to 14CO2. Anaplerotic carbon flux was assessed by labeling with NaH14CO3. The data showed that these metabolic fluxes displayed distinct patterns over culture time. However, phosphate depletion had relatively little impact on the various fluxes. Sucrose uptake was higher during the first six days of culture, followed by a decline, which was steeper in Pi-sufficient cells. Anaplerotic pathway flux was more important at day three and decreased thereafter. In contrast, the flux between sucrose and CO2 was at a maximum in the mid-log phase of the culture, with a peak at Day 6. Metabolization of [U-14C]sucrose into neutral, basic and acidic fractions was also unaffected by phosphate nutrition. Hence, the well-documented changes in central metabolism enzymes activities in response to Pi deficiency do not drastically modify metabolic fluxes, but rather result in the maintenance of the carbon fluxes that support respiration.

Role of SiPHR1 in the Response to Low Phosphate in Foxtail Millet via Comparative Transcriptomic and Co-Expression Network Analyses

Enhancing the absorption and utilization of phosphorus by crops is an important aim for ensuring food security worldwide. However, the gene regulatory network underlying phosphorus use in foxtail millet remains unclear. In this study, the molecular mechanism underlying low-phosphorus (LP) responsiveness in foxtail millet was evaluated using a comparative transcriptome analysis. LP reduced the chlorophyll content in shoots, increased the anthocyanin content in roots, and up-regulated purple acid phosphatase and phytase activities as well as antioxidant systems (CAT, POD, and SOD). Finally, 13 differentially expressed genes related to LP response were identified and verified using transcriptomic data and qRT-PCR. Two gene co-expression network modules related to phosphorus responsiveness were positively correlated with POD, CAT, and PAPs. Of these, SiPHR1, functionally annotated as PHOSPHATE STARVATION RESPONSE 1, was identified as an MYB transcription factor related to phosphate responsiveness. SiPHR1 overexpression in Arabidopsis significantly modified the root architecture. LP stress caused cellular, physiological, and phenotypic changes in seedlings. SiPHR1 functioned as a positive regulator by activating downstream genes related to LP tolerance. These results improve our understanding of the molecular mechanism underlying responsiveness to LP stress, thereby laying a theoretical foundation for the genetic modification and breeding of new LP-tolerant foxtail millet varieties.

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

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

Transcriptome analysis provides insights into the response of Lotus corniculatus roots to low-phosphorus stress

INTRODUCTION

A lack of soil phosphorus (P) is a principal factor restricting the normal growth of Lotus corniculatus in the karst area of Guizhou Province, China, but the response mechanism of L. corniculatus under low-phosphorus stress remains unclear.

METHODS

Therefore, we treated two selected L. corniculatus lines (low-P-intolerant line 08518 and low-P-tolerant line 01549) from 13 L. corniculatus lines with normal phosphorus (0.5 mmol/L KH₂PO₄, NP) and low phosphorus (0.005 mmol/L KH₂PO₄, LP) concentrations to study changes in morphological, physiological and transcriptome data under low-phosphorus stress.

RESULTS

The low-P-tolerant line 01549 exhibited better performance under low-phosphorus stress. Compared with the NP treatment, all root morphological indicators of the low-P-tolerant line 01549 increased, and those of the low-P-intolerant line 08518 decreased under low-P stress. Compared with the NP treatment, acid phosphatase (ACP), catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) activities, and the malondialdehyde (MDA), soluble sugar (SS), soluble protein (SP) and proline (Pro) contents of the two L. corniculatus lines increased under low-P stress. A transcriptome analysis of L. corniculatus showed that a total of 656 and 2243 differentially expressed genes (DEGs) were identified in line 01549 and line 08518, respectively. Meanwhile, the main pathways, such as carbohydrate metabolism, acid phosphatases, phosphate transporters and biosynthesis of secondary metabolites, as well as related genes were also screened by performing a KEGG enrichment analysis.

DISCUSSION

The findings provide an essential point of reference for studying the physiological and molecular mechanism of the response to low-P stress in L. corniculatus.

Genome-wide association analysis discovered new loci and candidate genes associated with low-phosphorus tolerance based on shoot mineral elements concentrations in soybean

Phosphorus (P) deficiency is one of the major limitations for soybean production. Moreover, it has been well reported P and other mineral elements function interdependently or antagonistically to control nutrients homeostasis in plants. Thus, it is urgently needed to understand the genetic mechanism of the accumulation of mineral elements in response to low-P stress. In this study, to identify single nucleotide polymorphisms (SNPs) and candidate genes controlling the accumulation of mineral elements suffering low-P stress in seedling stage of soybean plants, we measured concentrations of mineral elements, including P, Zn, Fe, Mn, Mg and Ca, in shoots of 211 soybean accessions under normal phosphorus (+P) and low phosphorus (-P) conditions in two hydroponic experiments. And genome-wide association study (GWAS) using high density NJAU 355K SoySNP array and concentrations of five of these mineral elements except P was performed. A total of 36 SNPs distributed on 13 chromosomes were identified to be significantly associated with low-P tolerance, and nine SNPs on chromosome 10 formed a SNP cluster. Meanwhile, the candidate gene GmFeB1 was found to serve as a negative regulator element involved in soybean P metabolism and the haplotype1 (Hap1) of GmFeB1 showed significantly higher shoot Fe concentration under -P condition than that of Hap2. In summary, we uncover 36 SNPs significantly associated with shoot mineral elements concentrations under different P conditions and a soybean low-P related gene GmFeB1, which will provide additional genetic information for soybean low-P tolerance and new gene resources for P-efficient soybean varieties breeding.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Relationships Among Arsenic-Related Traits, Including Rice Grain Arsenic Concentration and Straighthead Resistance, as Revealed by Genome-Wide Association

There is global concern that rice grains and foods can contain harmful amounts of arsenic (As), motivating breeders to produce cultivars that restrict As accumulation in grains to protect human health. Arsenic is also toxic to plants, with straighthead disorder (StHD), causing panicle sterility, being observed in rice. The genetic variation in StHD resistance suggests that plants have evolved mechanisms that reduce As toxicity, possibly via regulation of As uptake, transport, or detoxification/sequestration. Because these mechanisms could also underlie the wide (3- to 100-fold) differences in grain As concentration (grain-As) observed among diverse rice genotypes, it was hypothesized that some genes reduce both grain-As content and StHD susceptibility and may be detectable as co-located StDH and As quantitative trait loci (QTL). We used a machine-learning Bayesian network approach plus high-resolution genome-wide association study (GWAS) to identify QTL for grain-As and StHD resistance within the USDA Rice Minicore Collection (RMC). Arsenic enters roots through phosphorus (P) and silica (Si) transporters, As detoxification involves sulfur (S), and cell signaling to activate stress tolerance mechanisms is impacted by Si, calcium (Ca), and copper (Cu). Therefore, concentrations of Si, P, S, Ca, and Cu were included in this study to elucidate physiological mechanisms underlying grain-As and StHD QTL. Multiple QTL (from 9 to 33) were identified for each of the investigated As-associated traits. Although the QTL for StHD, Si, and grain-As did not overlap as heavily as our hypothesis predicted (4/33 StHD and 4/15 As QTL co-located), they do provide useful guidance to future research. Furthermore, these are the first StHD and Si QTL to be identified using high-density mapping, resulting in their being mapped to shorter, more precise genomic regions than previously reported QTL. The candidate genes identified provide guidance for future research, such as gene editing or mutation studies to further investigate the role of antioxidants and ROS scavenging to StHD resistance, as indicated by candidate genes around the commonly reported qStHD8-2 QTL. Other genes indicated for future study for improving grain-As and StHD include several multidrug and toxic compound extrusion (MATE) genes, F-box genes, and NIPs not documented to date to transport As.

Improving phosphate use efficiency in the aquatic crop watercress (Nasturtium officinale)

Watercress is a nutrient-dense leafy green crop, traditionally grown in aquatic outdoor systems and increasingly seen as well-suited for indoor hydroponic systems. However, there is concern that this crop has a detrimental impact on the environment through direct phosphate additions causing environmental pollution. Phosphate-based fertilisers are supplied to enhanced crop yield, but their use may contribute to eutrophication of waterways downstream of traditional watercress farms. One option is to develop a more phosphate use efficient (PUE) crop. This review identifies the key traits for this aquatic crop (the ideotype), for future selection, marker development and breeding. Traits identified as important for PUE are (i) increased root surface area through prolific root branching and adventitious root formation, (ii) aerenchyma formation and root hair growth. Functional genomic traits for improved PUE are (iii) efficacious phosphate remobilisation and scavenging strategies and (iv) the use of alternative metabolic pathways. Key genomic targets for this aquatic crop are identified as: PHT phosphate transporter genes, global transcriptional regulators such as those of the SPX family and genes involved in galactolipid and sulfolipid biosynthesis such as MGD2/3, PECP1, PSR2, PLDζ1/2 and SQD2. Breeding for enhanced PUE in watercress will be accelerated by improved molecular genetic resources such as a full reference genome sequence that is currently in development.

Integrative analysis of the metabolome and transcriptome reveal the phosphate deficiency response pathways of alfalfa

Understanding the mechanisms underlying the responses to inorganic phosphate (Pi) deficiency in alfalfa will help enhance Pi acquisition efficiency and the sustainable use of phosphorous resources. Integrated global metabolomic and transcriptomic analyses of mid-vegetative alfalfa seedlings under 12-day Pi deficiency were conducted. Limited seedling growth were found, including 13.24%, 16.85% and 33.36% decreases in height, root length and photosynthesis, and a 24.10% increase in root-to-shoot ratio on day 12. A total of 322 and 448 differentially abundant metabolites and 1199 and 1061 differentially expressed genes were identified in roots and shoots. Increased (>3.68-fold) inorganic phosphate transporter 1;4 and SPX proteins levels in the roots (>2.15-fold) and shoots (>2.50-fold) were related to Pi absorption and translocation. The levels of phospholipids and Pi-binding carbohydrates and nucleosides were decreased, while those of phosphatases and pyrophosphatases in whole seedlings were induced under reduced Pi. In addition, nitrogen assimilation was affected by inhibiting high-affinity nitrate transporters (NRT2.1 and NRT3.1), and nitrate reductase. Increased delphinidin-3-glucoside might contribute to the gray-green leaves induced by Pi limitation. Stress-induced MYB, WRKY and ERF transcription factors were identified. The responses of alfalfa to Pi deficiency were summarized as local systemic signaling pathways, including root growth, stress-related responses consisting of enzymatic and nonenzymatic systems, and hormone signaling and systemic signaling pathways including Pi recycling and Pi sensing in the whole plant, as well as Pi recovery, and nitrate and metal absorption in the roots. This study provides important information on the molecular mechanism of the response to Pi deficiency in alfalfa.

Loss of LaMATE impairs isoflavonoid release from cluster roots of phosphorus-deficient white lupin

White lupin (Lupinus albus L.) forms brush-like root structures called cluster roots under phosphorus-deficient conditions. Clusters secrete citrate and other organic compounds to mobilize sparingly soluble soil phosphates. In the context of aluminum toxicity tolerance mechanisms in other species, citrate is released via a subgroup of MATE/DTX proteins (multidrug and toxic compound extrusion/detoxification). White lupin contains 56 MATE/DTX genes. Many of these are closely related to gene orthologs with known substrates in other species. LaMATE is a marker gene for functional, mature clusters and is, together with its close homolog LaMATE3, a candidate for the citrate release. Both were highest expressed in mature clusters and when expressed in oocytes, induced inward-rectifying currents that were likely carried by endogenous channels. No citrate efflux was associated with LaMATE and LaMATE3 expression in oocytes. Furthermore, citrate secretion was largely unaffected in P-deficient composite mutant plants with genome-edited or RNAi-silenced LaMATE in roots. Moderately lower concentrations of citrate and malate in the root tissue and consequently less organic acid anion secretion and lower malate in the xylem sap were identified. Interestingly, however, less genistein was consistently found in mutant exudates, opening the possibility that LaMATE is involved in isoflavonoid release.

Identification of ABC transporter G subfamily in white lupin and functional characterization of L.albABGC29 in phosphorus use

BACKGROUND

White lupin (Lupinus albus) is a leguminous crop with elite adaptive ability in phosphorus-deficient soil and used as a model plant for studying phosphorus (P) use. However, the genetic basis of its adaptation to low P (LP) remains unclear. ATPase binding cassette (ABC) transports G subfamily play a crucial role in the transportation of biological molecules across the membrane. To date, identification of this subfamily has been analyzed in some plants, but no systematic analysis of these transporters in phosphorus acquisition is available for white lupin.

RESULTS

This study identified 66 ABCG gene family members in the white lupin genome using comprehensive approaches. Phylogenetic analysis of white lupin ABCG transporters revealed six subclades based on their counterparts in Arabidopsis, displaying distinct gene structure and motif distribution in each cluster. Influences of the whole genome duplication on the evolution of L.albABCGs were investigated in detail. Segmental duplications appear to be the major driving force for the expansion of ABCGs in white lupin. Analysis of the Ka/Ks ratios indicated that the paralogs of the L.albABCG subfamily members principally underwent purifying selection. However, it was found that L.albABCG29 was a result of both tandem and segmental duplications. Overexpression of L.albABCG29 in white lupin hairy root enhanced P accumulation in cluster root under LP and improved plant growth. Histochemical GUS staining indicated that L.albABCG29 expression increased under LP in white lupin roots. Further, overexpression of L.albABCG29 in rice significantly improved P use under combined soil drying and LP by improving root growth associated with increased rhizosheath formation.

CONCLUSION

Through systematic and comprehensive genome-wide bioinformatics analysis, including conserved domain, gene structures, chromosomal distribution, phylogenetic relationships, and gene duplication analysis, the L.albABCG subfamily was identified in white lupin, and L.albABCG29 characterized in detail. In summary, our results provide deep insight into the characterization of the L.albABCG subfamily and the role of L.albABCG29 in improving P use.

Transcription Factor GmWRKY46 Enhanced Phosphate Starvation Tolerance and Root Development in Transgenic Plants

Phosphorus (P) is one of the essential macronutrients, whose deficiency limits the growth and development of plants. In this study, we investigated the possible role of GmWRKY46 in the phosphate (Pi) starvation stress tolerance of soybean. GmWRKY46 belonged to the group III subfamily of the WRKY transcription factor family, which was localized in the nucleus and had transcriptional activator activity. GmWRKY46 could be strongly induced by Pi starvation, especially in soybean roots. Overexpression of GmWRKY46 significantly enhanced tolerance to Pi starvation and lateral root development in transgenic Arabidopsis. RNA-seq analysis showed that overexpression of GmWRKY46 led to change in many genes related to energy metabolisms, stress responses, and plant hormone signal transduction in transgenic Arabidopsis. Among these differential expression genes, we found that overexpression of AtAED1 alone could enhance the tolerance of transgenic Arabidopsis to Pi starvation. Y1H and ChIP-qPCR analyses showed that GmWRKY46 could directly bind to the W-box motif of the AtAED1 promoter in vitro and in vivo. Furthermore, results from intact soybean composite plants with GmWRKY46 overexpression showed that GmWRKY46 was involved in hairy roots development and subsequently affected plant growth and Pi uptake. These results provide a basis for the molecular genetic breeding of soybean tolerant to Pi starvation.

Truncated hemoglobin 2 modulates phosphorus deficiency response by controlling of gene expression in nitric oxide-dependent pathway in Chlamydomonas reinhardtii

Truncated hemoglobin 2 is involved in fine-tuning of PSR1-regulated gene expression during phosphorus deprivation. Truncated hemoglobins form a large family found in all domains of life. However, a majority of physiological functions of these proteins remain to be elucidated. In the model alga Chlamydomonas reinhardtii, macro-nutritional deprivation is known to elevate truncated hemoglobin 2 (THB2). This study investigated the role of THB2 in the regulation of a subset of phosphorus (P) limitation-responsive genes in cells suffering from P-deficiency. Underexpression of THB2 in amiTHB2 strains resulted in downregulation of a suite of P deprivation-induced genes encoding proteins with different subcellular location and functions (e.g., PHOX, LHCSR3.1, LHCSR3.2, PTB2, and PTB5). Moreover, our results provided primary evidence that the soluble guanylate cyclase 12 gene (CYG12) is a component of the P deprivation regulation. Furthermore, the transcription of PSR1 gene for the most critical regulator in the acclimation process under P restriction was repressed by nitric oxide (NO). Collectively, the results indicated a tight regulatory link between the THB2-controlled NO levels and PSR1-dependent induction of several P deprivation responsive genes with various roles in cells during P-limitation.

The role of strigolactones in P deficiency induced transcriptional changes in tomato roots

BACKGROUND

Phosphorus (P) is an essential macronutrient for plant growth and development. Upon P shortage, plant responds with massive reprogramming of transcription, the Phosphate Starvation Response (PSR). In parallel, the production of strigolactones (SLs)-a class of plant hormones that regulates plant development and rhizosphere signaling molecules-increases. It is unclear, however, what the functional link is between these two processes. In this study, using tomato as a model, RNAseq was used to evaluate the time-resolved changes in gene expression in the roots upon P starvation and, using a tomato CAROTENOID CLEAVAGE DIOXYGENASES 8 (CCD8) RNAi line, what the role of SLs is in this.

RESULTS

Gene ontology (GO)-term enrichment and KEGG analysis of the genes regulated by P starvation and P replenishment revealed that metabolism is an important component of the P starvation response that is aimed at P homeostasis, with large changes occurring in glyco-and galactolipid and carbohydrate metabolism, biosynthesis of secondary metabolites, including terpenoids and polyketides, glycan biosynthesis and metabolism, and amino acid metabolism. In the CCD8 RNAi line about 96% of the PSR genes was less affected than in wild-type (WT) tomato. For example, phospholipid biosynthesis was suppressed by P starvation, while the degradation of phospholipids and biosynthesis of substitute lipids such as sulfolipids and galactolipids were induced by P starvation. Around two thirds of the corresponding transcriptional changes depend on the presence of SLs. Other biosynthesis pathways are also reprogrammed under P starvation, such as phenylpropanoid and carotenoid biosynthesis, pantothenate and CoA, lysine and alkaloids, and this also partially depends on SLs. Additionally, some plant hormone biosynthetic pathways were affected by P starvation and also here, SLs are required for many of the changes (more than two thirds for Gibberellins and around one third for Abscisic acid) in the gene expression.

CONCLUSIONS

Our analysis shows that SLs are not just the end product of the PSR in plants (the signals secreted by plants into the rhizosphere), but also play a major role in the regulation of the PSR (as plant hormone).

Comparative transcriptome analysis reveals changes in gene expression in sea cucumber (Holothuria leucospilota) in response to acute temperature stress

Ambient temperature is an important abiotic factor that influences growth performance and physiological functions in sea cucumbers. To understand the molecular responses of the sea cucumber Holothuria leucospilota to acute temperature stress, we performed a de novo transcriptome analysis of body wall tissue from H. leucospilota exposed to 2 hoursh of acute heat (35 ± 1 °C) and cold stress (15 ± 1 °C). A total of 99,015 unigenes were obtained after assembly of the sequenced reads. Compared with a control group maintained at 25.0 ± 1 °C, 1169 differentially expressed unigenes (DEGs) were identified after heat stress, 781 were up-regulated and 388 were down-regulated. After cold stress, 1464 DEGs were identified; 900 were up-regulated and 564 were down-regulated. The annotation of DEGs revealed that heat shock proteins play important roles in protecting H. leucospilota from high temperature stress. Furthermore, KEGG pathway enrichment analysis showed that the categories: "Ribosome" and "Protein processing in endoplasmic reticulum" were strongly affected by heat stress. These two pathways are associated with biosynthesis and processing of proteins, and refolding of misfolded proteins. The lipid metabolism pathways "Sphingolipid metabolism" and "Ether lipid metabolism", were affected by cold stress. The RNA-Seq results for eight selected DEGs were verified the expression by quantitative real-time PCR analysis. Our results will improve the understanding of the molecular response mechanisms of H. leucospilota to ambient temperature stress.

Quantitative Control of Early Flowering in White Lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a pulse annual plant cultivated from the tropics to temperate regions for its high-protein grain as well as a cover crop or green manure. Wild populations are typically late flowering and have high vernalization requirements. Nevertheless, some early flowering and thermoneutral accessions were found in the Mediterranean basin. Recently, quantitative trait loci (QTLs) explaining flowering time variance were identified in bi-parental population mapping, however, phenotypic and genotypic diversity in the world collection has not been addressed yet. In this study, a diverse set of white lupin accessions (n = 160) was phenotyped for time to flowering in a controlled environment and genotyped with PCR-based markers (n = 50) tagging major QTLs and selected homologs of photoperiod and vernalization pathway genes. This survey highlighted quantitative control of flowering time in white lupin, providing statistically significant associations for all major QTLs and numerous regulatory genes, including white lupin homologs of CONSTANS, FLOWERING LOCUS T, FY, MOTHER OF FT AND TFL1, PHYTOCHROME INTERACTING FACTOR 4, SKI-INTERACTING PROTEIN 1, and VERNALIZATION INDEPENDENCE 3. This revealed the complexity of flowering control in white lupin, dispersed among numerous loci localized on several chromosomes, provided economic justification for future genome-wide association studies or genomic selection rather than relying on simple marker-assisted selection.

Probing into the Effects of Grapevine Leafroll-Associated Viruses on the Physiology, Fruit Quality and Gene Expression of Grapes

Grapevine leafroll is one of the most widespread and highly destructive grapevine diseases that is responsible for great economic losses to the grape and wine industries throughout the world. Six distinct viruses have been implicated in this disease complex. They belong to three genera, all in the family Closteroviridae. For the sake of convenience, these viruses are named as grapevine leafroll-associated viruses (GLRaV-1, -2, -3, -4, -7, and -13). However, their etiological role in the disease has yet to be established. Furthermore, how infections with each GLRaV induce the characteristic disease symptoms remains unresolved. Here, we first provide a brief overview on each of these GLRaVs with a focus on genome structure, expression strategies and gene functions, where available. We then provide a review on the effects of GLRaV infection on the physiology, fruit quality, fruit chemical composition, and gene expression of grapevine based on the limited information so far reported in the literature. We outline key methodologies that have been used to study how GLRaV infections alter gene expression in the grapevine host at the transcriptomic level. Finally, we present a working model as an initial attempt to explain how infections with GLRaVs lead to the characteristic symptoms of grapevine leafroll disease: leaf discoloration and downward rolling. It is our hope that this review will serve as a starting point for grapevine virology and the related research community to tackle this vastly important and yet virtually uncharted territory in virus-host interactions involving woody and perennial fruit crops.

Nitrogen and Phosphorus Interplay in Lupin Root Nodules and Cluster Roots

Nitrogen (N) and phosphorus (P) are two major plant nutrients, and their deficiencies often limit plant growth and crop yield. The uptakes of N or P affect each other, and consequently, understanding N-P interactions is fundamental. Their signaling mechanisms have been studied mostly separately, and integrating N-P interactive regulation is becoming the aim of some recent works. Lupins are singular plants, as, under N and P deficiencies, they are capable to develop new organs, the N2-fixing symbiotic nodules, and some species can also transform their root architecture to form cluster roots, hundreds of short rootlets that alter their metabolism to induce a high-affinity P transport system and enhance synthesis and secretion of organic acids, flavonoids, proteases, acid phosphatases, and proton efflux. These modifications lead to mobilization in the soil of, otherwise unavailable, P. White lupin (Lupinus albus) represents a model plant to study cluster roots and for understanding plant acclimation to nutrient deficiency. It tolerates simultaneous P and N deficiencies and also enhances uptake of additional nutrients. Here, we present the structural and functional modifications that occur in conditions of P and N deficiencies and lead to the organogenesis and altered metabolism of nodules and cluster roots. Some known N and P signaling mechanisms include different factors, including phytohormones and miRNAs. The combination of the individual N and P mechanisms uncovers interactive regulation pathways that concur in nodules and cluster roots. L. albus interlinks N and P recycling processes both in the plant itself and in nature.

PHO1 family members transport phosphate from infected nodule cells to bacteroids in Medicago truncatula

Legumes play an important role in the soil nitrogen availability via symbiotic nitrogen fixation (SNF). Phosphate (Pi) deficiency severely impacts SNF because of the high Pi requirement of symbiosis. Whereas PHT1 transporters are involved in Pi uptake into nodules, it is unknown how Pi is transferred from the plant infected cells to nitrogen-fixing bacteroids. We hypothesized that Medicago truncatula genes homologous to Arabidopsis PHO1, encoding a vascular apoplastic Pi exporter, are involved in Pi transfer to bacteroids. Among the seven MtPHO1 genes present in M. truncatula, we found that two genes, namely MtPHO1.1 and MtPHO1.2, were broadly expressed across the various nodule zones in addition to the root vascular system. Expressions of MtPHO1.1 and MtPHO1.2 in Nicotiana benthamiana mediated specific Pi export. Plants with nodule-specific downregulation of both MtPHO1.1 and MtPHO1.2 were generated by RNA interference (RNAi) to examine their roles in nodule Pi homeostasis. Nodules of RNAi plants had lower Pi content and a three-fold reduction in SNF, resulting in reduced shoot growth. Whereas the rate of 33Pi uptake into nodules of RNAi plants was similar to control, transfer of 33Pi from nodule cells into bacteroids was reduced and bacteroids activated their Pi-deficiency response. Our results implicate plant MtPHO1 genes in bacteroid Pi homeostasis and SNF via the transfer of Pi from nodule infected cells to bacteroids.

Phosphate or nitrate imbalance induces stronger molecular responses than combined nutrient deprivation in roots and leaves of chickpea plants

The negative effects of phosphate (Pi) and/or nitrate (NO₃^- ) fertilizers on the environment have raised an urgent need to develop crop varieties with higher Pi and/or nitrogen use efficiencies for cultivation in low-fertility soils. Achieving this goal depends upon research that focuses on the identification of genes involved in plant responses to Pi and/or NO₃^- starvation. Although plant responses to individual deficiency in either Pi (-Pi/+NO₃^- ) or NO₃^- (+Pi/-NO₃^- ) have been separately studied, our understanding of plant responses to combined Pi and NO₃^- deficiency (-Pi/-NO₃^- ) is still very limited. Using RNA-sequencing approach, transcriptome changes in the roots and leaves of chickpea cultivated under -Pi/+NO₃^- , +Pi/-NO₃^- or -Pi/-NO₃^- conditions were investigated in a comparative manner. -Pi/-NO₃^- treatment displayed lesser effect on expression changes of genes related to Pi or NO₃^- transport, signalling networks, lipid remodelling, nitrogen and Pi scavenging/remobilization/recycling, carbon metabolism and hormone metabolism than -Pi/+NO₃^- or +Pi/-NO₃^- treatments. Therefore, the plant response to -Pi/-NO₃^- is not simply an additive result of plant responses to -Pi/+NO₃^- and +Pi/-NO₃^- treatments. Our results indicate that nutrient imbalance is a stronger stimulus for molecular reprogramming than an overall deficiency.

Profiling of MicroRNAs and Their Targets in Roots and Shoots Reveals a Potential MiRNA-Mediated Interaction Network in Response to Phosphate Deficiency in the Forestry Tree Betula luminifera

Inorganic phosphate (Pi) is often lacking in natural and agro-climatic environments, which impedes the growth of economically important woody species. Plants have developed strategies to cope with low Pi (LP) availability. MicroRNAs (miRNAs) play important roles in responses to abiotic stresses, including nutrition stress, by regulating target gene expression. However, the miRNA-mediated regulation of these adaptive responses and their underlying coordinating signals are still poorly understood in forestry trees such as Betula luminifera. Transcriptomic libraries, small RNA (sRNA) libraries, and a mixed degradome cDNA library of B. luminifera roots and shoots treated under LP and normal conditions (CK) were constructed and sequenced using next-generation deep sequencing. A comprehensive B. luminifera transcriptome derived from its roots and shoots was constructed, and a total of 76,899 unigenes were generated. Analysis of the transcriptome identified 8,095 and 5,584 differentially expressed genes in roots and shoots, respectively, under LP conditions. sRNA sequencing analyses indicated that 66 and 60 miRNAs were differentially expressed in roots and shoots, respectively, under LP conditions. A total of 109 and 112 miRNA-target pairs were further validated in the roots and shoots, respectively, using degradome sequencing. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of differential miRNA targets indicated that the "ascorbate and aldarate metabolism" pathway responded to LP, suggesting miRNA-target pairs might participating in the removing of reactive oxidative species under LP stress. Moreover, a putative network of miRNA-target interactions involved in responses to LP stress in B. luminifera is proposed. Taken together, these findings provide useful information to decipher miRNA functions and establish a framework for exploring P signaling networks regulated by miRNAs in B. luminifera and other woody plants. It may provide new insights into the genetic engineering of high use efficiency of Pi in forestry trees.

Gene Expression Responses to Sequential Nutrient Deficiency Stresses in Soybean

Throughout the growing season, crops experience a multitude of short periods of various abiotic stresses. These stress events have long-term impacts on plant performance and yield. It is imperative to improve our understanding of the genes and biological processes underlying plant stress tolerance to mitigate end of season yield loss. The majority of studies examining transcriptional changes induced by stress focus on single stress events. Few studies have been performed in model or crop species to examine transcriptional responses of plants exposed to repeated or sequential stress exposure, which better reflect field conditions. In this study, we examine the transcriptional profile of soybean plants exposed to iron deficiency stress followed by phosphate deficiency stress (-Fe-Pi). Comparing this response to previous studies, we identified a core suite of genes conserved across all repeated stress exposures (-Fe-Pi, -Fe-Fe, -Pi-Pi). Additionally, we determined transcriptional response to sequential stress exposure (-Fe-Pi) involves genes usually associated with reproduction, not stress responses. These findings highlight the plasticity of the plant transcriptome and the complexity of unraveling stress response pathways.

Ca14-3-3 Interacts With CaWRKY58 to Positively Modulate Pepper Response to Low-Phosphorus Starvation

Low-phosphorus stress (LPS) and pathogen attack are two important stresses frequently experienced by plants in their natural habitats, but how plant respond to them coordinately remains under-investigated. Here, we demonstrate that CaWRKY58, a known negative regulator of the pepper (Capsicum annuum) response to attack by Ralstonia solanacearum, is upregulated by LPS. Virus-induced gene silencing (VIGS) and overexpression of CaWRKY58 in Nicotiana benthamiana plants in combination with chromatin immunoprecipitation (ChIP) and electrophoretic mobility shift assays (EMSA) demonstrated that CaWRKY58 positively regulates the response of pepper to LPS by directly targeting and regulating genes related to phosphorus-deficiency tolerance, including PHOSPHATE STARVATION RESPONSE1 (PHR1). Yeast two-hybrid assays revealed that CaWRKY58 interacts with a 14-3-3 protein (Ca14-3-3); this interaction was confirmed by pull-down, bimolecular fluorescence complementation (BiFC), and microscale thermophoresis (MST) assays. The interaction between Ca14-3-3 and CaWRKY58 enhanced the activation of PHR1 expression by CaWRKY58, but did not affect the expression of the immunity-related genes CaNPR1 and CaDEF1, which are negatively regulated by CaWRKY58 in pepper upon Ralstonia solanacearum inoculation. Collectively, our data indicate that CaWRKY58 negatively regulates immunity against Ralstonia solanacearum, but positively regulates tolerance to LPS and that Ca14-3-3 transcriptionally activates CaWRKY58 in response to LPS.

Transcriptional, metabolic, physiological and developmental responses of switchgrass to phosphorus limitation

Knowing how switchgrass (Panicum virgatum L.) responds and adapts to phosphorus (P)-limitation will aid efforts to optimize P acquisition and use in this species for sustainable biomass production. This integrative study investigated the impacts of mild, moderate, and severe P-stress on genome transcription and whole-plant metabolism, physiology and development in switchgrass. P-limitation reduced overall plant growth, increased root/shoot ratio, increased root branching at moderate P-stress, and decreased root diameter with increased density and length of root hairs at severe P-stress. RNA-seq analysis revealed thousands of genes that were differentially expressed under moderate and severe P-stress in roots and/or shoots compared to P-replete plants, with many stress-induced genes involved in transcriptional and other forms of regulation, primary and secondary metabolism, transport, and other processes involved in P-acquisition and homeostasis. Amongst the latter were multiple miRNA399 genes and putative targets of these. Metabolite profiling showed that levels of most sugars and sugar alcohols decreased with increasing P stress, while organic and amino acids increased under mild and moderate P-stress in shoots and roots, although this trend reversed under severe P-stress, especially in shoots.

Development of PCR-based markers and whole-genome selection model for anthracnose resistance in white lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a high-protein grain legume crop, grown since ancient Greece and Rome. Despite long domestication history, its cultivation remains limited, partly because of susceptibility to anthracnose. Only some late-flowering, bitter, low-yielding landraces from Ethiopian mountains displayed resistance to this devastating disease. The resistance is controlled by various genes, thereby complicating the breeding efforts. The objective of this study was developing tools for molecular tracking of Ethiopian resistance genes based on genotyping-by-sequencing (GBS) data, envisaging linkage mapping and genomic selection approaches. Twenty GBS markers from two major quantitative trait loci (QTLs), antr04_1/antr05_1 and antr04_2/antr05_2, were converted to PCR-based markers using assigned transcriptome sequences. Newly developed markers improved mapping resolution around both anthracnose resistance loci, providing more precise QTL estimation. PCR-based screening of diversified domesticated and primitive germplasm revealed the high specificity of two markers for the antr04_1/antr05_1 locus (TP222136 and TP47110) and one for the antr04_2/antr05_2 locus (TP338761), highlighted by simple matching coefficients of 0.96 and 0.89, respectively. Moreover, a genomic selection approach based on GBS data of a recombinant inbred line mapping population was assessed, providing an average predictive ability of 0.56. These tools can be used for preselection of candidate white lupin germplasm for anthracnose resistance assays.

Fulvic acid increases forage legume growth inducing preferential up-regulation of nodulation and signalling-related genes

The use of potential biostimulants is of broad interest in plant science for improving yields. The application of a humic derivative called fulvic acid (FA) may improve forage crop production. FA is an uncharacterized mixture of chemicals and, although it has been reported to increase growth parameters in many species including legumes, its mode of action remains unclear. Previous studies of the action of FA have lacked appropriate controls, and few have included field trials. Here we report yield increases due to FA application in three European Medicago sativa cultivars, in studies which include the appropriate nutritional controls which hitherto have not been used. No significant growth stimulation was seen after FA treatment in grass species in this study at the treatment rate tested. Direct application to bacteria increased Rhizobium growth and, in M. sativa trials, root nodulation was stimulated. RNA transcriptional analysis of FA-treated plants revealed up-regulation of many important early nodulation signalling genes after only 3 d. Experiments in plate, glasshouse, and field environments showed yield increases, providing substantial evidence for the use of FA to benefit M. sativa forage production.

De novo transcriptome assembly, functional annotation, and expression profiling of rye (Secale cereale L.) hybrids inoculated with ergot (Claviceps purpurea)

Rye is used as food, feed, and for bioenergy production and remain an essential grain crop for cool temperate zones in marginal soils. Ergot is known to cause severe problems in cross-pollinated rye by contamination of harvested grains. The molecular response of the underlying mechanisms of this disease is still poorly understood due to the complex infection pattern. RNA sequencing can provide astonishing details about the transcriptional landscape, hence we employed a transcriptomic approach to identify genes in the underlying mechanism of ergot infection in rye. In this study, we generated de novo assemblies from twelve biological samples of two rye hybrids with identified contrasting phenotypic responses to ergot infection. The final transcriptome of ergot susceptible (DH372) and moderately ergot resistant (Helltop) hybrids contain 208,690 and 192,116 contigs, respectively. By applying the BUSCO pipeline, we confirmed that these transcriptome assemblies contain more than 90% of gene representation of the available orthologue groups at Virdiplantae odb10. We employed a de novo assembled and the draft reference genome of rye to count the differentially expressed genes (DEGs) between the two hybrids with and without inoculation. The gene expression comparisons revealed that 228 genes were linked to ergot infection in both hybrids. The genome ontology enrichment analysis of DEGs associated them with metabolic processes, hydrolase activity, pectinesterase activity, cell wall modification, pollen development and pollen wall assembly. In addition, gene set enrichment analysis of DEGs linked them to cell wall modification and pectinesterase activity. These results suggest that a combination of different pathways, particularly cell wall modification and pectinesterase activity contribute to the underlying mechanism that might lead to resistance against ergot in rye. Our results may pave the way to select genetic material to improve resistance against ergot through better understanding of the mechanism of ergot infection at molecular level. Furthermore, the sequence data and de novo assemblies are valuable as scientific resources for future studies in rye.

Genome-Wide Dissection of the CRF Gene Family in Brassica napus Indicates that BnaCRF8s Specifically Regulate Root Architecture and Phosphate Homeostasis against Phosphate Fluctuation in Plants

Phosphorus (P) is an essential macronutrient required for plant growth and development. The involvement of cytokinin response factors (CRFs) in phosphate (Pi) homeostasis and lateral root (LR) initiation in Arabidopsis has been revealed. However, little is known in oil crops. Here, we performed genome-wide dissection of the CRF family in Brassica napus to identify 44 members, which were evolutionally classified into 6 subgroups. Among them, four BnaCRF8 genes were strongly upregulated by P deprivation, and were selected to be further investigated. Time course qRT-PCR analyses showed that four BnaCRF8 genes were enhanced dramatically after 12 h of P stress. Analyses of the subcellular localization in tobacco leaves indicated that BnaA7.CRF8 and BnaC2.CRF8 were localized in the nucleus. The expression of BnaCRF8 genes had constant negative effects on primary root growth and LR initiation and growth, and it reduced Pi acquisition and plant growth in Arabidopsis. Moreover, the expression of Pi homeostasis-related genes was modulated in BnaA7.CRF8 overexpression plants. These results suggest that BnaCRF8 genes might negatively regulate root architecture and plant growth through transcriptional modification of Pi homeostasis-related components. Overall, this study suggests that upregulation of BnaCRF8 genes might be a smart adaptive strategy to cope with continuous Pi deficiency in the environment.

Metabolic and transcriptomic analysis of two Cucurbita moschata germplasms throughout fruit development

BACKGROUND

Pumpkins (Cucurbita moschata; Cucurbitaceae) are valued for their fruits and seeds and are rich in nutrients. Carotenoids and sugar contents, as main feature of pumpkin pulp, are used to determine the fruit quality.

RESULTS

Two pumpkin germplasms, CMO-X and CMO-E, were analyzed regarding the essential quality traits such as dry weight, soluble solids, organic acids, carotenoids and sugar contents. For the comparison of fruit development in these two germplasms, fruit transcriptome was analyzed at 5 different developmental stages from 0 d to 40 d in a time course manner. Putative pathways for carotenoids biosynthesis and sucrose metabolism were developed in C. moschata fruit and homologs were identified for each key gene involved in the pathways. Gene expression data was found consistent with the accumulation of metabolites across developmental stages and also between two germplasms. PSY, PDS, ZEP, CRTISO and SUS, SPS, HK, FK were found highly correlated with the accumulation of carotenoids and sucrose metabolites, respectively, at different growth stages of C. moschata as shown by whole transcriptomic analysis. The results of qRT-PCR analysis further confirmed the association of these genes.

CONCLUSION

Developmental regulation of the genes associated with the metabolite accumulation can be considered as an important factor for the determination of C. moschata fruit quality. This research will facilitate the investigation of metabolic profiles in other cultivars.

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.

Improvement of alfalfa forage quality and management through the down-regulation of MsFTa1

Alfalfa (Medicago sativa L.) is one of the most important forage crops worldwide. As a perennial, alfalfa is cut several times each year. Farmers face a dilemma: if cut earlier, forage nutritive value is much higher but regrowth is affected and the longevity of the stand is severely compromised. On the other hand, if alfalfa is cut later at full flower, stands persist longer and more biomass may be harvested, but the nutritive value diminishes. Alfalfa is a strict long-day plant. We reasoned that by manipulating the response to photoperiod, we could delay flowering to improve forage quality and widen each harvesting window, facilitating management. With this aim, we functionally characterized the FLOWERING LOCUS T family of genes, represented by five members: MsFTa1, MsFTa2, MsFTb1, MsFTb2 and MsFTc. The expression of MsFTa1 correlated with photoperiodic flowering and its down-regulation led to severe delayed flowering. Altogether, with late flowering, low expression of MsFTa1 led to changes in plant architecture resulting in increased leaf to stem biomass ratios and forage digestibility. By manipulating photoperiodic flowering, we were able to improve the quality of alfalfa forage and management, which may allow farmers to cut alfalfa of high nutritive value without compromising stand persistence.

Transcriptome sequencing of the apricot (Prunus armeniaca L.) and identification of differentially expressed genes involved in drought stress

Apricot (Prunus armeniaca L.) is an important fruit crop that is widely planted throughout the world. But drought affects both yield and quality of apricot. In order to study the effects of long-term drought on the molecular and physiological mechanisms of apricot, we used transcriptome sequencing and measured physiological indices. First, 322 million high-quality clean reads were obtained, and 74,892 unigenes were generated for the transcriptome. Among the assembled unigenes, 18,671 simple sequence repeats (SSRs) and 5581 differentially expressed genes (DEGs) were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs revealed that starch and sucrose metabolism, plant-pathogen interaction and plant hormone signal transduction pathways are enriched. Additionally, we used quantitative real-time PCR (qRT-PCR) to confirm the RNA-seq results with 11 drought-related DEGs. Second, through the physiological analysis of apricot leaves under constant drought stress, and the results show the internal microstructure of apricot leaves changed to withstand drought stress. At the same time, plants exposed to long-term drought stress showed higher degree of membrane damage, which reduced photosynthesis in the damaged leaves. Our findings enrich the genome resources for apricot and refine our understanding of the molecular and physiological mechanisms of drought response in this fruit crop, providing insights into drought adaptation of the apricot.

The genome evolution and low-phosphorus adaptation in white lupin

White lupin (Lupinus albus) is a legume crop that develops cluster roots and has high phosphorus (P)-use efficiency (PUE) in low-P soils. Here, we assemble the genome of white lupin and find that it has evolved from a whole-genome triplication (WGT) event. We then decipher its diploid ancestral genome and reconstruct the three sub-genomes. Based on the results, we further reveal the sub-genome dominance and the genic expression of the different sub-genomes varying in relation to their transposable element (TE) density. The PUE genes in white lupin have been expanded through WGT as well as tandem and dispersed duplications. Furthermore, we characterize four main pathways for high PUE, which include carbon fixation, cluster root formation, soil-P remobilization, and cellular-P reuse. Among these, auxin modulation may be important for cluster root formation through involvement of potential genes LaABCG36s and LaABCG37s. These findings provide insights into the genome evolution and low-P adaptation of white lupin.

De novo transcriptome characterisation of two auxin-related genes associated with plant growth habit in Astragalus adsurgens Pall

Astragalus adsurgens Pall., a perennial legume native to China, is commonly used as a forage crop. And it has great value for sustainable development of grasslands in arid and semi arid regions. However, to date, little is known regarding the A. adsurgens genome, and no studies have determined whether it would be possible to improve the germplasm of A. adsurgens through genetic modification. In this study, we used an RNA-seq protocol to generate a de novo transcriptome including 151,516 unigenes of A. adsurgens. We compared the transcriptomes of A. adsurgens having different growth habits (prostrate/erect) and identified 14,133 single nucleotide polymorphism sites (SNP) in 8,139 unigenes. Differential expression gene (DEG) analysis suggested that 10,982 unigenes were up-regulated in the prostrate plant relative to the erect plant, while 10,607 unigenes were down-regulated. Of the 21,589 DEG, Unigene72782_All (LAX4) and CL12494.Contig3_All (TIR1), an auxin transporter gene and an auxin transport inhibitor gene, respectively, were predicted to influence the growth habit of A. adsurgens, which were verified by qRT-PCR in these phenotypes. These results suggest that auxin transport was more active in the prostrate plant than in the erect plant, resulting in asymmetric distribution of auxin that affects the growth habit of A. adsurgens. Overall, this study may provide a basis for future research on key genes in A. adsurgens and may deepen our understanding of the molecular mechanisms regulating plant growth habit.

Legume genomics and transcriptomics: From classic breeding to modern technologies

Legumes are essential and play a significant role in maintaining food standards and augmenting physiochemical soil properties through the biological nitrogen fixation process. Biotic and abiotic factors are the main factors limiting legume production. Classical breeding methodologies have been explored extensively about the problem of truncated yield in legumes but have not succeeded at the desired rate. Conventional breeding improved legume genotypes but with more resources and time. Recently, the invention of next-generation sequencing (NGS) and high-throughput methods for genotyping have opened new avenues for research and developments in legume studies. During the last decade, genome sequencing for many legume crops documented. Sequencing and re-sequencing of important legume species have made structural variation and functional genomics conceivable. NGS and other molecular techniques such as the development of markers; genotyping; high density genetic linkage maps; quantitative trait loci (QTLs) identification, expressed sequence tags (ESTs), single nucleotide polymorphisms (SNPs); and transcription factors incorporated into existing breeding technologies have made possible the accurate and accelerated delivery of information for researchers. The application of genome sequencing, RNA sequencing (transcriptome sequencing), and DNA sequencing (re-sequencing) provide considerable insights for legume development and improvement programs.

Moreover, RNA-Seq helps to characterize genes, including differentially expressed genes, and can be applied for functional genomics studies, especially when there is limited information available for the studied genomes. Genome-based crop development studies and the availability of genomics data as well as decision-making gears look be specific for breeding programs. This review mainly presents an overview of the path from classical breeding to new emerging genomics tools, which will trigger and accelerate genomics-assisted breeding for recognition of novel genes for yield and quality characters for sustainable legume crop production.

Common and specific responses to iron and phosphorus deficiencies in roots of apple tree (Malus × domestica)

Iron and phosphorus are abundant elements in soils but poorly available for plant nutrition. The availability of these two nutrients represents a major constraint for fruit tree cultivation such as apple (Malus × domestica) leading very often to a decrease of fruit productivity and quality worsening. Aim of this study was to characterize common and specific features of plant response to Fe and P deficiencies by ionomic, transcriptomic and exudation profiling of apple roots. Under P deficiency, the root release of oxalate and flavonoids increased. Genes encoding for transcription factors and transporters involved in the synthesis and release of root exudates were upregulated by P-deficient roots, as well as those directly related to P acquisition. In Fe-deficiency, plants showed an over-accumulation of P, Zn, Cu and Mn and induced the transcription of those genes involved in the mechanisms for the release of Fe-chelating compounds and Fe mobilization inside the plants. The intriguing modulation in roots of some transcription factors, might indicate that, in this condition, Fe homeostasis is regulated by a FIT-independent pathway. In the present work common and specific features of apple response to Fe and P deficiency has been reported. In particular, data indicate similar modulation of a. 230 genes, suggesting the occurrence of a crosstalk between the two nutritional responses involving the transcriptional regulation, shikimate pathway, and the root release of exudates.

The first genetic map for yellow lupin enables genetic dissection of adaptation traits in an orphan grain legume crop

BACKGROUND

Yellow lupin (Lupinus luteus L.) is a promising grain legume for productive and sustainable crop rotations. It has the advantages of high tolerance to soil acidity and excellent seed quality, but its current yield potential is poor, especially in low rainfall environments. Key adaptation traits such as phenology and enhanced stress tolerance are often complex and controlled by several genes. Genomic-enabled technologies may help to improve our basic understanding of these traits and to provide selective markers in breeding. However, in yellow lupin there are very limited genomic resources to support research and no published information is available on the genetic control of adaptation traits.

RESULTS

We aimed to address these deficiencies by developing the first linkage map for yellow lupin and conducting quantitative trait locus (QTL) analysis of yield under well-watered (WW) and water-deficit (WT) conditions. Two next-generation sequencing marker approaches - genotyping-by-sequencing (GBS) and Diversity Array Technology (DArT) sequencing - were employed to genotype a recombinant inbred line (RIL) population developed from a bi-parental cross between wild and domesticated parents. A total of 2,458 filtered single nucleotide polymorphism (SNP) and presence / absence variation (PAV) markers were used to develop a genetic map comprising 40 linkage groups, the first reported for this species. A number of significant QTLs controlling total biomass and 100-seed weight under two water (WW and WD) regimes were found on linkage groups YL-03, YL-09 and YL-26 that together explained 9 and 28% of total phenotypic variability. QTLs associated with length of the reproductive phase and time to flower were found on YL-01, YL-21, YL-35 and YL-40 that together explained a total of 12 and 44% of total phenotypic variation.

CONCLUSION

These genomic resources and the QTL information offer significant potential for use in marker-assisted selection in yellow lupin.

Development of gene-based molecular markers tagging low alkaloid pauper locus in white lupin (Lupinus albus L.)

White lupin (Lupinus albus L.) is a legume grain crop cultivated since ancient Greece and Egypt. Modern white lupin cultivars are appreciated as a source of protein with positive nutraceutical impact. However, white lupins produce anti-nutritional compounds, quinolizidine alkaloids, which provide bitter taste and have a negative influence on human health. During domestication of this species, several recessive alleles at unlinked loci controlling low alkaloid content were selected. One of these loci, pauper, was exploited worldwide providing numerous low-alkaloid cultivars. However, molecular tracking of pauper has been hampered due to the lack of diagnostic markers. In the present study, the synteny-based approach was harnessed to target pauper locus. Single-nucleotide polymorphisms flanking pauper locus on white lupin linkage map as well as candidate gene sequences elucidated from the narrow-leafed lupin (L. angustifolius L.) chromosome segment syntenic to the pauper linkage group region were transformed to PCR-based molecular markers. These markers were analyzed both in the mapping population and world germplasm collection. From fourteen markers screened, eleven were localized at a distance below 1.5 cM from this locus, including five co-segregating with pauper. The linkage of these markers was confirmed by high LOD values (up to 58.4). Validation performed in the set of 127 bitter and 23 sweet accessions evidenced high applicability of one marker, LAGI01_35805_F1_R1, for pauper locus selection, highlighted by the low ratio of false-positive scores (2.5%). LAGI01_35805 represents a homolog of L. angustifolius acyltransferase-like (LaAT) gene which might hypothetically participate in the alkaloid biosynthesis process in lupins.

Physiological and transcriptomic data highlight common features between iron and phosphorus acquisition mechanisms in white lupin roots

In agricultural soil, the bioavailability of iron (Fe) and phosphorus (P) is often below the plant's requirement causing nutritional deficiency in crops. Under P-limiting conditions, white lupin (Lupinus albus L.) activates mechanisms that promote P solubility in the soil through morphological, physiological and molecular adaptations. Similar changes occur also in Fe-deficient white lupin roots; however, no information is available on the molecular bases of the response. In the present work, responses to Fe and P deficiency and their reciprocal interactions were studied. Transcriptomic analyses indicated that white lupin roots upregulated Fe-responsive genes ascribable to Strategy-I response, this behaviour was mainly evident in cluster roots. The upregulation of some components of Fe-acquisition mechanism occurred also in P-deficient cluster roots. Concerning P acquisition, some P-responsive genes (as phosphate transporters and transcription factors) were upregulated by P deficiency as well by Fe deficiency. These data indicate a strong cross-connection between the responses activated under Fe or P deficiency in white lupin. The activation of Fe- and P-acquisition mechanisms might play a crucial role to enhance the plant's capability to mobilize both nutrients in the rhizosphere, especially P from its associated metal cations.

The LaCEP1 peptide modulates cluster root morphology in Lupinus albus

White lupin cluster roots are specialized brush-like root structures that are formed in some species under phosphorus (P)-deficient conditions. They intensely secrete protons and organic acid anions for solubilization and acquisition of sparingly soluble phosphates. Phytohormones and sucrose modulate cluster root number, but the molecular mechanisms of cluster root formation have been elusive. Here, a novel peptide phytohormone was identified that affects cluster root development. It belongs to the C-TERMINALLY-ENCODED PEPTIDE (CEP) family. Members of that family arrest root growth and modulate branching in model species. LaCEP1 was highly expressed in the pre-emergence zone of clusters. Over-expression of the gene encoding the LaCEP1 propeptide resulted in moderate inhibition of cluster root formation. The primary and lateral root lengths of lupin were little affected by the overexpression, but LaCEP1 reduced cluster rootlet and root hair elongation. Addition of a 15-mer core peptide derived from LaCEP1 similarly altered root morphology and modified cluster activity, suggesting that a core sequence of the propeptide is functionally sufficient. Stable overexpression in Arabidopsis confirmed the LaCEP1 function in root growth inhibition across species. Taken together, the root inhibitory effects of the LaCEP1 phytohormone suggest a role as of a regulatory module involved in cluster root development in white lupin.

Structural basis for the Target DNA recognition and binding by the MYB domain of phosphate starvation response 1

The phosphate starvation response 1 (PHR1) protein has a central role in mediating the response to phosphate starvation in plants. PHR1 is composed of a number of domains including a MYB domain involved with DNA binding and a coiled-coil domain proposed to be involved with dimer formation. PHR1 binds to the promoter of phosphate starvation-induced genes to control the levels of phosphate required for nutrition. Previous studies have shown that both the MYB domain and the coiled-coil domain of PHR1 are required for binding the target DNA. Here, we describe the crystal structure of the PHR1 MYB domain and two structures of its complex with the PHR1-binding DNA sequence (P1BS). Structural and isothermal titration calorimetry has been carried out showing that the MYB domain of PHR1 alone is sufficient for target DNA recognition and binding. Two copies of the PHR1 MYB domain bind to the same major groove of the P1BS DNA with few direct interactions between the individual MYB domains. In addition, the PHR1 MYB-P1BS DNA complex structures reveal amino acid residues involved in DNA recognition and binding. Mutagenesis of these residues results in lost or impaired ability of PHR1 MYB to bind to its target DNA. The results presented reveal the structural basis for DNA recognition by the PHR1 MYB domain and demonstrate that two PHR1 MYB domains attach to their P1BS DNA targeting sequence. DATABASE: Coordinates and structure factors have been deposited in the Protein Data Bank under accession codes 6J4K (PHR1 MYB), 6J4R (PHR1 MYB-R-P1BS), 6J5B (MYB-CC-R2-P1BS).

The Use of Lupin as a Source of Protein in Animal Feeding: Genomic Tools and Breeding Approaches

Livestock production in the European Union EU is highly dependent on imported soybean, exposing the livestock farming system to risks related to the global trade of soybean. Lupin species could be a realistic sustainable alternative source of protein for animal feeding. Lupinus is a very diverse genus with many species. However, only four of them—namely, L. albus, L. angustifolius, L. luteus and L. mutabilis—are cultivated. Their use in livestock farming systems has many advantages in relation to economic and environmental impact. Generally, lupin grains are characterized by high protein content, while their oil content is relatively low but of high quality. On the other hand, the presence of quinolizidine alkaloids and their specific carbohydrate composition are the main antinutritional factors that prevent their use in animal feeding. This research is mainly related to L. albus and to L. angustifolius, and to a lesser extent, to L. lauteus and L. mutabilis. The breeding efforts are mostly focused on yield stabilization, resistance to biotic and abiotic stresses, biochemical structure associated with seed quality and late maturing. Progress is made in improving lupin with respect to the seed quality, as well as the tolerance to biotic and abiotic stress. It has to be noted that modern cultivars, mostly of L. albus and L. angustifolius, contain low levels of alkaloids. However, for future breeding efforts, the implementation of marker-assisted selection and the available genomic tools is of great importance.

Analysis of the common bean (Phaseolus vulgaris L.) transcriptome regarding efficiency of phosphorus use

Common bean is a highly important food in tropical regions, where most production occurs on small farms with limited use of technology and, consequently, greater vulnerability to abiotic stresses such as nutritional stress. Usually phosphorus (P) is the most limiting nutrient for crop growth in these regions. The aim of this study was to characterize the gene expression profiles of the genotypes of common bean IAC Imperador (P-responsive) and DOR 364 (P-unresponsive) under different P concentrations using RNA-seq transcriptome sequencing technology. Plants were grown hydroponically, with application of two P concentrations (4.00 mg L^-1 restrictive level and 8.00 mg L^-1 control level). Differential expression analyses, annotation, and functional classification were performed comparing genotypes within each P rate administered and comparing each genotype response to the different P levels. Considering differential expression analyses within genotypes, IAC Imperador exhibited 1538 up-regulated genes under P restriction and 1679 up-regulated genes in the control, while DOR 364 exhibited 13 up-regulated genes in the control and only 2 up-regulated genes under P restriction, strongly corroborating P-unresponsiveness of this genotype. Genes related to phosphorus restriction were identified among the differentially expressed genes, including transcription factors such as WRKY, ERF, and MYB families, phosphatase related genes such as pyrophosphatase, acid phosphatase, and purple acid phosphatase, and phosphate transporters. The enrichment test for the P restriction treatment showed 123 enriched gene ontologies (GO) for IAC Imperador, while DOR 364 enriched only 24. Also, the enriched GO correlated with P metabolism, compound metabolic processes containing phosphate, nucleoside phosphate binding, phosphorylation, and also response to stresses. Thus, this study proved to be informative to phosphorus limitation in common bean showing global changes at transcript level.

Legume Cytosolic and Plastid Acetyl-Coenzyme-A Carboxylase Genes Differ by Evolutionary Patterns and Selection Pressure Schemes Acting before and after Whole-Genome Duplications

Acetyl-coenzyme A carboxylase (ACCase, E.C.6.4.1.2) catalyzes acetyl-coenzyme A carboxylation to malonyl coenzyme A. Plants possess two distinct ACCases differing by cellular compartment and function. Plastid ACCase contributes to de novo fatty acid synthesis, whereas cytosolic enzyme to the synthesis of very long chain fatty acids, phytoalexins, flavonoids, and anthocyanins. The narrow leafed lupin (Lupinus angustifolius L.) represents legumes, a plant family which evolved by whole-genome duplications (WGDs). The study aimed on the contribution of these WGDs to the multiplication of ACCase genes and their further evolutionary patterns. The molecular approach involved bacterial artificial chromosome (BAC) library screening, fluorescent in situ hybridization, linkage mapping, and BAC sequencing. In silico analysis encompassed sequence annotation, comparative mapping, selection pressure calculation, phylogenetic inference, and gene expression profiling. Among sequenced legumes, the highest number of ACCase genes was identified in lupin and soybean. The most abundant plastid ACCase subunit genes were accB. ACCase genes in legumes evolved by WGDs, evidenced by shared synteny and Bayesian phylogenetic inference. Transcriptional activity of almost all copies was confirmed. Gene duplicates were conserved by strong purifying selection, however, positive selection occurred in Arachis (accB2) and Lupinus (accC) lineages, putatively predating the WGD event(s). Early duplicated accA and accB genes underwent transcriptional sub-functionalization.

Genome Wide Transcriptome Analysis Reveals Complex Regulatory Mechanisms Underlying Phosphate Homeostasis in Soybean Nodules

Phosphorus (P) deficiency is a major limitation for legume crop production. Although overall adaptations of plant roots to P deficiency have been extensively studied, only fragmentary information is available in regard to root nodule responses to P deficiency. In this study, genome wide transcriptome analysis was conducted using RNA-seq analysis in soybean nodules grown under P-sufficient (500 μM KH₂PO₄) and P-deficient (25 μM KH₂PO₄) conditions to investigate molecular mechanisms underlying soybean (Glycine max) nodule adaptation to phosphate (Pi) starvation. Phosphorus deficiency significantly decreased soybean nodule growth and nitrogenase activity. Nodule Pi concentrations declined by 49% in response to P deficiency, but this was well below the 87% and 88% decreases observed in shoots and roots, respectively. Nodule transcript profiling revealed that a total of 2055 genes exhibited differential expression patterns between Pi sufficient and deficient conditions. A set of (differentially expressed genes) DEGs appeared to be involved in maintaining Pi homeostasis in soybean nodules, including eight Pi transporters (PTs), eight genes coding proteins containing the SYG1/PHO81/XPR1 domain (SPXs), and 16 purple acid phosphatases (PAPs). The results suggest that a complex transcriptional regulatory network participates in soybean nodule adaption to Pi starvation, most notable a Pi signaling pathway, are involved in maintaining Pi homeostasis in nodules.

Differential response of hexaploid and tetraploid wheat to interactive effects of elevated [CO2] and low phosphorus

Hexaploid wheat is more responsive than tetraploid to the interactive effects of elevated [CO₂] and low P in terms of carboxylate efflux, enzyme activity and gene expression (TaPT1 and TaPAP). Availability of mineral nutrients to plants under changing climate has become a serious challenge to food security and economic development. An understanding of how elevated [CO₂] influences phosphorus (P) acquisition processes at the whole-plant level would be critical in selecting cultivars as well as to maintain optimum yield in limited-P conditions. Wheat (Triticum aestivum and T. durum) grown hydroponically with sufficient and low P concentration were exposed to elevated and ambient [CO₂]. Improved dry matter partitioning towards root resulted in increased root-to-shoot ratio, root length, volume, surface area, root hair length and density at elevated [CO₂] with low P. Interaction of low P and [CO₂] induced activity of enzymes (phosphoenolpyruvate carboxylase, malate dehydrogenase and citrate synthase) in root tissue resulting in twofold increase in carboxylates and acid phosphatase exudation. Physiological absorption capacity of roots showed that plants alter their uptake kinetics by increasing affinity (low K_m) in response to elevated [CO₂] under low P supply. Increased relative expression of genes, purple acid phosphatase (TaPAP) and high-affinity Pi transporter (TaPT1) in roots induced by elevated [CO₂] and low P supported our physiological observations. Hexaploid wheat (PBW-396) being more responsive to elevated [CO₂] at low P supply as compared to tetraploid (PDW-233) necessitates the ploidy effect to be explored further which might be advantageous under changing climate.