Quantitative proteomics and transcriptomics reveal key metabolic processes associated with cotton fiber initiation

Transcriptome, Ectopic Expression and Genetic Population Analysis Identify Candidate Genes for Fiber Quality Improvement in Cotton

Comparative transcriptome analysis of fiber tissues between Gossypium barbadense and Gossypium hirsutum could reveal the molecular mechanisms underlying high-quality fiber formation and identify candidate genes for fiber quality improvement. In this study, 759 genes were found to be strongly upregulated at the elongation stage in G. barbadense, which showed four distinct expression patterns (I-IV). Among them, the 346 genes of group IV stood out in terms of the potential to promote fiber elongation, in which we finally identified 42 elongation-related candidate genes by comparative transcriptome analysis between G. barbadense and G. hirsutum. Subsequently, we overexpressed GbAAR3 and GbTWS1, two of the 42 candidate genes, in Arabidopsis plants and validated their roles in promoting cell elongation. At the secondary cell wall (SCW) biosynthesis stage, 2275 genes were upregulated and exhibited five different expression profiles (I-V) in G. barbadense. We highlighted the critical roles of the 647 genes of group IV in SCW biosynthesis and further picked out 48 SCW biosynthesis-related candidate genes by comparative transcriptome analysis. SNP molecular markers were then successfully developed to distinguish the SCW biosynthesis-related candidate genes from their G. hirsutum orthologs, and the genotyping and phenotyping of a BC₃F₅ population proved their potential in improving fiber strength and micronaire. Our results contribute to the better understanding of the fiber quality differences between G. barbadense and G. hirsutum and provide novel alternative genes for fiber quality improvement.

Transcriptomic analysis reveals the key role of histone deacetylation via mediating different phytohormone signalings in fiber initiation of cotton

Background

Histone deacetylation is one of the most important epigenetic modifications and plays diverse roles in plant development. However, the detailed functions and mechanisms of histone deacetylation in fiber development of cotton are still unclear. HDAC inhibitors (HDACi) have been commonly used to study the molecular mechanism underlying histone deacetylation or to facilitate disease therapy in humans through hindering the histone deacetylase catalytic activity. Trichostatin A (TSA)—the most widely used HDACi has been extensively employed to determine the role of histone deacetylation on different developmental stages of plants.

Results

Through in vitro culture of ovules, we observed that exogenous application of TSA was able to inhibit the fiber initiation development. Subsequently, we performed a transcriptomic analysis to reveal the underlying mechanisms. The data showed that TSA treatment resulted in 4209 differentially expressed genes, which were mostly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, and carbon metabolism pathways. The phytohormone signal transduction pathways harbor the most differentially expressed genes. Deeper studies showed that some genes promoting auxin, Gibberellic Acid (GA) signaling were down-regulated, while some genes facilitating Abscisic Acid (ABA) and inhibiting Jasmonic Acid (JA) signaling were up-regulated after the TSA treatments. Further analysis of plant hormone contents proved that TSA significantly promoted the accumulation of ABA, JA and GA₃.

Conclusions

Collectively, histone deacetylation can regulate some key genes involved in different phytohormone pathways, and consequently promoting the auxin, GA, and JA signaling, whereas repressing the ABA synthesis and signaling to improve the fiber cell initiation. Moreover, the genes associated with energy metabolism, phenylpropanoid, and glutathione metabolism were also regulated by histone deacetylation. The above results provided novel clues to illuminate the underlying mechanisms of epigenetic modifications as well as related different phytohormones in fiber cell differentiation, which is also very valuable for the molecular breeding of higher quality cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00840-4.

Proteome and physiological analyses reveal tobacco (Nicotiana tabacum) peroxidase 7 (POD 7) functions in responses to copper stress

Copper is a micronutrient essential for plant growth and development. However, Cu is also a heavy metal element that has deleterious impacts on plants when excessively accumulated in the environment. To understand the molecular mechanism underlying tobacco in response to Cu stress, iTRAQ based technology was used to identify differentially expressed proteins (DEPs) and important metabolic pathways in tobacco plants treated with excessive CuSO4. The results showed that 180 DEPs were detected between the treatment and control, among which 78 were upregulated and 102 were downregulated. These DEPs can be functionally divided into 65 categories and are closely related to metabolic pathways, carbon metabolism, secondary metabolite biosynthesis, biosynthesis of antibiotics, glyoxylate and dicarboxylate metabolism, and glycolysis/gluconeogenesis. Peroxidase7 was significantly upregulated and was selected and overexpressed in tobacco. Then, positive transgenic lines and wild type plants were exposed to a Cu stress environment. The results showed that Peroxidase7 transgenic tobacco plants exhibited enhanced Cu stress resistance with decreased malondialdehyde and Cu contents, and increased shoot dry weight, root length, secondary root number, SOD, POD and CAT activity. The present study suggests that the ROS scavenging mechanism is essential for tobacco plants in response to Cu stress and that Peroxidase7 functions in tobacco plant resistance to excessive Cu environment.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

Identification and Structure Analysis of KCS Family Genes Suggest Their Reponding to Regulate Fiber Development in Long-Staple Cotton Under Salt-Alkaline Stress

