A genetic and metabolic analysis revealed that cotton fiber cell development was retarded by flavonoid naringenin

GhLCYε-3 characterized as a lycopene cyclase gene responding to drought stress in cotton

Drought stress significantly affects crop productivity. Carotenoids are essential photosynthetic pigment for plants, bacteria, and algae, with signaling and antioxidant functions. Lutein is a crucial branch product in the carotenoid synthesis pathway, which effectively improves the stress tolerance of higher plants. lycopene cyclase, a central enzyme for lutein synthesis, holds great significance in regulating lutein production. This research establishes a correlation between lutein content and stress resistance by measuring the drought resistance and lutein content of various cotton materials. To identify which crucial genes are associated with lutein, the lycopene cyclase family (LCYs) was analyzed. The research found that LCYs form a highly conserved family divided into two subfamilies, LCY-ε (lycopene ε-cyclase) and LCY-β (lycopene β-cyclase). Most members of the LCY family contain photoresponsive elements and abscisic acid elements. qRT-PCR demonstrates showed that most genes responded positively to drought stress, and GhLCYε-3 was expressed significantly differently under drought stress. Virus-induced gene silencing (VIGS) assay showed that the content of GhLCYε-3 was significantly increased with MDA and PRO, and the contents of chlorophyll and lutein were significantly decreased in pYL156 plants. The decrease in GhLCYε-3 expression is speculated to lead to reduced lutein content in vivo, resulting in the accumulation of reactive oxygen species (ROS) and decreased drought tolerance. This research enriched the understanding of LCY gene family and lutein function, and provided a new reference for cotton planting in arid areas.

Synopsis

Lycopene cyclase plays an important role in enhancing the ability of scavenging ROS and drought resistance of plants.

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Synergistic interplay of redox homeostasis and polysaccharide synthesis promotes cotton fiber elongation

Cell polarity is the foundation of cell development and tissue morphogenesis. The investigation of polarized growth provides opportunities to gain profound insights into morphogenesis and tissue functionality in organisms. Currently, there are still many mysteries surrounding the mechanisms that regulate polarized cell growth. Cotton fiber cells serve as an excellent model for studying polarized growth, and provide important clues for unraveling the molecular mechanisms, signaling pathways, and regulatory networks of polarized growth. In this study, we characterized two functional genes, GhMDHAR1AT/DT and GhDHAR2AT/DT with predominant expression during fiber elongation. Loss of function of both genes contributed to a significant increase in fiber length. Transcriptomic data revealed up-regulated expression of antioxidant genes in CRISPR mutant lines, along with delayed expression of secondary wall-related genes and temporally prolonged expression of primary wall-related genes. Experimental evidence demonstrated that the increase in GSH content and glutathione peroxidase (GPX) enzyme activity led to enhanced total antioxidant capacity (T-AOC), resulting in reduced H2O2 levels, which contributed to the extension of fiber elongation stage in CRISPR mutant lines. Moreover, the increased polysaccharide synthesis in CRISPR mutant lines was found to provide an abundant supply of raw materials for fiber cell wall elongation, suggesting that synergistic interplay between redox homeostasis and polysaccharide synthesis in fiber cells may facilitate cell wall remodeling and fiber elongation. This study provides valuable insights for deciphering the mechanisms of cell polarized growth and improving cotton fiber quality.

Naturally colored cotton for wearable applications

Naturally colored cotton (NCC) offers an environmentally friendly fiber for textile applications. Processing white cotton fiber into textiles requires extensive energy, water, and chemicals, whereas processing of NCC skips the most polluting activity, scouring-bleaching and dyeing; therefore, NCC provides an avenue to minimize the harmful impacts of textile production. NCC varieties are suitable for organic agriculture since they are naturally insect and disease-resistant, salt and drought-tolerant. Various fiber shades, ranging from light green to tan and brown, are available in the cultivated NCC (Gossypium hirsutum L.) species. The pigments responsible for the color of brown cotton fiber are proanthocyanidins or their derivatives synthesized by the flavonoid pathway. Due to pigments, the NCC has excellent ultraviolet protection properties. Some brown cotton varieties exhibited superior thermal resistance of fiber that can be used to make fabrics with enhanced flame retardancy. Here, we review molecular mechanisms involved in the pigment production of brown cotton and challenges in breeding NCC varieties with a wide range of colors but without penalty in fiber quality. Also, we discuss opportunities for NCC with flame-retarding properties in textile applications.

The WRKY46-MYC2 module plays a critical role in E-2-hexenal-induced anti-herbivore responses by promoting flavonoid accumulation

Volatile organic compounds (VOCs) play key roles in plant-plant communication, especially in response to pest attack. E-2-hexenal is an important component of VOCs, but it is unclear whether it can induce endogenous plant resistance to insects. Here, we show that E-2-hexenal activates early signaling events in Arabidopsis (Arabidopsis thaliana) mesophyll cells, including an H2O2 burst at the plasma membrane, the directed flow of calcium ions, and an increase in cytosolic calcium concentration. Treatment of wild-type Arabidopsis plants with E-2-hexenal increases their resistance when challenged with the diamondback moth Plutella xylostella L., and this phenomenon is largely lost in the wrky46 mutant. Mechanistically, E-2-hexenal induces the expression of WRKY46 and MYC2, and the physical interaction of their encoded proteins was verified by yeast two-hybrid, firefly luciferase complementation imaging, and in vitro pull-down assays. The WRKY46-MYC2 complex directly binds to the promoter of RBOHD to promote its expression, as demonstrated by luciferase reporter, yeast one-hybrid, chromatin immunoprecipitation, and electrophoretic mobility shift assays. This module also positively regulates the expression of E-2-hexenal-induced naringenin biosynthesis genes (TT4 and CHIL) and the accumulation of total flavonoids, thereby modulating plant tolerance to insects. Together, our results highlight an important role for the WRKY46-MYC2 module in the E-2-hexenal-induced defense response of Arabidopsis, providing new insights into the mechanisms by which VOCs trigger plant defense responses.

Melatonin receptor, GhCAND2-D5 motivated responding to NaCl signaling in cotton

As a receptor for plant melatonin, CAND2/PMTR plays an important role in melatonin signaling. Most of the CANDs are membrane proteins and play indispensable roles in signal transduction. In this study, the CANDs from four cotton species were characterized, and the phylogenetic relationships, expression patterns, stress responses of cotton CANDs were analyzed by bioinformatics. Through the analysis of phylogenetic and protein structure, it was found that the CANDs in clade Ⅱ might function as cotton melatonin receptors, and most of the GhCANDs in clade Ⅱ were induced by melatonin. A putative cotton melatonin receptor, GhCAND2-D5, was functionally probed by gene silencing. The plants with silenced expression of this gene exhibited decreased salt tolerance. Protein interaction prediction identified that GhCAND2-D5 interacted with several membrane proteins and played an important role in melatonin signaling. This study provided a theoretical reference for further investigation of melatonin signaling in cotton.

Morphological and transcriptomic responses/acclimations of erect-type submerged macrophyte Hydrilla verticillata both at low-light exposure and light recovery phases

Light intensity is a determinant for submerged macrophytes. Little is known about their molecular responses to low-light exposure, despite more informative and responsive than morphological traits. For erect-type submerged macrophytes, the stem is more crucial relative to the leaf in acclimation to low-light stress, but receives less attention. We determined morphological and stem transcriptomic responses/acclimations of Hydrilla verticillata to extremely and mildly low light (7.2 and 36 μmol photons m-2 s-1, respectively), that is, EL and ML, with the radiation intensity of 180 μmol photons m-2 s-1 as the control. Low-light exposure continued for 9 days, followed by a 7-day recovery phase (180 μmol photons m-2 s-1). At the exposure phase, the low-light treatments, in particular the EL, decreased the relative growth ratio, but induced greater height and longer stem internode distance and epidermal cell. Such responses/acclimations continued into the recovery phase, despite more or less changes in the magnitude. Transcriptome showed that the photosynthetic system was inhibited at the exposure phase, but the macrophyte adjusted hormone synthesis relating to cell division and elongation. Moreover, the EL activated cell stress responses such as DNA repair. Following light recovery, the macrophyte exhibited a strong-light response, although energy metabolism enhanced. Especially, the EL enriched the pathways relating to anthocyanin synthesis at such phase, indicating an activation of photoprotective mechanism. Our findings suggest that negative influences of low light occur at both low-light exposure and recovery phases, but submerged macrophytes would acclimate to light environments. Transcriptome can show molecular basis of plant responses/acclimations, including but not limited to morphology. This study establishes a bridge connecting morphological and molecular responses/acclimations.

