Integrating transcriptomic and metabolomic analysis to understand natural leaf senescence in sunflower

Comprehensive Proteome and Acetylome Analysis of Needle Senescence in Larix gmelinii

Leaf senescence is essential for the growth and development of deciduous trees in the next season. Larix gmelinii, a deciduous coniferous tree, exhibits its most distinctive feature by turning yellow in the autumn and eventually shedding its leaves, resulting in significant changes in its appearance during the fall. Lysine acetylation plays an important role in diverse cellular processes; however, limited knowledge is available regarding acetylations in the needle senescence of L. gmelinii. In this study, the proteomics and acetylated modification omics of two phenotypic leaves, yellow and green (senescent and non-senescent) needles, were analyzed before autumn defoliation. In total, 5022 proteins and 4469 unique acetylation sites in 2414 lysine acylated proteins were identified, and this resulted in the discovery of 1335 differentially expressed proteins (DEPs) and 605 differentially expressed acetylated proteins (DAPs) in yellow versus green needles. There are significant differences between the proteome and acetylome; only 269 proteins were found to be DEP and DAP, of which 136 proteins were consistently expressed in both the DEP and DAP, 91 proteins were upregulated, and 45 proteins were down-regulated. The DEPs participate in the metabolism of starch and sucrose, while the DAPs are involved in glycolysis and the tricarboxylic acid cycle. Among them, DEPs underwent significant changes in glycolysis and citric acid cycling. Most of the enzymes involved in glycolysis and the citrate cycle were acetylated. DAPs were down-regulated in glycolysis and up-regulated in the citrate cycle. In all, the results of this study reveal the important role of lysine acetylation in the senescence of L. gmelinii needles and provide a new perspective for understanding the molecular mechanism of leaf senescence and tree seasonal growth.

Integrative Physiological, Transcriptome, and Proteome Analyses Provide Insights into the Photosynthetic Changes in Maize in a Maize-Peanut Intercropping System

Intercropping is a traditional and sustainable planting method that can make rational use of natural resources such as light, temperature, fertilizer, water, and CO2. Due to its efficient resource utilization, intercropping, in particular, maize and legume intercropping, is widespread around the world. However, the molecular details of these pathways remain largely unknown. In this study, physiological, transcriptome, and proteome analyses were compared between maize monocropping and maize-peanut intercropping. The results show that an intercropping system enhanced the ability of carbon fixation and carboxylation of maize leaves. Apparent quantum yield (AQY), the light-saturated net photosynthetic rate (LSPn), the light saturation point (LSP), and the light compensation point (LCP) were increased by 11.6%, 9.4%, 8.9%, and 32.1% in the intercropping system, respectively; carboxylation efficiency (CE), the CO2 saturation point (Cisat), the Rubisco maximum carboxylation rate (Vcmax), the maximum electron transfer rate (Jmax), and the triose phosphate utilization rate (TPU) were increased by 28.5%, 7.3%, 18.7%, 29.2%, and 17.0%, respectively; meanwhile, the CO2 compensation point (Γ) decreased by 22.6%. Moreover, the transcriptome analysis confirmed the presence of 588 differentially expressed genes (DEGs), and the numbers of up-regulated and down-regulated genes were 383 and 205, respectively. The DEGs were primarily concerned with ribosomes, plant hormone signal transduction, and photosynthesis. Furthermore, 549 differentially expressed proteins (DEPs) were identified in the maize leaves in both the maize monocropping and maize-peanut intercropping systems. Bioinformatics analysis revealed that 186 DEPs were related to 37 specific KEGG pathways in each of the two treatment groups. Based on the physiological, transcriptome, and proteome analyses, it was demonstrated that the photosynthetic characteristics in maize leaves can be improved by maize-peanut intercropping. This may be related to PS I, PS II, cytochrome b6f complex, ATP synthase, and photosynthetic CO2 fixation, which is caused by the improved CO2 carboxylation efficiency. Our results provide a more in-depth understanding of the high yield and high-efficiency mechanism in maize and peanut intercropping.

Integrated transcriptomics and miRNAomics provide insights into the complex multi-tiered regulatory networks associated with coleoptile senescence in rice

Coleoptile is the small conical, short-lived, sheath-like organ that safeguards the first leaf and shoot apex in cereals. It is also the first leaf-like organ to senesce that provides nutrition to the developing shoot and is, therefore, believed to play a crucial role in seedling establishment in rice and other grasses. Though histochemical studies have helped in understanding the pattern of cell death in senescing rice coleoptiles, genome-wide expression changes during coleoptile senescence have not yet been explored. With an aim to investigate the gene regulation underlying the coleoptile senescence (CS), we performed a combinatorial whole genome expression analysis by sequencing transcriptome and miRNAome of senescing coleoptiles. Transcriptome analysis revealed extensive reprogramming of 3439 genes belonging to several categories, the most prominent of which encoded for transporters, transcription factors (TFs), signaling components, cell wall organization enzymes, redox homeostasis, stress response and hormone metabolism. Small RNA sequencing identified 41 known and 21 novel miRNAs that were differentially expressed during CS. Comparison of gene expression and miRNA profiles generated for CS with publicly available leaf senescence (LS) datasets revealed that the two aging programs are remarkably distinct at molecular level in rice. Integration of expression data of transcriptome and miRNAome identified high confidence 140 miRNA-mRNA pairs forming 42 modules, thereby demonstrating multi-tiered regulation of CS. The present study has generated a comprehensive resource of the molecular networks that enrich our understanding of the fundamental pathways regulating coleoptile senescence in rice.

Function of Protein Kinases in Leaf Senescence of Plants

Leaf senescence is an evolutionarily acquired process and it is critical for plant fitness. During senescence, macromolecules and nutrients are disassembled and relocated to actively growing organs. Plant leaf senescence process can be triggered by developmental cues and environmental factors, proper regulation of this process is essential to improve crop yield. Protein kinases are enzymes that modify their substrates activities by changing the conformation, stability, and localization of those proteins, to play a crucial role in the leaf senescence process. Impressive progress has been made in understanding the role of different protein kinases in leaf senescence recently. This review focuses on the recent progresses in plant leaf senescence-related kinases. We summarize the current understanding of the function of kinases on senescence signal perception and transduction, to help us better understand how the orderly senescence degeneration process is regulated by kinases, and how the kinase functions in the intricate integration of environmental signals and leaf age information.

