Cloning and characterization of a calcium dependent protein kinase gene associated with cotton fiber development

Calcium-Dependent Protein Kinase GhCDPK16 Exerts a Positive Regulatory Role in Enhancing Drought Tolerance in Cotton

Cotton is essential for the textile industry as a primary source of natural fibers. However, environmental factors like drought present significant challenges to its cultivation, adversely affecting both production levels and fiber quality. Enhancing cotton's drought resilience has the potential to reduce yield losses and support the growth of cotton farming. In this study, the cotton calcium-dependent protein kinase GhCDPK16 was characterized, and the transcription level of GhCDPK16 was significantly upregulated under drought and various stress-related hormone treatments. Physiological analyses revealed that the overexpression of GhCDPK16 improved drought stress resistance in Arabidopsis by enhancing osmotic adjustment capacity and boosting antioxidant enzyme activities. In contrast, silencing GhCDPK16 in cotton resulted in increased dehydration compared with the control. Furthermore, reduced antioxidant enzyme activities and downregulation of ABA-related genes were observed in GhCDPK16-silenced plants. These findings not only enhanced our understanding of the biological functions of GhCDPK16 and the mechanisms underlying drought stress resistance but also underscored the considerable potential of GhCDPK16 in improving drought resilience in cotton.

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.

Calcium-dependent protein kinase GhCDPK4 plays a role in drought and abscisic acid stress responses

Drought is an important factor limiting the yield and quality of cotton. In the present study, the gene encoding the cotton calcium-dependent protein kinase GhCDPK4 was identified and characterized in the transcriptome of cotton under PEG-induced drought stress. In RT-qPCR experiments, GhCDPK4 expression was found to be up-regulated under drought and abscisic acid (ABA) stress. Under drought conditions, heterologous overexpression of GhCDPK4 in tobacco showed a better phenotypic status, higher antioxidant enzyme activity, and lower relative electrolyte leakage (REL) and malondialdehyde (MDA) content. Meanwhile, ghcdpk4-silenced cotton plants, which were extremely sensitive to drought, exhibited higher levels of O^2-，H₂O₂, and MDA contents compared to the control. Meanwhile, silenced lines showed impaired stomatal closure under drought stress, resulting in increased water loss from transpiration in silenced lines. GhCDPK4 expression was induced by ABA, and there are five ABA-responsive elements in its promoter. and C2-DOMAIN ABA-RELATED 4(CAR4, Gh_D09G1653) were found to interact and be co-expressed in the GhCDPK4 interaction network. Therefore, GhCDPK4 may reduce the extent of water loss and oxidative damage in cotton under drought by positively regulating ABA-controlled stomatal closure and reactive oxygen species (ROS) scavenging systems. This study demonstrates the great potential of GhCDPK4 in improving drought resistance in crops.

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.

Subgenome Bias and Temporal Postponement of Gene Expression Contributes to the Distinctions of Fiber Quality in Gossypium Species

As two cultivated widely allotetraploid cotton species, although Gossypium hirsutum and Gossypium barbadense evolved from the same ancestor, they differ in fiber quality; the molecular mechanism of that difference should be deeply studied. Here, we performed RNA-seq of fiber samples from four G. hirsutum and three G. barbadense cultivars to compare their gene expression patterns on multiple dimensions. We found that 15.90-37.96% of differentially expressed genes showed biased expression toward the A or D subgenome. In particular, interspecific biased expression was exhibited by a total of 330 and 486 gene pairs at 10 days post-anthesis (DPA) and 20 DPA, respectively. Moreover, 6791 genes demonstrated temporal differences in expression, including 346 genes predominantly expressed at 10 DPA in G. hirsutum (TM-1) but postponed to 20 DPA in G. barbadense (Hai7124), and 367 genes predominantly expressed at 20 DPA in TM-1 but postponed to 25 DPA in Hai7124. These postponed genes mainly participated in carbohydrate metabolism, lipid metabolism, plant hormone signal transduction, and starch and sucrose metabolism. In addition, most of the co-expression network and hub genes involved in fiber development showed asymmetric expression between TM-1 and Hai7124, like three hub genes detected at 10 DPA in TM-1 but not until 25 DPA in Hai7124. Our study provides new insights into interspecific expression bias and postponed expression of genes associated with fiber quality, which are mainly tied to asymmetric hub gene network. This work will facilitate further research aimed at understanding the mechanisms underlying cotton fiber improvement.

Advances about the Roles of Membranes in Cotton Fiber Development

Cotton fiber is an extremely elongated single cell derived from the ovule epidermis and is an ideal model for studying cell development. The plasma membrane is tremendously expanded and accompanied by the coordination of various physiological and biochemical activities on the membrane, one of the three major systems of a eukaryotic cell. This review compiles the recent progress and advances for the roles of the membrane in cotton fiber development: the functions of membrane lipids, especially the fatty acids, sphingolipids, and phytosterols; membrane channels, including aquaporins, the ATP-binding cassette (ABC) transporters, vacuolar invertase, and plasmodesmata; and the regulation mechanism of membrane proteins, such as membrane binding enzymes, annexins, and receptor-like kinases.

Genetic mapping and transcriptomic characterization of a new fuzzless-tufted cottonseed mutant