Integrated Transcriptome and Metabolome Analysis of Salinity Tolerance in Response to Foliar Application of β-Alanine in Cotton Seedlings

Salinity is amongst the serious abiotic stresses cotton plants face, impairing crop productivity. Foliar application of β-alanine is employed to improve salt tolerance in various crops, but the exact mechanism behind it is not yet completely understood. An advanced line SDS-01 of upland cotton Gossypium hirsutum L. was utilized to determine its salt tolerance. Foliar treatment with the β-alanine solution at different concentrations was applied to the seedlings stressed with 0.8% NaCl solution. On the 10th day of treatment, samples were collected for transcriptome and metabolome analyses. β-alanine solution at a concentration of 25 mM was found to be the best treatment with the lowest mortality rate and highest plant height and above-ground biomass under salt stress. Both differentially expressed genes and accumulated metabolites analyses showed improved tolerance of treated seedlings. The photosynthetic efficiency improved in seedlings due to higher expression of photosynthesis-antenna proteins and activation of hormones signal transduction after treatment with β-alanine. Highly expressed transcription factors observed were MYB, HD-ZIP, ARF, MYC, EREB, DELLA, ABF, H2A, H4, WRKY, and HK involved in the positive regulation of salinity tolerance in β-alanine-treated seedlings. Furthermore, compared to the control, the high accumulation of polyamines, coumarins, organic acids, and phenolic compounds in the β-alanine-treated seedlings helped regulate cellular antioxidant (glutathione and L-Cysteine) production. Hence, to improve salt tolerance and productivity in cotton, foliar application of β-alanine at the seedling stage can be a valuable management practice.

Profile of cotton flavonoids: Their composition and important roles in development and adaptation to adverse environments

Cotton is a commercial crop that is cultivated in more than 50 countries. The production of cotton has severely diminished in recent years owing to adverse environments. Thus, it is a high priority of the cotton industry to produce resistant cultivars to prevent diminished cotton yields and quality. Flavonoids comprise one of the most important groups of phenolic metabolites in plants. However, the advantage and biological roles of flavonoids in cotton have yet not been studied in depth. In this study, we performed a widely targeted metabolic study and identified 190 flavonoids in cotton leaves that span seven different classes with flavones and flavonols as the dominant groups. Furthermore, flavanone-3-hydroxylase was cloned and silenced to knock down flavonoid production. The results show that the inhibition of flavonoid biosynthesis affects the growth and development of cotton and causes semi-dwarfing in cotton seedlings. We also revealed that the flavonoids contribute to cotton defense against ultraviolet radiation and Verticillium dahliae. Moreover, we discuss the promising role of flavonoids in cotton development and defense against biotic and abiotic stresses. This study provides valuable information to study the variety and biological functions of flavonoids in cotton and will help to profile the advantages of flavonoids in cotton breeding.

GhWRKY41 forms a positive feedback regulation loop and increases cotton defence response against Verticillium dahliae by regulating phenylpropanoid metabolism

Despite the established significance of WRKY proteins and phenylpropanoid metabolism in plant immunity, how WRKY proteins modulate aspects of the phenylpropanoid pathway remains undetermined. To understand better the role of WRKY proteins in plant defence, we identified a cotton (Gossypium hirsutum) protein, GhWRKY41, that is, universally and rapidly induced in three disease-resistant cotton cultivars following inoculation with the plant pathogenic fungus, Verticillium dahliae. We show that overexpression of GhWRKY41 in transgenic cotton and Arabidopsis enhances resistance to V. dahliae, while knock-down increases cotton more susceptibility to the fungus. GhWRKY41 physically interacts with itself and directly activates its own transcription. A genome-wide chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq), in combination with RNA sequencing (RNA-seq) analyses, revealed that 43.1% of GhWRKY41-binding genes were up-regulated in cotton upon inoculation with V. dahliae, including several phenylpropanoid metabolism master switches, receptor kinases, and disease resistance-related proteins. We also show that GhWRKY41 homodimer directly activates the expression of GhC4H and Gh4CL, thereby modulating the accumulation of lignin and flavonoids. This finding expands our understanding of WRKY-WRKY protein interactions and provides important insights into the regulation of the phenylpropanoid pathway in plant immune responses by a WRKY protein.

Coordination of floral and fiber development in cotton (Gossypium) by hormone- and flavonoid-signalling associated regulatory miRNAs

Various plant development activities and stress responses are tightly regulated by various microRNAs (miRNA) and their target genes, or transcription factors in a spatiotemporal manner. Here, to exemplify how flowering-associated regulatory miRNAs synchronize their expression dynamics during floral and fiber development in cotton, constitutive expression diminution transgenic lines of auxin-signaling regulatory Gh-miR167 (35S-MIM167) were developed through target mimicry approach. 'Moderate' (58% to 80%)- and 'high' (> 80%)-Gh-miR167 diminution mimic lines showed dosage-dependent developmental deformities in anther development, pollen maturation, and fruit (= boll) formation. Cross pollination of 'moderate' 35S-MIM167 mimic lines with wild type (WT) plant partially restored boll formation and emergence of fiber initials on the ovule surface. Gh-miR167 diminution favored organ-specific transcription biases in miR159, miR166 as well as miR160, miR164, and miR172 along with their target genes during anther and petal development, respectively. Similarly, accumulative effect of percent Gh-miR167 diminution, cross regulation of its target ARF6/8 genes, and temporal mis-expression of hormone signaling- and flavonoid biosynthesis-associated regulatory miRNAs at early fiber initiation stage caused irregular fiber formation. Spatial and temporal transcription proportions of regulatory miRNAs were also found crucial for the execution of hormone- and flavonoid-dependent progression of floral and fiber development. These observations discover how assorted regulatory genetic circuits get organized in response to Gh-miR167 diminution and converge upon ensuing episodes of floral and fiber development in cotton.

Cloning and functional analysis of GhDFR1, a key gene of flavonoid synthesis pathway in naturally colored cotton

BACKGROUND

The naturally colored brown cotton fiber is the most widely used environmentally friendly textile material, which primarily contains proanthocyanidins and their derivatives. Many structural genes in the flavonoid synthesis pathway are known to improve the genetic resources of naturally colored cotton. Among them, DFR is a crucial late enzyme to synthesis both anthocyanins and proanthocyanidins in the plant flavonoid pathway.

METHODS

The protein sequences of GhDFRs were analyzed using bioinformatic tools. The expression levels of GhDFRs in various tissues and organs of upland cotton Zongxu1 (ZX1), were analyzed by quantitative real-time PCR, and the expression pattern of GhDFR1 during fiber development of white cotton and brown cotton was analyzed further. The function of GhDFR1 in NCC ZX1 was preliminarily analyzed by virus induced gene silencing (VIGS) technology.

RESULTS

Bioinformatic analysis revealed that GhDFRs sequences in upland cotton genome were extremely conserved. Furthermore, evolutionary tree analysis revealed that the functions of GhDFR1 and GhDFR2, and GhDFR3 and GhDFR4, presented different and shared some similarities. Our study showed GhDFR1 and GhDFR2 were specifically expressed in fibers, while GhDFR3 and GhDFR4 were specifically expressed in petals. GhDFR1 was exclusively expressed in brown cotton fiber at various stages of development and progressively increased with the growth of fiber, but the trend of expression in white cotton was quite the opposite. We silenced GhDFR1 expression in brown cotton fiber using VIGS technology, and observed the VIGS-interference plants. After reducing the expression level of GhDFR1, the period for significant GhDFR1 expression in the developing fibers changed, reducing the content of anthocyanins, and lightening the color of mature cotton fibers.

CONCLUSION

GhDFR1 was preferentially expressed in brown cotton during fiber development. The timing of GhDFR1 expression for flavonoid synthesis altered, resulting in anthocyanin contents reduced and the fiber color of the GhDFR1i lines lightened. These findings showed the role of GhDFR1 in fiber coloration of NCC and provided a new candidate for NCC genetic improvement.

Cotton 4-coumarate-CoA ligase 3 enhanced plant resistance to Verticillium dahliae by promoting jasmonic acid signaling-mediated vascular lignification and metabolic flux

