Genes expression analyses of sea-island cotton (Gossypium barbadense L.) during fiber development

GhL1L1 regulates the contents of unsaturated fatty acids by activating the expression of GhFAD2 genes in cotton

Edible oils with high unsaturated fatty acids, particularly oleic acid, are beneficial to human health. Cotton is one of the top five oil crops in the world, but the mechanism of high-quality oil synthesis and regulatory networks in cotton are largely unclear. Here, we identified Leafy cotyledon1-like 1 (GhL1L1), a NF-YB subfamily gene that is specifically expressed during somatic embryogenesis and seed maturation in cotton. Overexpression of GhL1L1 regulates the contents of unsaturated fatty acids in cotton, especially in the seeds, which is associated with altered expression of the cotton fatty acid biosynthesis-related genes. GhL1L1 synergistically enhanced the expression of GhFAD2-1A by binding to the G-box in its promoter, leading to an increase in the content of linoleic acid. Furthermore, this activation could be enhanced by GhNF-YC2 and GhNF-YA1 by form a transcriptional complex. Collectively, these results contribute to provide new insights into the molecular mechanism of oil biosynthesis in cotton and can facilitate genetic manipulation of cotton varieties with enhanced oil content.

Domestication over Speciation in Allopolyploid Cotton Species: A Stronger Transcriptomic Pull

Cotton has been domesticated independently four times for its fiber, but the genomic targets of selection during each domestication event are mostly unknown. Comparative analysis of the transcriptome during cotton fiber development in wild and cultivated materials holds promise for revealing how independent domestications led to the superficially similar modern cotton fiber phenotype in upland (G. hirsutum) and Pima (G. barbadense) cotton cultivars. Here we examined the fiber transcriptomes of both wild and domesticated G. hirsutum and G. barbadense to compare the effects of speciation versus domestication, performing differential gene expression analysis and coexpression network analysis at four developmental timepoints (5, 10, 15, or 20 days after flowering) spanning primary and secondary wall synthesis. These analyses revealed extensive differential expression between species, timepoints, domestication states, and particularly the intersection of domestication and species. Differential expression was higher when comparing domesticated accessions of the two species than between the wild, indicating that domestication had a greater impact on the transcriptome than speciation. Network analysis showed significant interspecific differences in coexpression network topology, module membership, and connectivity. Despite these differences, some modules or module functions were subject to parallel domestication in both species. Taken together, these results indicate that independent domestication led G. hirsutum and G. barbadense down unique pathways but that it also leveraged similar modules of coexpression to arrive at similar domesticated phenotypes.

Cloning and functional analysis of GhDFR1, a key gene of flavonoid synthesis pathway in naturally colored cotton

BACKGROUND

The naturally colored brown cotton fiber is the most widely used environmentally friendly textile material, which primarily contains proanthocyanidins and their derivatives. Many structural genes in the flavonoid synthesis pathway are known to improve the genetic resources of naturally colored cotton. Among them, DFR is a crucial late enzyme to synthesis both anthocyanins and proanthocyanidins in the plant flavonoid pathway.

METHODS

The protein sequences of GhDFRs were analyzed using bioinformatic tools. The expression levels of GhDFRs in various tissues and organs of upland cotton Zongxu1 (ZX1), were analyzed by quantitative real-time PCR, and the expression pattern of GhDFR1 during fiber development of white cotton and brown cotton was analyzed further. The function of GhDFR1 in NCC ZX1 was preliminarily analyzed by virus induced gene silencing (VIGS) technology.

RESULTS

Bioinformatic analysis revealed that GhDFRs sequences in upland cotton genome were extremely conserved. Furthermore, evolutionary tree analysis revealed that the functions of GhDFR1 and GhDFR2, and GhDFR3 and GhDFR4, presented different and shared some similarities. Our study showed GhDFR1 and GhDFR2 were specifically expressed in fibers, while GhDFR3 and GhDFR4 were specifically expressed in petals. GhDFR1 was exclusively expressed in brown cotton fiber at various stages of development and progressively increased with the growth of fiber, but the trend of expression in white cotton was quite the opposite. We silenced GhDFR1 expression in brown cotton fiber using VIGS technology, and observed the VIGS-interference plants. After reducing the expression level of GhDFR1, the period for significant GhDFR1 expression in the developing fibers changed, reducing the content of anthocyanins, and lightening the color of mature cotton fibers.

CONCLUSION

GhDFR1 was preferentially expressed in brown cotton during fiber development. The timing of GhDFR1 expression for flavonoid synthesis altered, resulting in anthocyanin contents reduced and the fiber color of the GhDFR1i lines lightened. These findings showed the role of GhDFR1 in fiber coloration of NCC and provided a new candidate for NCC genetic improvement.

Single-cell RNA-seq reveals fate determination control of an individual fibre cell initiation in cotton (Gossypium hirsutum)

Cotton fibre is a unicellular seed trichome, and lint fibre initials per seed as a factor determines fibre yield. However, the mechanisms controlling fibre initiation from ovule epidermis are not understood well enough. Here, with single-cell RNA sequencing (scRNA-seq), a total of 14 535 cells were identified from cotton ovule outer integument of Xu142_LF line at four developmental stages (1.5, 1, 0.5 days before anthesis and the day of anthesis). Three major cell types, fibre, non-fibre epidermis and outer pigment layer were identified and then verified by RNA in situ hybridization. A comparative analysis on scRNA-seq data between Xu142 and its fibreless mutant Xu142 fl further confirmed fibre cluster definition. The developmental trajectory of fibre cell was reconstructed, and fibre cell was identified differentiated at 1 day before anthesis. Gene regulatory networks at four stages revealed the spatiotemporal pattern of core transcription factors, and MYB25-like and HOX3 were demonstrated played key roles as commanders in fibre differentiation and tip-biased diffuse growth respectively. A model for early development of a single fibre cell was proposed here, which sheds light on further deciphering mechanism of plant trichome and the improvement of cotton fibre yield.

Analysis of transcriptome data and quantitative trait loci enables the identification of candidate genes responsible for fiber strength in Gossypium barbadense

Gossypium barbadense possesses a superior fiber quality because of its fiber length and strength. An in-depth analysis of the underlying genetic mechanism could aid in filling the gap in research regarding fiber strength and could provide helpful information for Gossypium barbadense breeding. Three quantitative trait loci related to fiber strength were identified from a Gossypium barbadense recombinant inbred line (PimaS-7 × 5917) for further analysis. RNA sequencing was performed in the fiber tissues of PimaS-7 × 5917 0–35 days postanthesis. Four specific modules closely related to the secondary wall-thickening stage were obtained using the weighted gene coexpression network analysis. In total, 55 genes were identified as differentially expressed from 4 specific modules. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes were used for enrichment analysis, and Gbar_D11G032910, Gbar_D08G020540, Gbar_D08G013370, Gbar_D11G033670, and Gbar_D11G029020 were found to regulate fiber strength by playing a role in the composition of structural constituents of cytoskeleton and microtubules during fiber development. Quantitative real-time PCR results confirmed the accuracy of the transcriptome data. This study provides a quick strategy for exploring candidate genes and provides new insights for improving fiber strength in cotton.

Data mining of transcriptional biomarkers at different cotton fiber developmental stages

Advancement of the gene expression study provides comprehensive information on pivotal genes at different cotton fiber development stages. For the betterment of cotton fiber yield and their quality, genetic improvement is a major target point for the cotton community. Therefore, various studies were carried out to understand the transcriptional machinery of fiber leading to the detailed integrative as well as innovative study. Through data mining and statistical approaches, we identified and validated the transcriptional biomarkers for staged specific differentiation of fiber. With the unique mapping read matrix of ~ 200 cotton transcriptome data and sequential statistical analysis, we identified several important genes that have a deciding and specific role in fiber cell commitment, initiation and elongation, or secondary cell wall synthesis stage. Based on the importance score and validation analysis, IQ domain 26, Aquaporin, Gibberellin regulated protein, methionine gamma lyase, alpha/beta hydrolases, and HAD-like superfamily have shown the specific and determining role for fiber developmental stages. These genes are represented as transcriptional biomarkers that provide a base for molecular characterization for cotton fiber development which will ultimately determine the high yield.

A 13-Lipoxygenase, GhLOX2, positively regulates cotton tolerance against Verticillium dahliae through JA-mediated pathway

