Identification and characterization of candidate genes for primary root length in Asiatic cotton (Gossypium arboreum L.).

KEY MESSAGE: One major gene controlling primary root length (PRL) in Gossypium arboreum is identified and this research provides a theoretical basis for root development for cotton. Primary root elongation is an essential process in plant root system structure. Here, we investigated the primary root length (PRL) of 215 diploid cotton (G. arboreum) accessions at 5, 8, 10, 15 days after sowing. A Genome-wide association study was performed for the PRL, resulting in 49 significant SNPs associated with 32 putative candidate genes. The SNP with the strongest signal (Chr07_8047530) could clearly distinguish the PRLs between accessions with two haplotypes. GamurG is the only gene that showed higher relative expression in the long PRL genotypes than the short PRL genotypes, which indicated it was the most likely candidate gene for regulating PRL. Moreover, the GamurG-silenced cotton seedlings showed a shorter PRL, while the GamurG-overexpressed Arabidopsis exhibited a significantly longer PRL. Our findings provide insight into the regulation mechanism of cotton root growth and will facilitate future breeding programs to optimize the root system structure in cotton.

Multistage temporal transcriptomic atlas unveils major contributor to reproductive phase in upland cotton.

Flowering is a major developmental transition in plants, but asynchronous flowering hinders the utilization of wild cotton relatives in breeding programs. We performed comparative transcriptomic profiling of early- and late-flowering Gossypium hirsutum genotypes to elucidate genetic factors influencing reproductive timing. Shoot apices were sampled from the photoperiod-sensitive landrace G. hirsutum purpurascens (GhP) and early-maturing variety ZhongMianSuo (ZMS) at five time points following the emergence of sympodial nodes. RNA-sequencing revealed extensive transcriptional differences during floral transition. Numerous flowering-associated genes exhibited genotype-specific expression, including FLOWERING LOCUS T (FT) homologs upregulated in ZMS. FT-interacting factors like SOC1 and CO-like also showed higher expression in ZMS, implicating florigen pathways in early flowering. Additionally, circadian clock and light signalling components were misregulated between varieties, suggesting altered photoperiod responses in GhP. Weighted co-expression network analysis specifically linked a module enriched for circadian-related genes to GhP's late flowering. Through an integrated transcriptome analysis, we defined a regulatory landscape of reproductive phase change in cotton. Differentially expressed genes related to photoperiod, circadian clock, and light signalling likely contribute to delayed flowering in wild cottons. Characterization of upstream flowering regulators will enable modifying photoperiod sensitivity and expand germplasm use for cotton improvement. This study provides candidate targets for elucidating interactive mechanisms that control cotton flowering time across diverse genotypes.

Silencing of GhORP_A02 enhances drought tolerance in Gossypium hirsutum.

BACKGROUND: ORP (Oxysterol-binding protein-related proteins) genes play a role in lipid metabolism, vesicular transferring and signaling, and non-vesicular sterol transport. However, no systematic identification and analysis of ORP genes have been reported in cotton. RESULT: In this study, we identified 14, 14, 7, and 7 ORP genes in G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. Phylogenetic analysis showed that all ORP genes could be classified into four groups. Gene structure and conserved motif analysis suggest that the function of this gene family was conserved. The Ka/Ks analysis showed that this gene family was exposed to purifying selection during evolution. Transcriptome data showed that four ORP genes, especially GhORP_A02, were induced by abiotic stress treatment. The cis-acting elements in the ORP promoters were responsive to phytohormones and various abiotic stresses. The silenced plants of GhORP_A02 were more sensitive to drought stress when compared to control. CONCLUSION: The major finding of this study shed light on the potential role of ORP genes in abiotic stress and provided a fundamental resource for further analysis in cotton.

Transcriptomic analysis reveals the beneficial effects of salt priming on enhancing defense responses in upland cotton under successive salt stress.