Lignins and their antimicrobial-related polymers cooperatively enhance plant resistance to pathogens. Several isoforms of 4-coumarate-coenzyme A ligases (4CLs) have been identified as indispensable enzymes involved in lignin and flavonoid biosynthetic pathways. However, their roles in plant-pathogen interaction are still poorly understood. This study uncovers the role of Gh4CL3 in cotton resistance to the vascular pathogen Verticillium dahliae. The cotton 4CL3-CRISPR/Cas9 mutant (CR4cl) exhibited high susceptibility to V. dahliae. This susceptibility was most probably due to the reduction in the total lignin content and the biosynthesis of several phenolic metabolites, e.g., rutin, catechin, scopoletin glucoside, and chlorogenic acid, along with jasmonic acid (JA) attenuation. These changes were coupled with a significant reduction in 4CL activity toward p-coumaric acid substrate, and it is likely that recombinant Gh4CL3 could specifically catalyze p-coumaric acid to form p-coumaroyl-coenzyme A. Thus, overexpression of Gh4CL3 (OE4CL) showed increasing 4CL activity that augmented phenolic precursors, cinnamic, p-coumaric, and sinapic acids, channeling into lignin and flavonoid biosyntheses and enhanced resistance to V. dahliae. Besides, Gh4CL3 overexpression activated JA signaling that instantly stimulated lignin deposition and metabolic flux in response to pathogen, which all established an efficient plant defense response system, and inhibited V. dahliae mycelium growth. Our results propose that Gh4CL3 acts as a positive regulator for cotton resistance against V. dahliae by promoting JA signaling-mediated enhanced cell wall rigidity and metabolic flux.

A glutathione S-transferase GhTT19 determines flower petal pigmentation via regulating anthocyanin accumulation in cotton

Anthocyanin accumulations in the flowers can improve seed production of hybrid lines, and produce higher commodity value in cotton fibre. However, the genetic mechanism underlying the anthocyanin pigmentation in cotton petals is poorly understood. Here, we showed that the red petal phenotype was introgressed from Gossypium bickii through recombination with the segment containing the R₃ ^bic region in the A07 chromosome of Gossypium hirsutum variety LR compared with the near-isogenic line of LW with white flower petals. The cyanidin-3-O-glucoside (Cy3G) was the major anthocyanin in red petals of cotton. A GhTT19 encoding a TT19-like GST was mapped to the R₃ ^bic site associated with red petals via map-based cloning, but GhTT19 homologue gene from the D genome was not expressed in G. hirsutum. Intriguingly, allelic variations in the promoters between GhTT19^LW and GhTT19^LR , rather than genic regions, were found as genetic causal of petal colour variations. GhTT19-GFP was found localized in both the endoplasmic reticulum and tonoplast for facilitating anthocyanin transport. An additional MYB binding element found only in the promoter of GhTT19^LR , but not in that of GhTT19^LW , enhanced its transactivation by the MYB activator GhPAP1. The transgenic analysis confirmed the function of GhTT19 in regulating the red flower phenotype in cotton. The essential light signalling component GhHY5 bonded to and activated the promoter of GhPAP1, and the GhHY5-GhPAP1 module together regulated GhTT19 expression to mediate the light-activation of petal anthocyanin pigmentation in cotton. This study provides new insights into the molecular mechanisms for anthocyanin accumulation and may lay a foundation for faster genetic improvement of cotton.

Comparative proteomic analysis identified proteins and the phenylpropanoid biosynthesis pathway involved in the response to ABA treatment in cotton fiber development

Abscisic acid (ABA) is a plant hormone that plays an important role in cotton fiber development. In this study, the physiological changes and proteomic profiles of cotton (Gossypium hirsutum) ovules were analyzed after 20 days of ABA or ABA inhibitor (ABAI) treatment. The results showed that compared to the control (CK), the fiber length was significantly decreased under ABA treatment and increased under ABAI treatment. Using a tandem mass tags-based quantitative technique, the proteomes of cotton ovules were comprehensively analyzed. A total of 7321 proteins were identified, of which 365 and 69 differentially accumulated proteins (DAPs) were identified in ABA versus CK and ABAI versus CK, respectively. Specifically, 345 and 20 DAPs were up- and down-regulated in the ABA group, and 65 and 4 DAPs were up- and down-regulated in the ABAI group, respectively. The DAPs in the ABA group were mainly enriched in the biosynthesis of secondary metabolites, phenylpropanoid biosynthesis and flavonoid secondary metabolism, whereas the DAPs in the ABAI group were mainly enriched in the indole alkaloid biosynthesis and phenylpropanoid biosynthesis pathways. Moreover, 9 proteins involved in phenylpropanoid biosynthesis were upregulated after ABA treatment, suggesting that this pathway might play important roles in the response to ABA, and 3 auxin-related proteins were upregulated, indicating that auxin might participate in the regulation of fiber development under ABAI treatment.

Transcriptomic analysis reveals the key role of histone deacetylation via mediating different phytohormone signalings in fiber initiation of cotton

Background

Histone deacetylation is one of the most important epigenetic modifications and plays diverse roles in plant development. However, the detailed functions and mechanisms of histone deacetylation in fiber development of cotton are still unclear. HDAC inhibitors (HDACi) have been commonly used to study the molecular mechanism underlying histone deacetylation or to facilitate disease therapy in humans through hindering the histone deacetylase catalytic activity. Trichostatin A (TSA)—the most widely used HDACi has been extensively employed to determine the role of histone deacetylation on different developmental stages of plants.

Results

Through in vitro culture of ovules, we observed that exogenous application of TSA was able to inhibit the fiber initiation development. Subsequently, we performed a transcriptomic analysis to reveal the underlying mechanisms. The data showed that TSA treatment resulted in 4209 differentially expressed genes, which were mostly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, photosynthesis, and carbon metabolism pathways. The phytohormone signal transduction pathways harbor the most differentially expressed genes. Deeper studies showed that some genes promoting auxin, Gibberellic Acid (GA) signaling were down-regulated, while some genes facilitating Abscisic Acid (ABA) and inhibiting Jasmonic Acid (JA) signaling were up-regulated after the TSA treatments. Further analysis of plant hormone contents proved that TSA significantly promoted the accumulation of ABA, JA and GA₃.

Conclusions

Collectively, histone deacetylation can regulate some key genes involved in different phytohormone pathways, and consequently promoting the auxin, GA, and JA signaling, whereas repressing the ABA synthesis and signaling to improve the fiber cell initiation. Moreover, the genes associated with energy metabolism, phenylpropanoid, and glutathione metabolism were also regulated by histone deacetylation. The above results provided novel clues to illuminate the underlying mechanisms of epigenetic modifications as well as related different phytohormones in fiber cell differentiation, which is also very valuable for the molecular breeding of higher quality cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00840-4.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

Impact of Novel Varietal and Regional Differences on Cotton Fiber Quality Characteristics

Modernization and the global fashion market demand continuous improvements in upland cotton cultivars (Gossypium hirustum L.) to meet these improved fiber characteristics (fiber length, fiber strength, micronaire) requirements. Researchers have centered their efforts on improved fiber quality; however, the efforts are not immediately supporting the textile sector. The daily mean and temperature amplitude fluctuation affect cotton yield and fiber characteristics. This study analyzed four newly developed cotton varieties in two cotton regions for fiber characteristics’ variations. It was observed that cotton fiber quality characteristics (fiber length, uniformity, strength, and micronaire) are impacted in diverse ways. Fiber quality is mainly affected by the genotype and environmental conditions, e.g., weather conditions, irrigation management, fertilization, and cultural practices. The Khanewal region had shown better fiber characteristics than the Multan region, whereas cotton variety CIM-785 had better fiber characteristics in both regions.

Transcriptome Co-expression Network and Metabolome Analysis Identifies Key Genes and Regulators of Proanthocyanidins Biosynthesis in Brown Cotton

