Search

Abstract Background Focusing on key indicators of drought resistance is highly important for quickly mining candidate genes related to drought resistance in cotton. Results In the present study, drought resistance was identified in drought resistance-related RIL populations during the flowering and boll stages, and multiple traits were evaluated; these traits included three key indicators: plant height (PH), single boll weight (SBW) and transpiration rate (Tr). Based on these three key indicators, three groups of extreme mixing pools were constructed for BSA-seq. Based on the mapping interval of each trait, a total of 6.27 Mb QTL intervals were selected on chromosomes A13 (3.2 Mb), A10 (2.45 Mb) and A07 (0.62 Mb) as the focus of this study. Based on the annotation information and qRT‒PCR analysis, three key genes that may be involved in the drought stress response of cotton were screened: GhF6'H1, Gh3AT1 and GhPER55. qRT‒PCR analysis of parental and extreme germplasm materials revealed that the expression of these genes changed significantly under drought stress. Cotton VIGS experiments verified the important impact of key genes on cotton drought resistance. Conclusions This study focused on the key indicators of drought resistance, laying the foundation for the rapid mining of drought-resistant candidate genes in cotton and providing genetic resources for directed molecular breeding of drought resistance in cotton.

Cotton seeds, as the main by-product of cotton, are not only an important raw material for edible oil and feed but also a source of biofuel. The quality of cotton seeds directly affects cotton planting and is closely related to the yield and fiber quality. However, the molecular mechanism governing cotton seed size remains largely unexplored. This study investigates the regulatory mechanisms of cotton seed size by focusing on two cotton genotypes, N10 and N12, which exhibit notable phenotypic variations across multiple environments. Developing seeds were sampled at various stages (5, 20, 30, and 35 DPA) and subjected to RNA-seq. Temporal pattern clustering and WGCNA on differentially expressed genes identified 413 candidate genes, including these related to sugar metabolism that were significantly enriched in transcriptional regulation. A genetic transformation experiment indicated that the overexpression of the GhUXS5 gene encoding UDP-glucuronate decarboxylase 5 significantly increased seed size, suggesting an important role of GhUXS5 in regulating cotton seed size. This discovery provides crucial insights into the molecular mechanisms controlling cotton seed size, helping to unravel the complex regulatory network and offering new strategies and targets for cotton breeding to enhance the economic value of cotton seeds and overall cotton yield.

Drought stress significantly affects the growth, development, and yield of cotton, triggering the response of multiple genes. Among them, ascorbate peroxidase (APX) is one of the important antioxidant enzymes in the metabolism of reactive oxygen species in plants, and APX enhances the ability of plants to resist oxidation, thus increasing plant stress tolerance. Therefore, enhancing the activity of APX in cells is crucial to improving plant stress resistance. Previous studies have isolated differentially expressed proteins under drought stress (GhAPX7) in drought-resistant (KK1543) and drought-sensitive (XLZ26) plants. Thus, this study analyzed the expression patterns of GhAPX7 in different cotton tissues to verify the drought resistance function of GhAPX7 and explore its regulatory pathways. GhAPX7 had the highest expression in cotton leaves, which significantly increased under drought stress, suggesting that GhAPX7 is essential for improving antioxidant capacity and enzyme activities in cotton. GhAPX7 silencing indirectly affects pronounced leaf yellowing and wilting in drought-resistant and drought-sensitive plants under drought stress. Malondialdehyde (MDA) content was significantly increased and chlorophyll and proline content and APX enzyme activity were generally decreased in silenced plants compared to the control. This result indicates that GhAPX7 may improve drought resistance by influencing the contents of MDA, chlorophyll, proline, and APX enzyme activity through increased expression levels. Transcriptome analysis revealed that the drought-related differentially expressed genes between the control and treated groups enriched plant hormone signal transduction, MAPK signaling, and plant–pathogen interaction pathways. Therefore, the decreased expression of GhAPX7 significantly affects the expression levels of genes in these three pathways, reducing drought resistance in plants. This study provides insights into the molecular mechanisms of GhAPX7 and its role in drought resistance and lays a foundation for further research on the molecular mechanisms of response to drought stress in cotton.

Abstract Background Flavonols play important roles in antioxidation and anticancer activities, longevity, and cardiovascular protection. Flavonol synthase (FLS) is a key enzyme for flavonol synthesis. Result Phenotypic, transcriptional and metabolic data were analyzed, which showed that there was a close relationship between salt stress and flavonoids, and flavonols were significantly upregulated under salt stress. Nine, seven, four, and four FLS genes were identified in Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum, and Gossypium raimondii, respectively. The results of subcellular localization showed that FLS existed in the nucleus and cytoplasmic. Through phylogenetic analysis, 24 FLS genes were divided into three subfamilies. The results of the RNA sequencing showed that the expression of GhFLS genes was mainly induced by salt, drought, low temperature, and heat stress. GhFLS promoter mainly comprised plant hormone response elements and abiotic stress elements, indicating that the GhFLS gene may play a key role in abiotic stress response. The proline contents of pYL156:GhFLS1 was reduced significantly compared to pYL156 under salt stress, thereby reducing the resistance of cotton to salt stress. Conclusion This study will lay a foundation for further study on the antioxidant regulation mechanism of the FLS gene under abiotic stress.

A-galactosidases (AGALs), the oligosaccharide (RFO) catabolic genes of the raffinose family, play crucial roles in plant growth and development and in adversity stress. They can break down the non-reducing terminal galactose residues of glycolipids and sugar chains. In this study, the whole genome of AGALs was analyzed. Bioinformatics analysis was conducted to analyze members of the AGAL family in Gossypium hirsutum, Gossypium arboreum, Gossypium barbadense, and Gossypium raimondii. Meanwhile, RT-qPCR was carried out to analyze the expression patterns of AGAL family members in different tissues of terrestrial cotton. It was found that a series of environmental factors stimulated the expression of the GhAGAL3 gene. The function of GhAGAL3 was verified through virus-induced gene silencing (VIGS). As a result, GhAGAL3 gene silencing resulted in milder wilting of seedlings than the controls, and a significant increase in the raffinose content in cotton, indicating that GhAGAL3 responded to NaCl stress. The increase in raffinose content improved the tolerance of cotton. Findings in this study lay an important foundation for further research on the role of the GhAGAL3 gene family in the molecular mechanism of abiotic stress resistance in cotton.

Background Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is the central enzyme of glycolysis and plays important regulatory roles in plant growth and development and responses to adverse stress conditions. However, studies on the characteristics and functions of cotton GAPDH family genes are still lacking. Methods In this study, genome-wide identification of the cotton GAPDH gene family was performed, and the phylogeny, gene structures, promoter progenitors and expression profiles of upland cotton GAPDH gene family members were explored by bioinformatics analysis to highlight potential functions. The functions of GhGAPDH9 in response to drought stress were initially validated based on RNA-seq, qRT‒PCR, VIGS techniques and overexpression laying a foundation for further studies on the functions of GAPDH genes. Results This study is the first systematic analysis of the cotton GAPDH gene family, which contains a total of 84 GAPDH genes, among which upland cotton contains 27 members. Quantitative, phylogenetic and covariance analyses of the genes revealed that the GAPDH gene family has been conserved during the evolution of cotton. Promoter analysis revealed that most cis-acting elements were related to MeJA and ABA. Based on the identified promoter cis-acting elements and RNA-seq data, it was hypothesized that Gh_GAPDH9, Gh_GAPDH11, Gh_GAPDH19 and Gh_GAPDH21 are involved in the response of cotton to abiotic stress. The expression levels of the Gh_GAPDH9 gene in two drought-resistant and two drought-sensitive materials were analyzed by qRT‒PCR and found to be high early in the treatment period in the drought-resistant material. The silencing of Gh_GAPDH9 based on virus-induced gene silencing (VIGS) technology resulted in significant leaf wilting or whole-plant dieback in silenced plants after drought stress compared to the control. The content of—malondialdehyde (MDA) in cotton leaves was significantly increased, and the content of proline (Pro) and chlorophyll (Chl) was reduced. In addition, the leaf wilting and dryness of transgenic lines under drought stress were lower than those of wild-type Arabidopsis, indicating that Gh_GAPDH9 is a positive regulator of drought resistance. In conclusion, our results demonstrate that GAPDH genes play an important role in the response of cotton to abiotic stresses and provide preliminary validation of the function of the Gh_GAPDH9 gene under drought stress. These findings provide an important theoretical basis for further studies on the function of the Gh_GAPDH9 gene and the molecular mechanism of the drought response in cotton.

