Comparative proteomic analysis reveals the mechanisms governing cotton fiber differentiation and initiation

GhMAC3e is involved in plant growth and defense response to Verticillium dahliae

KEY MESSAGE

GhMAC3e expression was induced by various stresses and hormones. GhMAC3e may regulate plant growth by influencing auxin distribution, and play important roles in Verticillium wilt resistance via mediating SA signaling. The MOS4-Associated Complex (MAC) is a highly conserved protein complex involved in pre-mRNA splicing and spliceosome assembly, which plays a vital role in plant immunity. It comprises key components such as MOS4, CDC5, and PRL1. MAC3A/B, as U-box E3 ubiquitin ligases, are crucial for various plant processes including development, stress responses, and disease resistance. However, their roles in cotton remain largely unknown. In this study, we first cloned the GhMAC3e gene from cotton and explored its biological functions by using virus-induced gene silencing (VIGS) in cotton and transgenic overexpression in Arabidopsis. The results showed that GhMAC3e is ubiquitously expressed in cotton tissues and could be induced by salt stress, Verticillium dahliae (VD) infection, PEG, ABA, ETH, GA3, MeJA, and SA. Silencing GhMAC3e retarded primary stem growth and reduced biomass of cotton coupled with the reduced auxin content in the petioles and veins. Silencing GhMAC3e up-regulated expression of cell growth-related genes GhXTH16 and Gh3.6, while down-regulated GhSAUR12 expression. Ectopic expression of GhMAC3e in Arabidopsis significantly enhanced its resistance to Verticillium wilt (VW) in terms of decreased pathogen biomass and lowered plant mortality. Overexpression of GhMAC3e dramatically upregulated AtPR1 by around 15 fold and more than 262 fold under basal and VD inoculation condition, respectively. This change was not associated with the expression of GhNPR1. In conclusion, GhMAC3e may not only regulate plant growth by influencing auxin distribution and growth-related gene expression, but also play important roles in VW resistance via mediating SA signaling independent of NPR1 transcription level.

Cell type-specific cytonuclear coevolution in three allopolyploid plant species

Cytonuclear disruption may accompany allopolyploid evolution as a consequence of the merger of different nuclear genomes in a cellular environment having only one set of progenitor organellar genomes. One path to reconcile potential cytonuclear mismatch is biased expression for maternal gene duplicates (homoeologs) encoding proteins that target to plastids and/or mitochondria. Assessment of this transcriptional form of cytonuclear coevolution at the level of individual cells or cell types remains unexplored. Using single-cell (sc-) and single-nucleus (sn-) RNAseq data from eight tissues in three allopolyploid species, we characterized cell type-specific variations of cytonuclear coevolutionary homoeologous expression and demonstrated the temporal dynamics of expression patterns across development stages during cotton fiber development. Our results provide unique insights into transcriptional cytonuclear coevolution in plant allopolyploids at the single-cell level.

Molecular Regulation of Cotton Fiber Development: A Review

Cotton (Gossypium spp.) is an economically important natural fiber crop. The quality of cotton fiber has a substantial effect on the quality of cotton textiles. The identification of cotton fiber development-related genes and exploration of their biological functions will not only enhance our understanding of the elongation and developmental mechanisms of cotton fibers but also provide insights that could aid the cultivation of new cotton varieties with improved fiber quality. Cotton fibers are single cells that have been differentiated from the ovule epidermis and serve as a model system for research on single-cell differentiation, growth, and fiber production. Genes and fiber formation mechanisms are examined in this review to shed new light on how important phytohormones, transcription factors, proteins, and genes linked to fiber development work together. Plant hormones, which occur in low quantities, play a critically important role in regulating cotton fiber development. Here, we review recent research that has greatly contributed to our understanding of the roles of different phytohormones in fiber development and regulation. We discuss the mechanisms by which phytohormones regulate the initiation and elongation of fiber cells in cotton, as well as the identification of genes involved in hormone biosynthetic and signaling pathways that regulate the initiation, elongation, and development of cotton fibers.

Identification of key proteins related to high-quality fiber in Upland cotton via proteomics analysis

The dynamics of cotton fiber elongation and microfibirl deposition orientation were delineated; advanced ethylene synthesis and redox reaction homeostasis may be crucial for high-quality fiber formation. Fiber length, strength, and fineness determine the use and commercial value of cotton fiber, but their underlying molecular mechanisms remain obscure. We compared the dynamic change trajectories of length, diameter and microfibril orientation angle of the fibers produced by an introgression line SY6167 which generates high-quality fibers even better than Sea island cotton with those of the common-quality fibers from TM-1 across 5 to 30 days post anthesis (DPA). The proteomes were profiled and compared at six representative time points using 2-DE and MS/MS. 14 proteins differentially expressed inside each of cotton line temporally and significantly different tween the two lines were identified. The dynamic change trajectories of fiber length and microfibril angle are close to "s" and reverse "s" growth curves, respectively. SY6167 and TM-1 fibers entered the logarithmic elongation phase simultaneously at 10 DPA, and SY6167 kept elongating logarithmically for 2 more days than TM-1. In parallel to logarithmic elongation, microfibril orientation angles dived sharply, and SY6167 declined faster for a shorter duration than TM-1. 53% of the identified proteins are related to redox homeostasis, and most of them are expressed at higher levels in SY6167 during logarithmic elongation. 1-Aminocyclopropane-1-Carboxylic Acid Oxidase (ACO) started to accumulate at 16 DPA in SY6167, and its encoding genes were highly expressed at this stage, with a much higher level than TM-1. These findings suggest high-quality fibers are associated with high expression of the proteins related to stress and redox homeostasis, the continuously elevated expression of ethylene synthesis ACO gene may play an essential role.

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.

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.

