Correlation of mRNA and protein abundance in the developing maize leaf

Multi-omics analysis of aggregative multicellularity

All organisms have to carefully regulate their gene expression, not least during development. mRNA levels are often used as proxy for protein output; however, this approach ignores post-transcriptional effects. In particular, mRNA-protein correlation remains elusive for organisms that exhibit aggregative rather than clonal multicellularity. We addressed this issue by generating a paired transcriptomics and proteomics time series during the transition from uni-to multicellular stage in the social ameba Dictyostelium discoideum. Our data reveals that mRNA and protein levels correlate highly during unicellular growth, but decrease when multicellular development is initiated. This accentuates that transcripts alone cannot accurately predict protein levels. The dataset provides a useful resource to study gene expression during aggregative multicellular development. Additionally, our study provides an example of how to analyze and visualize mRNA and protein levels, which should be broadly applicable to other organisms and conditions.

Gene expression response under thermal stress in two Hawaiian corals is dominated by ploidy and genotype

Transcriptome data are frequently used to investigate coral bleaching; however, the factors controlling gene expression in natural populations of these species are poorly understood. We studied two corals, Montipora capitata and Pocillopora acuta, that inhabit the sheltered Kāne'ohe Bay, Hawai'i. M. capitata colonies in the bay are outbreeding diploids, whereas P. acuta is a mixture of clonal diploids and triploids. Populations were sampled from six reefs and subjected to either control (no stress), thermal stress, pH stress, or combined pH and thermal stress treatments. RNA-seq data were generated to test two competing hypotheses: (1) gene expression is largely independent of genotype, reflecting a shared treatment-driven response (TDE) or, (2) genotype dominates gene expression, regardless of treatment (GDE). Our results strongly support the GDE model, even under severe stress. We suggest that post-transcriptional processes (e.g., control of translation, protein turnover) modify the signal from the transcriptome, and may underlie the observed differences in coral bleaching sensitivity via the downstream proteome and metabolome.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein coexpression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Maize multi-omics reveal leaf water status controlling of differential transcriptomes, proteomes and hormones as mechanisms of age-dependent osmotic stress response in leaves

Drought-induced osmotic stress severely affects the growth and yield of maize. However, the mechanisms underlying the different responses of young and old maize leaves to osmotic stress remain unclear. To gain a systematic understanding of age-related stress responses, we compared osmotic-stress-induced changes in maize leaves of different ages using multi-omics approaches. After short-term osmotic stress, old leaves suffered more severe water deficits than young leaves. The adjustments of transcriptomes, proteomes, and hormones in response to osmotic stress were more dynamic in old leaves. Metabolic activities, stress signaling pathways, and hormones (especially abscisic acid) responded to osmotic stress in an age-dependent manner. We identified multiple functional clusters of genes and proteins with potential roles in stress adaptation. Old leaves significantly accumulated stress proteins such as dehydrin, aquaporin, and chaperones to cope with osmotic stress, accompanied by senescence-like cellular events, whereas young leaves exhibited an effective water conservation strategy mainly by hydrolyzing transitory starch and increasing proline production. The stress responses of individual leaves are primarily determined by their intracellular water status, resulting in differential transcriptomes, proteomes, and hormones. This study extends our understanding of the mechanisms underlying plant responses to osmotic stress.

Detection of the Arabidopsis Proteome and Its Post-translational Modifications and the Nature of the Unobserved (Dark) Proteome in PeptideAtlas

This study describes a new release of the Arabidopsis thaliana PeptideAtlas proteomics resource (build 2023-10) providing protein sequence coverage, matched mass spectrometry (MS) spectra, selected post-translational modifications (PTMs), and metadata. 70 million MS/MS spectra were matched to the Araport11 annotation, identifying ∼0.6 million unique peptides and 18,267 proteins at the highest confidence level and 3396 lower confidence proteins, together representing 78.6% of the predicted proteome. Additional identified proteins not predicted in Araport11 should be considered for the next Arabidopsis genome annotation. This release identified 5198 phosphorylated proteins, 668 ubiquitinated proteins, 3050 N-terminally acetylated proteins, and 864 lysine-acetylated proteins and mapped their PTM sites. MS support was lacking for 21.4% (5896 proteins) of the predicted Araport11 proteome: the "dark" proteome. This dark proteome is highly enriched for E3 ligases, transcription factors, and for certain (e.g., CLE, IDA, PSY) but not other (e.g., THIONIN, CAP) signaling peptides families. A machine learning model trained on RNA expression data and protein properties predicts the probability that proteins will be detected. The model aids in discovery of proteins with short half-life (e.g., SIG1,3 and ERF-VII TFs) and for developing strategies to identify the missing proteins. PeptideAtlas is linked to TAIR, tracks in JBrowse, and several other community proteomics resources.

Implications of tRNA abundance on translation elongation across bovine tissues

Introduction: Translation is a crucial stage of gene expression. It may also act as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are believed to be codon-biased to support increased protein production, in which quickly translated codons correspond to highly abundant tRNAs. Synonymous SNPs, considered to be silent due to the degeneracy of the genetic code, may shift protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Methods: Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNAseq) and ribosome profiling across bovine tissues in order to investigate the relationship between tRNA expression and slowed translation. Results: Moreover, we have identified genes modulated at transcriptional and/or translational levels underlying tissue-specific biological processes. We have also successfully defined pausing sites that depict the regulatory information encoded within the open reading frame of transcripts, which could be related to translation rate and facilitate proper protein folding. This work offers an atlas of distinctive pausing sites across three bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

Compartmentation of photosynthesis gene expression in C4 maize depends on time of day

Compared with the ancestral C3 state, C4 photosynthesis occurs at higher rates with improved water and nitrogen use efficiencies. In both C3 and C4 plants, rates of photosynthesis increase with light intensity and are maximal around midday. We determined that in the absence of light or temperature fluctuations, photosynthesis in maize (Zea mays) peaks in the middle of the subjective photoperiod. To investigate the molecular processes associated with these temporal changes, we performed RNA sequencing of maize mesophyll and bundle sheath strands over a 24-h time course. Preferential expression of C4 cycle genes in these cell types was strongest between 6 and 10 h after dawn when rates of photosynthesis were highest. For the bundle sheath, DNA motif enrichment and gene coexpression analyses suggested members of the DNA binding with one finger (DOF) and MADS (MINICHROMOSOME MAINTENANCE FACTOR 1/AGAMOUS/DEFICIENS/Serum Response Factor)-domain transcription factor families mediate diurnal fluctuations in C4 gene expression, while trans-activation assays in planta confirmed their ability to activate promoter fragments from bundle sheath expressed genes. The work thus identifies transcriptional regulators and peaks in cell-specific C4 gene expression coincident with maximum rates of photosynthesis in the maize leaf at midday.

Changing turn-over rates regulate abundance of tryptophan, GS biosynthesis, IAA transport and photosynthesis proteins in Arabidopsis growth defense transitions

BACKGROUND

Shifts in dynamic equilibria of the abundance of cellular molecules in plant-pathogen interactions need further exploration. We induced PTI in optimally growing Arabidopsis thaliana seedlings for 16 h, returning them to growth conditions for another 16 h.

METHODS

Turn-over and abundance of 99 flg22 responding proteins were measured chronologically using a stable heavy nitrogen isotope partial labeling strategy and targeted liquid chromatography coupled to mass spectrometry (PRM LC-MS). These experiments were complemented by measurements of mRNA and phytohormone levels.