Verticillium wilt is a major limiting factor for sustainable production of cotton but the mechanism of controlling this disease is still poorly understood. Lipoxygenase (LOX)-derived oxylipins have been implicated in defense responses against diverse pathogens; however there is limited information about the functional characterization of LOXs in response to Verticillium dahliae infection. In this study, we report the characterization of a cotton LOX gene, GhLOX2, which phylogenetically clustered into 13-LOX subfamily and is closely related to Arabidopsis LOX2 gene. GhLOX2 was predominantly expressed in leaves and strongly induced following V. dahliae inoculation and treatment of methyl jasmonate (MeJA). RNAi-mediated knock-down of GhLOX2 enhanced cotton susceptibility to V. dahliae and was coupled with suppression of jasmonic acid (JA)-related genes both after inoculation with the cotton defoliating strain V991 or MeJA treatment. Interestingly, lignin contents, transcripts of lignin synthesis genes and H₂O₂ contents were also decreased in GhLOX2-silenced plants. This study suggests that GhLOX2 is involved in defense responses against infection of V. dahliae in cotton and supports that JA is one of the major defense hormones against this pathogen.

Overexpression of KcNHX1 gene confers tolerance to multiple abiotic stresses in Arabidopsis thaliana

Abiotic stresses such as drought, salinity, and heat affect plant growth and development. Karelinia caspica is a unique perennial herb that grows in desert area for a long time and has strong tolerance to environmental stresses. In order to explore the functions of the Na⁺/H⁺ antiporter gene from eremophyte K. caspica (KcNHX1) in the abiotic stress response of K. caspica and the underlying regulatory mechanisms, we constructed a vector overexpressing KcNHX1 and transformed it into Arabidopsis thaliana. The physiological results showed that the overexpression of KcNHX1 in A. thaliana not only enhanced the plant's tolerance to salt stress, but also enhanced its tolerance to drought and heat stress at the seedling stage. In addition, KcNHX1-overexpressing plants exhibited enhanced reproductive growth under high temperature, which was mediated by increased auxin accumulation. Taken together, our results indicate that KcNHX1 from an eremophyte can be used as a candidate gene to improve multiple stress tolerance in other plants.

Genome-wide analysis of CBL and CIPK family genes in cotton: conserved structures with divergent interactions and expression

Calcineurin B-like proteins (CBLs) interact with CBL-interacting protein kinases (CIPKs) to form complex molecular modules in response to diverse abiotic stresses. Although previous studies demonstrated that the CBL-CIPK networks play a crucial role in plants response to abiotic stresses, however, little is known about their functions in cotton. In the present study, a total of 22 GhCBL and 79 GhCIPK gene family members were identified in upland cotton (Gossypium hirsutum Linn). Synteny analysis revealed that most genes of GhCBL and GhCIPK exist in pairs between At sub-genome and Dt sub-genome. Interaction analysis between GhCBL and GhCIPK proteins by yeast two-hybrid (Y2H) suggested that the GhCBL-GhCIPK networks were complex, and exhibited functional redundancy in cotton. Quantitative expression analysis by public transcriptome datasets revealed that some GhCBL and GhCIPK genes are differentially expressed under abiotic stress treatments, and especially under drought stress. Our results not only contribute to understanding the structural features of GhCBL and GhCIPK genes but also provide the basis for in-depth functional studies of GhCBL-GhCIPK networks in stress response for plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at (doi:10.1007/s12298-021-00943-1).

Transcriptome analysis reveals differences in the mechanisms of fiber initiation and elongation between long- and short-fiber cotton (Gossypium hirsutum L.) lines

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Crosstalk between cytokinin and ethylene signaling pathways regulates leaf abscission in cotton in response to chemical defoliants

Abscission is a process that allows plants to shed tissues or organs via cell separation, and occurs throughout the life cycle. Removal of leaves through the use of chemical defoliants is very important for mechanical harvesting of cotton (Gossypium hirsutum). However, our knowledge of the molecular mechanisms of the defoliation response involved is limited. In this study, RNA-seq was conducted in order to profile the differentially expressed genes (DEGs) between cultivars X50 (sensitive to chemical defoliants) and X33 (relatively insensitive) at different time points after treatment with thidiazuron and ethephon (TE). A total of 2434 DEGs were identified between the two cultivars across the different time-points. Functional categories according to GO and KEGG analyses revealed that plant hormone signal transduction and zeatin biosynthesis were involved in the response to TE. Cytokinin oxidase/dehydrogenase (CKX) genes and ethylene-related genes were up-regulated following TE treatment, and were associated with increased level of ethylene, especially in cultivar X50. Down-regulation of GhCKX3 resulted in delayed defoliation and a reduced ethylene response. The results show that crosstalk between cytokinin and ethylene regulates cotton defoliation, and provide new insights into the molecular mechanisms underlying the mode of action of defoliants in cotton.

Genome-Wide Identification of Papain-Like Cysteine Proteases in Gossypium hirsutum and Functional Characterization in Response to Verticillium dahliae

Cotton, a natural fiber producing crop of huge importance, is often prone to attack of Verticillium dahliae. Papain-like cysteine proteases (PLCPs) constitute a large family in plants and were proposed to involve in plant defense against pathogen attack in a number of studies. However, there is no detailed characterization of PLCP genes in cotton against infection of V. dahliae. In this study, we carried out a genome-wide analysis in cotton and identified seventy-eight PLCPs, which were divided into nine subfamilies based on their evolution phylogeny: RD21 (responsive to desiccation 21), CEP (cysteine endopeptidase), XCP (xylem cysteine peptidase), XBCP3 (xylem bark cysteine peptidase 3), THI, SAG12 (senescence-associated gene 12), RD19 (responsive to desiccation 19), ALP (aleurain-like protease) and CTB (cathepsin B-like). Genes in each subfamily exhibit a similar structure and motif composition. The expression patterns of these genes in different organs were examined, and subfamily RD21 was the most abundant in these families. Expression profiles under abiotic stress showed that thirty-five PLCP genes were induced by multiple stresses. Further transcriptome analysis showed that sixteen PLCP genes were up-regulated in response to V. dahliae in cotton. Among those, GhRD21-7 showed a higher transcription level than most other PLCP genes. Additionally, over-expression of GhRD21-7 led to enhanced resistance and RNAi lines were more susceptible to V. dahliae in cotton. Our results provide valuable information for future functional genomic studies of PLCP gene family in cotton.

GhL1L1 affects cell fate specification by regulating GhPIN1-mediated auxin distribution

Auxin is as an efficient initiator and regulator of cell fate during somatic embryogenesis (SE), but the molecular mechanisms and regulating networks of this process are not well understood. In this report, we analysed SE process induced by Leafy cotyledon1-like 1 (GhL1L1), a NF-YB subfamily gene specifically expressed in embryonic tissues in cotton. We also identified the target gene of GhL1L1, and its role in auxin distribution and cell fate specification during embryonic development was analysed. Overexpression of GhL1L1 accelerated embryonic cell formation, associated with an increased concentration of IAA in embryogenic calluses (ECs) and in the shoot apical meristem, corresponding to altered expression of the auxin transport gene GhPIN1. By contrast, GhL1L1-deficient explants showed retarded embryonic cell formation, and the concentration of IAA was decreased in GhL1L1-deficient ECs. Disruption of auxin distribution accelerated the specification of embryonic cell fate together with regulation of GhPIN1. Furthermore, we showed that PHOSPHATASE 2AA2 (GhPP2AA2) was activated by GhL1L1 through targeting the G-box of its promoter, hence regulating the activity of GhPIN1 protein. Our results indicate that GhL1L1 functions as a key regulator in auxin distribution to regulate cell fate specification in cotton and contribute to the understanding of the complex process of SE in plant species.

Genome-wide identification of lipoxygenase gene family in cotton and functional characterization in response to abiotic stresses

BACKGROUND

Plant lipoxygenase (LOX) genes are members of the non-haeme iron-containing dioxygenase family that catalyze the oxidation of polyunsaturated fatty acids into functionally diverse oxylipins. The LOX family genes have been extensively studied under biotic and abiotic stresses, both in model and non-model plant species; however, information on their roles in cotton is still limited.

RESULTS

