Recent insights into cotton functional genomics: progress and future perspectives

Unveiling the power of MYB transcription factors: Master regulators of multi-stress responses and development in cotton

Plants, being immobile, are subject to environmental stresses more than other creatures, necessitating highly effective stress tolerance systems. Transcription factors (TFs) play a crucial role in the adaptation mechanism as they can be activated by diverse signals and ultimately control the expression of stress-responsive genes. One of the most prominent plant TFs family is MYB (myeloblastosis), which is involved in secondary metabolites, developmental mechanisms, biological processes, cellular architecture, metabolic pathways, and stress responses. Extensive research has been conducted on the involvement of MYB TFs in crops, while their role in cotton remains largely unexplored. We also utilized genome-wide data to discover potential 440 MYB genes and investigated their plausible roles in abiotic and biotic stress conditions, as well as in different tissues across diverse transcriptome databases. This review primarily summarized the structure and classification of MYB TFs biotic and abiotic stress tolerance and their role in secondary metabolism in different crops, especially in cotton. However, it intends to identify gaps in current knowledge and emphasize the need for further research to enhance our understanding of MYB roles in plants.

Yield-related quantitative trait loci identification and lint percentage hereditary dissection under salt stress in upland cotton

Salinity is frequently mentioned as a major constraint in worldwide agricultural production. Lint percentage (LP) is a crucial yield-component in cotton lint production. While the genetic factors affect cotton yield in saline soils are still unclear. Here, we employed a recombinant inbred line population in upland cotton (Gossypium hirsutum L.) and investigated the effects of salt stress on five yield and yield component traits, including seed cotton yield per plant, lint yield per plant, boll number per plant, boll weight, and LP. Between three datasets of salt stress (E1), normal growth (E2), and the difference values dataset of salt stress and normal conditions (D-value), 87, 82, and 55 quantitative trait loci (QTL) were detectable, respectively. In total, five QTL (qLY-Chr6-2, qBNP-Chr4-1, qBNP-Chr12-1, qBNP-Chr15-5, qLP-Chr19-2) detected in both in E1 and D-value were salt related QTL, and three stable QTL (qLP-Chr5-3, qLP-Chr13-1, qBW-Chr5-5) were detected both in E1 and E2 across 3 years. Silencing of nine genes within a stable QTL (qLP-Chr5-3) highly expressed in fiber developmental stages increased LP and decreased fiber length (FL), indicating that multiple minor-effect genes clustered on Chromosome 5 regulate LP and FL. Additionally, the difference in LP caused by Gh_A05G3226 is mainly in transcription level rather than in the sequence difference. Moreover, silencing of salt related gene (GhDAAT) within qBNP-Chr4-1 decreased salt tolerance in cotton. Our findings shed light on the regulatory mechanisms underlining cotton salt tolerance and fiber initiation.

Microbial influencers and cotton leaf curl disease (CLCuD) susceptibility: a network perspective

Biotic stresses, such as plant viruses, e.g., cotton leaf curl virus (CLCuV), can alter root-associated and leaf-associated microbial diversities in plants. There are complex ecological dynamics at play, with each microbe contributing to a multitude of biotic and abiotic interactions, thus deciding the stability of the plant's ecosystem in response to the disease. Deciphering these networks of interactions is a challenging task. The inferential research in microbiome is also at a nascent stage, often constrained by the underlying analytical assumptions and the limitations with respect to the depth of sequencing. There is also no real consensus on network-wide statistics to identify the influential microbial players in a network. Guided by the latest developments in network science, including recently published metrics such as Integrated View of Influence (IVI) and some other centrality measures, this study provides an exposé of the most influential nodes in the rhizospheric and phyllospheric microbial networks of the cotton leaf curl disease (CLCuD) susceptible, partially tolerant, and resistant cotton varieties. It is evident from our results that the CLCuD-resistant Gossypium arboreum possesses an equal share of keystone species, which helps it to withstand ecological pressures. In the resistant variety, the phyllosphere harbors the most influential nodes, whereas in the susceptible variety, they are present in the rhizosphere. Based on hubness score, spreading score, and IVI, the top 10 occurring keystone species in the FDH-228 (resistant) variety include Actinokineospora, Cohnella, Thermobacillus, Clostridium, Desulfofarcimen, and MDD-D21. Elusimicrobia, Clostridium-sensu-stricto_12, Candidatus woesebacteria, and Dyella were identified as the most influential nodes in the PFV-1 (partially tolerant) variety. In the PFV-2 (susceptible) variety, the keystone species were identified as Georginia, Nesterenkonia, Elusimicrobia MVP-88, Acetivibrio, Tepedisphaerales, Chelatococcus, Nitrosospira, and RCP2-54. This concept deciphers the diseased and healthy plant's response to viral disease, which may be microbially mediated.

Cotton SNARE complex component GhSYP121 regulates salicylic acid signaling during defense against Verticillium dahliae

Syntaxin of plant (SYP) plays a crucial role in SNARE-mediated membrane trafficking during endocytic and secretory pathways, contributing to the regulation and execution of plant immunity against pathogens. Verticillium wilt is among the most destructive fungal diseases affecting cotton worldwide. However, information regarding SYP family genes in cotton is scarce. Through genome-wide identification and transcriptome profiling, we identified GhSYP121, a Qa SNARE gene in Gossypium hirsutum. GhSYP121 is notably induced by Verticillium dahliae, the causal agent of Verticillium wilt in cotton, and acts as a negative regulator of defense against V. dahliae. This is evidenced by the reduced resistance of GhSYP121-deficient cotton and the increased susceptibility of GhSYP121-overexpressing lines. Furthermore, the activation of the salicylic acid (SA) pathway by V. dahliae is inversely correlated with the expression level of GhSYP121. GhSYP121 interacts with its cognate SNARE component, GhSNAP33, which is required for the penetration resistance against V. dahliae in cotton. Collectively, GhSYP121, as a member of the cotton SNARE complex, is involved in regulating the SA pathway during plant defense against V. dahliae. This finding enhances our understanding of the potential role of GhSYP121 in these distinct pathways that contribute to plant defense against V. dahliae infection.

A comprehensive review on Gossypium hirsutum resistance against cotton leaf curl virus

Cotton (Gossypium hirsutum L.) is a significant fiber crop. Being a major contributor to the textile industry requires continuous care and attention. Cotton is subjected to various biotic and abiotic constraints. Among these, biotic factors including cotton leaf curl virus (CLCuV) are dominant. CLCuV is a notorious disease of cotton and is acquired, carried, and transmitted by the whitefly (Bemisia tabaci). A cotton plant affected with CLCuV may show a wide range of symptoms such as yellowing of leaves, thickening of veins, upward or downward curling, formation of enations, and stunted growth. Though there are many efforts to protect the crop from CLCuV, long-term results are not yet obtained as CLCuV strains are capable of mutating and overcoming plant resistance. However, systemic-induced resistance using a gene-based approach remained effective until new virulent strains of CLCuV (like Cotton Leaf Curl Burewala Virus and others) came into existence. Disease control by biological means and the development of CLCuV-resistant cotton varieties are in progress. In this review, we first discussed in detail the evolution of cotton and CLCuV strains, the transmission mechanism of CLCuV, the genetic architecture of CLCuV vectors, and the use of pathogen and nonpathogen-based approaches to control CLCuD. Next, we delineate the uses of cutting-edge technologies like genome editing (with a special focus on CRISPR-Cas), next-generation technologies, and their application in cotton genomics and speed breeding to develop CLCuD resistant cotton germplasm in a short time. Finally, we delve into the current obstacles related to cotton genome editing and explore forthcoming pathways for enhancing precision in genome editing through the utilization of advanced genome editing technologies. These endeavors aim to enhance cotton's resilience against CLCuD.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Identification, evolution, and expression of GDSL-type Esterase/Lipase (GELP) gene family in three cotton species: a bioinformatic analysis