RESULTS

Changes in synthesis and degradation rate constants (Ks and Kd) regulated tryptophane and glucosinolate, IAA transport, and photosynthesis-associated protein (PAP) homeostasis in growth/PTI transitions independently of mRNA levels. Ks values increased after elicitation while protein and mRNA levels became uncorrelated. mRNA returned to pre-elicitation levels, yet protein abundance remained at PTI levels even 16 h after media exchange, indicating protein levels were robust and unresponsive to transition back to growth. The abundance of 23 PAPs including FERREDOXIN-NADP( +)-OXIDOREDUCTASE (FNR1) decreased 16 h after PAMP exposure, their depletion was nearly abolished in the myc234 mutant. FNR1 Kd increased as mRNA levels decreased early in PTI, its Ks decreased in prolonged PTI. FNR1 Kd was lower in myc234, mRNA levels decreased as in wild type.

CONCLUSIONS

Protein Kd and Ks values change in response to flg22 exposure and constitute an additional layer of protein abundance regulation in growth defense transitions next to changes in mRNA levels. Our results suggest photosystem remodeling in PTI to direct electron flow away from the photosynthetic carbon reaction towards ROS production as an active defense mechanism controlled post-transcriptionally and by MYC2 and homologs. Target proteins accumulated later and PAP and auxin/IAA depletion was repressed in myc234 indicating a positive effect of the transcription factors in the establishment of PTI.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein co-expression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3.6%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Highlights

At the level of individual genes, the abundance of transcripts and proteins is weakly correlated within a species ( ρ = 0.165). While the proteome is not imprinted by population structure, co-expression patterns recapitulate the cellular functional landscapeWild populations exhibit a higher abundance of respiration-related proteins compared to domesticated populationsLoci that influence protein abundance differ from those that impact transcript levels, likely because of protein turnover.

Targeted mutagenesis of the vacuolar H+ translocating pyrophosphatase gene reduces grain chalkiness in rice

Grain chalkiness is a major concern in rice production because it impacts milling yield and cooking quality, eventually reducing market value of the rice. A gene encoding vacuolar H+ translocating pyrophosphatase (V-PPase) is a major quantitative trait locus in indica rice, controlling grain chalkiness. Higher transcriptional activity of this gene is associated with increased chalk content. However, whether the suppression of V-PPase could reduce chalkiness is not clear. Furthermore, natural variation in the chalkiness of japonica rice has not been linked with V-PPase. Here, we describe promoter targeting of the japonica V-PPase allele that led to reduced grain chalkiness and the development of more translucent grains. Disruption of a putative GATA element by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 suppressed V-PPase activity, reduced grain chalkiness and impacted post-germination growth that could be rescued by the exogenous supply of sucrose. The mature grains of the targeted lines showed a much lower percentage of large or medium chalk. Interestingly, the targeted lines developed a significantly lower chalk under heat stress, a major inducer of grain chalk. Metabolomic analysis showed that pathways related to starch and sugar metabolism were affected in the developing grains of the targeted lines that correlated with higher inorganic pyrophosphate and starch contents and upregulation of starch biosynthesis genes. In summary, we show a biotechnology approach of reducing grain chalkiness in rice by downregulating the transcriptional activity of V-PPase that presumably leads to altered metabolic rates, including starch biosynthesis, resulting in more compact packing of starch granules and formation of translucent rice grains.

Seven plant capacities to adapt to abiotic stress

Abiotic stresses such as drought and heat continue to impact crop production in a warming world. This review distinguishes seven inherent capacities that enable plants to respond to abiotic stresses and continue growing, although at a reduced rate, to achieve a productive yield. These are the capacities to selectively take up essential resources, store them and supply them to different plant parts, generate the energy required for cellular functions, conduct repairs to maintain plant tissues, communicate between plant parts, manage existing structural assets in the face of changed circumstances, and shape-shift through development to be efficient in different environments. By illustration, we show how all seven plant capacities are important for reproductive success of major crop species during drought, salinity, temperature extremes, flooding, and nutrient stress. Confusion about the term 'oxidative stress' is explained. This allows us to focus on the strategies that enhance plant adaptation by identifying key responses that can be targets for plant breeding.

Effects of Turmeric Extract on Age-Related Skeletal Muscle Atrophy in Senescence-Accelerated Mice

Muscle atrophy is one of the main causes of sarcopenia—the age-related loss of skeletal muscle. In this study, we investigated the effect of turmeric (Curcuma longa) extract (TE) supplementation on age-related muscle atrophy in a senescence-accelerated mouse model and explored the underlying mechanisms. Twenty-six-week-old male, senescence-accelerated mouse resistant (SAMR) mice received the AIN-93G basal diet, while twenty-six-week-old male, senescence-accelerated mouse prone 8 (SAMP8) mice received the AIN-93G basal diet or a 2% TE powder-supplemented diet for ten weeks. Our findings revealed that TE supplementation showed certain effects on ameliorating the decrease in body weight, tibialis anterior weight, and mesenteric fat tissue weight in SAMP8 mice. TE improved gene expression in the glucocorticoid receptor-FoxO signaling pathway in skeletal muscle, including redd1, klf15, foxo1, murf1, and mafbx. Furthermore, TE might have the certain potential on improving the dynamic balance between anabolic and catabolic processes by inhibiting the binding of glucocorticoid receptor or FoxO1 to the glucocorticoid response element or FoxO-binding element in the MuRF1 promoter in skeletal muscle, thereby promoting muscle mass and strength, and preventing muscle atrophy and sarcopenia prevention. Moreover, TE may have reduced mitochondrial damage and maintained cell growth and division by downregulating the mRNA expression of the genes mfn2 and tsc2. Thus, the results indicated TE’s potential for preventing age-related muscle atrophy and sarcopenia.

Impact of Growth Rate on the Protein-mRNA Ratio in Pseudomonas aeruginosa

Our understanding of how bacterial pathogens colonize and persist during human infection has been hampered by the limited characterization of bacterial physiology during infection and a research bias toward in vitro, fast-growing bacteria. Recent research has begun to address these gaps in knowledge by directly quantifying bacterial mRNA levels during human infection, with the goal of assessing microbial community function at the infection site. However, mRNA levels are not always predictive of protein levels, which are the primary functional units of a cell. Here, we used carefully controlled chemostat experiments to examine the relationship between mRNA and protein levels across four growth rates in the bacterial pathogen Pseudomonas aeruginosa. We found a genome-wide positive correlation between mRNA and protein abundances across all growth rates, with genes required for P. aeruginosa viability having stronger correlations than nonessential genes. We developed a statistical method to identify genes whose mRNA abundances poorly predict protein abundances and calculated an RNA-to-protein (RTP) conversion factor to improve mRNA predictions of protein levels. The application of the RTP conversion factor to publicly available transcriptome data sets was highly robust, enabling the more accurate prediction of P. aeruginosa protein levels across strains and growth conditions. Finally, the RTP conversion factor was applied to P. aeruginosa human cystic fibrosis (CF) infection transcriptomes to provide greater insights into the functionality of this bacterium in the CF lung. This study addresses a critical problem in infection microbiology by providing a framework for enhancing the functional interpretation of bacterial human infection transcriptome data. IMPORTANCE Our understanding of bacterial physiology during human infection is limited by the difficulty in assessing bacterial function at the infection site. Recent studies have begun to address this question by quantifying bacterial mRNA levels in human-derived samples using transcriptomics. One challenge for these studies is the poor predictivity of mRNA for protein levels for some genes. Here, we addressed this challenge by measuring the transcriptomes and proteomes of P. aeruginosa grown at four growth rates. Our results revealed that the growth rate does not impact the genome-wide correlation of mRNA and protein levels. We used statistical methods to identify the genes for which mRNA and protein were poorly correlated and developed an RNA-to-protein (RTP) conversion factor that improved the predictivity of protein levels across strains and growth conditions. Our results provide new insights into mRNA-protein correlations and tools to enhance our understanding of bacterial physiology from transcriptome data.