Fiber mutants are unique and valuable resources for understanding the genetic and molecular mechanisms controlling initiation and development of cotton fibers that are extremely elongated single epidermal cells protruding from the seed coat of cottonseeds. In this study, we reported a new fuzzless-tufted cotton mutant (Gossypium hirsutum) and showed that fuzzless-tufted near-isogenic lines (NILs) had similar agronomic traits and a higher ginning efficiency compared to their recurrent parents with normal fuzzy seeds. Genetic analysis revealed that the mutant phenotype is determined by a single incomplete dominant locus, designated N5. The mutation was fine mapped to an approximately 250-kb interval containing 33 annotated genes using a combination of bulked segregant sequencing, SNP chip genotyping, and fine mapping. Comparative transcriptomic analysis using 0-6 days post-anthesis (dpa) ovules from NILs segregating for the phenotypes of fuzzless-tufted (mutant) and normal fuzzy cottonseeds (wild-type) uncovered candidate genes responsible for the mutant phenotype. It also revealed that the flanking region of the N5 locus is enriched with differentially expressed genes (DEGs) between the mutant and wild-type. Several of those DEGs are members of the gene families with demonstrated roles in cell initiation and elongation, such as calcium-dependent protein kinase and expansin. The transcriptome landscape of the mutant was significantly reprogrammed in the 6 dpa ovules and, to a less extent, in the 0 dpa ovules, but not in the 2 and 4 dpa ovules. At both 0 and 6 dpa, the reprogrammed mutant transcriptome was mainly associated with cell wall modifications and transmembrane transportation, while transcription factor activity was significantly altered in the 6 dpa mutant ovules. These results imply a similar molecular basis for initiation of lint and fuzz fibers despite certain differences.

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

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

Comparative phosphoproteomic analysis of BR-defective mutant reveals a key role of GhSK13 in regulating cotton fiber development

Brassinosteroid (BR), a steroid phytohormone, whose signaling transduction pathways include a series of phosphorylation and dephosphorylation events, and GSK3s are the main negative regulator kinases. BRs have been shown to play vital roles in cotton fiber elongation. However, the underlying mechanism is still elusive. In this study, fibers of a BR-defective mutant Pagoda 1 (pag1), and its corresponding wild-type (ZM24) were selected for a comparative global phosphoproteome analysis at critical developmental time points: fast-growing stage (10 days after pollination (DPA)) and secondary cell wall synthesis stage (20 DPA). Based on the substrate characteristics of GSK3, 900 potential substrates were identified. Their GO and KEGG annotation results suggest that BR functions in fiber development by regulating GhSKs (GSK3s of Gossypium hirsutum L.) involved microtubule cytoskeleton organization, and pathways of glucose, sucrose and lipid metabolism. Further experimental results revealed that among the GhSK members identified, GhSK13 not only plays a role in BR signaling pathway, but also functions in developing fiber by respectively interacting with an AP2-like ethylene-responsive factor GhAP2L, a nuclear transcription factor Gh_DNF_YB19, and a homeodomain zipper member GhHDZ5. Overall, our phosphoproteomic research advances the understanding of fiber development controlled by BR signal pathways especially through GhSKs, and also offers numbers of target proteins for improving cotton fiber quality.

Insights into Drought Stress Signaling in Plants and the Molecular Genetic Basis of Cotton Drought Tolerance

Drought stress restricts plant growth and development by altering metabolic activity and biological functions. However, plants have evolved several cellular and molecular mechanisms to overcome drought stress. Drought tolerance is a multiplex trait involving the activation of signaling mechanisms and differentially expressed molecular responses. Broadly, drought tolerance comprises two steps: stress sensing/signaling and activation of various parallel stress responses (including physiological, molecular, and biochemical mechanisms) in plants. At the cellular level, drought induces oxidative stress by overproduction of reactive oxygen species (ROS), ultimately causing the cell membrane to rupture and stimulating various stress signaling pathways (ROS, mitogen-activated-protein-kinase, Ca2+, and hormone-mediated signaling). Drought-induced transcription factors activation and abscisic acid concentration co-ordinate the stress signaling and responses in cotton. The key responses against drought stress, are root development, stomatal closure, photosynthesis, hormone production, and ROS scavenging. The genetic basis, quantitative trait loci and genes of cotton drought tolerance are presented as examples of genetic resources in plants. Sustainable genetic improvements could be achieved through functional genomic approaches and genome modification techniques such as the CRISPR/Cas9 system aid the characterization of genes, sorted out from stress-related candidate single nucleotide polymorphisms, quantitative trait loci, and genes. Exploration of the genetic basis for superior candidate genes linked to stress physiology can be facilitated by integrated functional genomic approaches. We propose a third-generation sequencing approach coupled with genome-wide studies and functional genomic tools, including a comparative sequenced data (transcriptomics, proteomics, and epigenomic) analysis, which offer a platform to identify and characterize novel genes. This will provide information for better understanding the complex stress cellular biology of plants.

Proteomic profiling of cotton fiber developmental transition from cell elongation to secondary wall deposition

Cotton fiber developmental transition from elongation to secondary cell wall biosynthesis is a critical growth shifting phase that affects fiber final length, strength, and other properties. Morphological dynamic analysis indicated an asynchronous fiber developmental pattern between two most important commercial cotton species, Gossypium hirsutum (Gh) and G. barbadense (Gb). Using isobaric tags for relative and absolute quantitation techniques, we examined the temporal changes of protein expression at three representative development periods (15-19, 19-23, and 23-27 dpa) in both species. Strikingly, a large proportion of differentially expressed proteins (DEPs) were identified at 19-23 dpa in Gh and at 23-27 dpa in Gb, corresponding to their fiber developmental transition timing. To better understand fiber transitional development, we comparatively analyzed those DEPs in 19-23 dpa of Gh vs. in 23-27 dpa of Gb, and noted that these cotton species indeed share fundamentally similar fiber developmental features under the biological processes. We also showed that there are limited overlaps in both specific upregulated and downregulated proteins between the two species, suggesting species-specific protein regulations in the development process. Proteomic profiling results revealed dynamic changes of several key proteins and biological processes that are potentially correlated with fiber developmental transition. During the transition, upregulated proteins are mainly involved in carbohydrate/energy metabolism, oxidation-reduction, cytoskeleton, protein turnover, Ca2+ signaling, etc., whereas important downregulated proteins are mostly involved in phenylpropanoid and flavonoid secondary metabolism pathways. The gene expressions of several changed proteins in this key stage were also examined by quantitative reverse transcription polymerase chain reaction. Overall, the present study provides accurate pictures of the regulatory networks of functional proteins during the fiber developmental transition, therefore highlighting candidate genes/proteins and related pathways for the cotton fiber improvement.

Identification of Loci and Candidate Genes Responsible for Fiber Length in Upland Cotton (Gossypium hirsutum L.) via Association Mapping and Linkage Analyses