BACKGROUND

GDSL esterase/lipases (GELPs) play important roles in plant growth, development, and response to biotic and abiotic stresses. Presently, an extensive and in-depth analysis of GELP family genes in cotton is still not clear enough, which greatly limits the further understanding of cotton GELP function and regulatory mechanism.

RESULTS

A total of 389 GELP family genes were identified in three cotton species of Gossypium hirsutum (193), G. arboreum (97), and G. raimondii (99). These GELPs could be classified into three groups and eight subgroups, with the GELPs in same group to have similar gene structures and conserved motifs. Evolutionary event analysis showed that the GELP family genes tend to be diversified at the spatial dimension and certain conservative at the time dimension, with a trend of potential continuous expansion in the future. The orthologous or paralogous GELPs among different genomes/subgenomes indicated the inheritance from genome-wide duplication during polyploidization, and the paralogous GELPs were derived from chromosomal segment duplication or tandem replication. GELP genes in the A/D subgenome underwent at least three large-scale replication events in the evolutionary process during the period of 0.6-3.2 MYA, with two large-scale evolutionary events between 0.6-1.8 MYA that were associated with tetraploidization, and the large-scale duplication between 2.6-9.1 MYA that occurred during diploidization. The cotton GELPs indicated diverse expression patterns in tissue development, ovule and fiber growth, and in response to biotic and abiotic stresses, combining the existing cis-elements in the promoter regions, suggesting the GELPs involvements of functions to be diversification and of the mechanisms to be a hormone-mediated manner.

CONCLUSIONS

Our results provide a systematic and comprehensive understanding the function and regulatory mechanism of cotton GELP family, and offer an effective reference for in-depth genetic improvement utilization of cotton GELPs.

Origin, evolution, and diversification of the wall-associated kinase gene family in plants

KEY MESSAGE

The study of the origin, evolution, and diversification of the wall-associated kinase gene family in plants facilitates their functional investigations in the future. Wall-associated kinases (WAKs) make up one subfamily of receptor-like kinases (RLKs), and function directly in plant cell elongation and responses to biotic and abiotic stresses. The biological functions of WAKs have been extensively characterized in angiosperms; however, the origin and evolutionary history of the WAK family in green plants remain unclear. Here, we performed a comprehensive analysis of the WAK family to reveal its origin, evolution, and diversification in green plants. In total, 1061 WAK genes were identified in 37 species from unicellular algae to multicellular plants, and the results showed that WAK genes probably originated before bryophyte differentiation and were widely distributed in land plants, especially angiosperms. The phylogeny indicated that the land plant WAKs gave rise to five clades and underwent lineage-specific expansion after species differentiation. Cis-acting elements and expression patterns analyses of WAK genes in Arabidopsis and rice demonstrated the functional diversity of WAK genes in these two species. Many gene gains and losses have occurred in angiosperms, leading to an increase in the number of gene copies. The evolutionary trajectory of the WAK family during polyploidization was uncovered using Gossypium species. Our results provide insights into the evolution of WAK genes in green plants, facilitating their functional investigations in the future.

Transgenerational effect of Afidopyropen on Bemisia tabaci Gennadius (Homoptera: Aleyrodidae)

Bemisia tabaci Gennadius (Homoptera: Aleyrodidae) is the most devastating insect-pest in cotton crop. It is vector of the cotton leaf curl virus (CLCV) and is responsible for huge losses to cotton industry in Pakistan and worldwide. It is mainly controlled by insecticides but the injudicious use of insecticides has resulted in insecticide resistance and population resurgence in addition to various harmful effects on the humans, non-target organisms and the environment. Transgenerational studies are very helpful to choose a best insecticidal option. In the current study, age-stage two-sex life table analysis was used to identify transgenerational effects of sublethal doses of afidopyropen. The adults of B. tabaci were treated with three concentrations of afidopyropen i.e., LC10, LC30 and LC50. The results indicated significant changes in the progeny i.e. the fecundity decreased in treated population; and female and male longevity of their progeny were more in control as compared to treated populations. Similarly, population parameters like intrinsic rate of growth (r), net reproductive rate (R0) and limiting rate of growth (λ) were significantly decreased in the treated adult progeny with values of 0.08-0.11, 4.85-7.46 and 1.09-1.12 per day, respectively. Based on the reduced biotic potential, afidopyropen can be suggested as an effective alternative option for the management of B. tabaci.

Complete chloroplast genome assembly and phylogenetic analysis of blackcurrant (Ribes nigrum), red and white currant (Ribes rubrum), and gooseberry (Ribes uva-crispa) provide new insights into the phylogeny of Grossulariaceae

Background

Blackcurrant (Ribes nigrum), red currant (R. rubrum), white currant (R. rubrum), and gooseberry (R. uva-crispa) belong to Grossulariaceae and are popular small-berry crops worldwide. The lack of genomic data has severely limited their systematic classification and molecular breeding.

Methods

The complete chloroplast (cp) genomes of these four taxa were assembled for the first time using MGI-DNBSEQ reads, and their genome structures, repeat elements and protein-coding genes were annotated. By genomic comparison of the present four and previous released five Ribes cp genomes, the genomic variations were identified. By phylogenetic analysis based on maximum-likelihood and Bayesian methods, the phylogeny of Grossulariaceae and the infrageneric relationships of the Ribes were revealed.

Results

The four cp genomes have lengths ranging from 157,450 to 157,802 bp and 131 shared genes. A total of 3,322 SNPs and 485 Indels were identified from the nine released Ribes cp genomes. Red currant and white currant have 100% identical cp genomes partially supporting the hypothesis that white currant (R. rubrum) is a fruit color variant of red currant (R. rubrum). The most polymorphic genic and intergenic region is ycf1 and trnT-psbD, respectively. The phylogenetic analysis demonstrated the monophyly of Grossulariaceae in Saxifragales and the paraphyletic relationship between Saxifragaceae and Grossulariaceae. Notably, the Grossularia subgenus is well nested within the Ribes subgenus and shows a paraphyletic relationship with the co-ancestor of Calobotrya and Coreosma sections, which challenges the dichotomous subclassification of the Ribes genus based on morphology (subgenus Ribes and subgenus Grossularia). These data, results, and insights lay a foundation for the phylogenetic research and breeding of Ribes species.

Comprehensive Evaluation and Transcriptome Analysis Reveal the Salt Tolerance Mechanism in Semi-Wild Cotton (Gossypium purpurascens)