Brown cotton fiber (BCF) is a unique raw material of naturally colored cotton (NCC). But characteristics of the regulatory gene network and metabolic components related to the proanthocyanidins biosynthesis pathway at various stages of its fiber development remain unclear. Here, the dynamic changes in proanthocyanidins biosynthesis components and transcripts in the BCF variety "Zong 1-61" and its white near-isogenic lines (NILs) "RT" were characterized at five fiber developmental stages (0, 5, 10, 15, and 20 days post-anthesis; DPA). Enrichment analysis of differentially expressed genes (DEGs), comparison of metabolome differences, and pathway enrichment analysis of a weighted gene correlation network analysis together revealed the dominant gene expression of flavonoid biosynthesis (FB), phenylpropanoid metabolisms, and some carbohydrate metabolisms at 15 or 20 DPA than white cotton. Eventually, 63 genes were identified from five modules putatively related to FB. Three R2R3-MYB and two bHLH transcription factors were predicted as the core genes. Further, GhANS, GhANR1, and GhUFGT2 were preliminarily regulated by GhMYB46, GhMYB6, and GhMYB3, respectively, according to yeast one-hybrid assays in vitro. Our findings provide an important transcriptional regulatory network of proanthocyanidins biosynthesis pathway and dynamic flavonoid metabolism profiles.

Dynamic Expression, Differential Regulation and Functional Diversity of the CNGC Family Genes in Cotton

Cyclic nucleotide-gated channels (CNGCs) constitute a family of non-selective cation channels that are primarily permeable to Ca2+ and activated by the direct binding of cyclic nucleotides (i.e., cAMP and cGMP) to mediate cellular signaling, both in animals and plants. Until now, our understanding of CNGCs in cotton (Gossypium spp.) remains poorly addressed. In the present study, we have identified 40, 41, 20, 20, and 20 CNGC genes in G. hirsutum, G. barbadense, G. herbaceum, G. arboreum, and G. raimondii, respectively, and demonstrated characteristics of the phylogenetic relationships, gene structures, chromosomal localization, gene duplication, and synteny. Further investigation of CNGC genes in G. hirsutum, named GhCNGC1-40, indicated that they are not only extensively expressed in various tissues and at different developmental stages, but also display diverse expression patterns in response to hormones (abscisic acid, salicylic acid, methyl jasmonate, ethylene), abiotic (salt stress) and biotic (Verticillium dahlia infection) stimuli, which conform with a variety of cis-acting regulatory elements residing in the promoter regions; moreover, a set of GhCNGCs are responsive to cAMP signaling during cotton fiber development. Protein–protein interactions supported the functional aspects of GhCNGCs in plant growth, development, and stress responses. Accordingly, the silencing of the homoeologous gene pair GhCNGC1&18 and GhCNGC12&31 impaired plant growth and development; however, GhCNGC1&18-silenced plants enhanced Verticillium wilt resistance and salt tolerance, whereas GhCNGC12&31-silenced plants had opposite effects. Together, these results unveiled the dynamic expression, differential regulation, and functional diversity of the CNGC family genes in cotton. The present work has laid the foundation for further studies and the utilization of CNGCs in cotton genetic improvement.

Identification of SNAT Family Genes Suggests GhSNAT3D Functional Reponse to Melatonin Synthesis Under Salinity Stress in Cotton

Serotonin N-acetyltransferase (SNAT) is a key enzyme in the biosynthesis of melatonin, and plays an important role in the regulation of melatonin synthesis. The study of SNAT is of great significance to understand the function of melatonin. In this study, we analyzed the structural characteristics, phylogenetic relationship, gene structure, expression pattern, evolutionary relationship and stress response of the members of the SNAT gene family in upland cotton through bioinformatics. A putative Serotonin n-acetyltransferase gene GhSNAT3D was identified, and preliminarily function of GhSNAT3D was verified by virus-induced gene silencing. We identified a total of 52 SNAT genes in the whole genome of G. hirsutum, and part of the GhSNATs were regulated by exogenous melatonin. The content of melatonin, antioxidant enzyme activity and Ca²⁺ content of GhSNAT3D gene silenced plants decreased, and the salt tolerance of GhSNAT3D gene silenced plants was reduced. Exogenous melatonin supplementation restored the salt tolerance of GhSNAT3D gene silenced plants. GhSNAT3D may interact with GhSNAT25D and ASMT to regulate melatonin synthesis. This study provided an important basis for further study on the regulation of melatonin in cotton against abiotic stress.

Dynamic 3D genome architecture of cotton fiber reveals subgenome-coordinated chromatin topology for 4-staged single-cell differentiation

Background

Despite remarkable advances in our knowledge of epigenetically mediated transcriptional programming of cell differentiation in plants, little is known about chromatin topology and its functional implications in this process.

Results

To interrogate its significance, we establish the dynamic three-dimensional (3D) genome architecture of the allotetraploid cotton fiber, representing a typical single cell undergoing staged development in plants. We show that the subgenome-relayed switching of the chromatin compartment from active to inactive is coupled with the silencing of developmentally repressed genes, pinpointing subgenome-coordinated contribution to fiber development. We identify 10,571 topologically associating domain-like (TAD-like) structures, of which 25.6% are specifically organized in different stages and 75.23% are subject to partition or fusion between two subgenomes. Notably, dissolution of intricate TAD-like structure cliques showing long-range interactions represents a prominent characteristic at the later developmental stage. Dynamic chromatin loops are found to mediate the rewiring of gene regulatory networks that exhibit a significant difference between the two subgenomes, implicating expression bias of homologous genes.

Conclusions

This study sheds light on the spatial-temporal asymmetric chromatin structures of two subgenomes in the cotton fiber and offers a new insight into the regulatory orchestration of cell differentiation in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02616-y.

Global expression dynamics and miRNA evolution profile govern floral/fiber architecture in the modern cotton (Gossypium)

Majority of differentially expressed miRNAs with functional attributes have been recruited independently and parallelly during allopolyploidy followed by the millennia of human selection of both domesticated G. hirsutum and G. barbadense. The genus Gossypium is a marvelous evolutionary model for studying allopolyploidy and morpho-evolution of long-spinnable fibers from the ancestral wild-fuzz. Many genes, transcription factors, and notably, the regulatory miRNAs essentially govern such remarkable modern fiber phenotypes. To comprehend the impact of allopolyploidy on the evolutionary selection of transcriptional dynamicity of key miRNAs, comparative transcriptome profiling of vegetative and fiber tissues of domesticated diploid G. arboreum (A₂) and allopolyploid cotton species G. hirsutum (AD₁), and G. barbadense (AD₂) identified > 300 differentially expressed miRNAs (DEmiRs) within or between corresponding tissues of A₂, AD₁ and AD₂ species. Up to 49% and 32% DEmiRs were up- and down-regulated at fiber initiation stage of AD₁ and AD₂ species, respectively, whereas 50% and 18% DEmiRs were up- and down-regulated at fiber elongation stage of both the allopolyploid species. Interestingly, A-subgenome-specific DEmiRs exhibit expression dominance in the allopolyploid genetic backgrounds. Comparative spatio-temporal expression analyses of AD₁ and AD₂ species discovered that a majority of DEmiRs were recruited independently under millennia of human selection during domestication. Functional annotations of these DEmiRs revealed selection of associated molecular functions such as hormone-signaling, calcium-signaling and reactive oxygen species (ROS) signaling during fiber initiation and elongation. To validate the functional attributes of annotated DEmiRs, we demonstrated for the first time that the target-mimicry-based constitutive diminution of auxin-signaling associated miR167 directly affected the differentiation of floral and fiber tissues of transgenic cotton. These results strongly suggested that the evolutionarily favored DEmiRs including miR167 were involved in the transcriptional regulation of numerous genes during cotton evolution for enhanced fiber-associated agronomic traits.

Melatonin Improves Cotton Salt Tolerance by Regulating ROS Scavenging System and Ca2 + Signal Transduction

As one of the cash crops, cotton is facing the threat of abiotic stress during its growth and development. It has been reported that melatonin is involved in plant defense against salt stress, but whether melatonin can improve cotton salt tolerance and its molecular mechanism remain unclear. We investigated the role of melatonin in cotton salt tolerance by silencing melatonin synthesis gene and exogenous melatonin application in upland cotton. In this study, applicating of melatonin can improve salt tolerance of cotton seedlings. The content of endogenous melatonin was different in cotton varieties with different salt tolerance. The inhibition of melatonin biosynthesis related genes and endogenous melatonin content in cotton resulted in the decrease of antioxidant enzyme activity, Ca²⁺ content and salt tolerance of cotton. To explore the protective mechanism of exogenous melatonin against salt stress by RNA-seq analysis. Melatonin played an important role in the resistance of cotton to salt stress, improved the salt tolerance of cotton by regulating antioxidant enzymes, transcription factors, plant hormones, signal molecules and Ca²⁺ signal transduction. This study proposed a regulatory network for melatonin to regulate cotton's response to salt stress, which provided a theoretical basis for improving cotton's salt tolerance.

