Background: The function of 4-coumarate-CoA ligases (4CL) under abiotic stresses has been studied in plants, however, limited is known about the 4CL genes in cotton (G. hirsutum L.) and their roles in response to drought stress. Results: We performed genome-wide identification of the 4CL genes in G. hirsutum and investigated the expression profiles of the identified genes in various cotton tissues and in response to stress conditions with an aim to identify 4CL gene(s) associated with drought tolerance. We identified 34 putative 4CL genes in G. hirsutum that were clustered into three classes. Genes of the same class usually share a similar gene structure and motif composition. Many cis-elements related to stress and phytohormone responses were found in the promoters of the Gh4CL genes. Of the 34 Gh4CL genes, 26 were induced by at least one abiotic stress and 10 (including Gh4CL7) were up-regulated under the polyethylene glycol (PEG) simulated drought stress conditions. Virus-induced gene silencing (VIGS) in cotton and overexpression (OE) in Arabidopsis thaliana were applied to investigate the biological function of Gh4CL7 in drought tolerance. The Gh4CL7-silencing cotton plants showed more sensitive to drought stress, probably due to decreased relative water content (RWC), chlorophyll content and antioxidative enzyme activity, increased stomatal aperture, and the contents of malondialdehyde (MDA) and hydrogen peroxide (H₂O₂). Arabidopsis lines overexpressing Gh4CL7, however, were more tolerant to drought treatment, which was associated with improved antioxidative enzyme activity, reduced accumulation of MDA and H₂O₂ and up-regulated stress-related genes under the drought stress conditions. In addition, compared to their respective controls, the Gh4CL7-silencing cotton plants and the Gh4CL7-overexpressing Arabidopsis lines had a ~ 20% reduction and a ~ 10% increase in lignin content, respectively. The expression levels of genes related to lignin biosynthesis, including PAL, CCoAOMT, COMT, CCR and CAD, were lower in Gh4CL7-silencing plants than in controls. Taken together, these results demonstrated that Gh4CL7 could positively respond to drought stress and therefore might be a candidate gene for improvement of drought tolerance in cotton. Conclusion: We characterized the 4CL gene family in upland cotton and revealed a role of Gh4CL7 in lignin biosynthesis and drought tolerance.

Green-colored fiber (GCF) is the unique raw material for naturally colored cotton textile but we know little about the pigmentation process in GCF. Here we compared transcriptomes and metabolomes of 12, 18 and 24 days post-anthesis (DPA) fibers from a green fiber cotton accession and its white-colored fiber (WCF) near-isogenic line. We found a total of 2047 non-redundant metabolites in GCF and WCF that were enriched in 80 pathways, including those of biosynthesis of phenylpropanoid, cutin, suberin, and wax. Most metabolites, particularly sinapaldehyde, of the phenylpropanoid pathway had a higher level in GCF than in WCF, consistent with the significant up-regulation of the genes responsible for biosynthesis of those metabolites. Weighted gene co-expression network analysis (WGCNA) of genes differentially expressed between GCF and WCF was used to uncover gene-modules co-expressed or associated with the accumulation of green pigments. Of the 16 gene-modules co-expressed with fiber color or time points, the blue module associated with G24 (i.e., GCF at 24 DPA) was of particular importance because a large proportion of its genes were significantly up-regulated at 24 DPA when fiber color was visually distinguishable between GCF and WCF. A total of 56 hub genes, including the two homoeologous Gh4CL4 that could act in green pigment biosynthesis, were identified among the genes of the blue module that are mainly involved in lipid metabolism, phenylpropanoid biosynthesis, RNA transcription, signaling, and transport. Our results provide novel insights into the mechanisms underlying pigmentation of green fibers and clues for developing cottons with stable green colored fibers.

The clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) gene editing system has been shown to be able to induce highly efficient mutagenesis in the targeted DNA of many plants, including cotton, and has become an important tool for investigation of gene function and crop improvement. Here, we developed a simple and easy to operate CRISPR/Cas9 system and demonstrated its high editing efficiency in cotton by targeting-ALARP, a gene encoding alanine-rich protein that is preferentially expressed in cotton fibers. Based on sequence analysis of the target site in the 10 transgenic cottons containing CRISPR/Cas9, we found that the mutation frequencies of GhALARP-A and GhALARP-D target sites were 71.4–100% and 92.9–100%, respectively. The most common editing event was deletion, but deletion together with large insertion was also observed. Mosaic mutation editing events were detected in most transgenic plants. No off-target mutation event was detected in any the 15 predicted sites analyzed. This study provided mutants for further study of the function of GhALARP in cotton fiber development. Our results further demonstrated the feasibility of use of CRISPR/Cas9 as a targeted mutagenesis tool in cotton, and provided an efficient tool for targeted mutagenesis and functional genomics in cotton.