Plant 3-ketoacyl-CoA synthase (KCS) gene family catalyzed a β ketoacyl-CoA synthase, which was the rate-limiting enzyme for the synthesis of very long chain fatty acids (VLCFAs). Gossypium barbadense was well-known not only for high-quality fiber, which was perceived as a cultivated species of Gossypium. In this study, a total of 131 KCS genes were identified in four cotton species, there were 38, 44, 26, 23 KCS genes in the G. barbadense, the G. hirsutum, the G. arboreum and G. raimondii, respectively. The gene structure and expression pattern were analyzed. GBKCS genes were divided into six subgroups, the chromosome distribution of members of the family were mapped. The prediction of cis-acting elements of the GBKCS gene promoters suggested that the GBKCS genes may be involved in hormone signaling, defense and the stress response. Collinearity analysis on the KCS genes of the four cotton species were formulated. Tandem duplication played an indispensable role in the evolution of the KCS gene family. Specific expression analysis of 20 GBKCS genes indicated that GBKCS gene were widely expressed in the first 25 days of fiber development. Among them, GBKCS3, GBKCS8, GBKCS20, GBKCS34 were expressed at a high level in the initial long-term level of the G. barbadense fiber. This study established a foundation to further understanding of the evolution of KCS genes and analyze the function of GBKCS genes.

RNA-Seq with a novel glabrous-ZM24fl reveals some key lncRNAs and the associated targets in fiber initiation of cotton

BACKGROUND

Cotton fiber is an important natural resource for textile industry and an excellent model for cell biology study. Application of glabrous mutant cotton and high-throughput sequencing facilitates the identification of key genes and pathways for fiber development and cell differentiation and elongation. LncRNA is a type of ncRNA with more than 200 nt in length and functions in the ways of chromatin modification, transcriptional and post-transcriptional modification, and so on. However, the detailed lncRNA and associated mechanisms for fiber initiation are still unclear in cotton.

RESULTS

In this study, we used a novel glabrous mutant ZM24fl, which is endowed with higher somatic embryogenesis, and functions as an ideal receptor for cotton genetic transformation. Combined with the high-throughput sequencing, fatty acid pathway and some transcription factors such as MYB, ERF and bHLH families were identified the important roles in fiber initiation; furthermore, 3,288 lncRNAs were identified, and some differentially expressed lncRNAs were also analyzed. From the comparisons of ZM24_0 DPA vs ZM24_-2 DPA and fl_0 DPA vs ZM24_0 DPA, one common lncRNA MSTRG 2723.1 was found that function upstream of fatty acid metabolism, MBY25-mediating pathway, and pectin metabolism to regulate fiber initiation. In addition, other lncRNAs MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1 were also showed potential important roles in fiber development; and the co-expression analysis between lncRNAs and targets showed the distinct models of different lncRNAs and complicated interaction between lncRNAs in fiber development of cotton.

CONCLUSIONS

From the above results, a key lncRNA MSTRG 2723.1 was identified that might mediate some key genes transcription of fatty acid metabolism, MYB25-mediating pathway, and pectin metabolism to regulate fiber initiation of ZM24 cultivar. Co-expression analysis implied that some other important lncRNAs (e.g., MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1) were also showed the different regulatory model and interaction between them, which proposes some valuable clues for the lncRNAs associated mechanisms in fiber development.

Integrative transcriptome and proteome analyses provide deep insights into the molecular mechanism of salt tolerance in Limonium bicolor

Integrative transcriptome and proteome analyses revealed many candidate members that may involve in salt secretion from salt glands in Limonium bicolor. Limonium bicolor, a typical recretohalophyte, protects itself from salt damage by excreting excess salt out of its cells through salt glands. Here, to provide an overview of the salt-tolerance mechanism of L. bicolor, we conducted integrative transcriptome and proteome analyses of this species under salt treatment. We identified numerous differentially expressed transcripts and proteins that may be related to the salt-tolerance mechanism of L. bicolor. By measuring the Na⁺ secretion rate, were found that this cation secretion rate of a single salt gland was significantly increased after high salinity treatment compared with that in control and then reached the maximum in a short time. Interestingly, transcripts and proteins involved in transmembrane transport of ions were differentially expressed in response to high salinity treatment, suggesting a number of genes and proteins they may play important roles in the salt-stress response. Correlation between differentially expressed transcript and protein profiles revealed several transcripts and proteins that may be responsible for salt tolerance, such as cellulose synthases and annexins. Our findings uncovered many candidate transcripts and proteins in response to the salt tolerance of L. bicolor, providing deep insights into the molecular mechanisms of this important process in recretohalophytes.

Comparative Metabolomics Analysis Reveals Sterols and Sphingolipids Play a Role in Cotton Fiber Cell Initiation