Fiber length (FL) is an important fiber quality trait in cotton. Although many fiber quality quantitative trait loci (QTL) responsible for FL have been identified, most cannot be applied to breeding programs, mainly due to unstable environments or large confidence intervals. In this study, we combined a genome-wide association study (GWAS) and linkage mapping to identify and validate high-quality QTLs responsible for FL. For the GWAS, we developed 93,250 high-quality single-nucleotide polymorphism (SNP) markers based on 355 accessions, and the FL was measured in eight different environments. For the linkage mapping, we constructed an F ₂ population from two extreme accessions. The high-density linkage maps spanned 3,848.29 cM, with an average marker interval of 1.41 cM. In total, 14 and 13 QTLs were identified in the association and linkage mapping analyses, respectively. Most importantly, a major QTL on chromosome D03 identified in both populations explained more than 10% of the phenotypic variation (PV). Furthermore, we found that a sucrose synthesis-related gene (Gh_D03G1338) was associated with FL in this QTL region. The RNA-seq data showed that Gh_D03G1338 was highly expressed during the fiber development stage, and the qRT-PCR analysis showed significant expression differences between the long fiber and short fiber varieties. These results suggest that Gh_D03G1338 may determine cotton fiber elongation by regulating the synthesis of sucrose. Favorable QTLs and candidate genes should be useful for increasing fiber quality in cotton breeding.

YC3.60-Based Imaging Analysis on Calcium Level in Cotton Cells

Calcium ion (Ca²⁺) is a core regulator of cell functions in response to many developmental and environmental stimuli. A hallmark for Ca²⁺ signaling is the change of this ion in cells. Fluorescent resonance energy transfer (FRET)-based Ca²⁺ sensors provide a powerful tool for qualitatively and quantitatively measuring cytosolic Ca²⁺ level. Using YC3.60, one of those sensors, we have imaged cytosolic changes of Ca²⁺ in cotton fibers during the initiation stage. In this chapter, the imaging method is described in detail. The description is not limited to fiber cells but also examples leaf trichomes and protoplasts of cotton.

Coping with drought: stress and adaptive mechanisms, and management through cultural and molecular alternatives in cotton as vital constituents for plant stress resilience and fitness

Increased levels of greenhouse gases in the atmosphere and associated climatic variability is primarily responsible for inducing heat waves, flooding and drought stress. Among these, water scarcity is a major limitation to crop productivity. Water stress can severely reduce crop yield and both the severity and duration of the stress are critical. Water availability is a key driver for sustainable cotton production and its limitations can adversely affect physiological and biochemical processes of plants, leading towards lint yield reduction. Adaptation of crop husbandry techniques suitable for cotton crop requires a sound understanding of environmental factors, influencing cotton lint yield and fiber quality. Various defense mechanisms e.g. maintenance of membrane stability, carbon fixation rate, hormone regulation, generation of antioxidants and induction of stress proteins have been found play a vital role in plant survival under moisture stress. Plant molecular breeding plays a functional role to ascertain superior genes for important traits and can offer breeder ready markers for developing ideotypes. This review highlights drought-induced damage to cotton plants at structural, physiological and molecular levels. It also discusses the opportunities for increasing drought tolerance in cotton either through modern gene editing technology like clustered regularly interspaced short palindromic repeat (CRISPR/Cas9), zinc finger nuclease, molecular breeding as well as through crop management, such as use of appropriate fertilization, growth regulator application and soil amendments.

GWAS Analysis and QTL Identification of Fiber Quality Traits and Yield Components in Upland Cotton Using Enriched High-Density SNP Markers

It is of great importance to identify quantitative trait loci (QTL) controlling fiber quality traits and yield components for future marker-assisted selection (MAS) and candidate gene function identifications. In this study, two kinds of traits in 231 F6:8 recombinant inbred lines (RILs), derived from an intraspecific cross between Xinluzao24, a cultivar with elite fiber quality, and Lumianyan28, a cultivar with wide adaptability and high yield potential, were measured in nine environments. This RIL population was genotyped by 122 SSR and 4729 SNP markers, which were also used to construct the genetic map. The map covered 2477.99 cM of hirsutum genome, with an average marker interval of 0.51 cM between adjacent markers. As a result, a total of 134 QTLs for fiber quality traits and 122 QTLs for yield components were detected, with 2.18-24.45 and 1.68-28.27% proportions of the phenotypic variance explained by each QTL, respectively. Among these QTLs, 57 were detected in at least two environments, named stable QTLs. A total of 209 and 139 quantitative trait nucleotides (QTNs) were associated with fiber quality traits and yield components by four multilocus genome-wide association studies methods, respectively. Among these QTNs, 74 were detected by at least two algorithms or in two environments. The candidate genes harbored by 57 stable QTLs were compared with the ones associated with QTN, and 35 common candidate genes were found. Among these common candidate genes, four were possibly "pleiotropic." This study provided important information for MAS and candidate gene functional studies.

Genome-wide classification, evolutionary analysis and gene expression patterns of the kinome in Gossypium

The protein kinase (PK, kinome) family is one of the largest families in plants and regulates almost all aspects of plant processes, including plant development and stress responses. Despite their important functions, comprehensive functional classification, evolutionary analysis and expression patterns of the cotton PK gene family has yet to be performed on PK genes. In this study, we identified the cotton kinomes in the Gossypium raimondii, Gossypium arboretum, Gossypium hirsutum and Gossypium barbadense genomes and classified them into 7 groups and 122-24 subfamilies using software HMMER v3.0 scanning and neighbor-joining (NJ) phylogenetic analysis. Some conserved exon-intron structures were identified not only in cotton species but also in primitive plants, ferns and moss, suggesting the significant function and ancient origination of these PK genes. Collinearity analysis revealed that 16.6 million years ago (Mya) cotton-specific whole genome duplication (WGD) events may have played a partial role in the expansion of the cotton kinomes, whereas tandem duplication (TD) events mainly contributed to the expansion of the cotton RLK group. Synteny analysis revealed that tetraploidization of G. hirsutum and G. barbadense contributed to the expansion of G. hirsutum and G. barbadense PKs. Global expression analysis of cotton PKs revealed stress-specific and fiber development-related expression patterns, suggesting that many cotton PKs might be involved in the regulation of the stress response and fiber development processes. This study provides foundational information for further studies on the evolution and molecular function of cotton PKs.