DUSPs, a diverse group of protein phosphatases, play a pivotal role in orchestrating cellular growth and development through intricate signaling pathways. Notably, they actively participate in the MAPK pathway, which governs crucial aspects of plant physiology, including growth regulation, disease resistance, pest resistance, and stress response. DUSP is a key enzyme, and it is the enzyme that limits the rate of cell metabolism. At present, complete understanding of the DUSP gene family in cotton and its specific roles in resistance to Verticillium wilt (VW) remains elusive. To address this knowledge gap, we conducted a comprehensive identification and analysis of four key cotton species: Gossypium arboreum, Gossypium barbadense, Gossypium hirsutum, and Gossypium raimondii. The results revealed the identification of a total of 120 DUSP genes in the four cotton varieties, which were categorized into six subgroups and randomly distributed at both ends of 26 chromosomes, predominantly localized within the nucleus. Our analysis demonstrated that closely related DUSP genes exhibited similarities in terms of the conserved motif composition and gene structure. A promoter analysis performed on the GhDUSP gene promoter revealed the presence of several cis-acting elements, which are associated with abiotic and biotic stress responses, as well as hormone signaling. A tissue expression pattern analysis demonstrated significant variations in GhDUSP gene expression under different stress conditions, with roots exhibiting the highest levels, followed by stems and leaves. In terms of tissue-specific detection, petals, leaves, stems, stamens, and receptacles exhibited higher expression levels of the GhDUSP gene. The gene expression analysis results for GhDUSPs under stress suggest that DUSP genes may have a crucial role in the cotton response to stress in cotton. Through Virus-Induced Gene Silencing (VIGS) experiments, the silencing of the target gene significantly reduced the resistance efficiency of disease-resistant varieties against Verticillium wilt (VW). Consequently, we conclude that GH_A11G3500-mediated bispecific phosphorylated genes may serve as key regulators in the resistance of G. hirsutum to Verticillium wilt (VW). This study presents a comprehensive structure designed to provide an in-depth understanding of the potential biological functions of cotton, providing a strong foundation for further research into molecular breeding and resistance to plant pathogens.

Cotton is a critical crop with massive economic implications worldwide. Verticillium wilt is a soil-borne ailment caused by Verticillium dahliae, which harms the growth and development of cotton. Therefore, investigating the genes associated with resistance to verticillium wilt is of particular significance. In this study, we identified the GhIQD1 gene through transcriptome analysis and experimentally characterized the role of the GhIQD1 gene in cotton against V. dahliae. The findings indicated that GhIQD1 acts as a calmodulin-binding protein. The expression of GhIQD1 was the highest in stems, and the expression level increased significantly following infection with V. dahliae. The expression in resistant cotton varieties was higher than in susceptible cotton varieties. Through overexpression of the GhIQD1 gene in tobacco, these transgenic plants exhibited improved resistance to V. dahliae. In contrast, by silencing the GhIQD1 gene in cotton through VIGS, the resistance to V. dahliae was reduced. Following inoculation, the leaves yellowed, and the disease index was higher. Transcriptome analysis of transgenic tobacco 72 h after inoculation indicated that overexpression of GhIQD1 increased the enrichment of the calmodulin pathway and stimulated the production of plant hormones alongside secondary metabolites. Consequently, we investigated the relationship between the GhIQD1 gene and plant disease-resistant hormones SA, JA, and ABA. In summary, this study uncovered the mechanism by which GhIQD1 conferred resistance to V. dahliae in cotton through positive regulation of JA and ABA, providing crucial information for further research on the adaptation of plants to pathogen invasion.

Abstract Background Abscisic acid (ABA) is an important stress hormone, the changes of abscisic acid content can alter plant tolerance to stress, abscisic acid is crucial for studying plant responses to abiotic stress. The abscisic acid aldehyde oxidase (AAO) plays a vital role in the final step in the synthesis of abscisic acid, therefore, understanding the function of AAO gene family is of great significance for plants to response to abiotic stresses. Result In this study, 6, 8, 4 and 4 AAO genes were identified in four cotton species. According to the structural characteristics of genes and the traits of phylogenetic tree, we divided the AAO gene family into 4 clades. Gene structure analysis showed that the AAO gene family was relatively conservative. The analysis of cis-elements showed that most AAO genes contained cis-elements related to light response and plant hormones. Tissue specificity analysis under NaHCO3 stress showed that GhAAO2 gene was differentially expressed in both roots and leaves. After GhAAO2 gene silencing, the degree of wilting of seedlings was lighter than that of the control group, indicating that GhAAO2 could respond to NaHCO3 stress. Conclusions In this study, the AAO gene family was analyzed by bioinformatics, the response of GhAAO gene to various abiotic stresses was preliminarily verified, and the function of the specifically expressed gene GhAAO2 was further verified. These findings provide valuable information for the study of potential candidate genes related to plant growth and stress.

Abstract Background Serine carboxypeptidase-like protein (SCPL) plays an important role in response to stress in plant. However, our knowledge of the function of the SCPL gene family is limited. Results In this study, a comprehensive and systematic analysis of SCPL gene family was conducted to explore the phylogeny and evolution of the SCPL gene in Gossypium hirsutum. The phenotype and molecular mechanism of silencing of the Gh_SCPL42 under Verticillium wilt stress was also studied. Our results showed that 96 SCPL genes were observed in genome of G. hirsutum, which distributed on 25 chromosomes and most of them were located in the nucleus. The phylogenetic tree analysis showed that members of SCPL gene family can be divided into three subgroups in G. hirsutum, which are relatively conservative in evolution. SCPL gene has a wide range of tissue expression types in G. hirsutum. Promoter analysis showed that the most cis-acting elements related to MeJA and ABA were contained. Through RNA-seq combined with genotyping, it was found that 11 GhSCPL genes not only had significant expression changes during Verticillium wilt stress but also had differential SNPs in the upstream, downstream, exonic or intronic regions. The expression of these 11 genes in the resistant (Zhongzhimian 2) and susceptible (Junmian 1) materials was further analyzed by qRT-PCR, it was found that 6 genes showed significant expression differences in the two materials. Among them, Gh_SCPL42 has the most obvious expression change. Furthermore, virus-induced gene silencing (VIGS) showed necrosis and yellowing of leaves and significantly higher disease severity index (DSI) and disease severity rate (DSR) values in VIGS plants than in control silenced Gh_SCPL42 plants. Moreover, the expression levels of genes related to the SA and JA pathways were significantly downregulated. These results show that Gh_SCPL42 might improve resistance to Verticillium wilt through the SA and JA pathways in G. hirsutum. Conclusion In conclusion, our findings indicated that Gh_SCPL42 gene plays an important role in resistance to Verticillium wilt in cotton. It was provided an important theoretical basis for further research on the function of SCPL gene family and the molecular mechanism of resistance to Verticillium wilt in cotton.

DVL is one of the small polypeptides which plays an important role in regulating plant growth and development, tissue differentiation, and organ formation in the process of coping with stress conditions. So far, there has been no comprehensive analysis of the expression profile and function of the cotton DVL gene. According to previous studies, a candidate gene related to the development of fuzz was screened, belonging to the DVL family, and was related to the development of trichomes in Arabidopsis thaliana. However, the comprehensive identification and systematic analysis of DVL in cotton have not been conducted. In this study, we employed bioinformatics approaches to conduct a novel analysis of the structural characteristics, phylogenetic tree, gene structure, expression pattern, evolutionary relationship, and selective pressure of the DVL gene family members in four cotton species. A total of 117 DVL genes were identified, including 39 members in G. hirsutum. Based on the phylogenetic analysis, the DVL protein sequences were categorized into five distinct subfamilies. Additionally, we successfully mapped these genes onto chromosomes and visually represented their gene structure information. Furthermore, we predicted the presence of cis-acting elements in DVL genes in G. hirsutum and characterized the repeat types of DVL genes in the four cotton species. Moreover, we computed the Ka/Ks ratio of homologous genes across the four cotton species and elucidated the selective pressure acting on these homologous genes. In addition, we described the expression patterns of the DVL gene family using RNA-seq data, verified the correlation between GhMDVL3 and fuzz development through VIGS technology, and found that some DVL genes may be involved in resistance to biotic and abiotic stress conditions through qRT-PCR technology. Furthermore, a potential interaction network was constructed by WGCNA, and our findings demonstrated the potential of GhM_A05G1032 to interact with numerous genes, thereby playing a crucial role in regulating fuzz development. This research significantly contributed to the comprehension of DVL genes in upland cotton, thereby establishing a solid basis for future investigations into the functional aspects of DVL genes in cotton.