Background: Cotton Verticillium wilt is one of the most devastating diseases for cotton production in the world. Although this diseases have been widely studied at the molecular level from pathogens, the molecular basis of V. dahliae interacted with cotton has not been well examined. Results: In this study, RNA-seq analysis was carried out on V. dahliae samples cultured by different root exudates from three cotton cultivars (a susceptible upland cotton cultivar, a tolerant upland cotton cultivar and a resistant island cotton cultivar) and water for 0 h, 6 h, 12 h, 24 h and 48 h. Statistical analysis of differentially expressed genes revealed that V. dahliae responded to all kinds of root exudates but more strongly to susceptible cultivar than to tolerant and resistant cultivars. Go analysis indicated that ‘hydrolase activity, hydrolyzing O-glycosyl compounds’ related genes were highly enriched in V. dahliae cultured by root exudates from susceptible cotton at early stage of interaction, suggesting genes related to this term were closely related to the pathogenicity of V. dahliae. Additionally, ‘transmembrane transport’, ‘coenzyme binding’, ‘NADP binding’, ‘cofactor binding’, ‘oxidoreductase activity’, ‘flavin adenine dinucleotide binding’, ‘extracellular region’ were commonly enriched in V. dahliae cultured by all kinds of root exudates at early stage of interaction (6 h and 12 h), suggesting that genes related to these terms were required for the initial steps of the roots infections. Conclusions: Based on the GO analysis results, the early stage of interaction (6 h and 12 h) were considered as the critical stage of V. dahliae-cotton interaction. Comparative transcriptomic analysis detected that 31 candidate genes response to root exudates from cotton cultivars with different level of V. dahliae resistance, 68 response to only susceptible cotton cultivar, and 26 genes required for development of V. dahliae. Collectively, these expression data have advanced our understanding of key molecular events in the V. dahliae interacted with cotton, and provided a framework for further functional studies of candidate genes to develop better control strategies for the cotton wilt disease.

Fatty acyl‐CoA reductases (FAR) are involved in plant wax synthesis and play important roles in plant growth and development. However, little information is available for cotton FAR genes. In this study, we carried out a genome‐wide identification of FAR genes in Upland cotton (Gossypium hirsutum L.) and found 10 GhFARs that form five pairs of homoeologs (GhFAR2A to GhFAR3.4D) distributed on six chromosomes. The 10 GhFARs were separated into two subfamilies. Most GhFARs showed tissue‐specific expression patterns, and at least one GhFAR of each pair of homoeologs was relatively highly expressed in at least one of the tissues investigated. GhFAR3.1 was highly expressed in leaves. The function of GhFAR3.1 in wax accumulation and drought tolerance was analyzed using virus‐induced gene silencing (VIGS). Silencing GhFAR3.1 reduced the total wax content and relative water content of leaves by over 60% and 13%, respectively, suggesting a role of GhFAR3.1 in wax synthesis and protection against water loss. Compared to the well‐watered conditions, drought stress induced significant accumulation of wax in leaves of wild‐type plants but not in leaves of GhFAR3.1 silenced plants, leading to less water holding capacity in GhFAR3.1 silenced plants and plant wilting. Silencing GhFAR3.1 had no effect on the expression levels of the wax biosynthesis pathway genes KAS, KCS, and LACS (upstream GhFAR3.1), but reduced the transcript level of its downstream gene WSD. Together, these results suggest that leaf wax content is important for water retention and drought tolerance and that GhFAR3.1 is essential for wax synthesis in cotton leaves. These results also provide the basis for further study on the molecular regulation mechanism of GhFARs in cotton development and surface lipid synthesis.

Upland cotton (Gossypium hirsutum L.) is one of the important fiber crops. Cotton flowers usually appear white (or cream-colored) without colored spots at the petal base, and turn pink on the next day after flowering. In this study, using a mutant showing pink petals with crimson spots at their base, we conducted comparative metabolome and transcriptome analyses to investigate the molecular mechanism of coloration in cotton flowers. Metabolic profiling showed that cyanidin-3-O-glucoside and glycosidic derivatives of pelargonidins and peonidins are the main pigments responsible for the coloration of the pink petals of the mutant. A total of 2443 genes differentially expressed (DEGs) between the white and pink petals were identified by RNA-sequencing. Many DEGs are structural genes and regulatory genes of the anthocyanin biosynthesis pathway. Among them, MYB21, UGT88F3, GSTF12, and VPS32.3 showed significant association with the accumulation of cyanidin-3-O-glucoside in the pink petals. Taken together, our study preliminarily revealed the metabolites responsible for the pink petals and the key genes regulating the biosynthesis and accumulation of anthocyanins in the pink petals. The results provide new insights into the biochemical and molecular mechanism underlying anthocyanin biosynthesis in upland cotton.