The roles of plant proteases and protease inhibitors in drought response: a review

Upon exposure to drought, plants undergo complex signal transduction events with concomitant changes in the expression of genes, proteins and metabolites. For example, proteomics studies continue to identify multitudes of drought-responsive proteins with diverse roles in drought adaptation. Among these are protein degradation processes that activate enzymes and signalling peptides, recycle nitrogen sources, and maintain protein turnover and homeostasis under stressful environments. Here, we review the differential expression and functional activities of plant protease and protease inhibitor proteins under drought stress, mainly focusing on comparative studies involving genotypes of contrasting drought phenotypes. We further explore studies of transgenic plants either overexpressing or repressing proteases or their inhibitors under drought conditions and discuss the potential roles of these transgenes in drought response. Overall, the review highlights the integral role of protein degradation during plant survival under water deficits, irrespective of the genotypes' level of drought resilience. However, drought-sensitive genotypes exhibit higher proteolytic activities, while drought-tolerant genotypes tend to protect proteins from degradation by expressing more protease inhibitors. In addition, transgenic plant biology studies implicate proteases and protease inhibitors in various other physiological functions under drought stress. These include the regulation of stomatal closure, maintenance of relative water content, phytohormonal signalling systems including abscisic acid (ABA) signalling, and the induction of ABA-related stress genes, all of which are essential for maintaining cellular homeostasis under water deficits. Therefore, more validation studies are required to explore the various functions of proteases and their inhibitors under water limitation and their contributions towards drought adaptation.

Cell-specific pathways recruited for symbiotic nodulation in the Medicago truncatula legume

Medicago truncatula is a model legume species that has been studied for decades to understand the symbiotic relationship between legumes and soil bacteria collectively named rhizobia. This symbiosis called nodulation is initiated in roots with the infection of root hair cells by the bacteria, as well as the initiation of nodule primordia from root cortical, endodermal, and pericycle cells, leading to the development of a new root organ, the nodule, where bacteria fix and assimilate the atmospheric dinitrogen for the benefit of the plant. Here, we report the isolation and use of the nuclei from mock and rhizobia-inoculated roots for the single nuclei RNA-seq (sNucRNA-seq) profiling to gain a deeper understanding of early responses to rhizobial infection in Medicago roots. A gene expression map of the Medicago root was generated, comprising 25 clusters, which were annotated as specific cell types using 119 Medicago marker genes and orthologs to Arabidopsis cell-type marker genes. A focus on root hair, cortex, endodermis, and pericycle cell types, showing the strongest differential regulation in response to a short-term (48 h) rhizobium inoculation, revealed not only known genes and functional pathways, validating the sNucRNA-seq approach, but also numerous novel genes and pathways, allowing a comprehensive analysis of early root symbiotic responses at a cell type-specific level.

Multi-Omics Approaches and Resources for Systems-Level Gene Function Prediction in the Plant Kingdom

In higher plants, the complexity of a system and the components within and among species are rapidly dissected by omics technologies. Multi-omics datasets are integrated to infer and enable a comprehensive understanding of the life processes of organisms of interest. Further, growing open-source datasets coupled with the emergence of high-performance computing and development of computational tools for biological sciences have assisted in silico functional prediction of unknown genes, proteins and metabolites, otherwise known as uncharacterized. The systems biology approach includes data collection and filtration, system modelling, experimentation and the establishment of new hypotheses for experimental validation. Informatics technologies add meaningful sense to the output generated by complex bioinformatics algorithms, which are now freely available in a user-friendly graphical user interface. These resources accentuate gene function prediction at a relatively minimal cost and effort. Herein, we present a comprehensive view of relevant approaches available for system-level gene function prediction in the plant kingdom. Together, the most recent applications and sought-after principles for gene mining are discussed to benefit the plant research community. A realistic tabulation of plant genomic resources is included for a less laborious and accurate candidate gene discovery in basic plant research and improvement strategies.

Identification of stably expressed reference genes for expression studies in Arabidopsis thaliana using mass spectrometry-based label-free quantification

Arabidopsis thaliana has been used regularly as a model plant in gene expression studies on transcriptional reprogramming upon pathogen infection, such as that by Pseudomonas syringae pv. tomato DC3000 (Pst DC3000), or when subjected to stress hormone treatments including jasmonic acid (JA), salicylic acid (SA), and abscisic acid (ABA). Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) has been extensively employed to quantitate these gene expression changes. However, the accuracy of the quantitation is largely dependent on the stability of the expressions of reference genes used for normalization. Recently, RNA sequencing (RNA-seq) has been widely used to mine stably expressed genes for use as references in RT-qPCR. However, the amplification step in RNA-seq creates an intrinsic bias against those genes with relatively low expression levels, and therefore does not provide an accurate quantification of all expressed genes. In this study, we employed mass spectrometry-based label-free quantification (LFQ) in proteomic analyses to identify those proteins with abundances unaffected by Pst DC3000 infection. We verified, using RT-qPCR, that the levels of their corresponding mRNAs were also unaffected by Pst DC3000 infection. Compared to commonly used reference genes for expression studies in A. thaliana upon Pst DC3000 infection, the candidate reference genes reported in this study generally have a higher expression stability. In addition, using RT-qPCR, we verified that the mRNAs of the candidate reference genes were stably expressed upon stress hormone treatments including JA, SA, and ABA. Results indicated that the candidate genes identified here had stable expressions upon these stresses and are suitable to be used as reference genes for RT-qPCR. Among the 18 candidate reference genes reported in this study, many of them had greater expression stability than the commonly used reference genes, such as ACT7, in previous studies. Here, besides proposing more appropriate reference genes for Arabidopsis expression studies, we also demonstrated the capacity of mass spectrometry-based LFQ to quantify protein abundance and the possibility to extend protein expression studies to the transcript level.

Suppressive Effects of Turmeric Extract on Muscle Atrophy in Dexamethasone-Treated Mice and Myotubes

Sarcopenia is the decline in skeletal muscle mass, strength, and functions, which decreases the quality of life in elderly people. This study investigated the suppressive effect of turmeric (Curcuma longa) extract (TE) on muscle atrophy in dexamethasone (DEX)-treated mice and C2C12 myotubes. DEX treatment significantly decreased the muscle weight and significantly increased Fbxo32 and Murf1 expression in mice, and these changes were suppressed by the supplementation of an AIN-93 based diet with 2% TE. A similar pattern was observed in FBXO32 and MuRF1 protein expression. In C2C12 myotubes, DEX treatment significantly increased FBXO32 and MuRF1 gene and protein expression, and these increases were significantly suppressed by TE supplementation at a concentration of 200 µg/mL. Furthermore, one of the five TE fractions, which were separated by high-performance liquid chromatography had a similar effect with TE supplementation. The present study proposes the suppressive effect of turmeric on sarcopenia.

