Co-localization and analysis of miR477b with fiber length quantitative trait loci in cotton.

Cotton is an important cash crop for the textile industry. However, the understanding of natural genetic variation of fiber elongation in relation to miRNA is lacking. A miRNA gene (miR477b) was found to co-localize with a previously mapped fiber length (FL) quantitative trait locus (QTL). The miR477b was differentially expressed during fiber elongation between two backcross inbred lines (BILs) differing in FL and its precursor sequences. Bioinformatics and qRT-PCR analysis were further used to analyse the miRNA genes, which could produce mature miR477b. Cotton plants with virus-induced gene silencing (VIGS) constructs to over-express the allele of miR477b from the BIL with longer fibers had significantly longer fibers as compared with negative control plants, while the VIGS plants with suppressed miRNA expression had significantly shorter fibers. The expression level of the target gene (DELLA) and related genes (RDL1 and EXPA1 for DELLA through HOX3 protein) in the two BILs and/or the VIGS plants were generally congruent, as expected. This report represents one of the first comprehensive studies to integrate QTL linkage mapping and physical mapping of small RNAs with both small and mRNA transcriptome analysis, followed by VIGS, to identify candidate small RNA genes affecting the natural variation of fiber elongation in cotton.

Genomic and co-expression network analyses reveal candidate genes for oil accumulation based on an introgression population in Upland cotton (Gossypium hirsutum).

KEY MESSAGE: Integrated QTL mapping and WGCNA condense the potential gene regulatory network involved in oil accumulation. A glycosyl hydrolases gene (GhHSD1) for oil biosynthesis was confirmed in Arabidopsis, which will provide useful knowledge to understand the functional mechanism of oil biosynthesis in cotton. Cotton is an economical source of edible oil for the food industry. The genetic mechanism that regulates oil biosynthesis in cottonseeds is essential for the genetic enhancement of oil content (OC). To explore the functional genomics of OC, this study utilized an interspecific backcross inbred line population to dissect the quantitative trait locus (QTL) interlinked with OC. In total, nine OC QTLs were identified, four of which were novel, and each QTL explained 3.62-34.73% of the phenotypic variation of OC. The comprehensive transcript profiling of developing cottonseeds revealed 3,646 core genes differentially expressed in both inbred parents. Functional enrichment analysis determined 43 genes were annotated with oil biosynthesis processes. Implementation of weighted gene co-expression network analysis showed that 803 differential genes had a significant correlation with the OC phenotype. Further integrated analysis identified seven important genes located in OC QTLs. Of which, the GhHSD1 gene located in stable QTL qOC-Dt3-1 exhibited the highest functional linkages with the other network genes. Phylogenetic analysis showed significant evolutionary differences in the HSD1 sequences between oilseed- and starch- crops. Furthermore, the overexpression of GhHSD1 in Arabidopsis yielded almost 6.78% higher seed oil. This study not only uncovers important genetic loci for oil accumulation in cottonseed, but also provides a set of new candidate genes that potentially influence the oil biosynthesis pathway in cottonseed.

EB1C forms dimer and interacts with protein phosphatase 2A (PP2A) to regulate fiber elongation in upland cotton (Gossypium hirsutum).

Cotton is the most economically important natural fiber crop grown in more than sixty-five countries of the world. Fiber length is the main factor affecting fiber quality, but the existing main varieties are short in length and cannot suit the higher demands of the textile industry. It is necessary to discover functional genes that enable fiber length improvement in cotton through molecular breeding. In this study, overexpression of GhEB1C in Arabidopsis thaliana significantly promotes trichomes, tap roots, and root hairs elongation. The molecular regulation of GhEB1C involves its interactions with itself and GhB'ETA, and the function of GhEB1C regulation mainly depends on the two cysteine residues located at the C-terminal. In particular, the function activity of GhEB1C protein triggered with the regulation of protein phosphatase 2A, while silencing of GhEB1C in cotton significantly influenced the fiber protrusions and elongation mechanisms., Further, influenced the expression of MYB-bHLH-WD40 complex, brassinosteroids, and jasmonic acid-related genes, which showed that transcriptional regulation of GhEB1C is indispensable for cotton fiber formation and elongation processes. Our study analyzed the brief molecular mechanism of GhEB1C regulation. Further elucidated that GhEB1C can be a potential target gene to improve cotton fiber length through transgenic breeding.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1.

BACKGROUND: Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear. RESULTS: In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR). CONCLUSIONS: Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