Breeding Next-Generation Naturally Colored Cotton

Use of naturally colored cotton (NCC) in the textile industry is limited due to its inferior fiber quality. Recent advances in understanding pigmentation in NCC fiber and the genetic relationship between pigmentation and fiber development may overcome this hindrance and can guide developing high fiber quality and yield NCC.

GhMYB4 downregulates lignin biosynthesis and enhances cotton resistance to Verticillium dahliae

GhMYB4 acts as a negative regulator in lignin biosynthesis, which results in alteration of cell wall integrity and activation of cotton defense response. Verticillium wilt of cotton (Gossypium hirsutum) caused by the soil-borne fungus Verticillium dahliae (V. dahliae) represents one of the most important constraints of cotton production worldwide. Mining of the genes involved in disease resistance and illuminating the molecular mechanisms that underlie this resistance is of great importance in cotton breeding programs. Defense-induced lignification in plants is necessary for innate immunity, and there are reports of a correlation between increased lignification and disease resistance. In this study, we present an example in cotton whereby plants with reduced lignin content also exhibit enhanced disease resistance. We identified a negative regulator of lignin synthesis, in cotton encoded in GhMYB4. Overexpression of GhMYB4 in cotton and Arabidopsis enhanced resistance to V. dahliae  with reduced lignin deposition. Moreover, GhMYB4 could bind the promoters of several genes involved in lignin synthesis, such as GhC4H-1, GhC4H-2, Gh4CL-4, and GhCAD-3, and impair their expression. The reduction of lignin content in GhMYB4-overexpressing cotton led to alterations of cell wall integrity (CWI) and released more oligogalacturonides (OGs) which may act as damage-associated molecular patterns (DAMPs) to stimulate plant defense responses. In support of this hypothesis, exogenous application with polygalacturonic acid (PGA) in cotton activated biosynthesis of jasmonic acid (JA) and JA-mediated defense against V. dahliae, similar to that described for cotton plants overexpressing GhMYB4. This study provides a new candidate gene for cotton disease-resistant breeding and an increased understanding of the relationship between lignin synthesis, OG release, and plant immunity.

Evolution, expression and functional analysis of cultivated allotetraploid cotton DIR genes

BACKGROUND

Dirigent (DIR) proteins mediate regioselectivity and stereoselectivity during lignan biosynthesis and are also involved in lignin, gossypol and pterocarpan biosynthesis. This gene family plays a vital role in enhancing stress resistance and in secondary cell-wall development, but systematical understanding is lacking in cotton.

RESULTS

In this study, 107 GbDIRs and 107 GhDIRs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. Most of these genes have a classical gene structure without intron and encode proteins containing a signal peptide. Phylogenetic analysis showed that cotton DIR genes were classified into four distinct subfamilies (a, b/d, e, and f). Of these groups, DIR-a and DIR-e were evolutionarily conserved, and segmental and tandem duplications contributed equally to their formation. In contrast, DIR-b/d mainly expanded by recent tandem duplications, accompanying with a number of gene clusters. With the rapid evolution, DIR-b/d-III was a Gossypium-specific clade involved in atropselective synthesis of gossypol. RNA-seq data highlighted GhDIRs in response to Verticillium dahliae infection and suggested that DIR gene family could confer Verticillium wilt resistance. We also identified candidate DIR genes related to fiber development in G. barbadense and G. hirsutum and revealed their differential expression. To further determine the involvement of DIR genes in fiber development, we overexpressed a fiber length-related gene GbDIR78 in Arabidopsis and validated its function in trichomes and hypocotyls.

CONCLUSIONS

These findings contribute novel insights towards the evolution of DIR gene family and provide valuable information for further understanding the roles of DIR genes in cotton fiber development as well as in stress responses.

Detection of candidate proteins in the indican biosynthetic pathway of Persicaria tinctoria (Polygonum tinctorium) using protein-protein interactions and transcriptome analyses

Persicaria tinctoria (Polygonum tinctorium) synthesizes indican (indoxyl-β-D-glucoside) as a specialized metabolite. Indican is synthesized in the cytosol of leaf cells from indoxyl and UDP-glucose by the catalysis of indoxyl-β-D-glucoside synthase (PtIGS), then transported into vacuoles. As a portion of PtIGS is found on the microsomal membrane, we assume that it is present on the ER membrane as a large complex involving other indican metabolism-related proteins. Based on this hypothesis, the existence of such a complex was investigated using two separate approaches: a protein-protein interaction assay and transcriptome analysis. We first performed a co-immunoprecipitation using the anti-PtIGS antibody and a pull-down assay using recombinant PtIGS, then identified the candidate proteins through MS/MS analysis. Secondly, we performed a transcriptome analysis to examine the differential gene expression between the first and the second leaves. The expressions of candidate genes detected by protein-protein interaction analyses were collated with transcriptome data and validated by quantitative reverse transcription polymerase chain reaction, showing that the expression of sucrose synthase and cytochrome P450 genes decreased in the second leaves compared with the first leaves. Furthermore, we detected several additional proteins, such as heat shock and cytoskeletal proteins, suggesting that PtIGS may form a large complex, a metabolon.

Characterization of floral morphoanatomy and identification of marker genes preferentially expressed during specific stages of cotton flower development

Characterization of anther and ovule developmental programs and expression analyses of stage-specific floral marker genes in Gossypium hirsutum allowed to build a comprehensive portrait of cotton flower development before fiber initiation. Gossypium hirsutum is the most important cotton species that is cultivated worldwide. Although cotton reproductive development is important for fiber production, since fiber is formed on the epidermis of mature ovules, cotton floral development remains poorly understood. Therefore, this work aims to characterize the cotton floral morphoanatomy by performing a detailed description of anther and ovule developmental programs and identifying stage-specific floral marker genes in G. hirsutum. Using light microscopy and scanning electron microscopy, we analyzed anther and ovule development during 11 stages of flower development. To better characterize the ovule development in cotton, we performed histochemical analyses to evaluate the accumulation of phenolic compounds, pectin, and sugar in ovule tissues. After identification of major hallmarks of floral development, three key stages were established in G. hirsutum floral development: in stage 1 (early-EF), sepal, petal, and stamen primordia were observed; in stage 2 (intermediate-IF), primordial ovules and anthers are present, and the differentiating archesporial cells were observed, marking the beginning of microsporogenesis; and in stage 6 (late-LF), flower buds presented initial anther tapetum degeneration and microspore were released from the tetrad, and nucellus and both inner and outer integuments are developing. We used transcriptome data of cotton EF, IF and LF stages to identify floral marker genes and evaluated their expression by real-time quantitative PCR (qPCR). Twelve marker genes were preferentially expressed in a stage-specific manner, including the putative homologs for AtLEAFY, AtAPETALA 3, AtAGAMOUS-LIKE 19 and AtMALE STERILITY 1, which are crucial for several aspects of reproductive development, such as flower organogenesis and anther and petal development. We also evaluated the expression profile of B-class MADS-box genes in G. hirsutum floral transcriptome (EF, IF, and LF). In addition, we performed a comparative analysis of developmental programs between Arabidopsis thaliana and G. hirsutum that considered major morphoanatomical and molecular processes of flower, anther, and ovule development. Our findings provide the first detailed analysis of cotton flower development.

Quantitative metabolome and transcriptome analysis reveals complex regulatory pathway underlying photoinduced fiber color formation in cotton

As an important plant single cell model and textile application materials, poorly known about fiber color formation in cotton, which is sensitively regulated by environmental signals. Our studies underline the importance of photo signal on sensitive fiber color formation and characterize fiber color early initiation (15 DPA) and late accumulated metabolites (45 DPA) in different lighting condition. The results revealed 236 differential metabolites between control and shading, of which phenylpropanoids metabolites accounted for 20%, including uncharacterized novel metabolites and pathways. Furthermore, the early initiation specific genes respond to the absence of light are highly correlated with phenylpropanoid metabolites related to pigmentation. The current study reveals the complex pathways involving early initiation regulation and late metabolic pathways. In addition, the collection composed of uncharacterized photoinduced metabolites and early initiation signaling/regulatory genes were identified, which are important resources for understanding fiber color formation. This report provides new insight into molecular regulatory and biochemical basis underlying photoinduced fiber color formation in cotton.

