Cotton (Gossypium hirsutum) VIRMA as an N6-Methyladenosine RNA Methylation Regulator Participates in Controlling Chloroplast-Dependent and Independent Leaf Development.

N6-methyladenosine (m6A) is one of the most abundant internal modifications of mRNA, which plays important roles in gene expression regulation, and plant growth and development. Vir-like m6A methyltransferase associated (VIRMA) serves as a scaffold for bridging the catalytic core components of the m6A methyltransferase complex. The role of VIRMA in regulating leaf development and its related mechanisms have not been reported. Here, we identified and characterized two upland cotton (Gossypium hirsutum) VIRMA genes, named as GhVIR-A and GhVIR-D, which share 98.5% identity with each other. GhVIR-A and GhVIR-D were ubiquitously expressed in different tissues and relatively higher expressed in leaves and main stem apexes (MSA). Knocking down the expression of GhVIR genes by the virus-induced gene silencing (VIGS) system influences leaf cell size, cell shape, and total cell numbers, thereby determining cotton leaf morphogenesis. The dot-blot assay and colorimetric experiment showed the ratio of m6A to A in mRNA is lower in leaves of GhVIR-VIGS plants compared with control plants. Messenger RNA (mRNA) high-throughput sequencing (RNA-seq) and a qRT-PCR experiment showed that GhVIRs regulate leaf development through influencing expression of some transcription factor genes, tubulin genes, and chloroplast genes including photosystem, carbon fixation, and ribosome assembly. Chloroplast structure, chlorophyll content, and photosynthetic efficiency were changed and unsuitable for leaf growth and development in GhVIR-VIGS plants compared with control plants. Taken together, our results demonstrate GhVIRs function in cotton leaf development by chloroplast dependent and independent pathways.

The cotton (Gossypium hirsutum) NAC transcription factor (FSN1) as a positive regulator participates in controlling secondary cell wall biosynthesis and modification of fibers.

Cotton (Gossypium hirsutum) fibers are the highly elongated and thickened single-cell trichomes on the seed epidermis. However, little is known about the molecular base of fiber cell wall thickening in detail. In this study, a cotton NAC transcription factor (GhFSN1) that is specifically expressed in secondary cell wall (SCW) thickening fibers was functionally characterized. The GhFSN1 transgenic cotton plants were generated to study how FSN1 regulates fiber SCW formation. Up-regulation of GhFSN1 expression in cotton resulted in an increase in SCW thickness of fibers but a decrease in fiber length. Transcriptomic analysis revealed that GhFSN1 activates or represses numerous downstream genes. GhFSN1 has the ability to form homodimers, binds to its promoter to activate itself, and might be degraded by the ubiquitin-mediated proteasome pathway. The direct targets of GhFSN1 include the fiber SCW-related GhDUF231L1, GhKNL1, GhMYBL1, GhGUT1 and GhIRX12 genes. GhFSN1 binds directly to a consensus sequence (GhNBS), (C/T)(C/G/T)TN(A/T)(G/T)(A/C/G)(A/G)(A/T/G)(A/T/G)AAG, which exists in the promoters of these SCW-related genes. Our data demonstrate that GhFSN1 acts as a positive regulator in controlling SCW formation of cotton fibers by activating its downstream SCW-related genes. Thus, these findings give us novel insights into comprehensive understanding of GhFSN1 function in fiber development.

A R2R3-MYB transcription factor that is specifically expressed in cotton (Gossypium hirsutum) fibers affects secondary cell wall biosynthesis and deposition in transgenic Arabidopsis.

Secondary cell wall (SCW) is an important industrial raw material for pulping, papermaking, construction, lumbering, textiles and potentially for biofuel production. The process of SCW thickening of cotton fibers lays down the cellulose that will constitute the bulk (up to 96%) of the fiber at maturity. In this study, a gene encoding a MYB-domain protein was identified in cotton (Gossypium hirsutum) and designated as GhMYBL1. Quantitative real-time polymerase chain reaction (RT-PCR) analysis revealed that GhMYBL1 was specifically expressed in cotton fibers at the stage of secondary wall deposition. Further analysis indicated that this protein is a R2R3-MYB transcription factor, and is targeted to the cell nucleus. Overexpression of GhMYBL1 in Arabidopsis affected the formation of SCW in the stem xylem of the transgenic plants. The enhanced SCW thickening also occurred in the interfascicular fibers, xylary fibers and vessels of the GhMYBL1-overexpression transgenic plants. The expression of secondary wall-associated genes, such as CesA4, CesA7, CesA8, PAL1, F5H and 4CL1, were upregulated, and consequently, cellulose and lignin biosynthesis were enhanced in the GhMYBL1 transgenic plants. These data suggested that GhMYBL1 may participate in modulating the process of secondary wall biosynthesis and deposition of cotton fibers.

