Search

Cotton is an important natural fiber crop and economic crop worldwide. The quality of cotton fiber directly determines the quality of cotton textiles. Identifying cotton fiber development-related genes and exploring their biological functions will not only help to better understand the elongation and development mechanisms of cotton fibers but also provide a theoretical basis for the cultivation of new cotton varieties with excellent fiber quality. In this study, RNA sequencing technology was used to construct transcriptome databases for different nonfiber tissues (root, leaf, anther and stigma) and fiber developmental stages (7 days post-anthesis (DPA), 14 DPA, and 26 DPA) of upland cotton Coker 312. The sizes of the seven transcriptome databases constructed ranged from 4.43 to 5.20 Gb, corresponding to approximately twice the genome size of Gossypium hirsutum (2.5 Gb). Among the obtained clean reads, 83.32% to 88.22% could be compared to the upland cotton TM-1 reference genome. By analyzing the differential gene expression profiles of the transcriptome libraries of fiber and nonfiber tissues, we obtained 1205, 1135 and 937 genes with significantly upregulated expression at 7 DPA, 14 DPA and 26 DPA, respectively, and 124, 179 and 213 genes with significantly downregulated expression. Subsequently, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analyses were performed, which revealed that these genes were mainly involved in catalytic activity, carbohydrate metabolism, the cell membrane and organelles, signal transduction and other functions and metabolic pathways. Through gene annotation analysis, many transcription factors and genes related to fiber development were screened. Thirty-six genes were randomly selected from the significantly upregulated genes in fiber, and expression profile analysis was performed using qRT-PCR. The results were highly consistent with the gene expression profile analyzed by RNA-seq, and all of the genes were specifically or predominantly expressed in fiber. Therefore, our RNA sequencing-based comparative transcriptome analysis will lay a foundation for future research to provide new genetic resources for the genetic engineering of improved cotton fiber quality and for cultivating new transgenic cotton germplasms for fiber quality improvement., (© 2021. The Author(s).)

Background: Numerous quantitative trait loci (QTLs) and candidate genes associated with yield-related traits have been identified in cotton by genome-wide association study (GWAS) analysis. However, most of the phenotypic data were from a single or few environments, and the stable loci remained to be validated under multiple field environments., Results: Here, 242 upland cotton accessions collected from different origins were continuously investigated for phenotypic data of four main yield components, including boll weight (BW) and lint percentage (LP) under 13 field environments, and boll number per plant (BN) and seed index (SI) under 11 environments. Correlation analysis revealed a positive correlation between BN and LP, BW and SI, while SI had a negative correlation with LP and BN. Genetic analysis indicated that LP had the highest heritability estimates of 94.97%, followed by 92.08% for SI, 86.09% for BW, and 72.92% for BN, indicating LP and SI were more suitable traits for genetic improvement. Based on 56,010 high-quality single nucleotide polymorphisms (SNPs) and GWAS analysis, a total of 95 non-redundant QTLs were identified, including 12 of BN, 23 of BW, 45 of LP, and 33 of SI, respectively. Of them, 10 pairs of homologous QTLs were detected between A and D sub-genomes. We also found that 15 co-located QTLs with more than two traits and 12 high-confidence QTLs were detected under more than six environments, respectively. Further, two NET genes (GH_A08G0716 and GH_A08G0783), located in a novel QTL hotspot (qtl24, qtl25 and qlt26) were predominately expressed in early fiber development stages, exhibited significant correlation with LP and SI. The GH_A07G1389 in the stable qtl19 region encoded a tetratricopeptide repeat (TPR)-like superfamily protein and was a homologous gene involved in short fiber mutant ligon lintless-y (Li y ), implying important roles in cotton yield., Conclusions: The present study provides a foundation for understanding the regulatory mechanisms of yield components and may enhance yield improvement through molecular breeding in cotton.

Key Message: Combining GWAS, QTL-seq and transcriptome sequencing detected basal defense-related genes showing gDNA sequence variation and expression difference in diverse cotton lines, which might be the molecular mechanisms of VW resistance in G. hirsutum. Verticillium wilt (VW), which is caused by the soil-borne fungus Verticillium dahliae, is a major disease in cotton (Gossypim hirsutum) worldwide. To facilitate the understanding of the genetic basis for VW resistance in cotton, a genome-wide association study (GWAS), QTL-seq and transcriptome sequencing were performed. The GWAS of VW resistance in a panel of 120 core elite cotton accessions using the Cotton 63K Illumina Infinium SNP array identified 5 QTL from 18 significant SNPs meeting the 5% false discovery rate threshold on 5 chromosomes. All QTL identified through GWAS were found to be overlapped with previously reported QTL. By combining GWAS, QTL-seq and transcriptome sequencing, we identified eight candidate genes showing both gDNA sequence variation and expression difference between resistant and susceptible lines, most related to transcription factors (TFs), flavonoid biosynthesis and those involving in the plant basal defense and broad-spectrum disease resistance. Ten KASP markers were successfully validated in diverse cotton lines and could be deployed in marker-assisted breeding to enhance VW resistance. These results supported our inference that the gDNA sequence variation or expression difference of those genes involving in the basal defense in diverse cotton lines might be the molecular mechanisms of VW resistance in G. hirsutum.

Background: In plants, 14-3-3 proteins, also called GENERAL REGULATORY FACTORs (GRFs), encoded by a large multigene family, are involved in protein-protein interactions and play crucial roles in various physiological processes. No genome-wide analysis of the GRF gene family has been performed in cotton, and their functions in flowering are largely unknown., Results: In this study, 17, 17, 31, and 17 GRF genes were identified in Gossypium herbaceum, G. arboreum, G. hirsutum, and G. raimondii, respectively, by genome-wide analyses and were designated as GheGRFs, GaGRFs, GhGRFs, and GrGRFs, respectively. A phylogenetic analysis revealed that these proteins were divided into ε and non-ε groups. Gene structural, motif composition, synteny, and duplicated gene analyses of the identified GRF genes provided insights into the evolution of this family in cotton. GhGRF genes exhibited diverse expression patterns in different tissues. Yeast two-hybrid and bimolecular fluorescence complementation assays showed that the GhGRFs interacted with the cotton FLOWERING LOCUS T homologue GhFT in the cytoplasm and nucleus, while they interacted with the basic leucine zipper transcription factor GhFD only in the nucleus. Virus-induced gene silencing in G. hirsutum and transgenic studies in Arabidopsis demonstrated that GhGRF3/6/9/15 repressed flowering and that GhGRF14 promoted flowering., Conclusions: Here, 82 GRF genes were identified in cotton, and their gene and protein features, classification, evolution, and expression patterns were comprehensively and systematically investigated. The GhGRF3/6/9/15 interacted with GhFT and GhFD to form florigen activation complexs that inhibited flowering. However, GhGRF14 interacted with GhFT and GhFD to form florigen activation complex that promoted flowering. The results provide a foundation for further studies on the regulatory mechanisms of flowering.

The CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION (ESR)-RELATED (CLE) gene family encodes a large number of polypeptide signaling molecules involved in the regulation of shoot apical meristem division and root and vascular bundle development in a variety of plants. CLE family genes encode important short peptide hormones; however, the functions of these signaling polypeptides in cotton remain largely unknown. In the current work, we studied the effects of the CLE family genes on growth and development in cotton. Based on the presence of a conserved CLE motif of 13 amino acids, 93 genes were characterized as GhCLE gene family members, and these were subcategorized into 7 groups. A preliminary analysis of the cotton CLE gene family indicated that the activity of its members tends to be conserved in terms of both the 13-residue conserved domain at the C-terminus and their subcellular localization pattern. Among the 14 tested genes, the ectopic overexpression of GhCLE5::GFP partially mimicked the phenotype of the clv3 mutant in Arabidopsis. GhCLE5 could affect the endogenous CLV3 in binding to the receptor complex, comprised of CLV1, CLV2, and CRN, in the yeast two-hybrid assay and split-luciferase assay. Silencing GhCLE5 in cotton caused a short seedling phenotype. Therefore, we concluded that the cotton GhCLE gene family is functionally conserved in apical shoot development regulation. These results indicate that CLE also plays roles in cotton development as a short peptide hormone.

Most commercially produced cotton cultivars have two types of fibers on the seed coat, short fuzz and long lint. Lint fiber is used in the textile industry, while fuzz is considered an undesirable trait. Both types of fibers are believed to be controlled by the same regulators; however, their mechanisms of actions are still obscure. Cotton fiber mutants provide an excellent system to study the genes that regulate fiber development. Here we described four uncharacterized and three previously reported cotton mutants with fuzzless seed phenotypes. To evaluate whether or not the genes previously associated with fuzzless seed phenotypes have mutations we sequenced whole genomic DNA of seven mutants and wild type varieties. We identified multiple polymorphic changes among the tested genes. Non-synonymous SNPs in the coding region of the MML3-A gene was common in the six mutant lines tested in this study, showing both dominant and recessive fuzzless phenotypes. We have mapped the locus of the causative mutation for one of the uncharacterized fuzzless lines using an F 2 population that originated from a cross between the dominant fuzzless mutant and a wild type. Further, we have clarified the current knowledge about the causative n 2 mutations by analyzing the sequence data and previously reported mapping data. The key genes and possible mechanisms of fiber differentiation are discussed in this study.

Background: The WOX (WUSCHEL-RELATED HOMEOBOX) gene family encodes a class of transcription factors that are unique to green plants, where they are involved in regulating the development of plant tissues and organs by determining cell fate. Although the importance of the WOX gene is well known, there are few studies describing their functions in cotton., Results: In this study, 32 WOX genes were found in Gossypium hirsutum. Phylogenetic analysis showed that WOX proteins of cotton can be divided into three clades: the ancient, intermediate, and WUS clades. The number of WOX proteins in the WUS clade was greater than the sum of the proteins in the other two clades. Our analysis revealed that 20 GhWOX genes are distributed on 16 cotton chromosomes and that duplication events are likely to have contributed to the expansion of the GhWOX family. All GhWOX genes have introns, and each GhWOX protein contains multiple motifs. RNA-seq data and real-time PCR showed that GhWOX13 gene subfamily is specifically expressed at a high level in cotton fibers. We also identified putative GA, NAA, and BR response elements in the promoter regions of the GhWOX13 genes and GhWOX13 transcripts were significantly induced by GA, NAA, and BR., Conclusions: Our data provides a useful resource for future studies on the functional roles of cotton WOX genes and shows that the GhWOX13 genes may influence cotton fiber development. Our results also provide an approach for identifying and characterizing WOX protein genes in other species.

Upland cotton (Gossypium hirsutum L.) race stocks may possess desirable traits for the genetic improvement of cotton. Quantitative trait locus (QTL) analysis can assist in uncovering new alleles from unadapted race stocks. In this study, three sets of chromosome segment introgression lines (ILs) were developed from three backcrosses (BC 3 ) between three race stocks, G. hirsutum races latifolium accs. TX-34 and TX-48 and punctatum acc. TX-114, as donor parents and Texas Marker-1 (TM-1) as the recurrent parent. Based on a total of 452 polymorphic simple sequence repeat (SSR) markers in BC 3 F 2 genotyping, 149, 150 and 184 ILs were obtained from TM-1 × TX-34, TM-1 × TX-48 and TM-1 × TX-114, respectively. The average introgressed chromosomal segment length was 12.7 cM, and the total genetic distance was 3268 cM covering approximately 73.4% of the Upland cotton genome. The BC 3 F 2 , BC 3 F 2:3 and BC 3 F 2:4 progeny, which produced the ILs, were evaluated for yield and fibre quality traits. A total of 128 QTLs were detected, each of which explained 1.6-13.0% of the phenotypic variation. Thirty-five common QTLs related to eight traits were detected. Six QTL clusters were found on five chromosomes. Thirty-eight QTLs were previously unreported, and they may be footprints of cotton domestication. Domestication or artificial selection by humans successfully eliminated most unfavourable QTLs (21/38); however, some favourable QTLs (17/38) are not present in modern cultivars, demonstrating the importance of race stocks for improving cotton cultivars. The 26 elite ILs developed could be used to improve the yield and fibre quality components simultaneously. These results provide information on desirable QTLs for cotton improvement.