A total of 64 putative LOX genes were identified in four cotton species (Gossypium (G. hirsutum, G. barbadense, G. arboreum, and G. raimondii)). In the phylogenetic tree, these genes were clustered into three categories (9-LOX, 13-LOX type I, and 13-LOX type II). Segmental duplication of putative LOX genes was observed between homologues from A2 to At and D5 to Dt hinting at allopolyploidy in cultivated tetraploid species (G. hirsutum and G. barbadense). The structure and motif composition of GhLOX genes appears to be relatively conserved among the subfamilies. Moreover, many cis-acting elements related to growth, stresses, and phytohormone signaling were found in the promoter regions of GhLOX genes. Gene expression analysis revealed that all GhLOX genes were induced in at least two tissues and the majority of GhLOX genes were up-regulated in response to heat and salinity stress. Specific expressions of some genes in response to exogenous phytohormones suggest their potential roles in regulating growth and stress responses. In addition, functional characterization of two candidate genes (GhLOX12 and GhLOX13) using virus induced gene silencing (VIGS) approach revealed their increased sensitivity to salinity stress in target gene-silenced cotton. Compared with controls, target gene-silenced plants showed significantly higher chlorophyll degradation, higher H2O2, malondialdehyde (MDA) and proline accumulation but significantly reduced superoxide dismutase (SOD) activity, suggesting their reduced ability to effectively degrade accumulated reactive oxygen species (ROS).

CONCLUSION

This genome-wide study provides a systematic analysis of the cotton LOX gene family using bioinformatics tools. Differential expression patterns of cotton LOX genes in different tissues and under various abiotic stress conditions provide insights towards understanding the potential functions of candidate genes. Beyond the findings reported here, our study provides a basis for further uncovering the biological roles of LOX genes in cotton development and adaptation to stress conditions.

Transcriptome, cytological and biochemical analysis of cytoplasmic male sterility and maintainer line in CMS-D8 cotton

Key message This research based on RNA-seq, biochemical, and cytological analyses sheds that ROS may serve as important signaling molecules of cytoplasmic male sterility in CMS-D8 cotton. To understand the mechanism of cytoplasmic male sterility in cotton (Gossypium hirsutum), transcriptomic, cytological, and biochemical analysis were performed between the cytoplasmic male sterility CMS-D8 line, Zhong41A, and its maintainer line Zhong41B. A total of 2335 differentially expressed genes (DEGs) were identified in the CMS line at three different stages of anther development. Bioinformatics analysis of these DEGs indicated their relationship to reactive oxygen species (ROS) homeostasis, including reduction-oxidation reactions and the metabolism of glutathione and ascorbate. At the same time, DEGs associated with tapetum development, especially the transition to secretory tapetum, were down-regulated in the CMS line. Biochemical analysis indicated that the ability of the CMS line to eliminate ROS was decreased, which led to the rapid release of H₂O₂. Cytological analysis revealed that the most crucial defect in the CMS line was the abnormal tapetum. All these results are consistent with the RNA sequencing data. On the basis of our findings, we propose that ROS act as signal molecules, which are released from mitochondria and transferred to the nucleus, triggering the formation of abnormal tapetum.

Transcriptomic repertoires depict the initiation of lint and fuzz fibres in cotton (Gossypium hirsutum L.)

Cotton fibre is an important natural fibre for the textile industry. The number of fibre initials determines the lint percentage, which is an important factor for cotton fibre yield. Although fibre development has been described by transcriptomic analysis, the mechanism by which the long noncoding RNA manipulates the initiation of lint and fuzz fibres remains unknown. In this study, three lines with different lint percentages were developed by crossing Xu142 with its fibreless mutant Xu142 fl. We collected the epidermal cells from the ovules with attached fibres at 0 and 5 days post anthesis (DPA) from Xu142, the fibreless mutant Xu142 fl and the three lint percent diversified lines for deep transcriptome sequencing. A total of 2641 novel genes, 35 802 long noncoding RNAs (lncRNAs) and 2262 circular RNAs (circRNAs) were identified, of which 645 lncRNAs were preferentially expressed in the fibreless mutant Xu142 fl and 651 lncRNAs were preferentially expressed in the fibre-attached lines. We demonstrated the functional roles of the three lncRNAs in fibre development via a virus-induced gene silencing (VIGS) system. Our results showed that silencing XLOC_545639 and XLOC_039050 in Xu142 fl increased the number of fibre initials on the ovules, but silencing XLOC_079089 in Xu142 resulted in a short fibre phenotype. This study established the transcriptomic repertoires in cotton fibre initiation and provided evidence for the potential functions of lncRNAs in fibre development.

Overexpression of rice thaumatin-like protein (Ostlp) gene in transgenic cassava results in enhanced tolerance to Colletotrichum gloeosporioides f. sp. manihotis

Cassava (Manihot esculenta Crantz) is the most important staple food for more than 300 million people in Africa, and anthracnose disease caused by Colletotrichum gloeosporioides f. sp. manihotis is the most destructive fungal disease affecting cassava production in sub-Saharan Africa. The main objective of this study was to improve anthracnose resistance in cassava through genetic engineering. Transgenic cassava plants harbouring rice thaumatin-like protein (Ostlp) gene, driven by the constitutive CaMV35S promoter, were generated using Agrobacterium-mediated transformation of friable embryogenic calli (FEC) of cultivar TMS 60444. Molecular analysis confirmed the presence, integration, copy number of the transgene all the independent transgenic events. Semi-quantitative RT-PCR confirmed high expression levels of Ostlp in six transgenic lines tested. The antifungal activity of the transgene against Colletotrichum gloeosporioides pathogen was evaluated using the leaves and stem cuttings bioassay. The results demonstrated significantly delayed disease development and reduced size of necrotic lesions in leaves and stem cuttings of all transgenic lines compared to the leaves and stem cuttingss of non-transgenic control plants. Therefore, constitutive overexpression of rice thaumatin-like protein in transgenic cassava confers enhanced tolerance to the fungal pathogen C. gloeosporioides f. sp. manihotis. These results can therefore serve as an initial step towards genetic engineering of farmer-preffered cassava cultivars for resistance to anthracnose disease.

The JASMONATE ZIM-Domain Gene Family Mediates JA Signaling and Stress Response in Cotton

JASMONATE ZIM-domain (JAZ) family proteins are involved in regulating diverse biological processes in plants. However, their functions have not been well characterized in cotton (Gossypium spp.). In the present study, 13, 15, 25 and 30 JAZ genes were identified in Gossypium arboretum, Gossypium raimondii, Gossypium barbadense and Gossypium hirsutum, respectively, based on gene homology. Selection and variation analyses showed that the single nucleotide polymorphism (SNP) density of GhJAZ genes in wild species was much higher than that in domesticated species. Expression pattern analysis showed that all the GhJAZ genes are expressed in at least one tissue and respond to one or more stress factors, as well as being induced by some phytohormones. The co-expression network indicated that GhJAZ genes might mediate multiple stress response pathways. Yeast two-hybrid (Y2H) experiments showed extensive interactions among GhJAZ proteins, including homo- and heteromeric interactions. Overexpressing one member of the JAZ gene family, GhJAZ2 (Gh_D06G0810), significantly enhanced sensitivity to salt stress in transgenic cotton. Transcriptome analysis indicated that GhJAZ2 regulates stress responses possibly by participating in α-linolenic acid metabolism and jasmonate signaling, and is involved in the repression of GhMYC2 regulated by GhJAZ2. Our data provide important clues for further elucidating the functions of JAZ genes in cotton.

Transcriptome Analysis Suggests That Chromosome Introgression Fragments from Sea Island Cotton (Gossypium barbadense) Increase Fiber Strength in Upland Cotton (Gossypium hirsutum)

As high-strength cotton fibers are critical components of high quality cotton, developing cotton cultivars with high-strength fibers as well as high yield is a top priority for cotton development. Recently, chromosome segment substitution lines (CSSLs) have been developed from high-yield Upland cotton (Gossypium hirsutum) crossed with high-quality Sea Island cotton (G. barbadense). Here, we constructed a CSSL population by crossing CCRI45, a high-yield Upland cotton cultivar, with Hai1, a Sea Island cotton cultivar with superior fiber quality. We then selected two CSSLs with significantly higher fiber strength than CCRI45 (MBI7747 and MBI7561), and one CSSL with lower fiber strength than CCRI45 (MBI7285), for further analysis. We sequenced all four transcriptomes at four different time points postanthesis, and clustered the 44,678 identified genes by function. We identified 2200 common differentially-expressed genes (DEGs): those that were found in both high quality CSSLs (MBI7747 and MBI7561), but not in the low quality CSSL (MBI7285). Many of these genes were associated with various metabolic pathways that affect fiber strength. Upregulated DEGs were associated with polysaccharide metabolic regulation, single-organism localization, cell wall organization, and biogenesis, while the downregulated DEGs were associated with microtubule regulation, the cellular response to stress, and the cell cycle. Further analyses indicated that three genes, XLOC_036333 [mannosyl-oligosaccharide-α-mannosidase (MNS1)], XLOC_029945 (FLA8), and XLOC_075372 (snakin-1), were potentially important for the regulation of cotton fiber strength. Our results suggest that these genes may be good candidates for future investigation of the molecular mechanisms of fiber strength formation and for the improvement of cotton fiber quality through molecular breeding.