Genome-wide and comparative phylogenetic analysis of senescence-associated NAC transcription factors in sunflower (Helianthus annuus)

BACKGROUND

Leaf senescence delay impacts positively in grain yield by maintaining the photosynthetic area during the reproductive stage and during grain filling. Therefore a comprehensive understanding of the gene families associated with leaf senescence is essential. NAC transcription factors (TF) form a large plant-specific gene family involved in regulating development, senescence, and responses to biotic and abiotic stresses. The main goal of this work was to identify sunflower NAC TF (HaNAC) and their association with senescence, studying their orthologous to understand possible functional relationships between genes of different species.

RESULTS

To clarify the orthologous relationships, we used an in-depth comparative study of four divergent taxa, in dicots and monocots, with completely sequenced genomes (Arabidopsis thaliana, Vitis vinifera, Musa acuminata and Oryza sativa). These orthologous groups provide a curated resource for large scale protein sequence annotation of NAC TF. From the 151 HaNAC genes detected in the latest version of the sunflower genome, 50 genes were associated with senescence traits. These genes showed significant differential expression in two contrasting lines according to an RNAseq assay. An assessment of overexpressing the Arabidopsis line for HaNAC001 (a gene of the same orthologous group of Arabidopsis thaliana ORE1) revealed that this line displayed a significantly higher number of senescent leaves and a pronounced change in development rate.

CONCLUSIONS

This finding suggests HaNAC001 as an interesting candidate to explore the molecular regulation of senescence in sunflower.

Genetic Improvement in Sunflower Breeding—Integrated Omics Approach

Foresight in climate change and the challenges ahead requires a systematic approach to sunflower breeding that will encompass all available technologies. There is a great scarcity of desirable genetic variation, which is in fact undiscovered because it has not been sufficiently researched as detection and designing favorable genetic variation largely depends on thorough genome sequencing through broad and deep resequencing. Basic exploration of genomes is insufficient to find insight about important physiological and molecular mechanisms unique to crops. That is why integrating information from genomics, epigenomics, transcriptomics, proteomics, metabolomics and phenomics enables a comprehensive understanding of the molecular mechanisms in the background of architecture of many important quantitative traits. Omics technologies offer novel possibilities for deciphering the complex pathways and molecular profiling through the level of systems biology and can provide important answers that can be utilized for more efficient breeding of sunflower. In this review, we present omics profiling approaches in order to address their possibilities and usefulness as a potential breeding tools in sunflower genetic improvement.

Identification of a New Giant Emrbryo Allele, and Integrated Transcriptomics and Metabolomics Analysis of Giant Embryo Development in Rice

Rice embryos are rich in high-quality protein, lipid, vitamins and minerals, representing the most important nutritional part of brown rice. However, the molecular mechanism of rice embryo development is poorly understood. In this study, two rice cultivars with contrasting embryo size (the giant embryo cultivar Dapeimi and the normal embryo cultivar 187R) were used to explore excellent genes controlling embryo size, and the developed near-isogenic lines (NILs) (NIL-D, which has the giant embryo phenotype, and its matching line, NIL-X) were used to explore transcript and metabolic properties in the earlier maturation stage of giant embryo development under natural conditions. The map-based cloning results demonstrated that Dapeimi is a novel allelic mutant of the rice GIANT EMBRYO (GE) gene, and the functional mutation site is a single cytosine deletion in the exon1. A total of 285 differentially accumulated metabolites (DAMs) and 677 differentially expressed genes (DEGs) were identified between NIL-D and NIL-X. The analysis of DAMs indicated that plants lacking GE mainly promoted energy metabolism, amino acid metabolism, and lipid metabolism pathways in the rice embryo. Pearson correlation coefficient showed that 300 pairs of gene-metabolites were highly correlated. Among them, OsZS_02G0528500 and OsZS_12G0013700 were considered to be key genes regulating L-Aspartic acid and L-Tryptophan content during rice giant embryo development, which are promising to be good candidate genes to improve rice nutrition. By analyzing rice embryo development through a combination of strategies, this research contributes to a greater understanding of the molecular mechanism of rice embryo development, and provides a theoretical foundation for breeding high-nutrition varieties.

Comparative Transcriptome Analysis of Salt Stress-Induced Leaf Senescence in Medicago truncatula

Leaves are the most critical portion of forage crops such as alfalfa (Medicago sativa). Leaf senescence caused by environmental stresses significantly impacts the biomass and quality of forages. To understand the molecular mechanisms and identify the key regulator of the salt stress-induced leaf senescence process, we conducted a simple and effective salt stress-induced leaf senescence assay in Medicago truncatula, which was followed by RNA-Seq analysis coupled with physiological and biochemical characterization. By comparing the observed expression data with that derived from dark-induced leaf senescence at different time points, we identified 3,001, 3,787, and 4,419 senescence-associated genes (SAGs) for salt stress-induced leaf senescence on day 2, 4, and 6, respectively. There were 1546 SAGs shared by dark and salt stress treatment across the three time points. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that the 1546 SAGs were mainly related to protein and amino acids metabolism, photosynthesis, chlorophyll metabolism, and hormone signaling during leaf senescence. Strikingly, many different transcription factors (TFs) families out of the 1546 SAGs, including NAC, bHLH, MYB, and ERF, were associated with salt stress-induced leaf senescence processes. Using the transient expression system in Nicotiana benthamiana, we verified that three functional NAC TF genes from the 1546 SAGs were related to leaf senescence. These results clarify SAGs under salt stress in M. truncatula and provide new insights and additional genetic resources for further forage crop breeding.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression patterns during somatic embryogenesis in Dimocarpus longan Lour