Elevated salinity significantly threatens cotton growth, particularly during the germination and seedling stages. The utilization of primitive species of Gossypium hirsutum, specifically Gossypium purpurascens, has the potential to facilitate the restoration of genetic diversity that has been depleted due to selective breeding in modern cultivars. This investigation evaluated 45 G. purpurascens varieties and a salt-tolerant cotton variety based on 34 morphological, physiological, and biochemical indicators and comprehensive salt tolerance index values. This study effectively identified a total of 19 salt-tolerant and two salt-resistant varieties. Furthermore, transcriptome sequencing of a salt-tolerant genotype (Nayanmian-2; NY2) and a salt-sensitive genotype (Sanshagaopao-2; GP2) revealed 2776, 6680, 4660, and 4174 differentially expressed genes (DEGs) under 0.5, 3, 12, and 24 h of salt stress. Gene ontology enrichment analysis indicated that the DEGs exhibited significant enrichment in biological processes like metabolic (GO:0008152) and cellular (GO:0009987) processes. MAPK signaling, plant-pathogen interaction, starch and sucrose metabolism, plant hormone signaling, photosynthesis, and fatty acid metabolism were identified as key KEGG pathways involved in salinity stress. Among the DEGs, including NAC, MYB, WRKY, ERF, bHLH, and bZIP, transcription factors, receptor-like kinases, and carbohydrate-active enzymes were crucial in salinity tolerance. Weighted gene co-expression network analysis (WGCNA) unveiled associations of salt-tolerant genotypes with flavonoid metabolism, carbon metabolism, and MAPK signaling pathways. Identifying nine hub genes (MYB4, MYB105, MYB36, bZIP19, bZIP43, FRS2 SMARCAL1, BBX21, F-box) across various intervals offered insights into the transcriptional regulation mechanism of salt tolerance in G. purpurascens. This study lays the groundwork for understanding the important pathways and gene networks in response to salt stress, thereby providing a foundation for enhancing salt tolerance in upland cotton.

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.

The Na+/H+ antiporter GbSOS1 interacts with SIP5 and regulates salt tolerance in Gossypium barbadense

Cotton is a globally cultivated economic crop and is a major source of natural fiber and edible oil. However, cotton production is severely affected by salt stress. Although Salt Overly Sensitive 1 (SOS1) is a well-studied Na⁺/H⁺ antiporter in multiple plant species, little is known about its function and regulatory mechanism in cotton. Here, we cloned a salt-induced SOS1 from sea-island cotton. Real-time quantitative PCR analysis revealed that GbSOS1 was induced by multiple stresses and phytohormones. Silencing GbSOS1 through virus-induced gene silencing significantly reduced cotton resistance to high Na⁺ but mildly affected Li⁺ tolerance. On the other hand, overexpression of GbSOS1 enhanced salt tolerance in yeast, Arabidopsis, and cotton largely due to the ability to maintain Na⁺ homeostasis in protoplasts. Yeast-two-hybrid assays and bimolecular fluorescence complementation identified a novel protein interacting with GbSOS1 on the plasma membrane, which we named SOS Interaction Protein 5 (SIP5). We found that the SIP5 gene encoded an unknown protein localized on the cell membrane. Silencing SIP5 significantly increased cotton tolerance to salt, exhibited by less wilting and plant death under salt stress. Our results revealed that GbSOS1 is crucial for cotton survival in saline soil, and SIP5 is a potentially negative regulator of SOS1-mediated salt tolerance in cotton. Overall, this study provides a theoretical basis for elucidating the molecular mechanism of SOS1, and a candidate gene for breeding salt-tolerant crops.

Systematic analysis and expression of Gossypium 2ODD superfamily highlight the roles of GhLDOXs responding to alkali and other abiotic stress in cotton

BACKGROUND

2-oxoglutarate-dependent dioxygenase (2ODD) is the second largest family of oxidases involved in various oxygenation/hydroxylation reactions in plants. Many members in the family regulate gene transcription, nucleic acid modification/repair and secondary metabolic synthesis. The 2ODD family genes also function in the formation of abundant flavonoids during anthocyanin synthesis, thereby modulating plant development and response to diverse stresses.

RESULTS

Totally, 379, 336, 205, and 204 2ODD genes were identified in G. barbadense (Gb), G. hirsutum (Gh), G. arboreum (Ga), and G. raimondii (Gb), respectively. The 336 2ODDs in G. hirsutum were divided into 15 subfamilies according to their putative functions. The structural features and functions of the 2ODD members in the same subfamily were similar and evolutionarily conserved. Tandem duplications and segmental duplications served essential roles in the large-scale expansion of the cotton 2ODD family. Ka/Ks values for most of the gene pairs were less than 1, indicating that 2ODD genes undergo strong purifying selection during evolution. Gh2ODDs might act in cotton responses to different abiotic stresses. GhLDOX3 and GhLDOX7, two members of the GhLDOX subfamily from Gh2ODDs, were significantly down-regulated in transcription under alkaline stress. Moreover, the expression of GhLDOX3 in leaves was significantly higher than that in other tissues. These results will provide valuable information for further understanding the evolution mechanisms and functions of the cotton 2ODD genes in the future.

CONCLUSIONS

Genome-wide identification, structure, and evolution and expression analysis of 2ODD genes in Gossypium were carried out. The 2ODDs were highly conserved during evolutionary. Most Gh2ODDs were involved in the regulation of cotton responses to multiple abiotic stresses including salt, drought, hot, cold and alkali.

Potassium ameliorates cotton (Gossypium hirsutum L.) fiber length by regulating osmotic and K+ /Na+ homeostasis under salt stress

Potassium (K) application can alleviate cotton salt stress, but the regulatory mechanisms affecting cotton fiber elongation and ion homeostasis are still unclear. A two-year field experiment was conducted to explore the effects of K on the osmolyte contents (soluble sugar, K⁺ content, and malate) and related enzyme activities during the fiber elongation of two cotton cultivars with contrasting salt sensitivity (CCRI-79; salt tolerant cultivar, and Simian 3; salt-sensitive cultivar). Three K application treatments (0, 150, and 300 kg K₂ O ha^-1 ) were applied at three soil salinity levels (low salinity, EC = 1.73 ± 0.05 dS m^-1 ; medium salinity, EC = 6.32 ± 0.10 dS m^-1 ; high salinity, EC = 10.84 ± 0.24 dS m^-1 ). K application improved fiber length and alleviated salt stress by increasing the maximum velocity of fiber elongation (V_max ). The increase rate of K on fiber length decreased with elevating salt stress, and the increase rate of K on V_max of CCRI-79 was greater than that of Simian3. K application can increase the enzyme activities (phosphoenolpyruvate carboxylase, PEPC, E.C. 4.1.1.31; pyrophosphatase, PPase, E.C. 3.6.1.1; and plasma membrane H⁺ -ATPase, PM H⁺ -ATPase, E.C. 3.6.3.14) as well as the content of osmolytes associated with the enzymes mentioned above. K increased the osmolyte contents under salt stress, and the increase in the K⁺ content of the fibers was much higher than that of soluble sugar and malate. The results of this study indicated K fertilizer application rates regulate the metabolism of osmolytes in cotton fiber development under salt stress, K⁺ is more critical to fiber elongation.

An overview of salinity stress, mechanism of salinity tolerance and strategies for its management in cotton