Overexpression of cotton Trihelix transcription factor GhGT-3b_A04 enhances resistance to Verticillium dahliae and affects plant growth in Arabidopsis thaliana.

Verticillium wilt is a soil-borne fungal disease that severely affects cotton fiber yield and quality. Herein, a cotton Trihelix family gene, GhGT-3b_A04, was strongly induced by the fungal pathogen Verticillium dahliae. Overexpression of the gene in Arabidopsis thaliana enhanced the plant's resistance to Verticillium wilt but inhibited the growth of rosette leaves. In addition, the primary root length, root hair number, and root hair length increased in GhGT-3b_A04-overexpressing plants. The density and length of trichomes on the rosette leaves also increased. GhGT-3b_A04 localized to the nucleus, and transcriptome analysis revealed that it induced gene expression for salicylic acid synthesis and signal transduction and activated gene expression for disease resistance. The gene expression for auxin signal transduction and trichome development was reduced in GhGT-3b_A04-overexpressing plants. Our results highlight important regulatory genes for Verticillium wilt resistance and cotton fiber quality improvement. The identification of GhGT-3b_A04 and other important regulatory genes can provide crucial reference information for future research on transgenic cotton breeding.

Polarization multiplexed active mode-locking optoelectronic oscillator for frequency tunable dual-band microwave pulse signals generation.

We propose and experimentally demonstrate a polarization multiplexed active mode-locking optoelectronic oscillator (AML-OEO) based on a single dual-polarization binary phase-shift keying (DP-BPSK) modulator for frequency tunable dual-band microwave pulse signal generation. In order to realize mode-locking, two single-tone signals whose frequency are integer multiple of the free spectrum range (FSR) of AML-OEO are applied as active modulation signals (AMSs) at the bias ports of the DP-BPSK modulator. By dividing the AML-OEO into two loops with polarization demultiplexing, both the carrier frequency and pulse repetition frequency (PRF) of the dual-band microwave pulses are independently adjustable. In the experiment, microwave pulses with different PRFs of 162.4 kHz, 324.8 kHz and 812 kHz are generated based on fundamental, second-order harmonic and fifth-order harmonic mode-locking, respectively. In addition, the carrier frequency tunability within 4∼10 GHz is verified by inserting a frequency tunable electrical filter. The phase noise of the generated pulse signal at 10 kHz offset is better than -125 dBc/Hz.

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton.

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton.

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Linkage and association analyses reveal that hub genes in energy-flow and lipid biosynthesis pathways form a cluster in upland cotton.

Upland cotton is an important allotetraploid crop that provides both natural fiber for the textile industry and edible vegetable oil for the food or feed industry. To better understand the genetic mechanism that regulates the biosynthesis of storage oil in cottonseed, we identified the genes harbored in the major quantitative trait loci/nucleotides (QTLs/QTNs) of kernel oil content (KOC) in cottonseed via both multiple linkage analyses and genome-wide association studies (GWAS). In 'CCRI70' RILs, six stable QTLs were simultaneously identified by linkage analysis of CHIP and SLAF-seq strategies. In '0-153' RILs, eight stable QTLs were detected by consensus linkage analysis integrating multiple strategies. In the natural panel, thirteen and eight loci were associated across multiple environments with two algorithms of GWAS. Within the confidence interval of a major common QTL on chromosome 3, six genes were identified as participating in the interaction network highly correlated with cottonseed KOC. Further observations of gene differential expression showed that four of the genes, LtnD, PGK, LPLAT1, and PAH2, formed hub genes and two of them, FER and RAV1, formed the key genes in the interaction network. Sequence variations in the coding regions of LtnD, FER, PGK, LPLAT1, and PAH2 genes may support their regulatory effects on oil accumulation in mature cottonseed. Taken together, clustering of the hub genes in the lipid biosynthesis interaction network provides new insights to understanding the mechanism of fatty acid biosynthesis and TAG assembly and to further genetic improvement projects for the KOC in cottonseeds.

Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton.