The NAM, ATAF1/2, and CUC2 form a huge plant-specific gene family of NAC TFs that are involved in the growth, development, and regulation of biotic and abiotic stress responses. Although the draft genome of longan (Dimocarpus longan Lour.) has been published, however the comprehensive data regarding the functions, evolution, and expression patterns of the NAC family are still unavailable. In this study, a comprehensive analysis of the NAC transcription factor family in longan was performed, and a total of 114 NAC genes were found. We investigated the NAC gene family exploring the phylogeny, domain conservation, intron/exon, motifs, cis-regulatory elements, protein-protein interaction, and expression profiles of RNA-seq samples in different tissues and early somatic embryogenesis of longan. Phylogenetic analysis showed that the genes with similar gene structure and motif distribution were clustered in the same group. Cis-element identification indicates the possible role of NAC genes in biological and physiological processes. Protein-protein interaction identified the DlNACs homologous with Arabidopsis proteins. We further investigated the expression pattern of DlNAC genes in different tissues (pulp, stem, large fruit, young fruit, and flower) during somatic embryogenesis at embryogenic callus (EC), incomplete compact pro-embryogenic cultures (ICpEC), and globular embryos (GE) stages. The qRT-PCR results showed that the DlNAC genes were expressed higher at EC and GE stage compared with ICpEC stage. In conclusion, our results provide insight into the evolution, diversity, and characterization of NAC genes in the longan and provide a base for understanding their biological roles and molecular mechanisms in plants.

Exploring sunflower responses to Sclerotinia head rot at early stages of infection using RNA-seq analysis

Sclerotinia head rot (SHR), caused by the necrotrophic fungus Sclerotinia sclerotiorum, is one of the most devastating sunflower crop diseases. Despite its worldwide occurrence, the genetic determinants of plant resistance are still largely unknown. Here, we investigated the Sclerotinia-sunflower pathosystem by analysing temporal changes in gene expression in one susceptible and two tolerant inbred lines (IL) inoculated with the pathogen under field conditions. Differential expression analysis showed little overlapping among ILs, suggesting genotype-specific control of cell defense responses possibly related to differences in disease resistance strategies. Functional enrichment assessments yielded a similar pattern. However, all three ILs altered the expression of genes involved in the cellular redox state and cell wall remodeling, in agreement with current knowledge about the initiation of plant immune responses. Remarkably, the over-representation of long non-coding RNAs (lncRNA) was another common feature among ILs. Our findings highlight the diversity of transcriptional responses to SHR within sunflower breeding lines and provide evidence of lncRNAs playing a significant role at early stages of defense.

Source-sink modifications affect leaf senescence and grain mass in wheat as revealed by proteomic analysis

BACKGROUND

The grain yield of cereals is determined by the synergistic interaction between source activity and sink capacity. However, source-sink interactions are far from being fully understood. Therefore, a field experiment was performed in wheat to investigate the responses of flag leaves and grains to sink/source manipulations.

RESULTS

Half-degraining delayed but partial defoliation enhanced leaf senescence. Sink/source manipulations influenced the content of reactive oxygen species in the flag leaf and the concentration of phytohormones, including cytokinins, indoleacetic 3-acid and jasmonic acid, in the flag leaves (LDef) and grains (GDef) in defoliated plants and flag leaves (LDG) and grain (GDG) in de-grained plants. Isobaric tag for relative and absolute quantitation (iTRAQ)-based quantitative proteomic analysis indicated that at 16 days after manipulation, a total of 97 and 59 differentially expressed proteins (DEPs) from various functional categories were observed in the LDG and LDef groups, respectively, compared with the control, and 115 and 121 DEPs were observed in the GDG and GDef groups, respectively. The gene ontology annotation terms of the DEPs mainly included carbon fixation, hydrogen peroxide catabolic process, chloroplast and cytoplasm, oxidoreductase activity and glutamate synthase activity in the flag leaves of manipulated plants and organonitrogen compound metabolic process, cytoplasm, vacuolar membrane, CoA carboxylase activity, starch synthase activity and nutrient reservoir activity in the grains of manipulated plants. KEGG pathway enrichment analysis revealed that photosynthesis, carbon, nitrogen and pyruvate metabolism and glycolysis/gluconeogenesis were the processes most affected by sink/source manipulations. Sink/source manipulations affected the activities of amylase and proteinases and, ultimately, changed the mass per grain.

CONCLUSIONS

Manipulations to change the sink/source ratio affect hormone levels; hydrolytic enzyme activities; metabolism of carbon, nitrogen and other main compounds; stress resistance; and leaf senescence and thus influence grain mass.

Integrating physiological and metabolites analysis to identify ethylene involvement in petal senescence in Tulipa gesneriana

Flower senescence is classified into ethylene-dependent and ethylene-independent manners and determines the flower longevity which is valuable for ornamental plants. However, the manner of petal senescence in tulip is still less defined. In this study, we characterized the physiological indexes in the process of petal senescence, as well as metabolic and ethylene responses in tulip cultivar 'American Dream', and further identified the role of ethylene biosynthesis genes TgACS by transgenic and transient assays. Primary metabolites profiling revealed that sugars, amino acids and organic acids preferentially accumulated in senescent petals. Additionally, senescence-associated genes were identified and significantly up-regulated, coupled with increased ROS contents, rapid water loss and accelerated cell membrane breakdown. Moreover, ethylene production was stimulated as evidenced by increasing in ACS activity and ethylene biosynthesis-related genes expression. Exogenous treatment of cutting flowers with 1-MCP or ethephon resulted in delayed or enhanced petal senescence, respectively. Transient down-regulation of TgACS by VIGS assay in tulip petals delayed senescence, while over-expressed TgACS1 in tobacco promoted leaf senescence. Taken together, this study provides evidences to certify ethylene roles and TgACS functions during flower senescence in tulip.

Transgenic expression of rice MYB102 (OsMYB102) delays leaf senescence and decreases abiotic stress tolerance in Arabidopsis thaliana