Non-target molecular network and putative genes of flavonoid biosynthesis in Erythrina velutina Willd., a Brazilian semiarid native woody plant

Erythrina velutina is a Brazilian native tree of the Caatinga (a unique semiarid biome). It is widely used in traditional medicine showing anti-inflammatory and central nervous system modulating activities. The species is a rich source of specialized metabolites, mostly alkaloids and flavonoids. To date, genomic information, biosynthesis, and regulation of flavonoids remain unknown in this woody plant. As part of a larger ongoing research goal to better understand specialized metabolism in plants inhabiting the harsh conditions of the Caatinga, the present study focused on this important class of bioactive phenolics. Leaves and seeds of plants growing in their natural habitat had their metabolic and proteomic profiles analyzed and integrated with transcriptome data. As a result, 96 metabolites (including 43 flavonoids) were annotated. Transcripts of the flavonoid pathway totaled 27, of which EvCHI, EvCHR, EvCHS, EvCYP75A and EvCYP75B1 were identified as putative main targets for modulating the accumulation of these metabolites. The highest correspondence of mRNA vs. protein was observed in the differentially expressed transcripts. In addition, 394 candidate transcripts encoding for transcription factors distributed among the bHLH, ERF, and MYB families were annotated. Based on interaction network analyses, several putative genes of the flavonoid pathway and transcription factors were related, particularly TFs of the MYB family. Expression patterns of transcripts involved in flavonoid biosynthesis and those involved in responses to biotic and abiotic stresses were discussed in detail. Overall, these findings provide a base for the understanding of molecular and metabolic responses in this medicinally important species. Moreover, the identification of key regulatory targets for future studies aiming at bioactive metabolite production will be facilitated.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

Transcriptome Analysis Reveals the Regulatory Networks of Cytokinin in Promoting Floral Feminization in Castanea henryi

Castanea henryi is a monoecious plant with a low female-to-male ratio, which limits its yield. The phytohormone cytokinin (CK) plays a crucial role in flower development, especially gynoecium development. Here, the feminizing effect of CK on the development of C. henryi was confirmed by the exogenous spraying of N-(2-chloro-4-pyridyl)-N'-phenylurea (CPPU). Spraying CPPU at 125 mg·L^-1 thrice changed the male catkin into a pure female catkin, whereas at 5 mg·L^-1 and 25 mg·L^-1, only a part of the male catkin was transformed into a female catkin. A comparative transcriptome analysis of male catkins subjected to CPPU was performed to study the mechanism of the role of CKs in sex differentiation. Using Pearson's correlation analysis between hormone content and hormone synthesis gene expression, four key genes, LOG1, LOG3, LOG7 and KO, were identified in the CK and GA synthesis pathways. Moreover, a hub gene in the crosstalk between JA and the other hormone signaling pathways, MYC2, was identified, and 15 flowering-related genes were significantly differentially expressed after CPPU treatment. These results suggest that CK interacts with other phytohormones to determine the sex of C. henryi, and CK may directly target floral organ recognition genes to control flower sex.

Advances and applications of single-cell omics technologies in plant research

Single-cell sequencing approaches reveal the intracellular dynamics of individual cells and answer biological questions with high-dimensional catalogs of millions of cells, including genomics, transcriptomics, chromatin accessibility, epigenomics, and proteomics data across species. These emerging yet thriving technologies have been fully embraced by the field of plant biology, with a constantly expanding portfolio of applications. Here, we introduce the current technical advances used for single-cell omics, especially single-cell genome and transcriptome sequencing. Firstly, we overview methods for protoplast and nucleus isolation and genome and transcriptome amplification. Subsequently, we use well-executed benchmarking studies to highlight advances made through the application of single-cell omics techniques. Looking forward, we offer a glimpse of additional hurdles and future opportunities that will introduce broad adoption of single-cell sequencing with revolutionary perspectives in plant biology.

To the proteome and beyond: advances in single-cell omics profiling for plant systems

Recent advances in single-cell proteomics for animal systems could be adapted for plants to increase our understanding of plant development, response to stimuli, and cell-to-cell signaling.

How Stress Affects Your Budget-Stress Impacts on Starch Metabolism

Starch is a polysaccharide that is stored to be used in different timescales. Transitory starch is used during nighttime when photosynthesis is unavailable. Long-term starch is stored to support vegetative or reproductive growth, reproduction, or stress responses. Starch is not just a reserve of energy for most plants but also has many other roles, such as promoting rapid stomatal opening, making osmoprotectants, cryoprotectants, scavengers of free radicals and signals, and reverting embolised vessels. Biotic and abiotic stress vary according to their nature, strength, duration, developmental stage of the plant, time of the day, and how gradually they develop. The impact of stress on starch metabolism depends on many factors: how the stress impacts the rate of photosynthesis, the affected organs, how the stress impacts carbon allocation, and the energy requirements involved in response to stress. Under abiotic stresses, starch degradation is usually activated, but starch accumulation may also be observed when growth is inhibited more than photosynthesis. Under biotic stresses, starch is usually accumulated, but the molecular mechanisms involved are largely unknown. In this mini-review, we explore what has been learned about starch metabolism and plant stress responses and discuss the current obstacles to fully understanding their interactions.

Ribosome profiling elucidates differential gene expression in bundle sheath and mesophyll cells in maize

The efficiencies offered by C4 photosynthesis have motivated efforts to understand its biochemical, genetic, and developmental basis. Reactions underlying C4 traits in most C4 plants are partitioned between two cell types, bundle sheath (BS), and mesophyll (M) cells. RNA-seq has been used to catalog differential gene expression in BS and M cells in maize (Zea mays) and several other C4 species. However, the contribution of translational control to maintaining the distinct proteomes of BS and M cells has not been addressed. In this study, we used ribosome profiling and RNA-seq to describe translatomes, translational efficiencies, and microRNA abundance in BS- and M-enriched fractions of maize seedling leaves. A conservative interpretation of our data revealed 182 genes exhibiting cell type-dependent differences in translational efficiency, 31 of which encode proteins with core roles in C4 photosynthesis. Our results suggest that non-AUG start codons are used preferentially in upstream open reading frames of BS cells, revealed mRNA sequence motifs that correlate with cell type-dependent translation, and identified potential translational regulators that are differentially expressed. In addition, our data expand the set of genes known to be differentially expressed in BS and M cells, including genes encoding transcription factors and microRNAs. These data add to the resources for understanding the evolutionary and developmental basis of C4 photosynthesis and for its engineering into C3 crops.

Genotype-dependent and heat-induced grain chalkiness in rice correlates with the expression patterns of starch biosynthesis genes

Starch biosynthesis is a complex process underlying grain chalkiness in rice in a genotype-dependent manner. Coordinated expression of starch biosynthesis genes is important for producing translucent rice grains, while disruption in this process leads to opaque or chalky grains. To better understand the dynamics of starch biosynthesis genes in grain chalkiness, six rice genotypes showing variable chalk levels were subjected to gene expression analysis during reproductive stages. In the chalky genotypes, peak expression of the large subunit genes of ADP-glucose pyrophosphorylase (AGPase), encoding the first key step in starch biosynthesis, occurred in the stages before grain filling commenced, creating a gap with the upregulation of starch synthase genes, granule bound starch synthase I (GBSSI) and starch synthase IIA (SSIIA). Whereas, in low-chalk genotypes, AGPase large subunit genes expressed at later stages, generally following the expression patterns of GBSSI and SSIIA. However, heat treatment altered the expression in a genotype-dependent manner that was accompanied by transformed grain morphology and increased chalkiness. The suppression of AGPase subunit genes during early grain filling stages was observed in the chalky genotypes or upon heat treatment, which could result in a limited pool of ADP-Glucose for synthesizing amylose and amylopectin, the major components of the starch. This suboptimal starch biosynthesis process could subsequently lead to inefficient grain filling and air pockets that contribute to chalkiness. In summary, this study suggests a mechanism of grain chalkiness based on the expression patterns of the starch biosynthesis genes in rice.

