Standardization of data processing and statistical analysis in comparative plant proteomics experiment

Protein research in millets: current status and way forward

MAIN CONCLUSION

Millets' protein studies are lagging behind those of major cereals. Current status and future insights into the investigation of millet proteins are discussed. Millets are important small-seeded cereals majorly grown and consumed by people in Asia and Africa and are considered crops of future food security. Although millets possess excellent climate resilience and nutrient supplementation properties, their research advancements have been lagging behind major cereals. Although considerable genomic resources have been developed in recent years, research on millet proteins and proteomes is currently limited, highlighting a need for further investigation in this area. This review provides the current status of protein research in millets and provides insights to understand protein responses for climate resilience and nutrient supplementation in millets. The reference proteome data is available for sorghum, foxtail millet, and proso millet to date; other millets, such as pearl millet, finger millet, barnyard millet, kodo millet, tef, and browntop millet, do not have any reference proteome data. Many studies were reported on stress-responsive protein identification in foxtail millet, with most studies on the identification of proteins under drought-stress conditions. Pearl millet has a few reports on protein identification under drought and saline stress. Finger millet is the only other millet to have a report on stress-responsive (drought) protein identification in the leaf. For protein localization studies, foxtail millet has a few reports. Sorghum has the highest number of 40 experimentally proven crystal structures, and other millets have fewer or no experimentally proven structures. Further proteomics studies will help dissect the specific proteins involved in climate resilience and nutrient supplementation and aid in breeding better crops to conserve food security.

Identification of Proteases and Protease Inhibitors in Seeds of the Recalcitrant Forest Tree Species Quercus ilex

Proteases and protease inhibitors have been identified in the recalcitrant species Quercus ilex using in silico and wet methods, with focus on those present in seeds during germination. In silico analyses showed that the Q. ilex transcriptome database contained 2,240 and 97 transcripts annotated as proteases and protease inhibitors, respectively. They belonged to the different families according to MEROPS, being the serine and metallo ones the most represented. The data were compared with those previously reported for other Quercus species, including Q. suber, Q. lobata, and Q. robur. Changes in proteases and protease inhibitors alongside seed germination in cotyledon and embryo axis tissues were assessed using proteomics and in vitro and in gel activity assays. Shotgun (LC-MSMS) analysis of embryo axes and cotyledons in nonviable (NV), mature (T1) and germinated (T3) seeds allowed the identification of 177 proteases and 12 protease inhibitors, mostly represented by serine and metallo types. Total protease activity, as determined by in vitro assays using azocasein as substrate, was higher in cotyledons than in embryo axes. There were not differences in activity among cotyledon samples, while embryo axis peaked at germinated T4 stage. Gel assays revealed the presence of protease activities in at least 10 resolved bands, in the Mr range of 60-260 kDa, being some of them common to cotyledons and embryo axes in either nonviable, mature, and germinated seeds. Bands showing quantitative or qualitative changes upon germination were observed in embryo axes but not in cotyledons at Mr values of 60-140 kDa. Proteomics shotgun analysis of the 10 bands with protease activity supported the results obtained in the overall proteome analysis, with 227 proteases and 3 protease inhibitors identified mostly represented by the serine, cysteine, and metallo families. The combined use of shotgun proteomics and protease activity measurements allowed the identification of tissue-specific (e.g., cysteine protease inhibitors in embryo axes of mature acorns) and stage-specific proteins (e.g., those associated with mobilization of storage proteins accumulated in T3 stage). Those proteins showing differences between nonviable and viable seeds could be related to viability, and those variables between mature and germinated could be associated with the germination process. These differences are observed mostly in embryo axes but not in cotyledons. Among them, those implicated in mobilization of reserve proteins, such as the cathepsin H cysteine protease and Clp proteases, and also the large number of subunits of the CNS and 26S proteasome complex differentially identified in embryos of the several stages suggests that protein degradation via CNS/26S plays a major role early in germination. Conversely, aspartic proteases such as nepenthesins were exclusively identified in NV seeds, so their presence could be used as indicator of nonviability.

