Effect of H2O2 on fiber initiation using fiber retardation initiation mutants in cotton (Gossypium hirsutum)

Comparative Transcriptomic Analysis of Gossypium hirsutum Fiber Development in Mutant Materials (xin w 139) Provides New Insights into Cotton Fiber Development

Cotton is the most widely planted fiber crop in the world, and improving cotton fiber quality has long been a research hotspot. The development of cotton fibers is a complex process that includes four consecutive and overlapping stages, and although many studies on cotton fiber development have been reported, most of the studies have been based on cultivars that are promoted in production or based on lines that are used in breeding. Here, we report a phenotypic evaluation of Gossypium hirsutum based on immature fiber mutant (xin w 139) and wild-type (Xin W 139) lines and a comparative transcriptomic study at seven time points during fiber development. The results of the two-year study showed that the fiber length, fiber strength, single-boll weight and lint percentage of xin w 139 were significantly lower than those of Xin W 139, and there were no significant differences in the other traits. Principal component analysis (PCA) and cluster analysis of the RNA-sequencing (RNA-seq) data revealed that these seven time points could be clearly divided into three different groups corresponding to the initiation, elongation and secondary cell wall (SCW) synthesis stages of fiber development, and the differences in fiber development between the two lines were mainly due to developmental differences after twenty days post anthesis (DPA). Differential expression analysis revealed a total of 5131 unique differentially expressed genes (DEGs), including 290 transcription factors (TFs), between the 2 lines. These DEGs were divided into five clusters. Each cluster functional category was annotated based on the KEGG database, and different clusters could describe different stages of fiber development. In addition, we constructed a gene regulatory network by weighted correlation network analysis (WGCNA) and identified 15 key genes that determined the differences in fiber development between the 2 lines. We also screened seven candidate genes related to cotton fiber development through comparative sequence analysis and qRT-PCR; these genes included three TFs (GH_A08G1821 (bHLH), GH_D05G3074 (Dof), and GH_D13G0161 (C3H)). These results provide a theoretical basis for obtaining an in-depth understanding of the molecular mechanism of cotton fiber development and provide new genetic resources for cotton fiber research.

Systematically and Comprehensively Understanding the Regulation of Cotton Fiber Initiation: A Review

Cotton fibers provide an important source of raw materials for the textile industry worldwide. Cotton fiber is a kind of single cell that differentiates from the epidermis of the ovule and provides a perfect research model for the differentiation and elongation of plant cells. Cotton fiber initiation is the first stage throughout the entire developmental process. The number of fiber cell initials on the seed ovule epidermis decides the final fiber yield. Thus, it is of great significance to clarify the mechanism underlying cotton fiber initiation. Fiber cell initiation is controlled by complex and interrelated regulatory networks. Plant phytohormones, transcription factors, sugar signals, small signal molecules, functional genes, non-coding RNAs, and histone modification play important roles during this process. Here, we not only summarize the different kinds of factors involved in fiber cell initiation but also discuss the mechanisms of these factors that act together to regulate cotton fiber initiation. Our aim is to synthesize a systematic and comprehensive review of different factors during fiber initiation that will provide the basics for further illustrating these mechanisms and offer theoretical guidance for improving fiber yield in future molecular breeding work.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

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.

Dynamic roles and intricate mechanisms of ethylene in epidermal hair development in Arabidopsis and cotton

Ethylene affects many aspects of plant growth and development, including root hairs and trichomes growth in Arabidopsis, as well as fiber development in cotton, though the underlying mechanism is unclear. In this article, we update the research progress associated with the main genes in ethylene biosynthesis and signaling pathway, and we propose a clear ethylene pathway based on genome-wide identification of homologues in cotton. Expression pattern analysis using transcriptome data revealed that some candidate genes may contribute to cotton fiber development through the ethylene pathway. Moreover, we systematically summarized the effects of ethylene on the development of epidermal hair and the underlying regulatory mechanisms in Arabidopsis. Based on the knowledge of ethylene-promoted cell differentiation, elongation, and development in different tissues or plants, we advised a possible regulatory network for cotton fiber development with ethylene as the hub. Importantly, we emphasized the roles of ethylene as an important node in regulating cotton vegetative growth, and stress resistance, and suggested utilizing multiple methods to subtly modify ethylene synthesis or signaling in a tissue or spatiotemporal-specific manner to clarify its exact effect on architecture, adaptability of the plant, and fiber development, paving the way for basic research and genetic improvement of the cotton crop.

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.