Fumonisin B1-Induced Changes in Cotton Fiber Elongation Revealed by Sphingolipidomics and Proteomics

Sphingolipids are essential biomolecules and membrane components, but their regulatory role in cotton fiber development is poorly understood. Here, we found that fumonisin B1 (FB1)—a sphingolipid synthesis inhibitor—could block fiber elongation severely. Using liquid chromatography tandem mass spectrometry (LC-MS/MS), we detected 95 sphingolipids that were altered by FB1 treatment; of these, 29 (mainly simple sphingolipids) were significantly increased, while 33 (mostly complex sphingolipids) were significantly decreased. A quantitative analysis of the global proteome, using an integrated quantitative approach with tandem mass tag (TMT) labeling and LC-MS/MS, indicated the upregulation of 633 and the downregulation of 672 proteins after FB1 treatment. Most differentially expressed proteins (DEPs) were involved in processes related to phenylpropanoid and flavonoid biosynthesis. In addition, up to 20 peroxidases (POD) were found to be upregulated, and POD activity was also increased by the inhibitor. To our knowledge, this is the first report on the effects of FB1 treatment on cotton fiber and ovule sphingolipidomics and proteomics. Our findings provide target metabolites and biological pathways for cotton fiber improvement.

Understanding the role of phytohormones in cotton fiber development through omic approaches; recent advances and future directions

Cotton is among the most important fiber crops for the textile-based industry, thanks to its cellulose-rich mature fibers. The fiber initiation and elongation are one of the best models for deciphering mechanisms of single-cell differentiation and growth, that also target of fiber development programs. During the last couple of decades, high yielding omics approaches (genomics, transcriptomics, and proteomics), have helped in the identification of several genes and gene products involved in fiber development along with functional relationship to phytohormones. For example, MYB transcription factor family and Sus gene family have been evidenced by controlling cotton fiber initiation. Most importantly, the biosynthesis, responses, and transporting of phytohormones is documented to participate in the initiation of cotton fibers. Herein, in this review, the reliable genetic evidence by manipulating the above genes in cotton have been summarized to describe the relationships among key phytohormones, transcription factors, proteins, and downstream fiber growth-related genes such as Sus. The effect of other important factors such as ROS, fatty acid metabolism, and actin (globular multi-functional proteins) over fiber development has also been discussed. The challenges and deficiencies in the research of cotton fiber development have been mentioned along with a future perspective to discover new crucial genes using multiple omics analysis.

Combined GWAS and eQTL analysis uncovers a genetic regulatory network orchestrating the initiation of secondary cell wall development in cotton

The cotton fibre serves as a valuable experimental system to study cell wall synthesis in plants, but our understanding of the genetic regulation of this process during fibre development remains limited. We performed a genome-wide association study (GWAS) and identified 28 genetic loci associated with fibre quality in allotetraploid cotton. To investigate the regulatory roles of these loci, we sequenced fibre transcriptomes of 251 cotton accessions and identified 15 330 expression quantitative trait loci (eQTL). Analysis of local eQTL and GWAS data prioritised 13 likely causal genes for differential fibre quality in a transcriptome-wide association study (TWAS). Characterisation of distal eQTL revealed unequal genetic regulation patterns between two subgenomes, highlighted by an eQTL hotspot (Hot216) that established a genome-wide genetic network regulating the expression of 962 genes. The primary regulatory role of Hot216, and specifically the gene encoding a KIP-related protein, was found to be the transcriptional regulation of genes responsible for cell wall synthesis, which contributes to fibre length by modulating the developmental transition from rapid cell elongation to secondary cell wall synthesis. This study uncovered the genetic regulation of fibre-cell development and revealed the molecular basis of the temporal modulation of secondary cell wall synthesis during plant cell elongation.

Map-based cloning of qBWT-c12 discovered brassinosteroid-mediated control of organ size in cotton

Assuring fiber yield stability is the primary objective for cotton breeders since the world population is on the rise, and the demand for cotton fiber is increasing every year. Thus, enhancing average cotton boll weight (BWT) could improve seed cotton production, and ultimately to increase cotton fiber yield. This study accomplished the map-based cloning of a novel boll weight regulating locus, qBWT-c12, in cotton. Bulk segregation analysis detected linked markers, aided in the detection of a stable BWT regulating locus, qBWT-c12, on Chr12 in a novel boll size mutant, BS41. Progeny evaluation confined the qBWT-c12 to a 0.89 cM interval between the AD-A12_07 and AD-FM_44 markers in recombinant derived F₃ and F₄ populations. Homology mapping detected a 40 bp insertion-deletion (InDel) site in the AD-FM_44 clone sequence situated +341 downstream of GhBRH1_A12, which showed complete linkage to the BWT phenotype. The suppressed expression of GhBRH1_A12 suggested its putative involvement during early boll development events in BS41. Although brassinosteroid (BR) biosynthesis and signaling pathway genes were up regulated in different tissues, but the organ growth was suppressed leading to dwarf plants, smaller leaves, and de-morphed smaller bolls in BS41. Thus, a disruption in the BR signal cascade is anticipated and could be related to lower GhBRH1_A12 expression in BS41.This study firstly reported the genetic dissection of boll size regulation of G. barbadense in G. hirsutum background using map-based cloning of a BWT regulating locus, qBWT-c12. Moreover, it also emphasized the putative role GhBRH1_A12 in regulating BR homeostasis and its potential to modulate plant growth and boll development in cotton.

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

BACKGROUND

Sulfotransferases (SOTs) (EC 2.8.2.-) play a crucial role in the sulphate conjugation reaction involved in plant growth, vigor, stress resistance and pathogen infection. SOTs in Arabidopsis have been carried out and divided into 8 groups. However, the systematic analysis and functional information of SOT family genes in cotton have rarely been reported.

RESULTS

According to the results of BLASTP and HMMER, we isolated 46, 46, 76 and 77 SOT genes in the genome G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. A total of 170 in 245 SOTs were further classified into four groups based on the orthologous relationships comparing with Arabidopsis, and tandem replication primarily contributed to the expansion of SOT gene family in G. hirsutum. Expression profiles of the GhSOT showed that most genes exhibited a high level of expression in the stem, leaf, and the initial stage of fiber development. The localization analysis indicated that GhSOT67 expressed in cytoplasm and located in stem and leaf tissue. Additionally, the expression of GhSOT67 were induced and the length of stem and leaf hairs were shortened after gene silencing mediated by Agrobacterium, compared with the blank and negative control plants.

CONCLUSIONS

Our findings indicated that SOT genes might be associated with fiber development in cotton and provided valuable information for further studies of SOT genes in Gossypium.

Melatonin enhances cotton immunity to Verticillium wilt via manipulating lignin and gossypol biosynthesis

Plants endure challenging environments in which they are constantly threatened by diverse pathogens. The soil-borne fungus Verticillium dahliae is a devastating pathogen affecting many plant species including cotton, in which it significantly reduces crop yield and fiber quality. Melatonin involvement in plant immunity to pathogens has been reported, but the mechanisms of melatonin-induced plant resistance are unclear. In this study, the role of melatonin in enhancing cotton resistance to V. dahliae was investigated. At the transcriptome level, exogenous melatonin increased the expression of genes in phenylpropanoid, mevalonate (MVA), and gossypol pathways after V. dahliae inoculation. As a result, lignin and gossypol, the products of these metabolic pathways, significantly increased. Silencing the serotonin N-acetyltransferase 1 (GhSNAT1) and caffeic acid O-methyltransferase (GhCOMT) melatonin biosynthesis genes compromised cotton resistance, with reduced lignin and gossypol levels after V. dahliae inoculation. Exogenous melatonin pre-treatment prior to V. dahliae inoculation restored the level of cotton resistance reduced by the above gene silencing effects. Melatonin levels were higher in resistant cotton cultivars than in susceptible cultivars after V. dahliae inoculation. The findings indicate that melatonin affects lignin and gossypol synthesis genes in phenylpropanoid, MVA, and gossypol pathways, thereby enhancing cotton resistance to V. dahliae.

Eriodictyol can modulate cellular auxin gradients to efficiently promote in vitro cotton fibre development

BACKGROUND

Flavonoids have essential roles in flower pigmentation, fibre development and disease resistance in cotton. Previous studies show that accumulation of naringenin in developing cotton fibres significantly affects fibre growth. This study focused on determining the effects of the flavonoids naringenin, dihydrokaempferol, dihydroquerectin and eriodictyol on fibre development in an in vitro system.