GaHD1, a candidate gene for the Gossypium arboreum SMA-4 mutant, promotes trichome and fiber initiation by cellular H2O2 and Ca2+ signals

Cotton fibers are initiated from the epidermal cells of the ovule before or on the day of anthesis. Gossypium arboreum SMA-4 mutant contains recessive mutation (sma-4(ha)) and has the phenotypes of fibreless seeds and glabrous stems. In this study, fine mapping and alternative splicing analysis indicated a nucleotide substitution (AG → AC) at splicing site in a homeodomain-leucine zipper IV family gene (GaHD1) might cause gene A3S (Alternative 3' splicing) mistake, suggested that GaHD1 was the candidate gene of sma-4(ha). Many genes related to the fiber initiation are identified to be differentially expressed in the mutant which could result in the blocked fiber initiation signals such as H₂O_2, or Ca in the mutant. Further comparative physiological analysis of H₂O₂ production and Ca²⁺ flux in the SMA-4 and wide type cotton confirmed that H₂O₂ and Ca were important fiber initiation signals and regulated by GaHD1. The in vitro ovule culture of the mutant with hormones recovered the fibered phenotype coupled with the restoration of these signals. Overexpressing of GaHD1 in Arabidopsis increased trichome densities on the sepal, leaf, and stem tissues while transient silencing of the GaHD1 gene in G. arboreum reduced the trichome densities. These phenotypes indicated that GaHD1 is the candidate gene of SMA-4 with a crucial role in acting upstream molecular switch of signal transductions for cotton trichome and fiber initiations.

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.

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 analysis reveals differences in the mechanisms of fiber initiation and elongation between long- and short-fiber cotton (Gossypium hirsutum L.) lines

Background

Improving the yield and fiber quality of upland cotton is a goal of plant breeders. However, increasing the yield and quality of cotton fibers is becoming more urgent. While the growing human population needs more cotton fiber, climate change is reducing the amount of land on which cotton can be planted, or making it difficult to ensure that water and other resources will be available in optimal quantities. The most logical means of improving yield and quality is understanding and manipulating the genes involved. Here, we used comparative transcriptomics to explore differences in gene expression between long- and short-fiber cotton lines to identify candidate genes useful for cotton improvement.

Results

Light and electron microscopy revealed that the initial fiber density was significantly greater in our short-fiber group (SFG) than in our long-fiber group (LFG). Compared with the SFG fibers, the LFG fibers were longer at all developmental stages. Comparison of the LFG and SFG transcriptomes revealed a total of 3538 differentially expressed genes (DEGs). Notably, at all three developmental stages examined, two expression patterns, consistently downregulated (profile 0) and consistently upregulated (profile 7), were identified, and both were significantly enriched in the SFG and LFG. Twenty-two DEGs known to be involved in fiber initiation were detected in profile 0, while 31 DEGs involved in fiber elongation were detected in profile 7. Functional annotation suggested that these DEGs, which included ERF1, TUA2, TUB1, and PER64, affect fiber elongation by participating in the ethylene response, microtubule synthesis, and/or the peroxidase (POD) catalytic pathway. qRT-PCR was used to confirm the RNA sequencing results for select genes.

Conclusions

A comparison of SFG and LFG transcription profiles revealed modest but important differences in gene expression between the groups. Notably, our results confirm those of previous studies suggesting that genes involved in ethylene, tubulin, and POD pathways play important roles in fiber development. The 22 consistently downregulated DEGs involved in fiber initiation and the 31 consistently upregulated genes involved in fiber elongation are seemingly good candidate genes for improving fiber initiation and elongation in cotton.

Electronic supplementary material

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

Fiber Quality Improvement in Upland Cotton (Gossypium hirsutum L.): Quantitative Trait Loci Mapping and Marker Assisted Selection Application

Genetic improvement in fiber quality is one of the main challenges for cotton breeders. Fiber quality traits are controlled by multiple genes and are classified as complex quantitative traits, with a negative relationship with yield potential, so the genetic gain is low in traditional genetic improvement by phenotypic selection. The availability of Gossypium genomic sequences facilitates the development of high-throughput molecular markers, quantitative trait loci (QTL) fine mapping and gene identification, which helps us to validate candidate genes and to use marker assisted selection (MAS) on fiber quality in breeding programs. Based on developments of high density linkage maps, QTLs fine mapping, marker selection and omics, we have performed trait dissection on fiber quality traits in diverse populations of upland cotton. QTL mapping combined with multi-omics approaches such as, RNA sequencing datasets to identify differentially expressed genes have benefited the improvement of fiber quality. In this review, we discuss the application of molecular markers, QTL mapping and MAS for fiber quality improvement in upland cotton.

PIN-formed protein, a door to reveal the mechanism for auxin-triggered initiation of cotton fiber

Cotton fibers are differentiated ovule epidermal cells that provide an ideal model to study cell differentiation and elongation. Establishment of auxin maximum in fiber cells is crucial for cotton-fiber protrusion from ovule surface. However, it is unclear where the auxin originates from and how the auxin accumulates in fiber cells. Our recent results indicate that the auxin is mainly imported from the outside of ovules, and transported to fiber cells through GhPIN (homolog of PIN-formed proteins in cotton) -mediated polar auxin transport, rather than in situ synthesis. Based on our finding in GhPINs, we discuss here briefly how auxin flow to fiber cells and auxin gradient in ovule epidermis is established mainly by GhPIN3a protein.

Proteomic analysis of common bean stem under drought stress using in-gel stable isotope labeling

Drought is an abiotic stress that strongly influences plant growth, development and productivity. Proteome changes in the stem of the drought-tolerant common bean (Phaseolus vulgaris L.) cultivar Tiber have were when the plants were exposed to drought. Five-week-old plants were subjected to water deficit by withholding irrigation for 7, 12 and 17days, whereas control plants were regularly irrigated. Relative water content (RWC) of leaves, as an indicator of the degree of cell and tissue hydration, showed the highest statistically significant differences between control and drought-stressed plants after 17days of treatment, where RWC remained at 90% for control and declined to 45% for stressed plants. Plants exposed to drought for 17days and control plants at the same developmental stage were included in quantitative proteomic analysis using in-gel stable isotope labeling of proteins in combination with mass spectrometry. The quantified proteins were grouped into several functional groups, mainly into energy metabolism, photosynthesis, proteolysis, protein synthesis and proteins related to defense and stress. 70kDa heat shock protein showed the greatest increase in abundance under drought of all the proteins, suggesting its role in protecting plants against stress by re-establishing normal protein conformations and thus cellular homeostasis. The abundance of proteins involved in protein synthesis also increased under drought stress, important for recovery of damaged proteins involved in the plant cell's metabolic activities. Other important proteins in this study were related to proteolysis and folding, which are necessary for maintaining proper cellular protein homeostasis. Taken together, these results reveal the complexity of pathways involved in the drought stress response in common bean stems and enable comparison with the results of proteomic analysis of leaves, thus providing important information to further understand the biochemical and molecular mechanisms of drought response in this important legume.

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.

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

To investigate the molecular mechanisms of fiber initiation in cotton (Gossypium spp.), an integrated approach combining transcriptome, iTRAQ-based proteome and genetic mapping was taken to compare the ovules of the Xuzhou 142 wild type (WT) with its fuzzless-lintless (fl) mutant at −3 and 0 day post-anthesis. A total of 1,953 mRNAs, 187 proteins, and 131 phosphoproteins were differentially expressed (DE) between WT and fl, and the levels of transcripts and their encoded proteins and phosphoproteins were highly congruent. A functional analysis suggested that the abundance of proteins were mainly involved in amino sugar, nucleotide sugar and fatty acid metabolism, one carbon pool for folate metabolism and flavonoid biosynthesis. qRT-PCR, Western blotting, and enzymatic assays were performed to confirm the regulation of these transcripts and proteins. A molecular mapping located the lintless gene li3 in the fl mutant on chromosome 26 for the first time. A further in-silico physical mapping of DE genes with sequence variations between fl and WT identified one and four candidate genes in the li3 and n2 regions, respectively. Taken together, the transcript abundance, phosphorylation status of proteins at the fiber initiation stage and candidate genes have provided insights into regulatory processes underlying cotton fiber initiation.

Proteomics of terpenoid biosynthesis and secretion in trichomes of higher plant species

Among the specialized (secondary) plant metabolites, terpenoids represent the most diverse family and are often involved in the defense against pathogens and herbivores. Terpenoids can be produced both constitutively and in response to the environment. At the front line of this defense strategy are the glandular trichomes, which are organs dedicated primarily to the production of specialized metabolites. Mass spectrometry-based proteomics is a powerful tool, which is very useful to investigate enzymes involved in metabolic pathways, such as the synthesis and secretion of terpenoids in glandular trichomes. Here we review the strategies used to investigate the specific roles of these particular organs from non-model plant species, mainly belonging to the Lamiaceae, Solanaceae, and Cannabaceae families. We discuss how proteomics helps to accurately pinpoint candidate proteins to be functionally characterized, and how technological progresses create opportunities for studying low-abundance proteins, such as the ones related to the synthesis and transport of specialized metabolites. 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.

Comparative Proteomic Analysis of Cotton Fiber Development and Protein Extraction Method Comparison in Late Stage Fibers

The distinct stages of cotton fiber development and maturation serve as a single-celled model for studying the molecular mechanisms of plant cell elongation, cell wall development and cellulose biosynthesis. However, this model system of plant cell development is compromised for proteomic studies due to a lack of an efficient protein extraction method during the later stages of fiber development, because of a recalcitrant cell wall and the presence of abundant phenolic compounds. Here, we compared the quality and quantities of proteins extracted from 25 dpa (days post anthesis) fiber with multiple protein extraction methods and present a comprehensive quantitative proteomic study of fiber development from 10 dpa to 25 dpa. Comparative analysis using a label-free quantification method revealed 287 differentially-expressed proteins in the 10 dpa to 25 dpa fiber developmental period. Proteins involved in cell wall metabolism and regulation, cytoskeleton development and carbohydrate metabolism among other functional categories in four fiber developmental stages were identified. Our studies provide protocols for protein extraction from maturing fiber tissues for mass spectrometry analysis and expand knowledge of the proteomic profile of cotton fiber development.

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.

Comparative Transcriptomics Reveals Jasmonic Acid-Associated Metabolism Related to Cotton Fiber Initiation

Analysis of mutants and gene expression patterns provides a powerful approach for investigating genes involved in key stages of plant fiber development. In this study, lintless-fuzzless XinWX and linted-fuzzless XinFLM with a single genetic locus difference for lint were used to identify differentially expressed genes. Scanning electron microscopy showed fiber initiation in XinFLM at 0 days post anthesis (DPA). Fiber transcriptional profiling of the lines at three initiation developmental stages (-1, 0, 1 DPA) was performed using an oligonucleotide microarray. Loop comparisons of the differentially expressed genes within and between the lines was carried out, and functional classification and enrichment analysis showed that gene expression patterns during fiber initiation were heavily associated with hormone metabolism, transcription factor regulation, lipid transport, and asparagine biosynthetic processes, as previously reported. Further, four members of the allene-oxide cyclase (AOC) family that function in jasmonate biosynthesis were parallel up-regulation in fiber initiation, especially at -1 DPA, compared to other tissues and organs in linted-fuzzed TM-1. Real time-quantitative PCR (RT-qPCR) analysis in different fiber mutant lines revealed that AOCs were up-regulated higher at -1 DPA in lintless-fuzzless than that in linted-fuzzless and linted-fuzzed materials, and transcription of the AOCs was increased under jasmonic acid (JA) treatment. Expression analysis of JA biosynthesis-associated genes between XinWX and XinFLM showed that they were up-regulated during fiber initiation in the fuzzless-lintless mutant. Taken together, jasmonic acid-associated metabolism was related to cotton fiber initiation. Parallel up-regulation of AOCs expression may be important for normal fiber initiation development, while overproduction of AOCs might disrupt normal fiber development.

iTRAQ-facilitated proteomic profiling of anthers from a photosensitive male sterile mutant and wild-type cotton (Gossypium hirsutum L.)

Male sterility is a common phenomenon in flowering plants, and it has been successfully developed in several crops by taking advantage of heterosis. Cotton (Gossypium hirsutum L.) is an important economic crop, used mainly for the production of textile fiber. Using a space mutation breeding technique, a novel photosensitive genetic male sterile mutant CCRI9106 was isolated from the wild-type upland cotton cultivar CCRI040029. To use CCRI9106 in cotton hybrid breeding, it is of great importance to study the molecular mechanisms of its male sterility. Here, histological and iTRAQ-facilitated proteomic analyses of anthers were performed to explore male sterility mechanisms of the mutant. Scanning and transmission electron microscopy of the anthers showed that the development of pollen wall in CCRI9106 was severely defective with a lack of exine formation. At the protein level, 6121 high-confidence proteins were identified and 325 of them showed differential expression patterns between mutant and wild-type anthers. The proteins up- or down-regulated in MT anthers were mainly involved in exine formation, protein degradation, calcium ion binding,etc. These findings provide valuable information on the proteins involved in anther and pollen development, and contribute to elucidate the mechanism of male sterility in upland cotton.

Cottonseed protein, oil, and mineral status in near-isogenic Gossypium hirsutum cotton lines expressing fuzzy/linted and fuzzless/linted seed phenotypes under field conditions

Cotton is an important crop in the world and is a major source of oil for human consumption and cotton meal for livestock. Cottonseed nutrition (seed composition: protein, oil, and minerals) determines the quality of seeds. Therefore, maintaining optimum levels of cottonseed nutrition is critical. Physiological and genetic mechanisms controlling the levels of these constituents in cottonseed are still largely unknown. Our previous research conducted under greenhouse conditions showed that seed and leaf nutrition differed between fuzzless and fuzzy seed isolines. Therefore, the objective of this research was to investigate the seed fuzz phenotype (trait) effects on seed protein, oil, N, C, S, and minerals in five sets of near-isogenic mutant cotton lines for seed fuzz in a 2-year experiment under field condition to evaluate the stability of the effect of the trait on seed nutrition. The isolines (genotypes) in each set differ for the seed fuzz trait (fuzzless/linted seed line, N lines, and fuzzy/linted seed line, F lines). Results showed that seed protein was higher in the fuzzy genotype in all sets, but seed oil was higher in fuzzless genotype in all sets. The concentrations of seed Ca and C were higher in all fuzzless genotypes, but N, S, B, Fe, and Zn were higher in most of the fuzzy genotypes. Generally, minerals were higher in leaves of F lines, suggesting the translocation of minerals from leaves to seeds was limited. The research demonstrated that fiber development could be involved in cottonseed composition. This may be due to the involvement of fiber development in carbon and nitrogen metabolism, and the mobility of nutrients from leaves (source) to seed (sink). This information is beneficial to breeders to consider fuzzless cottonseed for potential protein and oil use and select for higher oil or higher protein content, and to physiologists to further understand the mobility of minerals to increase the quality of cottonseed nutrition for food and feed.

Protein expression changes during cotton fiber elongation in response to drought stress and recovery

An investigation to better understand the molecular mechanism of cotton (Gossypium hirsutum L.) fiber elongation in response to drought stress and recovery was conducted using a comparative proteomics analysis. Cotton plants (cv. NuCOTN 33B) were subjected to water deprivation for 10 days followed by a recovery period (with watering) of 5 days. The temporal changes in total proteins in cotton fibers were examined using 2DE. The results revealed that 163 proteins are significantly drought responsive. MS analysis led to the identification of 132 differentially expressed proteins that include some known as well as some novel drought-responsive proteins. These drought responsive fiber proteins in NuCOTN 33B are associated with a variety of cellular functions, i.e. signal transduction, protein processing, redox homeostasis, cell wall modification, metabolisms of carbon, energy, lipid, lignin, and flavonoid. The results suggest that the enhancement of the perception of drought stress, a new balance of the metabolism of the biosynthesis of cell wall components and cytoskeleton homeostasis plays an important role in the response of cotton fibers to drought stress. Overall, the current study provides an overview of the molecular mechanism of drought response in cotton fiber cells.

Comparative proteomic and biochemical analyses reveal different molecular events occurring in the process of fiber initiation between wild-type allotetraploid cotton and its fuzzless-lintless mutant

To explore lint fiber initiation-related proteins in allotetraploid cotton (Gossypium hirsutum L.), a comparative proteomic analysis was performed between wild-type cotton (Xu-142) and its fuzzless-lintless mutant (Xu-142-fl) at five developmental time points for lint fiber initiation from -3 to +3 days post-anthesis (dpa). Using two-dimensional gel electrophoresis (2-DE) combined with mass spectrometry (MS) analyses, 91 differentially accumulated protein (DAP) species that are related to fiber initiation were successfully identified, of which 58 preferentially accumulated in the wild-type and 33 species in the fl mutant. These DAPs are involved in various cellular and metabolic processes, mainly including important energy/carbohydrate metabolism, redox homeostasis, amino acid and fatty acid biosynthesis, protein quality control, cytoskeleton dynamics, and anthocyanidin metabolism. Further physiological and biochemical experiments revealed dynamic changes in the carbohydrate flux and H₂O₂ levels in the cotton fiber initiation process. Compared with those in the fl mutant, the contents of glucose and fructose in wild-type ovules sharply increased after anthesis with a relatively higher rate of amino acid biosynthesis. The relative sugar starvation and lower rate of amino acid biosynthesis in the fl mutant ovules may impede the carbohydrate/energy supply and cell wall synthesis, which is consistent with the proteomic results. However, the H₂O₂ burst was only observed in the wild-type ovules on the day of anthesis. Cotton boll injection experiments in combination with electron microscope observation collectively indicated that H₂O₂ burst, which is negatively regulated by ascorbate peroxidases (APx), plays an important role in the fiber initiation process. Taken together, our study demonstrates a putative network of DAP species related to fiber initiation in cotton ovules and provides a foundation for future studies on the specific functions of these proteins in fiber development.

Proteomics profiling of fiber development and domestication in upland cotton (Gossypium hirsutum L.)

Comparative proteomic analyses were performed to detail the evolutionary consequences of strong directional selection for enhanced fiber traits in modern upland cotton (Gossypium hirsutum L.). Using two complementary proteomic approaches, 2-DE and iTRAQ LC-MS/MS, fiber proteomes were examined for four representative stages of fiber development. Approximately 1,000 protein features were characterized using each strategy, collectively resulting in the identification and functional categorization of 1,223 proteins. Unequal contributions of homoeologous proteins were detected for over a third of the fiber proteome, but overall expression was balanced with respect to the genome-of-origin in the allopolyploid G. hirsutum. About 30% of the proteins were differentially expressed during fiber development within wild and domesticated cotton. Notably, domestication was accompanied by a doubling of protein developmental dynamics for the period between 10 and 20 days following pollination. Expression levels of 240 iTRAQ proteins and 293 2-DE spots were altered by domestication, collectively representing multiple cellular and metabolic processes, including metabolism, energy, protein synthesis and destination, defense and stress response. Analyses of homoeolog-specific expression indicate that duplicated gene products in cotton fibers can be differently regulated in response to selection. These results demonstrate the power of proteomics for the analysis of crop domestication and phenotypic evolution.

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

An iTRAQ-based proteomics of ovules from the upland cotton species Gossypium hirsutum and its fuzzless-lintless mutant was performed, and finally 2729 proteins that preferentially accumulated at anthesis in wild-type ovules were identified. We confirmed that the gene expression levels of 2005 among these proteins also increased by performing an RNA sequencing transcriptomics. Expression of proteins involved in carboxylic acid metabolism, small-molecule metabolic processes, hormone regulation, and lipid metabolism was significantly enhanced in wild-type ovules. Quantitative real-time PCR verified the increased expression of 26 genes involved in these processes. Cotton 3-hydroxyacyl-CoA dehydratase (GhPAS2) catalyzing the third reaction of very long-chain fatty acid (VLCFA) biosynthesis, accumulated at anthesis in wild-type ovules. Heterogeneous expression of GhPAS2 restored viability to the Saccharomyces cerevisiae haploid psh1-deletion strain deficient in PAS2 activity. Application of VLCFA biosynthesis inhibitor acetochlor (2-chloro-N-[ethoxymethyl]-N-[2-ethyl-6-methyl-phenyl]-acetamide; ACE) and gibberellic acid to the unfertilized cotton ovules significantly suppressed fiber cell protrusion. In this study, the profiling of gene expression at both transcriptome and proteome levels provides new insights into cotton fiber cell initiation.

BIOLOGICAL SIGNIFICANCE

Cotton fiber initiation determines the ultimate number of fibers per ovule, thereby determining fiber yield. In total, 2729 proteins were preferentially accumulated in wild-type ovules at anthesis. The most up-regulated proteins were assigned to carboxylic acid metabolism, small-molecule metabolic processes, hormone regulation, and lipid metabolism. In consistence with these findings, we characterized GhPAS2 gene coding for the enzyme that catalyzes VLCFA production. VLCFA biosynthesis inhibitor, acetochlor, was shown to significantly suppress fiber initiation. This study provides a genome-scale transcriptomic and proteomic characterization of fiber initial cells, laying a solid basis for further investigation of the molecular processes governing fiber cell development.

Quantitative phosphoproteomic profiling of fiber differentiation and initiation in a fiberless mutant of cotton

Background

The cotton (Gossypium spp.) fiber cell is an important unicellular model for studying cell differentiation. There is evidence suggesting that phosphorylation is a critical post-translational modification involved in regulation of a wide range of cell activities. Nevertheless, the sites of phosphorylation in G. hirsutum and their regulatory roles in fiber cell initiation are largely unknown. In this study, we employed a mass spectrometry-based phosphoproteomics to conduct a global and site-specific phosphoproteome profiling between ovules of a fuzzless-lintless (fl) Upland cotton (G. hirsutum) mutant and its isogenic parental wild type (WT) at -3 and 0 days post-anthesis (DPA).

Results

A total of 830 phosphopeptides and 1,592 phosphorylation sites from 619 phosphoproteins were identified by iTRAQ (isobaric tags for relative and absolute quantitation). Of these, 76 phosphoproteins and 1,100 phosphorylation sites were identified for the first time after searching the P³DB public database using the BLAST program. Among the detected phosphopeptides, 69 were differentially expressed between the fl mutant and its WT in ovules at -3 and 0 DPA. An analysis using the Motif-X program uncovered 19 phosphorylation motifs, 8 of which were unique to cotton. A further metabolic pathway analysis revealed that the differentially phosphorylated proteins were involved in signal transduction, protein modification, carbohydrate metabolic processes, and cell cycle and cell proliferation.

Conclusions

Our phosphoproteomics-based research provides the first global overview of phosphorylation during cotton fiber initiation, and also offers a helpful dataset for elucidation of signaling networks in fiber development of G. hirsutum.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-466) contains supplementary material, which is available to authorized users.