Although efforts have been taken to exploit diversity for yield and quality improvements, limited progress on using beneficial alleles in domesticated and undomesticated cotton varieties is limited. Given the complexity and limited amount of genomic information since the completion of four cotton genomes, characterizing significant variations and haplotype block inheritance under artificial selection has been challenging. Here we sequenced Gossypium hirsutum L. cv CRI-12 (the cotton variety with the largest acreage in China), its parental cultivars, and progeny cultivars, which were bred by the different institutes in China. In total, 3.3 million SNPs were identified and 118, 126 and 176 genes were remarkably correlated with Verticillium wilt, salinity and drought tolerance in CRI-12, respectively. Transcriptome-wide analyses of gene expression, and functional annotations, have provided support for the identification of genes tied to these tolerances. We totally discovered 58 116 haplotype blocks, among which 23 752 may be inherited and 1029 may be recombined under artificial selection. This survey of genetic diversity identified loci that may have been subject to artificial selection and documented the haplotype block inheritance and recombination, shedding light on the genetic mechanism of artificial selection and guiding breeding efforts for the genetic improvement of cotton., (© 2018 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Main Conclusion: Cotton GaTOP6B is involved in cellular endoreduplication and a positive response to drought stress via promoting plant leaf and root growth. Drought is deemed as one of adverse conditions that could cause substantial reductions in crop yields worldwide. Since cotton exhibits a moderate-tolerant phenotype under water-deficit conditions, the plant could therefore be used to characterize potential new genes regulating drought tolerance in crop plants. In this work, GaTOP6B, encoding DNA topoisomerase VI subunit B, was identified in Asian cotton (Gossypium arboreum). Virus-induced gene silencing (VIGS) and overexpression (OE) were used to investigate the biological function of GaTOP6B in G. arboreum and Arabidopsis thaliana under drought stress. The GaTOP6B-silencing plants showed a reduced ploidy level, and displayed a compromised tolerance phenotype including lowered relative water content (RWC), decreased proline content and antioxidative enzyme activity, and an increased malondialdehyde (MDA) content under drought stress. GaTOP6B-overexpressing Arabidopsis lines, however, had increased ploidy levels, and were more tolerant to drought treatment, associated with improved RWC maintenance, higher proline accumulation, and reduced stomatal aperture under drought stress. Transcriptome analysis showed that genes involved in the processes like cell cycle, transcription and signal transduction, were substantially up-regulated in GaTOP6B-overexpressing Arabidopsis, promoting plant growth and development. More specifically, under drought stress, the genes involved in the biosynthesis of secondary metabolites such as phenylpropanoid, starch and sucrose were selectively enhanced to improve tolerance in plants. Taken together, the results demonstrated that GaTOP6B could coordinately regulate plant leaf and root growth via cellular endoreduplication, and positively respond to drought stress. Thus, GaTOP6B could be a competent candidate gene for improvement of drought tolerance in crop species.

Background: Auxin-induced genes regulate many aspects of plant growth and development. The Gretchen Hagen 3 (GH3) gene family, one of three major early auxin-responsive families, is ubiquitous in the plant kingdom and its members function as regulators in modulating hormonal homeostasis, and stress adaptations. Specific Auxin-amido synthetase activity of GH3 subfamily II genes is reported to reversibly inactivate or fully degrade excess auxin through the formation of amino acid conjugates. Despite these crucial roles, to date, genome-wide analysis of the GH3 gene family has not been reported in cotton., Results: We identified a total of 10 GH3 subfamily II genes in G. arboreum, 10 in G. raimondii, and 20 in G. hirsutum, respectively. Bioinformatic analysis showed that cotton GH3 genes are conserved with the established GH3s in plants. Expression pattern analysis based on RNA-seq data and qRT-PCR revealed that 20 GhGH3 genes were differentially expressed in a temporally and spatially specific manner, indicating their diverse functions in growth and development. We further summarized the organization of promoter regulatory elements and monitored their responsiveness to treatment with IAA (indole-3-acetic acid), SA (salicylic acid), GA (gibberellic acid) and BL (brassinolide) by qRT-PCR in roots and stems. These hormones seemed to regulate the expression of GH3 genes in both a positive and a negative manner while certain members likely have higher sensitivity to all four hormones. Further, we tested the expression of GhGH3 genes in the BR-deficient mutant pag1 and the corresponding wild-type (WT) of CCRI24. The altered expression reflected the true responsiveness to BL and further suggested possible reasons, at least in part, responsible for the dramatic dwarf and shriveled phenotypes of pag1., Conclusion: We comprehensively identified GH3 subfamily II genes in cotton. GhGH3s are differentially expressed in various tissues/organs/stages. Their response to IAA, SA, BL and GA and altered expression in pag1 suggest that some GhGH3 genes might be simultaneously involved in multiple hormone signaling pathways. Taken together, our results suggest that members of the GhGH3 gene family could be possible candidate genes for mechanistic study and applications in cotton fiber development in addition to the reconstruction of plant architecture.

Background: Basic helix-loop-helix/helix-loop-helix (bHLH/HLH) transcription factors play important roles in plant development. Many reports have suggested that bHLH/HLH proteins participate in brassinosteroid (BR) hormone signaling pathways to promote cell elongation. Cotton fibers are single-cells and derived from seed surface. To explore the roles of bHLH/HLH proteins in cotton fiber development progress by modulating BR signaling pathway, we performed a systematic analysis of the bHLH/HLH gene family in upland cotton (Gossypium hirsutum) genome., Results: In this study, we identified 437 bHLH/HLH genes in upland cotton (G. hirsutum) genome. Phylogenetic analysis revealed that GhbHLH/HLH proteins were split into twenty six clades in the tree. These GhbHLH/HLH genes are distributed unevenly in different chromosomes of cotton genome. Segmental duplication is the predominant gene duplication event and the major contributor for amplification of GhbHLH/HLH gene family. The GhbHLH/HLHs within the same group have conserved exon/intron pattern and their encoding proteins show conserved motif composition. Based on transcriptome data, we identified 77 GhbHLH/HLH candidates that are expressed at relatively high levels in cotton fibers. As adding exogenous BR (brassinolide, BL) or brassinazole (Brz, a BR biosynthesis inhibitor), expressions of these GhbHLH/HLH genes were up-regulated or down-regulated in cotton fibers. Furthermore, overexpression of GhbHLH282 (one of the BR-response genes) in Arabidopsis not only promoted the plant growth, but also changed plant response to BR signaling., Conclusion: Collectively, these data suggested that these GhbHLH/HLH genes may participate in BR signaling transduction during cotton fiber development. Thus, our results may provide a valuable reference data as the basis for further studying the roles of these bHLH/HLH genes in cotton fiber development.

Genes of the CENTRORADIALIS/TERMINAL FLOWER 1/SELF-PRUNING (CETS) family influence meristem identity by controlling the balance between indeterminate and determinate growth, thereby profoundly impacting plant architecture. Artificial selection during cotton (Gossypium hirsutum) domestication converted photoperiodic trees to the day-neutral shrubs widely cultivated today. To understand the regulation of cotton architecture and exploit these principles to enhance crop productivity, we characterized the CETS gene family from tetraploid cotton. We demonstrate that genes of the TERMINAL FLOWER 1 (TFL1)-like clade show different roles in regulating growth patterns. Cotton has five TFL1-like genes: SELF-PRUNING (GhSP) is a single gene whereas there are two TFL1-like and BROTHER OF FT (BFT)-like genes, and these duplications are specific to the cotton lineage. All genes of the cotton TFL1-like clade delay flowering when ectopically expressed in transgenic Arabidopsis, with the strongest phenotypes failing to produce functional flowers. GhSP, GhTFL1-L2, and GhBFT-L2 rescue the early flowering Attfl1-14 mutant phenotype, and the encoded polypeptides interact with a cotton FD protein. Heterologous promoter::GUS fusions illustrate differences in the regulation of these genes, suggesting that genes of the GhTFL1-like clade may not act redundantly. Characterizations of the GhCETS family provide strategies for nuanced control of plant growth.

Background: Plant architecture and the vegetative-reproductive transition have major impacts on the agronomic success of crop plants, but genetic mechanisms underlying these traits in cotton (Gossypium spp.) have not been identified., Results: We identify four natural mutations in GoCEN-D t associated with cluster fruiting (cl) and early maturity. The situ hybridization shows that GhCEN is preferentially expressed in cotton shoot apical meristems (SAM) of the main stem and axillary buds. Constitutive GhCEN-Dt overexpression suppresses the transition of the cotton vegetative apex to a reproductive shoot. Silencing GoCEN leads to early flowering and determinate growth, and in tetraploids causes the main stem to terminate in a floral bud, a novel phenotype that exemplifies co-adaptation of polyploid subgenomes and suggests new research and/or crop improvement approaches. Natural cl variations are enriched in cottons adapted to high latitudes with short frost-free periods, indicating that mutants of GoCEN have been strongly selected for early maturity., Conclusion: We show that the cotton gene GoCEN-Dt, a homolog of Antirrhinum CENTRORADIALIS, is responsible for determinate growth habit and cluster fruiting. Insight into the genetic control of branch and flower differentiation offers new approaches to develop early maturing cultivars of cotton and other crops with plant architecture appropriate for mechanical harvesting.

Key Message: A total of 62 SNPs associated with yield-related traits were identified by a GWAS. Based on significant SNPs, two candidate genes pleiotropically increase lint yield. Improved fibre yield is considered a constant goal of upland cotton (Gossypium hirsutum) breeding worldwide, but the understanding of the genetic basis controlling yield-related traits remains limited. To better decipher the molecular mechanism underlying these traits, we conducted a genome-wide association study to determine candidate loci associated with six yield-related traits in a population of 719 upland cotton germplasm accessions; to accomplish this, we used 10,511 single-nucleotide polymorphisms (SNPs) genotyped by an Illumina CottonSNP63K array. Six traits, including the boll number, boll weight, lint percentage, fruit branch number, seed index and lint index, were assessed in multiple environments; large variation in all phenotypes was detected across accessions. We identified 62 SNP loci that were significantly associated with different traits on chromosomes A07, D03, D05, D09, D10 and D12. A total of 689 candidate genes were screened, and 27 of them contained at least one significant SNP. Furthermore, two genes (Gh_D03G1064 and Gh_D12G2354) that pleiotropically increase lint yield were identified. These identified SNPs and candidate genes provide important insights into the genetic control underlying high yields in G. hirsutum, ultimately facilitating breeding programmes of high-yielding cotton.

Key Message: Thirty significant associations between 22 SNPs and five plant architecture component traits in Chinese upland cotton were identified via GWAS. Four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits. A candidate gene, Gh_D03G0922, might be responsible for plant height in upland cotton. A compact plant architecture is increasingly required for mechanized harvesting processes in China. Therefore, cotton plant architecture is an important trait, and its components, such as plant height, fruit branch length and fruit branch angle, affect the suitability of a cultivar for mechanized harvesting. To determine the genetic basis of cotton plant architecture, a genome-wide association study (GWAS) was performed using a panel composed of 355 accessions and 93,250 single nucleotide polymorphisms (SNPs) identified using the specific-locus amplified fragment sequencing method. Thirty significant associations between 22 SNPs and five plant architecture component traits were identified via GWAS. Most importantly, four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits, and these SNPs were harbored in one linkage disequilibrium block. Furthermore, 21 candidate genes for plant architecture were predicted in a 0.95-Mb region including the four peak SNPs. One of these genes (Gh_D03G0922) was near the significant SNP D03_31584163 (8.40 kb), and its Arabidopsis homologs contain MADS-box domains that might be involved in plant growth and development. qRT-PCR showed that the expression of Gh_D03G0922 was upregulated in the apical buds and young leaves of the short and compact cotton varieties, and virus-induced gene silencing (VIGS) proved that the silenced plants exhibited increased PH. These results indicate that Gh_D03G0922 is likely the candidate gene for PH in cotton. The genetic variations and candidate genes identified in this study lay a foundation for cultivating moderately short and compact varieties in future Chinese cotton-breeding programs.

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

Xyloglucan is a matrix polysaccharide found in the cell walls of all land plants. In growing cells, xyloglucan is thought to connect cellulose microfibrils and regulate their separation during wall extension. Ligon lintless-2 (Li 2 ) is a monogenic dominant cotton fiber mutation that causes extreme reduction in lint fiber length with no pleiotropic effects on vegetative growth. Li 2 represents an excellent model system to study fiber elongation. To understand the role of xyloglucan in cotton fiber elongation we used the short fiber mutant Li 2 and its near isogenic wild type for analysis of xyloglucan content and expression of xyloglucan-related genes in developing fibers. Accumulation of xyloglucan was significantly higher in Li 2 developing fibers than in wild type. Genes encoding enzymes for nine family members of xyloglucan biosynthesis were identified in the draft Gossypium hirsutum genome. RNAseq analysis revealed that most differentially expressed xyloglucan-related genes were down-regulated in Li 2 fiber cells. RT-qPCR analysis revealed that the peak of expression for the majority of xyloglucan-related genes in wild type developing fibers was 5-16days post anthesis (DPA) compared to 1-3 DPA in Li 2 fibers. Thus, our results suggest that early activation of xyloglucan-related genes and down regulation of xyloglucan degradation genes during the elongation phase lead to elevated accumulation of xyloglucan that restricts elongation of fiber cells in Li 2 ., (Copyright © 2017. Published by Elsevier B.V.)

Key Message: Expression analysis of the AG -subfamily members from G. hirsutum during flower and fruit development. Reproductive development in cotton, including the fruit and fiber formation, is a complex process; it involves the coordinated action of gene expression regulators, and it is highly influenced by plant hormones. Several studies have reported the identification and expression of the transcription factor family MADS-box members in cotton ovules and fibers; however, their roles are still elusive during the reproductive development in cotton. In this study, we evaluated the expression profiles of five MADS-box genes (GhMADS3, GhMADS4, GhMADS5, GhMADS6 and GhMADS7) belonging to the AGAMOUS-subfamily in Gossypium hirsutum. Phylogenetic and protein sequence analyses were performed using diploid (G. arboreum, G. raimondii) and tetraploid (G. barbadense, G. hirsutum) cotton genomes, as well as the AG-subfamily members from Arabidopsis thaliana, Petunia hybrida and Antirrhinum majus. qPCR analysis showed that the AG-subfamily genes had high expression during flower and fruit development in G. hirsutum. In situ hybridization analysis also substantiates the involvement of AG-subfamily members on reproductive tissues of G. hirsutum, including ovule and ovary. The effect of plant hormones on AG-subfamily genes expression was verified in cotton fruits treated with gibberellin, auxin and brassinosteroid. All the genes were significantly regulated in response to auxin, whereas only GhMADS3, GhMADS4 and GhMADS7 genes were also regulated by brassinosteroid treatment. In addition, we have investigated the GhMADS3 and GhMADS4 overexpression effects in Arabidopsis plants. Interestingly, the transgenic plants from both cotton AG-like genes in Arabidopsis significantly altered the fruit size compared to the control plants. This alteration suggests that cotton AG-like genes might act regulating fruit formation. Our results demonstrate that members of the AG-subfamily in G. hirsutum present a conserved expression profile during flower development, but also demonstrate their expression during fruit development and in response to phytohormones.

COBRA-Like (COBL) genes, which encode a plant-specific glycosylphosphatidylinositol (GPI) anchored protein, have been proven to be key regulators in the orientation of cell expansion and cellulose crystallinity status. Genome-wide analysis has been performed in A. thaliana, O. sativa, Z. mays and S. lycopersicum, but little in Gossypium. Here we identified 19, 18 and 33 candidate COBL genes from three sequenced cotton species, diploid cotton G. raimondii, G. arboreum and tetraploid cotton G. hirsutum acc. TM-1, respectively. These COBL members were anchored onto 10 chromosomes in G. raimondii and could be divided into two subgroups. Expression patterns of COBL genes showed highly developmental and spatial regulation in G. hirsutum acc. TM-1. Of them, GhCOBL9 and GhCOBL13 were preferentially expressed at the secondary cell wall stage of fiber development and had significantly co-upregulated expression with cellulose synthase genes GhCESA4, GhCESA7 and GhCESA8. Besides, GhCOBL9 Dt and GhCOBL13 Dt were co-localized with previously reported cotton fiber quality quantitative trait loci (QTLs) and the favorable allele types of GhCOBL9 Dt had significantly positive correlations with fiber quality traits, indicating that these two genes might play an important role in fiber development.

For efficient spinning and superior fabric production, long fiber length is a desired trait for cotton production. To unveil the molecular basis of the cotton fiber length regulation, a short fiber mutant, Ligon lintless-1 (Li1), is selected to compare with its corresponding wild type (WT). Li1 is a monogenic dominant cotton mutant causing extremely short fibers (<6mm) on mature seeds with visible pleiotropic effects on vegetative growth and development. In this research, we compared the transcriptome of fiber bearing ovules at 1 DPA, 3 DPA, 8 DPA and leaf between Li1 mutant and WT. A total of 7,852 differentially expressed genes (DEGs) were detected in ovules and leaves, which mainly participated in sugar, secondary metabolite and lipid metabolism pathways based on KEGG analysis. The common DEGs at 1 DPA and 3 DPA were involved in the responses to endogenous stimulus, signal transduction and long-chain fatty acid biosynthesis. For 3 DPA, 8 DPA and leaf, the common DEGs were involved in the responses to auxin and receptor kinases related pathway. Further analysis showed that 37 genes involved in very-long-chain fatty acid biosynthesis were suppressed in Li1 mutant during fiber fast elongation development. Most of the DEGs involved in cell wall metabolism, such cellulose synthase, expansin family, and glycosyl hydrolase were differentially expressed at 3 DPA and 8 DPA. Our results provide new insights into the mechanisms of fiber elongation, and offer novel genes as potential objects for fiber length improvement.

Cotton, an important commercial crop, is cultivated for its natural fibers, and requires an adequate supply of soil nutrients, including phosphorus, for its growth. Soil phosporus exists primarily in insoluble forms. We isolated a mitochondrial malate dehydrogenase (MDH) gene, designated as GhmMDH1, from Gossypium hirsutum L. to assess its effect in enhancing P availability and absorption. An enzyme kinetic assay showed that the recombinant GhmMDH1 possesses the capacity to catalyze the interconversion of oxaloacetate and malate. The malate contents in the roots, leaves and root exudates was significantly higher in GhmMDH1-overexpressing plants and lower in knockdown plants compared with the wild-type control. Knockdown of GhmMDH1 gene resulted in increased respiration rate and reduced biomass whilst overexpression of GhmMDH1 gave rise to decreased respiration rate and higher biomass in the transgenic plants. When cultured in medium containing only insoluble phosphorus, Al-phosphorus, Fe-phosphorus, or Ca-phosphorus, GhmMDH1-overexpressing plants produced significantly longer roots and had a higher biomass and P content than WT plants, however, knockdown plants showed the opposite results for these traits. Collectively, our results show that GhmMDH1 is involved in plant and root growth under phosphorus deficiency conditions in cotton, owing to its functions in leaf respiration and P acquisition.

The growth regulator mepiquat chloride (MC) is globally used in cotton (Gossypium hirsutum L.) canopy manipulation to avoid excess growth and yield loss. However, little information is available as to whether the modification of plant architecture by MC is related to alterations in gibberellic acid (GA) metabolism and signaling. Here, the role of GA metabolism and signaling was investigated in cotton seedlings treated with MC. The MC significantly decreased endogenous GA3 and GA4 levels in the elongating internode, which inhibited cell elongation by downregulating GhEXP and GhXTH2, and then reducing plant height. Biosynthetic and metabolic genes of GA were markedly suppressed within 2-10d of MC treatment, which also downregulated the expression of DELLA-like genes. A remarkable feedback regulation was observed at the early stage of MC treatment when GA biosynthetic and metabolic genes expression was evidently upregulated. Mepiquat chloride action was controlled by temporal translocation and spatial accumulation which regulated GA biosynthesis and signal expression for maintaining GA homeostasis. The results suggested that MC application could reduce endogenous GA levels in cotton through controlled GA biosynthetic and metabolic genes expression, which might inhibit cell elongation, thereby shortening the internode and reducing plant height., (Copyright © 2014 Elsevier Ireland Ltd. All rights reserved.)

Cotton is an important economic crop worldwide. Due to its long growth cycle, large genome size and recalcitrance to stable transformation, traditional methods for the analysis of gene function in this crop are difficult and labour intensive. Here, we report a cotton leaf crumple virus (CLCrV)-based vector and its application in gene function analysis through virus-induced gene silencing (VIGS) and overexpression of microRNAs (miRNAs), small tandem target mimic (STTM) and artificial miRNA (amiRNA) in cotton via an Agrobacterium-mediated infiltration approach. Using this system, we were able to efficiently silence two endogenous genes, magnesium chelatase subunit I (CHLI) and elongation factor-1α (EF-1α), in Gossypium species and the Bacillus thuringiensis cry1A gene in transgenic cotton. Furthermore, our results show that this vector can be used to ectopically express endogenous miR156 in G. hirsutum, causing a reduction in miR156-targeted RNA transcripts resulting in the development of abnormal leaf phenotypes. Ectopic expression of miR165/166 STTM with this vector led to downward curling and crumpled leaves, and a significant increase in the miR165/166 target mRNAs. This versatile system is easy to use and can provide more uniform and persistent gene silencing in cotton, thereby providing a powerful approach for gene discovery in cotton., (© 2014 Society for Experimental Biology, Association of Applied Biologists and John Wiley & Sons Ltd.)

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., (© 2014 The Authors. New Phytologist © 2014 New Phytologist Trust.)

In this study, a cDNA, GhMADS14, encoding a typical MADS-box protein with 223 amino acids was isolated from a cotton cDNA library. Fluorescent microscopy indicated that the GhMADS14 protein was localized in the cell nucleus. GhMADS14 was specifically expressed in the elongating fibers, and its expression was gradually enhanced at early stages of fiber elongation and reached its peak in 9-10 DPA fibers. Overexpression of GhMADS14 in Arabidopsis hindered plant growth. Measurement and statistical analysis revealed that hypocotyl length of GhMADS14 transgenic seedlings was significantly reduced, and the height of the mature transgenic plants was remarkably less than that of the wild type. Furthermore, expression of GA 20-oxidase (AtGA20ox1 and AtGA20ox2) and GA 3-oxidase (AtGA3ox1 and AtGA3ox2) genes was remarkably reduced, whereas AtGA2ox1 and AtGA2ox8 were dramatically up-regulated in the transgenic plants, compared with the wild type. These results suggested that overexpression of GhMADS14 in Arabidopsis may alter expression levels of the genes related to GA biosynthetic and metabolic pathways, resulting in the reduction of endogenous GA amounts in cells. As a result, the transgenic plants grew slowly and display a GA-deficient phenotype., (Copyright © 2013 Elsevier Masson SAS. All rights reserved.)

Key Message: A major leaf shape locus (L) was mapped with molecular markers and genomically targeted to a small region in the D-genome of cotton. By using expression analysis and candidate gene mapping, two LMI1 -like genes are identified as possible candidates for leaf shape trait in cotton. Leaf shape in cotton is an important trait that influences yield, flowering rates, disease resistance, lint trash, and the efficacy of foliar chemical application. The leaves of okra leaf cotton display a significantly enhanced lobing pattern, as well as ectopic outgrowths along the lobe margins when compared with normal leaf cotton. These phenotypes are the hallmark characteristics of mutations in various known modifiers of leaf shape that culminate in the mis/over-expression of Class I KNOX genes. To better understand the molecular and genetic processes underlying leaf shape in cotton, a normal leaf accession (PI607650) was crossed to an okra leaf breeding line (NC05AZ21). An F2 population of 236 individuals confirmed the incompletely dominant single gene nature of the okra leaf shape trait in Gossypium hirsutum L. Molecular mapping with simple sequence repeat markers localized the leaf shape gene to 5.4 cM interval in the distal region of the short arm of chromosome 15. Orthologous mapping of the closely linked markers with the sequenced diploid D-genome (Gossypium raimondii) tentatively resolved the leaf shape locus to a small genomic region. RT-PCR-based expression analysis and candidate gene mapping indicated that the okra leaf shape gene (L (o) ) in cotton might be an upstream regulator of Class I KNOX genes. The linked molecular markers and delineated genomic region in the sequenced diploid D-genome will assist in the future high-resolution mapping and map-based cloning of the leaf shape gene in cotton.

MYB transcription factors are essential for cotton fiber development. We isolated the R2R3-MYB gene GbMYB2 from Gossypium barbadense. RNA in situ hybridization analysis showed that GbMYB2 is mainly expressed in the outer integuments of cotton ovules and in elongating fibers. GbMYB2 expression increased throughout fiber initiation and during the elongation stage. The expression level of GbMYB2 was higher in the Gossypium hirsutum cultivar Xu142 than in the Xu142 (fl) mutant. Overexpression of GbMYB2 in Arabidopsis caused thicker leaf trichomes and longer roots to develop due to the activation of trichome development-related genes such as GL2. These results indicate that GbMYB2 is an R2R3-MYB gene that is involved in fiber development., (Copyright © 2013 Elsevier Masson SAS. All rights reserved.)

A subtractive cDNA library from cotton buds was constructed to prospect for differentially expressed genes related to early bud development. A library was constructed and 768 cDNA sequences were obtained, comprising 168 clusters, with 126 contigs and 42 singlets. Both the Gossypium as well as Arabidopsis databases were utilized for the in silico analysis, since some genes identified in cotton have not yet been studied for functionality, although they have homology with genes from other species. The transcriptome revealed a large number of transcripts, some of them with unknown function, and others related to pollen development, pollen tubes, ovules, and fibers at different stages. The most populated contig was identified as fiber from 0-10 days after anthesis, with 12 reads. The success and novelty rates generated from the library were 67 and 51%, respectively. The information obtained here will provide a framework for research on functional cotton genomics.

Gossypium hirsutum and G. barbadense are two cultivated tetraploid cotton species with differences in fibre quality. The fibre of G. barbadense is longer, stronger and finer than that of G. hirsutum. To isolate genes expressed differently between the two species during fibre development, cDNA-SRAP (sequence-related amplified polymorphism) was applied. This technique was used to analyse genes at different stages of fibre development in G. hirsutum cv. Emian22 and G. barbadense acc. 3-79, the parents of our interspecific mapping population. A total of 4096 SRAP primer combinations were used to screen polymorphism between the DNA of the parents, and 275 highly polymorphic primers were picked out to analyse DNA and RNA from leaves and fibres at different developmental stages of the parents. A total of 168 DNA fragments were isolated from gels and sequenced: 54, 30, 38 and 41 from fibres of 5, 10, 15 and 20 days post-anthesis, respectively, and five from multi stages. To genetically map these sequences, 104 sequence-specific primers were developed and were used to screened polymorphism between the mapping parents. Finally, six markers were mapped on six chromosomes of our backbone interspecific genetic map. This work can give us a primary knowledge of differences in mechanism of fibre development between G. hirsutum and G. barbadense.

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.

A cotton fiber cDNA, GhRING1, and its corresponding gene have been cloned and characterized. The GhRING1 gene encodes a RING-type ubiquitin ligase (E3) containing 338 amino acids (aa). The GhRING1 protein contains a RING finger motif with conserved cysteine and histine residues at the C-terminus, and is classified as a C(3)H(2)C(3)-type RING protein. Blast searches show that GhRING1 has the highest homology to At3g19950, a zinc finger family protein from Arabidopsis. Real time RT-PCR analysis indicates that the GhRING1 gene is expressed in cotton fibers in a developmental manner. The transcript level of GhRING1 gene reaches a maximum in elongating fibers at 15 days post-anthesis (DPA). In vitro auto-ubiquitination assays using wheat germ extract and a reconstitution system demonstrate that GhRING1 has the ubiquitin E3 ligase activity. The histochemical GUS assay was performed to analyze tissue specificity of the GhRING1 and At3g19950 promoters in transgenic Arabidopsis plants. The GUS assay shows that the promoter of At3g19950 is highly activated in leaves, roots, trichomes, and also in anthers and stigma of flowers. In contrast, the GUS expression directed by the GhRING1 promoter is only located at stipules and anthers. The expression pattern of GhRING1 suggests that protein ubiquitination and turnover may be involved in transition to different stages of cotton fiber development.

Mitogen-activated protein kinase (MAPK) cascades play a pivotal role in environmental responses and developmental processes in plants. Previous researches mainly focus on the MAPKs in groups A and B, and little is known on group C. In this study, we isolated and characterized GhMPK7, which is a novel gene from cotton belonging to the group C MAPK. RNA blot analysis indicated that GhMPK7 transcript was induced by pathogen infection and multiple defense-related signal molecules. Transgenic Nicotina benthamiana overexpressing GhMPK7 displayed significant resistance to fungus Colletotrichum nicotianae and virus PVY, and the transcript levels of SA pathway genes were more rapidly and strongly induced. Furthermore, the transgenic N. benthamiana showed reduced ROS-mediated injuries by upregulating expression of oxidative stress-related genes. Interestingly, the transgenic plants germinated earlier and grew faster in comparison to wild-type plants. beta-glucuronidase activity driven by the GhMPK7 promoter was detected in the apical meristem at the vegetative stage, and it was enhanced by treatments with signal molecules and phytohormones. These results suggest that GhMPK7 might play an important role in SA-regulated broad-spectrum resistance to pathogen infection, and that it is also involved in regulation of plant growth and development.

Cotton fiber elongation, largely achieved by cell wall loosening, is an important stage during cotton fiber development. In this present research, a fiber preferential cDNA encoding a pectate lyase (PEL) which could exclusively degrade the de-esterified pectin was isolated from a cotton (Gossypium hirsutum) fiber cDNA library. Subsequently, the corresponding PEL genes were isolated from four different cotton species and characterized. In vitro enzyme assays indicated that GhPEL really exhibited cleavage-activity against de-esterified pectin. The temporal-spatial expression analyses revealed that the GhPEL gene was preferentially expressed in fibers at 10 days-post anthesis (DPA). Antisense GhPEL transgenic cotton plants were generated by Agrobacterium-mediated transformation. Six homozygous lines, each with one or two copies of the transgene inserted as determined by southern blot analysis of the NPTII gene, were selected for further functional analysis. The GhPEL expression during fiber elongation in these transgenic lines was significantly suppressed in various degrees. Furthermore, the reduction of GhPEL enzymatic activity by decreasing GhPEL transcripts severely affected the degradation of de-esterified pectin in primary cell walls of transgenic cotton fibers, which consequently blocked cell wall loosening in early fiber development. Ultimately, the fiber elongation of all these transgenic lines was repressed. These results suggested that GhPEL may play an important role in the process of normal fiber elongation in cotton.

Multiple cellular pathways have been shown to be involved during fiber initiation and elongation stages in the cultivated allotetraploid cotton (Gossypium hirsutum). The cell wall enzymes xyloglucan endotransglycosylase/hydrolases (XTH) have been reported to be associated with the biosynthesis of the cell wall and the growth of cotton fibers, probably regulating the plasticity of the primary cell wall. Among various cotton fiber cDNAs found to be preferentially expressed in cotton fibers, a xyloglucan endotransglycosylase (XTH) cDNA was significantly up-regulated during the elongation stage of cotton fiber development. In the present study, we isolated and characterized genomic clones encoding cotton XTH from cultivated cotton (Gossypium hirsutum) and its diploid progenitors (Gossypium arboreum and Gossypium raimondii), designated GhXTH1-1, GhXTH1-2, GaXTH1 and GrXTH, respectively. In addition, we isolated and characterized, by in silico methods, the putative promoter of XTH1 from Gossypium hirsutum. Sequence analysis revealed more than 50% homology to XTH's at the protein level. DNA gel blot hybridization indicated that at least two copies of GhXTH1 are present in Gossypium hirsutum whereas the diploid progenitor species Gossypium arboreum and Gossypium raimondii has only a single copy. Quantitative real-time PCR and high-resolution melting experiments indicated that in Gossypium hirsutum cultivars, in cotton fibers during early stages of fiber elongation specifically expressing only the GhXTH1-1 gene and expression levels of GhXTH1-1 in fibers varies among cultivars differing in fiber percentage and fiber length.

The gene GhCPK1 encoding a calcium dependent protein kinase was identified from cotton. Transcripts of GhCPK1 accumulated primarily in the elongating fiber, and Arabidopsis plants transformed with GhCPK1 promoter-GUS construct exhibited GUS activity mainly in the developing trichomes, roots, young leaves and sepals. In the bombarded onion epidermal cells, GhCPK1-GFP fusion proteins showed a subcellular distribution in the plasma membrane. In vitro assays indicated that GhCPK1 was a functional calcium-dependent kinase able to undergo autophosphorylation and phosphorylation of the known substrate histone III-S. Together, these results suggest that GhCPK1 may play a role in the calcium signaling events associated with fiber elongation.

One of the impediments in the genetic improvement of cotton fiber is the paucity of information about genes associated with fiber development. Availability of chromosome arm substitution line CS-B22sh (chromosome 22 short arm substitution from 3-79 (Gossypium barbadense) into a TM-1 (Gossypium hirsutum) background) provides a novel opportunity to study fiber-associated genes because previous studies revealed this line was associated with some superior fiber quality traits compared to TM-1. We used an integrated approach of suppression subtractive hybridization (SSH), microarray, and real-time reverse transcription-polymerase chain reaction (RT-PCR) technologies to identify the potential genes associated with fiber development. Utilizing mRNAs from 15 days post-anthesis (dpa) fibers, we constructed a SSH cDNA library with chromosome substitution line CS-B22sh as the tester and TM-1 as the driver. The SSH cDNA library was screened using microarrays. Microarray analysis showed that 36 genes were differentially expressed in CS-B22sh 15-dpa fiber compared to TM-1 as confirmed by real time RT-PCR. These genes include two beta-tubulins, an actin, a putative kinesin light chain, a cellulose synthase, glycosyl hydrolase family 3 protein, pyruvate decarboxylase, glycoside hydrolase family 5, GDP-mannose pyrophosphorylase, dynamin-like protein, annexin, and a number of genes involved in signal transduction, and protein, nucleic acid, and lipid metabolisms. To our knowledge, this is the first report on identification of differentially expressed fiber-associated genes in a cotton chromosomal substitution line.

The possible role of methylation in the performance of heterosis has been analyzed in many crops. To further study this possibility, we investigated both the differences in cytosine methylation patterns between cotton heterotic hybrid/nonheterotic hybrids and their parental lines and the change in methylation level from seedling stage to flowering stage by using the methylation-sensitive amplified polymorphism (MSAP) method. The results showed that the number of demethylation loci in highly heterotic hybrids was greater that in lowly heterotic hybrids, and the level of DNA cytosine methylation in cotton at the seedling stage is higher than that at the flowering stage. The altered methylation patterns at low-copy genomic regions can be confirmed by DNA gel blot analysis. A total of 39 fragments that showed different methylation patterns were cloned and sequenced. The methylation status of these genes was modified differentially in hybrid and parents, suggesting that these genes might play a role in the performance of heterosis.

cDNA amplified fragment length polymorphism (cDNA-AFLP) analysis was used to investigate the differentially expressed genes between sterile and fertile plants of ms5ms6 double-recessive genic male sterility (GMS) two-type line cotton (Gossypium hirsutum L.) at different stages, i.e., sporogenous cell stage, pollen mother cell (PMC) stage, and pollen grain stage. Seventeen differentially expressed fragments were identified. Functional analysis indicated that their corresponding genes may participate in the processes of signal transduction, transcription, energy metabolism, and plant cell wall development. Northern blot demonstrated the credibility of the result of cDNA-AFLP. A sterility restorer factor-like gene, which only expressed in fertile anther and was notably homologous to T cytoplasm male sterility restorer factor 2 of maize (Zea mays L.), was identified in this research.

A novel cDNA encoding DRE-binding transcription factor, designated GhDBP3, was cloned from Gossypium hirsutum. This protein was classified into A-4 group of DREB subfamily based on multiple sequence alignment and phylogenetic characterization. Semiquantitative RT-PCR showed that GhDBP3 was expressed in the leaves, cotyledons, roots and stems of 2-week-old cotton seedlings under non-stress conditions and was greatly induced in the cotton cotyledons by drought, NaCl, low temperature and ABA treatment. EMSA revealed that GhDBP3 was able to bind to the DRE cis-element in vitro. Transient assay using the particle bombardment method showed that GhDBP3 was a transcriptional activator, capable of activating expression of a reporter gene driven by the LEA D113 promoter containing a DRE like sequence in tobacco cells. Our results indicate that GhDBP3 could be a new member of DRE-binding transcription factor family and may play an important role in response to ABA and environmental stresses.

Cotton fibres are a subset of single epidermal cells that elongate from the seed coat to produce the long cellulose strands or lint used for spinning into yarn. To identify genes that might regulate lint fibre initiation, expression profiles of 0 days post-anthesis (dpa) whole ovules from six reduced fibre or fibreless mutants were compared with wild-type linted cotton using cDNA microarrays. Numerous clones were differentially expressed, but when only those genes that are normally expressed in the ovule outer integument (where fibres develop) were considered, just 13 different cDNA clones were down-regulated in some or all of the mutants. These included: a Myb transcription factor (GhMyb25) similar to the Antirrhinum Myb AmMIXTA, a putative homeodomain protein (related to Arabidopsis ATML1), a cyclin D gene, some previously identified fibre-expressed structural and metabolic genes, such as lipid transfer protein, alpha-expansin and sucrose synthase, as well as some unknown genes. Laser capture microdissection and reverse transcription-PCR were used to show that both the GhMyb25 and the homeodomain gene were predominantly ovule specific and were up-regulated on the day of anthesis in fibre initials relative to adjacent non-fibre ovule epidermal cells. Their spatial and temporal expression pattern therefore coincided with the time and location of fibre initiation. Constitutive overexpression of GhMyb25 in transgenic tobacco resulted in an increase in branched long-stalked leaf trichomes. The involvement of cell cycle genes prompted DNA content measurements that indicated that fibre initials, like leaf trichomes, undergo DNA endoreduplication. Cotton fibre initiation therefore has some parallels with leaf trichome development, although the detailed molecular mechanisms are clearly different.

Single-celled cotton fiber (Gossypium hirsutum) provides a unique experimental system to study cell elongation. To investigate the role of the actin cytoskeleton during fiber development, 15 G. hirsutum ACTIN (GhACT) cDNA clones were characterized. RNA gel blot and real-time RT-PCR analysis revealed that GhACT genes are differentially expressed in different tissues and can be classified into four groups. One group, represented by GhACT1, is expressed predominantly in fiber cells and was studied in detail. A 0.8-kb GhACT1 promoter sufficient to confirm its fiber-specific expression was identified. RNA interference of GhACT1 caused significant reduction of its mRNA and protein levels and disrupted the actin cytoskeleton network in fibers. No defined actin network was observed in these fibers and, consequently, fiber elongation was inhibited. Our results suggested that GhACT1 plays an important role in fiber elongation but not fiber initiation.

Plant Rac proteins belong to an important group of signal switches anchoring on membranes, involved in various physiological processes including cell polar growth, synthesis of secondary wall, resistance response and hormone signaling. In the attempt to elucidate the molecular mechanism of initiation and elongation of cotton fiber, two cotton Rac protein genes, designated as GhRacA and GhRacB, were amplified from elongating fibers and cloned. It was demonstrated that, the cDNA of GhRacA contained 959 bp and encoded a putative polypepetide of 211 aa, while GhRacB was 920 bp in length, encoding a predicted protein of 195 aa. These two cotton Rac proteins, GhRacA and GhRacB, contained conserved regions involved in GTP/GDP binding and activation, an effector region and a polybasic region. A conserved prenylation site CSIL was found in GhRacB, while no apparent prenylation site was discovered in GhRacA. Sequence comparisons showed that GhRacA and GhRacB were two novel Rac proteins from cotton. The expression patterns of GhRacA and GhRacB was analyzed by RT-PCR. It was demonstrated that these two Rac protein genes were both expressed in root, hypocotyls, stem, leaf and fibers, and the highest level of transcripts was to accumulate in the fibers at the stage of initiation and elongation, suggesting that the two Rac genes, GhRacA and GhRacB, might play an important role in the early stage of fiber development.

The fibers of cotton (Gossypium hirsutum) are single-cell trichomes that undergo rapid and synchronous elongation. Cortical microtubules provide spatial information necessary for the alignment of cellulose microfibrils that confine and regulate cell elongation. We used gene-specific probes to investigate alpha-tubulin transcript levels in elongating cotton fibers. Two discrete patterns of transcript accumulation were observed. Whereas transcripts of alpha-tubulin genes GhTua2/3 and GhTua4 increased in abundance from 10 to 20 d post anthesis (DPA), GhTua1 and GhTua5 transcripts were abundant only through to 14 DPA, and dropped significantly at 16 DPA with the onset of secondary wall synthesis. This is the first report, to our knowledge, of gene-specific changes in tubulin transcript levels during the development of a terminally differentiated plant cell. The decrease in abundance of GhTua1 and GhTua5 transcripts was correlated with pronounced changes in cell wall structure, suggesting that alpha-tubulin isoforms may be functionally distinct in elongating fiber cells. Although total alpha-tubulin transcript levels were much higher in fiber than several other tissues, including the hypocotyl and pollen, none of the alpha-tubulins was specific to fiber cells.

A gene (G9) expressed during late microsporogenesis in cotton (Gossypium hirsutum L.) was isolated. Sequence analysis of the cDNA (1.3 kb) as well as the gene (2.6 kb) revealed an open reading frame of 1233 bases encoding a protein of 43.9 kDa. The coding region of the gene is interrupted by three introns. Northern analysis of the RNA from developing anthers showed that the transcripts appear 12 days before anthesis and that the maximal concentration of RNA occurs in pollen on the day of anthesis. This pattern of gene expression suggests functions in post-anthesis events. Sequence comparisons with other known plant genes indicated that G9 is homologous to polygalacturonases. The G9 promoter conferred tissue and temporal specificity of beta-glucuronidase (GUS) expression in transgenic tobacco plants. Thus, the G9 promoter can be used to drive gene expression in homologous as well as heterologous plants in a tissue-specific manner.