RESULTS

20 μM eriodictyol treatment produced a maximum fibre growth, in terms of fibre length and total fibre units. To gain insight into the associated transcriptional regulatory networks, RNA-seq analysis was performed on eriodictyol-treated elongated fibres, and computational analysis of differentially expressed genes revealed that carbohydrate metabolism and phytohormone signaling pathways were differentially modulated. Eriodictyol treatment also promoted the biosynthesis of quercetin and dihydroquerectin in ovules and elongating fibres through enhanced expression of genes encoding chalcone isomerase, chalcone synthase and flavanone 3-hydroxylase. In addition, auxin biosynthesis and signaling pathway genes were differentially expressed in eriodictyol-driven in vitro fibre elongation. In absence of auxin, eriodictyol predominantly enhanced fibre growth when the localized auxin gradient was disrupted by the auxin transport inhibitor, triiodobenzoic acid.

CONCLUSION

Eriodictyol was found to significantly enhance fibre development through accumulating and maintaining the temporal auxin gradient in developing unicellular cotton fibres.

Overexpression of a Novel Cytochrome P450 Promotes Flavonoid Biosynthesis and Osmotic Stress Tolerance in Transgenic Arabidopsis

Flavonoids are mainly associated with growth, development, and responses to diverse abiotic stresses in plants. A growing amount of data have demonstrated the biosynthesis of flavonoids through multienzyme complexes of which the membrane-bounded cytochrome P450 supergene family shares a crucial part. However, the explicit regulation mechanism of Cytochrome P450s related to flavonoid biosynthesis largely remains elusive. In the present study, we reported the identification of a stress-tolerant flavonoid biosynthetic CtCYP82G24 gene from Carthamus tinctorius. The transient transformation of CtCYP82G24 determined the subcellular localization to the cytosol. Heterologously expressed CtCYP82G24 was effective to catalyze the substrate-specific conversion, promoting the de novo biosynthesis of flavonoids in vitro. Furthermore, a qRT-PCR assay and the accumulation of metabolites demonstrated that the expression of CtCYP82G24 was effectively induced by Polyethylene glycol stress in transgenic Arabidopsis. In addition, the overexpression of CtCYP82G24 could also trigger expression levels of several other flavonoid biosynthetic genes in transgenic plants. Taken together, our findings suggest that CtCYP82G24 overexpression plays a decisive regulatory role in PEG-induced osmotic stress tolerance and alleviates flavonoid accumulation in transgenic Arabidopsis.

Transcriptome Profiling of the Elongating Internode of Cotton (Gossypium hirsutum L.) Seedlings in Response to Mepiquat Chloride

The plant growth retardant mepiquat chloride (MC) has been extensively used to produce compact plant canopies and increase yield in cotton (Gossypium hirsutum L.). Previous studies showed that MC reduced plant height and internode length by inhibiting GA biosynthesis and cell elongation. However, whether there are other molecular mechanisms underlying MC-induced growth retardation has remained largely unknown. In the present study, we conducted histological, transcriptomic, and phytohormone analyses of the second elongating internodes of cotton seedlings treated with MC. Histological analysis revealed that the MC shortened the internodes through suppressing both cell division and cell elongation. Consistent with the observed phenotype, many genes related to cell growth were significantly downregulated by MC. Transcriptome profiling showed that the expression of genes related not only to GA, but also to auxin, brassinosteroid (BR), and ethylene metabolism and signaling was remarkably suppressed, whereas that of genes related to cytokinin (CK) and abscisic acid (ABA) metabolism was induced by MC. Consistent with the expression pattern, significant decrease of endogenous GA, auxin, and BR content, but an increase in CK content was observed after MC treatment. Most of these hormone related genes displayed opposite regulation pattern by exogenous GA₃ treatment compared to MC; simultaneous application of MC and GA₃ could alleviate the genes expression changes induced by MC treatment, indicating MC does not directly affect other plant hormones, but through the inhibition of the GA biosynthesis. In addition, the expression of genes related to secondary metabolism and many transcription factors (TFs) were differentially regulated by MC. In summary, we confirmed the important role of GA in MC-induced growth inhibition of cotton, and further found that other hormones were also involved in this process in a GA-dependent manner. This study provides novel insights into the molecular mechanism underlying the MC-mediated inhibition of internode elongation in cotton seedlings.

GbSOBIR1 confers Verticillium wilt resistance by phosphorylating the transcriptional factor GbbHLH171 in Gossypium barbadense

Receptor-like kinases (RLKs) are important components of plant innate immunity. Although recent studies have revealed that the RLK suppressor of BIR1-1 (SOBIR1) can interact with multiple receptor-like proteins and is required for resistance against fungal pathogens, how the signal is transduced and triggers immune responses remains enigmatic. In this study, we identified a defence-related RLK from Gossypium barbadense (designated GbSOBIR1) and investigated its functional mechanism. Expression of the GbSOBIR1 gene is ubiquitous in cotton plants and is induced by Verticillium dahliae inoculation. Knock-down of GbSOBIR1 by virus-induced gene silencing resulted in attenuated resistance of cotton plants to V. dahliae, while heterologous overexpression of GbSOBIR1 in Arabidopsis improves resistance. We also found that the kinase region of GbSOBIR1 interacts with a basic helix-loop-helix (bHLH) transcription factor identified as GbbHLH171 in a yeast-two-hybrid screen. GbbHLH171 could interact with and be phosphorylated by GbSOBIR1 in vitro and in vivo and contributes positively to the resistance of cotton against V. dahliae. Furthermore, we found that this phosphorylation is essential to the transcriptional activity and functional role of GbbHLH171. We also show by spectrometric analysis and site-directed mutagenesis that Ser413 is the GbSOBIR1-mediated phosphorylation site of GbbHLH171. These results demonstrate that GbSOBIR1 interacts with GbbHLH171 and plays a critical role in cotton resistance to V. dahliae.

Functional Differentiation of Duplicated Flavonoid 3-O-Glycosyltransferases in the Flavonol and Anthocyanin Biosynthesis of Freesia hybrida

Flavonols and anthocyanins are two widely distributed groups of flavonoids that occurred apart during plant evolution and biosynthesized by shared specific enzymes involved in flavonoid metabolism. UDP-glucose, flavonoid 3-O-glycosyltransferase (UF3GT), is one of the common enzymes which could catalyze the glycosylation of both flavonol and anthocyanidin aglycons simultaneously in vitro. However, whether and how UF3GT paralogous genes function diversely at the biochemical and transcriptional levels are largely unknown. Recently, Fh3GT1 was identified to be a member of UF3GTs in Freesia hybrida. However, its expression patterns and enzymatic characteristics could not coincide well with flavonol accumulation. In an attempt to characterize other flavonoids, especially flavonol glycosyltransferase genes in Freesia, three closest candidate UFGT genes-Fh3GT2, Fh3GT3, and Fh3GT4-were mined from the Freesia transcriptomic database and isolated from the flowers of the widely distributed Freesia cultivar, Red River^®. Based on bioinformatic analysis and enzymatic assays, Fh3GT2 turned out to be another bona fide glycosyltransferase gene. Biochemical analysis further proved that Fh3GT2 preferentially glucosylated kaempferol while Fh3GT1 controlled the glucosylation of quercetin and anthocyanidins. In addition, transfection assays demonstrated that Fh3GT2 could be mainly activated by the flavonol regulator FhMYBF1 or the anthocyanin regulator FhPAP1, whereas Fh3GT1 could only be activated by FhPAP1. These findings suggested that Fh3GTs might have functionally diverged in flavonoid biosynthesis at both the biochemical and transcriptional levels.

Up-regulation of GhTT2-3A in cotton fibres during secondary wall thickening results in brown fibres with improved quality