Recent insights into cotton functional genomics: progress and future perspectives

Functional genomics has transformed from futuristic concept to well-established scientific discipline during the last decade. Cotton functional genomics promise to enhance the understanding of fundamental plant biology to systematically exploit genetic resources for the improvement of cotton fibre quality and yield, as well as utilization of genetic information for germplasm improvement. However, determining the cotton gene functions is a much more challenging task, which has not progressed at a rapid pace. This article presents a comprehensive overview of the recent tools and resources available with the major advances in cotton functional genomics to develop elite cotton genotypes. This effort ultimately helps to filter a subset of genes that can be used to assemble a final list of candidate genes that could be employed in future novel cotton breeding programme. We argue that next stage of cotton functional genomics requires the draft genomes refinement, re-sequencing broad diversity panels with the development of high-throughput functional genomics tools and integrating multidisciplinary approaches in upcoming cotton improvement programmes.

Calcium-dependent protein kinases in cotton: insights into early plant responses to salt stress

BACKGROUND

Soil salinization is one of the major environmental constraints to plant growth and agricultural production worldwide. Signaling components involving calcium (Ca²⁺) and the downstream calcium-dependent protein kinases (CPKs) play key roles in the perception and transduction of stress signals. However, the study of CPKs in cotton and their functions in response to salt stress remain unexplored.

RESULTS

A total of 98 predicted CPKs were identified from upland cotton (Gossypium hirsutum L. 'TM-1'), and phylogenetic analyses classified them into four groups. Gene family distribution studies have revealed the substantial impacts of the genome duplication events to the total number of GhCPKs. Transcriptome analyses showed a wide distribution of CPKs' expression among different organs. A total of 19 CPKs were selected for their rapid responses to salt stress at the transcriptional level, most of which were also incduced by the thylene-releasing chemical ethephon, suggesting a partal overlap of the salinity and ethylene responses. Silencing of 4 of the 19 CPKs (GhCPK8, GhCPK38, GhCPK54, and GhCPK55) severely compromised the basal cotton resistance to salt stress.

CONCLUSIONS

Our genome-wide expression analysis of CPK genes from up-land cotton suggests that CPKs are involved in multiple developmental responses as well as the response to different abiotic stresses. A cluster of the cotton CPKs was shown to participate in the early signaling events in cotton responses to salt stress. Our results provide significant insights on functional analysis of CPKs in cotton, especially in the context of cotton adaptions to salt stress.

Interaction between calcium and potassium modulates elongation rate in cotton fiber cells

Calcium (Ca2+) is necessary for fiber cell development in cotton (Gossypium hirsutum), both as a cell wall structural component and for environmental signaling responses. It is also known that potassium (K+) plays a critical role in cotton fiber cell elongation. However, it is unclear whether Ca2+ integrates its activities with K+ to regulate fiber elongation. Here, we report the novel discovery that Ca2+ deficiency, when integrated with K+ signaling, promotes fiber elongation. Using inductively coupled plasma-mass spectrometry (ICP-MS), we determined dynamic profiles of the ionome in ovules and fibers at different developmental stages, and found that a high accumulation of macro-elements, but not Ca2+, was associated with longer fibers. Using an in vitro ovule culture system, we found that under Ca2+-deficient conditions, sufficient K+ (52 mM) rapidly induced ovule and fiber browning, while reduced K+ (2 or 27 mM) not only suppressed tissue browning but also altered fiber elongation. Reduced K+ also enhanced reactive oxygen species scavenging ability and maintained abscisic acid and jasmonic acid levels, which in turn compensated for Ca2+ deficiency. Ca2+ deficiency combined with reduced K+ (0 mM Ca2+ and 27 mM K+) produced longer fibers in cultured ovules, due to cell wall loosening by phytosulfokine (PSK), expansin (EXP), and xyloglucan endotransglycosylase/hydrolase (XTH), and an increase of the K+ content of fiber cells. Using transgenic cotton, we showed that the CBL-INTERACTING PROTEIN KINASE 6 (GhCIPK6) gene mediates the uptake of K+ under Ca2+-deficient conditions. This study establishes a new link between Ca2+, K+, and fiber elongation.

In vivo imaging of Ca2+ accumulation during cotton fiber initiation using fluorescent indicator YC3.60

Non-tip-focused Ca ²⁺ gradient indicated by genetically expressing a FRET-based calcium sensor YC3.60 was established in spherical expanding cotton fibers, which is vital for cotton fiber initiation. Cotton fiber is a single cell elongated from ovule epidermis. It is not only the most important natural fiber used in the textile industry but also an ideal model for studying cell differentiation and elongation. Before linear cell growth, cotton fibers undergo spherical expansion at the beginning of initiation. Ca²⁺, as an important secondary messenger, plays a central role in polarized cell growth including cotton fiber elongation. However, the role of Ca²⁺ in fiber initiation is far from well understood. In this paper, through ovule culture we demonstrate that Ca²⁺ is crucial for fiber initiation. Using transgenic cotton expressing the fluorescent Ca²⁺ indicator YC3.60, we show cellular and intracellular distribution of Ca²⁺ in cotton ovule epidermis and fiber cells. In the initiating fiber cell, Ca²⁺ accumulated mainly at the base of the cell, while in the fast elongating cell, the Ca²⁺ was enriched in the tip region. This cellular distribution of Ca²⁺ reported by YC3.60 was confirmed by the staining with a Ca²⁺-sensitive dye fluo-3/AM. Compared to the fluorescent dye staining, the YC3.60 system can reveal more detailed information on the intracellular distribution without photobleaching. Taken together, our data suggest that Ca²⁺ plays an important role in spherical expansion of cotton fiber initials.

MicroRNA expression profiles during cotton (Gossypium hirsutum L) fiber early development