Priming-mediated stress tolerance in plants stimulates defense mechanisms and enables plants to cope with future stresses. Seed priming has been proven effective for tolerance against abiotic stresses; however, underlying genetic mechanisms are still unknown. We aimed to assess upland cotton genotypes and their transcriptional behaviors under salt priming and successive induced salt stress. We pre-selected 16 genotypes based on previous studies and performed morpho-physiological characterization, from which we selected three genotypes, representing different tolerance levels, for transcriptomic analysis. We subjected these genotypes to four different treatments: salt priming (P0), salt priming with salinity dose at 3-true-leaf stage (PD), salinity dose at 3-true-leaf stage without salt priming (0D), and control (CK). Although the three genotypes displayed distinct expression patterns, we identified common differentially expressed genes (DEGs) under PD enriched in pathways related to transferase activity, terpene synthase activity, lipid biosynthesis, and regulation of acquired resistance, indicating the beneficial role of salt priming in enhancing salt stress resistance. Moreover, the number of unique DEGs associated with G. hirsutum purpurascens was significantly higher compared to other genotypes. Coexpression network analysis identified 16 hub genes involved in cell wall biogenesis, glucan metabolic processes, and ribosomal RNA binding. Functional characterization of XTH6 (XYLOGLUCAN ENDOTRANSGLUCOSYLASE/HYDROLASE) using virus-induced gene silencing revealed that suppressing its expression improves plant growth under salt stress. Overall, findings provide insights into the regulation of candidate genes in response to salt stress and the beneficial effects of salt priming on enhancing defense responses in upland cotton.

Integrating Genome-wide association and whole transcriptome analysis to reveal genetic control of leaf traits in Gossypium arboreum L.

Leaves are important organs for crop photosynthesis and transpiration, and their morphological characteristics can directly reflect the growth state of plants. Accurate measurement of leaf traits and mining molecular markers are of great significance to the study of cotton growth. Here, we performed a Genome-wide association study on 7 leaf traits in 213 Asian cotton accessions. 32 significant SNPs and 44 genes were identified. A field experiment showed significant difference in leaf hair and leaf area between DPL971 and its natural mutant DPL972. We also compared the leaf transcriptome difference between DPL971 and DPL972, and found a batch of differentially expressed genes and non-coding RNAs (including lncRNAs, microRNAs, and circRNAs). After integrating the GWAS and transcriptome results, we finally selected two coding genes (Ga03G2383 and Ga05G3412) and two microRNAs (hbr-miR156, unconservative_Chr03_contig343_2364) as the candidate for leaf traits. Those findings will provide important genomic resources for cotton leaf improvement breeding.

Genome-wide association study for seedling biomass-related traits in Gossypium arboreum L.

BACKGROUND: Seedling stage plant biomass is usually used as an auxiliary trait to study plant growth and development or stress adversities. However, few molecular markers and candidate genes of seedling biomass-related traits were found in cotton. RESULT: Here, we collected 215 Gossypium arboreum accessions, and investigated 11 seedling biomass-related traits including the fresh weight, dry weight, water content, and root shoot ratio. A genome-wide association study (GWAS) utilizing 142,5003 high-quality SNPs identified 83 significant associations and 69 putative candidate genes. Furthermore, the transcriptome profile of the candidate genes emphasized higher expression of Ga03G1298, Ga09G2054, Ga10G1342, Ga11G0096, and Ga11G2490 in four representative cotton accessions. The relative expression levels of those five genes were further verified by qRT-PCR. CONCLUSIONS: The significant SNPs, candidate genes identified in this study are expected to lay a foundation for studying the molecular mechanism for early biomass development and related traits in Asian cotton.

Deltapine 15 contributes to the genomic architecture of modern upland cotton cultivars.

Foundation parents play a critical role in the genetic constituents of the derived genotypes. Deltapine-15 (DLP-15), introduced to China in 1950, is one of the most commonly used parents for early breeding programs in China. However, the formation and inheritance patterns of genomic constituents have not been studied. Therefore, this study aimed at understanding and exploring the genomic architecture of 146 DLP-15 derived cultivars with a common foundation parent DLP-15. Population structure based on sequencing data clustered genotypes into two groups (G1 and G2) supported by principal component analysis. Further exploration led to the identification of Chr-A08 with significantly differentiated regions between two groups. Moreover, we identified genome-wide identity by descent (IBD) segments (840 segments) to understand the genomic inheritance pattern in DLP-15 derived cultivars, spanning the 20-95 Mb region on Chr-A08. Interestingly, Chr-A08 depicted a unique inheritance pattern from DLP-15 to its derived cultivars. IBD-segment-based haplotype analysis suggested significant differences among the two groups. Phenotypic trait association with DLP-derived haplotypes concerning Chr-A08 suggested a significant increase in yield and fiber quality. Furthermore, distinguished IBD segments overlapped with previously reported QTLs concerning fiber yield and quality. Our results systematically identified genomic signatures transmitted from the foundation parent DLP-15 to its derived cultivars and provided a basis for further exploiting excellent haplotypes associated with DLP-15.

