Integrative transcriptome, proteome, phosphoproteome and genetic mapping reveals new aspects in a fiberless mutant of cotton

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 BC3F5 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.

Transcriptome profiling for developmental stages Protaetia brevitarsis seulensis with focus on wing development and metamorphosis

A white-spotted flower chafer Protaetia brevitarsis seulensis widely distributed in Asian countries is traditionally used in oriental medicine. This study explored gene expression abundance with respect to wing development and metamorphosis in P. b. seulensis based on the large-scale RNA-seq data. The transcriptome assembly consists of 23,551 high-quality transcripts which are approximately 96.7% covered. We found 265 wing development genes, 19 metamorphosis genes, and 1,314 candidates. Of the 1,598 genes, 1,594 are included exclusively in cluster 4 with similar gene co-expression patterns. The network centrality analyses showed that wing development- and metamorphosis-related genes have a high degree of betweenness centrality and are expressed most highly in eggs, moderately in pupa and adults, and lowest in larva. This study provides some meaningful clues for elucidating the genetic modulation mechanism of wing development and metamorphosis in P. b. seulensis.

Integration of transcriptomic and proteomic analyses reveals several levels of metabolic regulation in the excess starch and early senescent leaf mutant lses1 in rice

BACKGROUND

The normal metabolism of transitory starch in leaves plays an important role in ensuring photosynthesis, delaying senescence and maintaining high yield in crops. OsCKI1 (casein kinase I1) plays crucial regulatory roles in multiple important physiological processes, including root development, hormonal signaling and low temperature-treatment adaptive growth in rice; however, its potential role in regulating temporary starch metabolism or premature leaf senescence remains unclear. To reveal the molecular regulatory mechanism of OsCKI1 in rice leaves, physiological, transcriptomic and proteomic analyses of leaves of osckI1 allele mutant lses1 (leaf starch excess and senescence 1) and its wild-type varieties (WT) were performed.

RESULTS

Phenotypic identification and physiological measurements showed that the lses1 mutant exhibited starch excess in the leaves and an obvious leaf tip withering phenotype as well as high ROS and MDA contents, low chlorophyll content and protective enzyme activities compared to WT. The correlation analyses between protein and mRNA abundance are weak or limited. However, the changes of several important genes related to carbohydrate metabolism and apoptosis at the mRNA and protein levels were consistent. The protein-protein interaction (PPI) network might play accessory roles in promoting premature senescence of lses1 leaves. Comprehensive transcriptomic and proteomic analysis indicated that multiple key genes/proteins related to starch and sugar metabolism, apoptosis and ABA signaling exhibited significant differential expression. Abnormal increase in temporary starch was highly correlated with the expression of starch biosynthesis-related genes, which might be the main factor that causes premature leaf senescence and changes in multiple metabolic levels in leaves of lses1. In addition, four proteins associated with ABA accumulation and signaling, and three CKI potential target proteins related to starch biosynthesis were up-regulated in the lses1 mutant, suggesting that LSES1 may affect temporary starch accumulation and premature leaf senescence through phosphorylation crosstalk ABA signaling and starch anabolic pathways.

CONCLUSION

The current study established the high correlation between the changes in physiological characteristics and mRNA and protein expression profiles in lses1 leaves, and emphasized the positive effect of excessive starch on accelerating premature leaf senescence. The expression patterns of genes/proteins related to starch biosynthesis and ABA signaling were analyzed via transcriptomes and proteomes, which provided a novel direction and research basis for the subsequent exploration of the regulation mechanism of temporary starch and apoptosis via LSES1/OsCKI1 in rice.

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.

Comparative Dynamic Transcriptome Reveals the Delayed Secondary-Cell-Wall Thickening Results in Altered Lint Percentage and Fiber Elongation in a Chromosomal Segment Substitution Line of Cotton (Gossypium hirsutum L.)

Lint percentage (LP) is an important yield component in cotton that is usually affected by initial fiber number and cell wall thickness. To explore how fiber cell wall development affects LP, phenotypic identification and dynamic transcriptome analysis were conducted using a single segment substitution line of chromosome 15 (SL15) that harbors a major quantitative trait locus (QTL) for LP. Compared to its recurrent parent LMY22, SL15 did not differ in initial fiber number, but the fiber cell wall thickness and single-fiber weight decreased significantly, altering LP. The comparative transcriptome profiles revealed that the secondary cell wall (SCW) development phase of SL15 was relatively delayed. Meanwhile, the expression of genes related to cell expansion decreased more slightly in SL15 with fiber development, resulting in relatively higher expression at SL15_25D than at LMY22_25D. SCW development-related genes, such as GhNACs and GhMYBs, in the putative NAC-MYB-CESA network differentially expressed at SL15_25D, along with the lower expression of CESA6, CSLC12, and CSLA2. The substituted chromosomal interval was further investigated, and found 6 of 146 candidate genes were differentially expressed in all four cell development periods including 10, 15, 20 and 25 DPA. Genetic variation and co-expression analysis showed that GH_D01G0052, GH_D01G0099, GH_D01G0100, and GH_D01G0140 may be important candidate genes associated with qLP-C15-1. Our results provide novel insights into cell wall development and its relationship with LP, which is beneficial for lint yield and fiber quality improvement.

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.