Salinity stress is one of the primary threats to agricultural crops resulting in impaired crop growth and development. Although cotton is considered as reasonably salt tolerant, it is sensitive to salt stress at some critical stages like germination, flowering, boll formation, resulting in reduced biomass and fiber production. The mechanism of partial ion exclusion (exclusion of Na+ and/or Cl-) in cotton appears to be responsible for the pattern of uptake and accumulation of harmful ions (Na+ and Cl) in tissues of plants exposed to saline conditions. Maintaining high tissue K+/Na+ and Ca2+/Na+ ratios has been proposed as a key selection factor for salt tolerance in cotton. The key adaptation mechanism in cotton under salt stress is excessive sodium exclusion or compartmentation. Among the cultivated species of cotton, Egyptian cotton (Gossypium barbadense L.) exhibit better salt tolerance with good fiber quality traits as compared to most cultivated cotton and it can be used to improve five quality traits and transfer salt tolerance into Upland or American cotton (Gossypium hirsutum L.) by interspecific introgression. Cotton genetic studies on salt tolerance revealed that the majority of growth, yield, and fiber traits are genetically determined, and controlled by quantitative trait loci (QTLs). Molecular markers linked to genes or QTLs affecting key traits have been identified, and they could be utilized as an indirect selection criterion to enhance breeding efficiency through marker-assisted selection (MAS). Transfer of genes for compatible solute, which are an important aspect of ion compartmentation, into salt-sensitive species is, theoretically, a simple strategy to improve tolerance. The expression of particular stress-related genes is involved in plant adaptation to environmental stressors. As a result, enhancing tolerance to salt stress can be achieved by marker assisted selection added with modern gene editing tools can boost the breeding strategies that defend and uphold the structure and function of cellular components. The intent of this review was to recapitulate the advancements in salt screening methods, tolerant germplasm sources and their inheritance, biochemical, morpho-physiological, and molecular characteristics, transgenic approaches, and QTLs for salt tolerance in cotton.

Role of WRKY Transcription Factors in Regulation of Abiotic Stress Responses in Cotton

Environmental factors are the major constraints in sustainable agriculture. WRKY proteins are a large family of transcription factors (TFs) that regulate various developmental processes and stress responses in plants, including cotton. On the basis of Gossypium raimondii genome sequencing, WRKY TFs have been identified in cotton and characterized for their functions in abiotic stress responses. WRKY members of cotton play a significant role in the regulation of abiotic stresses, i.e., drought, salt, and extreme temperatures. These TFs either activate or repress various signaling pathways such as abscisic acid, jasmonic acid, salicylic acid, mitogen-activated protein kinases (MAPK), and the scavenging of reactive oxygen species. WRKY-associated genes in cotton have been genetically engineered in Arabidopsis, Nicotiana, and Gossypium successfully, which subsequently enhanced tolerance in corresponding plants against abiotic stresses. Although a few review reports are available for WRKY TFs, there is no critical report available on the WRKY TFs of cotton. Hereby, the role of cotton WRKY TFs in environmental stress responses is studied to enhance the understanding of abiotic stress response and further improve in cotton plants.

Comparative Metabolome and Transcriptome Analysis of Anthocyanin Biosynthesis in White and Pink Petals of Cotton (Gossypium hirsutum L.)

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.

Comparative Ionomics and Metabolic Responses and Adaptive Strategies of Cotton to Salt and Alkali Stress

Soil salinization and alkalization severely inhibit agriculture. However, the response mechanisms of cotton to salt stress or alkali stress are unclear. Ionomics and metabolomics were used to investigate salt and alkali stresses in cotton roots and leaves. Compared with the control, salt-treated and alkali-treated cotton plants showed 51.8 and 53.0% decreases in biomass, respectively. Under salt stress, the concentration of N decreased in roots but increased in leaves, and the concentrations of P and K increased in roots but decreased in leaves. Salt stress inhibited Ca, B, N, and Fe uptake and Mg, K, P, S, and Cu transport, but promoted Mo, Mn, Zn, Mg, K, P, S, and Cu uptake and Mo, Mn, Zn, B, N, and Fe transport. Under alkali stress, the concentrations of N and P in roots and leaves decreased, while the concentrations of K in roots and leaves increased. Alkali stress inhibited P, Ca, S, N, Fe, and Zn uptake and N, P, Mg and B transport, but promoted K, Mn, Cu, Mo, Mg, and B uptake and K, Mn, Cu, Mo, Fe, and Zn transport. Under salt stress in the leaves, 93 metabolites increased, mainly organic acids, amino acids, and sugars, increased in abundance, while 6 decreased. In the roots, 72 metabolites increased, mainly amino acids, organic acids, and sugars, while 18 decreased. Under alkali stress, in the leaves, 96 metabolites increased, including organic acids, amino acids, and sugars, 83 metabolites decreased, including organic acids, amino acids, and sugars; In the roots, 108 metabolites increased, including organic acids, amino acids, and sugars. 83 metabolites decreased, including organic acids and amino acids. Under salt stress, cotton adapts to osmotic stress through the accumulation of organic acids, amino acids and sugars, while under alkali stress, osmoregulation was achieved via inorganic ion accumulation. Under salt stress, significant metabolic pathways in the leaves and roots were associated with amino acid and organic acid metabolism, sugar metabolism was mainly used as a source of energy, while under alkali stress, the pathways in the leaves were related to amino acid and linoleic acid metabolism, β-Oxidation, TCA cycle, and glycolysis were enhanced to provide the energy needed for life activities. Enhancing organic acid accumulation and metabolism in the roots is the key response mechanism of cotton to alkalinity.

The activation of gene expression and alternative splicing in the formation and evolution of allopolyploid Brassica napus

Allopolyploids contain two or more sets of subgenomes. To establish a compatible relationship between subgenomes, a series of gene expression changes occurred in allopolyploids. What evolutionary changes of transcripts have taken place in Brassica napus during the early establishment and subsequent evolution was a fascinating scientific question. Here, we study this issue using a set of materials (natural, resynthesized B. napus and their progenitors/parents) by long-read RNA sequencing technology. The results showed that more genes were up-regulated in resynthesized B. napus compared with its two parents, and more up-regulated expressed genes were observed in natural B. napus compared with resynthesized B. napus. The presence of up-regulation genes in organism may help it adapt to the influence of "genomic shock" and cope with natural environment. Isoforms are produced from precursor mRNAs by alternative splicing (AS) events, and more than 60% of novel isoforms were identified in all materials, which could improve the reference genome information of B. napus. We found that the isoform numbers, the number of genes potentially involved in AS and alternative polyadenylation increased in B. napus after evolution, which may involve in the adaptation of plants to natural environment. In addition, all identified isoforms were functional annotated by searching 7 databases. In general, this study could improve our overall understanding of the full-length transcriptome of B. napus, and help us recognize the significant gene expression changes and isoform abundance changes occurred in allopolyploid B. napus during evolution.

Determination of hormoligosis of organophosphate insecticides against Phenacoccus solenopsis

Abstract Cotton mealybug is a highly invasive pest of agricultural crops worldwide. Major agriculturists most rely on the use of insecticides for the control of pesticides. So, the indiscriminate use of insecticides leads to resistance development in recent years. For this purpose, an experiment was conducted using different concentrations of the three insecticides (profenfos chlorpyrifos and triazophos) to check the hormoligosis effects against cotton mealybug (CMB) in laboratory conditions. Investigation of variations for % mortality of adults of CMB after three days revealed that all treatments had statistically significant (P ˂ 0.05). The highest mortality was observed at the highest concentrations of profenofos 2.4% (38.55%). After 7 days, all the treatments were significant with difference in means (P ˂ 0.05). The highest mortality was recorded at the highest dilution of pesticide profenofos 2.4% (77.11%). The values of fecundity and longevity exposed a valid difference among treatments (P ˂ 0.05). Maximum fecundity was observed at the concentration 2.4% (181.41%) and longevity showed (38.46%). The highest mortality was observed at a concentration of triazophos 4% (27.98%). For chlorpyriphos the highest mortality was examined at concentration 4% (24.79%). The fecundity showed a statistically significant difference for different concentrations of triazophos and chlorpyriphos (P ˂ 0.05). The results of the recent study provide valuable information regarding the selection of insecticides and hormoligosis effects. The study can be helpful in the implications of integrated pest management of P. solenopsis.