Cotton is not only the most important fiber crop but also the fifth most important oilseed crop in the world because of its oil-rich seeds as a byproduct of fiber production. By comparative transcriptome analysis between two germplasms with diverse oil accumulation, we reveal pieces of the gene expression network involved in the process of oil synthesis in cottonseeds. Approximately, 197.16 Gb of raw data from 30 RNA sequencing samples with 3 biological replicates were generated. Comparison of the high-oil and low-oil transcriptomes enabled the identification of 7682 differentially expressed genes (DEGs). Based on gene expression profiles relevant to triacylglycerol (TAG) biosynthesis, we proposed that the Kennedy pathway (diacylglycerol acyltransferase-catalyzed diacylglycerol to TAG) is the main pathway for oil production, rather than the phospholipid diacylglycerol acyltransferase-mediated pathway. Using weighted gene co-expression network analysis, 5312 DEGs were obtained and classified into 14 co-expression modules, including the MEblack module containing 10 genes involved in lipid metabolism. Among the DEGs in the MEblack module, GhCYSD1 was identified as a potential key player in oil biosynthesis. The overexpression of GhCYSD1 in yeast resulted in increased oil content and altered fatty acid composition. This study may not only shed more light on the underlying molecular mechanism of oil accumulation in cottonseed oil, but also provide a set of new gene for potential enhancement of oil content in cottonseeds.

Genetics, Breeding and Genetic Engineering to Improve Cottonseed Oil and Protein: A Review.

Upland cotton (Gossypium hirsutum) is the world's leading fiber crop and one of the most important oilseed crops. Genetic improvement of cotton has primarily focused on fiber yield and quality. However, there is an increased interest and demand for enhanced cottonseed traits, including protein, oil, fatty acids, and amino acids for broad food, feed and biofuel applications. As a byproduct of cotton production, cottonseed is an important source of edible oil in many countries and could also be a vital source of protein for human consumption. The focus of cotton breeding on high yield and better fiber quality has substantially reduced the natural genetic variation available for effective cottonseed quality improvement within Upland cotton. However, genetic variation in cottonseed oil and protein content exists within the genus of Gossypium and cultivated cotton. A plethora of genes and quantitative trait loci (QTLs) (associated with cottonseed oil, fatty acids, protein and amino acids) have been identified, providing important information for genetic improvement of cottonseed quality. Genetic engineering in cotton through RNA interference and insertions of additional genes of other genetic sources, in addition to the more recent development of genome editing technology has achieved considerable progress in altering the relative levels of protein, oil, fatty acid profile, and amino acids composition in cottonseed for enhanced nutritional value and expanded industrial applications. The objective of this review is to summarize and discuss the cottonseed oil biosynthetic pathway and major genes involved, genetic basis of cottonseed oil and protein content, genetic engineering, genome editing through CRISPR/Cas9, and QTLs associated with quantity and quality enhancement of cottonseed oil and protein.

Quantitative Trait Locus Analysis and Identification of Candidate Genes Affecting Seed Size and Shape in an Interspecific Backcross Inbred Line Population of Gossypium hirsutum × Gossypium barbadense.

Seed size and shape are key agronomic traits affecting seedcotton yield and seed quality in cotton (Gossypium spp.). However, the genetic mechanisms that regulate the seed physical traits in cotton are largely unknown. In this study, an interspecific backcross inbred line (BIL) population of 250 BC1F7 lines, derived from the recurrent parent Upland CRI36 (Gossypium hirsutum) and Hai7124 (Gossypium barbadense), was used to investigate the genetic basis of cotton seed physical traits via quantitative trait locus (QTL) mapping and candidate gene identification. The BILs were tested in five environments, measuring eight seed size and shape-related traits, including 100-kernel weight, kernel length width and their ratio, kernel area, kernel girth, kernel diameter, and kernel roundness. Based on 7,709 single nucleotide polymorphic (SNP) markers, a total of 49 QTLs were detected and each explained 2.91-35.01% of the phenotypic variation, including nine stable QTLs mapped in at least three environments. Based on pathway enrichment, gene annotation, genome sequence, and expression analysis, five genes encoding starch synthase 4, transcription factor PIF7 and MYC4, ubiquitin-conjugating enzyme E27, and THO complex subunit 4A were identified as candidate genes that might be associated with seed size and shape. Our research provides valuable information to improve seed physical traits in cotton breeding.

Genome-Wide Analysis of the GW2-Like Genes in Gossypium and Functional Characterization of the Seed Size Effect of GhGW2-2D.