Correlation of Cry1Ac mRNA and protein abundance in transgenic Gossypium hirsutum plant

Transcription and translation in eukaryotes are distinct processes of the molecular cascade leading to protein production from genetic material. However, establishing correlation between mRNA expression and protein abundance, the end results of the two processes of central dogma, remains a challenge. For transgenic plants, such correlation between mRNA and protein expression serves as a guide to design the transgene, in particular the choices of promoter and codon usage to ensure stable expression of the target protein in relevant tissues under various stress conditions. To elucidate level of mRNA-protein correlation in a commercial transgenic cotton plant Gossypium hirsutum, Bollgard II^® (MON15985), we present the results of Cry1Ac protein expression correlating with corresponding mRNA levels. Protein was quantitated using a home-grown validated ELISA assay with a monoclonal-polyclonal antibody pair, whereas mRNA level was detected by a real-time quantitative PCR assay using standardized reference genes. Our results indicate that protein and mRNA levels are highly correlated in the leaves, but not in squares and stem. The correlations seem to be consistent between young and mature leaves and increase over time of harvesting of samples from months 1-3. These findings demonstrate that transcript level measurement could serve as a proxy to protein abundance for this commercially important cotton species, particularly for leaf tissues which are the most vulnerable organs to cotton bollworms and other pathogens.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02828-2.

Cellular and transcriptomic analyses reveal two-staged chloroplast biogenesis underpinning photosynthesis build-up in the wheat leaf

BACKGROUND

The developmental gradient in monocot leaves has been exploited to uncover leaf developmental gene expression programs and chloroplast biogenesis processes. However, the relationship between the two is barely understood, which limits the value of transcriptome data to understand the process of chloroplast development.

RESULTS

Taking advantage of the developmental gradient in the bread wheat leaf, we provide a simultaneous quantitative analysis for the development of mesophyll cells and of chloroplasts as a cellular compartment. This allows us to generate the first biologically-informed gene expression map of this leaf, with the entire developmental gradient from meristematic to fully differentiated cells captured. We show that the first phase of plastid development begins with organelle proliferation, which extends well beyond cell proliferation, and continues with the establishment and then the build-up of the plastid genetic machinery. The second phase is marked by the development of photosynthetic chloroplasts which occupy the available cellular space. Using a network reconstruction algorithm, we predict that known chloroplast gene expression regulators are differentially involved across those developmental stages.

CONCLUSIONS

Our analysis generates both the first wheat leaf transcriptional map and one of the most comprehensive descriptions to date of the developmental history of chloroplasts in higher plants. It reveals functionally distinct plastid and chloroplast development stages, identifies processes occurring in each of them, and highlights our very limited knowledge of the earliest drivers of plastid biogenesis, while providing a basis for their future identification.

Transcription initiation as a control point in plastid gene expression

The extensive processing and protein-assisted stabilization of transcripts have been taken as evidence for a viewpoint that the control of gene expression had shifted entirely in evolution from transcriptional in the bacterial endosymbiont to posttranscriptional in the plastid. This suggestion is however at odds with many observations on plastid gene transcription. Chloroplasts of flowering plants and mosses contain two or more RNA polymerases with distinct promoter preference and division of labor for the coordinated synthesis of plastid RNAs. Plant and algal plastids further possess multiple nonredundant sigma factors that function as transcription initiation factors. The controlled accumulation of plastid sigma factors and modification of their activity by sigma-binding proteins and phosphorylation constitute additional transcriptional regulatory strategies. Plant and algal plastids also contain dedicated one- or two-component transcriptional regulators. Transcription initiation thus continues to form a critical control point at which varied developmental and environmental signals intersect with plastid gene expression.

Comparison of protein and peptide fractionation approaches in protein identification and quantification from Saccharomyces cerevisiae

Shotgun proteomics is a high-throughput technology which has been developed with the aim of investigating the maximum number of proteins in cells in a given experiment. However, protein discovery and data generation vary in depth and coverage when different technical strategies are selected. In this study, three different sample preparation approaches, and peptide or protein fractionation methods, were applied to identify and quantify proteins from log-phase yeast lysate: sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), filter-aided sample preparation coupled with gas phase fractionation (FASP-GPF), and FASP - high pH reversed phase fractionation (HpH). Fractions were initially analyzed and compared using nanoflow liquid chromatography - tandem mass spectrometry (nanoLC-MS/MS) employing data dependent acquisition on a linear ion trap instrument. The number of fractions and analytical replicates was adjusted so that each experiment used a similar amount of mass spectrometric instrument time. A second set of experiments was performed, comparing FASP-GPF, SDS-PAGE and FASP-HpH using a Q Exactive Orbitrap mass spectrometer. Compared with results from the linear ion trap mass spectrometer, the use of a Q Exactive Orbitrap mass spectrometer enabled a substantial increase in protein identifications, and an even greater increase in peptide identifications. This shows that the main advantage of using the higher resolution instrument is in increased proteome coverage. A total of 1035, 1357 and 2134 proteins were separately identified by FASP-GPF, SDS-PAGE and FASP-HpH. Combining results from the Orbitrap experiments, there were a total of 2269 proteins found, with 94% of them identified using the FASP-HpH method. Therefore, the FASP-HpH method is the optimal choice among these approaches, when applied to this type of sample.

Reshaping of the Arabidopsis thaliana Proteome Landscape and Co-regulation of Proteins in Development and Immunity

Proteome remodeling is a fundamental adaptive response, and proteins in complexes and functionally related proteins are often co-expressed. Using a deep sampling strategy we define core proteomes of Arabidopsis thaliana tissues with around 10 000 proteins per tissue, and absolutely quantify (copy numbers per cell) nearly 16 000 proteins throughout the plant lifecycle. A proteome-wide survey of global post-translational modification revealed amino acid exchanges pointing to potential conservation of translational infidelity in eukaryotes. Correlation analysis of protein abundance uncovered potentially new tissue- and age-specific roles of entire signaling modules regulating transcription in photosynthesis, seed development, and senescence and abscission. Among others, the data suggest a potential function of RD26 and other NAC transcription factors in seed development related to desiccation tolerance as well as a possible function of cysteine-rich receptor-like kinases (CRKs) as ROS sensors in senescence. All of the components of ribosome biogenesis factor (RBF) complexes were found to be co-expressed in a tissue- and age-specific manner, indicating functional promiscuity in the assembly of these less-studied protein complexes in Arabidopsis.Furthermore, we characterized detailed proteome remodeling in basal immunity by treating Arabidopsis seeldings with flg22. Through simultaneously monitoring phytohormone and transcript changes upon flg22 treatment, we obtained strong evidence of suppression of jasmonate (JA) and JA-isoleucine (JA-Ile) levels by deconjugation and hydroxylation by IAA-ALA RESISTANT3 (IAR3) and JASMONATE-INDUCED OXYGENASE 2 (JOX2), respectively, under the control of JASMONATE INSENSITIVE 1 (MYC2), suggesting an unrecognized role of a new JA regulatory switch in pattern-triggered immunity. Taken together, the datasets generated in this study present extensive coverage of the Arabidopsis proteome in various biological scenarios, providing a rich resource available to the whole plant science community.