Cotton fiber is a seed trichome that protrudes from the outer epidermis of cotton ovule on the day of anthesis (0 day past anthesis, 0 DPA). The initial number and timing of fiber cells are closely related to fiber yield and quality. However, the mechanism underlying fiber initiation is still unclear. Here, we detected and compared the contents and compositions of sphingolipids and sterols in 0 DPA ovules of Xuzhou142 lintless-fuzzless mutants (Xufl) and Xinxiangxiaoji lintless-fuzzless mutants (Xinfl) and upland cotton wild-type Xuzhou142 (XuFL). Nine classes of sphingolipids and sixty-six sphingolipid molecular species were detected in wild-type and mutants. Compared with the wild type, the contents of Sphingosine-1-phosphate (S1P), Sphingosine (Sph), Glucosylceramide (GluCer), and Glycosyl-inositol-phospho-ceramides (GIPC) were decreased in the mutants, while the contents of Ceramide (Cer) were increased. Detail, the contents of two Cer molecular species, d18:1/22:0 and d18:1/24:0, and two Phyto-Cer molecular species, t18:0/22:0 and t18:0/h22:1 were significantly increased, while the contents of all GluCer and GIPC molecular species were decreased. Consistent with this result, the expression levels of seven genes involved in GluCer and GIPC synthesis were decreased in the mutants. Furthermore, exogenous application of a specific inhibitor of GluCer synthase, PDMP (1-phenyl-2-decanoylamino-3-morpholino-1-propanol), in ovule culture system, significantly inhibited the initiation of cotton fiber cells. In addition, five sterols and four sterol esters were detected in wild-type and mutant ovules. Compared with the wild type, the contents of total sterol were not significantly changed. While the contents of stigmasterol and campesterol were significantly increased, the contents of cholesterol were significantly decreased, and the contents of total sterol esters were significantly increased. In particular, the contents of campesterol esters and stigmasterol esters increased significantly in the two mutants. Consistently, the expression levels of some sterol synthase genes and sterol ester synthase genes were also changed in the two mutants. These results suggested that sphingolipids and sterols might have some roles in the initiation of fiber cells. Our results provided a novel insight into the regulatory mechanism of fiber cell initiation.

Integrative iTRAQ-based proteomic and transcriptomic analysis reveals the accumulation patterns of key metabolites associated with oil quality during seed ripening of Camellia oleifera

Camellia oleifera (C. oleifera) is one of the four major woody oil-bearing crops in the world and has relatively high ecological, economic, and medicinal value. Its seeds undergo a series of complex physiological and biochemical changes during ripening, which is mainly manifested as the accumulation and transformation of certain metabolites closely related to oil quality, especially flavonoids and fatty acids. To obtain new insights into the underlying molecular mechanisms, a parallel analysis of the transcriptome and proteome profiles of C. oleifera seeds at different maturity levels was conducted using RNA sequencing (RNA-seq) and isobaric tags for relative and absolute quantification (iTRAQ) complemented with gas chromatography-mass spectrometry (GC-MS) data. A total of 16,530 transcripts and 1228 proteins were recognized with significant differential abundances in pairwise comparisons of samples at various developmental stages. Among these, 317 were coexpressed with a poor correlation, and most were involved in metabolic processes, including fatty acid metabolism, α-linolenic acid metabolism, and glutathione metabolism. In addition, the content of total flavonoids decreased gradually with seed maturity, and the levels of fatty acids generally peaked at the fat accumulation stage; these results basically agreed with the regulation patterns of genes or proteins in the corresponding pathways. The expression levels of proteins annotated as upstream candidates of phenylalanine ammonia-lyase (PAL) and chalcone synthase (CHS) as well as their cognate transcripts were positively correlated with the variation in the flavonoid content, while shikimate O-hydroxycinnamoyltransferase (HCT)-encoding genes had the opposite pattern. The increase in the abundance of proteins and mRNAs corresponding to alcohol dehydrogenase (ADH) was associated with a reduction in linoleic acid synthesis. Using weighted gene coexpression network analysis (WGCNA), we further identified six unique modules related to flavonoid, oil, and fatty acid anabolism that contained hub genes or proteins similar to transcription factors (TFs), such as MADS intervening keratin-like and C-terminal (MIKC_MADS), type-B authentic response regulator (ARR-B), and basic helix-loop-helix (bHLH). Finally, based on the known metabolic pathways and WGCNA combined with the correlation analysis, five coexpressed transcripts and proteins composed of cinnamyl-alcohol dehydrogenases (CADs), caffeic acid 3-O-methyltransferase (COMT), flavonol synthase (FLS), and 4-coumarate: CoA ligase (4CL) were screened out. With this exploratory multiomics dataset, our results presented a dynamic picture regarding the maturation process of C. oleifera seeds on Hainan Island, not only revealing the temporal specific expression of key candidate genes and proteins but also providing a scientific basis for the genetic improvement of this tree species.

Integrated transcriptome and proteome analyses identify novel regulatory network of nucleus pulposus cells in intervertebral disc degeneration

Background

Degeneration of intervertebral disc is a major cause of lower back pain and neck pain. Studies have tried to unveil the regulatory network using either transcriptomic or proteomic analysis. However, neither have fully elucidated the exact mechanism of degeneration process. Since post-transcriptional regulation may affect gene expression by modulating the translational process of mRNA to protein product, a combined transcriptomic and proteomic study may provide more insight into the key regulatory network of Intervertebral disc degeneration.

Methods

In order to obtain the proteomic and transcriptomic data, we performed label-free proteome analysis on freshly isolated nucleus pulposus cells and obtained transcriptome profiling data from the Gene Expression Omnibus repository. To identify the key regulatory network of intervertebral disc degeneration in nucleus pulposus cells, we performed bioinformatic analyses and established a protein-RNA interacting network. To validate the candidate genes, we performed in vitro experimentation and immunochemistry labeling to identify their potential function during nucleus pulposus degeneration.

Results