Comparative transcriptome analysis of cotton fiber development of Upland cotton (Gossypium hirsutum) and Chromosome Segment Substitution Lines from G. hirsutum × G. barbadense

Background

How to develop new cotton varieties possessing high yield traits of Upland cotton and superior fiber quality traits of Sea Island cotton remains a key task for cotton breeders and researchers. While multiple attempts bring in little significant progresses, the development of Chromosome Segment Substitution Lines (CSSLs) from Gossypium barbadense in G. hirsutum background provided ideal materials for aforementioned breeding purposes in upland cotton improvement. Based on the excellent fiber performance and relatively clear chromosome substitution segments information identified by Simple Sequence Repeat (SSR) markers, two CSSLs, MBI9915 and MBI9749, together with the recurrent parent CCRI36 were chosen to conduct transcriptome sequencing during the development stages of fiber elongation and Secondary Cell Wall (SCW) synthesis (from 10DPA and 28DPA), aiming at revealing the mechanism of fiber development and the potential contribution of chromosome substitution segments from Sea Island cotton to fiber development of Upland cotton.

Results

In total, 15 RNA-seq libraries were constructed and sequenced separately, generating 705.433 million clean reads with mean GC content of 45.13% and average Q30 of 90.26%. Through multiple comparisons between libraries, 1801 differentially expressed genes (DEGs) were identified, of which the 902 up-regulated DEGs were mainly involved in cell wall organization and response to oxidative stress and auxin, while the 898 down-regulated ones participated in translation, regulation of transcription, DNA-templated and cytoplasmic translation based on GO annotation and KEGG enrichment analysis. Subsequently, STEM software was performed to explicate the temporal expression pattern of DEGs. Two peroxidases and four flavonoid pathway-related genes were identified in the “oxidation-reduction process”, which could play a role in fiber development and quality formation. Finally, the reliability of RNA-seq data was validated by quantitative real-time PCR of randomly selected 20 genes.

Conclusions

The present report focuses on the similarities and differences of transcriptome profiles between the two CSSLs and the recurrent parent CCRI36 and provides novel insights into the molecular mechanism of fiber development, and into further exploration of the feasible contribution of G. barbadense substitution segments to fiber quality formation, which will lay solid foundation for simultaneously improving fiber yield and quality of upland cotton through CSSLs.

Electronic supplementary material

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

Identification of candidate genes for fiber length quantitative trait loci through RNA-Seq and linkage and physical mapping in cotton

BACKGROUND

Cotton (Gossypium spp.) fibers are single-celled elongated trichomes, the molecular aspects of genetic variation in fiber length (FL) among genotypes are currently unknown. In this study, two backcross inbred lines (BILs), i.e., NMGA-062 ("Long") and NMGA-105 ("Short") with 32.1 vs. 27.2 mm in FL, respectively, were chosen to perform RNA-Seq on developing fibers at 10 days post anthesis (DPA). The two BILs differed in 4 quantitative trait loci (QTL) for FL and were developed from backcrosses between G. hirsutum as the recurrent parent and G. barbadense.

RESULTS

In total, 51.7 and 54.3 million reads were obtained and assembled to 49,508 and 49,448 transcripts in the two genotypes, respectively. Of 1551 differentially expressed genes (DEGs) between the two BILs, 678 were up-regulated and 873 down-regulated in "Long"; and 703 SNPs were identified in 339 DEGs. Further physical mapping showed that 8 DEGs were co-localized with the 4 FL QTL identified in the BIL population containing the two BILs. Four SNP markers in 3 DEGs that showed significant correlations with FL were developed. Among the three candidate genes encoding for proline-rich protein, D-cysteine desulfhydrase, and thaumatin-like protein, a SNP of thaumatin-like protein gene showed consistent correlations with FL across all testing environments.

CONCLUSIONS

This study represents one of the first investigations of positional candidate gene approach of QTL in cotton in integrating transcriptome and SNP identification based on RNA-Seq with linkage and physical mapping of QTL and genes, which will facilitate eventual cloning and identification of genes responsible for FL QTL. The candidate genes may serve as the foundation for further in-depth studies of the molecular mechanism of natural variation in fiber elongation.

A high-efficiency CRISPR/Cas9 system for targeted mutagenesis in Cotton (Gossypium hirsutum L.)

The complex allotetraploid genome is one of major challenges in cotton for repressing gene expression. Developing site-specific DNA mutation is the long-term dream for cotton breeding scientists. The clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system is emerging as a robust biotechnology for targeted-DNA mutation. In this study, two sgRNAs, GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2, were designed in the identical genomic regions of GhMYB25-like A and GhMYB25-like D, which were encoded by cotton A subgenome and the D subgenome, respectively, was assembled to direct Cas9-mediated allotetraploid cotton genome editing. High proportion (14.2–21.4%) CRISPR/Cas9-induced specific truncation events, either from GhMYB25-like A DNA site or from GhMYB25-like D DNA site, were detected in 50% examined transgenic cotton through PCR amplification assay and sequencing analyses. Sequencing results also demonstrated that 100% and 98.8% mutation frequency were occurred on GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2 target site respectively. The off-target effect was evaluated by sequencing two putative off-target sites, which have 3 and 1 mismatched nucleotides with GhMYB25-like-sgRNA1 and GhMYB25-like-sgRNA2, respectively; all the examined samples were not detected any off-target-caused mutation events. Thus, these results demonstrated that CRISPR/Cas9 is qualified for generating DNA level mutations on allotetraploid cotton genome with high-efficiency and high-specificity.

Identification of QTL for Fiber Quality and Yield Traits Using Two Immortalized Backcross Populations in Upland Cotton

Two immortalized backcross populations (DHBCF₁s and JMBCF₁s) were developed using a recombinant inbred line (RIL) population crossed with the two parents DH962 and Jimian5 (as the males), respectively. The fiber quality and yield component traits of the two backcross populations were phenotyped at four environments (two locations, two years). One hundred seventy-eight quantitative trait loci (QTL) were detected including 76 for fiber qualities and 102 for yield components, explaining 4.08–17.79% of the phenotypic variation (PV). Among the 178 QTL, 22 stable QTL were detected in more than one environment or population. A stable QTL, qFL-c10-1, was detected in the previous F₂ population, a RIL population in 3 environments and the current two BCF₁ populations in this study, explaining 5.79–37.09% of the PV. Additionally, 117 and 110 main-effect QTL (M-QTL) and 47 and 191 digenic epistatic QTL (E-QTL) were detected in the DHBCF₁s and JMBCF₁s populations, respectively. The effect of digenic epistasis played a more important role on lint percentage, fiber length and fiber strength. These results obtained in the present study provided more resources to obtain stable QTL, confirming the authenticity and reliability of the QTL for molecular marker-assisted selection breeding and QTL cloning.

GhJAZ2 negatively regulates cotton fiber initiation by interacting with the R2R3-MYB transcription factor GhMYB25-like

Jasmonic acid (JA) signaling has been well studied in Arabidopsis. Most reports focus on the role of JA in biological pathways, such as stress resistance, trichome initiation and anthocyanin accumulation. The JASMONATE ZIM-DOMAIN (JAZ) protein is one of the important repressors in the JA signaling pathway. Previous studies showed that JA functions in fiber initiation and elongation, but little is known about how JAZ genes function in fiber development. In this study, a cotton JAZ protein (GhJAZ2) containing a highly conserved TIFY motif and a C-terminal Jas domain was identified, and its function during cotton fiber development was analysed. Gene expression analysis showed that GhJAZ2 was preferentially expressed in the root, hypocotyl, flower and ovule 1 day before anthesis. Overexpression of GhJAZ2 inhibited both lint and fuzz fiber initiation, and reduced the fiber length. Yeast two-hybrid assays showed that GhJAZ2 interacted with the R2R3-MYB transcription factors GhMYB25-like and GhGL1, the bHLH transcription factor GhMYC2, the WD repeat protein GhWD40 and the unknown protein GhJI1. Among these transcription factors, previous studies showed that downregulation of GhMYB25-like leads to a fiberless phenotype in cotton seeds. Molecular and genetic evidence showed that the GhJAZ2 protein suppresses fiber initiation in the overexpressing lines by interacting with GhMYB25-like and suppressing GhMYB25-like activity. Our results suggested that GhJAZ2 functions as a primary transcription repressor during lint and fuzz fiber initiation by interacting with GhMYB25-like, GhGL1, GhMYC2, GhWD40 and GhJI1 to regulate the JA signaling pathway.