A fasciclin-like arabinogalactan protein, GhFLA1, is involved in fiber initiation and elongation of cotton.

Arabinogalactan proteins (AGPs) are involved in many aspects of plant development. In this study, biochemical and genetic approaches demonstrated that AGPs are abundant in developing fibers and may be involved in fiber initiation and elongation. To further investigate the role of AGPs during fiber development, a fasciclin-like arabinogalactan protein gene (GhFLA1) was identified in cotton (Gossypium hirsutum). Overexpression of GhFLA1 in cotton promoted fiber elongation, leading to an increase in fiber length. In contrast, suppression of GhFLA1 expression in cotton slowed down fiber initiation and elongation. As a result, the mature fibers of the transgenic plants were significantly shorter than those of the wild type. In addition, expression levels of GhFLAs and the genes related to primary cell wall biosynthesis were remarkably enhanced in the GhFLA1 overexpression transgenic fibers, whereas the transcripts of these genes were dramatically reduced in the fibers of GhFLA1 RNA interference plants. An immunostaining assay indicated that both AGP composition and primary cell wall composition were changed in the transgenic fibers. The levels of glucose, arabinose, and galactose were also altered in the primary cell wall of the transgenic fibers compared with those of the wild type. Together, our results suggested that GhFLA1 may function in fiber initiation and elongation by affecting AGP composition and the integrity of the primary cell wall matrix.

Cotton KNL1, encoding a class II KNOX transcription factor, is involved in regulation of fibre development.

In this study, the GhKNL1 (KNOTTED1-LIKE) gene, encoding a classical class II KNOX protein was identified in cotton (Gossypium hirsutum). GhKNL1 was preferentially expressed in developing fibres at the stage of secondary cell wall (SCW) biosynthesis. GhKNL1 was localized in the cell nucleus, and could interact with GhOFP4, as well as AtOFP1, AtOFP4, and AtMYB75. However, GhKNL1 lacked transcriptional activation activity. Dominant repression of GhKNL1 affected fibre development of cotton. The expression levels of genes related to fibre elongation and SCW biosynthesis were altered in transgenic fibres of cotton. As a result, transgenic cotton plants produced aberrant, shrunken, and collapsed fibre cells. Length and cell-wall thickness of fibres of transgenic cotton plants were significantly reduced compared with the wild type. Furthermore, overexpression and dominant repression of GhKNL1 in Arabidopsis resulted in a reduction in interfascicular fibre cell-wall thickening of basal stems of transgenic plants. Complementation revealed that GhKNL1 rescued the defective phenotype of Arabidopsis knat7 mutant in some extent. These data suggest that GhKNL1, as a transcription factor, participates in regulating fibre development of cotton.

Cotton photosynthesis-related PSAK1 protein is involved in plant response to aphid attack.

It is believed that hundreds of genes, including photosynthesis-related genes, are typically involved in plant response to aphid feeding. Up to now, however, it is little known on the relationship between the photosynthesis-related genes and plant response to herbivores. In this study, we identified a cotton photosynthesis-related gene (GhPSAK1) which belongs to PSI-PSAK family and encodes a putative protein of 162 amino acids. RT-PCR analysis revealed that GhPSAK1 transcripts in leaves were increased at 12-24 h, but decreased at 48-72 h after cotton aphid attack or wounding induction. Choice assay and no-choice assay demonstrated that overexpression of GhPSAK1 in Arabidopsis improved plant tolerance to green peach aphids (Myzus persicae). The defense response genes related to salicylic acid signaling pathway were enhanced in the GhPSAK1 overexpressing transgenic plants. In addition, the callose amount in transgenic Arabidopsis leaves was more than that of wild type. Contents of the soluble sugars and total amino acids were also altered in leaves of transgenic Arabidopsis plants. Activities of superoxide dismutase and peroxidase in transgenic leaves were higher than those of wild type. These results suggested that GhPSAK1 may be involved in regulation of cotton response and tolerance to aphid attack.