MYB-type transcription factors (TFs) play important roles in plant growth and development, and in the rapid responses to unfavorable environmental conditions. We recently reported the isolation and characterization of a rice (Oryza sativa) MYB TF, OsMYB102, which is involved in the regulation of leaf senescence by downregulating abscisic acid (ABA) biosynthesis and the downstream signaling response. Based on the similarities of their sequences and expression patterns, OsMYB102 appears to be a homolog of the Arabidopsis thaliana AtMYB44 TF. Since AtMYB44 is a key regulator of leaf senescence and abiotic stress responses, it is important to examine whether AtMYB44 homologs in other plants also act similarly. Here, we generated transgenic Arabidopsis plants expressing OsMYB102 (OsMYB102-OX). The OsMYB102-OX plants showed a delayed senescence phenotype during dark incubation and were more susceptible to salt and drought stresses, considerably similar to Arabidopsis plants overexpressing AtMYB44. Real-time quantitative PCR (RT-qPCR) revealed that, in addition to known senescence-associated genes, genes encoding the ABA catabolic enzymes AtCYP707A3 and AtCYP707A4 were also significantly upregulated in OsMYB102-OX, leading to a significant decrease in ABA accumulation. Furthermore, protoplast transient expression and chromatin immunoprecipitation assays revealed that OsMYB102 directly activated AtCYP707A3 expression. Based on our findings, it is probable that the regulatory functions of AtMYB44 homologs in plants are highly conserved and they have vital roles in leaf senescence and the abiotic stress responses.

Dissecting the Regulatory Network of Leaf Premature Senescence in Maize (Zea mays L.) Using Transcriptome Analysis of ZmELS5 Mutant

Leaf premature senescence largely determines maize (Zea mays L.) grain yield and quality. A natural recessive premature-senescence mutant was selected from the breeding population, and near-isogenic lines were constructed using Jing24 as the recurrent parent. In the near-isogenic lines, the dominant homozygous material was wild-type (WT), and the recessive material of early leaf senescence was the premature-senescence-type ZmELS5. To identify major genes and regulatory mechanisms involved in leaf senescence, a transcriptome analysis of the ZmELS5 and WT near-isogenic lines (NILs) was performed. A total of 8796 differentially expressed transcripts were identified between ZmELS5 and WT, including 3811 up-regulated and 4985 down-regulated transcripts. By combining gene ontology, Kyoto Encyclopedia of Genes and Genomes, gene set, and transcription factor enrichment analyses, key differentially expressed genes were screened. The senescence regulatory network was predicted based on these key differentially expressed genes, which indicated that the senescence process is mainly regulated by bHLH, WRKY, and AP2/EREBP family transcription factors, leading to the accumulations of jasmonic acid and ethylene. This causes stress responses and reductions in the chlorophyll a/b-binding protein activity level. Then, decreased ATP synthase activity leads to increased photosystem II photodamage, ultimately leading to leaf senescence.

Asparagine Synthesis During Tobacco Leaf Curing

Senescence is a genetically controlled mechanism that modifies leaf chemistry. This involves significant changes in the accumulation of carbon- and nitrogen-containing compounds, including asparagine through the activity of asparagine synthetases. These enzymes are required for nitrogen re-assimilation and remobilization in plants; however, their mechanisms are not fully understood. Here, we report how leaf curing—a senescence-induced process that allows tobacco leaves to dry out—modifies the asparagine metabolism. We show that leaf curing strongly alters the concentration of the four main amino acids, asparagine, glutamine, aspartate, and glutamate. We demonstrate that detached tobacco leaf or stalk curing has a different impact on the expression of asparagine synthetase genes and accumulation of asparagine. Additionally, we characterize the main asparagine synthetases involved in the production of asparagine during curing. The expression of ASN1 and ASN5 genes is upregulated during curing. The ASN1-RNAi and ASN5-RNAi tobacco plant lines display significant alterations in the accumulation of asparagine, glutamine, and aspartate relative to wild-type plants. These results support the idea that ASN1 and ASN5 are key regulators of asparagine metabolism during leaf curing.

Transcriptomic analysis reveals flavonoid biosynthesis of Syringa oblata Lindl. in response to different light intensity

BACKGROUND

Hazy weather significantly increase air pollution and affect light intensity which may also affect medicinal plants growth. Syringa oblata Lindl. (S. oblata), an effective anti-biofilm medicinal plants, is also vulnerable to changes in plant photoperiods and other abiotic stress responses. Rutin, one of the flavonoids, is the main bioactive ingredient in S. oblata that inhibits Streptococcus suis biofilm formation. Thus, the present study aims to explore the biosynthesis and molecular basis of flavonoids in S. oblata in response to different light intensity.

RESULTS

In this study, it was shown that compared with natural (Z0) and 25% ~ 35% (Z2) light intensities, the rutin content of S. oblata under 50% ~ 60% (Z1) light intensity increased significantly. In addition, an integrated analysis of metabolome and transcriptome was performed using light intensity stress conditions from two kinds of light intensities which S. oblata was subjected to: Z0 and Z1. The results revealed that differential metabolites and genes were mainly related to the flavonoid biosynthetic pathway. We found out that 13 putative structural genes and a transcription factor bHLH were significantly up-regulated in Z1. Among them, integration analysis showed that 3 putative structural genes including 4CL1, CYP73A and CYP75B1 significantly up-regulated the rutin biosynthesis, suggesting that these putative genes may be involved in regulating the flavonoid biosynthetic pathway, thereby making them key target genes in the whole metabolic process.

CONCLUSIONS

The present study provided helpful information to search for the novel putative genes that are potential targets for S. oblata in response to light intensity.

Exploring gene networks in two sunflower lines with contrasting leaf senescence phenotype using a system biology approach

BACKGROUND

Leaf senescence is a complex process, controlled by multiple genetic and environmental variables. In sunflower, leaf senescence is triggered abruptly following anthesis thereby limiting the capacity of plants to keep their green leaf area during grain filling, which subsequently has a strong impact on crop yield. Recently, we performed a selection of contrasting sunflower inbred lines for the progress of leaf senescence through a physiological, cytological and molecular approach. Here we present a large scale transcriptomic analysis using RNA-seq and its integration with metabolic profiles for two contrasting sunflower inbred lines, R453 and B481-6 (early and delayed senescence respectively), with the aim of identifying metabolic pathways associated to leaf senescence.

RESULTS