Fibre elongation requires normal redox homeostasis modulated by cytosolic ascorbate peroxidase in cotton (Gossypium hirsutum)

High-quality cotton fibre equates to a more comfortable textile. Fibre length is an important index of fibre quality. Hydrogen peroxide (H2O2) acts as a signalling molecule in the regulation of fibre elongation. Results from in vitro ovule culture suggest that the alteration of fibre cell H2O2 levels affects fibre development. Ascorbate peroxidase (APX) is an important reactive oxygen species (ROS) scavenging enzyme, and we found that GhAPX1AT/DT encoded one member of the previously unrealized group of cytosolic APXs (cAPXs) that were preferentially expressed during the fibre elongation stage. Transgenic cottons with up- and down-regulation of GhAPX1AT/DT were generated to control fibre endogenous levels of H2O2 Suppression of all cAPX (IAO) resulted in a 3.5-fold increase in H2O2 level in fibres and oxidative stress, which significantly suppressed fibre elongation. The fibre length of transgenic lines with over-expression or specific down-regulation of GhAPX1AT/DT did not show any obvious change. However, the fibres in the over-expression lines exhibited higher tolerance to oxidative stress. Differentially expressed genes (DEGs) in fibres at 10 days post-anthesis (DPA) of IAO lines identified by RNA-seq were related to redox homeostasis, signalling pathways, stress responses and cell wall synthesis, and the DEGs that were up-regulated in IAO lines were also up-regulated in the 10 DPA and 20 DPA fibres of wild cotton compared with domesticated cotton. These results suggest that optimal H2O2 levels and redox state regulated by cytosolic APX are key mechanisms regulating fibre elongation, and dysregulation of the increase in H2O2 induces oxidative stress and results in shorter fibres by initiating secondary cell wall-related gene expression.

Identification of Multiple Stress Responsive Genes by Sequencing a Normalized cDNA Library from Sea-Land Cotton (Gossypium barbadense L.)

Background

Plants often face multiple stresses including drought, extreme temperature, salinity, nutrition deficiency and biotic stresses during growth and development. All the stresses result in a series of physiological and metabolic reactions and then generate reversible inhibition of metabolism and growth and can cause seriously irreversible damage, even death. At each stage of cotton growth, environmental stress conditions pose devastating threats to plant growth and development, especially yield and quality. Due to the complex stress conditions and unclear molecular mechanisms of stress response, there is an urgent need to explore the mechanisms of cotton response against abiotic stresses.

Methodology and Principal Findings

A normalized cDNA library was constructed using Gossypium barbadense Hai-7124 treated with different stress conditions (heat, cold, salt, drought, potassium and phosphorus deficit and Verticillium dahliae infection). Random sequencing of this library generated 6,047 high-quality expressed sequence tags (ESTs). The ESTs were clustered and assembled into 3,135 uniESTs, composed of 2,497 contigs and 638 singletons. The blastx results demonstrated 2,746 unigenes showing significant similarity to known genes, 74 uniESTs displaying significant similarity to genes of predicted proteins, and 315 uniESTs remain uncharacterized. Functional classification unveiled the abundance of uniESTs in binding, catalytic activity, and structural molecule activity. Annotations of the uniESTs by the plant transcription factor database (PlantTFDB) and Plant Stress Protein Database (PSPDB) disclosed that transcription factors and stress-related genes were enriched in the current library. The expression of some transcription factors and specific stress-related genes were verified by RT-PCR under various stress conditions.

Conclusions/Significance

Annotation results showed that a huge number of genes respond to stress in our study, such as MYB-related, C2H2, FAR1, bHLH, bZIP, MADS, and mTERF. These results will improve our knowledge of stress tolerance in cotton. In addition, they are also helpful in discovering candidate genes related to stress tolerance. The publicly available ESTs from G. barbadense are a valuable genomic resource that will facilitate further molecular study and breeding of stress-tolerant cotton.

GbEXPATR, a species-specific expansin, enhances cotton fibre elongation through cell wall restructuring

Cotton provides us the most important natural fibre. High fibre quality is the major goal of cotton breeding, and introducing genes conferring longer, finer and stronger fibre from Gossypium barbadense to Gossypium hirsutum is an important breeding strategy. We previously analysed the G. barbadense fibre development mechanism by gene expression profiling and found two homoeologous fibre-specific α-expansins from G. barbadense, GbEXPA2 and GbEXPATR. GbEXPA2 (from the DT genome) is a classical α-expansin, while its homoeolog, GbEXPATR (AT genome), encodes a truncated protein lacking the normal C-terminal polysaccharide-binding domain of other α-expansins and is specifically expressed in G. barbadense. Silencing EXPA in G. hirsutum induced shorter fibres with thicker cell walls. GbEXPA2 overexpression in G. hirsutum had no effect on mature fibre length, but produced fibres with a slightly thicker wall and increased crystalline cellulose content. Interestingly, GbEXPATR overexpression resulted in longer, finer and stronger fibres coupled with significantly thinner cell walls. The longer and thinner fibre was associated with lower expression of a number of secondary wall-associated genes, especially chitinase-like genes, and walls with lower cellulose levels but higher noncellulosic polysaccharides which advocated that a delay in the transition to secondary wall synthesis might be responsible for better fibre. In conclusion, we propose that α-expansins play a critical role in fibre development by loosening the cell wall; furthermore, a truncated form, GbEXPATR, has a more dramatic effect through reorganizing secondary wall synthesis and metabolism and should be a candidate gene for developing G. hirsutum cultivars with superior fibre quality.

A cotton fiber-preferential promoter, PGbEXPA2, is regulated by GA and ABA in Arabidopsis

PGbEXPA2 (Promoter of GbEXPA2 ) was preferentially and strongly expressed during cotton fiber development, and the 461-bp PGbEXPA2 fragment was essential for responding to exogenous GA and ABA in Arabidopsis. Cotton fibers are highly elongated single-cell, unbranched and non-glandular seed trichomes. Previous studies have reported that the transcript level of GbEXPA2 is significantly up-regulated during fiber cell elongation, suggesting that GbEXPA2 has an important function in fiber development. In this study, the promoter of GbEXPA2 (839 bp) from the D(T) sub-genome was isolated from Gossypium barbadense 3-79. Consistent with the expression pattern of GbEXPA2, the promoter PGbEXPA2 was able to express GUS to high levels in elongating fibers, but not in the root, stem, or leaf. In Arabidopsis, GUS activity was only found in the rosette leaf trichomes and rosette leaf vascular tissue, indicating that the transcription factors which bind to PGbEXPA2 in the leaf trichomes of transgenic Arabidopsis were similar to those found in cotton fiber. A deletion analysis of PGbEXPA2 revealed that a 461-bp fragment was sufficient to drive GUS expression in cotton fibers and Arabidopsis rosette leaf trichomes. Exogenous phytohormonal treatments on transgenic Arabidopsis with different promoter lengths (P-839, P-705, P-588 and P-461) showed that GUS activity in Arabidopsis trichomes could be strongly up-regulated by GA and, in contrast, down-regulated by ABA.

Defective pollen wall contributes to male sterility in the male sterile line 1355A of cotton

To understand the mechanisms of male sterility in cotton (Gossypium spp.), combined histological, biochemical and transcription analysis using RNA-Seq was carried out in the anther of the single-gene recessive genic male sterility system of male sterile line 1355A and male fertile line 1355B, which are near-isogenic lines (NILs) differing only in the fertility trait. A total of 2,446 differentially expressed genes were identified between the anthers of 1355AB lines, at three different stages of development. Cluster analysis and functional assignment of differentially expressed genes revealed differences in transcription associated with pollen wall and anther development, including the metabolism of fatty acids, glucose, pectin and cellulose. Histological and biochemical analysis revealed that a major cellular defect in the 1355A was a thicker nexine, consistent with the RNA-seq data, and further gene expression studies implicated differences in fatty acids synthesis and metabolism. This study provides insight into the phenotypic characteristics and gene regulatory network of the genic male sterile line 1355A in upland cotton.