Looking into 'hair tonics' for cotton fiber initiation

Cotton fiber is the most important source of cellulose for the global textile industry. These hair-like single-celled trichomes develop from ovule epidermis. They are classified into long spinnable lint and short fuzz. A key objective in the cotton industry is to breed elite cultivars with fuzzless seeds carrying high lint yield. Molecular basis underlying lint and fuzz initiation remains obscure. Recent studies indicate fiber initiation is under the control of MYB-bHLH-WDR (MBW) transcription factor complex. Based on molecular genetic studies and gene expression patterns linking fiber phenotypes, we propose that specific but different sets of MBW genes are required to precisely regulate the initiation of the lint and fuzz fibers. Emerging evidence further points to sugar signaling as a 'hair-tonic' to boost fiber initiation through interaction with MBW complex and auxin signaling. An integrative model is provided as a conceptual framework for future studies to dissect the molecular network responsible for cotton fiber 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.

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.

Genome-wide identification of the GhARF gene family reveals that GhARF2 and GhARF18 are involved in cotton fibre cell initiation

Auxin signalling plays an essential role in regulating plant development. Auxin response factors (ARFs), which are critical components of auxin signalling, modulate the expression of early auxin-responsive genes by binding to auxin response factor elements (AuxREs). However, there has been no comprehensive characterization of this gene family in cotton. Here, we identified 56 GhARF genes in the assembled Gossypium hirsutum genome. This gene family was divided into 17 subfamilies, and 44 members of them were distributed across 21 chromosomes. GhARF6 and GhARF11 subfamily genes were predominantly expressed in vegetative tissues, whereas GhARF2 and GhARF18 subfamily genes were highly expressed during seed fibre cell initiation. GhARF2-1 and GhARF18-1 were exclusively expressed in trichomes, organs similar to cotton seed fibre cells, and overexpression of these genes in Arabidopsis enhances trichome initiation. Comparative transcriptome analysis combined with AuxRE prediction revealed 11 transcription factors as potential target genes of GhARF2 and GhARF18. Six of these genes were significantly expressed during seed fibre cell initiation and were bound by GhARF2-1 and GhARF18-1 in yeast one-hybrid assays. Our results suggest that GhARF2 and GhARF18 genes may be key regulators of cotton seed fibre initiation by regulating the expression of several transcription factor genes. This study deepens our understanding of auxin-mediated initiation of cotton seed fibre cells and helps us in breeding better cotton varieties in the future.

Inhibition of Heat Shock proteins HSP90 and HSP70 induce oxidative stress, suppressing cotton fiber development

Cotton fiber is a specialized unicellular structure useful for the study of cellular differentiation and development. Heat shock proteins (HSPs) have been shown to be involved in various developmental processes. Microarray data analysis of five Gossypium hirsutum genotypes revealed high transcript levels of GhHSP90 and GhHSP70 genes at different stages of fiber development, indicating their importance in the process. Further, we identified 26 and 55 members of HSP90 and HSP70 gene families in G. hirsutum. The treatment of specific inhibitors novobiocin (Nov; HSP90) and pifithrin/2-phenylethynesulfonamide (Pif; HSP70) in in-vitro cultured ovules resulted in a fewer number of fiber initials and retardation in fiber elongation. The molecular chaperone assay using bacterially expressed recombinant GhHSP90-7 and GhHSP70-8 proteins further confirmed the specificity of inhibitors. HSP inhibition disturbs the H₂O₂ balance that leads to the generation of oxidative stress, which consequently results in autophagy in the epidermal layer of the cotton ovule. Transmission electron microscopy (TEM) of inhibitor-treated ovule also corroborates autophagosome formation along with disrupted mitochondrial cristae. The perturbations in transcript profile of HSP inhibited ovules show differential regulation of different stress and fiber development-related genes and pathways. Altogether, our results indicate that HSP90 and HSP70 families play a crucial role in cotton fiber differentiation and development by maintaining cellular homeostasis.