Cotton proteomics for deciphering the mechanism of environment stress response and fiber development

Cotton fiber is considered as the backbone of the textile industry. The productivity of cotton crop is severely hampered by the occurrence of pathogens, pests, and various environmental factors. Nevertheless, cotton plant has developed sophisticated mechanisms to respond to environment stresses to avoid detrimental effects on its growth and development. Therefore, understanding the mechanisms of cotton fiber development and environment stress response is of considerable interest for designing agriculture breeding strategies to ensure sustainable productivity. The application of proteomics technologies to advance our knowledge in cotton fiber development and abiotic/biotic stress tolerance has increased dramatically in the last 5years as evidenced by the large amount of publications in this area. This review summarizes the work which has been reported for cotton proteomics and evaluates the findings in context of the approaches that are widely employed with the aim to generate novel insight useful for cotton improvement.

BIOLOGICAL SIGNIFICANCE

Cotton (Gossypium spp.) is considered as the foremost commercially important fiber crop grown all over the world and is deemed as the backbone of the textile industry. Cotton is also an important source of edible oil seed and a nutrient-rich food crop as cottonseed contains high-quality protein and oil. The growth and productivity of cotton crop are often hampered by various biotic stress factors, such as insect pests and pathogens. In addition, cotton plants are frequently subjected to unavoidable environmental factors that cause abiotic stress, such as salt, heat and drought. Proteomic techniques provide one of the best options for understanding the gene function and phenotypic changes during cotton fiber development and stress response. This review first summarizes the work which has been reported for cotton proteomics about cotton fiber development and abiotic/biotic stress tolerance, and also evaluates the findings in context of the approaches that are widely employed during last 5years, with the aim to generate novel insight useful for cotton improvement. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Identification of differentially accumulating pistil proteins associated with self-incompatibility of non-heading Chinese cabbage