The label-free proteome analysis identified altogether 656 proteins, and 503 of which were differentially expressed between nucleus pulposus cells from degenerated or normal disc cells. Using the existing nucleus pulposus transcriptomic profiling data, we integrated the proteomic and transcriptomic data of nucleus pulposus cells, and established a protein-RNA interacting network to show the combined regulatory network of intervertebral disc degeneration. In the network, we found 9 genes showed significant changes, and 6 of which (CHI3L1, KRT19, COL6A2, DPT, TNFAIP6 and COL11A2) showed concordant changes in both protein and mRNA level. Further functional analysis showed these candidates can significantly affect the degeneration of the nucleus pulposus cell when altering their expression.

Conclusions

This study is the first to use combined analysis of proteomic and transcriptomic profiling data to identify novel regulatory network of nucleus pulposus cells in intervertebral disc degeneration. Our established protein-RNA interacting network demonstrated novel regulatory mechanisms and key genes that may play vital roles in the pathogenesis of intervertebral disc degeneration.

Physiological and quantitative proteomic analysis of NtPRX63-overexpressing tobacco plants revealed that NtPRX63 functions in plant salt resistance

High salinity is harmful to crop yield and productivity. Peroxidases (PRXs) play crucial roles in H2O2 scavenging. In our previous study, PRX63 significantly upregulated in tobacco plants under salt stress. Thus, in order to understand the function of PRX63 in tobacco salt response, we overexpressed this gene in tobacco (Nicotiana tabacum L.), investigated the morphological, physiological and proteomic profiles of NtPRX63-overexpressing tobacco transgenic lines and wild type. The results showed that, compared with the wild type, the transgenic tobacco plants presented enhanced salt tolerance and displayed lower ROS (reactive oxygen species), malondialdehyde (MDA) and Na+ contents; higher biomass, potassium content, soluble sugar content, and peroxidase activity; and higher expression levels of NtSOD, NtPOD and NtCAT. Protein abundance analysis revealed 123 differentially expressed proteins between the transgenic and wild-type plants. These proteins were functionally classified into 18 categories and are involved in 41 metabolic pathways. Furthermore, among the 123 proteins, eight proteins involved in the ROS-scavenging system, 12 involved in photosynthesis and energy metabolism processes, two stress response proteins, one signal transduction protein and one disulfide isomerase were significantly upregulated. Furthermore, three novel proteins that may be involved in the plant salt response were also identified. The results of our study indicate that an enhanced ROS-scavenging ability, together with the expression of proteins related to energy mobilization and the stress response, functions in the confirmed salt resistance of transgenic tobacco plants. Our data provide valuable information for research on the function of NtPRX63 in tobacco in response to abiotic stress.

Progress of ethylene action mechanism and its application on plant type formation in crops

The plant hormone ethylene exerts a huge influence in the whole life cycle of plants, especially stress-resistance responses. With the development of functional genomics, that the action mechanism of ethylene takes part in mediated plant architecture has been clarified gradually, such as plant roots, stems, leaves, fiber elongation and so on. Accordingly, the application of ethylene on crops chemical control and genetic improvement is greatly expanded. From the view of ethylene mediated plant architecture in crops, here reviewed advances in ethylene biosynthesis and signal transduction pathway, stress-resistance responses and the yield potential enhance of crops in recently 20 years. On these grounds, the objectives of this paper were to provide scientific reference and a useful clue for the crop creation of ideal plant type.

Comparative Proteomic Analysis of Molecular Differences between Leaves of Wild-Type Upland Cotton and Its Fuzzless-Lintless Mutant

Fuzzless-lintless mutant (fl) ovules of upland cotton have been used to investigate cotton fiber development for decades. However, the molecular differences of green tissues between fl and wild-type (WT) cotton were barely reported. Here, we found that gossypol content, the most important secondary metabolite of cotton leaves, was higher in Gossypium hirsutum L. cv Xuzhou-142 (Xu142) WT than in fl. Then, we performed comparative proteomic analysis of the leaves from Xu142 WT and its fl. A total of 4506 proteins were identified, of which 103 and 164 appeared to be WT- and fl-specific, respectively. In the 4239 common-expressed proteins, 80 and 74 were preferentially accumulated in WT and fl, respectively. Pathway enrichment analysis and protein–protein interaction network analysis of both variety-specific and differential abundant proteins showed that secondary metabolism and chloroplast-related pathways were significantly enriched. Quantitative real-time PCR confirmed that the expression levels of 12 out of 16 selected genes from representative pathways were consistent with their protein accumulation patterns. Further analyses showed that the content of chlorophyll a in WT, but not chlorophyll b, was significantly increased compared to fl. This work provides the leaf proteome profiles of Xu142 and its fl mutant, indicating the necessity of further investigation of molecular differences between WT and fl leaves.

Transcriptome profiling of Gossypium arboreum during fiber initiation and the genome-wide identification of trihelix transcription factors

Cotton fiber initiation is the first step in fiber development, and it determines the yield. Here, genome-wide transcriptome profiling of Gossypium arboreum was performed to determine the molecular basis of cotton fiber initiation. A comparison of the transcriptomes of fiber-bearing ovules at -0.5, 0, 0.5, 1, 1.5, 2, 2.5 and 3 d post-anthesis detected 12,049 differentially expressed genes that mainly participated in ribosome, carbon metabolism and amino acid biosynthesis pathways. Genes encoding alcohol dehydrogenase 1 and hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase, involving in fatty acid degradation and flavonoid biosynthesis, were enriched. Furthermore, 1049 differentially expressed transcription factors were identified. Among these, 17 were trihelix family transcription factors, which play important roles in plant development and responses to biotic and abiotic stresses. In total, 52 full-length trihelix genes, named as GaGTs, were identified in G. arboreum and located in 12 of the 13 cotton chromosomes. Transcriptomic data and a quantitative real-time PCR analysis indicated that several GaGTs were significantly induced during fiber initiation in G. arboreum. Thus, the genome-wide comprehensive analysis of gene expression in G. arboreum fiber initiation will serve as a useful resource for unraveling the functions of specific genes. The phylogenetic relationships and expression analyses of the G. arboreum trihelix genes established a solid foundation for future comprehensive functional analyses of the GaGTs.