Fourth generation biofuel from genetically modified algal biomass: Challenges and future directions

Genetic engineering applications in the field of biofuel are rapidly expanding due to their potential to boost biomass productivity while lowering its cost and enhancing its quality. Recently, fourth-generation biofuel (FGB), which is biofuel obtained from genetically modified (GM) algae biomass, has gained considerable attention from academic and industrial communities. However, replacing fossil resources with FGB is still beset with many challenges. Most notably, technical aspects of genetic modification operations need to be more fully articulated and elaborated. However, relatively little attention has been paid to GM algal biomass. There is a limited number of reviews on the progress and challenges faced in the algal genetics of FGB. Therefore, the present review aims to fill this gap in the literature by recapitulating the findings of recent studies and achievements on safe and efficient genetic manipulation in the production of FGB. Then, the essential issues and parameters related to genome editing in algal strains are highlighted. Finally, the main challenges to FGB pertaining to the diffusion risk and regulatory frameworks are addressed. This review concluded that the technical and biosafety aspects of FGB, as well as the complexity and diversity of the related regulations, legitimacy concerns, and health and environmental risks, are among the most important challenges that require a strong commitment at the national/international levels to reach a global consensus.

Using Multiplexed CRISPR/Cas9 for Suppression of Cotton Leaf Curl Virus

In recent decades, Pakistan has suffered a decline in cotton production due to several factors, including insect pests, cotton leaf curl disease (CLCuD), and multiple abiotic stresses. CLCuD is a highly damaging plant disease that seriously limits cotton production in Pakistan. Recently, genome editing through CRISPR/Cas9 has revolutionized plant biology, especially to develop immunity in plants against viral diseases. Here we demonstrate multiplex CRISPR/Cas-mediated genome editing against CLCuD using transient transformation in N. benthamiana plants and cotton seedlings. The genomic sequences of cotton leaf curl viruses (CLCuVs) were obtained from NCBI and the guide RNA (gRNA) were designed to target three regions in the viral genome using CRISPR MultiTargeter. The gRNAs were cloned in pHSE401/pKSE401 containing Cas9 and confirmed through colony PCR, restriction analysis, and sequencing. Confirmed constructs were moved into Agrobacterium and subsequently used for transformation. Agroinfilteration in N. benthamiana revealed delayed symptoms (3–5 days) with improved resistance against CLCuD. In addition, viral titer was also low (20–40%) in infected plants co-infiltrated with Cas9-gRNA, compared to control plants (infected with virus only). Similar results were obtained in cotton seedlings. The results of transient expression in N. benthamiana and cotton seedlings demonstrate the potential of multiplex CRISPR/Cas to develop resistance against CLCuD. Five transgenic plants developed from three experiments showed resistance (60−70%) to CLCuV, out of which two were selected best during evaluation and screening. The technology will help breeding CLCuD-resistant cotton varieties for sustainable cotton production.

Using Multiplexed CRISPR/Cas9 for Suppression of Cotton Leaf Curl Virus

In recent decades, Pakistan has suffered a decline in cotton production due to several factors, including insect pests, cotton leaf curl disease (CLCuD), and multiple abiotic stresses. CLCuD is a highly damaging plant disease that seriously limits cotton production in Pakistan. Recently, genome editing through CRISPR/Cas9 has revolutionized plant biology, especially to develop immunity in plants against viral diseases. Here we demonstrate multiplex CRISPR/Cas-mediated genome editing against CLCuD using transient transformation in N. benthamiana plants and cotton seedlings. The genomic sequences of cotton leaf curl viruses (CLCuVs) were obtained from NCBI and the guide RNA (gRNA) were designed to target three regions in the viral genome using CRISPR MultiTargeter. The gRNAs were cloned in pHSE401/pKSE401 containing Cas9 and confirmed through colony PCR, restriction analysis, and sequencing. Confirmed constructs were moved into Agrobacterium and subsequently used for transformation. Agroinfilteration in N. benthamiana revealed delayed symptoms (3–5 days) with improved resistance against CLCuD. In addition, viral titer was also low (20–40%) in infected plants co-infiltrated with Cas9-gRNA, compared to control plants (infected with virus only). Similar results were obtained in cotton seedlings. The results of transient expression in N. benthamiana and cotton seedlings demonstrate the potential of multiplex CRISPR/Cas to develop resistance against CLCuD. Five transgenic plants developed from three experiments showed resistance (60−70%) to CLCuV, out of which two were selected best during evaluation and screening. The technology will help breeding CLCuD-resistant cotton varieties for sustainable cotton production.

Multi-Omics-Based Identification and Functional Characterization of Gh_A06G1257 Proves Its Potential Role in Drought Stress Tolerance in Gossypium hirsutum

Cotton is one of the most important fiber crops globally. Despite this, various abiotic stresses, including drought, cause yield losses. We used transcriptome profiles to investigate the co-expression patterns of gene networks associated with drought stress tolerance. We identified three gene modules containing 3,567 genes highly associated with drought stress tolerance. Within these modules, we identified 13 hub genes based on intramodular significance, for further validation. The yellow module has five hub genes (Gh_A07G0563, Gh_D05G0221, Gh_A05G3716, Gh_D12G1438, and Gh_D05G0697), the brown module contains three hub genes belonging to the aldehyde dehydrogenase (ALDH) gene family (Gh_A06G1257, Gh_A06G1256, and Gh_D06G1578), and the pink module has five hub genes (Gh_A02G1616, Gh_D12G2599, Gh_D07G2232, Gh_A02G0527, and Gh_D07G0629). Based on RT-qPCR results, the Gh_A06G1257 gene has the highest expression under drought stress in different plant tissues and it might be the true candidate gene linked to drought stress tolerance in cotton. Silencing of Gh_A06G1257 in cotton leaves conferred significant sensitivity in response to drought stress treatments. Overexpression of Gh_A06G1257 in Arabidopsis also confirms its role in drought stress tolerance. L-valine, Glutaric acid, L-proline, L-Glutamic acid, and L-Tryptophan were found to be the most significant metabolites playing roles in drought stress tolerance. These findings add significantly to existing knowledge of drought stress tolerance mechanisms in cotton.

Gene Expression Correlation Analysis Reveals MYC-NAC Regulatory Network in Cotton Pigment Gland Development

Plant NAC (NAM, ATAF1/2, and CUC2) family is involved in various development processes including Programmed Cell Death (PCD) associated development. However, the relationship between NAC family and PCD-associated cotton pigment gland development is largely unknown. In this study, we identified 150, 153 and 299 NAC genes in newly updated genome sequences of G. arboreum, G. raimondii and G. hirsutum, respectively. All NAC genes were divided into 8 groups by the phylogenetic analysis and most of them were conserved during cotton evolution. Using the vital regulator of gland formation GhMYC2-like as bait, expression correlation analysis screened out 6 NAC genes which were low-expressed in glandless cotton and high-expressed in glanded cotton. These 6 NAC genes acted downstream of GhMYC2-like and were induced by MeJA. Silencing CGF1(Cotton Gland Formation1), another MYC-coding gene, caused almost glandless phenotype and down-regulated expression of GhMYC2-like and the 6 NAC genes, indicating a MYC-NAC regulatory network in gland development. In addition, predicted regulatory mechanism showed that the 6 NAC genes were possibly regulated by light, various phytohormones and transcription factors as well as miRNAs. The interaction network and DNA binding sites of the 6 NAC transcription factors were also predicted. These results laid the foundation for further study of gland-related genes and gland development regulatory network.