Cotton WRKY1 mediates the plant defense-to-development transition during infection of cotton by Verticillium dahliae by activating JASMONATE ZIM-DOMAIN1 expression

Plants have evolved an elaborate signaling network to ensure an appropriate level of immune response to meet the differing demands of developmental processes. Previous research has demonstrated that DELLA proteins physically interact with JASMONATE ZIM-DOMAIN1 (JAZ1) and dynamically regulate the interaction of the gibberellin (GA) and jasmonate (JA) signaling pathways. However, whether and how the JAZ1-DELLA regulatory node is regulated at the transcriptional level in plants under normal growth conditions or during pathogen infection is not known. Here, we demonstrate multiple functions of cotton (Gossypium barbadense) GbWRKY1 in the plant defense response and during development. Although GbWRKY1 expression is induced rapidly by methyl jasmonate and infection by Verticillium dahliae, our results show that GbWRKY1 is a negative regulator of the JA-mediated defense response and plant resistance to the pathogens Botrytis cinerea and V. dahliae. Under normal growth conditions, GbWRKY1-overexpressing lines displayed GA-associated phenotypes, including organ elongation and early flowering, coupled with the down-regulation of the putative targets of DELLA. We show that the GA-related phenotypes of GbWRKY1-overexpressing plants depend on the constitutive expression of Gossypium hirsutum GhJAZ1. We also show that GhJAZ1 can be transactivated by GbWRKY1 through TGAC core sequences, and the adjacent sequences of this binding site are essential for binding specificity and affinity to GbWRKY1, as revealed by dual-luciferase reporter assays and electrophoretic mobility shift assays. In summary, our data suggest that GbWRKY1 is a critical regulator mediating the plant defense-to-development transition during V. dahliae infection by activating JAZ1 expression.

Comparison of the transcriptome between two cotton lines of different fiber color and quality

To understand the mechanism of fiber development and pigmentation formation, the mRNAs of two cotton lines were sequenced: line Z128 (light brown fiber) was a selected mutant from line Z263 (dark brown fiber). The primary walls of the fiber cell in both Z263 and Z128 contain pigments; more pigments were laid in the lumen of the fiber cell in Z263 compared with that in Z128. However, Z263 contained less cellulose than Z128. A total of 71,895 unigenes were generated: 13,278 (20.26%) unigenes were defined as differentially expressed genes (DEGs) by comparing the library of Z128 with that of Z263; 5,345 (8.16%) unigenes were up-regulated and 7,933 (12.10%) unigenes were down-regulated. qRT-PCR and comparative transcriptional analysis demonstrated that the pigmentation formation in brown cotton fiber was possibly the consequence of an interaction between oxidized tannins and glycosylated anthocyanins. Furthermore, our results showed the pigmentation related genes not only regulated the fiber color but also influenced the fiber quality at the fiber elongation stage (10 DPA). The highly expressed flavonoid gene in the fiber elongation stage could be related to the fiber quality. DEGs analyses also revealed that transcript levels of some fiber development genes (Ca²⁺/CaM, reactive oxygen, ethylene and sucrose phosphate synthase) varied dramatically between these two cotton lines.

Cotton fiber elongation network revealed by expression profiling of longer fiber lines introgressed with different Gossypium barbadense chromosome segments

Background

Cotton fiber, a highly elongated, thickened single cell of the seed epidermis, is a powerful cell wall research model. Fiber length, largely determined during the elongation stage, is a key property of fiber quality. Several studies using expressed sequence tags and microarray analysis have identified transcripts that accumulate preferentially during fiber elongation. To further show the mechanism of fiber elongation, we used Digital Gene Expression Tag Profiling to compare transcriptome data from longer fiber chromosome introgressed lines (CSILs) containing segments of various Gossypium barbadense chromosomes with data from its recurrent parent TM-1 during fiber elongation (from 5 DPA to 20 DPA).

Results

A large number of differentially expressed genes (DEGs) involved in carbohydrate, fatty acid and secondary metabolism, particularly cell wall biosynthesis, were highly upregulated during the fiber elongation stage, as determined by functional enrichment and pathway analysis. Furthermore, DEGs related to hormone responses and transcription factors showed upregulated expression levels in the CSILs. Moreover, metabolic and regulatory network analysis indicated that the same pathways were differentially altered, and distinct pathways exhibited altered gene expression, in the CSILs. Interestingly, mining of upregulated DEGs in the introgressed segments of these CSILs based on D-genome sequence data showed that these lines were enriched in glucuronosyltransferase, inositol-1, 4, 5-trisphosphate 3-kinase and desulfoglucosinolate sulfotransferase activity. These results were similar to the results of transcriptome analysis.

Conclusions

This report provides an integrative network about the molecular mechanisms controlling fiber length, which are mainly tied to carbohydrate metabolism, cell wall biosynthesis, fatty acid metabolism, secondary metabolism, hormone responses and Transcription factors. The results of this study provide new insights into the critical factors associated with cell elongation and will facilitate further research aimed at understanding the mechanisms underlying cotton fiber elongation.

Electronic supplementary material

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

GhH2A12, a replication-dependent histone H2A gene from Gossypium hirsutum, is negatively involved in the development of cotton fiber cells

GhH2A12 was preferentially expressed at the initiation and early elongation stage of cotton fiber development, and overexpression of GhH2A12 caused retardation of fiber initiation and produced shorter fibers. Histone H2A is a component of eukaryotic chromatin whose function has not been studied in cotton. We have isolated an H2A gene encoding 156 amino acids, named GhH2A12. Like other plant histone H2As, GhH2A12 contains a typical SPKK motif in the carboxy-terminal and a plant-unique peptide-binding A/T-rich DNA region, and it was localized to the nucleus. GhH2A12 was preferentially expressed at the initiation and early elongation stage of cotton fiber, from 0 to 5 days post anthesis and the transcript level declined rapidly when the fiber entered the fast elongation stage, suggesting that GhH2A12 was involved in fiber differentiation. Therefore, GhH2A12 overexpression and RNAi transgenic cotton lines were developed via Agrobacterium tumefaciens-mediated transformation. Overexpression of GhH2A12 caused retardation of fiber initiation and produced shorter fibers and lower lint percentages. Moreover, the overexpressors showed negative effects on seedling growth, and the leaf emergence was delayed compared to wild type. However, no significant change in the GhH2A12 suppression line was observed. Coupled with retardation of fiber initiation, upregulation of GhH2A12 downregulated the expression of genes involved in cell-cycle performance. These results suggest that GhH2A12 might regulate fiber differentiation via regulating the cell cycle-related genes.

A peptide hormone gene, GhPSK promotes fibre elongation and contributes to longer and finer cotton fibre

Cotton fibres, the single-celled trichomes derived from the ovule epidermis, provide the most important natural material for the global textile industry. A number of studies have demonstrated that regulating endogenous hormone levels through transgenic approaches can improve cotton fibre qualities. Phytosulfokine-α (PSK-α) is a novel peptide hormone in plants that is involved in regulating cell proliferation and elongation. However, its potential applications in crop genetic improvement have not been evaluated. In this study, we describe how exogenous PSK-α application promotes cotton fibre cell elongation in vitro. Chlorate, an effective inhibitor of peptide sulfation, suppressed fibre elongation in ovule culture. Exogenously applied PSK-α partly restored the chlorate-induced suppression. A putative PSK gene (GhPSK) was cloned from Gossypium hirsutum. Expression pattern analysis revealed that GhPSK is preferentially expressed in rapidly elongating fibre cells (5-20 days postanthesis). Overexpression of GhPSK in cotton increased the endogenous PSK-α level and promoted cotton fibre cell elongation, resulting in longer and finer fibres. Further results from electrophysiological and physiological analyses suggest that GhPSK affects fibre development through regulation of K(+) efflux. Digital gene expression (DGE) profile analysis of GhPSK overexpression lines indicates that PSK signalling may regulate the respiratory electron-transport chain and reactive oxygen species to affect cotton fibre development. These results imply that peptide hormones are involved in cotton fibre growth and suggest a new strategy for the biotechnological improvement of cotton fibre quality.

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.

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.

Comparative transcriptome analysis between somatic embryos (SEs) and zygotic embryos in cotton: evidence for stress response functions in SE development