Genome-wide dissection of hybridization for fiber quality- and yield-related traits in upland cotton.

An evaluation of combining ability can facilitate the selection of suitable parents and superior F1 hybrids for hybrid cotton breeding, although the molecular genetic basis of combining ability has not been fully characterized. In the present study, 282 female parents were crossed with four male parents in accordance with the North Carolina II mating scheme to generate 1128 hybrids. The parental lines were genotyped based on restriction site-associated DNA sequencing and 306 814 filtered single nucleotide polymorphisms were used for genome-wide association analysis involving the phenotypes, general combining ability (GCA) values, and specific combining ability values of eight fiber quality- and yield-related traits. The main results were: (i) all parents could be clustered into five subgroups based on population structure analyses and the GCA performance of the female parents had significant differences between subgroups; (ii) 20 accessions with a top 5% GCA value for more than one trait were identified as elite parents for hybrid cotton breeding; (iii) 120 significant single nucleotide polymorphisms, clustered into 66 quantitative trait loci, such as the previously reported Gh_A07G1769 and GhHOX3 genes, were found to be significantly associated with GCA; and (iv) identified quantitative trait loci for GCA had a cumulative effect on GCA of the accessions. Overall, our results suggest that pyramiding the favorable loci for GCA may improve the efficiency of hybrid cotton breeding.

Assessment of heterosis based on parental genetic distance estimated with SSR and SNP markers in upland cotton (Gossypium hirsutum L.).

BACKGROUND: Heterosis has been extensively utilized in different crops and made a significant contribution to global food security. Genetic distance (GD) is one of the valuable criteria for selecting parents in hybrid breeding. The objectives of this study were to estimate the GD between parents using both simple sequence repeat (SSR) markers and single nucleotide polymorphism (SNP) markers and to investigate the efficiency of the prediction of hybrid performance based on GD. The experiment comprised of four male parents, 282 female parents and 1128 F1, derived from NCII mating scheme. The hybrids, their parents and two check cultivars were evaluated for two years. Performance of F1, mid-parent heterosis (MPH), and best parent heterosis (BPH) were evaluated for ten agronomic and fiber quality traits, including plant height, boll weight, boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index. RESULTS: Heterosis was observed in all hybrids and, the traits like plant height, boll number, boll weight and lint percentage exhibited higher heterosis than the fiber quality traits. Correlations were significant between parental and F1 performances. The F1 performances between three hybrid sets (Elite×Elite, Exotic×Elite, and Historic×Elite) showed significant differences in eight traits, including boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index. The correlation of the GD assessed by both SSR and SNP markers was significantly positive. The cluster analysis based on GD results estimated using SNP showed that all the female parents divided into five groups and the F1 performance between these five groups showed significant differences in four traits, including lint percentage, micronaire, fiber strength, and fiber elongation ratio. The correlation between GD and F1 performance, MPH and BPH were significant for lint percentage and micronaire. CONCLUSIONS: Our results suggested that GD between parents could be helpful in heterosis prediction for certain traits. This study reveals that molecular marker analysis can serve as a basis for assigning germplasm into heterotic groups and to provide guidelines for parental selection in hybrid cotton breeding.

GhGASA10-1 promotes the cell elongation in fiber development through the phytohormones IAA-induced.

BACKGROUND: Cotton is an important cash crop. The fiber length has always been a hot spot, but multi-factor control of fiber quality makes it complex to understand its genetic basis. Previous reports suggested that OsGASR9 promotes germination, width, and thickness by GAs in rice, while the overexpression of AtGASA10 leads to reduced silique length, which is likely to reduce cell wall expansion. Therefore, this study aimed to explore the function of GhGASA10 in cotton fibers development. RESULTS: To explore the molecular mechanisms underlying fiber elongation regulation concerning GhGASA10-1, we revealed an evolutionary basis, gene structure, and expression. Our results emphasized the conservative nature of GASA family with its origin in lower fern plants S. moellendorffii. GhGASA10-1 was localized in the cell membrane, which may synthesize and transport secreted proteins to the cell wall. Besides, GhGASA10-1 promoted seedling germination and root extension in transgenic Arabidopsis, indicating that GhGASA10-1 promotes cell elongation. Interestingly, GhGASA10-1 was upregulated by IAA at fiber elongation stages. CONCLUSION: We propose that GhGASA10-1 may promote fiber elongation by regulating the synthesis of cellulose induced by IAA, to lay the foundation for future research on the regulation networks of GASA10-1 in cotton fiber development.