Brown cotton fibres are the most widely used naturally coloured raw materials for the eco-friendly textile industry. Previous studies have indicated that brown fibre pigments belong to proanthocyanidins (PAs) or their derivatives, and fibre coloration is negatively associated with cotton productivity and fibre quality. To date, the molecular basis controlling the biosynthesis and accumulation of brown pigments in cotton fibres is largely unknown. In this study, based on expressional and transgenic analyses of cotton homologs of ArabidopsisPA regulator TRANSPARENT TESTA 2 (TT2) and fine-mapping of the cotton dark-brown fibre gene (Lc1), we show that a TT2 homolog, GhTT2-3A, controls PA biosynthesis and brown pigmentation in cotton fibres. We observed that GhTT2-3A activated GhbHLH130D, a homolog of ArabidopsisTT8, which in turn synergistically acted with GhTT2-3A to activate downstream PA structural genes and PA synthesis and accumulation in cotton fibres. Furthermore, the up-regulation of GhTT2-3A in fibres at the secondary wall-thickening stage resulted in brown mature fibres, and fibre quality and lint percentage were comparable to that of the white-fibre control. The findings of this study reveal the regulatory mechanism controlling brown pigmentation in cotton fibres and demonstrate a promising biotechnological strategy to break the negative linkage between coloration and fibre quality and/or productivity.

Long noncoding RNAs involve in resistance to Verticillium dahliae, a fungal disease in cotton

Long noncoding RNAs (lncRNAs) have several known functions in plant development, but their possible roles in responding to plant disease remain largely unresolved. In this study, we described a comprehensive disease-responding lncRNA profiles in defence against a cotton fungal disease Verticillium dahliae. We further revealed the conserved and specific characters of disease-responding process between two cotton species. Conservatively for two cotton species, we found the expression dominance of induced lncRNAs in the Dt subgenome, indicating a biased induction pattern in the co-existing subgenomes of allotetraploid cotton. Comparative analysis of lncRNA expression and their proposed functions in resistant Gossypium barbadense cv. '7124' versus susceptible Gossypium hirsutum cv. 'YZ1' revealed their distinct disease response mechanisms. Species-specific (LS) lncRNAs containing more SNPs displayed a fiercer inducing level postinfection than the species-conserved (core) lncRNAs. Gene Ontology enrichment of LS lncRNAs and core lncRNAs indicates distinct roles in the process of biotic stimulus. Further functional analysis showed that two core lncRNAs, GhlncNAT-ANX2- and GhlncNAT-RLP7-silenced seedlings, displayed an enhanced resistance towards V. dahliae and Botrytis cinerea, possibly associated with the increased expression of LOX1 and LOX2. This study represents the first characterization of lncRNAs involved in resistance to fungal disease and provides new clues to elucidate cotton disease response mechanism.

Flavonoid biosynthesis controls fiber color in naturally colored cotton

The existence of only natural brown and green cotton fibers (BCF and GCF, respectively), as well as poor fiber quality, limits the use of naturally colored cotton (Gossypium hirsutum L.). A better understanding of fiber pigment regulation is needed to surmount these obstacles. In this work, transcriptome analysis and quantitative reverse transcription PCR revealed that 13 and 9 phenylpropanoid (metabolic) pathway genes were enriched during pigment synthesis, while the differential expression of phenylpropanoid (metabolic) and flavonoid metabolic pathway genes occurred among BCF, GCF, and white cotton fibers (WCF). Silencing the chalcone flavanone isomerase gene in a BCF line resulted in three fiber phenotypes among offspring of the RNAi lines: BCF, almost WCF, and GCF. The lines with almost WCF suppressed chalcone flavanone isomerase, while the lines with GCF highly expressed the glucosyl transferase (3GT) gene. Overexpression of the Gh3GT or Arabidopsis thaliana 3GT gene in BCF lines resulted in GCF. Additionally, the phenylpropanoid and flavonoid metabolites of BCF and GCF were significantly higher than those of WCF as assessed by a metabolomics analysis. Thus, the flavonoid biosynthetic pathway controls both brown and green pigmentation processes. Like natural colored fibers, the transgenic colored fibers were weaker and shorter than WCF. This study shows the potential of flavonoid pathway modifications to alter cotton fibers’ color and quality.

Inhibition of Heat Shock proteins HSP90 and HSP70 induce oxidative stress, suppressing cotton fiber development

Cotton fiber is a specialized unicellular structure useful for the study of cellular differentiation and development. Heat shock proteins (HSPs) have been shown to be involved in various developmental processes. Microarray data analysis of five Gossypium hirsutum genotypes revealed high transcript levels of GhHSP90 and GhHSP70 genes at different stages of fiber development, indicating their importance in the process. Further, we identified 26 and 55 members of HSP90 and HSP70 gene families in G. hirsutum. The treatment of specific inhibitors novobiocin (Nov; HSP90) and pifithrin/2-phenylethynesulfonamide (Pif; HSP70) in in-vitro cultured ovules resulted in a fewer number of fiber initials and retardation in fiber elongation. The molecular chaperone assay using bacterially expressed recombinant GhHSP90-7 and GhHSP70-8 proteins further confirmed the specificity of inhibitors. HSP inhibition disturbs the H₂O₂ balance that leads to the generation of oxidative stress, which consequently results in autophagy in the epidermal layer of the cotton ovule. Transmission electron microscopy (TEM) of inhibitor-treated ovule also corroborates autophagosome formation along with disrupted mitochondrial cristae. The perturbations in transcript profile of HSP inhibited ovules show differential regulation of different stress and fiber development-related genes and pathways. Altogether, our results indicate that HSP90 and HSP70 families play a crucial role in cotton fiber differentiation and development by maintaining cellular homeostasis.

Linkage and association mapping reveals the genetic basis of brown fibre (Gossypium hirsutum)

Brown fibre cotton is an environmental-friendly resource that plays a key role in the textile industry. However, the fibre quality and yield of natural brown cotton are poor, and fundamental research on brown cotton is relatively scarce. To understand the genetic basis of brown fibre cotton, we constructed linkage and association populations to systematically examine brown fibre accessions. We fine-mapped the brown fibre region, Lc1 , and dissected it into 2 loci, qBF-A07-1 and qBF-A07-2. The qBF-A07-1 locus mediates the initiation of brown fibre production, whereas the shade of the brown fibre is affected by the interaction between qBF-A07-1 and qBF-A07-2. Gh_A07G2341 and Gh_A07G0100 were identified as candidate genes for qBF-A07-1 and qBF-A07-2, respectively. Haploid analysis of the signals significantly associated with these two loci showed that most tetraploid modern brown cotton accessions exhibit the introgression signature of Gossypium barbadense. We identified 10 quantitative trait loci (QTLs) for fibre yield and 19 QTLs for fibre quality through a genome-wide association study (GWAS) and found that qBF-A07-2 negatively affects fibre yield and quality through an epistatic interaction with qBF-A07-1. This study sheds light on the genetics of fibre colour and lint-related traits in brown fibre cotton, which will guide the elite cultivars breeding of brown fibre cotton.

The GhmiR157a-GhSPL10 regulatory module controls initial cellular dedifferentiation and callus proliferation in cotton by modulating ethylene-mediated flavonoid biosynthesis

MicroRNAs (miRNAs) modulate many biological processes through inactivation of specific mRNA targets such as those encoding transcription factors. A delicate spatial/temporal balance between specific miRNAs and their targets is central to achieving the appropriate biological outcomes. Somatic embryogenesis in cotton (Gossypium hirsutum), which goes through initial cellular dedifferentiation, callus proliferation, and somatic embryo development, is of great importance for both fundamental research and biotechnological applications. In this study, we characterize the function of the GhmiR157a-GhSPL10 miRNA-transcription factor module during somatic embryogenesis in cotton. We show that overexpression of GhSPL10, a target of GhmiR157a, increases free auxin and ethylene content and expression of associated signaling pathways, activates the flavonoid biosynthesis pathway, and promotes initial cellular dedifferentiation and callus proliferation. Inhibition of expression of the flavonoid synthesis gene F3H in GhSPL10 overexpression lines (35S:rSPL10-7) blocked callus initiation, while exogenous application of several types of flavonol promoted callus proliferation, associated with cell cycle-related gene expression. Inhibition of ethylene synthesis by aminoethoxyvinylglycine treatment in the 35S:rSPL10-7 line severely inhibited callus initiation, while activation of ethylene signaling through 1-aminocyclopropane 1-carboxylic acid treatment, EIN2 overexpression, or inhibition of the ethylene negative regulator CTR1 by RNA interference promoted flavonoid-related gene expression and flavonol accumulation. These results show that an up-regulation of ethylene signaling and the activation of flavonoid biosynthesis in GhSPL10 overexpression lines were associated with initial cellular dedifferentiation and callus proliferation. Our results demonstrate the importance of a GhmiR157a-GhSPL10 gene module in regulating somatic embryogenesis via hormonal and flavonoid pathways.