As a product of asexual reproduction in plants, the somatic embryo (SE) differentiates into a new plantlet via a zygotic embryogenesis-like process. Here, we present the phenotypic and cellular differences between SEs and zygotic embryos (ZEs) revealed by histological section scanning using three parallel development stages of the two types of embryos of cotton (Gossypium hirsutum cv. YZ1), including globular, torpedo and cotyledonary-stages. To identify the molecular characteristics of SE development in cotton, the digital gene expression system was used to profile the genes active during SE and ZE development. A total of 4242 differentially expressed genes (DEGs) were identified in at least one developmental stage. Expression pattern and functional classification analysis based on these DEGs reveals that SE development exhibits a transcriptional activation of stress responses. RT-PCR analysis further confirmed enhanced expression levels of stress-related genes in SEs than in ZEs. Experimental stress treatment, induced by NaCl and ABA, accelerated SE development and increased the transcription of genes related to stress response, in parallel with decelerated proliferation of embryogenic calluses under stress treatment. Our data reveal that SE development involves the activation of stress responses, which we suggest may regulate the balance between cell proliferation and differentiation. These results provide new insight into the molecular mechanisms of SE development and suggest strategies that can be used for regulating the developmental processes of somatic embryogenesis.

The flavonoid pathway regulates the petal colors of cotton flower

Although biochemists and geneticists have studied the cotton flower for more than one century, little is known about the molecular mechanisms underlying the dramatic color change that occurs during its short developmental life following blooming. Through the analysis of world cotton germplasms, we found that all of the flowers underwent color changes post-anthesis, but there is a diverse array of petal colors among cotton species, with cream, yellow and red colors dominating the color scheme. Genetic and biochemical analyses indicated that both the original cream and red colors and the color changes post-anthesis were related to flavonoid content. The anthocyanin content and the expression of biosynthesis genes were both increased from blooming to one day post-anthesis (DPA) when the flower was withering and undergoing abscission. Our results indicated that the color changes and flavonoid biosynthesis of cotton flowers were precisely controlled and genetically regulated. In addition, flavonol synthase (FLS) genes involved in flavonol biosynthesis showed specific expression at 11 am when the flowers were fully opened. The anthocyanidin reductase (ANR) genes, which are responsible for proanthocyanidins biosynthesis, showed the highest expression at 6 pm on 0 DPA, when the flowers were withered. Light showed primary, moderate and little effects on flavonol, anthocyanin and proanthocyanidin biosynthesis, respectively. Flavonol biosynthesis was in response to light exposure, while anthocyanin biosynthesis was involved in flower color changes. Further expression analysis of flavonoid genes in flowers of wild type and a flavanone 3-hydroxylase (F3H) silenced line showed that the development of cotton flower color was controlled by a complex interaction between genes and light. These results present novel information regarding flavonoids metabolism and flower development.

GhCAX3 gene, a novel Ca(2+)/H(+) exchanger from cotton, confers regulation of cold response and ABA induced signal transduction

As a second messenger, Ca(2+) plays a major role in cold induced transduction via stimulus-specific increases in [Ca(2+)]cyt, which is called calcium signature. During this process, CAXs (Ca(2+)/H(+) exchangers) play critical role. For the first time, a putative Ca(2+)/H(+) exchanger GhCAX3 gene from upland cotton (Gossypium hirsutum cv. 'YZ-1') was isolated and characterized. It was highly expressed in all tissues of cotton except roots and fibers. This gene may act as a regulator in cotton's response to abiotic stresses as it could be up-regulated by Ca(2+), NaCl, ABA and cold stress. Similar to other CAXs, it was proved that GhCAX3 also had Ca(2+) transport activity and the N-terminal regulatory region (NRR) through yeast complementation assay. Over-expression of GhCAX3 in tobacco showed less sensitivity to ABA during seed germination and seedling stages, and the phenotypic difference between wild type (WT) and transgenic plants was more significant when the NRR was truncated. Furthermore, GhCAX3 conferred cold tolerance in yeast as well as in tobacco seedlings based on physiological and molecular studies. However, transgenic plant seeds showed more sensitivity to cold stress compared to WT during seed germination, especially when expressed in N-terminal truncated version. Finally, the extent of sensitivity in transgenic lines was more severe than that in WT line under sodium tungstate treatment (an ABA repressor), indicating that ABA could alleviate cold sensitivity of GhCAX3 seeds, especially in short of its NRR. Meanwhile, we also found that overexpression of GhCAX3 could enhance some cold and ABA responsive marker genes. Taken together, these results suggested that GhCAX3 plays important roles in the cross-talk of ABA and cold signal transduction, and compared to full-length of GhCAX3, the absence of NRR could enhance the tolerance or sensitivity to cold stress, depending on seedling's developmental stages.

A genetic and metabolic analysis revealed that cotton fiber cell development was retarded by flavonoid naringenin

The cotton (Gossypium spp.) fiber is a unique elongated cell that is useful for investigating cell differentiation. Previous studies have demonstrated the importance of factors such as sugar metabolism, the cytoskeleton, and hormones, which are commonly known to be involved in plant cell development, while the secondary metabolites have been less regarded. By mining public data and comparing analyses of fiber from two cotton species (Gossypium hirsutum and Gossypium barbadense), we found that the flavonoid metabolism is active in early fiber cell development. Different flavonoids exhibited distinct effects on fiber development during ovule culture; among them, naringenin (NAR) could significantly retard fiber development. NAR is a substrate of flavanone 3-hydroxylase (F3H), and silencing the F3H gene significantly increased the NAR content of fiber cells. Fiber development was suppressed following F3H silencing, but the overexpression of F3H caused no obvious effects. Significant retardation of fiber growth was observed after the introduction of the F3H-RNA interference segment into the high-flavonoid brown fiber G. hirsutum T586 line by cross. A greater accumulation of NAR as well as much shorter fibers were also observed in the BC1 generation plants. These results suggest that NAR is negatively associated with fiber development and that the metabolism mediated by F3H is important in fiber development, thus highlighting that flavonoid metabolism represents a novel pathway with the potential for cotton fiber improvement.

Deep sequencing reveals differences in the transcriptional landscapes of fibers from two cultivated species of cotton

Cotton (Gossypium) fiber is the most prevalent natural product used in the textile industry. The two major cultivated species, G. hirsutum (Gh) and G. barbadense (Gb), are allotetraploids with contrasting fiber quality properties. To better understand the molecular basis for their fiber differences, EST pyrosequencing was used to document the fiber transcriptomes at two key development stages, 10 days post anthesis (dpa), representing the peak of fiber elongation, and 22 dpa, representing the transition to secondary cell wall synthesis. The 617,000 high quality reads (89% of the total 692,000 reads) from 4 libraries were assembled into 46,072 unigenes, comprising 38,297 contigs and 7,775 singletons. Functional annotation of the unigenes together with comparative digital gene expression (DGE) revealed a diverse set of functions and processes that were partly linked to specific fiber stages. Globally, 2,770 contigs (7%) showed differential expression (>2-fold) between 10 and 22 dpa (irrespective of genotype), with 70% more highly expressed at 10 dpa, while 2,248 (6%) were differentially expressed between the genotypes (irrespective of stage). The most significant genes with differential DGE at 10 dpa included expansins and lipid transfer proteins (higher in Gb), while at 22 dpa tubulins, cellulose, and sucrose synthases showed higher expression in Gb. DGE was compared with expression data of 10 dpa-old fibers from Affymetrix microarrays. Among 543 contigs showing differential expression on both platforms, 74% were consistent in being either over-expressed in Gh (242 genes) or in Gb (161 genes). Furthermore, the unigene set served to identify 339 new SSRs and close to 21,000 inter-genotypic SNPs. Subsets of 88 SSRs and 48 SNPs were validated through mapping and added 65 new loci to a RIL genetic map. The new set of fiber ESTs and the gene-based markers complement existing available resources useful in basic and applied research for crop improvement in cotton.

GbTCP, a cotton TCP transcription factor, confers fibre elongation and root hair development by a complex regulating system

As the most important natural raw material for textile industry, cotton fibres are an excellent model for studying single-cell development. Although expression profiling and functional genomics have provided some data, the mechanism of fibre development is still not well known. A class I TCP transcription factor (designated GbTCP), encoding 344 amino acids, was isolated from the normalized cDNA library of sea-island cotton fibre (from -2 to 25 days post anthesis). GbTCP was preferentially expressed in the elongating cotton fibre from 5 to 15 days post anthesis. Some expression was also observed in stems, apical buds, and petals. RNAi silencing of GbTCP produced shorter fibre, a reduced lint percentage, and a lower fibre quality than the wild-type plants. Overexpression of GbTCP enhanced root hair initiation and elongation in Arabidopsis and regulated branching. Solexa sequencing and Affymetrix GeneChip analysis indicated that GbTCP positively regulates the level of jasmonic acid (JA) and, as a result, activates downstream genes (reactive oxygen species, calcium signalling, ethylene biosynthesis and response, and several NAC and WRKY transcription factors) necessary for elongation of fibres and root hairs. JA content analysis in cotton also confirmed that GbTCP has a profound effect on JA biosynthesis. In vitro ovule culture showed that an appropriate concentration of JA promoted fibre elongation. The results suggest that GbTCP is an important transcription factor for fibre and root hair development by regulating JA biosynthesis and response and other pathways, including reactive oxygen species, calcium channel and ethylene signalling.

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.