Non-heading Chinese cabbage (Brassica campestris L. ssp. chinensis Makino), an important vegetable crop in China, exhibits a typical sporophytic self-incompatibility (SI) system. To better understand the mechanism of SI response and identify potential candidate proteins involved in the SI system of this vegetable crop, the proteomic approach was taken to identify differential accumulating pistil proteins. Pistils were collected at 0 h and 2 h after self-pollination at anthesis in self-incompatible and compatible lines of non-heading Chinese cabbage, and total proteins were extracted and separated by two-dimensional gel electrophoresis (2-DE). A total of 25 protein spots that displayed differential abundance were identified by matrix-assisted laser desorption/ionisation-time of flight mass spectrometry (MALDI-TOF/TOF MS) and peptide mass fingerprinting (PMF). Among them, 22 protein spots were confidently established. The mRNA levels of the corresponding genes were detected by quantitative RT-PCR. The 22 identified protein spots are involved in energy metabolism (four), protein biosynthesis (three), photosynthesis (six), stress response and defence (five), and protein degradation (four). Among these potential candidate proteins, UDP-sugar pyrophosphorylase could be involved in sucrose degradation to influence pollen germination and growth. Glutathione S-transferases could be involved in pollen maturation, and affect pollen fertility. Senescence-associated cysteine protease, which is related to programmed cell death, could be mainly related to self pollen recognition of non-heading Chinese cabbage. The study will contribute to further investigations of molecular mechanism of sporophytic SI in Brassicaceae.