The role of microRNAs (miRNAs) during cotton fiber development remains unclear. Here, a total of 54 miRNAs belonging to 39 families were selected to characterize miRNA regulatory mechanism in eight different fiber development stages in upland cotton cv BM-1. Among 54 miRNAs, 18 miRNAs were involved in cotton fiber initiation and eight miRNAs were related to fiber elongation and secondary wall biosynthesis. Additionally, 3,576 protein-coding genes were candidate target genes of these miRNAs, which are potentially involved in cotton fiber development. We also investigated the regulatory network of miRNAs and corresponding targets in fiber initiation and elongation, and secondary wall formation. Our Gene Ontology-based term classification and KEGG-based pathway enrichment analyses showed that the miRNA targets covered 220 biological processes, 67 molecular functions, 45 cellular components, and 10 KEGG pathways. Three of ten KEGG pathways were involved in lignan synthesis, cell elongation, and fatty acid biosynthesis, all of which have important roles in fiber development. Overall, our study shows the potential regulatory roles of miRNAs in cotton fiber development and the importance of miRNAs in regulating different cell types. This is helpful to design miRNA-based biotechnology for improving fiber quality and yield.

Drought coping strategies in cotton: increased crop per drop

The growth and yield of many crops, including cotton, are affected by water deficit. Cotton has evolved drought specific as well as general morpho-physiological, biochemical and molecular responses to drought stress, which are discussed in this review. The key physiological responses against drought stress in cotton, including stomata closing, root development, cellular adaptations, photosynthesis, abscisic acid (ABA) and jasmonic acid (JA) production and reactive oxygen species (ROS) scavenging, have been identified by researchers. Drought stress induces the expression of stress-related transcription factors and genes, such as ROS scavenging, ABA or mitogen-activated protein kinases (MAPK) signalling genes, which activate various drought-related pathways to induce tolerance in the plant. It is crucial to elucidate and induce drought-tolerant traits via quantitative trait loci (QTL) analysis, transgenic approaches and exogenous application of substances. The current review article highlights the natural as well as engineered drought tolerance strategies in cotton.

Comprehensive Analysis of the CDPK-SnRK Superfamily Genes in Chinese Cabbage and Its Evolutionary Implications in Plants

The CDPK-SnRK (calcium-dependent protein kinase/Snf1-related protein kinase) gene superfamily plays important roles in signaling pathways for disease resistance and various stress responses, as indicated by emerging evidence. In this study, we constructed comparative analyses of gene structure, retention, expansion, whole-genome duplication (WGD) and expression patterns of CDPK-SnRK genes in Brassica rapa and their evolution in plants. A total of 49 BrCPKs, 14 BrCRKs, 3 BrPPCKs, 5 BrPEPRKs, and 56 BrSnRKs were identified in B. rapa. All BrCDPK-SnRK proteins had highly conserved kinase domains. By statistical analysis of the number of CDPK-SnRK genes in each species, we found that the expansion of the CDPK-SnRK gene family started from angiosperms. Segmental duplication played a predominant role in CDPK-SnRK gene expansion. The analysis showed that PEPRK was more preferentially retained than other subfamilies and that CPK was retained similarly to SnRK. Among the CPKs and SnRKs, CPKIII and SnRK1 genes were more preferentially retained than other groups. CRK was closest to CPK, which may share a common evolutionary origin. In addition, we identified 196 CPK genes and 252 SnRK genes in 6 species, and their different expansion and evolution types were discovered. Furthermore, the expression of BrCDPK-SnRK genes is dynamic in different tissues as well as in response to abiotic stresses, demonstrating their important roles in development in B. rapa. In summary, this study provides genome-wide insight into the evolutionary history and mechanisms of CDPK-SnRK genes following whole-genome triplication in B. rapa.

Solanum tuberosum StCDPK1 is regulated by miR390 at the posttranscriptional level and phosphorylates the auxin efflux carrier StPIN4 in vitro, a potential downstream target in potato development

Among many factors that regulate potato tuberization, calcium and calcium-dependent protein kinases (CDPKs) play an important role. CDPK activity increases at the onset of tuber formation with StCDPK1 expression being strongly induced in swollen stolons. However, not much is known about the transcriptional and posttranscriptional regulation of StCDPK1 or its downstream targets in potato development. To elucidate further, we analyzed its expression in different tissues and stages of the life cycle. Histochemical analysis of StCDPK1::GUS (β-glucuronidase) plants demonstrated that StCDPK1 is strongly associated with the vascular system in stems, roots, during stolon to tuber transition, and in tuber sprouts. In agreement with the observed GUS profile, we found specific cis-acting elements in StCDPK1 promoter. In silico analysis predicted miR390 to be a putative posttranscriptional regulator of StCDPK1. Quantitative real time-polymerase chain reaction (qRT-PCR) analysis showed ubiquitous expression of StCDPK1 in different tissues which correlated well with Western blot data except in leaves. On the contrary, miR390 expression exhibited an inverse pattern in leaves and tuber eyes suggesting a possible regulation of StCDPK1 by miR390. This was further confirmed by Agrobacterium co-infiltration assays. In addition, in vitro assays showed that recombinant StCDPK1-6xHis was able to phosphorylate the hydrophilic loop of the auxin efflux carrier StPIN4. Altogether, these results indicate that StCDPK1 expression is varied in a tissue-specific manner having significant expression in vasculature and in tuber eyes; is regulated by miR390 at posttranscriptional level and suggest that StPIN4 could be one of its downstream targets revealing the overall role of this kinase in potato development.

GhCaM7-like, a calcium sensor gene, influences cotton fiber elongation and biomass production