Molecular characterization of cotton C-repeat/dehydration-responsive element binding factor genes that are involved in response to cold stress.

Low temperature, drought and salinity are major abiotic stresses that influence survival, productivity and geographical distribution of many important crops across the globe. The C-repeat/dehydration-responsive element binding transcription factors (CBF/DREB) are important proteins involved in response to abiotic stresses in plants. In this study, twenty-one CBF genes were identified in cotton (Gossypium hirsutum) by bioinformatic approach. The twenty-one CBF genes (named as GhCBF1--GhCBF21) were characterized to encode proteins that share high similarity with those plant cold stress-related CBF proteins, which contain the classic AP2 domain of 58 amino acid residues. Phylogenetic analysis revealed that the isolated cotton CBF genes can be classified into 4 groups: GhCBF I, GhCBF II, GhCBF III and GhCBF IV. RT-PCR analysis indicated that GhCBF genes were up-regulated in cotton plants under cold stress. Furthermore, four GhCBF genes were up-regulated in cotton under salinity and drought treatments. Our data provided valuable information for further exploring the roles of the CBF genes in cotton development and in response to cold stress.

Genome-wide identification and analysis of DEAD-box RNA helicases in Gossypium hirsutum.

DEAD-box RNA helicases, a prominent subfamily within the RNA helicase superfamily 2 (SF2), play crucial roles in the growth, development, and abiotic stress responses of plants. This study identifies 146 DEAD-box RNA helicase genes (GhDEADs) and categorizes them into four Clades (Clade A-D) through phylogenetic analysis. Promoter analysis reveals cis-acting elements linked to plant responses to light, methyl jasmonate (MeJA), abscisic acid (ABA), low temperature, and drought. RNA-seq data demonstrate that Clade C GhDEADs exhibit elevated and ubiquitous expression across different tissues, validating their connection to leaf development through real-time quantitative polymerase chain reaction (RT-qPCR) analysis. Notably, over half of GhDEADs display up-regulation in the leaves of virus-induced gene silencing (VIGS) plants of GhVIR-A/D (members of m6A methyltransferase complex, which regulate leaf morphogenesis). In conclusion, this study offers a comprehensive insight into GhDEADs, emphasizing their potential involvement in leaf development.

Cotton PRP5 gene encoding a proline-rich protein is involved in fiber development.

Proline-rich proteins contribute to cell wall structure of specific cell types and are involved in plant growth and development. In this study, a fiber-specific gene, GhPRP5, encoding a proline-rich protein was functionally characterized in cotton. GhPRP5 promoter directed GUS expression only in trichomes of both transgenic Arabidopsis and tobacco plants. The transgenic Arabidopsis plants with overexpressing GhPRP5 displayed reduced cell growth, resulting in smaller cell size and consequently plant dwarfs, in comparison with wild type plants. In contrast, knock-down of GhPRP5 expression by RNA interference in cotton enhanced fiber development. The fiber length of transgenic cotton plants was longer than that of wild type. In addition, some genes involved in fiber elongation and wall biosynthesis of cotton were up-regulated or down-regulated in the transgenic cotton plants owing to suppression of GhPRP5. Collectively, these data suggested that GhPRP5 protein as a negative regulator participates in modulating fiber development of cotton.

The Functions of an NAC Transcription Factor, GhNAC2-A06, in Cotton Response to Drought Stress.

Drought stress imposes severe constraints on crop growth and yield. The NAC transcription factors (TF) play a pivotal role in regulating plant stress responses. However, the biological functions and regulatory mechanisms of many cotton NACs have not been explored. In this study, we report the cloning and characterization of GhNAC2-A06, a gene encoding a typical cotton NAC TF. The expression of GhNAC2-A06 was induced by PEG treatment, drought stress, and ABA treatment. Furthermore, we investigated its function using the virus-induced gene silencing (VIGS) method. GhNAC2-A06 silenced plants exhibited a poorer growth status under drought stress conditions compared to the controls. The GhNAC2-A06 silenced cotton plants had a lower leaf relative water and chlorophyll content and a higher MDA content compared to the controls under the drought treatment. The levels of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) enzyme activity in the GhNAC2-A06 silenced plants were found to be lower compared to the controls when exposed to drought stress. Additionally, the downregulation of the drought stress-related genes, GhSAP12-D07, GhNCED1-A01, GhLEA14-A11, GhZAT10-D02, GhPROT2-A05, GhABF3-A03, GhABF2-D05, GhSAP3-D07, and GhCPK1-D04, was observed in the GhNAC2-A06 silenced cotton. Together, our research reveals that GhNAC2-A06 plays a role in the reaction of cotton to drought stress by affecting the expression of genes related to drought stress. The data obtained from this study lay the theoretical foundation for further in-depth research on the biological function and regulatory mechanisms of GhNAC2-A06.