Proteomics and plant biology: contributions to date and a look towards the next decade

INTRODUCTION

This review presents the view of the author, that is opinionable and even speculative, on the field of proteomics, its application to plant biology knowledge, and translation to biotechnology. Written in a more academic than scientific style, it is based on past original and review articles by the author´s group, and those published by leading scientists in the last two years.

AREAS COVERED

Starting with a general definition and references to historical milestones, it covers sections devoted to the different platforms employed, the plant biology discourse in the protein language, challenges and future prospects, ending with the author opinion.

EXPERT OPINION

In 25 years, five proteomics platform generations have appeared. We are now moving from proteomics to Systems Biology. While feasible with model organisms, proteomics of orphan species remains challenging. Proteomics, even in its simplest approach, sheds light on plant biological processes, central dogma, and molecular bases of phenotypes of interest, and it can be translated to areas such as food traceability and allergen detection. Proteomics should be validated and optimized to each experimental system, objectives, and hypothesis. It has limitations, artifacts, and biases. We should not blindly accept proteomics data and just create a list of proteins, networks, and avoid speculative biological interpretations. From the hundred to thousand proteins identified and quantified, it is important to obtain a focus and validate some of them, otherwise it is merely. We are starting to have the protein pieces, so let, from now, build the proteomics and biological puzzle.

Protein Interaction Networks: Functional and Statistical Approaches

The evolution of next-generation sequencing and high-throughput technologies has created new opportunities and challenges in data science. Currently, a classic proteomics analysis can be complemented by going a step beyond the individual analysis of the proteome by using integrative approaches. These integrations can be focused either on inferring relationships among proteins themselves, with other molecular levels, phenotype, or even environmental data, giving the researcher new tools to extract and determine the most relevant information in biological terms. Furthermore, it is also important the employ of visualization methods that allow a correct and deep interpretation of data.To carry out these analyses, several bioinformatics and biostatistical tools are required. In this chapter, different workflows that enable the creation of interaction networks are proposed. Resulting networks reduce the complexity of original datasets, depicting complex statistical relationships (through PLS analysis and variants), functional networks (STRING, shinyGO), and a combination of both approaches. Recently developed methods for integrating different omics levels, such as coinertial analyses or DIABLO, are also described. Finally, the use of Cytoscape or Gephi was described for the representation and mining of the different networks.This approach constitutes a new way of acquiring a deeper knowledge of the function of proteins, such as the search for specific connections of each group to identify differentially connected modules, which may reflect involved protein complexes and key pathways.

Proteometabolomic characterization of apical bud maturation in Pinus pinaster

Bud maturation is a physiological process that implies a set of morphophysiological changes that lead to the transition of growth patterns from young to mature. This transition defines tree growth and architecture, and in consequence traits such as biomass production and wood quality. In Pinus pinaster Aiton, a conifer of great timber value, bud maturation is closely related to polycyclism (multiple growth periods per year). This process causes a lack of apical dominance, and consequently increased branching that reduces its timber quality and value. However, despite its importance, little is known about bud maturation. In this work, proteomics and metabolomics were employed to study apical and basal sections of young and mature buds in P. pinaster. Proteins and metabolites in samples were described and quantified using (n)UPLC-LTQ-Orbitrap. The datasets were analyzed employing an integrative statistical approach, which allowed the determination of the interactions between proteins and metabolites and the different bud sections and ages. Specific dynamics of proteins and metabolites such as histones H3 and H4, ribosomal proteins L15 and L12, chaperonin TCP1, 14-3-3 protein gamma, gibberellins A1, A3 and A8, strigolactones and abscisic acid, involved in epigenetic regulation, proteome remodeling, hormonal signaling and abiotic stress pathways showed their potential role during bud maturation. Candidates and pathways were validated employing interaction databases and targeted transcriptomics. These results increase our understanding of the molecular processes behind bud maturation, a key step towards improving timber production and natural pine forests management in a future scenario of climate change. However, further studies are necessary using different P. pinaster populations that show contrasting wood quality and stress tolerance in order to generalize the results.