IntroductionStarch metabolism is involved in the stress response. Starch synthase (SS) is the key enzyme in plant starch synthesis, which plays an indispensable role in the conversion of pyrophosphoric acid to starch. However, the SS gene family in cotton has not been comprehensively identified and systematically analyzed.ResultIn our study, a total of 76 SS genes were identified from four cotton genomes and divided into five subfamilies through phylogenetic analysis. Genetic structure analysis proved that SS genes from the same subfamily had similar genetic structure and conserved sequences. A cis-element analysis of the SS gene promoter showed that it mainly contains light response elements, plant hormone response elements, and abiotic stress elements, which indicated that the SS gene played key roles not only in starch synthesis but also in abiotic stress response. Furthermore, we also conducted a gene interaction network for SS proteins. Silencing GhSS9 expression decreased the resistance of cotton to drought stress. These findings suggested that SS genes could be related to drought stress in cotton, which provided theoretical support for further research on the regulation mechanism of SS genes on abiotic starch synthesis and sugar levels.

Nitrate transporters (NRTs) are crucial for the uptake, use, and storage of nitrogen by plants. In this study, 42 members of the GhNRT2 (Nitrate Transporter 2 family) were found in the four different cotton species. The conserved domains, phylogenetic relationships, physicochemical properties, subcellular localization, conserved motifs, gene structure, cis-acting elements, and promoter region expression patterns of these 42 members were analyzed. The findings confirmed that members of the NRT2 family behaved typically, and subcellular localization tests confirmed that they were hydrophobic proteins that were mostly located on the cytoplasmic membrane. The NRT2 family of genes with A.thaliana and rice underwent phylogenetic analysis, and the results revealed that GhNRT2 could be divided into three groups. The same taxa also shared similar gene structure and motif distribution. The composition of cis-acting elements suggests that most of the expression of GhNRT2 may be related to plant hormones, abiotic stress, and photoreactions. The GhNRT2 gene was highly expressed, mainly in roots. Drought, salt, and extreme temperature stress showed that GhNRT2 gene expression was significantly up-regulated or down-regulated, indicating that it may be involved in the stress response of cotton. In general, the genes of the NRT2 family of cotton were comprehensively analyzed, and their potential nitrogen uptake and utilization functions in cotton were preliminarily predicted. Additionally, we provide an experimental basis for the adverse stress conditions in which they may function.

Cotton is an important fiber crop. The cotton fiber is an extremely long trichome that develops from the epidermis of an ovule. The trichome is a general and multi-function plant organ, and trichome birefringence-like (TBL) genes are related to trichome development. At the genome-wide scale, we identified TBLs in four cotton species, comprising two cultivated tetraploids (Gossypium hirsutum and G. barbadense) and two ancestral diploids (G. arboreum and G. raimondii). Phylogenetic analysis showed that the TBL genes clustered into six groups. We focused on GH_D02G1759 in group IV because it was located in a lint percentage-related quantitative trait locus. In addition, we used transcriptome profiling to characterize the role of TBLs in group IV in fiber development. The overexpression of GH_D02G1759 in Arabidopsis thaliana resulted in more trichomes on the stems, thereby confirming its function in fiber development. Moreover, the potential interaction network was constructed based on the co-expression network, and it was found that GH_D02G1759 may interact with several genes to regulate fiber development. These findings expand our knowledge of TBL family members and provide new insights for cotton molecular breeding.

Fiber quality traits, especially fiber strength, length, and micronaire (FS, FL, and FM), have been recognized as critical fiber attributes in the textile industry, while the lint percentage (LP) was an important indicator to evaluate the cotton lint yield. So far, the genetic mechanism behind the formation of these traits is still unclear. Quantitative trait loci (QTL) identification and candidate gene validation provide an effective methodology to uncover the genetic and molecular basis of FL, FS, FM, and LP. A previous study identified three important QTL/QTL cluster loci, harboring at least one of the above traits on chromosomes A01, A07, and D12 via a recombinant inbred line (RIL) population derived from a cross of Lumianyan28 (L28) × Xinluzao24 (X24). A secondary segregating population (F2) was developed from a cross between L28 and an RIL, RIL40 (L28 × RIL40). Based on the population, genetic linkage maps of the previous QTL cluster intervals on A01 (6.70–10.15 Mb), A07 (85.48–93.43 Mb), and D12 (0.40–1.43 Mb) were constructed, which span 12.25, 15.90, and 5.56 cM, with 2, 14, and 4 simple sequence repeat (SSR) and insertion/deletion (Indel) markers, respectively. QTLs of FL, FS, FM, and LP on these three intervals were verified by composite interval mapping (CIM) using WinQTL Cartographer 2.5 software via phenotyping of F2 and its derived F2:3 populations. The results validated the previous primary QTL identification of FL, FS, FM, and LP. Analysis of the RNA-seq data of the developing fibers of L28 and RIL40 at 10, 20, and 30 days post anthesis (DPA) identified seven differentially expressed genes (DEGs) as potential candidate genes. qRT-PCR verified that five of them were consistent with the RNA-seq result. These genes may be involved in regulating fiber development, leading to the formation of FL, FS, FM, and LP. This study provides an experimental foundation for further exploration of these functional genes to dissect the genetic mechanism of cotton fiber development.