Multi-omics analysis of cellular pathways involved in different rapid growth stages of moso bamboo

Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau) is a rapidly growing grass of industrial and ecological importance. However, the molecular mechanisms of its remarkable growth are not well understood. In this study, we investigated the early-stage growth of moso bamboo shoots and defined three different growth stages based on histological and biochemical analyses, namely, starting of cell division (SD), rapid division (RD) and rapid elongation (RE). Further analyses on potentially relevant cellular pathways in these growth stages using multi-omics approaches such as transcriptomics and proteomics revealed the involvement of multiple cellular pathways, including DNA replication, repair and ribosome biogenesis. A total of 8045 differentially expressed genes (DEGs) and 1053 differentially expressed proteins (DEPs) were identified in our analyses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of detected DEGs identified several key biological pathways such as phytohormone metabolism, signal transduction, cell wall development and carbohydrate metabolism. The comparative analysis of proteins displayed that a total of 213 DEPs corresponded with DEGs and 3 significant expression profiles that could be promoting the fast growth of bamboo internodes. Moreover, protein-protein interaction network prediction analysis is suggestive of the involvement of five major proteins of signal transduction, DNA synthesis and RNA transcription, and may act as key elements responsible for the rapid shoot growth. Our work exploits multi-omics and bioinformatic approaches to unfurl the complexity of molecular networks involved in the rapid growth of moso bamboo and opens up questions related to the interactions between the functions played by individual molecular pathway.

Transgene behavior in Zea mays L. crosses across different genetic backgrounds: Segregation patterns, cry1Ab transgene expression, insecticidal protein concentration and bioactivity against insect pests

Brazil and South Africa, countries with economies in transition, are characterized by a dual agrarian structure co-occurring, sometimes, alongside in the same region. Large-scale commercial farming produces crops for export to global markets on the one hand, and small-scale farming, on the other hand, provides for subsistence and local markets. In both systems, maize (Zea mays) is a key crop for these two countries. For the commercial system, maize is a commodity crop while for the small-scale system it is the prime staple crop. In commercial systems, farmers predominantly grow genetically modified (GM) hybrid maize. In small-scale systems, farmers rely on open pollinated varieties (OPVs) and/or landraces and are largely dependent on seed saving systems. The aim of this study was to understand the relationship between transgene expression rates, the resulting concentrations of the transgene product (Bt protein) and its bioactivity in insect pests following transgene flow from GM hybrid maize into non-genetically modified, non-GM near-isogenic maize hybrid (ISO) and OPVs. We modeled segregation patterns and measured cry1Ab transgene expression (mRNA quantification), Cry1Ab protein concentration and pest survival. Two groups of F1, F2 crosses and backcrosses with GM, ISO and OPV maize varieties from Brazil and South Africa were used. Bioassays with the larvae of two lepidopteran maize pest species, Helicoverpa armigera and Spodoptera littoralis, were carried out. Overall, the cry1Ab transgene outcrossed effectively into the genetic backgrounds tested. The cry1Ab transgene was stably expressed in both ISO and OPV genetic backgrounds. Transgene introgression led to consistent, though highly variable, concentrations of Cry1Ab toxins that were similar to those observed in GM parental maize. Most crosses, but not all, suggested the expected Mendelian segregation pattern. Transgene expression rates were significantly higher than expected from purely Mendelian segregation in the South African crosses. In South African materials, ISO and OPV crosses produced significantly lower Cry1Ab concentrations compared to the GM parental maize. The Cry1Ab toxins from crosses were bioactive and induced mortality rates of ≥92.19% in H. armigera and ≥40.63% in S. littoralis after a period of only 4 days. However, no correlations were observed between the quantitation of mRNA for cry1Ab and the corresponding Cry1Ab protein concentrations, nor between the Cry1Ab concentrations and insect mortality rates across different genetic backgrounds. We therefore suggest that while transcription of the cry1Ab transgene reliably determines the presence of Cry1Ab protein, mRNA levels do not reflect, by themselves, the end Cry1Ab protein concentrations found in the plant. Because predictably high Cry1Ab concentrations are a key prerequisite for effective insect resistance management (IRM) programs, this observation raises questions about the effectiveness of such programs in scenarios with complex crop genetic backgrounds. On the other hand, confirmed bioactivity in all crosses should be expected to impact small farmer's selection behavior, unknowingly favoring the insecticidal trait. This may lead to a fixation of the trait in the wider population, and may influence the genetic diversity of varieties maintained by small-scale farmers.

Regulation of Carbon Metabolism by Environmental Conditions: A Perspective From Diatoms and Other Chromalveolates

Diatoms belong to a major, diverse and species-rich eukaryotic clade, the Heterokonta, within the polyphyletic chromalveolates. They evolved as a result of secondary endosymbiosis with one or more Plantae ancestors, but their precise evolutionary history is enigmatic. Nevertheless, this has conferred them with unique structural and biochemical properties that have allowed them to flourish in a wide range of different environments and cope with highly variable conditions. We review the effect of pH, light and dark, and CO2 concentration on the regulation of carbon uptake and assimilation. We discuss the regulation of the Calvin-Benson-Bassham cycle, glycolysis, lipid synthesis, and carbohydrate synthesis at the level of gene transcripts (transcriptomics), proteins (proteomics) and enzyme activity. In contrast to Viridiplantae where redox regulation of metabolic enzymes is important, it appears to be less common in diatoms, based on the current evidence, but regulation at the transcriptional level seems to be widespread. The role of post-translational modifications such as phosphorylation, glutathionylation, etc., and of protein-protein interactions, has been overlooked and should be investigated further. Diatoms and other chromalveolates are understudied compared to the Viridiplantae, especially given their ecological importance, but we believe that the ever-growing number of sequenced genomes combined with proteomics, metabolomics, enzyme measurements, and the application of novel techniques will provide a better understanding of how this important group of algae maintain their productivity under changing conditions.

trans-Zeatin-N-glucosides have biological activity in Arabidopsis thaliana

Cytokinin is an indispensable phytohormone responsible for physiological processes ranging from root development to leaf senescence. The term "cytokinin" refers to several dozen adenine-derived compounds occurring naturally in plants. Cytokinins (CKs) can be divided into various classes and forms; base forms are generally considered to be active while highly abundant cytokinin-N-glucosides (CKNGs), composed of a CK base irreversibly conjugated to a glucose molecule, are considered inactive. However, results from early CK studies suggest CKNGs do not always lack activity despite the perpetuation over several decades in the literature that they are inactive. Here we show that exogenous application of trans-Zeatin-N-glucosides (tZNGs, a specific class of CKNGs) to Arabidopsis results in CK response comparable to the application of an active CK base. These results are most apparent in senescence assays where both a CK base (tZ) and tZNGs (tZ7G, tZ9G) delay senescence in cotyledons. Further experiments involving root growth and shoot regeneration revealed tZNGs do not always have the same effects as tZ, and have largely distinct effects on the transcriptome and proteome. These data are in contrast to previous reports of CKNGs being inactive and raise questions about the function of these compounds as well as their mechanism of action.

Redox homeostasis in the growth zone of the rice leaf plays a key role in cold tolerance