Large-fragment insertion activates gene GaFZ (Ga08G0121) and is associated with the fuzz and trichome reduction in cotton (Gossypium arboreum).

Cotton seeds are typically covered by lint and fuzz fibres. Natural 'fuzzless' mutants are an ideal model system for identifying genes that regulate cell initiation and elongation. Here, using a genome-wide association study (GWAS), we identified a ~ 6.2 kb insertion, larINDELFZ , located at the end of chromosome 8, composed of a ~ 5.0 kb repetitive sequence and a ~ 1.2 kb fragment translocated from chromosome 12 in fuzzless Gossypium arboreum. The presence of larINDELFZ was associated with a fuzzless seed and reduced trichome phenotypes in G. arboreum. This distant insertion was predicted to be an enhancer, located ~ 18 kb upstream of the dominant-repressor GaFZ (Ga08G0121). Ectopic overexpression of GaFZ in Arabidopsis thaliana and G. hirsutum suggested that GaFZ negatively modulates fuzz and trichome development. Co-expression and interaction analyses demonstrated that GaFZ might impact fuzz fibre/trichome development by repressing the expression of genes in the very-long-chain fatty acid elongation pathway. Thus, we identified a novel regulator of fibre/trichome development while providing insights into the importance of noncoding sequences in cotton.

A copy number variant at the HPDA-D12 locus confers compact plant architecture in cotton.

Improving yield is a primary mission for cotton (Gossypium hirsutum) breeders; development of cultivars with suitable architecture for high planting density (HPDA) can increase yield per unit area. We characterized a natural cotton mutant, AiSheng98 (AS98), which exhibits shorter height, shorter branch length, and more acute branch angle than wild-type. A copy number variant at the HPDA locus on Chromosome D12 (HPDA-D12), encoding a dehydration-responsive element-binding (DREB) transcription factor, GhDREB1B, strongly affects plant architecture in the AS98 mutant. We found an association between a tandem duplication of a c. 13.5 kb segment in HPDA-D12 and elevated GhDREB1B expression resulting in the AS98 mutant phenotype. GhDREB1B overexpression confers a significant decrease in plant height and branch length, and reduced branch angle. Our results suggest that fine-tuning GhDREB1B expression may be a viable engineering strategy for modification of plant architecture favorable to high planting density in cotton.

Identification of Raf-Like Kinases B Subfamily Genes in Gossypium Species Revealed GhRAF42 Enhanced Salt Tolerance in Cotton.

Salinity is a critical abiotic factor that significantly reduces agricultural production. Cotton is an important fiber crop and a pioneer on saline soil, hence genetic architecture that underpins salt tolerance should be thoroughly investigated. The Raf-like kinase B-subfamily (RAF) genes were discovered to regulate the salt stress response in cotton plants. However, understanding the RAFs in cotton, such as Enhanced Disease Resistance 1 and Constitutive Triple Response 1 kinase, remains a mystery. This study obtained 29, 28, 56, and 54 RAF genes from G. arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively. The RAF gene family described allopolyploidy and hybridization events in allotetraploid cotton evolutionary connections. Ka/Ks analysis advocates that cotton evolution was subjected to an intense purifying selection of the RAF gene family. Interestingly, integrated analysis of synteny and gene collinearity suggested dispersed and segmental duplication events involved in the extension of RAFs in cotton. Transcriptome studies, functional validation, and virus-induced gene silencing on salt treatments revealed that GhRAF42 is engaged in salt tolerance in upland cotton. This research might lead to a better understanding of the role of RAFs in plants and the identification of suitable candidate salt-tolerant genes for cotton breeding.

Casparian strip membrane domain proteins in Gossypium arboreum: genome-wide identification and negative regulation of lateral root growth.