Effects of fuzzless cottonseed phenotype on cottonseed nutrient composition in near isogenic cotton (Gossypium hirsutum L.) mutant lines under well-watered and water stress conditions

There is no information available on the effect of fuzzless seed trait on cottonseed nutrient composition (minerals, N, S, protein, and oil) under drought stress. The objective of this research was to investigate the effect of the fuzzless seed trait on cottonseed nutrients using five sets of near-isogenic lines (NILs). Each set consists of two lines that share the same genetic background, but differ in seed fuzziness (fuzzy, F; fuzzless, N). The near isogenic lines will enable us to compare the effect of the trait without confounding the genotypic background effects. We hypothesized that since the fuzzless trait involved in fiber initiation development, and was reported to be involved in biochemical, molecular, and genetic processes, this trait may also alter cottonseed nutrient composition. Results showed that NIL sets accumulated different levels of minerals in seeds and leaves, and the fuzzless trait (N) in most of the lines altered seed and leaf mineral accumulations when compared with fuzzy lines (F) or the control line. For example, K, P, Mg, Cu, and Na concentrations in seeds were higher in MD N and STV N than in their equivalent MD F and STV F lines. Leaf concentrations of Ca, K, Mg, S, B, Cu, and Fe in MD N lines were higher than MD F line. Lower levels of nutrients in seeds and leaves were observed under water stress conditions, especially Ca, Mg, N, and B in seeds.Generally and with few exceptions, seed protein was higher in fuzzy lines than in fuzzless lines; however, seed oil was higher in fuzzless lines than in fuzzy lines. Our research demonstrated that fuzzless trait altered the composition and level of nutrients in seed and leaves in well watered and water stressed plants. Differences in protein and oil between fuzzy and fuzzless seeds may indicate alteration in nitrogen and carbon fixation and metabolism. The differential accumulation of seed nutrients in this germplasm could be used by cotton breeders to select for higher cottonseed quality.