Genetics and evolution of MIXTA genes regulating cotton lint fiber development

Cotton, with cellulose-enriched mature fibers, is the largest source of natural textiles. Through a map-based cloning strategy, we isolated an industrially important lint fiber development gene (Li₃ ) that encodes an MYB-MIXTA-like transcription factor (MML) on chromosome D12 (GhMML4_D12). Virus-induced gene silencing or decreasing the expression of the GhMML4_D12 gene in n₂ NSM plants resulted in a significant reduction in epidermal cell prominence and lint fiber production. GhMML4_D12 is arranged in tandem with GhMML3, another MIXTA gene responsible for fuzz fiber development. These two very closely related MIXTA genes direct fiber initiation production in two specialized cell forms: lint and fuzz fibers. They may control the same metabolic pathways in different cell types. The MIXTAs expanded in Malvaceae during their evolution and produced a Malvaceae-specific family that regulates epidermal cell differentiation, different from the gene family that regulates leaf hair trichome development. Cotton has developed a unique transcriptional regulatory network for fiber development. Characterization of target genes regulating fiber production has provided insights into the molecular mechanisms underlying cotton fiber development and has allowed the use of genetic engineering to increase lint yield by inducing more epidermal cells to develop into lint rather than fuzz fibers.

Transcript profiling of genes expressed during fibre development in diploid cotton (Gossypium arboreum L.)

Background

Cotton fibre is a single cell and it is one of the best platforms for unraveling the genes express during various stages of fibre development. There are reports devoted to comparative transcriptome study on fiber cell initiation and elongation in tetraploid cultivated cotton. However, in the present investigation, comparative transcriptome study was made in diploid cultivated cotton using isogenic fuzzy-lintless (Fl) and normal fuzzy linted (FL) lines belong to Gossypium arboreum, diploid species at two stages, 0 and 10 dpa (days post anthesis), using Affymetrix cotton GeneChip genome array.

Result

Scanning electron microscopy (SEM) analysis uncovered the occurrence of few fibre cell initials in the Fl line as compared to many in Normal FL at −2 and 0 dpa. However, at 10 dpa there were no fibre cells found elongated in Fl but many elongated cells were found in FL line. Up-regulation of transcription factors, AP2-EREBP, C2H2, C3H, HB and WRKY was observed at 0 dpa whereas in 10 dpa transcription factors, AP2-EREBP, AUX/IAA, bHLH, C2H2, C3H, HB, MYB, NAC, Orphans, PLATZ and WRKY were found down regulated in Fl line. These transcription factors were mainly involved in metabolic pathways such as phytohormone signaling, energy metabolism of cell, fatty acid metabolism, secondary metabolism and other signaling pathways and are related directly or indirectly in fiber development. Quantitative real-time PCR was performed to check fold up or down-regulation of these genes and transcription factors (TFs) down regulated in mutants as compared to normal at 0 and 10 dpa.

Conclusion

This study elucidates that the up-regulation of transcription factors like AP2-EREBP, C2H2, C3H, HB, WRKY and phytohormone signaling genes at 0 dpa and their down-regulation at the 10 dpa might have constrain the fibre elongation in fuzzy-lintless line. Along with this the down-regulation of genes involved in synthesis of VLCFA chain, transcripts necessary for energy and cell wall metabolism, EXPANSINs, arabinogalactan proteins (AGPs), tubulin might also be the probable reason for reduced growth of fibres in the Fl. Plant receptor-like kinases (RLKs), Leucine Rich Repeats) LRR- family protein and signal transduction coding for mitogen-activated protein kinase (MAPK) cascade, have been engaged in coordination of cell elongation and SCW biosynthesis, down-regulation of these might loss the function leads to reduced fibre growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4066-y) contains supplementary material, which is available to authorized users.

Cotton cytosolic pyruvate kinase GhPK6 participates in fast fiber elongation regulation in a ROS-mediated manner