GeLC-Orbitrap/MS and 2-DE-MALDI-TOF/TOF comparative proteomics analysis of seed cotyledons from the non-orthodox Quercus ilex tree species

Gel electrophoresis-based and shotgun approaches are the most employed proteomic platforms in plant biology research, with the latter replacing the former in the last years. We have compared 2-DE-MALDI-TOF/TOF and GeLC-Orbitrap/MS analyses using the same protein extracts from Quercus ilex cotyledons at different development stages. The results obtained (ProteomeXchange available data, PXD020603) showed that both platforms were complementary, showing common and specific proteins identified in each case, but leading to similar biological conclusions. Protein analysis identified 562 spots in gel-based (292 variables) and 2409 proteins in shotgun (560 variables), that were detected with both platforms and represent common key pathways related to maturation and germination. The main differences concern hormone metabolism, storage and late embryogenesis abundant proteins. Deeper proteome coverage was obtained with the shotgun approach, with a greater number of metabolic pathways represented, as gibberellin biosynthesis, not observed in the gel-based analysis. Nevertheless, several storage proteins, highly abundant in cotyledons and well represented in gel-based platform were not identified using the shotgun platform. These results support that when analyzing any plant biological process, the use of both platforms is complementary rather than redundant, that favors an in-depth proteomic analysis and a more confident biological interpretation of the data obtained.

In-depth analysis of the Quercus suber metabolome under drought stress and recovery reveals potential key metabolic players

Cork oak (Quercus suber L.) is a species of ecological, social and economic importance in the Mediterranean region. Given its xerophytic adaptability, the study of cork oak's response to drought stress conditions may provide important data in the global scenario of climate change. The mechanisms behind cork oak's adaptation to drought conditions can inform the design and development of tools to better manage this species under the changing climate patterns. Metabolomics is one of the most promising omics layers to capture a snapshot of a particular physiological state and to identify putative biomarkers of stress tolerance. Drastic changes were observed in the leaf metabolome of Q. suber between the different experimental conditions, namely at the beginning of the drought stress treatment, after one month under drought and post rehydration. All experimental treatments were analyzed through sPLS to inspect for global changes and stress and rehydration responses were analyzed independently for specific alterations. This allowed a more in-depth study and a search for biomarkers specific to a given hydric treatment. The metabolome analyses showed changes in both primary and secondary metabolism, but highlighted the role of secondary metabolism. In addition, a compound-specific response was observed in stress and rehydration. Key compounds such as L-phenylalanine and epigallocatechin 3-gallate were identified in relation to early drought response, terpenoid leonuridine and the flavonoid glycoside (-)-epicatechin-3'-O-glucuronide in long-term drought response, and flavone isoscoparine was identified in relation to the recovery process. The results here obtained provide novel insights into the biology of cork oak, highlighting pathways and metabolites potentially involved in the response of this species during drought and recovery that may be essential for its adaptation to long periods of drought. It is expected that this knowledge can encourage further functional studies in order to validate potential biomarkers of drought and recovery that maybe used to support decision-making in cork oak breeding programs.

Dissecting the Seed Maturation and Germination Processes in the Non-Orthodox Quercus ilex Species Based on Protein Signatures as Revealed by 2-DE Coupled to MALDI-TOF/TOF Proteomics Strategy