Glycoproteome of elongating cotton fiber cells

Cotton ovule epidermal cell differentiation into long fibers primarily depends on wall-oriented processes such as loosening, elongation, remodeling, and maturation. Such processes are governed by cell wall bound structural proteins and interacting carbohydrate active enzymes. Glycosylation plays a major role in the structural, functional, and localization aspects of the cell wall and extracellular destined proteins. Elucidating the glycoproteome of fiber cells would reflect its wall composition as well as compartmental requirement, which must be system specific. Following complementary proteomic approaches, we have identified 334 unique proteins comprising structural and regulatory families. Glycopeptide-based enrichment followed by deglycosylation with PNGase F and A revealed 92 unique peptides containing 106 formerly N-linked glycosylated sites from 67 unique proteins. Our results showed that structural proteins like arabinogalactans and carbohydrate active enzymes were relatively more abundant and showed stage- and isoform-specific expression patterns in the differentiating fiber cell. Furthermore, our data also revealed the presence of heterogeneous and novel forms of structural and regulatory glycoproteins. Comparative analysis with other plant glycoproteomes highlighted the unique composition of the fiber glycoproteome. The present study provides the first insight into the identity, abundance, diversity, and composition of the glycoproteome within single celled cotton fibers. The elucidated composition also indirectly provides clues about unicellular compartmental requirements underlying single cell differentiation.