Cytokinin inhibits cotton fiber initiation by disrupting PIN3a-mediated asymmetric accumulation of auxin in the ovule epidermis

Auxin-dependent cell expansion is crucial for initiation of fiber cells in cotton (Gossypium hirsutum), which ultimately determines fiber yield and quality. However, the regulation of this process is far from being well understood. In this study, we demonstrate an antagonistic effect between cytokinin (CK) and auxin on cotton fiber initiation. In vitro and in planta experiments indicate that enhanced CK levels can reduce auxin accumulation in the ovule integument, which may account for the defects in the fiberless mutant xu142fl. In turn, supplementation with auxin can recover fiber growth of CK-treated ovules and mutant ovules. We further found that GhPIN3a is a key auxin transporter for fiber-cell initiation and is polarly localized to the plasma membranes of non-fiber cells, but not to those of fiber cells. This polar localization allows auxin to be transported within the ovule integument while specifically accumulating in fiber cells. We show that CKs antagonize the promotive effect of auxin on fiber cell initiation by undermining asymmetric accumulation of auxin in the ovule epidermis through down-regulation of GhPIN3a and disturbance of the polar localization of the protein.

An improved protein extraction method applied to cotton leaves is compatible with 2-DE and LC-MS

Background

Two-dimensional electrophoresis (2-DE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) are widely used in plant proteomics research. However, these two techniques cannot be simultaneously satisfied by traditional protein extraction methods when investigate cotton leaf proteome.

Results

Here, we evaluated the efficiency of three different protein extraction methods for 2-DE and LC-MS/MS analyses of total proteins obtained from cotton leaves. The protein yield of the borax/PVPP/phenol (BPP) method (0.14%) was significantly lower than the yields of the trichloroacetic acid/acetone (TCA) precipitation method (1.42%) and optimized TCA combined with BPP (TCA-B) method (0.47%). The BPP method was failed to get a clear 2-DE electrophoretogram. Fifty pairs of protein spots were randomly selected from the 2-DE gels of TCA- and TCA-B-extracted proteins for identification by MALDI TOF/TOF, and the results of 42 pairs were consistent. High-throughput proteomic analysis showed that 6339, 9282 and 9697 unique proteins were identified from the total cotton leaf proteins extracted by the TCA, BPP and TCA-B methods, respectively. Gene Ontology (GO) analysis revealed that the proteins specifically identified by TCA method were primarily distributed in the plasma membrane, while BPP and TCA-B methods specific proteins distributed in the cytosol, indicating the sub-cellular preference of different protein extraction methods. Further, ATP-dependent zinc metalloprotease FTSH 8 could be observed in the 2-DE gels of TCA and TCA-B methods, and could only be detected in the LC-MS/MS results of the BPP and TCA-B methods, showing that TCA-B method might be the optimized choice for both 2-DE and LC-MS/MS.

Conclusion

Our data provided an improved TCA-B method for protein extraction that is compatible with 2-DE and LC-MS/MS for cotton leaves and similar plant tissues which is rich in polysaccharides and polyphenols.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5658-5) contains supplementary material, which is available to authorized users.

Quantitative phosphoproteomics of lectin receptor-like kinase VI.4 dependent abscisic acid response in Arabidopsis thaliana

Lectin receptor-like kinases (LecRKs) play important roles in the responses to adverse environment stress. Abscisic acid (ABA) is a plant hormone involved in plant growth, development and adverse environmental stress responses. Although some studies of ABA response LecRK genes have been reported, the molecular mechanisms of LecRKs regulation of downstream pathways under ABA induction are not well understood. The present study showed that LecRK-VI.4 responded to ABA and negatively regulated stomatal closure. Here, a quantitative phosphoproteomics approach based on mass spectrometry was employed to study the roles of LecRK-VI.4 in the ABA signaling pathway. Metal oxide affinity beads and C18 chromatography were used for phosphopeptide enrichment and separation. The isobaric tags for relative and absolute quantitation were used for profiling the phosphoproteome of mutant lecrk-vi.4-1 and wild-type Col-0 Arabidopsis under normal growth conditions or ABA treatments. In total, 475 unique phosphopeptides were quantified, including 81 phosphopeptides related to LecRK-VI.4 regulation. Gene ontology, protein-protein interaction and motif analysis were performed. The bioinformatics data showed that phosphorylated proteins regulated by LecRK-VI.4 had close relations with factors of stomatal function, which included aquaporin activity, H⁺ pump activity and the Ca²⁺ concentration in the cytoplasm. These data have expanded our understanding of how LecRK-VI.4 regulates ABA-mediated stomatal movements.

Physiological and quantitative proteomic analyses unraveling potassium deficiency stress response in alligator weed (Alternanthera philoxeroides L.) root