Efficient CRISPR/Cas9 mediated Pooled-sgRNAs assembly accelerates targeting multiple genes related to male sterility in cotton

BACKGROUND

Upland cotton (Gossypium hirsutum), harboring a complex allotetraploid genome, consists of A and D sub-genomes. Every gene has multiple copies with high sequence similarity that makes genetic, genomic and functional analyses extremely challenging. The recent accessibility of CRISPR/Cas9 tool provides the ability to modify targeted locus efficiently in various complicated plant genomes. However, current cotton transformation method targeting one gene requires a complicated, long and laborious regeneration process. Hence, optimizing strategy that targeting multiple genes is of great value in cotton functional genomics and genetic engineering.

RESULTS

To target multiple genes in a single experiment, 112 plant development-related genes were knocked out via optimized CRISPR/Cas9 system. We optimized the key steps of pooled sgRNAs assembly method by which 116 sgRNAs pooled together into 4 groups (each group consisted of 29 sgRNAs). Each group of sgRNAs was compiled in one PCR reaction which subsequently went through one round of vector construction, transformation, sgRNAs identification and also one round of genetic transformation. Through the genetic transformation mediated Agrobacterium, we successfully generated more than 800 plants. For mutants identification, Next Generation Sequencing technology has been used and results showed that all generated plants were positive and all targeted genes were covered. Interestingly, among all the transgenic plants, 85% harbored a single sgRNA insertion, 9% two insertions, 3% three different sgRNAs insertions, 2.5% mutated sgRNAs. These plants with different targeted sgRNAs exhibited numerous combinations of phenotypes in plant flowering tissues.

CONCLUSION

All targeted genes were successfully edited with high specificity. Our pooled sgRNAs assembly offers a simple, fast and efficient method/strategy to target multiple genes in one time and surely accelerated the study of genes function in cotton.

The Cotton Wall-Associated Kinase GhWAK7A Mediates Responses to Fungal Wilt Pathogens by Complexing with the Chitin Sensory Receptors

Plant receptor-like kinases (RLKs) are important players in response to pathogen infections. Verticillium and Fusarium wilts, caused by Verticillium dahliae (Vd) and Fusarium oxysporum f. sp vasinfectum (Fov), respectively, are among the most devastating diseases in cotton (Gossypium spp). To understand the cotton response to these soil-borne fungal pathogens, we performed a genome-wide in silico characterization and functional screen of diverse RLKs for their involvement in cotton wilt diseases. We identified Gossypium hirsutum GhWAK7A, a wall-associated kinase, that positively regulates cotton response to both Vd and Fov infections. Chitin, the major constituent of the fungal cell wall, is perceived by lysin-motif-containing RLKs (LYKs/CERK1), leading to the activation of plant defense against fungal pathogens. A conserved chitin sensing and signaling system is present in cotton, including chitin-induced GhLYK5-GhCERK1 dimerization and phosphorylation, and contributes to cotton defense against Vd and Fov Importantly, GhWAK7A directly interacts with both GhLYK5 and GhCERK1 and promotes chitin-induced GhLYK5-GhCERK1 dimerization. GhWAK7A phosphorylates GhLYK5, which itself does not have kinase activity, but requires phosphorylation for its function. Consequently, GhWAK7A plays a crucial role in chitin-induced responses. Thus, our data reveal GhWAK7A as an important component in cotton response to fungal wilt pathogens by complexing with the chitin receptors.

Characterization and Stress Response of the JmjC Domain-Containing Histone Demethylase Gene Family in the Allotetraploid Cotton Species Gossypium hirsutum

Histone modification is an important epigenetic modification that controls gene transcriptional regulation in eukaryotes. Histone methylation is accomplished by histone methyltransferase and can occur on two amino acid residues, arginine and lysine. JumonjiC (JmjC) domain-containing histone demethylase regulates gene transcription and chromatin structure by changing the methylation state of the lysine residue site and plays an important role in plant growth and development. In this study, we carried out genome-wide identification and comprehensive analysis of JmjC genes in the allotetraploid cotton species Gossypium hirsutum. In total, 50 JmjC genes were identified and in G. hirsutum, and 25 JmjC genes were identified in its two diploid progenitors, G. arboreum and G. raimondii, respectively. Phylogenetic analysis divided these JmjC genes into five subfamilies. A collinearity analysis of the two subgenomes of G. hirsutum and the genomes of G. arboreum and G. raimondii uncovered a one-to-one relationship between homologous genes of the JmjC gene family. Most homologs in the JmjC gene family between A and D subgenomes of G. hirsutum have similar exon-intron structures, which indicated that JmjC family genes were conserved after the polyploidization. All G. hirsutumJmjC genes were found to have a typical JmjC domain, and some genes also possess other special domains important for their function. Analysis of promoter regions revealed that cis-acting elements, such as those related to hormone and abiotic stress response, were enriched in G. hirsutum JmjC genes. According to a reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis, most G. hirsutumJmjC genes had high abundance expression at developmental stages of fibers, suggesting that they might participate in cotton fiber development. In addition, some G. hirsutumJmjC genes were found to have different degrees of response to cold or osmotic stress, thus indicating their potential role in these types of abiotic stress response. Our results provide useful information for understanding the evolutionary history and biological function of JmjC genes in cotton.

Application of compound material alleviates saline and alkaline stress in cotton leaves through regulation of the transcriptome

BACKGROUND

Soil salinization and alkalinization are the main factors that affect the agricultural productivity. Evaluating the persistence of the compound material applied in field soils is an important part of the regulation of the responses of cotton to saline and alkaline stresses.

RESULT

To determine the molecular effects of compound material on the cotton's responses to saline stress and alkaline stress, cotton was planted in the salinized soil (NaCl 8 g kg^- 1) and alkalized soil (Na₂CO₃ 8 g kg^- 1) after application of the compound material, and ion content, physiological characteristics, and transcription of new cotton leaves at flowering and boll-forming stage were analyzed. The results showed that compared with saline stress, alkaline stress increased the contents of Na⁺, K⁺, SOD, and MDA in leaves. The application of the compound material reduced the content of Na⁺ but increased the K⁺/Na⁺ ratio, the activities of SOD, POD, and CAT, and REC. Transcriptome analysis revealed that after the application of the compound material, the Na⁺/H⁺ exchanger gene in cotton leaves was down-regulated, while the K⁺ transporter, K⁺ channel, and POD genes were up-regulated. Besides, the down-regulation of genes related to lignin synthesis in phenylalanine biosynthesis pathway had a close relationship with the ion content and physiological characteristics in leaves. The quantitative analysis with PCR proved the reliability of the results of RNA sequencing.

CONCLUSION

These findings suggest that the compound material alleviated saline stress and alkaline stress on cotton leaves by regulating candidate genes in key biological pathways, which improves our understanding of the molecular mechanism of the compound material regulating the responses of cotton to saline stress and alkaline stress.

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.

Identification of NHXs in Gossypium species and the positive role of GhNHX1 in salt tolerance

BACKGROUND

Plant Na⁺/H⁺ antiporters (NHXs) are membrane-localized proteins that maintain cellular Na⁺/K⁺ and pH homeostasis. Considerable evidence highlighted the critical roles of NHX family in plant development and salt response; however, NHXs in cotton are rarely studied.

RESULTS