Cotton is one of the most economically important crops worldwide. Seed size is a vital trait for plants connected with yield and germination. GW2 encodes a RING_Ubox E3 ubiquitin ligase that controls seed development by affecting cell growth. Here, are few reports on GW2-like genes in cotton, and the function of GW2 in cotton is poorly understood. In the present study, a genome-wide analysis identified 6 and 3 GW2-like genes in each of the two cultivated tetraploids (Gossypium hirsutum and G. barbadense) and each of their diploid ancestral species (G. arboreum, G. raimondii), respectively. GhGW2-2D has the same functional domain and high sequence similarity with AtDA2 in Arabidopsis. Overexpression of GhGW2-2D in Arabidopsis significantly reduced seed and seedling size, suggesting GhGW2-2D is a potential target for regulating cotton seed size. These results provided information on the genetic and molecular basis of GW2-like genes in cotton, thus establishing a foundation for functional studies of cotton seeds.

Mapping of dynamic QTLs for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a backcross inbred line population of Upland cotton.

KEY MESSAGE: A backcross inbred line population of cotton was evaluated for Fusarium wilt race 4 resistance at different days after inoculation (DAI). Both constitutively expressed and developmentally regulated QTLs were detected. The soil-borne fungus Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) causes Fusarium wilt including seedling mortality in cotton. A backcross inbred line (BIL) population of 181 lines, derived from a bi-parental cross of moderately resistant non-recurrent Hai 7124 (Gossypium barbadense) and recurrent parent CCRI 36 (G. hirsutum), was evaluated under temperature-controlled conditions for FOV4 resistance with artificial inoculations. Based on three replicated tests evaluated at 7, 14, 21, and 28 days after inoculation (DAI), only 2-5 BILs showed lower disease severity ratings (DSR) than the parents while 22-50 BILs were more susceptible, indicating transgressive segregation toward susceptibility. Although DSR were overall congruent between DAI, there were many BILs displaying different responses to FOV4 across DAI. Genetic mapping using 7709 SNP markers identified 42 unique QTLs for four evaluation parameters- disease incidence (DI), DSR, mortality rate (MR), and area under disease progress curve (AUDPC), including 26 for two or more parameters. All five QTLs for AUDPC were co-localized with QTLs for DI, DSR, and/or MR at one or two DAI, indicating the unnecessary use of AUDPC in QTL mapping for FOV4 resistance. Those common QTLs explained the significant positive associations between parameters observed. Ten common QTLs with negative or positive additive effects were detected between DAI. DAI-specific and consistent QTLs were detected between DAI in cotton for the first time, suggesting the existence of both constitutively expressed and developmentally regulated QTLs for FOV4 resistance and the importance of evaluating genetic populations for FOV4 resistance at different growth stages.

Quantitative Trait Locus Analysis and Identification of Candidate Genes for Micronaire in an Interspecific Backcross Inbred Line Population of Gossypium hirsutum × Gossypium barbadense.

Cotton is the most important fiber crop and provides indispensable natural fibers for the textile industry. Micronaire (MIC) is determined by fiber fineness and maturity and is an important component of fiber quality. Gossypium barbadense L. possesses long, strong and fine fibers, while upland cotton (Gossypium hirsutum L.) is high yielding with high MIC and widely cultivated worldwide. To identify quantitative trait loci (QTLs) and candidate genes for MIC in G. barbadense, a population of 250 backcross inbred lines (BILs), developed from an interspecific cross of upland cotton CRI36 × Egyptian cotton (G. barbadense) Hai7124, was evaluated in 9 replicated field tests. Based on a high-density genetic map with 7709 genotyping-by-sequencing (GBS)-based single-nucleotide polymorphism (SNP) markers, 25 MIC QTLs were identified, including 12 previously described QTLs and 13 new QTLs. Importantly, two stable MIC QTLs (qMIC-D03-2 on D03 and qMIC-D08-1 on D08) were identified. Of a total of 338 genes identified within the two QTL regions, eight candidate genes with differential expression between TM-1 and Hai7124 were identified. Our research provides valuable information for improving MIC in cotton breeding.

Improving the accuracy of distance measurements based on beat frequency detection of a dual-frequency optoelectronic oscillator.

A high-accuracy distance measurements (DMs) approach based on beat frequency detection of a dual-frequency optoelectronic oscillator (OEO) is proposed and experimentally demonstrated. The dual-frequency OEO is formed with a single electro-optical modulator and a common length of energy storage fiber, and the beat frequency of the two oscillation signals can be directly achieved after a photodetector. Since the environmental disturbance has the same influence on the lengths of the two loops corresponding to the two oscillation frequencies, the environmental disturbance errors can be greatly reduced by beat frequency detection. In the experiment, we achieved a sensitivity of 49.9375 kHz/mm and a measurement error of ±15µm. The frequency stability of the beat-frequency signal was 8.57 times higher than the oscillation signal of the measurement loop.