Physiological and iTRAQ based proteomic analysis provided new insights into potassium deficiency stress response in alligator weed root. Alligator weed (Alternanthera philoxeroides) has a strong ability to adapt to potassium deficiency (LK) stress. Proteomic changes in response to this stress are largely unknown in alligator weed. In this study, we investigated physiological and molecular mechanisms under LK using isobaric tags for relative and absolute quantitation to characterize proteome-level changes in this plant. First, root physiology, 2, 3, 5-Triphenyl-trazolium chloride (TTC) assay and peroxidase activity were significantly altered after 10 and 15 days of LK treatment. The comparative proteomic analysis suggested a total of 375 proteins were differential abundance proteins. The proteomic results were verified by western blot assays and quantitative real-time PCR. Correlation analysis of transcription and proteomics suggested protein processing in the endoplasmic reticulum, endocytosis, and spliceosome pathways were significantly enriched. The protein responsible for energy metabolism, signal sensing and transduction and protein degradation played crucial roles in this stress. Twelve ubiquitin pathway related proteins were identified in our study, among them 11 proteins were up-regulated. All protein ubiquitination of lysine using pan antibodies were also increased after LK treatment. Our study provide a valuable insights of molecular mechanism underlying LK stress response in alligator weed roots and afford a vital basis to further study potassium nutrition molecular breeding of other plant species.

Plant Systems Biology at the Single-Cell Level

Our understanding of plant biology is increasingly being built upon studies using 'omics and system biology approaches performed at the level of the entire plant, organ, or tissue. Although these approaches open new avenues to better understand plant biology, they suffer from the cellular complexity of the analyzed sample. Recent methodological advances now allow plant scientists to overcome this limitation and enable biological analyses of single-cells or single-cell-types. Coupled with the development of bioinformatics and functional genomics resources, these studies provide opportunities for high-resolution systems analyses of plant phenomena. In this review, we describe the recent advances, current challenges, and future directions in exploring the biology of single-cells and single-cell-types to enhance our understanding of plant biology as a system.

Proteomic analysis of phytase transgenic and non-transgenic maize seeds

Proteomics has become a powerful technique for investigating unintended effects in genetically modified crops. In this study, we performed a comparative proteomics of the seeds of phytase-transgenic (PT) and non-transgenic (NT) maize using 2-DE and iTRAQ techniques. A total of 148 differentially expressed proteins (DEPs), including 106 down-regulated and 42 up-regulated proteins in PT, were identified. Of these proteins, 32 were identified through 2-DE and 116 were generated by iTRAQ. It is noteworthy that only three proteins could be detected via both iTRAQ and 2-DE, and most of the identified DEPs were not newly produced proteins but proteins with altered abundance. These results indicated that many DEPs could be detected in the proteome of PT maize seeds and the corresponding wild type after overexpression of the target gene, but the changes in these proteins were not substantial. Functional classification revealed many DEPs involved in posttranscriptional modifications and some ribosomal proteins and heat-shock proteins that may generate adaptive effects in response to the insertion of exogenous genes. Protein-protein interaction analysis demonstrated that the detected interacting proteins were mainly ribosomal proteins and heat-shock proteins. Our data provided new information on such unintended effects through a proteomic analysis of maize seeds.

Comparative proteomic analysis of autotetraploid and diploid Paulownia tomentosa reveals proteins associated with superior photosynthetic characteristics and stress adaptability in autotetraploid Paulownia

To enlarge the germplasm resource of Paulownia plants, we used colchicine to induce autotetraploid Paulownia tomentosa, as reported previously. Compared with its diploid progenitor, autotetraploid P. tomentosa exhibits better photosynthetic characteristics and higher stress resistance. However, the underlying mechanism for its predominant characteristics has not been determined at the proteome level. In this study, isobaric tag for relative and absolute quantitation coupled with liquid chromatography-tandem mass spectrometry was employed to compare proteomic changes between autotetraploid and diploid P. tomentosa. A total of 1427 proteins were identified in our study, of which 130 proteins were differentially expressed between autotetraploid and diploid P. tomentosa. Functional analysis of differentially expressed proteins revealed that photosynthesis-related proteins and stress-responsive proteins were significantly enriched among the differentially expressed proteins, suggesting they may be responsible for the photosynthetic characteristics and stress adaptability of autotetraploid P. tomentosa. The correlation analysis between transcriptome and proteome data revealed that only 15 (11.5%) of the differentially expressed proteins had corresponding differentially expressed unigenes between diploid and autotetraploid P. tomentosa. These results indicated that there was a limited correlation between the differentially expressed proteins and the previously reported differentially expressed unigenes. This work provides new clues to better understand the superior traits in autotetraploid P. tomentosa and lays a theoretical foundation for developing Paulownia breeding strategies in the future.

iTRAQ-based proteomic analysis reveals key proteins affecting muscle growth and lipid deposition in pigs