Separation of integument and nucellar tissues from cotton ovules (Gossypium hirsutum L.) for both high- and low-throughput molecular applications

Here we present a quick and low-cost method to separate the different layers of tissue from the ovules and young seeds of cotton (Gossypium hirsutum L.) for use in high- and low-throughput molecular applications. This method is performed at room temperature using standard laboratory equipment and does not require embedding of the samples, time-consuming fixation, or micro-sectioning procedures. We show that the three main tissues can be efficiently separated from isolated ovules collected on the day of anthesis. RNA and genomic DNA extracted from tissues separated by this method are of good quality and suitable for a variety of molecular applications to study the early stages of cotton seed and fiber development.

Functional genomics of fuzzless-lintless mutant of Gossypium hirsutum L. cv. MCU5 reveal key genes and pathways involved in cotton fibre initiation and elongation

Background

Fuzzless-lintless cotton mutants are considered to be the ideal material to understand the molecular mechanisms involved in fibre cell development. Although there are few reports on transcriptome and proteome analyses in cotton at fibre initiation and elongation stages, there is no comprehensive comparative transcriptome analysis of fibre-bearing and fuzzless-lintless cotton ovules covering fibre initiation to secondary cell wall (SCW) synthesis stages. In the present study, a comparative transcriptome analysis was carried out using G. hirsutum L. cv. MCU5 wild-type (WT) and it’s near isogenic fuzzless-lintless (fl) mutant at fibre initiation (0 dpa/days post anthesis), elongation (5, 10 and 15 dpa) and SCW synthesis (20 dpa) stages.