We analysed the cellular and molecular changes in the leaf growth zone of tolerant and sensitive rice varieties in response to suboptimal temperatures. Cold reduced the final leaf length by 35% and 51% in tolerant and sensitive varieties, respectively. Tolerant lines exhibited a smaller reduction of the leaf elongation rate and greater compensation by an increased duration of leaf growth. Kinematic analysis showed that cold reduced cell production in the meristem and the expansion rate in the elongation zone, but the latter was compensated for by a doubling of the duration of cell expansion. We performed iTRAQ proteome analysis on proliferating and expanding parts of the leaf growth zone. We identified 559 and 542 proteins, of which 163 and 210 were differentially expressed between zones, and 96 and 68 between treatments, in the tolerant and sensitive lines, respectively. The categories protein biosynthesis and redox homeostasis were significantly overrepresented in the up-regulated proteins. We therefore measured redox metabolites and enzyme activities in the leaf growth zone, demonstrating that tolerance of rice lines to suboptimal temperatures correlates with the ability to up-regulate enzymatic antioxidants in the meristem and non-enzymatic antioxidants in the elongation zone.

Tissue-specific gene expression and protein abundance patterns are associated with fractionation bias in maize

BACKGROUND

Maize experienced a whole-genome duplication event approximately 5 to 12 million years ago. Because this event occurred after speciation from sorghum, the pre-duplication subgenomes can be partially reconstructed by mapping syntenic regions to the sorghum chromosomes. During evolution, maize has had uneven gene loss between each ancient subgenome. Fractionation and divergence between these genomes continue today, constantly changing genetic make-up and phenotypes and influencing agronomic traits.

RESULTS

Here we regenerate the subgenome reconstructions for the most recent maize reference genome assembly. Based on both expression and abundance data for homeologous gene pairs across multiple tissues, we observed functional divergence of genes across subgenomes. Although the genes in the larger maize subgenome are often expressing more highly than their homeologs in the smaller subgenome, we observed cases where homeolog expression dominance switches in different tissues. We demonstrate for the first time that protein abundances are higher in the larger subgenome, but they also show tissue-specific dominance, a pattern similar to RNA expression dominance. We also find that pollen expression is uniquely decoupled from protein abundance.

CONCLUSION

Our study shows that the larger subgenome has a greater range of functional assignments and that there is a relative lack of overlap between the subgenomes in terms of gene functions than would be suggested by similar patterns of gene expression and protein abundance. Our study also revealed that some reactions are catalyzed uniquely by the larger and smaller subgenomes. The tissue-specific, nonequivalent expression-level dominance pattern observed here implies a change in regulatory control which favors differentiated selective pressure on the retained duplicates leading to eventual change in gene functions.

2-DE-based proteomic analysis of protein changes associated with etiolated mesocotyl growth in Zea mays

Background

The mesocotyl connects the coleoptilar node and the basal part of the seminal root of maize (Zea mays) seedling. The mesocotyl pushes the shoot of the seedling out of the soil during seed germination; thus, its growth is highly related to deep-sowing tolerance. Although many studies on the maize mesocotyl have been carried out at physiological and molecular levels, the proteomic changes associated with cellular and physiological activities during mesocotyl growth are still unknown.

Results

In the present study, the maize hybrid Zhengdan 958 was used to study mesocotyl growth and accompanying protein changes. The dark-grown etiolated mesocotyls exhibited a slow-fast-slow feature, with significant changes in the levels of indole-3-acetic acid (IAA) and cellulose and the activity of peroxidase (POD). In particular, POD activity increased with mesocotyl growth, showing higher activity at the mature (lower) end of the mesocotyl. For the proteomic analysis, soluble proteins were extracted from etiolated mesocotyls dark-grown for 48 h, 84 h, and 132 h, corresponding to the initial, rapid, and slow growth periods, respectively, and subjected to separation by two-dimensional gel electrophoresis (2-DE). As a result, 88 differentially abundant proteins (DAPs) were identified using MALDI-TOF-TOF analysis. At 48 h, most DAPs were stress proteins, heat shock proteins and storage proteins; at 84 h, oxidation/reduction proteins, carbohydrate biogenesis-related proteins and cytoskeleton-related proteins were highly accumulated; at 132 h, the most striking DAPs were those involved in the synthesis and modification of the cell wall and the biogenesis of carbohydrates. Gene ontology (GO) analysis showed that changes in the abundance and proportion of DAPs were consistent with cellular and physiological activities and biological processes during mesocotyl growth. The accumulation of nine DAPs of interest was verified by immunoblotting and RT-qPCR.

Conclusions

The present study revealed that the protein patterns in 2-D gels differed greatly with mesocotyl growth. At different growth periods, a specific set of DAPs participate in various biological processes and underlie the cellular and physiological activities of the mesocotyl. These results contributed to the understanding of mesocotyl growth and the cultivation of maize lines with deep-sowing tolerance.

Investigation and development of maize fused network analysis with multi-omics

Maize is a critically important staple crop in the whole world, which has contributed to both economic security and food in planting areas. The main target for researchers and breeding is the improvement of maize quality and yield. The use of computational biology methods combined with multi-omics for selecting biomolecules of interest for maize breeding has been receiving more attention. Moreover, the rapid growth of high-throughput sequencing data provides the opportunity to explore biomolecules of interest at the molecular level in maize. Furthermore, we constructed weighted networks for each of the omics and then integrated them into a final fused weighted network based on a nonlinear combination method. We also analyzed the final fused network and mined the orphan nodes, some of which were shown to be transcription factors that played a key role in maize development. This study could help to improve maize production via insights at the multi-omics level and provide a new perspective for maize researchers. All related data have been released at http://lab.malab.cn/∼jj/maize.htm.

Extensive post-transcriptional buffering of gene expression in the response to severe oxidative stress in baker's yeast

Cells responds to diverse stimuli by changing the levels of specific effector proteins. These changes are usually examined using high throughput RNA sequencing data (RNA-Seq); transcriptional regulation is generally assumed to directly influence protein abundances. However, the correlation between RNA-Seq and proteomics data is in general quite limited owing to differences in protein stability and translational regulation. Here we perform RNA-Seq, ribosome profiling and proteomics analyses in baker's yeast cells grown in rich media and oxidative stress conditions to examine gene expression regulation at various levels. With the exception of a small set of genes involved in the maintenance of the redox state, which are regulated at the transcriptional level, modulation of protein expression is largely driven by changes in the relative ribosome density across conditions. The majority of shifts in mRNA abundance are compensated by changes in the opposite direction in the number of translating ribosomes and are predicted to result in no net change at the protein level. We also identify a subset of mRNAs which is likely to undergo specific translational repression during stress and which includes cell cycle control genes. The study suggests that post-transcriptional buffering of gene expression may be more common than previously anticipated.

Variable protein homeostasis in housekeeping and non-housekeeping pathways under mycotoxins stress

Transcript levels are the primary factor determining protein levels, but for the majority of genes, fold changes in transcript levels are larger than the corresponding changes in protein levels, a phenomenon that is termed “protein homeostasis”. However, this phenomenon is not well characterized in the context of environmental changes. In this study, we sequenced the entire transcriptome and proteome of chicken primary hepatocytes administered three mycotoxin treatments Aflatoxin B₁ (AFB₁), Ochoratoxin A (OTA) and Zearalenone (ZEN). Each mycotoxin induced unique set of differential expressed transcripts and proteins, suggesting variable cytotoxicity and biochemical action in cell. We found a weak positive correlation between transcript and protein changes, and the transcript changes were higher than the protein changes. Furthermore, we observed pathway-specific protein homeostasis pattern under mycotoxin stress. Specifically, the “Metabolism”, “Transcription” and “Translation” pathways (housekeeping pathways) showed lower fold changes in protein/mRNA levels than non-housekeeping pathways such as “Cell growth and death” and “Immune system”. Protein molecular weight had a weak negative effect on protein production, and this effect was stronger for non-housekeeping pathways. Overall, we hypothesize housekeeping pathways maintain stable protein production for baseline cellular function, whereas non-housekeeping pathways is associated with the fitness response to environmental stress.