Unlike orthodox species, seed recalcitrance is poorly understood, especially at the molecular level. In this regard, seed maturation and germination were studied in the non-orthodox Quercus ilex by using a proteomics strategy based on two-dimensional gel electrophoresis coupled to matrix-assisted laser desorption ionization/time of flight (2-DE-MALDI-TOF).Cotyledons and embryo/radicle were sampled at different developmental stages, including early (M1–M3), middle (M4–M7), and late (M8–M9) seed maturation, and early (G1–G3) and late (G4–G5) germination. Samples corresponding to non-germinating, inviable, seeds were also included. Protein extracts were subjected to 2-dimensional gel electrophoresis (2-DE) and changes in the protein profiles were analyzed. Identified variable proteins were grouped according to their function, being the energy, carbohydrate, lipid, and amino acid metabolisms, together with protein fate, redox homeostasis, and response to stress are the most represented groups. Beyond the visual aspect, morphometry, weight, and water content, each stage had a specific protein signature. Clear tendencies for the different protein groups throughout the maturation and germination stages were observed for, respectively, cotyledon and the embryo axis. Proteins related to metabolism, translation, legumins, proteases, proteasome, and those stress related were less abundant in non-germinating seeds, it related to the loss of viability. Cotyledons were enriched with reserve proteins and protein-degrading enzymes, while the embryo axis was enriched with proteins of cell defense and rescue, including heat-shock proteins (HSPs) and antioxidants. The peaks of enzyme proteins occurred at the middle stages (M6–M7) in cotyledons and at late ones (M8–M9) in the embryo axis. Unlike orthodox seeds, proteins associated with glycolysis, tricarboxylic acid cycle, carbohydrate, amino acid and lipid metabolism are present at high levels in the mature seed and were maintained throughout the germination stages. The lack of desiccation tolerance in Q. ilex seeds may be associated with the repression of some genes, late embryogenesis abundant proteins being one of the candidates.

Proteomics, Holm Oak (Quercus ilex L.) and Other Recalcitrant and Orphan Forest Tree Species: How do They See Each Other?

Proteomics has had a big impact on plant biology, considered as a valuable tool for several forest species, such as Quercus, Pines, Poplars, and Eucalyptus. This review assesses the potential and limitations of the proteomics approaches and is focused on Quercus ilex as a model species and other forest tree species. Proteomics has been used with Q. ilex since 2003 with the main aim of examining natural variability, developmental processes, and responses to biotic and abiotic stresses as in other species of the genus Quercus or Pinus. As with the progress in techniques in proteomics in other plant species, the research in Q. ilex moved from 2-DE based strategy to the latest gel-free shotgun workflows. Experimental design, protein extraction, mass spectrometric analysis, confidence levels of qualitative and quantitative proteomics data, and their interpretation are a true challenge with relation to forest tree species due to their extreme orphan and recalcitrant (non-orthodox) nature. Implementing a systems biology approach, it is time to validate proteomics data using complementary techniques and integrate it with the -omics and classical approaches. The full potential of the protein field in plant research is quite far from being entirely exploited. However, despite the methodological limitations present in proteomics, there is no doubt that this discipline has contributed to deeper knowledge of plant biology and, currently, is increasingly employed for translational purposes.

Salt-stress induced proteomic changes of two contrasting alfalfa cultivars during germination stage

BACKGROUND

Alfalfa (Medicago sativa L.), the primary forage crop throughout the world, is sensitive to salt stress during the germination stage. To investigate the response of alfalfa to salt stress, a comprehensive proteomic analysis was performed comparing alfalfa cultivars that differ in salinity tolerance in the early seedling.

RESULTS

Five cultivars were examined for salt tolerance, and the most salt-tolerant cultivar, ZhongmuNo.3, and the most salt-sensitive cultivar, Daxiyang, were compared in terms of their physiological and proteomic responses. The two alfalfa cultivars seeds were exposed to 0 mmolL^-1 or 200 mmolL^-1 NaCl salt treatment for 10 days. Salt stress significantly reduced young seedling growth and the cotyledons' chlorophyll content; meanwhile, it increased the cotyledons' H₂ O₂ and malondialdehyde (MDA) levels, all of which were less adversely affected in ZhongmuNo.3 than in Daxiyang. A total of 51 spots (24 and 27 protein spots in the salt-sensitive and salt-tolerant cultivars, respectively) were identified as significantly differentially expressed using two-dimensional electrophoresis analysis. The proteins that were associated with salt tolerance included antioxidants/detoxifying enzymes, molecular chaperones, energy metabolic enzymes, a secondary metabolic enzyme, and pathogenesis-related proteins.