BACKGROUND: Root systems are critical for plant growth and development. The Casparian strip in root systems is involved in stress resistance and maintaining homeostasis. Casparian strip membrane domain proteins (CASPs) are responsible for the formation of Casparian strips. RESULTS: To investigate the function of CASPs in cotton, we identified and characterized 48, 54, 91 and 94 CASPs from Gossypium arboreum, Gossypium raimondii, Gossypium barbadense and Gossypium hirsutum, respectively, at the genome-wide level. However, only 29 common homologous CASP genes were detected in the four Gossypium species. A collinearity analysis revealed that whole genome duplication (WGD) was the primary reason for the expansion of the genes of the CASP family in the four cotton species. However, dispersed duplication could also contribute to the expansion of the GaCASPs gene family in the ancestors of G. arboreum. Phylogenetic analysis was used to cluster a total of 85 CASP genes from G. arboreum and Arabidopsis into six distinct groups, while the genetic structure and motifs of CASPs were conserved in the same group. Most GaCASPs were expressed in diverse tissues, with the exception of that five GaCASPs (Ga08G0113, Ga08G0114, Ga08G0116, Ga08G0117 and Ga08G0118) that were highly expressed in root tissues. Analyses of the tissue and subcellular localization suggested that GaCASP27 genes (Ga08G0117) are membrane protein genes located in the root. In the GaCASP27 silenced plants and the Arabidopsis mutants, the lateral root number significantly increased. Furthermore, GaMYB36, which is related to root development was found to regulate lateral root growth by targeting GaCASP27. CONCLUSIONS: This study provides a fundamental understanding of the CASP gene family in cotton and demonstrates the regulatory role of GaCASP27 on lateral root growth and development.

Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum).

KEY MESSAGE: Three extensive eco-haplotypes associated with population differentiation and environmental adaptability in Upland cotton were identified, with A06_85658585, A08_43734499 and A06_113104285 considered the eco-loci for environmental adaptability. Population divergence is suggested to be the primary force driving the evolution of environmental adaptability in various species. Chromosome inversion increases reproductive isolation between subspecies and accelerates population divergence to adapt to new environments. Although modern cultivated Upland cotton (Gossypium hirsutum L.) has spread worldwide, the noticeable phenotypic differences still existed among cultivars grown in different areas. In recent years, the long-distance migration of cotton cultivation areas throughout China has demanded that breeders better understand the genetic basis of environmental adaptability in Upland cotton. Here, we integrated the genotypes of 419 diverse accessions, long-term environment-associated variables (EAVs) and environment-associated traits (EATs) to evaluate subgroup differentiation and identify adaptive loci in Upland cotton. Two highly divergent genomic regions were found on chromosomes A06 and A08, which likely caused by extensive chromosome inversions. The subgroups could be geographically classified based on distinct haplotypes in the divergent regions. A genome-wide association study (GWAS) also confirmed that loci located in these regions were significantly associated with environmental adaptability in Upland cotton. Our study first revealed the cause of population divergence in Upland cotton, as well as the consequences of variation in its environmental adaptability. These findings provide new insights into the genetic basis of environmental adaptability in Upland cotton, which could accelerate the development of molecular markers for adaptation to climate change in future cotton breeding.

Identification and characterization analysis of sulfotransferases (SOTs) gene family in cotton (Gossypium) and its involvement in fiber development.

BACKGROUND: Sulfotransferases (SOTs) (EC 2.8.2.-) play a crucial role in the sulphate conjugation reaction involved in plant growth, vigor, stress resistance and pathogen infection. SOTs in Arabidopsis have been carried out and divided into 8 groups. However, the systematic analysis and functional information of SOT family genes in cotton have rarely been reported. RESULTS: According to the results of BLASTP and HMMER, we isolated 46, 46, 76 and 77 SOT genes in the genome G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. A total of 170 in 245 SOTs were further classified into four groups based on the orthologous relationships comparing with Arabidopsis, and tandem replication primarily contributed to the expansion of SOT gene family in G. hirsutum. Expression profiles of the GhSOT showed that most genes exhibited a high level of expression in the stem, leaf, and the initial stage of fiber development. The localization analysis indicated that GhSOT67 expressed in cytoplasm and located in stem and leaf tissue. Additionally, the expression of GhSOT67 were induced and the length of stem and leaf hairs were shortened after gene silencing mediated by Agrobacterium, compared with the blank and negative control plants. CONCLUSIONS: Our findings indicated that SOT genes might be associated with fiber development in cotton and provided valuable information for further studies of SOT genes in Gossypium.

Identifying favorable alleles for improving key agronomic traits in upland cotton.