Results

Scanning electron microscopy study revealed the delay in the initiation of fibre cells and lack of any further development after 2 dpa in the fl mutant. Transcriptome analysis showed major down regulation of transcripts (90%) at fibre initiation and early elongation (5 dpa) stages in the fl mutant. Majority of the down regulated transcripts at fibre initiation stage in the fl mutant represent calcium and phytohormone mediated signal transduction pathways, biosynthesis of auxin and ethylene and stress responsive transcription factors (TFs). Further, transcripts involved in carbohydrate and lipid metabolisms, mitochondrial electron transport system (mETS) and cell wall loosening and elongation were highly down-regulated at fibre elongation stage (5–15 dpa) in the fl mutant. In addition, cellulose synthases and sucrose synthase C were down-regulated at SCW biosynthesis stage (15–20 dpa). Interestingly, some of the transcripts (~50%) involved in phytohormone signalling and stress responsive transcription factors that were up-regulated at fibre initiation stage in the WT were found to be up-regulated at much later stage (15 dpa) in fl mutant.

Conclusions

Comparative transcriptome analysis of WT and its near isogenic fl mutant revealed key genes and pathways involved at various stages of fibre development. Our data implicated the significant role of mitochondria mediated energy metabolism during fibre elongation process. The delayed expression of genes involved in phytohormone signalling and stress responsive TFs in the fl mutant suggests the need for a coordinated expression of regulatory mechanisms in fibre cell initiation and differentiation.

Genomes for jeans: cotton genomics for engineering superior fiber

Twenty years ago, scientists predicted that better understanding of fiber development would lead to novel ways to engineer superior cotton fiber. Advances in genetic resources, DNA markers, DNA sequence information, and gene expression data have indeed provided new insights into fiber initiation, elongation and maturation. Many exciting applications of this knowledge offer the potential to select better cotton genotypes more effectively in mainstream breeding programs or engineer genotypes with improved agronomic and/or quality traits. Here, we discuss recent progress in understanding genes involved in fiber development, and their regulation and manipulation to engineer improved fibers. Better understanding of quantitative trait loci/gene interactions that influence fiber quality and yield may help to tailor superior cotton genotypes to diverse environments.