Cotton BCP genes encoding putative blue copper-binding proteins are functionally expressed in fiber development and involved in response to high-salinity and heavy metal stresses.

Copper is vitally required for plants at low concentrations but extremely toxic for plants at elevated concentrations. Plants have evolved a series of mechanisms to prevent the consequences of the excess or deficit of copper. These mechanisms require copper-interacting proteins involved in copper trafficking. Blue copper-binding proteins (BCPs) are a class of copper proteins containing one blue copper-binding domain binding a single type I copper. To investigate the role of BCPs in plant development and in response to stresses, we isolated nine cDNAs encoding the putative blue copper-binding proteins (GhBCPs) from cotton (Gossypium hirsutum). Meanwhile, four corresponding genes (including GhBCP1-GhBCP4), which contain a single intron inserted in their conserved position, were isolated from cotton genome. Quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) analysis indicated that the nine GhBCP genes are differentially expressed in cotton tissues. Among them, GhBCP1 and GhBCP4 were predominantly expressed in fibers, while the transcripts of GhBCP2 and GhBCP3 were accumulated at relatively high levels in fibers. These four genes were strongly expressed in early fiber elongation, but dramatically declined with further fiber development. In addition, these GhBCP genes were upregulated in fibers by Cu(2+) , Zn(2+) , high-salinity and drought stresses, but downregulated in fibers by Al(3+) treatment. Overexpression of GhBCP1 and GhBCP4 in yeast (Schizosaccharomyces pombe) significantly increased the cell growth rate under Cu(2+) , Zn(2+) and high-salinity stresses. These results suggested that these GhBCPs may participate in the regulation of fiber development and in response to high-salinity and heavy metal stresses in cotton.

GhHyPRP4, a cotton gene encoding putative hybrid proline-rich protein, is preferentially expressed in leaves and involved in plant response to cold stress.

Plant hybrid proline-rich proteins (HyPRPs) usually consist of an N-terminal signal peptide, a central proline-rich domain, and a conserved eight-cysteine motif C-terminal domain. In this study, one gene (designated as GhHyPRP4) encoding putative HyPRP was isolated from cotton cDNA library. Northern blot and quantitative reverse transcriptase-polymerase chain reaction analyses revealed that GhHyPRP4 was preferentially expressed in leaves. Under cold stress, GhHyPRP4 expression was significantly up-regulated in leaves of cotton seedlings. Using the genome walking approach, a promoter fragment of GhHyPRP4 gene was isolated from cotton genome. GUS (ß-glucuronidase) gene driven by GhHyPRP4 promoter was specifically expressed in leaves and cotyledons of the transgenic Arabidopsis thaliana. Furthermore, GUS expression in leaves was remarkably induced by cold stress. Overexpression of GhHyPRP4 in yeast (Schizosaccharomyces pombe) significantly enhanced the cell survival rate upon treatment under -20°C for 60 h. These data suggested that GhHyPRP4 may be involved in plant response to cold stress during seedling development of cotton.

Family-Wide Evaluation of Multiple C2 Domain and Transmembrane Region Protein in Gossypium hirsutum.

Multiple C2 domain and transmembrane region proteins (MCTPs) are a group of evolutionarily conserved proteins and show emerging roles in mediating protein trafficking and signaling transduction. Although, several studies showed that MCTPs play important roles during plant growth and development, their biological functions in cotton remain largely unknown. Here, we identify and characterize 33 GhMCTP genes from upland cotton (Gossypium hirsutum) and reveal the diverse expression patterns of GhMCTPs in various tissues. We also find that GhMCTP7, GhMCTP12, and GhMCTP17 are highly expressed in the main stem apex, suggesting their possible roles in shoot development. Through analyzing different cotton species, we discover plant heights are closely related to the expression levels of GhMCTP7, GhMCTP12, and GhMCTP17. Furthermore, we silence the expression of GhMCTP genes using virus-induced gene silencing (VIGS) system in cotton and find that GhMCTP7, GhMCTP12, and GhMCTP17 play an essential role in shoot meristem development. GhMCTPs interact with GhKNAT1 and GhKNAT2 and regulate meristem development through integrating multiple signal pathways. Taken together, our results demonstrate functional redundancy of GhMCTPs in cotton shoot meristem development and provide a valuable resource to further study various functions of GhMCTPs in plant growth and development.