BACKGROUND: Gossypium hirsutum L. is grown worldwide and is the largest source of natural fiber crop. We focus on exploring the favorable alleles (FAs) for upland cotton varieties improvement, and further understanding the history of accessions selection and acumination of favorable allele during breeding. RESULTS: The genetic basis of phenotypic variation has been studied. But the accumulation of favorable alleles in cotton breeding history in unknown, and potential favorable alleles to enhance key agronomic traits in the future cotton varieties have not yet been identified. Therefore, 419 upland cotton accessions were screened, representing a diversity of phenotypic variations of 7362 G. hirsutum, and 15 major traits were investigated in 6 environments. These accessions were categorized into 3 periods (early, medium, and modern) according to breeding history. All accessions were divided into two major groups using 299 polymorphic microsatellite markers: G1 (high fiber yield and quality, late maturity) and G2 (low fiber yield and quality, early maturity). The proportion of G1 genotype gradually increased from early to modern breeding periods. Furthermore, 21 markers (71 alleles) were significantly associated (-log P > 4) with 15 agronomic traits in multiple environments. Seventeen alleles were identified as FAs; these alleles accumulated more in the modern period than in other periods, consistent with their phenotypic variation trends in breeding history. Our results demonstrate that the favorable alleles accumulated through breeding effects, especially for common favorable alleles. However, the potential elite accessions could be rapidly screened by rare favorable alleles. CONCLUSION: In our study, genetic variation and genome-wide associations for 419 upland cotton accessions were analyzed. Two favorable allele types were identified during three breeding periods, providing important information for yield/quality improvement of upland cotton germplasm.

Association Analysis of Salt Tolerance in Asiatic cotton (Gossypium arboretum) with SNP Markers.

Salinity is not only a major environmental factor which limits plant growth and productivity, but it has also become a worldwide problem. However, little is known about the genetic basis underlying salt tolerance in cotton. This study was carried out to identify marker-trait association signals of seven salt-tolerance-related traits and one salt tolerance index using association analysis for 215 accessions of Asiatic cotton. According to a comprehensive index of salt tolerance (CIST), 215 accessions were mainly categorized into four groups, and 11 accessions with high salinity tolerance were selected for breeding. Genome-wide association studies (GWAS) revealed nine SNP rich regions significantly associated with relative fresh weight (RFW), relative stem length (RSL), relative water content (RRWC) and CIST. The nine SNP rich regions analysis revealed 143 polymorphisms that distributed 40 candidate genes and significantly associated with salt tolerance. Notably, two SNP rich regions on chromosome 7 were found to be significantly associated with two salinity related traits, RFW and RSL, by the threshold of -log10P ≥ 6.0, and two candidate genes (Cotton_A_37775 and Cotton_A_35901) related to two key SNPs (Ca7_33607751 and Ca7_77004962) were possibly associated with salt tolerance in G. arboreum. These can provide fundamental information which will be useful for future molecular breeding of cotton, in order to release novel salt tolerant cultivars.

Dissection of complicate genetic architecture and breeding perspective of cottonseed traits by genome-wide association study.

BACKGROUND: Cottonseed is one of the most important raw materials for plant protein, oil and alternative biofuel for diesel engines. Understanding the complex genetic basis of cottonseed traits is requisite for achieving efficient genetic improvement of the traits. However, it is not yet clear about their genetic architecture in genomic level. GWAS has been an effective way to explore genetic basis of quantitative traits in human and many crops. This study aims to dissect genetic mechanism seven cottonseed traits by a GWAS for genetic improvement. RESULTS: A genome-wide association study (GWAS) based on a full gene model with gene effects as fixed and gene-environment interaction as random, was conducted for protein, oil and 5 fatty acids using 316 accessions and ~ 390 K SNPs. Totally, 124 significant quantitative trait SNPs (QTSs), consisting of 16, 21, 87 for protein, oil and fatty acids (palmitic, linoleic, oleic, myristic, stearic), respectively, were identified and the broad-sense heritability was estimated from 71.62 to 93.43%; no QTS-environment interaction was detected for the protein, the palmitic and the oleic contents; the protein content was predominantly controlled by epistatic effects accounting for 65.18% of the total variation, but the oil content and the fatty acids except the palmitic were mainly determined by gene main effects and no epistasis was detected for the myristic and the stearic. Prediction of superior pure line and hybrid revealed the potential of the QTSs in the improvement of cottonseed traits, and the hybrid could achieve higher or lower genetic values compared with pure lines. CONCLUSIONS: This study revealed complex genetic architecture of seven cottonseed traits at whole genome-wide by mixed linear model approach; the identified genetic variants and estimated genetic component effects of gene, gene-gene and gene-environment interaction provide cotton geneticist or breeders new knowledge on the genetic mechanism of the traits and the potential molecular breeding design strategy.