Phosphorylation-mediated signalling in flowering: prospects and retrospects of phosphoproteomics in crops

Protein phosphorylation is a major post-translational modification, regulating protein function, stability, and subcellular localization. To date, annotated phosphorylation data are available mainly for model organisms and humans, despite the economic importance of crop species and their large kinomes. Our understanding of the phospho-regulation of flowering in relation to the biology and interaction between the pollen and pistil is still significantly lagging, limiting our knowledge on kinase signalling and its potential applications to crop production. To address this gap, we bring together relevant literature that were previously disconnected to present an overview of the roles of phosphoproteomic signalling pathways in modulating molecular and cellular regulation within specific tissues at different morphological stages of flowering. This review is intended to stimulate research, with the potential to increase crop productivity by providing a platform for novel molecular tools.

Understanding the role of phytohormones in cotton fiber development through omic approaches; recent advances and future directions

Cotton is among the most important fiber crops for the textile-based industry, thanks to its cellulose-rich mature fibers. The fiber initiation and elongation are one of the best models for deciphering mechanisms of single-cell differentiation and growth, that also target of fiber development programs. During the last couple of decades, high yielding omics approaches (genomics, transcriptomics, and proteomics), have helped in the identification of several genes and gene products involved in fiber development along with functional relationship to phytohormones. For example, MYB transcription factor family and Sus gene family have been evidenced by controlling cotton fiber initiation. Most importantly, the biosynthesis, responses, and transporting of phytohormones is documented to participate in the initiation of cotton fibers. Herein, in this review, the reliable genetic evidence by manipulating the above genes in cotton have been summarized to describe the relationships among key phytohormones, transcription factors, proteins, and downstream fiber growth-related genes such as Sus. The effect of other important factors such as ROS, fatty acid metabolism, and actin (globular multi-functional proteins) over fiber development has also been discussed. The challenges and deficiencies in the research of cotton fiber development have been mentioned along with a future perspective to discover new crucial genes using multiple omics analysis.

Comparative acetylome analysis of wild-type and fuzzless-lintless mutant ovules of upland cotton (Gossypium hirsutum Cv. Xu142) unveils differential protein acetylation may regulate fiber development

Protein acetylation (KAC) is a significant post-translational modification, which plays an essential role in the regulation of growth and development. Unfortunately, related studies are inadequately available in angiosperms, and to date, there is no report providing insight on the role of protein acetylation in cotton fiber development. Therefore, we first compared the lysine-acetylation proteome (acetylome) of upland cotton ovules in the early fiber development stages by using wild-type as well as its fuzzless-lintless mutant to identify the role of KAC in the fiber development. A total of 1696 proteins with 2754 acetylation sites identified with the different levels of acetylation belonging to separate subcellular compartments suggesting a large number of proteins differentially acetylated in two cotton cultivars. About 80% of the sites were predicted to localize in the cytoplasm, chloroplast, and mitochondria. Seventeen significantly enriched acetylation motifs were identified. Serine and threonine and cysteine located downstream and upstream to KAC sites. KEGG pathway enrichment analysis indicated oxidative phosphorylation, fatty acid, ribosome and protein, and folate biosynthesis pathways enriched significantly. To our knowledge, this is the first report of comparative acetylome analysis to compare the wild-type as well as its fuzzless-lintless mutant acetylome data to identify the differentially acetylated proteins, which may play a significant role in cotton fiber development.

Complementary iTRAQ Proteomic and Transcriptomic Analyses of Leaves in Tea Plant ( Camellia sinensis L.) with Different Maturity and Regulatory Network of Flavonoid Biosynthesis

The quality of tea is highly related with the maturity of the fresh tea leaves at harvest. The present study investigated the proteomic and transcriptomic profiles of tea leaves with different maturity, using iTRAQ and RNA-seq technologies. A total of 4455 proteins and 27 930 unigenes were identified, with functional enrichment analyses of GO categorization and KEGG annotation. The compositions of flavonoids (catechins and flavonols) in tea leaves were determined. The total content of flavonoids decreased with leaf maturity, in accordance with the protein regulation patterns of shikimate, phenylpropanoid, and flavonoid pathways. The abundance of ANR had a positive correlation with epi-catechin content, while LAR abundance was positively related with catechin content ( P < 0.05). The biosynthetic network of flavonoid biosynthesis was discussed in combination with photosynthesis, primary metabolism, and transcription factors. Bud had the lowest activities of photosynthesis and carbon fixation but the highest flavonoid biosynthesis ability in opposite to mature leaf. SUS-INV switch might be an important joint for carbon flow shifting into the follow-up biochemical syntheses. This work provided a comprehensive overview on the functional protein profile changes of tea leaves at different growing stages and also proposed a research direction regarding the correlations between primary metabolism and flavonoid biosynthesis.