Cotton cytosolic pyruvate kinase GhPK6 is preferentially expressed in the late stage of fiber elongation process, transgenic experiments indicated that its expression level was negatively correlated to cell expansion rate. Pyruvate kinase (PK) plays vital regulatory roles in rapid cell growth in mammals. However, the function of PK in plant cell growth remains unclear. In allotetraploid upland cotton (Gossypium hirsutum L.), a total of 33 PK genes are encoded by the genome. Analysis of the transcriptome data indicated that only two cytosolic PK genes, GhPK6 and its duplicated gene GhPK26, are preferentially expressed in elongating cotton fiber cells. RT-qPCR and western blot analyses revealed that the expression of GhPK6 was negatively correlated with fiber elongation rate, which well explains the observed sharp increase of cytosolic PK activity at the end of fast fiber elongation process. Furthermore, virus-induced gene silencing of GhPK6 in cotton plants resulted in increased fiber cell elongation and reduced reactive oxygen species (ROS) accumulation. On the contrary, Arabidopsis plants ectopically expressing GhPK6 exhibited ROS-mediated growth inhibition, whereas the addition of ROS scavenging reagents could partly rescue this inhibition. These data collectively suggested that GhPK6 might play an important role in regulating cotton fiber elongation in a ROS-dependent inhibition manner.

Characterization of in vivo phosphorylation modification of differentially accumulated proteins in cotton fiber-initiation process

Initiation of cotton fiber from ovule epidermal cells determines the ultimate number of fibers per cotton ovule, making it one of the restriction factors of cotton fiber yield. Previous comparative proteomics studies have collectively revealed 162 important differentially accumulated proteins (DAPs) in cotton fiber-initiation process, however, whether and how post-translational modifications, especially phosphorylation modification, regulate the expression and function of the DAPs are still unclear. Here we reported the successful identification of 17 phosphopeptides from 16 phosphoproteins out of the 162 DAPs using the integrated bioinformatics analyses of peptide mass fingerprinting data and targeted MS/MS identification method. In-depth analyses indicated that 15 of the 17 phosphorylation sites were novel phosphorylation sites first identified in plants, whereas 6 of the 16 phosphoproteins were found to be the phosphorylated isoforms of 6 proteins. The phosphorylation-regulated dynamic protein network derived from this study not only expanded our understanding of the cotton fiber-initiation process, but also provided a valuable resource for future functional studies of the phosphoproteins.

AtRAV1 and AtRAV2 overexpression in cotton increases fiber length differentially under drought stress and delays flowering

There is a longstanding problem of an inverse relationship between cotton fiber qualities versus high yields. To better understand drought stress signaling and adaptation in cotton (Gossypium hirsutum) fiber development, we expressed the Arabidopsis transcription factors RELATED_TO_ABA-INSENSITIVE3/VIVIPAROUS1/(RAV1) and AtRAV2, which encode APETALA2-Basic3 domain proteins shown to repress transcription of FLOWERING_LOCUS_T (FT) and to promote stomatal opening cell-autonomously. In three years of field trials, we show that AtRAV1 and AtRAV2-overexpressing cotton had ∼5% significantly longer fibers with only marginal decreases in yields under well-watered or drought stress conditions that resulted in 40-60% yield penalties and 3-7% fiber length penalties in control plants. The longer transgenic fibers from drought-stressed transgenics could be spun into yarn which was measurably stronger and more uniform than that from well-watered control fibers. The transgenic AtRAV1 and AtRAV2 lines flowered later and retained bolls at higher nodes, which correlated with repression of endogenous GhFT-Like (FTL) transcript accumulation. Elevated expression early in development of ovules was observed for GhRAV2L, GhMYB25-Like (MYB25L) involved in fiber initiation, and GhMYB2 and GhMYB25 involved in fiber elongation. Altered expression of RAVs controlling critical nodes in developmental and environmental signaling hierarchies has the potential for phenotypic modification of crops.

Metabolomic and transcriptomic insights into how cotton fiber transitions to secondary wall synthesis, represses lignification, and prolongs elongation

Background