GbPDF1 is involved in cotton fiber initiation via the core cis-element HDZIP2ATATHB2

Cotton (Gossypium spp.) fiber cells are seed trichomes derived from the epidermal layer of the cotton seed coat. The molecular components responsible for regulating fiber cell differentiation have not been fully elucidated. A cotton PROTODERMAL FACTOR1 gene (GbPDF1) was found to be expressed preferentially during fiber initiation and early elongation, with highest accumulation in fiber cells 5 d post anthesis. PDF1 silencing caused retardation of fiber initiation and produced shorter fibers and lower lint percentage compared with the wild type, indicating that the gene is required for cotton fiber development. Further analysis showed that a higher accumulation of hydrogen peroxide occurred in the RNA interference transgenic cotton lines. Meanwhile, the expression of several genes related to ethylene and pectin synthesis or sugar transport during cotton fiber growth was found to be significantly reduced in the PDF1-suppressed cotton. Three proteins interacting with GbPDF1 in yeast and in planta might involve cellular signaling or metabolism. GbPDF1 promoter::GUS constructs in transgenic cotton were predominantly expressed in the epidermis of ovules and developing fibers. Progressive deletions of the GbPDF1 promoter showed that a 236-bp promoter fragment was sufficient for basal GbPDF1 transcription in cotton. Mutation of putative regulatory sequences showed that HDZIP2ATATHB2, an element within the fragment, was essential for PGbPDF1-1 expression. The binding activity between this cis-element and nuclear extracts from fiber-bearing cotton ovules at 5 d post anthesis was specific. We conclude that GbPDF1 plays a critical role together with interaction partners in hydrogen peroxide homeostasis and steady biosynthesis of ethylene and pectin during fiber development via the core cis-element HDZIP2ATATHB2.

Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry

The incompatible pathosystem between resistant cotton (Gossypium barbadense cv. 7124) and Verticillium dahliae strain V991 was used to study the cotton transcriptome changes after pathogen inoculation by RNA-Seq. Of 32,774 genes detected by mapping the tags to assembly cotton contigs, 3442 defence-responsive genes were identified. Gene cluster analyses and functional assignments of differentially expressed genes indicated a significant transcriptional complexity. Quantitative real-time PCR (qPCR) was performed on selected genes with different expression levels and functional assignments to demonstrate the utility of RNA-Seq for gene expression profiles during the cotton defence response. Detailed elucidation of responses of leucine-rich repeat receptor-like kinases (LRR-RLKs), phytohormone signalling-related genes, and transcription factors described the interplay of signals that allowed the plant to fine-tune defence responses. On the basis of global gene regulation of phenylpropanoid metabolism-related genes, phenylpropanoid metabolism was deduced to be involved in the cotton defence response. A closer look at the expression of these genes, enzyme activity, and lignin levels revealed differences between resistant and susceptible cotton plants. Both types of plants showed an increased level of expression of lignin synthesis-related genes and increased phenylalanine-ammonia lyase (PAL) and peroxidase (POD) enzyme activity after inoculation with V. dahliae, but the increase was greater and faster in the resistant line. Histochemical analysis of lignin revealed that the resistant cotton not only retains its vascular structure, but also accumulates high levels of lignin. Furthermore, quantitative analysis demonstrated increased lignification and cross-linking of lignin in resistant cotton stems. Overall, a critical role for lignin was believed to contribute to the resistance of cotton to disease.

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.

Spatiotemporal manipulation of auxin biosynthesis in cotton ovule epidermal cells enhances fiber yield and quality

The capacity of conventional breeding to simultaneously improve the yield and quality of cotton fiber is limited. The accumulation of the plant hormone indole-3-acetic acid (IAA) in cotton fiber initials prompted us to investigate the effects of genetically engineering increased IAA levels in the ovule epidermis. Targeted expression of the IAA biosynthetic gene iaaM, driven by the promoter of the petunia MADS box gene Floral Binding protein 7 (FBP7), increased IAA levels in the epidermis of cotton ovules at the fiber initiation stage. This substantially increased the number of lint fibers, an effect that was confirmed in a 4-year field trial. The lint percentage of the transgenic cotton, an important component of fiber yield, was consistently higher in our transgenic plants than in nontransgenic controls, resulting in a >15% increase in lint yield. Fiber fineness was also notably improved.

Genome structure of cotton revealed by a genome-wide SSR genetic map constructed from a BC1 population between gossypium hirsutum and G. barbadense

Background

Cotton, with a large genome, is an important crop throughout the world. A high-density genetic linkage map is the prerequisite for cotton genetics and breeding. A genetic map based on simple polymerase chain reaction markers will be efficient for marker-assisted breeding in cotton, and markers from transcribed sequences have more chance to target genes related to traits. To construct a genome-wide, functional marker-based genetic linkage map in cotton, we isolated and mapped expressed sequence tag-simple sequence repeats (EST-SSRs) from cotton ESTs derived from the A₁, D₅, (AD)₁, and (AD)₂ genome.

Results

A total of 3177 new EST-SSRs developed in our laboratory and other newly released SSRs were used to enrich our interspecific BC₁ genetic linkage map. A total of 547 loci and 911 loci were obtained from our EST-SSRs and the newly released SSRs, respectively. The 1458 loci together with our previously published data were used to construct an updated genetic linkage map. The final map included 2316 loci on the 26 cotton chromosomes, 4418.9 cM in total length and 1.91 cM in average distance between adjacent markers. To our knowledge, this map is one of the three most dense linkage maps in cotton. Twenty-one segregation distortion regions (SDRs) were found in this map; three segregation distorted chromosomes, Chr02, Chr16, and Chr18, were identified with 99.9% of distorted markers segregating toward the heterozygous allele. Functional analysis of SSR sequences showed that 1633 loci of this map (70.6%) were transcribed loci and 1332 loci (57.5%) were translated loci.

Conclusions

This map lays groundwork for further genetic analyses of important quantitative traits, marker-assisted selection, and genome organization architecture in cotton as well as for comparative genomics between cotton and other species. The segregation distorted chromosomes can be a guide to identify segregation distortion loci in cotton. The annotation of SSR sequences identified frequent and rare gene ontology items on each chromosome, which is helpful to discover functions of cotton chromosomes.

Near-isogenic cotton germplasm lines that differ in fiber-bundle strength have temporal differences in fiber gene expression patterns as revealed by comparative high-throughput profiling

Gene expression profiles of developing cotton (Gossypium hirsutum L.) fibers from two near-isogenic lines (NILs) that differ in fiber-bundle strength, short-fiber content, and in fewer than two genetic loci were compared using an oligonucleotide microarray. Fiber gene expression was compared at five time points spanning fiber elongation and secondary cell wall (SCW) biosynthesis. Fiber samples were collected from field plots in a randomized, complete block design, with three spatially distinct biological replications for each NIL at each time point. Microarray hybridizations were performed in a loop experimental design that allowed comparisons of fiber gene expression profiles as a function of time between the two NILs. Overall, developmental expression patterns revealed by the microarray experiment agreed with previously reported cotton fiber gene expression patterns for specific genes. Additionally, genes expressed coordinately with the onset of SCW biosynthesis in cotton fiber correlated with gene expression patterns of other SCW-producing plant tissues. Functional classification and enrichment analysis of differentially expressed genes between the two NILs revealed that genes associated with SCW biosynthesis were significantly up-regulated in fibers of the high-fiber quality line at the transition stage of cotton fiber development. For independent corroboration of the microarray results, 15 genes were selected for quantitative reverse transcription PCR analysis of fiber gene expression. These analyses, conducted over multiple field years, confirmed the temporal difference in fiber gene expression between the two NILs. We hypothesize that the loci conferring temporal differences in fiber gene expression between the NILs are important regulatory sequences that offer the potential for more targeted manipulation of cotton fiber quality.