Sulfotransferases (SOTs) (EC 2.8.2.-) are sulfate regulatory proteins in a variety of organisms that have been previously shown to be involved in regulating a variety of physiological and biological processes, such as growth, development, adaptation to land, stomatal closure, drought tolerance, and response to pathogen infection. However, there is a lack of comprehensive identification and systematic analysis of SOT in cotton, especially in G. barbadense. In this study, we used bioinformatics methods to analyze the structural characteristics, phylogenetic relationships, gene structure, expression patterns, evolutionary relationships, selection pressure and stress response of SOT gene family members in G. barbadense. In this study, a total of 241 SOT genes were identified in four cotton species, among which 74 SOT gene members were found in G. barbadense. According to the phylogenetic tree, 241 SOT protein sequences were divided into five distinct subfamilies. We also mapped the physical locations of these genes on chromosomes and visualized the structural information of SOT genes in G. barbadense. We also predicted the cis-acting elements of the SOT gene in G. barbadense, and we identified the repetitive types and collinearity analysis of SOT genes in four cotton species. We calculated the Ka/Ks ratio between homologous gene pairs to elucidate the selective pressure between SOT genes. Transcriptome data were used to explore the expression patterns of SOT genes, and then qRT-PCR was used to detect the expression patterns of GBSOT4, GBSOT17 and GBSOT33 under FOV stress. WGCNA (weighted gene co-expression network analysis) showed that GB_A01G0479 (GBSOT4) belonged to the MEblue module, which may regulate the resistance mechanism of G. barbadense to FOV through plant hormones, signal transduction and glutathione metabolism. In addition, we conducted a VIGS (virus-induced gene silencing) experiment on GBSOT4, and the results showed that after FOV inoculation, the plants with a silenced target gene had more serious leaf wilting, drying and cracking than the control group, and the disease index of the plants with the silenced target gene was significantly higher than that of the control group. This suggests that GBSOT4 may be involved in protecting the production of G. barbadense from FOV infection. Subsequent metabolomics analysis showed that some flavonoid metabolites, such as Eupatorin-5-methylether (3′-hydroxy-5,6,7,4′-tetramethoxyflavone, were accumulated in cotton plants in response to FOV infection.

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

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.

Background Cotton is an economically important crop in China, and drought has seriously affected cotton production. Understanding genetic variation, genotype ×environment interactions, and the associations between these traits is critical for developing improved cotton varieties with high drought tolerance. Methods To screen ideal drought-resistant cotton germplasm lines and excellent genotypes, the yield traits of 103 cotton germplasm lines were analyzed. Cotton resource material was planted under normal watering and water deficit conditions for three consecutive years. The yield traits under normal irrigation and water stress conditions were measured, and then five screening indicators were calculated based on the cotton yield per plant under the two water treatments to determine the ideal genotype and most accurate identification indicators. Results The results of correlation analysis and principal component analysis showed that the geometric mean productivity (GMP), mean productivity (MP), and stress tolerance index (STI) were significantly positively correlated with yield under water stress and could be used to distinguish genotypes with high drought tolerance. Among the experimental germplasm lines, some had higher STI and GMP values, indicating their higher drought tolerance. This result indicates that best linear unbiased prediction (BLUP) analysis of the STI and GMP under drought stress can effectively improve screening for drought tolerance in cotton germplasm lines. The results from the screening index, three-dimensional map, and genotype ×environment (GGE) biplots were consistent with the above results. We determined that CQJ-5, Xin lu zao 45, Bellsno, Zhong R 2016 and ND 359-5 are drought-tolerant genotypes that can be used to breed drought-tolerant germplasm lines that produce high and stable yields.

The subtilisin-like proteases (SBTs) are a large family of serine peptidases that are unique to plants. Previous studies have shown that SBTs are associated with developmental processes and environmental responses. However, comprehensive identification and systematic analysis of the SBT family have not been conducted in cotton. We used bioinformatics methods to analyze the structural characteristics, phylogenetic relationships, gene structures, expression modes, evolutionary relationships, selection pressures and stress responses of SBT gene family members in upland cotton. In this study, we identified 120 and 112 SBTs in the tetraploid cotton species G. hirsutum and G. barbadense, while 67 and 69 SBTs were identified in the diploid species G. arboreum and G. raimondii, respectively; these SBTs were divided into five distinct subfamilies. We identified the SBT gene GhSBT27A, and explore its function through virus-induced gene silencing and transmission electron microscopy. These results suggested that the GhSBT27A gene was involved in the response to drought stress. These results lay a foundation for further study on the drought stress mechanism of cotton.

Plant virus-mediated sgRNA delivery and expression have great advantages; sgRNA expression can rapidly expand and accumulate along with virus replication and movement, resulting in efficient gene editing efficiency. In this study, a VIGE system based on cotton leaf crumple virus (CLCrV) was established using cotton overexpressing Cas9 (Cas9-OE) as the VIGE receptor. CLCrV-mediated VIGE could not only target and knock out the GhMAPKKK2, GhCLA1 and GhPDS genes subgroup A and D genome sequences but also achieve double mutation of GhCLA1 and GhPDS genes at the same time. These results verified the effectiveness and efficiency of this system. In addition, the off-target effect assay demonstrated that the CLCrV-mediated VIGE system not only has high gene editing efficiency but also high gene editing specificity in cotton. We further explored whether the FT-sgRNA strategy could transport sgRNA to cotton apical meristem (SAM) over long distances to avoid using tissue culture to obtain stable genetic mutants. The results showed that the sgRNA fused with FT mRNA at the 5’ end could also efficiently achieve targeted editing of endogenous genes in cotton, but it was difficult to detect heritable mutant progeny. The above results showed that the CLCrV-mediated VIGE system provided an accurate and rapid validation tool for screening effective sgRNAs in cotton.

Abstract Worldwide, Verticillium wilt is among the major harmful diseases in cotton production, causing substantial reduction in yields. While this disease has been extensively researched at the molecular level of the pathogen, the molecular basis of V. dahliae host response association is yet to be thoroughly investigated. In this study, RNA-seq analysis was carried out on V. dahliae infected two Gossypium hirsutum L. cultivars, Xinluzao-36 (susceptible) and Zhongzhimian-2 (disease resistant) for 0 h, 24 h, 72 h and 120 h time intervals. Statistical analysis revealed that V. dahliae infection elicited differentially expressed gene responses in the two cotton varieties, but more intensely in the susceptible cultivar than in the resistant cultivars. Data analysis revealed 4241 differentially expressed genes (DEGs) in the LT variety across the three treatment timepoints whereas 7657 in differentially expressed genes (DEGs) in the Vd592 variety across the three treatment timepoints. Six genes were randomly selected for qPCR validation of the RNA-Seq data. Numerous genes encompassed in disease resistance and defense mechanisms were identified. Further, RNA-Seq dataset was utilized in construction of the weighted gene co-expression network and 11 hub genes were identified, that encode for different proteins associated with lignin and immune response, Auxin response factor, cell wall and vascular development, microtubule, Ascorbate transporter, Serine/threonine kinase and Immunity and drought were identified. This significant research will aid in advancing crucial knowledge on virus-host interactions and identify key genes intricate in G. hirsutum L. resistance to V. dahliae infection.

Background PIN proteins are an important class of auxin polar transport proteins that play an important regulatory role in plant growth and development. However, their characteristics and functions have not been identified in Gossypium barbadense. Methods PIN family genes were identified in the cotton species G. barbadense, Gossypium hirsutum, Gossypium raimondii, and Gossypium arboreum, and detailed bioinformatics analyses were conducted to explore the roles of these genes in G. barbadense using transcriptome data and quantitative reverse-transcription polymerase chain reaction (qRT-PCR) technology. Functional verification of the genes was performed using virus-induced gene silencing (VIGS) technology. Results A total of 138 PIN family genes were identified in the four cotton species; the genes were divided into seven subgroups. GbPIN gene family members were widely distributed on 20 different chromosomes, and most had repeated duplication events. Transcriptome analysis showed that some genes had differential expression patterns in different stages of fiber development. According to ‘PimaS-7’ and ‘5917’ transcript component association analysis, the transcription of five genes was directly related to endogenous auxin content in cotton fibers. qRT-PCR analysis showed that the GbPIN7 gene was routinely expressed during fiber development, and there were significant differences among materials. Transient silencing of the GbPIN7 gene by VIGS led to significantly higher cotton plant growth rates and significantly lower endogenous auxin content in leaves and stems. This study provides comprehensive analyses of the roles of PIN family genes in G. barbadense and their expression during cotton fiber development. Our results will form a basis for further PIN auxin transporter research.

Abstract Background An evolutionary model using diploid and allotetraploid cotton species identified 80 % of non-coding transcripts in allotetraploid cotton as being uniquely activated in comparison with its diploid ancestors. The function of the lncRNAs activated in allotetraploid cotton remain largely unknown. Results We employed transcriptome analysis to examine the relationship between the lncRNAs and mRNAs of protein coding genes (PCGs) in cotton leaf tissue under abiotic stresses. LncRNA expression was preferentially associated with that of the flanking PCGs. Selected highly-expressed lncRNA candidates (n = 111) were subjected to a functional screening pilot test in which virus-induced gene silencing was integrated with abiotic stress treatment. From this low-throughput screen, we obtained candidate lncRNAs relating to plant height and tolerance to drought and other abiotic stresses. Conclusions Low-throughput screen is an effective method to find functional lncRNA for further study. LncRNAs were more active in abiotic stresses than PCG expression, especially temperature stress. LncRNA XLOC107738 may take a cis-regulatory role in response to environmental stimuli. The degree to which lncRNAs are constitutively expressed may impact expression patterns and functions on the individual gene level rather than in genome-wide aggregate.

Abstract Background The Raffinose synthetase (RAFS) genes superfamily is critical for the synthesis of raffinose, which accumulates in plant leaves under abiotic stress. However, it remains unclear whether RAFS contributes to resistance to abiotic stress in plants, specifically in the Gossypium species. Results In this study, we identified 74 RAFS genes from G. hirsutum, G. barbadense, G. arboreum and G. raimondii by using a series of bioinformatic methods. Phylogenetic analysis showed that the RAFS gene family in the four Gossypium species could be divided into four major clades; the relatively uniform distribution of the gene number in each species ranged from 12 to 25 based on species ploidy, most likely resulting from an ancient whole-genome polyploidization. Gene motif analysis showed that the RAFS gene structure was relatively conservative. Promoter analysis for cis-regulatory elements showed that some RAFS genes might be regulated by gibberellins and abscisic acid, which might influence their expression levels. Moreover, we further examined the functions of RAFS under cold, heat, salt and drought stress conditions, based on the expression profile and co-expression network of RAFS genes in Gossypium species. Transcriptome analysis suggested that RAFS genes in clade III are highly expressed in organs such as seed, root, cotyledon, ovule and fiber, and under abiotic stress in particular, indicating the involvement of genes belonging to clade III in resistance to abiotic stress. Gene co-expressed network analysis showed that GhRFS2A-GhRFS6A, GhRFS6D, GhRFS7D and GhRFS8A-GhRFS11A were key genes, with high expression levels under salt, drought, cold and heat stress. Conclusion The findings may provide insights into the evolutionary relationships and expression patterns of RAFS genes in Gossypium species and a theoretical basis for the identification of stress resistance materials in cotton.

Abstract Background Domain of unknown function 668 (DUF668) may play a crucial role in the plant growth and developmental response to adverse stress. However, our knowledge of the function of the DUF668 gene family is limited. Results Our study was conducted based on the DUF668 gene family identified from cotton genome sequencing. Phylogenetic analysis showed that the DUF668 family genes can be classified into four subgroups in cotton. We identified 32 DUF668 genes, which are distributed on 17 chromosomes and most of them located in the nucleus of Gossypium hirsutum. Gene structure and motif analyses revealed that the members of the DUF668 gene family can be clustered in G. hirsutum into two broad groups, which are relatively evolutionarily conserved. Transcriptome data analysis showed that the GhDUF668 genes are differentially expressed in different tissues under various stresses (cold, heat, drought, salt, and Verticillium dahliae), and expression is generally increased in roots and stems. Promoter and expression analyses indicated that Gh_DUF668–05, Gh_DUF668–08, Gh_DUF668–11, Gh_DUF668–23 and Gh_DUF668–28 in G. hirsutum might have evolved resistance to adverse stress. Additionally, qRT-PCR revealed that these 5 genes in four cotton lines, KK1543 (drought resistant), Xinluzao 26 (drought sensitive), Zhongzhimian 2 (disease resistant) and Simian 3 (susceptible), under drought and Verticillium wilt stress were all significantly induced. Roots had the highest expression of these 5 genes before and after the treatment. Among them, the expression levels of Gh_DUF668–08 and Gh_DUF668–23 increased sharply at 6 h and reached a maximum at 12 h under biotic and abiotic stress, which showed that they might be involved in the process of adverse stress resistance in cotton. Conclusion The significant changes in GhDUF668 expression in the roots after adverse stress indicate that GhDUF668 is likely to increase plant resistance to stress. This study provides an important theoretical basis for further research on the function of the DUF668 gene family and the molecular mechanism of adverse stress resistance in cotton.

Abstract 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 (Liy), 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.

The increasing water scarcity associated with environmental change brings significant negative impacts to the growth of cotton plants, whereby it is urgent to enhance plant tolerance to drought. Here, we overexpressed the com58276 gene isolated from the desert plant Caragana korshinskii in cotton plants. We obtained three OE plants and demonstrated that com58276 confers drought tolerance in cotton after subjecting transgenic seeds and plants to drought. RNA-seq revealed the mechanisms of the possible anti-stress response, and that the overexpression of com58276 does not affect growth and fiber content in OE cotton plants. The function of com58276 is conserved across species, improving the tolerance of cotton to salt and low temperature, and demonstrating its applicability to improve plant resistance to environmental change.

Abstract Background Cotton (Gossypium hirsutum) is considered a fairly salt tolerant crop however, salinity can still cause significant economic losses by affecting the yield and deteriorating the fiber quality. We studied a salt-tolerant upland cotton cultivar under temporal salt stress to unfold the salt tolerance molecular mechanisms. Biochemical response to salt stress (400 mM) was measured at 0 h, 3 h, 12 h, 24 h and 48 h post stress intervals and single-molecule long-read sequencing technology from Pacific Biosciences (PacBio) combined with the unique molecular identifiers approach was used to identify differentially expressed genes (DEG). Results Antioxidant enzymes including, catalase (CAT), peroxidase (POD), superoxide dismutase (SOD) were found significantly induced under temporal salt stress, suggesting that reactive oxygen species scavenging antioxidant machinery is an essential component of salt tolerance mechanism in cotton. We identified a wealth of novel transcripts based on the PacBio long reads sequencing approach. Prolonged salt stress duration induces high number of DEGs. Significant numbers of DEGs were found under key terms related to stress pathways such as “response to oxidative stress”, “response to salt stress”, “response to water deprivation”, “cation transport”, “metal ion transport”, “superoxide dismutase”, and “reductase”. Key DEGs related to hormone (abscisic acid, ethylene and jasmonic acid) biosynthesis, ion homeostasis (CBL-interacting serine/threonine-protein kinase genes, calcium-binding proteins, potassium transporter genes, potassium channel genes, sodium/hydrogen exchanger or antiporter genes), antioxidant activity (POD, SOD, CAT, glutathione reductase), transcription factors (myeloblastosis, WRKY, Apetala 2) and cell wall modification were found highly active in response to salt stress in cotton. Expression fold change of these DEGs showed both positive and negative responses, highlighting the complex nature of salt stress tolerance mechanisms in cotton. Conclusion Collectively, this study provides a good insight into the regulatory mechanism under salt stress in cotton and lays the foundation for further improvement of salt stress tolerance.

Fusarium wilt caused by Fusarium oxysporum f. sp. vasinfectum (FOV) is one of the most destructive diseases in cotton (Gossypium spp.) production, and use of resistant cultivars is the most cost-effective method managing the disease. To understand the genetic basis of cotton resistance to FOV race 7 (FOV7), this study evaluated a recombinant inbred line (RIL) population of 110 lines of G. barbadense from a cross between susceptible Xinhai 14 and resistant 06-146 in eight tests and constructed a high-density genetic linkage map with resequencing-based 933,845 single-nucleotide polymorphism (SNP) markers covering a total genetic distance of 2483.17 cM. Nine quantitative trait loci (QTLs) for FOV7 resistance were identified, including qFOV7-D03-1 on chromosome D03 in two tests. Through a comparative analysis of gene expression and DNA sequence for predicted genes within the QTL region between the two parents and selected lines inoculated with FOV7, GB_D03G0217 encoding for a calmodulin (CaM)-like (CML) protein was identified as a candidate gene. A further analysis confirmed that the expression of GB_D03G0217 was suppressed, leading to increased disease severity in plants of the resistant parent with virus induced gene silencing (VIGS).

Detached organs or differentiated tissues could form a mass of pluripotent cells termed as callus on an auxin-rich medium, the underlying molecular mechanism of which remains elusive in cotton. LATERAL ORGAN BOUNDARIES DOMAIN (LBD) transcription factor is a key regulator of plant cell totipotency/pluripotency, and a number of cotton GhLBDs with high-level differential expression during the callus induction process have been identified. Their overexpression in cotton calli fostered promotions in and callus induction without exogenous auxin. Expression analysis and histological observation using paraffin sectioning suggested that the first 72 h on culture is a key time point for callus initiation, whereby the GhLBDs showed high transcript abundance and enlarged calli that were rapidly developed from procambium and cambium. GhLBDs’ expression level could be precisely modulated by the gradient concentrations of exogenous auxin, whereas auxin transport inhibitor 2,3,5-triiodobenzoic acid could severely inhibit its expression. The LBD-mediated callus formation was also dependent on the expression levels of GhLBDs. Further, a β-estradiol-inducible promoter pER8 was used to drive GhLBD18-1 expression, which led to rapid callus proliferation, suggesting that pER8/GhLBD18-1 could be used as a selectable marker system to replace the existing antibiotic/herbicide-resistance selectable markers in plant transformation. Our study provides new insights for callus initiation regulatory mechanism and strategies for improving transformation efficiency in cotton.

Abstract Background Sea island cotton (Gossypium barbadense) has markedly superior high quality fibers, which plays an important role in the textile industry and acts as a donor for upland cotton (G. hirsutum) fiber quality improvement. The genetic characteristics analysis and the identification of key genes will be helpful to understand the mechanism of fiber development and breeding utilization in sea island cotton. Results In this study, 279 sea island cotton accessions were collected from different origins for genotyping and phenotyping fiber quality traits. A set of 6303 high quality single nucleotide polymorphisms (SNPs) were obtained by high-density CottonSNP80K array. The population characteristics showed that the sea island cotton accessions had wide genetic diversity and were clustered into three groups, with Group1 closely related to Menoufi, an original sea island cotton landrace, and Group2 and Group3 related to widely introduced accessions from Egypt, USA and Former Soviet Union. Further, we used 249 accessions and evaluated five fiber quality traits under normal and salt environments over 2 years. Except for fiber uniformity (FU), fiber length (FL) and fiber elongation (FE) were significantly decreased in salt conditions, while fiber strength (FS) and fiber micronaire (MIC) were increased. Based on 6303 SNPs and genome-wide association study (GWAS) analysis, a total of 34 stable quantitative trait loci (QTLs) were identified for the five fiber quality traits with 25 detected simultaneously under normal and salt environments. Gene Ontology (GO) analysis indicated that candidate genes in the 25 overlapped QTLs were enriched mostly in “cellular and biological process”. In addition, “xylem development” and “response to hormone” pathways were also found. Haplotype analyses found that GB_A03G0335 encoding an E3 ubiquitin-protein ligase in QTL TM6004 had SNP variation (A/C) in gene region, was significantly correlated with FL, FS, FU, and FE, implying a crucial role in fiber quality. Conclusions The present study provides a foundation for genetic diversity of sea island cotton accessions and will contribute to fiber quality improvement in breeding practice.

Abstract Background Plant height (PH) and fruit branch number (FBN) are important traits for improving yield and mechanical harvesting of cotton. In order to identify genes of PH and FBN in cotton germplasms to develop superior cultivars, quantitative trait loci (QTLs) for these traits were detected based on the phenotypic evaluation data in nine environments across four locations and 4 years and a previously reported genetic linkage map of an recombinant inbred line (RIL) population of upland cotton. Results In total, 53 QTLs of PH and FBN, were identified on 21 chromosomes of the cotton genome except chromosomes c02, c09-c11, and c22. For PH, 27 QTLs explaining 3.81%–8.54% proportions of phenotypic variance were identified on 18 chromosomes except c02, c08-c12, c15, and c22. For FBN, 26 QTLs explaining 3.23%–11.00% proportions of phenotypic variance were identified on 16 chromosomes except c02-c03, c06, c09-c11, c17, c22-c23, and c25. Eight QTLs were simultaneously identified in at least two environments. Three QTL clusters containing seven QTLs were identified on three chromosomes (c01, c18 and c21). Eleven QTLs were the same as previously reported ones, while the rest were newly identified. Conclusions The QTLs and QTL clusters identified in the current study will be helpful to further understand the genetic mechanism of PH and FBN development of cotton and will enhance the development of excellent cultivars for mechanical managements in cotton production.

Abstract Background Salinity is one of the most significant environmental factors limiting the productivity of cotton. However, the key genetic components responsible for the reduction in cotton yield in saline-alkali soils are still unclear. Results Here, we evaluated three main components of lint yield, single boll weight (SBW), lint percentage (LP) and boll number per plant (BNPP), across 316 G. hirsutum accessions under four salt conditions over two years. Phenotypic analysis indicated that LP was unchanged under different salt conditions, however BNPP decreased significantly and SBW increased slightly under high salt conditions. Based on 57,413 high-quality single nucleotide polymorphisms (SNPs) and genome-wide association study (GWAS) analysis, a total of 42, 91 and 25 stable quantitative trait loci (QTLs) were identified for SBW, LP and BNPP, respectively. Phenotypic and QTL analysis suggested that there was little correlation among the three traits. For LP, 8 stable QTLs were detected simultaneously in four different salt conditions, while fewer repeated QTLs for SBW or BNPP were identified. Gene Ontology (GO) analysis indicated that their regulatory mechanisms were also quite different. Via transcriptome profile data, we detected that 10 genes from the 8 stable LP QTLs were predominantly expressed during fiber development. Further, haplotype analyses found that a MYB gene (GhMYB103), with the two SNP variations in cis-regulatory and coding regions, was significantly correlated with lint percentage, implying a crucial role in lint yield. We also identified that 40 candidate genes from BNPP QTLs were salt-inducible. Genes related to carbohydrate metabolism and cell structure maintenance were rich in plants grown in high salt conditions, while genes related to ion transport were active in plants grown in low salt conditions, implying different regulatory mechanisms for BNPP at high and low salt conditions. Conclusions This study provides a foundation for elucidating cotton salt tolerance mechanisms and contributes gene resources for developing upland cotton varieties with high yields and salt stress tolerance.

Upland cotton (Gossypium hirsutum) is widely planted around the world for its natural fiber, and producing high-quality fiber is essential for the textile industry. CCRI70 is a hybrid cotton plant harboring superior yield and fiber quality, whose recombinant inbred line (RIL) population was developed from two upland cotton varieties (sGK156 and 901-001) and were used here to investigate the source of high-quality related alleles. Based on the material of the whole population, a high-density genetic map was constructed using specific locus-amplified fragment sequencing (SLAF-seq). It contained 24,425 single nucleotide polymorphism (SNP) markers, spanning a distance of 4,850.47 centimorgans (cM) over 26 chromosomes with an average marker interval of 0.20 cM. In evaluating three fiber quality traits in nine environments to detect multiple environments stable quantitative trait loci (QTLs), we found 289 QTLs, of which 36 of them were stable QTLs and 18 were novel. Based on the transcriptome analysis for two parents and two RILs, 24,941 unique differentially expressed genes (DEGs) were identified, 473 of which were promising genes. For the fiber strength (FS) QTLs, 320 DEGs were identified, suggesting that pectin synthesis, phenylpropanoid biosynthesis, and plant hormone signaling pathways could influence FS, and several transcription factors may regulate fiber development, such as GAE6, C4H, OMT1, AFR18, EIN3, bZIP44, and GAI. Notably, the marker D13_56413025 in qFS-chr18-4 provides a potential basis for enhancing fiber quality of upland cotton via marker-assisted breeding and gene cloning of important fiber quality traits.

Cotton is a globally important cash crop. Verticillium wilt (VW) is commonly known as “cancer” of cotton and causes serious loss of yield and fiber quality in cotton production around the world. Here, we performed a BSA-seq analysis using an F2:3 segregation population to identify the candidate loci involved in VW resistance. Two QTLs (qvw-D05-1 and qvw-D05-2) related to VW resistance in cotton were identified using two resistant/susceptible bulks from the F2 segregation population constructed by crossing the resistant cultivar ZZM2 with the susceptible cultivar J11. A total of 30stop-lost SNPs and 42 stop-gained SNPs, which included 17 genes, were screened in the qvw-D05-2 region by SnpEff analysis. Further analysis of the transcriptome data and qRT-PCR revealed that the expression level of Ghir_D05G037630 (designated as GhDRP) varied significantly at certain time points after infection with V. dahliae. The virus-induced gene silencing of GhDRP resulted in higher susceptibility of the plants to V. dahliae than the control, suggesting that GhDRP is involved in the resistance to V. dahlia infection. This study provides a method for rapid mining of quantitative trait loci and screening of candidate genes, as well as enriches the genomic information and gene resources for the molecular breeding of disease resistance in cotton.

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.

The APETALA2 (AP2)/ethylene response factor plays vital functions in response to environmental stimulus. The ethylene response factor (ERF) subfamily B3 group belongs to the AP2/ERF superfamily and contains a single AP2/ERF domain. Phylogenetic analysis of the ERF subfamily B3 group genes from Arabdiposis thaliana, Gossypium arboreum, Gossypium hirsutum, and Gossypium raimondii made it possible to divide them into three groups and showed that the ERF subfamily B3 group genes are conserved in cotton. Collinearity analysis identified172 orthologous/paralogous gene pairs between G. arboreum and G. hirsutum; 178 between G. hirsutum and G. raimondii; and 1,392 in G. hirsutum. The GhERF subfamily B3 group gene family experienced massive gene family expansion through either segmental or whole genome duplication events, with most genes showing signature compatible with the action of purifying selection during evolution. Most G. hirsutum ERF subfamily B3 group genes are responsive to salt stress. GhERF13.12 transgenic Arabidopsis showed enhanced salt stress tolerance and exhibited regulation of related biochemical parameters and enhanced expression of genes participating in ABA signaling, proline biosynthesis, and ROS scavenging. In addition, the silencing of the GhERF13.12 gene leads to increased sensitivity to salt stress in cotton. These results indicate that the ERF subfamily B3 group had remained conserved during evolution and that GhERF13.12 induces salt stress tolerance in Arabidopsis and cotton.

Protein palmitoylation is an ability of the frame of the cell marker protein is one of the most notable reversible changes after translation. However, studies on protein palmitoylation in cotton have not yet been performed. In our current research, the PAT gene family was systematically identified and bioinformatically analyzed in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, and 211 PAT genes were authenticated and classified into six subfamilies. Sixty-nine PAT genes were identified in upland cotton, mainly at the ends of its the 26 chromosomes of upland cotton. The majority of these genes are located in the nucleus of the plant. Gene structure analysis revealed that each member encodes a protein that which contains at least one DHHC structural domain. Cis-acting element analysis indicated that GhPATs genes are mainly involved in hormone production, light response and stress response. Gene expression pattern analysis indicated that most GhPATs genes were differentially expressed upon induction by pathogenic bacteria, drought, salt, hot and cold stresses, and some GhPATs could be induced by multiple abiotic stresses simultaneously. GhPATs genes showed different expression patterns in tissue-specific assays and were found to be preferentially expressed in roots, followed by expression in stems and leaves. Virus-induced gene silencing (VIGS) experiments showed that cotton was significantly less resistant to Verticillium dahliae when GhPAT27 was silenced. We conclude that the GhPAT27 gene, which mediates S-palmitoylation acetylation, may be involved in the regulation of upland cotton resistance to Verticillium wilt (VW). Overall, this work has provided a fundamental framework for understanding the latent capabilities of GhPATs and a solid foundation for molecular breeding and plant pathogen resistance in cotton.

Drought is one of the main abiotic stresses that seriously influences cotton production. Many indicators can be used to evaluate cotton drought tolerance, but the key indicators remain to be determined. The objective of this study was to identify effective cotton drought tolerance indicators from 19 indices, including morphology, photosynthesis, physiology, and yield-related indices, and to evaluate the yield potential of 104 cotton varieties under both normal and drought-stress field conditions. Combined with principal component analysis (PCA) and a regression analysis method, the results showed that the top five PCs among the 19, with eigenvalues > 1, contributed 65.52, 63.59, and 65.90% of the total variability during 2016 to 2018, respectively, which included plant height (PH), effective fruit branch number (EFBN), single boll weight (SBW), transpiration rate (Tr) and chlorophyll (Chl). Therefore, the indicator dimension decreased from 19 to 5. A comparison of the 19 indicators with the 5 identified indicators through PCA and a combined regression analysis found that the results of the final cluster of drought tolerance on 104 cotton varieties were basically consistent. The results indicated that these five traits could be used in combination to screen cotton varieties or lines for drought tolerance in cotton breeding programs, and Zhong R2016 and Xin lu zao 45 exhibited high drought tolerance and can be selected as superior parents for good yield performance under drought stress.

Heat shock transcription factors (HSFs) are involved in environmental stress response and plant development, such as heat stress and flowering development. According to the structural characteristics of the HSF gene family, HSF genes were classified into three major types (HSFA, HSFB, and HSFC) in plants. Using conserved domains of HSF genes, we identified 621 HSF genes among 13 cotton genomes, consisting of eight diploid and five tetraploid genomes. Phylogenetic analysis indicated that HSF genes among 13 cotton genomes were grouped into two different clusters: one cluster contained all HSF genes of HSFA and HSFC, and the other cluster contained all HSF genes of HSFB. Comparative analysis of HSF genes in Arabidopsis thaliana, Gossypium herbaceum (A1), Gossypium arboreum (A2), Gossypium raimondii (D5), and Gossypium hirsutum (AD1) genomes demonstrated that four HSF genes were inherited from a common ancestor, A0, of all existing cotton A genomes. Members of the HSF gene family in G. herbaceum (A1) genome indicated a significant loss compared with those in G. arboretum (A2) and G. hirsutum (AD1) A genomes. However, HSF genes in G. raimondii (D5) showed relative loss compared with those in G. hirsutum (AD1) D genome. Analysis of tandem duplication (TD) events of HSF genes revealed that protein-coding genes among different cotton genomes have experienced TD events, but only the two-gene tandem array was detected in Gossypium thurberi (D1) genome. The expression analysis of HSF genes in G. hirsutum (AD1) and Gossypium barbadense (AD2) genomes indicated that the expressed HSF genes were divided into two different groups, respectively, and the expressed HSF orthologous genes between the two genomes showed totally different expression patterns despite the implementation of the same abiotic stresses. This work will provide novel insights for the study of evolutionary history and expression characterization of HSF genes in different cotton genomes and a widespread application model for the study of HSF gene families in plants.

Anthocyanidin reductase (ANR) is an important regulator of flavonoid metabolism, and proanthocyanidins, the secondary metabolites of flavonoids, play an important role in the response of plants to pathogenic stress. Therefore, in this study, the expression analysis of the ANR gene family of Gossypium barbadense after inoculation with Fusarium oxysporum f. sp. vasinfectum (FOV) was performed at different time points. It was found that Gb_ANR-47 showed significant differences in the disease-resistant cultivar 06-146 and the susceptible cultivar Xinhai 14, as well as in the highest root expression. It was found that the expression of Gb_ANR-47 in the resistant cultivar was significantly higher than that in the susceptible cultivar by MeJA and SA, and different amounts of methyl jasmonate (MeJA) and salicylic acid (SA) response elements were found in the promoter region of Gb_ANR-47. After silencing GbANR-47 in 06-146 material by VIGS technology, its resistance to FOV decreased significantly. The disease severity index (DSI) was significantly increased, and the anthocyanin content was significantly decreased in silenced plants, compared to controls. Our findings suggest that GbANR-47 is a positive regulator of FOV resistance in Gossypium barbadense. The research results provide an important theoretical basis for in-depth analysis of the molecular mechanism of GbANR-47 and improving the anti-FOV of Gossypium barbadense.

Ethylene plays a pivotal role in plant stress resistance and 1-aminocyclopropane-1-carboxylic acid synthase (ACS) is the rate-limiting enzyme in ethylene biosynthesis. Upland cotton (Gossypium hirsutum L.) is the most important natural fiber crop, but the function of ACS in response to abiotic stress has rarely been reported in this plant. We identified 18 GaACS, 18 GrACS, and 35 GhACS genes in Gossypiumarboreum, Gossypium raimondii and Gossypiumhirsutum, respectively, that were classified as types I, II, III, or IV. Collinearity analysis showed that the GhACS genes were expanded from diploid cotton by the whole-genome-duplication. Multiple alignments showed that the C-terminals of the GhACS proteins were conserved, whereas the N-terminals of GhACS10 and GhACS12 were different from the N-terminals of AtACS10 and AtACS12, probably diverging during evolution. Most type II ACS genes were hardly expressed, whereas GhACS10/GhACS12 were expressed in many tissues and in response to abiotic stress; for example, they were highly and hardly expressed at the early stages of cold and heat exposure, respectively. The GhACS genes showed different expression profiles in response to cold, heat, drought, and salt stress by quantitative PCR analysis, which indicate the potential roles of them when encountering the various adverse conditions, and provide insights into GhACS functions in cotton’s adaptation to abiotic stress.

Abstract 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.

Abstract Background Gossypium barbadense (Sea Island, Egyptian or Pima cotton) cotton has high fiber quality, however, few studies have investigated the genetic basis of its traits using molecular markers. Genome complexity reduction approaches such as genotyping-by-sequencing have been utilized to develop abundant markers for the construction of high-density genetic maps to locate quantitative trait loci (QTLs). Results The Chinese G. barbadense cultivar 5917 and American Pima S-7 were used to develop a recombinant inbred line (RIL) population with 143 lines. The 143 RILs together with their parents were tested in three replicated field tests for lint yield traits (boll weight and lint percentage) and fiber quality traits (fiber length, fiber elongation, fiber strength, fiber uniformity and micronaire) and then genotyped using GBS to develop single-nucleotide polymorphism (SNP) markers. A high-density genetic map with 26 linkage groups (LGs) was constructed using 3557 GBS SNPs spanning a total genetic distance of 3076.23 cM at an average density of 1.09 cM between adjacent markers. A total of 42 QTLs were identified, including 24 QTLs on 12 LGs for fiber quality and 18 QTLs on 7 LGs for lint yield traits, with LG1 (9 QTLs), LG10 (7 QTLs) and LG14 (6 QTLs) carrying more QTLs. Common QTLs for the same traits and overlapping QTLs for different traits were detected. Each individual QTLs explained 0.97 to 20.7% of the phenotypic variation. Conclusions This study represents one of the first genetic mapping studies on the fiber quality and lint yield traits in a RIL population of G. barbadense using GBS-SNPs. The results provide important information for the subsequent fine mapping of QTLs and the prediction of candidate genes towards map-based cloning and marker-assisted selection in cotton.

Extreme high temperatures are threatening cotton production around the world due to the intensification of global warming. To cope with high-temperature stress, heat-tolerant cotton cultivars have been bred, but the heat-tolerant mechanism remains unclear. This study selected heat-tolerant (‘Xinluzao36′) and heat-sensitive (‘Che61-72′) cultivars of cotton treated with high-temperature stress as plant materials and performed comparative nanopore sequencing transcriptome analysis to reveal the potential heat-tolerant mechanism of cotton. Results showed that 120,605 nonredundant sequences were generated from the raw reads, and 78,601 genes were annotated. Differentially expressed gene (DEG) analysis showed that a total of 19,600 DEGs were screened; the DEGs involved in the ribosome, heat shock proteins, auxin and ethylene signaling transduction, and photosynthesis pathways may be attributed to the heat tolerance of the heat-tolerant cotton cultivar. This study also predicted a total of 5118 long non-coding RNAs (lncRNAs)and 24,462 corresponding target genes. Analysis of the target genes revealed that the expression of some ribosomal, heat shock, auxin and ethylene signaling transduction-related and photosynthetic proteins may be regulated by lncRNAs and further participate in the heat tolerance of cotton. This study deepens our understandings of the heat tolerance of cotton.

It is of great importance to identify quantitative trait loci (QTL) controlling fiber quality traits and yield components for future marker-assisted selection (MAS) and candidate gene function identifications. In this study, two kinds of traits in 231 F6:8 recombinant inbred lines (RILs), derived from an intraspecific cross between Xinluzao24, a cultivar with elite fiber quality, and Lumianyan28, a cultivar with wide adaptability and high yield potential, were measured in nine environments. This RIL population was genotyped by 122 SSR and 4729 SNP markers, which were also used to construct the genetic map. The map covered 2477.99 cM of hirsutum genome, with an average marker interval of 0.51 cM between adjacent markers. As a result, a total of 134 QTLs for fiber quality traits and 122 QTLs for yield components were detected, with 2.18–24.45 and 1.68–28.27% proportions of the phenotypic variance explained by each QTL, respectively. Among these QTLs, 57 were detected in at least two environments, named stable QTLs. A total of 209 and 139 quantitative trait nucleotides (QTNs) were associated with fiber quality traits and yield components by four multilocus genome-wide association studies methods, respectively. Among these QTNs, 74 were detected by at least two algorithms or in two environments. The candidate genes harbored by 57 stable QTLs were compared with the ones associated with QTN, and 35 common candidate genes were found. Among these common candidate genes, four were possibly “pleiotropic.” This study provided important information for MAS and candidate gene functional studies.

In plant cells, vacuolar H+-ATPases (V-ATPases) are responsible for deacidification of the cytosol and energisation of the secondary transport processes across the tonoplast. A number of V-ATPase subunit genes have been demonstrated to be involved in the regulation of the plant response to water deficit. However, there are no reports on the role of V-ATPase subunit A (VHA-A) in dehydration tolerance of cotton. In this study, cotton GhVHA-A gene was functionally characterized, especially with regard to its role in dehydration stress tolerance. Expression analysis showed that GhVHA-A was differentially expressed in various cotton organs and was induced by dehydration, low temperature, high salinity, and abscisic acid treatment in leaves. We also report that GhVHA-A improve dehydration tolerance in transgenic tobacco and cotton. Virus-induced gene silencing of GhVHA-A decreased the tolerance of cotton plantlets to dehydration stress. Silencing GhVHA-A decreased chlorophyll content and antioxidant enzyme activities and increased malondialdehyde (MDA) content in cotton under dehydration stress. However, transgenic tobacco expressing GhVHA-A exhibited enhanced dehydration resistance, resulting in reduced leaf water loss, higher average root length, and lower MDA levels under dehydration stress. Meanwhile, overexpression of GhVHA-A in tobacco conferred water deficit tolerance by enhancing osmotic adjustment (proline) and the activities of the antioxidant enzymes superoxide dismutase and peroxidase, thereby enhancing reactive oxygen species detoxification. These results suggest that GhVHA-A plays an important role in conferring resistance to dehydration stress. Our results have identified GhVHA-A as a candidate gene for improving dehydration tolerance in plants.

Members of the YABBY gene family, a small plant-specific family of genes, have been proposed to function in specifying abaxial cell fate. Although to date little has been learned about cotton YABBY genes, completion of the cotton genome enables a comprehensive genome-wide analysis of YABBY genes in cotton. Here, a total of 12, 12, and 23 YABBY genes were identified in Gossypium arboreum (2n = 26, A2), G. raimondii (2n = 26, D5), and G. hirsutum (2n = 4x = 52, [AD]t), respectively. Sequence analysis showed that the N-terminal zinc-finger and C-terminal YABBY domains in YABBY proteins are highly conserved among cotton, Arabidopsis, and rice. Eighty-five genes from eight sequenced species naturally clustered into five groups, and the YAB2-like group could be divided into three sub-groups, indicating that YABBYs are highly conserved among the examined species. Orthologs from the At and Dt sub-genomes (where “t” indicates tetraploid) showed good collinearity, indicating that YABBY loci are highly conserved between these two sub-genomes. Whole-genome duplication was the primary cause of upland cotton YABBY gene expansion, segmental duplication played important roles in YABBY gene expansion within the At and Dt sub-genomes, and the YAB5-like group was mainly generated by segmental duplication. The long-terminal repeat retroelements Copia and Gypsy were identified as major transposable elements accompanying the appearance of duplicated YABBY genes, suggesting that transposable element expansion might be involved in gene duplication. Selection pressure analyses using PAML revealed that relaxed purifying selection might be the main impetus during evolution of YABBY genes in the examined species. Furthermore, exon/intron pattern and motif analyses indicated that genes within the same group were significantly conserved between Arabidopsis and cotton. In addition, the expression patterns in different tissues suggest that YABBY proteins may play roles in ovule development because YABBYs are highly expressed in ovules. The expression pattern of YABBY genes showed that approximately half of the YABBYs were down-regulated under different stress treatments. Collectively, our results represent a comprehensive genome-wide study of the YABBY gene family, which should be helpful in further detailed studies on the gene function and evolution of YABBY genes in cotton.

Height is an important trait related to plant architecture and yield potential in bread wheat (Triticum aestivum L.). We previously identified a major quantitative trait locus QPH.caas-6A flanked by simple sequence repeat markers Xbarc103 and Xwmc256 that reduced height by 8.0–10.4%. Here QPH.caas-6A, designated as Rht24, was confirmed using recombinant inbred lines (RILs) derived from a Jingdong 8/Aikang 58 cross. The target sequences of Xbarc103 and Xwmc256 were used as queries to BLAST against International Wheat Genome Sequence Consortium database and hit a super scaffold of approximately 208 Mb. Based on gene annotation of the scaffold, three gene-specific markers were developed to genotype the RILs, and Rht24 was narrowed to a 1.85 cM interval between TaAP2 and TaFAR. In addition, three single nucleotide polymorphism (SNP) markers linked to Rht24 were identified from SNP chip-based screening in combination with bulked segregant analysis. The allelic efficacy of Rht24 was validated in 242 elite wheat varieties using TaAP2 and TaFAR markers. These showed a significant association between genotypes and plant height. Rht24 reduced plant height by an average of 6.0–7.9 cm across environments and were significantly associated with an increased TGW of 2.0–3.4 g. The findings indicate that Rht24 is a common dwarfing gene in wheat breeding, and TaAP2 and TaFAR can be used for marker-assisted selection.