Calcium signaling regulates many developmental processes in plants. Calmodulin (CaM) is one of the most conserved calcium sensors and has a flexible conformation in eukaryotes. The molecular functions of CaM are unknown in cotton, which is a major source of natural fiber. In this study, a Gossypium hirsutum L.CaM7-like gene was isolated from upland cotton. Bioinformatics analysis indicated that the GhCaM7-like gene was highly conserved as compared with Arabidopsis AtCaM7. The GhCaM7-like gene showed a high expression level in elongating fibers. Expression of β-glucuronidase was observed in trichomes on the stem, leaf and root in transgenic Arabidopsis plants of a PRO_GhCaM7-like:GUS fusion. Silencing of the GhCaM7-like gene resulted in decreased fiber length, but also caused reduction in stem height, leaf dimensions, seed length and 100-seed weight, in comparison with those of the control. Reduced expression of the GhCaM7-like gene caused decreased Ca²⁺ influx in cells of the leaf hypodermis and stem apex, and down-regulation of GhIQD1 (IQ67-domain containing protein), GhAnn2 (Annexins) and GhEXP2 (Expansin). These results indicate that the GhCaM7-like gene plays a vital role in calcium signaling pathways, and may regulate cotton fiber elongation and biomass production by affecting Ca²⁺ signatures and downstream signaling pathways of CaM.

Know where your clients are: subcellular localization and targets of calcium-dependent protein kinases

Calcium-dependent protein kinases (CDPKs) are at the forefront of decoding transient Ca(2+) signals into physiological responses. They play a pivotal role in many aspects of plant life starting from pollen tube growth, during plant development, and in stress response to senescence and cell death. At the cellular level, Ca(2+) signals have a distinct, narrow distribution, thus requiring a conjoined localization of the decoders. Accordingly, most CDPKs have a distinct subcellular distribution which enables them to 'sense' the local Ca(2+) concentration and to interact specifically with their targets. Here we present a comprehensive overview of identified CDPK targets and discuss them in the context of kinase-substrate specificity and subcellular distribution of the CDPKs. This is particularly relevant for calcium-mediated phosphorylation where different CDPKs, as well as other kinases, were frequently reported to be involved in the regulation of the same target. However, often these studies were not performed in an in situ context. Thus, considering the specific expression patterns, distinct subcellular distribution, and different Ca(2+) affinities of CDPKs will narrow down the number of potential CDPKs for one given target. A number of aspects still remain unresolved, giving rise to pending questions for future research.

Detection and validation of one stable fiber strength QTL on c9 in tetraploid cotton

Fiber strength is an essential trait of fiber property in cotton, and it is quantitatively inherited. Identification of stable quantitative trait loci (QTL) contributing to fiber strength would provide the key basis for marker-assisted selection (MAS) in cotton breeding. In this study, four interspecific hybridization populations were established with a common G. barbadense parent Pima 90-53 and two G. hirsutum parents (CCRI 8 and Handan 208), each of which had fiber strength characteristic. Based on the phenotypic data of fiber strength from seven environments, a stable QTL, qFS-c9-1, was detected and validated on c9 in a marker interval between SSR markers NAU2395 and NAU1092. The QTL explaining 14.4-17.9 % of the phenotypic variation was firstly detected in two populations (CCRI 8 × Pima 90-53, BC1F1 and BC1F2) and its derived lines in four environments. And it accounting for 12.1-14.8 % of the phenotypic variation was further confirmed in two populations (Handan 208 × Pima 90-53, BC1F1, and F2) under one environment. In silico mapping using three sequenced cotton genomes indicated that homologous genes, anchored by NAU2395 and NAU1092, were aligned to the G. arboreum genome within a physical distance between 81.10 Mbps and 87.07 Mbps. In that interval, several genes were confirmed in literatures to associate with fiber development. Among these genes, seven genes were further selected for an expression analysis through fiber development transcriptome database, revealing unique expression patterns across different stages of fiber development between CCRI 8 and Pima 90-53. The genes underlying qFS-c9-1 were favorable to fine mapping and cloning. The current study results provided valuable evidence for mapping stable QTL of fiber strength utilizing multiple populations and environments, as well as map-based cloning the candidate gene underlying the QTL for future prospective research directions.

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.

Down-regulating annexin gene GhAnn2 inhibits cotton fiber elongation and decreases Ca2+ influx at the cell apex

Cotton fiber is a single cell that differentiates from the ovule epidermis and undergoes synchronous elongation with high secretion and growth rate. Apart from economic importance, cotton fiber provides an excellent single-celled model for studying mechanisms of cell-growth. Annexins are Ca(2+)- and phospholipid-binding proteins that have been reported to be localized in multiple cellular compartments and involved in control of vesicle secretions. Although several annexins have been found to be highly expressed in elongating cotton fibers, their functional roles in fiber development remain unknown. Here, 14 annexin family members were identified from the fully sequenced diploid G. raimondii (D5 genome), half of which were expressed in fibers of the cultivated tetraploid species G. hirsutum (cv. YZ1). Among them, GhAnn2 from the D genome of the tetraploid species displayed high expression level in elongating fiber. The expression of GhAnn2 could be induced by some phytohormones that play important roles in fiber elongation, such as IAA and GA3. RNAi-mediated down-regulation of GhAnn2 inhibited fiber elongation and secondary cell wall synthesis, resulting in shorter and thinner mature fibers in the transgenic plants. Measurement with non-invasive scanning ion-selective electrode revealed that the rate of Ca(2+) influx from extracellular to intracellular was decreased at the fiber cell apex of GhAnn2 silencing lines, in comparison to that in the wild type. These results indicate that GhAnn2 may regulate fiber development through modulating Ca(2+) fluxes and signaling.

PAG1, a cotton brassinosteroid catabolism gene, modulates fiber elongation

Cotton (Gossypium hirsutum) is the major source of natural textile fibers. Brassinosteroids (BRs) play crucial roles in regulating fiber development. The molecular mechanisms of BRs in regulating fiber elongation, however, are poorly understood. pagoda1 (pag1) was identified via an activation tagging genetic screen and characterized by genome walking and brassinolide (BL) supplementation. RNA-Seq analysis was employed to elucidate the mechanisms of PAG1 in regulating fiber development. pag1 exhibited dwarfism and reduced fiber length due to significant inhibition of cell elongation and expansion. BL treatment rescued its growth and fiber elongation. PAG1 encodes a homolog of Arabidopsis CYP734A1 that inactivates BRs via C-26 hydroxylation. RNA-Seq analyses showed that the constitutive expression of PAG1 downregulated the expression of genes involved in very-long-chain fatty acids (VLCFA) biosynthesis, ethylene-mediated signaling, response to cadmium, cell wall development, cytoskeleton organization and cell growth. Our results demonstrate that PAG1 plays crucial roles in regulating fiber development via controlling the level of endogenous bioactive BRs, which may affect ethylene signaling cascade by mediating VLCFA. Therefore, BR may be a critical regulator of fiber elongation, a role which may in turn be linked to effects on VLCFA biosynthesis, ethylene and cadmium signaling, cell wall- and cytoskeleton-related gene expression.

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.