The comprehensive and systematic comparative study of NHXs in three Gossypium species was performed. We identified 12, 12, and 23 putative NHX proteins from G. arboreum, G. raimondii, and G. hirsutum, respectively. Phylogenetic study revealed that repeated polyploidization of Gossypium spp. contributed to the expansion of NHX family. Gene structure analysis showed that cotton NHXs contain many introns, which will lead to alternative splicing and help plants to adapt to high salt concentrations in soil. The expression changes of NHXs indicate the possible differences in the roles of distinct NHXs in salt response. GhNHX1 was proved to be located in the vacuolar system and intensively induced by salt stress in cotton. Silencing of GhNHX1 resulted in enhanced sensitivity of cotton seedlings to high salt concentrations, which suggests that GhNHX1 positively regulates cotton tolerance to salt stress.

CONCLUSION

We characterized the gene structure, phylogenetic relationship, chromosomal location, and expression pattern of NHX genes from G. arboreum, G. raimondii, and G. hirsutum. Our findings indicated that the cotton NHX genes are regulated meticulously and differently at the transcription level with possible alternative splicing. The tolerance of plants to salt stress may rely on the expression level of a particular NHX, rather than the number of NHXs in the genome. This study could provide significant insights into the function of plant NHXs, as well as propose promising candidate genes for breeding salt-resistant cotton cultivars.

Insights into Drought Stress Signaling in Plants and the Molecular Genetic Basis of Cotton Drought Tolerance

Drought stress restricts plant growth and development by altering metabolic activity and biological functions. However, plants have evolved several cellular and molecular mechanisms to overcome drought stress. Drought tolerance is a multiplex trait involving the activation of signaling mechanisms and differentially expressed molecular responses. Broadly, drought tolerance comprises two steps: stress sensing/signaling and activation of various parallel stress responses (including physiological, molecular, and biochemical mechanisms) in plants. At the cellular level, drought induces oxidative stress by overproduction of reactive oxygen species (ROS), ultimately causing the cell membrane to rupture and stimulating various stress signaling pathways (ROS, mitogen-activated-protein-kinase, Ca2+, and hormone-mediated signaling). Drought-induced transcription factors activation and abscisic acid concentration co-ordinate the stress signaling and responses in cotton. The key responses against drought stress, are root development, stomatal closure, photosynthesis, hormone production, and ROS scavenging. The genetic basis, quantitative trait loci and genes of cotton drought tolerance are presented as examples of genetic resources in plants. Sustainable genetic improvements could be achieved through functional genomic approaches and genome modification techniques such as the CRISPR/Cas9 system aid the characterization of genes, sorted out from stress-related candidate single nucleotide polymorphisms, quantitative trait loci, and genes. Exploration of the genetic basis for superior candidate genes linked to stress physiology can be facilitated by integrated functional genomic approaches. We propose a third-generation sequencing approach coupled with genome-wide studies and functional genomic tools, including a comparative sequenced data (transcriptomics, proteomics, and epigenomic) analysis, which offer a platform to identify and characterize novel genes. This will provide information for better understanding the complex stress cellular biology of plants.

Flavonoid accumulation in spontaneous cotton mutant results in red coloration and enhanced disease resistance

Cotton, the leading natural fiber, is cultivated worldwide, but its production is seriously threatened by pathogens. Accordingly, the selection of resistant cultivars has become a key priority of cotton breeding programs. In this study, a spontaneous mutant with red coloration (S156) and a control cultivar (S78) were used as experimental materials for a comparative analysis. Metabolomic analysis revealed the enrichment of flavonoids in S156 leaves compared with S78 leaves, and transcriptomic analysis revealed the upregulated expression of flavonoid biosynthesis genes in S156 leaves relative to S78 leaves. In addition, the red mutant showed a significantly increase in resistance to Verticillium dahliae, a fungal pathogen that poses a major threat to cotton production. The pathogen invasion process was suppressed in the red cotton cultivar. This study reveals the mechanism underlying the red coloration of S156 cotton and indicates the great potential of red cotton in pathogen- and insect-resistant breeding of cotton.

Overexpression of the Key Virulence Factor 1,3-1,4-β-d-Glucanase in the Endophytic Bacterium Bacillus halotolerans Y6 To Improve Verticillium Resistance in Cotton

Verticillium wilt, caused by Verticillium dahliae, results in a dramatic loss of cotton yields in China. There is great potential for biocontrol to manage this destructive crop disease. In this study, we obtained the endophytic bacterium Bacillus halotolerans Y6 from Verticillium wilt-resistant cotton Gossypium barbadense Xinhai15; this bacterium possesses strong antagonistic abilities that inhibit V. dahliae spore germination and mycelial growth. The results of the enzyme activity assay, heterologous expression, and gene knockdown showed that the key virulence factor of Y6 for antagonizing V. dahliae was β -glucanase Bgy6. To facilitate field tests of biological control, we constructed the homologous Bgy6-overexpression strain OY6. Compared with the wild-type Y6 strain, the β-glucanase activity of OY6 was increased by 91.79%, and the inhibition rate of OY6 against V. dahliae V991 exceeded 96.7%. Moreover, the spores of V. dahliae V991 treated with OY6 showed more mucus and larger holes on the surface, as observed by scanning electron microscopy. Potting test results illustrated that both OY6 and Y6 could improve the resistance of upland cotton to Verticillium wilt. With the inoculation of V. dahliae V991 for 45 days, the disease index of G. hirsutum TM-1 treated with OY6 was only 8.33, which was significantly lower than that in plants treated with the wild-type strain Y6 (17.86) or the controls without bacteria (35.94). Our research provides a new idea for the control of Verticillium wilt in upland cotton via transforming endophytic bacteria of Verticillium wilt-resistant cotton and proposes a new solution to prevent and control Verticillium wilt.

Evolutionary Conservation and Divergence of Genes Encoding 3-Hydroxy-3-methylglutaryl Coenzyme A Synthase in the Allotetraploid Cotton Species Gossypium hirsutum

Polyploidization is important for the speciation and subsequent evolution of many plant species. Analyses of the duplicated genes produced via polyploidization events may clarify the origin and evolution of gene families. During terpene biosynthesis, 3-hydroxy-3-methylglutaryl coenzyme A synthase (HMGS) functions as a key enzyme in the mevalonate pathway. In this study, we first identified a total of 53 HMGS genes in 23 land plant species, while no HMGS genes were detected in three green algae species. The phylogenetic analysis suggested that plant HMGS genes may have originated from a common ancestral gene before clustering in different branches during the divergence of plant lineages. Then, we detected six HMGS genes in the allotetraploid cotton species (Gossypium hirsutum), which was twice that of the two diploid cotton species (Gossypium raimondii and Gossypium arboreum). The comparison of gene structures and phylogenetic analysis of HMGS genes revealed conserved evolution during polyploidization in Gossypium. Moreover, the expression patterns indicated that six GhHMGS genes were expressed in all tested tissues, with most genes considerably expressed in the roots, and they were responsive to various phytohormone treatments and abiotic stresses. The sequence and expression divergence of duplicated genes in G. hirsutum implied the sub-functionalization of GhHMGS1A and GhHMGS1D as well as GhHMGS3A and GhHMGS3D, whereas it implied the pseudogenization of GhHMGS2A and GhHMGS2D. Collectively, our study unraveled the evolutionary history of HMGS genes in green plants and from diploid to allotetraploid in cotton and illustrated the different evolutionary fates of duplicated HMGS genes resulting from polyploidization.