Transcriptome analysis reveals genes potentially related to high fiber strength in a Gossypium hirsutum line IL9 with Gossypium mustelinum introgression.

Cotton (Gossypium L.) is the most important fiber crop worldwide. Here, transcriptome analysis was conducted on developing fibers of a G. mustelinum introgression line, IL9, and its recurrent parent, PD94042, at 17 and 21 days post-anthesis (dpa). Differentially expressed genes (DEGs) of PD94042 and IL9 were identified. Gene Ontology (GO) enrichment analysis showed that the annotated DEGs were rich in two main biological processes and two main molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis likewise showed that the annotated DEGs were mainly enriched in metabolic pathways and biosynthesis of secondary metabolites. In total, 52 DEGs were selected as candidate genes based on comparison of the DEGs and GO function annotation information. Quantitative real-time PCR (RT-qPCR) analysis results for 12 randomly selected DEGs were consistent with transcriptome analysis. SNP identification based on G. mustelinum chromatin segment introgression showed that 394 SNPs were identified in 268 DEGs, and two genes with known functions were identified within fiber strength quantitative trait loci (QTL) regions or near the confidence intervals. We identified 52 key genes potentially related to high fiber strength in a G. mustelinum introgression line and provided significant insights into the study of cotton fiber quality improvement.

Genome-Wide Analysis of DA1-Like Genes in Gossypium and Functional Characterization of GhDA1-1A Controlling Seed Size.

Cotton (Gossypium spp.) is an economically important crop grown for natural fiber and seed oil production. DA1 is a ubiquitin receptor that determines final seed and organ size by restricting the period of cell proliferation. In the present study, we identified 7 DA1-like genes each in cultivated tetraploid (AADD) G. hirsutum and G. barbadense, and 4 and 3 DA1-like genes in their ancestral diploid G. arboreum (A2A2) and G. raimondii (D5D5), respectively. The 7 GhDA1 genes were confirmed to be distributed on four At and three Dt subgenome chromosomes in G. hirsutum. GhDA1-1A showed a high sequence similarity to AtDA1 in Arabidopsis, and they possessed the same functional domains, suggesting conserved functions. The overexpression of GhDA1-1A R301K in Arabidopsis significantly increased seed size and seed weight, indicating that GhDA1-1A is a promising target for cotton improvement. This study provides information on the molecular evolutionary properties of DA1-like genes in cotton, which will be useful for the genetic improvement of cotton.

The cellulose synthase (CesA) gene family in four Gossypium species: phylogenetics, sequence variation and gene expression in relation to fiber quality in Upland cotton.

Cellulose synthases (CesAs) are multi-subunit enzymes found on the plasma membrane of plant cells and play a pivotal role in cellulose production. The cotton fiber is mainly composed of cellulose, and the genetic relationships between CesA genes and cotton fiber yield and quality are not fully understood. Through a phylogenetic analysis, the CesA gene family in diploid Gossypium arboreum and Gossypium raimondii, as well as tetraploid Gossypium hirsutum ('TM-1') and Gossypium barbadense ('Hai-7124' and '3-79'), was divided into 6 groups and 15 sub-groups, with each group containing two to five homologous genes. Most CesA genes in the four species are highly collinear. Among the five cotton genomes, 440 and 1929 single nucleotide polymorphisms (SNPs) in the CesA gene family were identified in exons and introns, respectively, including 174 SNPs resulting in amino acid changes. In total, 484 homeologous SNPs between the A and D genomes were identified in diploids, while 142 SNPs were detected between the two tetraploids, with 32 and 82 SNPs existing within G. hirsutum and G. barbadense, respectively. Additionally, 74 quantitative trait loci near 18 GhCesA genes were associated with fiber quality. One to four GhCesA genes were differentially expressed (DE) in ovules at 0 and 3 days post anthesis (DPA) between two backcross inbred lines having different fiber lengths, but no DE genes were identified between these lines in developing fibers at 10 DPA. Twenty-seven SNPs in above DE CesA genes were detected among seven cotton lines, including one SNP in Ghi_A08G03061 that was detected in four G. hirsutum genotypes. This study provides the first comprehensive characterization of the cotton CesA gene family, which may play important roles in determining cotton fiber quality.