Role of GhHDA5 in H3K9 deacetylation and fiber initiation in Gossypium hirsutum

Cotton fibers are single-celled trichomes that initiate from the epidermal cells of the ovules at or before anthesis. Here, we identified that the histone deacetylase (HDAC) activity is essential for proper cotton fiber initiation. We further identified 15 HDACs homoeologs in each of the A- and D-subgenomes of Gossypium hirsutum. Few of these HDAC homoeologs expressed preferentially during the early stages of fiber development [-1, 0 and 6 days post-anthesis (DPA)]. Among them, GhHDA5 expressed significantly at the time of fiber initiation (-1 and 0 DPA). The in vitro assay for HDAC activity indicated that GhHDA5 primarily deacetylates H3K9 acetylation marks. Moreover, the reduced expression of GhHDA5 also suppresses fiber initiation and lint yield in the RNA interference (RNAi) lines. The 0 DPA ovules of GhHDA5RNAi lines also showed alterations in reactive oxygen species homeostasis and elevated autophagic cell death in the developing fibers. The differentially expressed genes (DEGs) identified through RNA-seq of RNAi line (DEP12) and their pathway analysis showed that GhHDA5 modulates expression of many stress and development-related genes involved in fiber development. The reduced expression of GhHDA5 in the RNAi lines also resulted in H3K9 hyper-acetylation on the promoter region of few DEGs assessed by chromatin immunoprecipitation assay. The positively co-expressed genes with GhHDA5 showed cumulative higher expression during fiber initiation, and gene ontology annotation suggests their involvement in fiber development. Furthermore, the predicted protein interaction network in the positively co-expressed genes indicates HDA5 modulates fiber initiation-specific gene expression through a complex involving reported repressors.

Integration of proteomic and transcriptomic profiles reveals multiple levels of genetic regulation of salt tolerance in cotton

BACKGROUND

Salinity is a major abiotic stress that limits upland cotton growth and reduces fibre production worldwide. To reveal genetic regulation via transcript and protein levels after salt stress, we comprehensively analysed the global changes in mRNA, miRNA, and protein profiles in response to salt stress in two contrasting salt-tolerant cotton genotypes.

RESULTS

In the current study, proteomic and mRNA-seq data were combined to reveal that some genes are differentially expressed at both the proteomic and mRNA levels. However, we observed no significant change in mRNA corresponding to most of the strongly differentially abundant proteins. This finding may have resulted from global changes in alternative splicing events and miRNA levels under salt stress conditions. Evidence was provided indicating that several salt stress-responsive proteins can alter miRNAs and modulate alternative splicing events in upland cotton. The results of the stringent screening of the mRNA-seq and proteomic data between the salt-tolerant and salt-sensitive genotypes identified 63 and 85 candidate genes/proteins related to salt tolerance after 4 and 24 h of salt stress, respectively, between the tolerant and sensitive genotype. Finally, we predicted an interaction network comprising 158 genes/proteins and then discovered that two main clusters in the network were composed of ATP synthase (CotAD_74681) and cytochrome oxidase (CotAD_46197) in mitochondria. The results revealed that mitochondria, as important organelles involved in energy metabolism, play an essential role in the synthesis of resistance proteins during the process of salt exposure.

CONCLUSION

We provided a plausible schematic for the systematic salt tolerance model; this schematic reveals multiple levels of gene regulation in response to salt stress in cotton and provides a list of salt tolerance-related genes/proteins. The information here will facilitate candidate gene discovery and molecular marker development for salt tolerance breeding in cotton.

Proteomic Analysis of Differences in Fiber Development between Wild and Cultivated Gossypium hirsutum L

Upland cotton (Gossypium hirsutum L.) is one of the world's most important fiber crops, accounting for more than 90% of all cotton production. While their wild progenitors have relatively short and coarse, often tan-colored fibers, modern cotton cultivars possess longer, finer, stronger, and whiter fiber. In this study, the wild and cultivated cottons (YU-3 and TM-1) selected show significant differences on fibers at 10 days postanthesis (DPA), 20 DPA, and mature stages at the morphological level. To explore the effects of domestication, reveal molecular mechanisms underlying these phenotypic differences, and better inform our efforts to further enhance cotton fiber quality, isobaric tags for relative and absolute protein quantification-facilitated proteomic methods were performed on developing fibers. There were 6990 proteins identified; among them, 336 were defined as differentially expressed proteins between fibers of wild versus domesticated cotton. The down- or up-regulated proteins in wild cotton were involved in phenylpropanoid biosynthesis, zeatin biosynthesis, fatty acid elongation, and other processes. Association analysis between transcriptome and proteome showed positive correlations between transcripts and proteins at both 10 DPA and 20 DPA. Differences in proteomics have been verified at the mRNA level by quantitative real-time polymerase chain reaction and have been validated at the physiological and biochemical levels by POD (peroxidase) activity assays and ZA (zeatin) content estimates. This work corroborates the major pathways involved in cotton fiber development and demonstrates that POD activity and zeatin content have a great potential related to fiber elongation and thickening.

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.