Growth rate and meat quality, two economically important traits in pigs, are controlled by multiple genes and biological pathways. In the present study, we performed a proteomic analysis of longissimus dorsi muscle from six-month-old pigs from two Chinese native mini-type breeds (TP and DSP) and two introduced western breeds (YY and LL) using isobaric tag for relative and absolute quantification (iTRAQ). In total, 4,815 peptides corresponding to 969 proteins were detected. Comparison of expression patterns between TP-DSP and YY-LL revealed 288 differentially expressed proteins (DEPs), of which 169 were up-regulated and 119 were down-regulated. Functional annotation suggested that 28 DEPs were related to muscle growth and 15 to lipid deposition. Protein interaction network predictions indicated that differences in muscle growth and muscle fibre between TP-DSP and YY-LL groups were regulated by ALDOC, ENO3, PGK1, PGK2, TNNT1, TNNT3, TPM1, TPM2, TPM3, MYL3, MYH4, and TNNC2, whereas differences in lipid deposition ability were regulated by LPL, APOA1, APOC3, ACADM, FABP3, ACADVL, ACAA2, ACAT1, HADH, and PECI. Twelve DEPs were analysed using parallel reaction monitoring to confirm the reliability of the iTRAQ analysis. Our findings provide new insights into key proteins involved in muscle growth and lipid deposition in the pig.

Two-Dimensional Gel Electrophoresis-Based Proteomic Analysis Reveals N-terminal Truncation of the Hsc70 Protein in Cotton Fibers In Vivo

On two-dimensional electrophoresis gels, six protein spots from cotton ovules and fibers were identified as heat shock cognate 70 kD protein (Hsc70). Three spots corresponded to an experimental molecular weight (MW) of 70 kD (spots 1, 2 and 3), and the remaining three spots corresponded to an experimental MW slightly greater than 45 kD (spots 4, 5 and 6). Protein spots 1, 2 and 3 were abundant on gels of 0-day (the day of anthesis) wild-type (WT) ovules, 0-day fuzzless-lintless mutant ovules and 10-day WT ovules but absent from gels of 10-day WT fibers. Three individual transcripts encoding these six protein spots were obtained by using rapid amplification of cDNA ends (RACE). Edman degradation and western blotting confirmed that the three 45 kD Hsc70 protein spots had the same N-terminal, which started from the T271 amino acid in the intact Hsc70 protein. Furthermore, quadrupole time-of-flight mass spectrometry analysis identified a methylation modification on the arginine at position 475 for protein spots 4 and 5. Our data demonstrate that site-specific in vivo N-terminal truncation of the Hsc70 protein was particularly prevalent in cotton fibers, indicating that post-translational regulation might play an important role in cotton fiber development.

Comparative proteomic analyses of Asian cotton ovules with attached fibers in the early stages of fiber elongation process

BACKGROUND

Plenty of proteomic studies were performed to characterize the allotetraploid upland cotton fiber elongation process, whereas little is known about the elongating diploid cotton fiber proteome.

METHODS

In this study, we used a two-dimensional electrophoresis-based comparative proteomic approach to profile dynamic proteomes of diploid Asian cotton ovules with attached fibers in the early stages of fiber elongation process. One-way ANOVA and Student-Newman-Keuls test were used to find the differentially displayed protein (DDP) spots.

RESULTS

A total of 55 protein spots were found having different abundance ranging from 1 to 9 days post-anthesis (DPA) in a two-day interval. These 55 DDP spots were all successfully identified using high-resolution mass spectrometric analyses. Gene ontology analyses revealed that proteoforms involved in energy/carbohydrate metabolism, redox homeostasis, and protein metabolism are the most abundant. In addition, orthologues of the 13 DDP spots were also found in differential proteome of allotetraploid elongating cotton fibers, suggesting their possible essential roles in fiber elongation process.

CONCLUSIONS

Our results not only revealed the dynamic proteome change of diploid Asian cotton fiber and ovule during early stages of fiber elongation process but also provided valuable resource for future studies on the molecular mechanism how the polyploidization improves the trait of fiber length.

Ovarian proteomic study reveals the possible molecular mechanism for hyperprolificacy of Small Tail Han sheep

Small Tail Han sheep is a widely bred farm animal in China which has attracted lots of attention due to their high prolificacy and year-round estrus. However, the molecular mechanism of its fecundity remains unrevealed. The FecB gene polymorphism has been found to be associated with the ovulation rate and litter size of sheep. In the present study, we constructed an iTRAQ-based quantitative proteomics analysis to compare the ovarian proteomes of FecB+FecB+ genotype Small Tail Han sheep ewes (Han ++), FecB^BFecB^B Han ewes (Han BB) and Dorset ewes (Dorset). Hundreds of differentially expressed proteins between each two groups were identified; GO and KEGG pathway analysis indicated that the expressions of those proteins involved in ribosome assembly, protein translation and mTOR pathway between Dorset and both Han groups were highly different. Between Han ++ and Han BB groups, higher level of protein expressions were related to mitochondrial oxidation functions such as oxidoreductase activity, cytochrome-c oxidase activity and electron carrier activity. This was identified in Han BB group, which may contribute to the elevated ovulation rate of Han BB ewes. In conclusion, our work provided a prospective understanding of the molecular mechanism for high prolificacy of Small Tail Han sheep.

Polyploidy and the proteome