Genome-wide survey and expression analysis of calcium-dependent protein kinase in Gossypium raimondii

Calcium-dependent protein kinases (CDPKs) are one of the largest protein kinases in plants and participate in different physiological processes through regulating downstream components of calcium signaling pathways. In this study, 41 CDPK genes, from GrCPK1 to GrCPK41, were identified in the genome of the diploid cotton, Gossypium raimondii. The phylogenetic analysis indicated that all these genes were divided into four subgroups and members within the same subgroup shared conserved exon-intron structures. The expansion of GrCPKs family in G. raimondii was due to the segmental duplication events, and the analysis of Ka/Ks ratios implied that the duplicated GrCPKs had mainly undergone strong purifying selection pressure with limited functional divergence. The cold-responsive elements in promoter regions were detected in the majority of GrCPKs. The expression analysis of 11 selected genes showed that GrCPKs exhibited tissue-specific expression patterns and the expression of GrCPKs were induced or repressed by cold treatment. These observations would lay an important foundation for functional and evolutionary analysis of CDPK gene family in Gossypium species.

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.

Genome-wide analysis of small RNAs reveals eight fiber elongation-related and 257 novel microRNAs in elongating cotton fiber cells

BACKGROUND

MicroRNAs (miRNAs) and other types of small regulatory RNAs play critical roles in the regulation of gene expression at the post-transcriptional level in plants. Cotton is one of the most economically important crops, but little is known about the roles of miRNAs during cotton fiber elongation.

RESULTS

Here, we combined high-throughput sequencing with computational analysis to identify small RNAs (sRNAs) related to cotton fiber elongation in Gossypium hirsutum L. (G. hirsutum). The sequence analysis confirmed the expression of 79 known miRNA families in elongating fiber cells and identified 257 novel miRNAs, primarily derived from corresponding specific loci in the Gossypium raimondii Ulbr. (G. raimondii) genome. Furthermore, a comparison of the miRNAomes revealed that 46 miRNA families were differentially expressed throughout the elongation period. Importantly, the predicted and experimentally validated targets of eight miRNAs were associated with fiber elongation, with obvious functional relationships with calcium and auxin signal transduction, fatty acid metabolism, anthocyanin synthesis and the xylem tissue differentiation. Moreover, one tasiRNA was also identified, and its target, ARF4, was experimentally validated in vivo.

CONCLUSION

This study not only facilitated the discovery of 257 novel low-abundance miRNAs in elongating cotton fiber cells but also revealed a potential regulatory network of nine sRNAs important for fiber elongation. The identification and characterization of miRNAs in elongating cotton fiber cells might promote the further study of fiber miRNA regulation mechanisms and provide insight into the importance of miRNAs in cotton.

A cotton annexin protein AnxGb6 regulates fiber elongation through its interaction with actin 1

Annexins are assumed to be involved in regulating cotton fiber elongation, but direct evidence remains to be presented. Here we cloned six Annexin genes (AnxGb) abundantly expressed in fiber from sea-island cotton (G. barbadense). qRT-PCR results indicated that all six G. barbadense annexin genes were expressed in elongating cotton fibers, while only the expression of AnxGb6 was cotton fiber-specific. Yeast two hybridization and BiFC analysis revealed that AnxGb6 homodimer interacted with a cotton fiber specific actin GbAct1. Ectopic-expressed AnxGb6 in Arabidopsis enhanced its root elongation without increasing the root cell number. Ectopic AnxGb6 expression resulted in more F-actin accumulation in the basal part of the root cell elongation zone. Analysis of AnxGb6 expression in three cotton genotypes with different fiber length confirmed that AnxGb6 expression was correlated to cotton fiber length, especially fiber elongation rate. Our results demonstrated that AnxGb6 was important for fiber elongation by potentially providing a domain for F-actin organization.

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.

Transcriptome analysis reveals critical genes and key pathways for early cotton fiber elongation in Ligon lintless-1 mutant

Fiber length is a key determinant of cotton yield and quality. Using a monogenic dominant cotton mutant Ligon lintless-1 with extremely short fibers, we employed microarray technology and quantitative real time PCR to compare transcriptomes of Li(1) and the normal wild-type TM-1, the results showed that only a few genes differentially expressed in 0 days postanthesis (DPA) ovules and 3 DPA fibers, whereas 577 transcripts differentially expressed in 6 DPA fibers. 6 DPA is probably a key phase determining fiber elongation. Gene ontology analyses showed such processes as response to stimulus, signal transduction, and lipid metabolism were readjusted by the mutant gene. Pathway studio analysis indicated that auxin signaling and sugar signaling pathways play major roles in modulation of early fiber elongation. This work provides new insight into the mechanisms of fiber development, and offers novel genes as potential objects for genetic manipulation to achieve improvement of fiber properties.

How cotton fibers elongate: a tale of linear cell-growth mode

Cotton fibers (cotton lint) are single-celled trichomes that differentiate from the ovule epidermis. Unidirectional and fast-growing cells generally expand at the dome-shaped apical zone (tip-growth mode); however, previous studies suggest that elongating fiber cells expand via a diffuse-growth mode. Tip-localized Ca(2+) gradient and active secretary vesicle trafficking are two important phenomena of tip-growth. Recently, a high Ca(2+) gradient is found in the cytoplasm of fast-elongating cotton fiber cells near the growing tip. Several protein coding genes participating in vesicle coating and transport are highly expressed in elongating fiber cells. Taken together with the observation that ethylene acts as a positive regulator for cotton fiber and several Arabidopsis tissues that are known to elongate via tip growth prompted us to propose a linear-growth mode for similar cell types.