Characterization of Proteome Variation During Modern Maize Breeding

The success of modern maize breeding has been demonstrated by remarkable increases in productivity with tremendous modification of agricultural phenotypes over the last century. Although the underlying genetic changes of the maize adaptation from tropical to temperate regions have been extensively studied, our knowledge is limited regarding the accordance of protein and mRNA expression levels accompanying such adaptation. Here we conducted an integrative analysis of proteomic and transcriptomic changes in a maize association panel. The minimum extent of correlation between protein and RNA levels suggests that variation in mRNA expression is often not indicative of protein expression at a population scale. This is corroborated by the observation that mRNA- and protein-based coexpression networks are relatively independent of each other, and many pQTLs arise without the presence of corresponding eQTLs. Importantly, compared with transcriptome, the subtypes categorized by the proteome show a markedly high accuracy to resemble the genomic subpopulation. These findings suggest that proteome evolved under a greater evolutionary constraint than transcriptome during maize adaptation from tropical to temperate regions. Overall, the integrated multi-omics analysis provides a functional context to interpret gene expression variation during modern maize breeding.

The controversies of silicon's role in plant biology

Contents Summary 67 I. Introduction 68 II. Silicon transport in plants: to absorb or not to absorb 69 III. The role of silicon in plants: not just a matter of semantics 71 IV. Silicon and biotic stress: beyond mechanical barriers and defense priming 76 V. Silicon and abiotic stress: a proliferation of proposed mechanisms 78 VI. The apoplastic obstruction hypothesis: a working model 79 VII. Perspectives and conclusions 80 Acknowledgements 81 References 81 SUMMARY: Silicon (Si) is not classified as an essential plant nutrient, and yet numerous reports have shown its beneficial effects in a variety of species and environmental circumstances. This has created much confusion in the scientific community with respect to its biological roles. Here, we link molecular and phenotypic data to better classify Si transport, and critically summarize the current state of understanding of the roles of Si in higher plants. We argue that much of the empirical evidence, in particular that derived from recent functional genomics, is at odds with many of the mechanistic assertions surrounding Si's role. In essence, these data do not support reports that Si affects a wide range of molecular-genetic, biochemical and physiological processes. A major reinterpretation of Si's role is therefore needed, which is critical to guide future studies and inform agricultural practice. We propose a working model, which we term the 'apoplastic obstruction hypothesis', which attempts to unify the various observations on Si's beneficial influences on plant growth and yield. This model argues for a fundamental role of Si as an extracellular prophylactic agent against biotic and abiotic stresses (as opposed to an active cellular agent), with important cascading effects on plant form and function.

Coregulation of ribosomal RNA with hundreds of genes contributes to phenotypic variation

Ribosomal repeats occupy 5% of a plant genome, yet there has been little study of their diversity in the modern age of genomics. Ribosomal copy number and expression variation present an opportunity to tap a novel source of diversity. In the present study, we estimated the ribosomal DNA (rDNA) copy number and ribosomal RNA (rRNA) expression for a population of maize inbred lines and investigated the potential role of rDNA and rRNA dosage in regulating global gene expression. Extensive variation was found in both ribosomal DNA copy number and ribosomal RNA expression among maize inbred lines. However, rRNA abundance was not consistent with the copy number of the rDNA. We have not found that the rDNA gene dosage has a regulatory role in gene expression; however, thousands of genes are identified to be coregulated with rRNA expression, including genes participating in ribosome biogenesis and other functionally relevant pathways. We further investigated the potential roles of copy number and the expression level of rDNA on agronomic traits and found that both correlated with flowering time but through different regulatory mechanisms. This comprehensive analysis suggested that rRNA expression variation is a valuable source of functional diversity that affects gene expression variation and field-based phenotypic changes.

Impact of drought and salt stress on the biosynthesis of alkaloids in Chelidonium majus L

When plants are exposed to various stress situations, their alkaloid concentration frequently is enhanced. This well-known phenomenon is presumably due to a passively enhanced rate of biosynthesis, caused by greatly elevated concentrations of NADPH in stressed plants. Here, we used Chelidonium majus L. plants, which accumulate high concentrations of dihydrocoptisine in their leaves, to study the impact of drought and salt stress on the biosynthesis and accumulation of alkaloids. In comparison to well-watered controls, in the transcriptome of the gene encoding the key enzyme in alkaloid biosynthesis, stylopine synthase, is enhanced in stressed C. majus plants. If we presuppose that increased transcript levels correlate with increased enzymatic activity of the gene products, these data indicate, for the first time, that stress-related increases in alkaloid concentration might not only be caused by the well-known stress-related passive shift, but may also be due to an enhancement of enzymatic capacity.

Post-transcriptional regulation of the oxidative stress response in plants

Due to their sessile lifestyle, plants can be exposed to several kinds of stresses that will increase the production of reactive oxygen species (ROS), such as hydrogen peroxide, singlet oxygen, and hydroxyl radicals, in the plant cells and activate several signaling pathways that cause alterations in the cellular metabolism. Nevertheless, when ROS production outreaches a certain level, oxidative damage to nucleic acids, lipids, metabolites, and proteins will occur, finally leading to cell death. Until now, the most comprehensive and detailed readout of oxidative stress responses is undoubtedly obtained at the transcriptome level. However, transcript levels often do not correlate with the corresponding protein levels. Indeed, together with transcriptional regulations, post-transcriptional, translational, and/or post-translational regulations will shape the active proteome. Here, we review the current knowledge on the post-transcriptional gene regulation during the oxidative stress responses in planta.

Stress responsive mitochondrial proteins in Arabidopsis thaliana

In the last decade plant mitochondria have emerged as a target, sensor and initiator of signalling cascades to a variety of stress and adverse growth conditions. A combination of various 'omic profiling approaches combined with forward and reverse genetic studies have defined how mitochondria respond to stress and the signalling pathways and regulators of these responses. Reactive oxygen species (ROS)-dependent and -independent pathways, specific metabolites, complex I dysfunction, and the mitochondrial unfolded protein response (UPR) pathway have been proposed to date. These pathways are regulated by kinases (sucrose non-fermenting response like kinase; cyclin dependent protein kinase E 1) and transcription factors from the abscisic acid-related, WRKY and NAC families. A number of independent studies have revealed that these mitochondrial signalling pathways interact with a variety of phytohormone signalling pathways. While this represents significant progress in the last decade there are more pathways to be uncovered. Post-transcriptional/translational regulation is also a likely determinant of the mitochondrial stress response. Unbiased analyses of the expression of genes encoding mitochondrial proteins in a variety of stress conditions reveal a modular network exerting a high degree of anterograde control. As abiotic and biotic stresses have significant impact on the yield of important crops such as rice, wheat and barley we will give an outlook of how knowledge gained in Arabidopsis may help to increase crop production and how emerging technologies may contribute.