Transcriptome analysis reveals differentially expressed ERF transcription factors associated with salt response in cotton

Salinity is a major abiotic stress limiting plant growth and development that has caused severe damage to yield and quality of cotton fiber. Uncovering the mechanisms of response to salt stress is important in breeding salt-tolerant cotton varieties. Transcriptome analysis identified 2356 differentially expressed genes in cotton under salt stress, of which 9.4% were predicted transcription factors (TFs). Approximately 17.6% (39 out of 222) of the differentially expressed TFs belonged to the ethylene response factor (ERF) family. Expression pattern analysis showed significant changes in these ERFs during salt stress. Moreover, the number of down-regulated ERFs was more than that of the up-regulated ERFs. Two of the ERFs, GhERF4L and GhERF54L, showed increased (12-15 times) expression after 12 h of salt treatment. Silencing of GhERF4L and GhERF54L significantly reduced salt tolerance of cotton seedlings, indicating their role in regulating cotton response to salt stress. This study revealed the essential role of ERF transcription factors in the salt response mechanism of plants, and provided important genetic resources for breeding salt-tolerant cotton.

Coping with drought: stress and adaptive mechanisms, and management through cultural and molecular alternatives in cotton as vital constituents for plant stress resilience and fitness

Increased levels of greenhouse gases in the atmosphere and associated climatic variability is primarily responsible for inducing heat waves, flooding and drought stress. Among these, water scarcity is a major limitation to crop productivity. Water stress can severely reduce crop yield and both the severity and duration of the stress are critical. Water availability is a key driver for sustainable cotton production and its limitations can adversely affect physiological and biochemical processes of plants, leading towards lint yield reduction. Adaptation of crop husbandry techniques suitable for cotton crop requires a sound understanding of environmental factors, influencing cotton lint yield and fiber quality. Various defense mechanisms e.g. maintenance of membrane stability, carbon fixation rate, hormone regulation, generation of antioxidants and induction of stress proteins have been found play a vital role in plant survival under moisture stress. Plant molecular breeding plays a functional role to ascertain superior genes for important traits and can offer breeder ready markers for developing ideotypes. This review highlights drought-induced damage to cotton plants at structural, physiological and molecular levels. It also discusses the opportunities for increasing drought tolerance in cotton either through modern gene editing technology like clustered regularly interspaced short palindromic repeat (CRISPR/Cas9), zinc finger nuclease, molecular breeding as well as through crop management, such as use of appropriate fertilization, growth regulator application and soil amendments.

Genome-Wide Identification and Comparative Analysis of the 3-Hydroxy-3-methylglutaryl Coenzyme A Reductase (HMGR) Gene Family in Gossypium

Terpenes are the largest and most diverse class of secondary metabolites in plants and play a very important role in plant adaptation to environment. 3-Hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-limiting enzyme in the process of terpene biosynthesis in the cytosol. Previous study found the HMGR genes underwent gene expansion in Gossypium raimondii, but the characteristics and evolution of the HMGR gene family in Gossypium genus are unclear. In this study, genome-wide identification and comparative study of HMGR gene family were carried out in three Gossypium species with genome sequences, i.e., G. raimondii, Gossypium arboreum, and Gossypium hirsutum. In total, nine, nine and 18 HMGR genes were identified in G. raimondii, G. arboreum, and G. hirsutum, respectively. The results indicated that the HMGR genes underwent gene expansion and a unique gene cluster containing four HMGR genes was found in all the three Gossypium species. The phylogenetic analysis suggested that the expansion of HMGR genes had occurred in their common ancestor. There was a pseudogene that had a 10-bp deletion resulting in a frameshift mutation and could not be translated into functional proteins in G. arboreum and the A-subgenome of G. hirsutum. The expression profiles of the two pseudogenes showed that they had tissue-specific expression. Additionally, the expression pattern of the pseudogene in the A-subgenome of G. hirsutum was similar to its paralogous gene in the D-subgenome of G. hirsutum. Our results provide useful information for understanding cytosolic terpene biosynthesis in Gossypium species.

Genome-Wide Association Studies Reveal Genetic Variation and Candidate Genes of Drought Stress Related Traits in Cotton (Gossypium hirsutum L.)

Cotton is an important industrial crop worldwide and upland cotton (Gossypium hirsutum L.) is most widely cultivated in the world. Due to ever-increasing water deficit, drought stress brings a major threat to cotton production. Thus, it is important to reveal the genetic basis under drought stress and develop drought tolerant cotton cultivars. To address this issue, in present study, 319 upland cotton accessions were genotyped by 55,060 single nucleotide polymorphisms (SNPs) from high-density CottonSNP80K array and phenotyped nine drought tolerance related traits. The two datasets were used to identify quantitative trait nucleotides (QTNs) for the above nine traits using multi-locus random-SNP-effect mixed linear model method. As a result, a total of 20 QTNs distributed on 16 chromosomes were found to be significantly associated with six drought tolerance related traits. Of the 1,326 genes around the 20 QTNs, 205 were induced after drought stress treatment, and 46 were further mapped to Gene ontology (GO) term "response to stress." Taken genome-wide association study (GWAS) analysis, RNA-seq data and qRT-PCR verification, four genes, RD2 encoding a response to desiccation 2 protein, HAT22 encoding a homeobox-leucine zipper protein, PIP2 encoding a plasma membrane intrinsic protein 2, and PP2C encoding a protein phosphatase 2C, were proposed to be potentially important for drought tolerance in cotton. These results will deepen our understanding of the genetic basis of drought stress tolerance in cotton and provide candidate markers to accelerate the development of drought-tolerant cotton cultivars.

GhSNAP33, a t-SNARE Protein From Gossypium hirsutum, Mediates Resistance to Verticillium dahliae Infection and Tolerance to Drought Stress

Soluble N-ethylmaleimide-sensitive fusion protein attachment protein receptor (SNARE) proteins mediate membrane fusion and deliver cargo to specific cellular locations through vesicle trafficking. Synaptosome-associated protein of 25 kDa (SNAP25) is a target membrane SNARE that drives exocytosis by fusing plasma and vesicular membranes. In this study, we isolated GhSNAP33, a gene from cotton (Gossypium hirsutum), encoding a SNAP25-type protein containing glutamine (Q)b- and Qc-SNARE motifs connected by a linker. GhSNAP33 expression was induced by H₂O₂, salicylic acid, abscisic acid, and polyethylene glycol 6000 treatment and Verticillium dahliae inoculation. Ectopic expression of GhSNAP33 enhanced the tolerance of yeast cells to oxidative and osmotic stresses. Virus-induced gene silencing of GhSNAP33 induced spontaneous cell death and reactive oxygen species accumulation in true leaves at a later stage of cotton development. GhSNAP33-deficient cotton was susceptible to V. dahliae infection, which resulted in severe wilt on leaves, an elevated disease index, enhanced vascular browning and thylose accumulation. Conversely, Arabidopsis plants overexpressing GhSNAP33 showed significant resistance to V. dahliae, with reduced disease index and fungal biomass and elevated expression of PR1 and PR5. Leaves from GhSNAP33-transgenic plants showed increased callose deposition and reduced mycelia growth. Moreover, GhSNAP33 overexpression enhanced drought tolerance in Arabidopsis, accompanied with reduced water loss rate and enhanced expression of DERB2A and RD29A during dehydration. Thus, GhSNAP33 positively mediates plant defense against stress conditions and V. dahliae infection, rendering it a candidate for the generation of stress-resistant engineered cotton.