Gene expression profiles revealed a higher number of differentially expressed genes, as well as, higher expression levels in R453, providing evidence for early activation of the senescence program in this line. Metabolic pathways associated with sugars and nutrient recycling were differentially regulated between the lines. Additionally, we identified transcription factors acting as hubs in the co-expression networks; some previously reported as senescence-associated genes in model species but many are novel candidate genes.

CONCLUSIONS

Understanding the onset and the progress of the senescence process in crops and the identification of these new candidate genes will likely prove highly useful for different management strategies to mitigate the impact of senescence on crop yield. Functional characterization of candidate genes will help to develop molecular tools for biotechnological applications in breeding crop yield.

Transcription Factors Associated with Leaf Senescence in Crops

Leaf senescence is a complex mechanism controlled by multiple genetic and environmental variables. Different crops present a delay in leaf senescence with an important impact on grain yield trough the maintenance of the photosynthetic leaf area during the reproductive stage. Additionally, because of the temporal gap between the onset and phenotypic detection of the senescence process, candidate genes are key tools to enable the early detection of this process. In this sense and given the importance of some transcription factors as hub genes in senescence pathways, we present a comprehensive review on senescence-associated transcription factors, in model plant species and in agronomic relevant crops. This review will contribute to the knowledge of leaf senescence process in crops, thus providing a valuable tool to assist molecular crop breeding.

Maize expressing the sunflower transcription factor HaHB11 has improved productivity in controlled and field conditions

HaHB11 is a sunflower transcription factor from the homeodomain-leucine zipper I family. Transgenic Arabidopsis plants expressing HaHB11 had larger rosettes and improved seed yield. In this work maize plants from hybrid HiII were transformed with 35S:HaHB11, ZmUBI:HaHB11 and ProHaHB11:HaHB11 and then backcrossed to B73 to obtain a more homozygous inbred phenotype. Transgene expression levels were stable at least during three generations. Greenhouse-grown HaHB11 transgenic lines had larger leaf area and delayed senescence than controls, together with increased total biomass (up to 25%) and seed yield (up to 28%). Field trials conducted with T2 and T4 generations indicated that enhanced leaf area (up to 18%), stem diameter (up to 28%) and total biomass (up to 40%) as well as delayed leaf senescence were maintained among transgenic individuals when upscaling from pots in the greenhouse to communal plants in the field. The T4 field-grown transgenic generation had increased light interception and radiation use efficiency as well as seed yield (43-47% for events driven by the 35S promoter). Results suggest that HaHB11 is a promising tool for crop improvement because differential traits observed in the Arabidopsis model plant were preserved in a crop like maize independently of growth conditions and backcross level.

Wheat drought tolerance in the field is predicted by amino acid responses to glasshouse-imposed drought

Water limits crop productivity, so selecting for a minimal yield gap in drier environments is critical to mitigate against climate change and land-use pressure. We investigated the responses of relative water content (RWC), stomatal conductance, chlorophyll content, and metabolites in flag leaves of commercial wheat (Triticum aestivum L.) cultivars to three drought treatments in the glasshouse and in field environments. We observed strong genetic associations between glasshouse-based RWC, metabolites, and yield gap-based drought tolerance (YDT; the ratio of yield in water-limited versus well-watered conditions) across 18 field environments spanning sites and seasons. Critically, RWC response to glasshouse drought was strongly associated with both YDT (r2=0.85, P<8E-6) and RWC under field drought (r2=0.77, P<0.05). Moreover, multiple regression analyses revealed that 98% of genetic YDT variance was explained by drought responses of four metabolites: serine, asparagine, methionine, and lysine (R2=0.98; P<0.01). Fitted coefficients suggested that, for given levels of serine and asparagine, stronger methionine and lysine accumulation was associated with higher YDT. Collectively, our results demonstrate that high-throughput, targeted metabolic phenotyping of glasshouse-grown plants may be an effective tool for selection of wheat cultivars with high field-derived YDT.

Metabolite signatures of grasspea suspension-cultured cells illustrate the complexity of dehydration response

This represents the first report deciphering the dehydration response of suspension-cultured cells of a crop species, highlighting unique and shared pathways, and adaptive mechanisms via profiling of 330 metabolites. Grasspea, being a hardy legume, is an ideal model system to study stress tolerance mechanisms in plants. In this study, we investigated the dehydration-responsive metabolome in grasspea suspension-cultured cells (SCCs) to identify the unique and shared metabolites crucial in imparting dehydration tolerance. To reveal the dehydration-induced metabolite signatures, SCCs of grasspea were exposed to 10% PEG, followed by metabolomic profiling. Chromatographic separation by HPLC coupled with MRM-MS led to the identification of 330 metabolites, designated dehydration-responsive metabolites (DRMs), which belonged to 28 varied functional classes. The metabolome was found to be constituted by carboxylic acids (17%), amino acids (13.5%), flavonoids (10.9%) and plant growth regulators (10%), among others. Pathway enrichment analysis revealed predominance of metabolites involved in phytohormone biosynthesis, secondary metabolism and osmotic adjustment. Exogenous application of DRMs, arbutin and acetylcholine, displayed improved physiological status in stress-resilient grasspea as well as hypersensitive pea, while administration of lauric acid imparted detrimental effects. This represents the first report on stress-induced metabolomic landscape of a crop species via a suspension culture system, which would provide new insights into the molecular mechanism of stress responses and adaptation in crop species.

De novo transcriptome assembly and co-expression network analysis of Cynanchum thesioides: Identification of genes involved in resistance to drought stress

Cynanchum thesioides are upright, xerophytic shrubs that are widely distributed in arid and semi-arid areas of China, North Korea, Mongolia and Siberia. To date, little is known about the molecular mechanisms of drought resistance in C. thesioides. To better understand drought resistance, we used transcriptome analysis and Illumina sequencing technology on C. thesioides, to identify drought-responsive genes. Using de novo assembly 55,268 unigenes were identified from 207.58 Gb of clean data. Amongst these, 36,265 were annotated with gene descriptions, conserved domains, gene ontology terms and metabolic pathways. The sequencing results showed that genes that were differentially expressed (DEGs) under drought stress were enriched in pathways such as carbon metabolism, starch and sucrose metabolism, amino acid biosynthesis, phenylpropanoid biosynthesis and plant hormone signal transduction. Moreover, many functional genes were up-regulated under severe drought stress to enhance tolerance. Weighted gene co-expression network analysis showed that there were key hub genes related to drought stress. Hundreds of candidate genes were identified under severe drought stress, including transcriptional factors such as MYB, G2-like, ERF, C2H2, NAC, NF-X1, GRF, HD-ZIP, HB-other, HSF, C3H, GRAS, WRKY, bHLH and Trihelix. These data are a valuable resource for further investigation into the molecular mechanism for drought stress in C. thesioides and will facilitate exploration of drought resistance genes.