Three cotton genes preferentially expressed in flower tissues encode actin-depolymerizing factors which are involved in F-actin dynamics in cells.

To investigate whether the high expression levels of actin-depolymerizing factor genes are related to pollen development, three GhADF genes (cDNAs) were isolated and characterized in cotton. Among them, GhADF6 and GhADF8 were preferentially expressed in petals, whereas GhADF7 displayed the highest level of expression in anthers, revealing its anther specificity. The GhADF7 transcripts in anthers reached its peak value at flowering, suggesting that its expression is developmentally-regulated in anthers. The GhADF7 gene including the promoter region was isolated from the cotton genome. To demonstrate the specificity of the GhADF7 promoter, the 5'-flanking region, including the promoter and 5'-untranslated region, was fused with the GUS gene. Histochemical assays demonstrated that the GhADF7:GUS gene was specifically expressed in pollen grains. When pollen grains germinated, very strong GUS staining was detected in the elongating pollen tube. Furthermore, overexpression of GhADF7 gene in Arabidopsis thaliana reduced the viable pollen grains and, consequently, transgenic plants were partially male-sterile. Overexpression of GhADF7 in fission yeast (Schizosaccharomyces pombe) altered the balance of actin depolymerization and polymerization, leading to the defective cytokinesis and multinucleate formation in the cells. Given all the above results together, it is proposed that the GhADF7 gene may play an important role in pollen development and germination.

Expression and localization of GhH6L, a putative classical arabinogalactan protein in cotton (Gossypium hirsutum).

Arabinogalactan proteins (AGPs) are a large family of highly glycosylated of hydroxyproline-rich glycoproteins that play important roles in plant growth, development, and signal transduction. A cDNA encoding a putative classical AGP named GhH6L was isolated from cotton fiber cDNA libraries, and the deduced protein contains 17 copies of repetitive motif of X-Y-proline-proline-proline (where X is serine or alanine and Y is threonine or serine). Northern blotting analysis and quantitative RT-PCR results showed that it was preferentially expressed in 10 days post-anthesis (dpa) fibers and was also developmentally regulated. A promoter fragment was isolated from cotton (Gossypium hirsutum) by genome walking PCR. Expression of beta-glucuronidase (GUS) gene under the GhH6L promoter was examined in the transgenic Arabidopsis plants; only petiole and pedicel were stained, no staining was detected in other tissues. Subcellular localization indicated that GhH6L was localized to the plasma membrane and in the cytoplasm. These data further our understanding of GhH6L as well as shed light on functional insight to GhH6L in cotton.

Molecular Basis of Natural Variation in Photoperiodic Flowering Responses.

Plants exhibit different flowering behaviors in response to variable photoperiods across a wide geographical range. Here, we identify MYC3, a bHLH transcription factor, and its cis-element form the long-sought regulatory module responsible for cis-regulatory changes at the florigen gene FLOWERING LOCUS T (FT) that mediate natural variation in photoperiodic flowering responses in Arabidopsis. MYC3 is stabilized by DELLAs in the gibberellin pathway to suppress FT through binding the ACGGAT motif and antagonizing CONSTANS (CO) activation. Changing photoperiods modulate the relative abundance of MYC3 and CO, thus determining either of them as the predominant regulator for FT expression under different day lengths. Cis-regulatory changes in the MYC3 binding site at FT are associated with natural variation in day-length requirement for flowering in Arabidopsis accessions. Our findings reveal that environmental and developmental signals converge at MYC3 suppression of FT, an elementary event underlying natural variation in photoperiodic flowering responses.

GhEIN3, a cotton (Gossypium hirsutum) homologue of AtEIN3, is involved in regulation of plant salinity tolerance.