Arabidopsis calcium-dependent protein kinase CPK10 functions in abscisic acid- and Ca2+-mediated stomatal regulation in response to drought stress

Plant calcium-dependent protein kinases (CDPKs) may function as calcium sensors and play important roles in the regulation of plant growth and development and in plant responses to biotic and abiotic stresses. The Arabidopsis (Arabidopsis thaliana) genome encodes 34 CDPKs, and most of them have not been functionally characterized. Here, we report the functional characterization of CPK10 in Arabidopsis response to drought stress. The cpk10 mutant, a T-DNA insertion mutant for the Arabidopsis CPK10 gene, showed a much more sensitive phenotype to drought stress compared with wild-type plants, while the CPK10 overexpression lines displayed enhanced tolerance to drought stress. Induction of stomatal closure and inhibition of stomatal opening by abscisic acid (ABA) and Ca(2+) were impaired in the cpk10 mutants. Using yeast two-hybrid methods, a heat shock protein, HSP1, was identified as a CPK10-interacting protein. The interaction between CPK10 and HSP1 was further confirmed by pull-down and bimolecular fluorescence complementation assays. The HSP1 knockout mutant (hsp1) plants showed a similar sensitive phenotype under drought stress as the cpk10 mutant plants and were similarly less sensitive to ABA and Ca(2+) in regulation of stomatal movements. Electrophysiological experiments showed that ABA and Ca(2+) inhibition of the inward K(+) currents in stomatal guard cells were impaired in the cpk10 and hsp1 mutants. All presented data demonstrate that CPK10, possibly by interacting with HSP1, plays important roles in ABA- and Ca(2+)-mediated regulation of stomatal movements.

Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant

Cotton fiber as a single-celled trichome is a biological model system for studying cell differentiation and elongation. However, the complexity of its gene expression and regulatory mechanism allows only marginal progress. Here, we report the high-throughput tag-sequencing (Tag-seq) analysis using Solexa Genome Analyzer platform on transcriptome of -2 to 1 (fiber initiation, stage I) and 2-8 (fiber elongation, stage II) days post anthesis (DPA) cotton (Gossypium hirsutum) ovules (wild type: WT; Xuzhou 142 and its mutant: fuzzless/lintless or flM, in the same background). To this end, we sequenced 3.5-3.8 million tags representing 0.7-1.0 million unique transcripts for each library (WT1, WT2, M1, and M2). After removal of low quality tags, we obtained a total of 2,973,104, 3,139,306, 2,943,654, and 3,392,103 clean sequences that corresponded to 357,852, 280,787, 372,952, and 382,503 distinct tags for WT1, WT2, M1, and M2, respectively. All clean tags were aligned to the publicly available cotton transcript database (TIGR, http://www.tigr.org). About 15% of the distinct tags were uniquely mapped to the reference genes, and 31.4% of existing genes were matched by tags. The tag mapping to the database sequences generated 23,854, 24,442, 23,497, and 19,957 annotated genes for WT1, WT2, M1, and M2 libraries, respectively. Analyses of differentially expressed genes revealed the substantial changes in gene type and abundance between the wild type and mutant libraries. Among the 20 most differentially expressed genes in WT1/M1 and WT2/M2 libraries were cellulose synthase, phosphatase, and dehydrogenase, all of which are involved in the fiber cell development. Overall, the deep-sequencing analyses demonstrate the high degree of transcriptional complexity in early developing fibers and represent a major improvement over the microarrays for analyzing transcriptional changes on a large scale.

Proteomic identification of differentially expressed proteins in the Ligon lintless mutant of upland cotton (Gossypium hirsutum L.)

Cotton fiber is an ideal model for studying plant cell elongation. To date, the underlying mechanisms controlling fiber elongation remain unclear due to their high complexity. In this study, a comparative proteomic analysis between a short-lint fiber mutant (Ligon lintless, Li(1)) and its wild-type was performed to identify fiber elongation-related proteins. By 2-DE combined with local EST database-assisted MS/MS analysis, 81 differentially expressed proteins assigned to different functional categories were identified from Li(1) fibers, of which 54 were down-regulated and 27 were up-regulated. Several novel aspects regarding cotton fiber elongation can be illustrated from our data. First, over half of the down-regulated proteins were newly identified at the protein level, which is mainly involved in protein folding and stabilization, nucleocytoplasmic transport, signal transduction, and vesicular-mediated transport. Second, a number of cytoskeleton-related proteins showed a remarkable decrease in protein abundance in the Li(1) fibers. Accordingly, the architecture of actin cytoskeleton was severely deformed and the microtubule organization was moderately altered, accompanied with dramatic disruption of vesicle trafficking. Third, the expression of several proteins involved in unfolded protein response (UPR) was activated in Li(1) fibers, indicating that the deficiency of fiber cell elongation was related to ER stress. Collectively, these findings significantly advanced our understanding of the mechanisms associated with cotton fiber elongation.

Breaking the code: Ca2+ sensors in plant signalling

Ca2+ ions play a vital role as second messengers in plant cells during various developmental processes and in response to environmental stimuli. Plants have evolved a diversity of unique proteins that bind Ca2+ using the evolutionarily conserved EF-hand motif. The currently held hypothesis is that these proteins function as Ca2+ sensors by undergoing conformational changes in response to Ca2+-binding that facilitate their regulation of target proteins and thereby co-ordinate various signalling pathways. The three main classes of these EF-hand Ca2+sensors in plants are CaMs [calmodulins; including CMLs (CaM-like proteins)], CDPKs (calcium-dependent protein kinases) and CBLs (calcineurin B-like proteins). In the plant species examined to date, each of these classes is represented by a large family of proteins, most of which have not been characterized biochemically and whose physiological roles remain unclear. In the present review, we discuss recent advances in research on CaMs and CMLs, CDPKs and CBLs, and we attempt to integrate the current knowledge on the different sensor classes into common physiological themes.