CONCLUSIONS

Under salt stress, ZhongmuNo.3 possessed a higher capacity for reactive oxygen species (ROS) scavenging, a more abundant energy supply, and stronger photosynthesis than the salt-sensitive cultivar Daxiyang, and these physiological processes may be the primary contributors to salt tolerance in ZhongmuNo.3. These advanced proteome data expand our knowledge of the physiology of the response of alfalfa to salt stress, providing a potentially valuable foundation for molecular breeding to enhance salt tolerance. © 2018 Society of Chemical Industry.

Metabolome Integrated Analysis of High-Temperature Response in Pinus radiata

The integrative omics approach is crucial to identify the molecular mechanisms underlying high-temperature response in non-model species. Based on future scenarios of heat increase, Pinus radiata plants were exposed to a temperature of 40°C for a period of 5 days, including recovered plants (30 days after last exposure to 40°C) in the analysis. The analysis of the metabolome using complementary mass spectrometry techniques (GC-MS and LC-Orbitrap-MS) allowed the reliable quantification of 2,287 metabolites. The analysis of identified metabolites and highlighter metabolic pathways across heat time exposure reveal the dynamism of the metabolome in relation to high-temperature response in P. radiata, identifying the existence of a turning point (on day 3) at which P. radiata plants changed from an initial stress response program (shorter-term response) to an acclimation one (longer-term response). Furthermore, the integration of metabolome and physiological measurements, which cover from the photosynthetic state to hormonal profile, suggests a complex metabolic pathway interaction network related to heat-stress response. Cytokinins (CKs), fatty acid metabolism and flavonoid and terpenoid biosynthesis were revealed as the most important pathways involved in heat-stress response in P. radiata, with zeatin riboside (ZR) and isopentenyl adenosine (iPA) as the key hormones coordinating these multiple and complex interactions. On the other hand, the integrative approach allowed elucidation of crucial metabolic mechanisms involved in heat response in P. radiata, as well as the identification of thermotolerance metabolic biomarkers (L-phenylalanine, hexadecanoic acid, and dihydromyricetin), crucial metabolites which can reschedule the metabolic strategy to adapt to high temperature.

System-wide analysis of short-term response to high temperature in Pinus radiata

Pinus radiata seedlings, the most widely planted pine species in the world, were exposed to temperatures within a range mimicking future scenarios based on current models of heat increase. The short-term heat response in P. radiata was studied in detail by exploring the metabolome, proteome and targeted transcriptome. The use of complementary mass spectrometry techniques, GC-MS and LC-Orbitrap-MS, together with novel bioinformatics tools allowed the reliable quantification of 2,075 metabolites and 901 protein groups. Integrative analyses of different functional levels and plant physiological status revealed a complex molecular interaction network of positive and negative correlations between proteins and metabolites involved in short-term heat response, including three main physiological functions as: 1) A hormone subnetwork, where fatty acids, flavonoids and hormones presented a key role; 2) An oxidoreductase subnetwork, including several dehydrogenase and peroxidase proteins; and 3) A heat shock protein subnetwork, with numerous proteins that contain a HSP20 domain, all of which were overexpressed at the transcriptional level. Integrated analysis pinpointed the basic mechanisms underlying the short-term physiological reaction of P. radiata during heat response. This approach was feasible in forest species and unmasked two novel candidate biomarkers of heat resistance, PHO1 and TRANSCRIPTION FACTOR APFI, and a MITOCHONDRIAL SMALL HEAT SHOCK PROTEIN, for use in future breeding programs.

System-level network analysis of nitrogen starvation and recovery in Chlamydomonas reinhardtii reveals potential new targets for increased lipid accumulation

BACKGROUND