Rice transcription factor OsMYB102 delays leaf senescence by down-regulating abscisic acid accumulation and signaling

MYB-type transcription factors (TFs) play important roles in plant growth and development, and in the responses to several abiotic stresses. In rice (Oryza sativa), the roles of MYB-related TFs in leaf senescence are not well documented. Here, we examined rice MYB TF gene OsMYB102 and found that an OsMYB102 T-DNA activation-tagged line (termed osmyb102-D), which constitutively expresses OsMYB102 under the control of four tandem repeats of the 35S promoter, and OsMYB102-overexpressing transgenic lines (35S:OsMYB102 and 35S:GFP-OsMYB102) maintain green leaves much longer than the wild-type under natural, dark-induced, and abscisic acid (ABA)-induced senescence conditions. Moreover, an osmyb102 knockout mutant showed an accelerated senescence phenotype under dark-induced and ABA-induced leaf senescence conditions. Microarray analysis showed that a variety of senescence-associated genes (SAGs) were down-regulated in the osmyb102-D line. Further studies demonstrated that overexpression of OsMYB102 controls the expression of SAGs, including genes associated with ABA degradation and ABA signaling (OsABF4, OsNAP, and OsCYP707A6), under dark-induced senescence conditions. OsMYB102 inhibits ABA accumulation by directly activating the transcription of OsCYP707A6, which encodes the ABA catabolic enzyme ABSCISIC ACID 8'-HYDROXYLASE. OsMYB102 also indirectly represses ABA-responsive genes, such as OsABF4 and OsNAP. Collectively, these results demonstrate that OsMYB102 plays a critical role in leaf senescence by down-regulating ABA accumulation and ABA signaling responses.

Opposite fates of the purine metabolite allantoin under water and nitrogen limitations in bread wheat

Degradation of nitrogen-rich purines is tightly and oppositely regulated under drought and low nitrogen supply in bread wheat. Allantoin is a key target metabolite for improving nitrogen homeostasis under stress. The metabolite allantoin is an intermediate of the catabolism of purines (components of nucleotides) and is known for its housekeeping role in nitrogen (N) recycling and also for its function in N transport and storage in nodulated legumes. Allantoin was also shown to differentially accumulate upon abiotic stress in a range of plant species but little is known about its role in cereals. To address this, purine catabolic pathway genes were identified in hexaploid bread wheat and their chromosomal location was experimentally validated. A comparative study of two Australian bread wheat genotypes revealed a highly significant increase of allantoin (up to 29-fold) under drought. In contrast, allantoin significantly decreased (up to 22-fold) in response to N deficiency. The observed changes were accompanied by transcriptional adjustment of key purine catabolic genes, suggesting that the recycling of purine-derived N is tightly regulated under stress. We propose opposite fates of allantoin in plants under stress: the accumulation of allantoin under drought circumvents its degradation to ammonium (NH₄⁺) thereby preventing N losses. On the other hand, under N deficiency, increasing the NH₄⁺ liberated via allantoin catabolism contributes towards the maintenance of N homeostasis.

Effect of post-silking drought stress on the expression profiles of genes involved in carbon and nitrogen metabolism during leaf senescence in maize (Zea mays L.)

Drought stress during reproductive growth stages greatly affects the growth and productivity of maize plants. To better understand the metabolic regulation during post-silking drought (PD) stress, an RNA sequencing (RNA-Seq) analysis was performed at the late stage of leaf senescence in maize. Physiological measurements showed that PD stress reduced both leaf carbon and nitrogen levels. A total of 4013 differentially expressed genes (DEGs) were found based on RNA-Seq analysis, 115 of which were identified to be involved in photosynthesis and in the metabolism of sucrose, starch, and amino acids. Among these DEGs, 14 genes involved in photosynthesis were down-regulated. The genes coding for sucrose and pectin synthesis were up-regulated under PD stress. The two genes of asparagine synthetase (ZmAS3 and ZmAS4), which are responsible for nitrogen remobilization in leaves, were also significantly induced by the drought treatment. The expression profiles of these genes involved in carbon and nitrogen metabolism suggests their regulatory roles during drought-induced leaf senescence.

Integration of multi-omics techniques and physiological phenotyping within a holistic phenomics approach to study senescence in model and crop plants

The study of senescence in plants is complicated by diverse levels of temporal and spatial dynamics as well as the impact of external biotic and abiotic factors and crop plant management. Whereas the molecular mechanisms involved in developmentally regulated leaf senescence are very well understood, in particular in the annual model plant species Arabidopsis, senescence of other organs such as the flower, fruit, and root is much less studied as well as senescence in perennials such as trees. This review addresses the need for the integration of multi-omics techniques and physiological phenotyping into holistic phenomics approaches to dissect the complex phenomenon of senescence. That became feasible through major advances in the establishment of various, complementary 'omics' technologies. Such an interdisciplinary approach will also need to consider knowledge from the animal field, in particular in relation to novel regulators such as small, non-coding RNAs, epigenetic control and telomere length. Such a characterization of phenotypes via the acquisition of high-dimensional datasets within a systems biology approach will allow us to systematically characterize the various programmes governing senescence beyond leaf senescence in Arabidopsis and to elucidate the underlying molecular processes. Such a multi-omics approach is expected to also spur the application of results from model plants to agriculture and their verification for sustainable and environmentally friendly improvement of crop plant stress resilience and productivity and contribute to improvements based on postharvest physiology for the food industry and the benefit of its customers.

Sunflower Hybrid Breeding: From Markers to Genomic Selection