Although major advances have been made during the past 20 years in our understanding of the genetic and genomic consequences of polyploidy, our knowledge of polyploidy and the proteome is in its infancy. One of our goals is to stimulate additional study, particularly broad-scale proteomic analyses of polyploids and their progenitors. Although it may be too early to generalize regarding the extent to which transcriptomic data are predictive of the proteome of polyploids, it is clear that the proteome does not always reflect the transcriptome. Despite limited data, important observations on the proteomes of polyploids are emerging. In some cases, proteomic profiles show qualitatively and/or quantitatively non-additive patterns, and proteomic novelty has been observed. Allopolyploids generally combine the parental contributions, but there is evidence of parental dominance of one contributing genome in some allopolyploids. Autopolyploids are typically qualitatively identical to but quantitatively different from their parents. There is also evidence of parental legacy at the proteomic level. Proteomes clearly provide insights into the consequences of genomic merger and doubling beyond what is obtained from genomic and/or transcriptomic data. Translating proteomic changes in polyploids to differences in morphology and physiology remains the holy grail of polyploidy--this daunting task of linking genotype to proteome to phenotype should emerge as a focus of polyploidy research in the next decade. This article is part of a Special Issue entitled: Plant Proteomics--a bridge between fundamental processes and crop production, edited by Dr. Hans-Peter Mock.

Two-dimensional gel electrophoresis-based analysis provides global insights into the cotton ovule and fiber proteomes

Proteomic analysis of upland cotton was performed to profile the global detectable proteomes of ovules and fibers using two-dimensional electrophoresis (2DE). A total of 1,203 independent protein spots were collected from representative 2DE gels, which were digested with trypsin and identified by matrix-assisted laser desorption and ionization-time-offlight/ time-of-flight (MALDI-TOF/TOF) mass spectrometry. The mass spectrometry or tandem mass spectrometry (MS or MS/MS) data were then searched against a local database constructed from Gossypium hirsutum genome sequences, resulting in successful identification of 975 protein spots (411 for ovules and 564 for fibers). Functional annotation analysis of the 975 identified proteins revealed that ovule-specific proteins were mainly enriched in functions related to fatty acid elongation, sulfur amino acid metabolism and post-replication repair, while fiber-specific proteins were enriched in functions related to root hair elongation, galactose metabolism and D-xylose metabolic processes. Further annotation analysis of the most abundant protein spots showed that 28.96% of the total proteins in the ovule were mainly located in the Golgi apparatus, endoplasmic reticulum, mitochondrion and ribosome, whereas in fibers, 27.02% of the total proteins were located in the cytoskeleton, nuclear envelope and cell wall. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses of the ovule-specific protein spots P61, P93 and P198 and fiber-specific protein spots 230, 477 and 511 were performed to validate the proteomics data. Protein-protein interaction network analyses revealed very different network cluster patterns between ovules and fibers. This work provides the largest protein identification dataset of 2DE-detectable proteins in cotton ovules and fibers and indicates potentially important roles of tissue-specific proteins, thus providing insights into the cotton ovule and fiber proteomes on a global scale.

Quantitative Proteomic Analysis Reveals Molecular Adaptations in the Hippocampal Synaptic Active Zone of Chronic Mild Stress-Unsusceptible Rats

BACKGROUND

While stressful events are recognized as an important cause of major depressive disorder, some individuals exposed to life stressors maintain normal psychological functioning. The molecular mechanism(s) underlying this phenomenon remain unclear. Abnormal transmission and plasticity of hippocampal synapses have been implied to play a key role in the pathoetiology of major depressive disorder.

METHODS

A chronic mild stress protocol was applied to separate susceptible and unsusceptible rat subpopulations. Proteomic analysis using an isobaric tag for relative and absolute quantitation coupled with tandem mass spectrometry was performed to identify differential proteins in enriched hippocampal synaptic junction preparations.

RESULTS

A total of 4318 proteins were quantified, and 89 membrane proteins were present in differential amounts. Of these, SynaptomeDB identified 81 (91%) having a synapse-specific localization. The unbiased profiles identified several candidate proteins within the synaptic junction that may be associated with stress vulnerability or insusceptibility. Subsequent functional categorization revealed that protein systems particularly involved in membrane trafficking at the synaptic active zone exhibited a positive strain as potential molecular adaptations in the unsusceptible rats. Moreover, through STRING and immunoblotting analysis, membrane-associated GTP-bound Rab3a and Munc18-1 appear to coregulate syntaxin-1/SNAP25/VAMP2 assembly at the hippocampal presynaptic active zone of unsusceptible rats, facilitating SNARE-mediated membrane fusion and neurotransmitter release, and may be part of a stress-protection mechanism in actively maintaining an emotional homeostasis.

CONCLUSIONS

The present results support the concept that there is a range of potential protein adaptations in the hippocampal synaptic active zone of unsusceptible rats, revealing new investigative targets that may contribute to a better understanding of stress insusceptibility.

The Proteomic Response of Arabidopsis thaliana to Cadmium Sulfide Quantum Dots, and Its Correlation with the Transcriptomic Response

A fuller understanding of the interaction between plants and engineered nanomaterials is of topical relevance because the latter are beginning to find applications in agriculture and the food industry. There is a growing need to establish objective safety criteria for their use. The recognition of two independent Arabidopsis thaliana mutants displaying a greater level of tolerance than the wild type plant to exposure to cadmium sulfide quantum dots (CdS QDs) has offered the opportunity to characterize the tolerance response at the physiological, transcriptomic, and proteomic levels. Here, a proteomics-based comparison confirmed the conclusions drawn from an earlier transcriptomic analysis that the two mutants responded to CdS QD exposure differently both to the wild type and to each other. Just over half of the proteomic changes mirrored documented changes at the level of gene transcription, but a substantial number of transcript/gene product pairs were altered in the opposite direction. An interpretation of the discrepancies is given, along with some considerations regarding the use and significance of -omics when monitoring the potential toxicity of ENMs for health and environment.