Nitrogen starvation is known to cause drastic alterations in physiology and metabolism leading to the accumulation of lipid bodies in many microalgae, and it thus presents an important alternative for biofuel production. However, despite the importance of this process, the molecular mechanisms that mediate the metabolic remodeling induced by N starvation and especially by stress recovery are still poorly understood, and new candidates for bioengineering are needed to make this process useful for biofuel production.

RESULTS

We have studied the molecular changes involved in the adaptive mechanisms to N starvation and full recovery of the vegetative cells in the microalga Chlamydomonas reinhardtii during a four-day time course. High throughput mass spectrometry was employed to integrate the proteome and the metabolome with physiological changes. N starvation led to an accumulation of oil bodies and reduced Fv/Fm.. Distinct enzymes potentially participating in the carbon-concentrating mechanism (CAH7, CAH8, PEPC1) are strongly accumulated. The membrane composition is changed, as indicated by quantitative lipid profiles. A reprogramming of protein biosynthesis was observed by increased levels of cytosolic ribosomes, while chloroplastidic were dramatically reduced. Readdition of N led to, the identification of early responsive proteins mediating stress recovery, indicating their key role in regaining and sustaining normal vegetative growth. Analysis of the data with multivariate correlation analysis, Granger causality, and sparse partial least square (sPLS) provided a functional network perspective of the molecular processes. Cell growth and N metabolism were clearly linked by the branched chain amino acids, suggesting an important role in this stress. Lipid accumulation was also tightly correlated to the COP II protein, involved in vesicle and lysosome coating, and a major lipid droplet protein. This protein, together with other key proteins mediating signal transduction and adaption (BRI1, snRKs), constitute a series of new metabolic and regulatory targets.

CONCLUSIONS

This work not only provides new insights and corrects previous models by analyzing a complex dataset, but also increases our biochemical understanding of the adaptive mechanisms to N starvation in Chlamydomonas, pointing to new bioengineering targets for increased lipid accumulation, a key step for a sustainable and profitable microalgae-based biofuel production.

Systemic cold stress adaptation of Chlamydomonas reinhardtii

Chlamydomonas reinhardtii is one of the most important model organisms nowadays phylogenetically situated between higher plants and animals (Merchant et al. 2007). Stress adaptation of this unicellular model algae is in the focus because of its relevance to biomass and biofuel production. Here, we have studied cold stress adaptation of C. reinhardtii hitherto not described for this algae whereas intensively studied in higher plants. Toward this goal, high throughput mass spectrometry was employed to integrate proteome, metabolome, physiological and cell-morphological changes during a time-course from 0 to 120 h. These data were complemented with RT-qPCR for target genes involved in central metabolism, signaling, and lipid biosynthesis. Using this approach dynamics in central metabolism were linked to cold-stress dependent sugar and autophagy pathways as well as novel genes in C. reinhardtii such as CKIN1, CKIN2 and a hitherto functionally not annotated protein named CKIN3. Cold stress affected extensively the physiology and the organization of the cell. Gluconeogenesis and starch biosynthesis pathways are activated leading to a pronounced starch and sugar accumulation. Quantitative lipid profiles indicate a sharp decrease in the lipophilic fraction and an increase in polyunsaturated fatty acids suggesting this as a mechanism of maintaining membrane fluidity. The proteome is completely remodeled during cold stress: specific candidates of the ribosome and the spliceosome indicate altered biosynthesis and degradation of proteins important for adaptation to low temperatures. Specific proteasome degradation may be mediated by the observed cold-specific changes in the ubiquitinylation system. Sparse partial least squares regression analysis was applied for protein correlation network analysis using proteins as predictors and Fv/Fm, FW, total lipids, and starch as responses. We applied also Granger causality analysis and revealed correlations between proteins and metabolites otherwise not detectable. Twenty percent of the proteins responsive to cold are uncharacterized proteins. This presents a considerable resource for new discoveries in cold stress biology in alga and plants.