In sunflower, molecular markers for simple traits as, e.g., fertility restoration, high oleic acid content, herbicide tolerance or resistances to Plasmopara halstedii, Puccinia helianthi, or Orobanche cumana have been successfully used in marker-assisted breeding programs for years. However, agronomically important complex quantitative traits like yield, heterosis, drought tolerance, oil content or selection for disease resistance, e.g., against Sclerotinia sclerotiorum have been challenging and will require genome-wide approaches. Plant genetic resources for sunflower are being collected and conserved worldwide that represent valuable resources to study complex traits. Sunflower association panels provide the basis for genome-wide association studies, overcoming disadvantages of biparental populations. Advances in technologies and the availability of the sunflower genome sequence made novel approaches on the whole genome level possible. Genotype-by-sequencing, and whole genome sequencing based on next generation sequencing technologies facilitated the production of large amounts of SNP markers for high density maps as well as SNP arrays and allowed genome-wide association studies and genomic selection in sunflower. Genome wide or candidate gene based association studies have been performed for traits like branching, flowering time, resistance to Sclerotinia head and stalk rot. First steps in genomic selection with regard to hybrid performance and hybrid oil content have shown that genomic selection can successfully address complex quantitative traits in sunflower and will help to speed up sunflower breeding programs in the future. To make sunflower more competitive toward other oil crops higher levels of resistance against pathogens and better yield performance are required. In addition, optimizing plant architecture toward a more complex growth type for higher plant densities has the potential to considerably increase yields per hectare. Integrative approaches combining omic technologies (genomics, transcriptomics, proteomics, metabolomics and phenomics) using bioinformatic tools will facilitate the identification of target genes and markers for complex traits and will give a better insight into the mechanisms behind the traits.

Initiation, Progression, and Genetic Manipulation of Leaf Senescence

As a representative form of plant senescence, leaf senescence has received the most attention during the last two decades. In this chapter we summarize the initiation of leaf senescence by various internal and external signals, the progression of senescence including switches in gene expression, as well as changes at the biochemical and cellular levels during leaf senescence. Impacts of leaf senescence in agriculture and genetic approaches that have been used in manipulating leaf senescence of crop plants are discussed.

Intergrative metabolomic and transcriptomic analyses unveil nutrient remobilization events in leaf senescence of tobacco

Leaf senescence in plants is a coordinated process that involves remobilization of nutrients from senescing leaves to sink tissues. The molecular events associated with nutrient remobilization are however not well understood. In this study the tobacco system with a source-sink relationship between different leaf positions was used in analyzing the spatiotemporal changes of 76 metabolites from leaves at 3 different stalk positions and 8 developmental stages. The metabolomic data was then compared with RNA-seq data from the same samples to analyze the activities of the metabolic pathways that are important for nutrient remobilization. Integrative analyses on metabolites accumulation and expression changes of enzyme-encoding genes in corresponding metabolic pathways indicated a significant up-regulation of the tricarboxylic acid cycle and related metabolism of sugars, amino acids and fatty acids, suggesting the importance of energy metabolism during leaf senescence. Other changes of the metabolism during tobacco leaf senescence include increased activities of the GS/GOGAT cycle which is responsible for nitrogen recycling, and increased accumulation of nicotine. The results also suggested that a number of compounds seemed to be transported from senescing leaves at lower positions to sink leaves at upper positions. Some of these metabolites could play a role in nutrient remobilization.

Integration of transcriptomic and metabolic data reveals hub transcription factors involved in drought stress response in sunflower (Helianthus annuus L.)

By integration of transcriptional and metabolic profiles we identified pathways and hubs transcription factors regulated during drought conditions in sunflower, useful for applications in molecular and/or biotechnological breeding. Drought is one of the most important environmental stresses that effects crop productivity in many agricultural regions. Sunflower is tolerant to drought conditions but the mechanisms involved in this tolerance remain unclear at the molecular level. The aim of this study was to characterize and integrate transcriptional and metabolic pathways related to drought stress in sunflower plants, by using a system biology approach. Our results showed a delay in plant senescence with an increase in the expression level of photosynthesis related genes as well as higher levels of sugars, osmoprotectant amino acids and ionic nutrients under drought conditions. In addition, we identified transcription factors that were upregulated during drought conditions and that may act as hubs in the transcriptional network. Many of these transcription factors belong to families implicated in the drought response in model species. The integration of transcriptomic and metabolomic data in this study, together with physiological measurements, has improved our understanding of the biological responses during droughts and contributes to elucidate the molecular mechanisms involved under this environmental condition. These findings will provide useful biotechnological tools to improve stress tolerance while maintaining crop yield under restricted water availability.

Metabolic Investigation of Phelipanche aegyptiaca Reveals Significant Changes during Developmental Stages and in Its Different Organs

Phelipanche aegyptiaca Pers. is a root holoparasitic plant considered to be among the most destructive agricultural weeds worldwide. In order to gain more knowledge about the metabolic profile of the parasite during its developmental stages, we carried out primary metabolic and lipid profiling using GC-MS analysis. In addition, the levels of amino acids that incorporate into proteins, total protein in the albumin fraction, nitrogen, reduced sugars, and phenols were determined. For the assays, the whole plants from the four developmental stages-tubercle, pre-emergent shoot, post-emergent shoot, and mature flowering plants-were taken. Thirty-five metabolites out of 66 differed significantly between the various developmental stages. The results have shown that the first three developmental stages were distinguished in their profiles, but the latter two did not differ from the mature stage. Yet, 46% of the metabolites detected did not change significantly during the developmental stages. This is unlike other studies of non-parasitic plants showing that their metabolic levels tend to alter significantly during development. This implies that the parasite can control the levels of these metabolites. We further studied the metabolic nature of five organs (adventitious roots, lower and upper shoot, floral buds, and flowers) in mature plants. Similar to non-parasitic plants, the parasite exhibited significant differences between the vegetative and reproductive organs. Compared to other organs, floral buds had higher levels of free amino acids and total nitrogen, whereas flowers accumulated higher levels of simple sugars such as sucrose, and the putative precursors for nectar synthesis, color, and volatiles. This suggests that the reproductive organs have the ability to accumulate metabolites that are required for the production of seeds and as a source of energy for the reproductive processes. The data contribute to our knowledge about the metabolic behavior of parasites that rely on their host for its basic nutrients.