The morphogenesis of single-celled cotton fiber includes extreme elongation and staged cell wall differentiation. Designing strategies for improving cotton fiber for textiles and other uses relies on uncovering the related regulatory mechanisms. In this research we compared the transcriptomes and metabolomes of two Gossypium genotypes, Gossypium barbadense cv Phytogen 800 and G. hirsutum cv Deltapine 90. When grown in parallel, the two types of fiber developed similarly except for prolonged fiber elongation in the G. barbadense cultivar. The data were collected from isolated fibers between 10 to 28 days post anthesis (DPA) representing: primary wall synthesis to support elongation; transitional cell wall remodeling; and secondary wall cellulose synthesis, which was accompanied by continuing elongation only in G. barbadense fiber.

Results

Of 206 identified fiber metabolites, 205 were held in common between the two genotypes. Approximately 38,000 transcripts were expressed in the fiber of each genotype, and these were mapped to the reference set and interpreted by homology to known genes. The developmental changes in the transcriptomes and the metabolomes were compared within and across genotypes with several novel implications. Transitional cell wall remodeling is a distinct stable developmental stage lasting at least four days (18 to 21 DPA). Expression of selected cell wall related transcripts was similar between genotypes, but cellulose synthase gene expression patterns were more complex than expected. Lignification was transcriptionally repressed in both genotypes. Oxidative stress was lower in the fiber of G. barbadense cv Phytogen 800 as compared to G. hirsutum cv Deltapine 90. Correspondingly, the G. barbadense cultivar had enhanced capacity for management of reactive oxygen species during its prolonged elongation period, as indicated by a 138-fold increase in ascorbate concentration at 28 DPA.

Conclusions

The parallel data on deep-sequencing transcriptomics and non-targeted metabolomics for two genotypes of single-celled cotton fiber showed that a discrete developmental stage of transitional cell wall remodeling occurs before secondary wall cellulose synthesis begins. The data showed how lignification can be transcriptionally repressed during secondary cell wall synthesis, and they implicated enhanced capacity to manage reactive oxygen species through the ascorbate-glutathione cycle as a positive contributor to fiber length.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1708-9) contains supplementary material, which is available to authorized users.

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.

Phytohormonal networks promote differentiation of fiber initials on pre-anthesis cotton ovules grown in vitro and in planta

The number of cotton (Gossypium sp.) ovule epidermal cells differentiating into fiber initials is an important factor affecting cotton yield and fiber quality. Despite extensive efforts in determining the molecular mechanisms regulating fiber initial differentiation, only a few genes responsible for fiber initial differentiation have been discovered. To identify putative genes directly involved in the fiber initiation process, we used a cotton ovule culture technique that controls the timing of fiber initial differentiation by exogenous phytohormone application in combination with comparative expression analyses between wild type and three fiberless mutants. The addition of exogenous auxin and gibberellins to pre-anthesis wild type ovules that did not have visible fiber initials increased the expression of genes affecting auxin, ethylene, ABA and jasmonic acid signaling pathways within 1 h after treatment. Most transcripts expressed differentially by the phytohormone treatment in vitro were also differentially expressed in the ovules of wild type and fiberless mutants that were grown in planta. In addition to MYB25-like, a gene that was previously shown to be associated with the differentiation of fiber initials, several other differentially expressed genes, including auxin/indole-3-acetic acid (AUX/IAA) involved in auxin signaling, ACC oxidase involved in ethylene biosynthesis, and abscisic acid (ABA) 8'-hydroxylase an enzyme that controls the rate of ABA catabolism, were co-regulated in the pre-anthesis ovules of both wild type and fiberless mutants. These results support the hypothesis that phytohormonal signaling networks regulate the temporal expression of genes responsible for differentiation of cotton fiber initials in vitro and in planta.

The calcium sensor GhCaM7 promotes cotton fiber elongation by modulating reactive oxygen species (ROS) production