Ethylene insensitive 3 (EIN3), a key transcription factor in ethylene signal transduction, play important roles in plant stress signaling pathways. In this study, we isolated and characterized an EIN3-like gene from cotton (Gossypium hirsutum), designated as GhEIN3. GhEIN3 is highly expressed in vegetative tissues, and its expression is induced by 1-aminocyclopropane-1-carboxylic acid (ACC) and NaCl. Ectopic expression of GhEIN3 in Arabidopsis elevated plants' response to ethylene, which exhibit smaller leaves, more root hairs, shorter roots and hypocotyls. The germination rate, survival rate and root length of GhEIN3 transgenic plants were significantly improved compared to wild type under salt stress. GhEIN3 transgenic plants accumulated less H2O2 and malondialdehyde (MDA), while higher superoxide dismutase (SOD) and peroxidase (POD) activities were detected under salt stress. In addition, expression of several genes related to reactive oxygen species (ROS) pathway and ABA signaling pathway was increased in the GhEIN3 transgenic plants under salt stress. In contrast, virus-induced gene silencing (VIGS) of GhEIN3 in cotton enhanced the sensitivity of transgenic plants to salt stress, accumulating higher H2O2 and MDA and lower SOD and POD activities compared to control plants. Collectively, our results revealed that GhEIN3 might be involved in the regulation of plant response to salt stress by regulating ABA and ROS pathway during plant growth and development.

Characterization of 19 novel cotton FLA genes and their expression profiling in fiber development and in response to phytohormones and salt stress.

Fasciclin-like arabinogalactan proteins (FLAs), a subclass of arabinogalactan proteins (AGPs), are usually involved in cell development in plants. To investigate the expression profiling as well as the role of FLA genes in fiber development, 19 GhFLA genes (cDNAs) were isolated from cotton (Gossypium hirsutum). Among them, 15 are predicted to be glycosylphosphatidylinositol anchored to the plasma membranes. The isolated cotton FLAs could be divided into four groups. Real-time quantitative reverse transcriptase polymerase chain reaction results indicated that the GhFLA genes are differentially expressed in cotton tissues. Three genes (GhFLA1/2/4) were specifically or predominantly expressed in 10 days post-anthesis fibers, and the transcripts of the other four genes (GhFLA6/14/15/18) were accumulated at relatively high levels in cotton fibers. Furthermore, expressions of the GhFLA genes are regulated in fiber development and in response to phytohormones and NaCl. The identification of cotton FLAs will facilitate the study of their roles in cotton fiber development and cell wall biogenesis.

Genome-wide identification and characterization of JAZ gene family in upland cotton (Gossypium hirsutum).

Plant JAZ (Jasmonate ZIM-domain) proteins play versatile roles in multiple aspects of plant development and defense. However, little is known about the JAZ family in allotetraploid upland cotton (Gossypium hirsutum) so far. In this study, 30 non-redundant JAZ genes were identified in upland cotton through genome-wide screening. Phylogenetic analysis revealed that the 30 proteins in cotton JAZ family are further divided into five groups (I - V), and members in the same group share highly conserved motif structures. Subcellular localization assay demonstrated that GhJAZ proteins are localized in the cell nucleus. Quantitative RT-PCR analysis indicated that GhJAZs display different expression patterns in cotton tissues, and most of them could be induced by Jasmonic (JA). Furthermore, some GhJAZ genes are preferentially expressed in cotton ovules and fibers, and showed differential expression in ovules of wild type cotton and fiberless mutant (fl) during fiber initiation. GhJAZ proteins could interact with each other to form homodimer or heterodimer, and they also interacted with some JA signaling regulators and the proteins involved in cotton fiber initiation. Collectively, our data suggested that some GhJAZ proteins may play important roles in cotton fiber initiation and development by regulating JA signaling as well as some fiber-related proteins.

Genome-wide identification and characterization of TCP transcription factor genes in upland cotton (Gossypium hirsutum).

TCP proteins are plant-specific transcription factors (TFs), and perform a variety of physiological functions in plant growth and development. In this study, 74 non-redundant TCP genes were identified in upland cotton (Gossypium hirsutum L.) genome. Cotton TCP family can be classified into two classes (class I and class II) that can be further divided into 11 types (groups) based on their motif composition. Quantitative RT-PCR analysis indicated that GhTCPs display different expression patterns in cotton tissues. The majority of these genes are preferentially or specifically expressed in cotton leaves, while some GhTCP genes are highly expressed in initiating fibers and/or elongating fibers of cotton. Yeast two-hybrid results indicated that GhTCPs can interact with each other to form homodimers or heterodimers. In addition, GhTCP14a and GhTCP22 can interact with some transcription factors which are involved in fiber development. These results lay solid foundation for further study on the functions of TCP genes during cotton fiber development.