Excessive nitrogen application dampens antioxidant capacity and grain filling in wheat as revealed by metabolic and physiological analyses

In this study, field-grown wheat (Triticum aestivum L.) was treated with normal (Nn) and excessive (Ne) levels of fertilizer N. Results showed that Ne depressed the activity of superoxide dismutase and peroxidase and increased the accumulation of reactive oxygen species (ROS) and malondialdehyde. The normalized difference vegetation index (NDVI) was higher under Ne at anthesis and medium milk but similar at the early dough stage and significantly lower at the hard dough stage than that under Nn. The metabolomics analysis of the leaf responses to Ne during grain filling showed 99 metabolites that were different between Ne and Nn treatments, including phenolic and flavonoid compounds, amino acids, organic acids and lipids, which are primarily involved in ROS scavenging, N metabolism, heat stress adaptation and disease resistance. Organic carbon (C) and total N contents were affected by the Ne treatment, with lower C/N ratios developing after medium milk. Ultimately, grain yields decreased with Ne. Based on these data, compared with the normal N fertilizer treatment, we concluded that excessive N application decreased the ability to scavenge ROS, increased lipid peroxidation and caused significant metabolic changes disturbing N metabolism, secondary metabolism and lipid metabolism, which led to reduced grain filling in wheat.

A sunflower WRKY transcription factor stimulates the mobilization of seed-stored reserves during germination and post-germination growth

The sunflower transcription factor HaWRKY10 stimulates reserves mobilization in Arabidopsis. Gene expression and enzymes activity assays indicated that lipolysis and gluconeogenesis were increased. Microarray results suggested a parallelism in sunflower. Germinating oilseeds converts stored lipids into sugars, and thereafter in metabolic energy that is used in seedling growth and establishment. During germination, the induced lipolysis linked to the glyoxylate pathway and gluconeogenesis produces sucrose, which is then transported to the embryo and driven through catabolic routes. Herein, we report that the sunflower transcription factor HaWRKY10 regulates carbon partitioning by reducing carbohydrate catabolism and increasing lipolysis and gluconeogenesis. HaWRKY10 was regulated by abscisic acid and gibberellins in the embryo leaves 48 h after seed imbibition and highly expressed during sunflower seed germination and seedling growth, concomitantly with lipid mobilization. Sunflower leaf disks overexpressing HaWRKY10 showed repressed expression of genes related to sucrose cleavage and glycolysis compared with controls. Moreover, HaWRKY10 constitutive expression in Arabidopsis seeds produced higher decrease in lipid reserves, whereas starch and sucrose were more preserved compared with wild type. Gene transcripts abundance and enzyme activities involved in stored lipid mobilization and gluconeogenesis increased more in transgenic than in wild type seeds 36 h after imbibition, whereas the negative regulator of lipid mobilization, ABI4, was repressed. Altogether, the results point out a functional parallelism between tissues and plant species, and reveal HaWRKY10 as a positive regulator of storage reserve mobilization in sunflower.

Toward Systems Understanding of Leaf Senescence: An Integrated Multi-Omics Perspective on Leaf Senescence Research

Leaf senescence is a complex but tightly regulated developmental process involving a coordinated sequence of multiple molecular events, which ultimately leads to death of the leaf. Efforts to understand the mechanistic principles underlying leaf senescence have been largely made by transcriptomic, proteomic, and metabolomic studies over the past decade. This review focuses on recent milestones in leaf senescence research obtained using multi-omics technologies, as well as future endeavors toward systems understanding of leaf senescence processes. In particular, we discuss recent advances in understanding molecular events during leaf senescence through genome-wide transcriptome analyses in Arabidopsis. We also describe comparative transcriptome analyses used to unveil the commonality and diversity of regulatory mechanisms governing leaf senescence in the plant kingdom. Finally, we provide current illustrations of epigenomic, proteomic, and metabolomic landscapes of leaf senescence. We envisage that integration of multi-omics leaf senescence data will enable us to address unresolved questions regarding leaf senescence, including determining the molecular principles that coordinate concurrent and ordered changes in biological events during leaf senescence.

Network and biosignature analysis for the integration of transcriptomic and metabolomic data to characterize leaf senescence process in sunflower

Background

In recent years, high throughput technologies have led to an increase of datasets from omics disciplines allowing the understanding of the complex regulatory networks associated with biological processes. Leaf senescence is a complex mechanism controlled by multiple genetic and environmental variables, which has a strong impact on crop yield. Transcription factors (TFs) are key proteins in the regulation of gene expression, regulating different signaling pathways; their function is crucial for triggering and/or regulating different aspects of the leaf senescence process. The study of TF interactions and their integration with metabolic profiles under different developmental conditions, especially for a non-model organism such as sunflower, will open new insights into the details of gene regulation of leaf senescence.

Results

Weighted Gene Correlation Network Analysis (WGCNA) and BioSignature Discoverer (BioSD, Gnosis Data Analysis, Heraklion, Greece) were used to integrate transcriptomic and metabolomic data. WGCNA allowed the detection of 10 metabolites and 13 TFs whereas BioSD allowed the detection of 1 metabolite and 6 TFs as potential biomarkers. The comparative analysis demonstrated that three transcription factors were detected through both methodologies, highlighting them as potentially robust biomarkers associated with leaf senescence in sunflower.

Conclusions

The complementary use of network and BioSignature Discoverer analysis of transcriptomic and metabolomic data provided a useful tool for identifying candidate genes and metabolites which may have a role during the triggering and development of the leaf senescence process. The WGCNA tool allowed us to design and test a hypothetical network in order to infer relationships across selected transcription factor and metabolite candidate biomarkers involved in leaf senescence, whereas BioSignature Discoverer selected transcripts and metabolites which discriminate between different ages of sunflower plants. The methodology presented here would help to elucidate and predict novel networks and potential biomarkers of leaf senescence in sunflower.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1045-2) contains supplementary material, which is available to authorized users.