Fiber elongation is the key determinant of fiber quality and output in cotton (Gossypium hirsutum). Although expression profiling and functional genomics provide some data, the mechanism of fiber development is still not well understood. Here, a gene encoding a calcium sensor, GhCaM7, was isolated based on its high expression level relative to other GhCaMs in fiber cells at the fast elongation stage. The level of expression of GhCaM7 in the wild-type and the fuzzless/lintless mutant correspond to the presence and absence, respectively, of fiber initials. Overexpressing GhCaM7 promotes early fiber elongation, whereas GhCaM7 suppression by RNAi delays fiber initiation and inhibits fiber elongation. Reactive oxygen species (ROS) play important roles in early fiber development. ROS induced by exogenous hydrogen peroxide (H2 O2 ) and Ca(2+) starvation promotes early fiber elongation. GhCaM7 overexpression fiber cells show increased ROS concentrations compared with the wild-type, while GhCaM7 RNAi fiber cells have reduced concentrations. Furthermore, we show that H2 O2 enhances Ca(2+) influx into the fiber and feedback-regulates the expression of GhCaM7. We conclude that GhCaM7, Ca(2+) and ROS are three important regulators involved in early fiber elongation. GhCaM7 might modulate ROS production and act as a molecular link between Ca(2+) and ROS signal pathways in early fiber development.

Epidermal cell differentiation in cotton mediated by the homeodomain leucine zipper gene, GhHD-1

Gossypium hirsutum L. (cotton) fibres are specialized trichomes a few centimetres in length that grow from the seed coat. Few genes directly involved in the differentiation of these epidermal cells have been identified. These include GhMYB25-like and GhMYB25, two related MYB transcription factors that regulate fibre cell initiation and expansion. We have also identified a putative homeodomain leucine zipper (HD-ZIP) transcription factor, GhHD-1, expressed in trichomes and early fibres that might play a role in cotton fibre initiation. Here, we characterize GhHD-1 homoeologues from tetraploid G. hirsutum and show, using reporter constructs and quantitative real-time PCR (qRT-PCR), that they are expressed predominantly in epidermal tissues during early fibre development, and in other tissues bearing epidermal trichomes. Silencing of GhHD-1 reduced trichome formation and delayed the timing of fibre initiation. Constitutive overexpression of GhHD-1 increased the number of fibres initiating on the seed, but did not affect leaf trichomes. Expression of GhHD-1 in cotton silenced for different fibre MYBs suggest that in ovules it acts downstream of GhMYB25-like, but is unaffected in GhMYB25- or GhMYB109-silenced plants. Microarray analysis of silencing and overexpression lines of GhHD-1 indicated that it potentially regulates the levels of ethylene and reactive oxidation species (ROS) through a WRKY transcription factor and calcium-signalling pathway genes to activate downstream genes necessary for cell expansion and elongation.

Generation, annotation and analysis of first large-scale expressed sequence tags from developing fiber of Gossypium barbadense L

Background

Cotton fiber is the world's leading natural fiber used in the manufacture of textiles. Gossypium is also the model plant in the study of polyploidization, evolution, cell elongation, cell wall development, and cellulose biosynthesis. G. barbadense L. is an ideal candidate for providing new genetic variations useful to improve fiber quality for its superior properties. However, little is known about fiber development mechanisms of G. barbadense and only a few molecular resources are available in GenBank.

Methodology and Principal Findings

In total, 10,979 high-quality expressed sequence tags (ESTs) were generated from a normalized fiber cDNA library of G. barbadense. The ESTs were clustered and assembled into 5852 unigenes, consisting of 1492 contigs and 4360 singletons. The blastx result showed 2165 unigenes with significant similarity to known genes and 2687 unigenes with significant similarity to genes of predicted proteins. Functional classification revealed that unigenes were abundant in the functions of binding, catalytic activity, and metabolic pathways of carbohydrate, amino acid, energy, and lipids. The function motif/domain-related cytoskeleton and redox homeostasis were enriched. Among the 5852 unigenes, 282 and 736 unigenes were identified as potential cell wall biosynthesis and transcription factors, respectively. Furthermore, the relationships among cotton species or between cotton and other model plant systems were analyzed. Some putative species-specific unigenes of G. barbadense were highlighted.

Conclusions/Significance

The ESTs generated in this study are from the first large-scale EST project for G. barbadense and significantly enhance the number of G. barbadense ESTs in public databases. This knowledge will contribute to cotton improvements by studying fiber development mechanisms of G. barbadense, establishing a breeding program using marker-assisted selection, and discovering candidate genes related to important agronomic traits of cotton through oligonucleotide array. Our work will also provide important resources for comparative genomics, polyploidization, and genome evolution among Gossypium species.