Genetic basis for glandular trichome formation in cotton

Terpene synthases GhTPS6 and GhTPS47 participate in resistance to Verticillium dahliae in upland cotton

Terpene synthases (TPSs) are enzymes responsible for catalyzing the production of diverse terpenes, the largest class of secondary metabolites in plants. Here, we identified 107 TPS gene loci encompassing 92 full-length TPS genes in upland cotton (Gossypium hirsutum L.). Phylogenetic analysis showed they were divided into six subfamilies. Segmental duplication and tandem duplication events contributed greatly to the expansion of TPS gene family, particularly the TPS-a and TPS-b subfamilies. Expression profile analysis screened out that GhTPSs may mediate the interaction between cotton and Verticillium dahliae. Three-dimensional structures and subcellular localizations of the two selected GhTPSs, GhTPS6 and GhTPS47, which belong to the TPS-a subfamily, demonstrated similarity in protein structures and nucleus and cytoplasm localization. Virus-induced gene silencing (VIGS) of the two GhTPSs yielded plants characterized by increased wilting and chlorosis, more severe vascular browning, and higher disease index than control plants. Additionally, knockdown of GhTPS6 and GhTPS47 led to the down-regulation of cotton terpene synthesis following V. dahliae infection, indicating that these two genes may positively regulate resistance to V. dahliae through the modulation of disease-resistant terpene biosynthesis. Overall, our study represents a comprehensive analysis of the G. hirsutum TPS gene family, revealing their potential roles in defense responses against Verticillium wilt.

GhHAM regulates GoPGF-dependent gland development and contributes to broad-spectrum pest resistance in cotton

Cotton is a globally cultivated crop, producing 87% of the natural fiber used in the global textile industry. The pigment glands, unique to cotton and its relatives, serve as a defense structure against pests and pathogens. However, the molecular mechanism underlying gland formation and the specific role of pigment glands in cotton's pest defense are still not well understood. In this study, we cloned a gland-related transcription factor GhHAM and generated the GhHAM knockout mutant using CRISPR/Cas9. Phenotypic observations, transcriptome analysis, and promoter-binding experiments revealed that GhHAM binds to the promoter of GoPGF, regulating pigment gland formation in cotton's multiple organs via the GoPGF-GhJUB1 module. The knockout of GhHAM significantly reduced gossypol production and increased cotton's susceptibility to pests in the field. Feeding assays demonstrated that more than 80% of the cotton bollworm larvae preferred ghham over the wild type. Furthermore, the ghham mutants displayed shorter cell length and decreased gibberellins (GA) production in the stem. Exogenous application of GA3 restored stem cell elongation but not gland formation, thereby indicating that GhHAM controls gland morphogenesis independently of GA. Our study sheds light on the functional differentiation of HAM proteins among plant species, highlights the significant role of pigment glands in influencing pest feeding preference, and provides a theoretical basis for breeding pest-resistant cotton varieties to address the challenges posed by frequent outbreaks of pests.

Regulation of Glandular Size and Phytoalexin Biosynthesis by a Negative Feedback Loop in Cotton

Plants have evolved diverse defense mechanisms encompassing physical and chemical barriers. Cotton pigment glands are known for containing various defense metabolites, but the precise regulation of gland size to modulate defense compound levels remains enigmatic. Here, it is discovered that the VQ domain-containing protein JAVL negatively regulates pigment gland size and the biosynthesis of defense compounds, while the MYC2-like transcription factor GoPGF has the opposite effect. Notably, GoPGF directly activates the expression of JAVL, whereas JAVL suppresses GoPGF transcription, establishing a negative feedback loop that maintains the expression homeostasis between GoPGF and JAVL. Furthermore, it is observed that JAVL negatively regulates jasmonate levels by inhibiting the expression of jasmonate biosynthetic genes and interacting with GoPGF to attenuate its activation effects, thereby maintaining homeostatic regulation of jasmonate levels. The increased expression ratio of GoPGF to JAVL leads to enlarged pigment glands and elevated jasmonates and defense compounds, enhancing insect and pathogen resistance in cotton. These findings unveil a new mechanism for regulating gland size and secondary metabolites biosynthesis, providing innovative strategies for strengthening plant defense.

Deep learning-based quantification and transcriptomic profiling reveal a methyl jasmonate-mediated glandular trichome formation pathway in Cannabis sativa

Cannabis glandular trichomes (GTs) are economically and biotechnologically important structures that have a remarkable morphology and capacity to produce, store, and secrete diverse classes of secondary metabolites. However, our understanding of the developmental changes and the underlying molecular processes involved in cannabis GT development is limited. In this study, we developed Cannabis Glandular Trichome Detection Model (CGTDM), a deep learning-based model capable of differentiating and quantifying three types of cannabis GTs with a high degree of efficiency and accuracy. By profiling at eight different time points, we captured dynamic changes in gene expression, phenotypes, and metabolic processes associated with GT development. By integrating weighted gene co-expression network analysis with CGTDM measurements, we established correlations between phenotypic variations in GT traits and the global transcriptome profiles across the developmental gradient. Notably, we identified a module containing methyl jasmonate (MeJA)-responsive genes that significantly correlated with stalked GT density and cannabinoid content during development, suggesting the existence of a MeJA-mediated GT formation pathway. Our findings were further supported by the successful promotion of GT development in cannabis through exogenous MeJA treatment. Importantly, we have identified CsMYC4 as a key transcription factor that positively regulates GT formation via MeJA signaling in cannabis. These findings provide novel tools for GT detection and counting, as well as valuable information for understanding the molecular regulatory mechanism of GT formation, which has the potential to facilitate the molecular breeding, targeted engineering, informed harvest timing, and manipulation of cannabinoid production.

Molecular regulation of oil gland development and biosynthesis of essential oils in Citrus spp

Secretory structures in terrestrial plants serve as reservoirs for a variety of secondary metabolites. Among these, the secretory cavity of the Rutaceae family is notable for containing essential oils with a wide range of applications. However, the molecular basis underlying secretory cavity development is unknown. Here, we reveal a molecular framework for Citrus oil gland formation. Using genetic mapping and genome editing, we demonstrated that this process requires LATE MERISTEM IDENTITY1 (LMI1), a key regulator of leaf serration. A conserved GCC box element of the LMI1 promoter recruits DORNROSCHEN-like (DRNL) for transcriptional activation. This DRNL-LMI1 cascade triggers MYC5 activation, facilitating the development of oil glands and the biosynthesis of essential oils. Our findings spotlight cis-regulatory divergence within leaf shape genes, propelling novel functional tissue formation.

Overexpression of cucumber CYP82D47 enhances resistance to powdery mildew and Fusarium oxysporum f. sp. cucumerinum

Cytochrome P450s are a large family of protein-encoding genes in plant genomes, many of which have not yet been comprehensively characterized. Here, a novel P450 gene, CYP82D47, was isolated and functionally characterized from cucumber (Cucumis sativus L.). Quantitative real-time reverse-transcription polymerase chain reaction analysis revealed that CYP82D47 expression was triggered by salicylic acid (SA) and ethephon (ETH). Expression analysis revealed a correlation between CYP82D47 transcript levels and plant defense responses against powdery mildew (PM) and Fusarium oxysporum f. sp. cucumerinum (Foc). Although no significant differences were observed in disease resistance between CYP82D47-RNAi and wild-type cucumber, overexpression (OE) of CYP82D47 enhanced PM and Foc resistance in cucumber. Furthermore, the expression levels of SA-related genes (PR1, PR2, PR4, and PR5) increased in CYP82D47-overexpressing plants 7 days post fungal inoculation. The levels of ETH-related genes (EIN3 and EBF2) were similarly upregulated. The observed enhanced resistance was associated with the upregulation of SA/ETH-signaling-dependent defense genes. These findings indicate the crucial role of CYP82D47 in pathogen defense in cucumber. CYP82D47-overexpressing cucumber plants exhibited heightened susceptibility to both diseases. The study results offer important insights that could aid in the development of disease-resistant cucumber cultivars and elucidate the molecular mechanisms associated with the functions of CYP82D47.

Applications of Virus-Induced Gene Silencing in Cotton

Virus-induced gene silencing (VIGS) is an RNA-mediated reverse genetics technique that has become an effective tool to investigate gene function in plants. Cotton is one of the most important economic crops globally. In the past decade, VIGS has been successfully applied in cotton functional genomic studies, including those examining abiotic and biotic stress responses and vegetative and reproductive development. This article summarizes the traditional vectors used in the cotton VIGS system, the visible markers used for endogenous gene silencing, the applications of VIGS in cotton functional genomics, and the limitations of VIGS and how they can be addressed in cotton.

The dual role of GoPGF reveals that the pigment glands are synthetic sites of gossypol in aerial parts of cotton

Gossypol and the related terpenoids are stored in the pigment gland to protect cotton plants from biotic stresses, but little is known about the synthetic sites of these metabolites. Here, we showed that GoPGF, a key gene regulating gland formation, was expressed in gland cells and roots. The chromatin immunoprecipitation sequencing (ChIP-seq) analysis demonstrated that GoPGF targets GhJUB1 to regulate gland morphogenesis. RNA-sequencing (RNA-seq) showed high accumulation of gossypol biosynthetic genes in gland cells. Moreover, integrated analysis of the ChIP-seq and RNA-seq data revealed that GoPGF binds to the promoter of several gossypol biosynthetic genes. The cotton callus overexpressing GoPGF had dramatically increased the gossypol levels, indicating that GoPGF can directly activate the biosynthesis of gossypol. In addition, the gopgf mutant analysis revealed the existence of both GoPGF-dependent and -independent regulation of gossypol production in cotton roots. Our study revealed that the pigment glands are synthetic sites of gossypol in aerial parts of cotton and that GoPGF plays a dual role in regulating gland morphogenesis and gossypol biosynthesis. The study provides new insights for exploring the complex relationship between glands and the metabolites they store in cotton and other plant species.

QTL Mapping of Trichome Traits and Analysis of Candidate Genes in Leaves of Wheat (Triticum aestivum L.)

Trichome plays an important role in heat dissipation, cold resistance, water absorption, protection of leaves from mechanical damage, and direct exposure to ultraviolet rays. It also plays an important role in the photosynthesis, transpiration, and respiration of plants. However, the genetic basis of trichome traits is not fully understood in wheat. In this study, wheat DH population (Hanxuan 10 × Lumai 14) was used to map quantitative trait loci (QTL) for trichome traits in different parts of flag leaf at 10 days after anther with growing in Zhao County, Hebei Province, and Taigu County, Shanxi Province, respectively. The results showed that trichome density (TD) was leaf center > leaf tip > leaf base and near vein > middle > edge, respectively, in both environments. The trichome length (TL) was leaf tip > leaf center > leaf base and edge > middle > near vein. Significant phenotypic positive correlations were observed between the trichome-related traits of different parts. A total of 83 QTLs for trichome-related traits were mapped onto 18 chromosomes, and each one accounted for 2.41 to 27.99% of the phenotypic variations. Two QTL hotspots were detected in two marker intervals: AX-95232910~AX-95658735 on 3A and AX-94850949~AX-109507404 on 7D. Six possible candidate genes (TraesCS3A02G406000, TraesCS3A02G414900, TraesCS3A02G440900, TraesCS7D02G145200, TraesCS7D02G149200, and TraesCS7D02G152400) for trichome-related traits of wheat leaves were screened out according to their predicted expression levels in wheat leaves. The expression of these genes may be induced by a variety of abiotic stresses. The results provide the basis for further validation and functional characterization of the candidate genes.

Single-cell RNA sequencing reveals a hierarchical transcriptional regulatory network of terpenoid biosynthesis in cotton secretory glandular cells

Plants can synthesize a wide range of terpenoids in response to various environmental cues. However, the specific regulatory mechanisms governing terpenoid biosynthesis at the cellular level remain largely elusive. In this study, we employed single-cell RNA sequencing to comprehensively characterize the transcriptome profile of cotton leaves and established a hierarchical transcriptional network regulating cell-specific terpenoid production. We observed substantial expression levels of genes associated with the biosynthesis of both volatile terpenes (such as β-caryophyllene and β-myrcene) and non-volatile gossypol-type terpenoids in secretory glandular cells. Moreover, two novel transcription factors, namely GoHSFA4a and GoNAC42, are identified to function downstream of the Gossypium PIGMENT GLAND FORMATION genes. Both transcription factors could directly regulate the expression of terpenoid biosynthetic genes in secretory glandular cells in response to developmental and environmental stimuli. For convenient retrieval of the single-cell RNA sequencing data generated in this study, we developed a user-friendly web server . Our findings not only offer valuable insights into the precise regulation of terpenoid biosynthesis genes in cotton leaves but also provide potential targets for cotton breeding endeavors.

Characterization and fine mapping of a yellow leaf gene regulating chlorophyll biosynthesis and chloroplast development in cotton (Gossypium arboreum)

Chlorophyll biosynthesis and chloroplast development are essential for photosynthesis and plant growth. Gossypium arboreum, a valuable source of genetic variation for cotton improvement, remains poorly studied for the mechanisms regulating chlorophyll biosynthesis and chloroplast development. Here we created a G. arboreum etiolated leaf and stuntedness (els) mutant that displayed a distinct yellow color of leaves, bracts and stems throughout the whole growth, where chlorophyll accumulation in leaves was reduced and chloroplast development was delayed. The GaCHLH gene, which encodes the H subunit of magnesium chelatase (Mg-chelatase), was screened by MutMap and KASP analysis. Compared to GaCHLH, the gene Gachlh of the mutant had a single nucleotide transition (G to A) at 1549 bp, which causes the substitution of a glycine (G) by a serine (S) at the 517th amino acid, resulting in an abnormal secondary structure of the Gachlh protein. GaCHLH-silenced SXY1 and ZM24 plants exhibited a lower GaCHLH expression level, a lower chlorophyll content, and the yellow-leaf phenotype. Gachlh expression affected the expression of key genes in the tetrapyrrole pathway. GaCHLH and Gachlh were located in the chloroplasts and that alteration of the mutation site did not affect the final target position. The BiFC assay result indicated that Gachlh could not bind to GaCHLD properly, which prevented the assembly of Mg-chelatase and thus led to the failure of chlorophyll synthesis. In this study, the Gachlh gene of G. arboreum els was finely localized and identified for the first time, providing new insights into the chlorophyll biosynthesis pathway in cotton.

GoSTR, a negative modulator of stem trichome formation in cotton

Trichomes, the outward projection of plant epidermal tissue, provide an effective defense against stress and insect pests. Although numerous genes have been identified to be involved in trichome development, the molecular mechanism for trichome cell fate determination is not well enunciated. Here, we reported GoSTR functions as a master repressor for stem trichome formation, which was isolated by map-based cloning based on a large F2 segregating population derived from a cross between TM-1 (pubescent stem) and J220 (smooth stem). Sequence alignment revealed a critical G-to-T point mutation in GoSTR's coding region that converted codon 2 from GCA (Alanine) to TCA (Serine). This mutation occurred between the majority of Gossypium hirsutum with pubescent stem (GG-haplotype) and G. barbadense with glabrous stem (TT-haplotype). Silencing of GoSTR in J220 and Hai7124 via virus-induced gene silencing resulted in the pubescent stems but no visible change in leaf trichomes, suggesting stem trichomes and leaf trichomes are genetically distinct. Yeast two-hybrid assay and luciferase complementation imaging assay showed GoSTR interacts with GoHD1 and GoHOX3, two key regulators of trichome development. Comparative transcriptomic analysis further indicated that many transcription factors such as GhMYB109, GhTTG1, and GhMYC1/GhDEL65 which function as positive regulators of trichomes were significantly upregulated in the stem from the GoSTR-silencing plant. Taken together, these results indicate that GoSTR functions as an essential negative modulator of stem trichomes and its transcripts will greatly repress trichome cell differentiation and growth. This study provided valuable insights for plant epidermal hair initiation and differentiation research.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Genome-wide identification and expression analysis of terpene synthases in Gossypium species in response to gossypol biosynthesis

Cottonseed is an invaluable resource, providing protein, oil, and abundant minerals that significantly contribute to the well-being and nutritional needs of both humans and livestock. However, cottonseed also contains a toxic substance called gossypol, a secondary metabolite in Gossypium species that plays an important role in cotton plant development and self-protection. Herein, genome-wide analysis and characterization of the terpene synthase (TPS) gene family identified 304 TPS genes in Gossypium. Bioinformatics analysis revealed that the gene family was grouped into six subgroups TPS-a, TPS-b, TPS-c, TPS-e, TPS-f, and TPS-g. Whole-genome, segmental, and tandem duplication contributed to the evolution of TPS genes. According to the analysis of selection pressure, it was predicted that TPS genes experience predominantly negative selection, with positive selection occurring subsequently. RT-qPCR analysis in TM-1 and CRI-12 lines revealed GhTPS48 gene as the candidate gene for silencing experiments. To summarize, comprehensive genome-wide studies, RT-qPCR, and gene silencing experiments have collectively demonstrated the involvement of the TPS gene family in the biosynthesis of gossypol in cotton.

AaSEPALLATA1 integrates jasmonate and light-regulated glandular secretory trichome initiation in Artemisia annua

Glandular secretory trichomes (GSTs) can secrete and store a variety of specific metabolites. By increasing GST density, valuable metabolites can be enhanced in terms of productivity. However, the comprehensive and detailed regulatory network of GST initiation still needs further investigation. By screening a complementary DNA library derived from young leaves of Artemisia annua, we identified a MADS-box transcription factor, AaSEPALLATA1 (AaSEP1), that positively regulates GST initiation. Overexpression of AaSEP1 in A. annua substantially increased GST density and artemisinin content. The HOMEODOMAIN PROTEIN 1 (AaHD1)-AaMYB16 regulatory network regulates GST initiation via the jasmonate (JA) signaling pathway. In this study, AaSEP1 enhanced the function of AaHD1 activation on downstream GST initiation gene GLANDULAR TRICHOME-SPECIFIC WRKY 2 (AaGSW2) through interaction with AaMYB16. Moreover, AaSEP1 interacted with the JA ZIM-domain 8 (AaJAZ8) and served as an important factor in JA-mediated GST initiation. We also found that AaSEP1 interacted with CONSTITUTIVE PHOTOMORPHOGENIC 1 (AaCOP1), a major repressor of light signaling. In this study, we identified a MADS-box transcription factor that is induced by JA and light signaling and that promotes the initiation of GST in A. annua.

Comparative Transcriptome Analysis Revealed Key Genes Regulating Gossypol Synthesis in Tetraploid Cultivated Cotton

Tetraploid cultivated cotton (Gossypium spp.) produces cottonseeds rich in protein and oil. Gossypol and related terpenoids, stored in the pigment glands of cottonseeds, are toxic to human beings and monogastric animals. However, a comprehensive understanding of the genetic basis of gossypol and gland formation is still lacking. We performed a comprehensive transcriptome analysis of four glanded versus two glandless tetraploid cultivars distributed in Gossypium hirsutum and Gossypium barbadense. A weighted gene co-expression network analysis (WGCNA) based on 431 common differentially expressed genes (DEGs) uncovered a candidate module that was strongly associated with the reduction in or disappearance of gossypol and pigment glands. Further, the co-expression network helped us to focus on 29 hub genes, which played key roles in the regulation of related genes in the candidate module. The present study contributes to our understanding of the genetic basis of gossypol and gland formation and serves as a rich potential source for breeding cotton cultivars with gossypol-rich plants and gossypol-free cottonseed, which is beneficial for improving food safety, environmental protection, and economic gains of tetraploid cultivated cotton.

Genome-wide identification of terpenoid synthase family genes in Gossypium hirsutum and functional dissection of its subfamily cadinene synthase A in gossypol synthesis

Plant terpenoid synthase (TPS) family genes participate in metabolite synthesis, hormones, gossypol, etc. Here, we genome-widely identified TPS family genes in 12 land plant species. Four hundred and thirty TPS-related genes were divided into seven subfamilies. The TPS-c in Bryophytes was suggested to be the earliest subfamily, followed by the TPS-e/f and TPS-h presence in ferns. TPS-a, the largest number of genes, was derived from monocotyledonous and dicotyledonous plants. Collinearity analysis showed that 38 out of the 76 TPS genes in G. hirsutum were collinear within G. arboreum and G. raimondii. Twenty-one GhTPS-a genes belong to the cadinene synthase (GhCDN) subfamily and were divided into five groups, A, B, C, D, and E. The special cis-elements in the promoters of 12 GhCDN-A genes suggested that the JA and ethylene signaling pathways may be involved in their expression regulation. When 12 GhCDN-A genes were simultaneously silenced through virus-induced gene silencing, the glandular color of GhCDN-A-silenced plants was lighter than that of the control, supported by a gossypol content decrease based on HPLC testing, suggesting that GhCDN-A subgroup genes participate in gossypol synthesis. According to RNA-seq analysis, gossypol synthesis-related genes and disease-resistant genes in the glandular variety exhibited upregulated expression compared to the glandless variety, whereas hormone signaling-related genes were downregulated. All in all, these results revealed plant TPS gene evolution rules and dissected the TPS subfamily, GhCDN-A, function in gossypol synthesis in cotton.

Identification of MYB gene family and functional analysis of GhMYB4 in cotton (Gossypium spp.)

Myeloblastosis (MYB) transcription factors (TFs) form a large gene family involved in a variety of biological processes in plants. Little is known about their roles in the development of cotton pigment glands. In this study, 646 MYB members were identified in Gossypium hirsutum genome and phylogenetic classification was analyzed. Evolution analysis revealed assymetric evolution of GhMYBs during polyploidization and sequence divergence of MYBs in G. hirustum was preferentially happend in D sub-genome. WGCNA (weighted gene co-expression network analysis) showed that four modules had potential relationship with gland development or gossypol biosynthesis in cotton. Eight differentially expressed GhMYB genes were identified by screening transcriptome data of three pairs of glanded and glandless cotton lines. Of these, four were selected as candidate genes for cotton pigment gland formation or gossypol biosynthesis by qRT-PCR assay. Silencing of GH_A11G1361 (GhMYB4) downregulated expression of multiple genes in gossypol biosynthesis pathway, indicating it could be involved in gossypol biosynthesis. The potential protein interaction network suggests that several MYBs may have indirect interaction with GhMYC2-like, a key regulator of pigment gland formation. Our study was the systematic analysis of MYB genes in cotton pigment gland development, providing candidate genes for further study on the roles of cotton MYB genes in pigment gland formation, gossypol biosynthesis and future crop plant improvement.

Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of pigment glands in Gossypium bickii

Comprehensive utilization of cottonseeds is limited by the presence of pigment glands and its inclusion gossypol. The ideal cotton has glandless seeds but a glanded plant, a trait found in only a few Australian wild cotton species, including Gossypium bickii. Introgression of this trait into cultivated species has proved to be difficult. Understanding the biological processes toward pigment gland morphogenesis and the associated underlying molecular mechanisms will facilitate breeding of cultivated cotton varieties with the trait of glandless seeds and glanded plant. In this study, single-cell RNA sequencing (scRNA-seq) was performed on 12 222 protoplasts isolated from cotyledons of germinating G. bickii seeds 48 h after imbibition. Clustered into 14 distinct clusters unsupervisedly, these cells could be grouped into eight cell populations with the assistance of known cell marker genes. The pigment gland cells were well separated from others and could be separated into pigment gland parenchyma cells, secretory cells, and apoptotic cells. By integrating the pigment gland cell developmental trajectory, transcription factor regulatory networks, and core transcription factor functional validation, we established a model for pigment gland formation. In this model, light and gibberellin were verified to promote the formation of pigment glands. In addition, three novel genes, GbiERF114 (ETHYLENE RESPONSE FACTOR 114), GbiZAT11 (ZINC FINGER OF ARABIDOPSIS THALIANA 11), and GbiNTL9 (NAC TRANSCRIPTION FACTOR-LIKE 9), were found to affect pigment gland formation. Collectively, these findings provide new insights into pigment gland morphogenesis and lay the cornerstone for future cotton scRNA-seq investigations.

Dirigent gene editing of gossypol enantiomers for toxicity-depleted cotton seeds

Axial chirality of biaryls can generate varied bioactivities. Gossypol is a binaphthyl compound made by cotton plants. Of its two axially chiral isomers, (-)-gossypol is the bioactive form in mammals and has antispermatogenic activity, and its accumulation in cotton seeds poses health concerns. Here we identified two extracellular dirigent proteins (DIRs) from Gossypium hirsutum, GhDIR5 and GhDIR6, which impart the hemigossypol oxidative coupling into (-)- and (+)-gossypol, respectively. To reduce cotton seed toxicity, we disrupted GhDIR5 by genome editing, which eliminated (-)-gossypol but had no effects on other phytoalexins, including (+)-gossypol, that provide pest resistance. Reciprocal mutagenesis identified three residues responsible for enantioselectivity. The (-)-gossypol-forming DIRs emerged later than their enantiocomplementary counterparts, from tandem gene duplications that occurred shortly after the cotton genus diverged. Our study offers insight into how plants control enantiomeric ratios and how to selectively modify the chemical spectra of cotton plants and thereby improve crop quality.

Identification and Functional Characterization of CsMYCs in Cucumber Glandular Trichome Development

Glandular trichomes (GTs), specialized structures formed by the differentiation of plant epidermal cells, are known to play important roles in the resistance of plants to external biotic and abiotic stresses. These structures are capable of storing and secreting secondary metabolites, which often have important agricultural and medicinal values. In order to better understand the molecular developmental mechanisms of GTs, studies have been conducted in a variety of crops, including tomato (Solanum lycopersicum), sweetworm (Artemisia annua), and cotton (Gossypium hirsutum). The MYC transcription factor of the basic helix-loop-helix (bHLH) transcription factor family has been found to play an important role in GT development. In this study, a total of 13 cucumber MYC transcription factors were identified in the cucumber (Cucumis sativus L.) genome. After performing phylogenetic analyses and conserved motifs on the 13 CsMYCs in comparison to previously reported MYC transcription factors that regulate trichome development, seven candidate MYC transcription factors were selected. Through virus-induced gene silencing (VIGS), CsMYC2 is found to negatively regulate GT formation while CsMYC4, CsMYC5, CsMYC6, CsMYC7, and CsMYC8 are found to positively regulate GT formation. Furthermore, the two master effector genes, CsMYC2 and CsMYC7, are observed to have similar expression patterns indicating that they co-regulate the balance of GT development in an antagonistic way.

Transcriptome analysis reveals the effect of grafting on gossypol biosynthesis and gland formation in cotton

BACKGROUND

Gossypol is a unique secondary metabolite and sesquiterpene in cotton, which is mainly synthesized in the root system of cotton and exhibits many biological activities. Previous research found that grafting affected the density of pigment glands and the gossypol content in cotton.

RESULTS

This study performed a transcriptome analysis on cotton rootstocks and scions of four grafting methods. The gene expression of mutual grafting and self-grafting was compared to explore the potential genes involved in gossypol biosynthesis. A total of six differentially expressed enzymes were found in the main pathway of gossypol synthesis-sesquiterpene and triterpene biosynthesis (map00909): lupeol synthase (LUP1, EC:5.4.99.41), beta-amyrin synthase (LUP2, EC:5.4.99.39), squalene monooxygenase (SQLE, EC:1.14.14.17), squalene synthase (FDFT1, EC:2.5.1.21), (-)-germacrene D synthase (GERD, EC:4.2.3.75), ( +)-delta-cadinene synthase (CADS, EC:4.2.3.13). By comparing the results of the gossypol content and the density of the pigment gland, we speculated that these six enzymes might affect the biosynthesis of gossypol. It was verified by qRT-PCR analysis that grafting could influence gene expression of scion and stock. After suppressing the expression of the LUP1, FDFT1, and CAD genes by VIGS technology, the gossypol content in plants was significantly down-regulated.

CONCLUSIONS

These results indicate the potential molecular mechanism of gossypol synthesis during the grafting process and provide a theoretical foundation for further research on gossypol biosynthesis.

A male-sterile mutant with necrosis-like dark spots on anthers was generated in cotton

Although conventional hybrid breeding has paved the way for improving cotton production and other properties, it is undoubtedly time and labor consuming, while the cultivation of male sterile line can fix the problem. Here, we induced male sterile mutants by simultaneously editing three cotton EXCESS MICROSPOROCYTES1 (GhEMS1) genes by CRISPR/Cas9. Notably, the GhEMS1 genes are homologous to AtEMS1 genes, which inhibit the production of middle layer and tapetum cells as well, leading to male sterility in cotton. Interestingly, there are necrosis-like dark spots on the surface of the anthers of GhEMS1s mutants, which is different from AtEMS1 mutant whose anther surface is clean and smooth, suggesting that the function of EMS1 gene has not been uncovered yet. Moreover, we have detected mutations in GhEMS1 genes from T₀ to T₃ mutant plants, which had necrosis-like dark spots as well, indicating that the mutation of the three GhEMS1 genes could be stably inherited. Dynamic transcriptomes showed plant hormone pathway and anther development genetic network were differential expression in mutant and wild-type anthers. And the lower level of IAA content in the mutant anthers than that in the wild type at four anther developmental stages may be the reason for the male sterility. This study not only facilitates the exploration of the basic research of cotton male sterile lines, but also provides germplasms for accelerating the hybrid breeding in cotton.

A highly accumulated secretory protein from cotton bollworm interacts with basic helix-loop-helix transcription factors to dampen plant defense

Caterpillar oral secretion (OS) contains active molecules that modulate plant defense signaling. We isolated an effector-like protein (Highly Accumulated Secretory Protein 1, HAS1) from cotton bollworm (Helicoverpa armigera) that is the most highly accumulated secretory protein of the nondigestive components in OS and belongs to venom R-like protein. Elimination of HAS1 by plant-mediated RNA interference reduced the suppression of OS on the defense response in plants. Plants expressing HAS1 are more susceptible to insect herbivory accompanied by the reduced expressions of multiple defense genes. HAS1 binds to the basic helix-loop-helix (bHLH) transcription factors, including GoPGF involved in pigmented gland formation and defense compounds biosynthesis in cotton and MYC3/MYC4 the main regulators in jasmonate (JA) signaling in Arabidopsis. The binding activity is required for HAS1 to inhibit the activation of bHLHs on plant defense gene expressions. Together with our previous study that another venom R-like protein HARP1 in cotton bollworm OS blocks JA signaling by interacting with JASMONATE-ZIM-domain repressors, we conclude that the venom R-like proteins in OS interfere with plant defense in a dual suppression manner. Considering the venom proteins in parasitic wasp assault the immune system of its host animal, our investigation reveals their conserved function in carnivorous and herbivorous insects.

Duplicate mutations of GhCYP450 lead to the production of ms5m6 male sterile line in cotton

The duplicated male sterile genes ms₅m₆ in cotton were map-based cloned and validated by the virus-induced gene silencing assays. Duplicate mutations of the GhCYP450 gene encoding a cytochrome P450 protein are responsible for the male sterility in cotton. The utilization of male sterility in cotton plays a vital role in improving yield and fiber quality. A complete male sterile line (ms₅ms₆) has been extensively used to develop hybrid cotton worldwide. Using Zhongkang-A (ZK-A) developed by transferring Bt and ms₅ms₆ genes into the commercial cultivar Zhongmiansuo 12, the duplicate genes were map-based cloned and confirmed via the virus-induced gene silencing (VIGS) assays. The duplicate mutations of GhCYP450 genes encoding a cytochrome P450 protein were responsible for producing male sterility in ms₅ms₆ in cotton. Sequence alignment showed that GhCYP450-Dt in ZK-A differed in two critical aspects from the fertile wild-type TM-1: GhCYP450-Dt has three amino acid (D98E, E168K, G198R) changes in the coding region and a 7-bp (GGAAAAA) insertion in the promoter domain; GhCYP450-At appears to be premature termination of GhCYP450 translation. Further morphological observation and cytological examination of GhCYP450-silenced plants induced by VIGS exhibited shorter filaments and no mature pollen grains. These results indicate that GhCYP450 is essential for pollen exine formation and pollen development for male fertility. Investigating the mechanisms of ms₅ms₆ male sterility will deepen our understanding of the development and utilization of heterosis.

High-resolution sequencing of nine elite upland cotton cultivars uncovers genic variations and breeding improvement targets

Structural variations (SVs) are critical factors affecting genome evolution and important traits. However, identification results and functional analyses of SVs in upland cotton are rare. Here, based on the genetic relationships, breeding history and cumulative planting area of upland cotton in China, nine predominant cultivars from the past 60 years (1950s-2010s) were selected for long read sequencing to uncover genic variations and breeding improvement targets for this crop. Based on the ZM24 reference genome, 0.88-1.47 × 10⁴ SVs per cultivar were identified, and an SV set was constructed. SVs affected the expression of a large number of genes during fiber elongation, and a transposable element insertion resulted in the glandless phenotype in upland cotton. Six widespread inversions were identified based on nine draft genomes and high-throughput chromosome conformation capture data. Multiple haplotype blocks that were always associated with aggregated SVs were demonstrated to play a pivotal role in the agronomic traits of upland cotton and drove its adaptation to the northern planting region. Exotic introgression was the source of these haplotype blocks and increased the genetic diversity of upland cotton. Our results enrich the genome resources of upland cotton, and the identified SVs will promote genetic and breeding research in cotton.

Development of the engineered "glanded plant and glandless seed" cotton

After fiber, cottonseed is the second most important by-product of cotton production. However, high concentrations of toxic free gossypol deposited in the glands of the cottonseed greatly hamper its effective usage as food or feed. Here, we developed a cotton line with edible cottonseed by specifically silencing the endogenous expression of GoPGF in the seeds, which led to a glandless phenotype with an ultra-low gossypol content in the seeds and nearly normal gossypol in other parts of the plants. This engineered cotton maintains normal resistance to insect pests, but the gossypol content in the seeds dropped by 98%, and thus, it can be consumed directly as food. The trait of a low gossypol content in the cottonseeds was stable and heritable, while the protein, oil content, and fiber yield or quality were nearly unchanged compared to the transgenic receptor W0. In addition, comparative transcriptome analysis showed that down-regulated genes in the ovules of the glandless cotton were enriched in terpenoid biosynthesis, indicating the underlying relationship between gland formation and gossypol biosynthesis. These results pave the way for the comprehensive utilization of cotton as a fiber, oil, and feed crop in the future.

Gossypol and related compounds are produced and accumulate in the aboveground parts of the cotton plant, independent of roots as the source

Use of Ultra-low gossypol cottonseed event as a scion in a graft combination confirmed that roots are not a source of terpenoids in the aboveground parts of a cotton plant. Gossypol and related terpenoids, derived from the same basic biosynthetic pathway, are present in the numerous lysigenous glands in the aboveground parts of a cotton plant. Roots, with sparse presence of such glands, do produce significant amount of gossypol and a different set of terpenoids. These compounds serve a defensive function against various pests and pathogens. This investigation was undertaken to examine whether gossypol produced in the roots can replenish the gossypol content of the cottonseed-glands that are largely devoid of this terpenoid in a genetically engineered event. Graft unions between a scion derived from the RNAi-based, Ultra-low gossypol cottonseed (ULGCS) event, TAM66274, and a rootstock derived from wild-type parental genotype, Coker 312 (Coker), were compared with various other grafts that served as controls. The results showed that the seeds developing within the scion of test grafts (ULGCS/Coker) continued to maintain the ultra-low gossypol levels found in the TAM66274 seeds. Molecular analyses confirmed that while the key gene involved in gland development showed normal activity in the developing embryos in the scion, two genes encoding the enzymes involved in gossypol biosynthesis were suppressed. Thus, the gene expression data confirmed the results obtained from biochemical measurements and collectively demonstrated that roots are not a source of gossypol for the aboveground parts of the cotton plant. These findings, combined with the results from previous investigations, support the assertion that gossypol and related terpenoids are produced in a highly localized manner in various organs of the cotton plant and are retained therein.

Carpel-specific down-regulation of GhCKXs in cotton significantly enhances seed and fiber yield

Cytokinin is considered to be an important driver of seed yield. To increase the yield of cotton while avoiding the negative consequences caused by constitutive overproduction of cytokinin, we down-regulated specifically the carpel genes for cytokinin oxidase/dehydrogenase (CKX), a key negative regulator of cytokinin levels, in transgenic cotton. The carpel-specific down-regulation of CKXs significantly enhanced cytokinin levels in the carpels. The elevated cytokinin promoted the expression of carpel- and ovule-development-associated genes, GhSTK2, GhAG1, and GhSHP, boosting ovule formation and thus producing more seeds in the ovary. Field experiments showed that the carpel-specific increase of cytokinin significantly increased both seed yield and fiber yield of cotton, without resulting in detrimental phenotypes. Our study details the regulatory mechanism of cytokinin signaling for seed development, and provides an effective and feasible strategy for yield improvement of seed crops.

Natural variation of GhSI7 increases seed index in cotton

qSI07.1, a major QTL for seed index in cotton, was fine-mapped to a 17.45-kb region, and the candidate gene GhSI7 was verified in transgenic plants. Improving production to meet human needs is a vital objective in cotton breeding. The yield-related trait seed index is a complex quantitative trait, but few candidate genes for seed index have been characterized. Here, a major QTL for seed index qSI07.1 was fine-mapped to a 17.45-kb region by linkage analysis and substitutional mapping. Only GhSI7, encoding the transcriptional regulator STERILE APETALA, was contained in the candidate region. Association test and genetic analysis indicated that an 845-bp-deletion in its intron was responsible for the seed index variation. Origin analysis revealed that this variation was unique in Gossypium hirsutum and originated from race accessions. Overexpression of GhSI7 (haplotype 2) significantly increased the seed index and organ size in cotton plants. Our findings provided a diagnostic marker for breeding and selecting cotton varieties with high seed index, and laid a foundation for further studies to understand the molecular mechanism of cotton seed morphogenesis.

Comparative transcriptome analysis reveals key pathways and genes involved in trichome development in tea plant (Camellia sinensis)

Trichomes, which develop from epidermal cells, are considered one of the important characteristics of the tea plant [Camellia sinensis (L.) O. Kuntze]. Many nutritional and metabolomic studies have indicated the important contributions of trichomes to tea products quality. However, understanding the regulation of trichome formation at the molecular level remains elusive in tea plants. Herein, we present a genome-wide comparative transcriptome analysis between the hairless Chuyeqi (CYQ) with fewer trichomes and the hairy Budiaomao (BDM) with more trichomes tea plant genotypes, toward the identification of biological processes and functional gene activities that occur during trichome development. In the present study, trichomes in both cultivars CYQ and BDM were unicellular, unbranched, straight, and soft-structured. The density of trichomes was the highest in the bud and tender leaf periods. Further, using the high-throughput sequencing method, we identified 48,856 unigenes, of which 31,574 were differentially expressed. In an analysis of 208 differentially expressed genes (DEGs) encoding transcription factors (TFs), five may involve in trichome development. In addition, on the basis of the Gene Ontology (GO) annotation and the weighted gene co-expression network analysis (WGCNA) results, we screened several DEGs that may contribute to trichome growth, including 66 DEGs related to plant resistance genes (PRGs), 172 DEGs related to cell wall biosynthesis pathway, 29 DEGs related to cell cycle pathway, and 45 DEGs related to cytoskeleton biosynthesis. Collectively, this study provided high-quality RNA-seq information to improve our understanding of the molecular regulatory mechanism of trichome development and lay a foundation for additional trichome studies in tea plants.

Visible gland constantly traces virus-induced gene silencing in cotton

A virus-induced gene silencing (VIGS) system was established to induce endogenous target gene silencing by post-transcriptional gene silencing (PTGS), which is a powerful tool for gene function analysis in plants. Compared with stable transgenic plant via Agrobacterium-mediated gene transformation, phenotypes after gene knockdown can be obtained rapidly, effectively, and high-throughput through VIGS system. This approach has been successfully applied to explore unknown gene functions involved in plant growth and development, physiological metabolism, and biotic and abiotic stresses in various plants. In this system, GhCLA1 was used as a general control, however, silencing of this gene leads to leaf albino, wilting, and plant death ultimately. As such, it cannot indicate the efficiency of target gene silencing throughout the whole plant growth period. To address this question, in this study, we developed a novel marker gene, Gossypium PIGMENT GLAND FORMATION GENE (GoPGF), as the control to trace the efficiency of gene silencing in the infected tissues. GoPGF has been proved a key gene in gland forming. Suppression of GoPGF does not affect the normal growth and development of cotton. The number of gland altered related to the expression level of GoPGF gene. So it is a good marker that be used to trace the whole growth stages of plant. Moreover, we further developed a method of friction inoculation to enhance and extend the efficiency of VIGS, which facilitates the analysis of gene function in both the vegetative stage and reproductive stage. This improved VIGS technology will be a powerful tool for the rapid functional identification of unknown genes in genomes.

GhMYC2 activates cytochrome P450 gene CYP71BE79 to regulate gossypol biosynthesis in cotton

GhMYC2 regulates the gossypol biosynthesis pathway in cotton through activation of the expression of gossypol synthesis gene CYP71BE79, CDNC, CYP706B1, DH1, and CYP82D113. Cotton is one of the main cash crops globally. Cottonseed contains fiber, fat, protein, and starch, and has important economic value. However, gossypol in cottonseed seriously affects the development and utilization of cottonseed. Nonetheless, gossypol has great application potential in agriculture, medicine, and industry. Therefore, it is very important to study gossypol biosynthesis and its upstream regulatory pathways. It has been reported that the content of gossypol in hairy roots of cotton is regulated through jasmonic acid signaling; however, the specific molecular mechanism has not been revealed yet. We found that the expression of basic helix-loop-helix family transcription factor GhMYC2 was significantly upregulated after exogenous administration of methyl jasmonate to cotton seedlings, and the content of gossypol changed significantly with the variation of GhMYC2 expression. Further studies revealed that GhMYC2 could specifically bind to the G-Box in the promoter region of CDNC, CYP706B1, DH1, CYP82D113, CYP71BE79 to activate its expression and regulate gossypol synthesis, and its activation of CYP71BE79 promoter was inhibited by GhJAZ2. Not only that GhMYC2 could also interact with GoPGF. In this work, the molecular mechanisms of gossypol biosynthesis regulated by GhMYC2 were analyzed. The results provide a theoretical basis for cultivating new varieties of low-gossypol or high-gossypol cotton and creating excellent germplasm resources.

Combined morphological and multi-omics analyses to reveal the developmental mechanism of Zanthoxylum bungeanum prickles

Zanthoxylum bungeanum Maxim. as an important economic forest, its epidermis bears prickles which complicate the harvesting process and increase the labor costs. To explore the developmental mechanism of prickles, three varieties of Zanthoxylum bungeanum (PZB, SZB, GSZB) were selected for morphological and multi-omics analyses. The absorption spectra of prickles and stems were detected using Fourier-transform infrared spectroscopy (FTIR), and they were found different at 1617, 1110, 3319, and 1999 cm^-1. The morphology of prickles and stems were observed using light microscopy and transmission electron microscopy (TEM). The growth direction of cells on the prickle side and stem side were perpendicular to each other, and there was a resembling abscission zone (RAZ) between them. The vacuolar deposits of prickle cells were much more than stem cells, indicating that the lignification degree of prickles was higher than stems. In addition, 9 candidate genes (ZbYABBY2, ZbYABBY1, ZbYABBY5, ZbWRKY, ZbLOG5, ZbAZG2, ZbGh16, ZbIAA33, and ZbGh16X1) were screened out and validated base on transcriptome and qRT-PCA. As well as, 30 key metabolites were found related to prickle development base on metabolome analysis. Among them, 4-hydroxy-2-oxopentanoate, trans-2-hydroxy-cinnamate, trans-cinnamate, polyhydroxy-fatty acid, 10,16-dihydroxypalmitate, cinnamic acid were related to the biosynthesis of cutin, suberine and wax. Indole-3-acetate, tryptamine, anthranilate, fromylanthranilate, N6-(delta2-isopentenyl)-adenine were related to plant hormone signal transduction. Generally, this is the first study to reveal the developmental mechanism of prickles. The results of this study lay the foundation for the breeding of non-prickle Zanthoxylum bungeanum.

Transcriptomic and functional analysis provides molecular insights into multicellular trichome development

Trichomes, the hair-like structures located on aerial parts of most vascular plants, are associated with a wide array of biological processes and affect the economic value of certain species. The processes involved in unicellular trichome formation have been well-studied in Arabidopsis (Arabidopsis thaliana). However, our understanding of the morphological changes and the underlying molecular processes involved in multicellular trichome development is limited. Here, we studied the dynamic developmental processes involved in glandular and nonglandular multicellular trichome formation in cucumber (Cucumis sativus L.) and divided these processes into five sequential stages. To gain insights into the underlying mechanisms of multicellular trichome formation, we performed a time-course transcriptome analysis using RNA-sequencing analysis. A total of 711 multicellular trichome-related genes were screened and a model for multicellular trichome formation was developed. The transcriptome and co-expression datasets were validated by reverse transcription-quantitative PCR and in situ hybridization. In addition, virus-induced gene silencing analysis revealed that CsHOMEOBOX3 (CsHOX3) and CsbHLH1 are involved in nonglandular trichome elongation and glandular trichome formation, respectively, which corresponds with the transcriptome data. This study presents a transcriptome atlas that provides insights into the molecular processes involved in multicellular trichome formation in cucumber and can be an important resource for future functional studies.

An Overview of Cotton Gland Development and Its Transcriptional Regulation

Cotton refers to species in the genus Gossypium that bear spinnable seed coat fibers. A total of 50 species in the genus Gossypium have been described to date. Of these, only four species, viz. Gossypium, hirsutum, G. barbadense, G. arboretum, and G. herbaceum are cultivated; the rest are wild. The black dot-like structures on the surfaces of cotton organs or tissues, such as the leaves, stem, calyx, bracts, and boll surface, are called gossypol glands or pigment glands, which store terpenoid aldehydes, including gossypol. The cotton (Gossypium hirsutum) pigment gland is a distinctive structure that stores gossypol and its derivatives. It provides an ideal system for studying cell differentiation and organogenesis. However, only a few genes involved in the process of gland formation have been identified to date, and the molecular mechanisms underlying gland initiation remain unclear. The terpenoid aldehydes in the lysigenous glands of Gossypium species are important secondary phytoalexins (with gossypol being the most important) and one of the main defenses of plants against pests and diseases. Here, we review recent research on the development of gossypol glands in Gossypium species, the regulation of the terpenoid aldehyde biosynthesis pathway, discoveries from genetic engineering studies, and future research directions.

Identification and stress function verification of the HAK/KUP/KT family in Gossypium hirsutum

The potassium transporter family HAK/KUP/KT is a large group of proteins that are important in plant potassium transport and play a crucial role in plant growth and development. The members of the family play an important role in the response of plants to abiotic stress by maintaining osmotic balance. However, the function of the family in cotton is unclear. In this study, whole genome identification and characterization of the HAK/KUP/KT family from upland cotton (Gossypium hirsutum) were carried out. Bioinformatics methods were used to identify HAK/KUP/KT family members from the G. hirsutum genome and to analyse the physical and chemical properties, basic characteristics, phylogeny, chromosome location and expression of HAK/KUP/KT family members. A total of 41 HAK/KUP/KT family members were identified in the G. hirsutum genome. Phylogenetic analysis grouped these genes into four clusters (I, II, III, IV), containing 6, 10, 3 and 22 genes, respectively. Chromosomal distribution, gene structure and conserved motif analyses of the 41 GhHAK genes were subsequently performed. The RNA-seq data and qRT-PCR results showed that the family had a wide range of tissue expression patterns, and they responded to certain drought stresses. Through expression analysis, seven HAK/KUP/KT genes involved in drought stress were screened, and four genes with obvious phenotypes under drought stress were obtained by VIGS verification, which laid a theoretical foundation for the function of the cotton HAK/KUP/KT family.

GhCDNC and GhCYP706B1 genes mediate gossypol biosynthesis in upland cotton

BACKGROUND

In cotton the identification and characterization of natural defense is a cost-effective, sustainable, and environment-friendly strategy to combat cotton pests. The secondary metabolites traits in cotton plant i.e., toxic gossypol glands play significant role for development and self-defense mechanism. To utilize gossypol in breeding implements, the understanding of gossypol initiation biosynthesis genes has vital importance at reproductive organ development stages.

METHODS

Cotton germplasm of 100 genotypes screened visually based on gossypol glandedness and a core set of ten genotypes was developed. Further three genotypes FH-330 (high glanding), F-280 (low glanding) and IRMA-197 (glandless) were used for determining the transcript abundance of twelve gossypol biosynthesis genes.

RESULTS

Out of 100, germplasm categorized as (76) high glanding, (22) medium glanding), one genotype for each (low glanding) and (glandless) category. Real-time qPCR analysis revealed varied expression patterns among selected three genotypes. Out of twelve, three genes CYP706B1, CDNC and 2ODD-1 had strong expression levels in all tested tissues in high glanded genotype, while, slight or no expression of these genes was recorded in low glanding and glandless genotype, respectively. The shell of developing boll (10, 20, 30 DPA), and developing embryo (20, 30 DPA) showed substantially medium to maximum expression, respectively while high to medium expression was recorded in sepals and leaf tissue.

CONCLUSIONS

Our study demonstrated that CYP706B1, CDNC and 2ODD-1 are the most promising genes involved in gossypol biosynthesis. Developing boll shell, developing embryo, leaf and sepal also have significant ability to synthesize gossypol. This will provide scientists a way to manipulate gossypol contents in economically important organs of cotton plant for targeted breeding.

Transcriptional Landscape of Cotton Fiber Development and Its Alliance With Fiber-Associated Traits

Cotton fiber development is still an intriguing question to understand fiber commitment and development. At different fiber developmental stages, many genes change their expression pattern and have a pivotal role in fiber quality and yield. Recently, numerous studies have been conducted for transcriptional regulation of fiber, and raw data were deposited to the public repository for comprehensive integrative analysis. Here, we remapped > 380 cotton RNAseq data with uniform mapping strategies that span ∼400 fold coverage to the genome. We identified stage-specific features related to fiber cell commitment, initiation, elongation, and Secondary Cell Wall (SCW) synthesis and their putative cis-regulatory elements for the specific regulation in fiber development. We also mined Exclusively Expressed Transcripts (EETs) that were positively selected during cotton fiber evolution and domestication. Furthermore, the expression of EETs was validated in 100 cotton genotypes through the nCounter assay and correlated with different fiber-related traits. Thus, our data mining study reveals several important features related to cotton fiber development and improvement, which were consolidated in the "CottonExpress-omics" database.

Comparative transcriptome analysis reveals the regulation network for fiber strength in cotton

OBJECTIVE

Determine the effect of secondary cell wall (SCW) thickness and microcrystalline cellulose content (MCC) on mature fiber strength (FS) and reveal through comparative transcriptome analysis the molecular regulation network governing FS in cotton.

RESULTS

Transmission electron microscope (TEM) analysis of two parent varieties, Prema with elite FS and 86-1 with weak fiber, revealed significant difference in the SCW but not in MCC. Transcriptome analysis revealed that genes differentially expressed during SCW thickening (20 DPA) are highly related to FS; in particular, up-regulated genes such as UDPG, CESA2, and NAC83 were important in SCW thickening, likely contributing to higher FS. GO and KEGG enrichment analysis revealed the common up-regulated genes to be enriched in carbon metabolism and terms relating to the cell wall.

CONCLUSIONS

We developed two recombinant inbred lines with elite FS, selected from the filial generation of Prema and 86-1. By comparing transcriptomic data, we revealed the gene expression network governing SCW thickness in mature fiber. Our results provide solid insights into the relationship of the SCW and FS.

Cottonseed extracts regulate gene expression in human colon cancer cells

Cotton plant provides economically important fiber and cottonseed, but cottonseed contributes 20% of the crop value. Cottonseed value could be increased by providing high value bioactive compounds and polyphenolic extracts aimed at improving nutrition and preventing diseases because plant polyphenol extracts have been used as medicinal remedy for various diseases. The objective of this study was to investigate the effects of cottonseed extracts on cell viability and gene expression in human colon cancer cells. COLO 225 cells were treated with ethanol extracts from glanded and glandless cottonseed followed by MTT and qPCR assays. Cottonseed extracts showed minor effects on cell viability. qPCR assay analyzed 55 mRNAs involved in several pathways including DGAT, GLUT, TTP, IL, gossypol-regulated and TTP-mediated pathways. Using BCL2 mRNA as the internal reference, qPCR analysis showed minor effects of ethanol extracts from glanded seed coat and kernel and glandless seed coat on mRNA levels in the cells. However, glandless seed kernel extract significantly reduced mRNA levels of many genes involved in glucose transport, lipid biosynthesis and inflammation. The inhibitory effects of glandless kernel extract on gene expression may provide a useful opportunity for improving nutrition and healthcare associated with colon cancer. This in turn may provide the potential of increasing cottonseed value by using ethanol extract as a nutrition/health intervention agent.

Optimizing the Protein Fluorescence Reporting System for Somatic Embryogenesis Regeneration Screening and Visual Labeling of Functional Genes in Cotton

Protein fluorescence reporting systems are of crucial importance to in-depth life science research, providing systematic labeling tools for visualization of microscopic biological activities in vivo and revolutionizing basic research. Cotton somatic cell regeneration efficiency is low, causing difficulty in cotton transformation. It is conducive to screening transgenic somatic embryo using the fluorescence reporting system. However, available fluorescence labeling systems in cotton are currently limited. To optimize the fluorescence reporting system of cotton with an expanded range of available fluorescent proteins, we selected 11 fluorescent proteins covering red, green, yellow, and cyan fluorescence colors and expressed them in cotton. Besides mRuby2 and G3GFP, the other nine fluorescent proteins (mCherry, tdTomato, sfGFP, Clover, EYFP, YPet, mVenus, mCerulean, and ECFP) were stably and intensely expressed in transgenic callus and embryo, and inherited in different cotton organs derive from the screened embryo. In addition, transgenic cotton expressing tdTomato appears pink under white light, not only for callus and embryo tissues but also various organs of mature plants, providing a visual marker in the cotton genetic transformation process, accelerating the evaluation of transgenic events. Further, we constructed transgenic cotton expressing mCherry-labeled organelle markers in vivo that cover seven specific subcellular compartments: plasma membrane, endoplasmic reticulum, tonoplast, mitochondrion, plastid, Golgi apparatus, and peroxisome. We also provide a simple and highly efficient strategy to quickly determine the subcellular localization of uncharacterized proteins in cotton cells using organelle markers. Lastly, we built the first cotton stomatal fluorescence reporting system using stomata-specific expression promoters (ProKST1, ProGbSLSP, and ProGC1) to drive Clover expression. The optimized fluorescence labeling system for transgenic somatic embryo screening and functional gene labeling in this study offers the potential to accelerating somatic cell regeneration efficiency and the in vivo monitoring of diverse cellular processes in cotton.

Genome-Wide Analysis of the Cytochrome P450 Gene Family Involved in Salt Tolerance in Gossypium hirsutum

Plant cytochrome P450 (P450) participates in a wide range of biosynthetic reactions and targets a variety of biological molecules. These reactions lead to various fatty acid conjugates, plant hormones, secondary metabolites, lignin, and various defensive compounds. In our previous research, transcriptome analysis was performed on the salt-tolerant upland cotton "Tongyan No. 1." Many differentially expressed genes (DEGs) belong to the P450 family, and their domains occur widely in plants. In this current research, P450 genes were identified in Gossypium hirsutum with the aid of bioinformatics methods for investigating phylogenetic relations, gene structure, cis-elements, chromosomal localization, and collinearity within a genome. qRT-PCR was conducted to analyze P450 gene expression patterns under salt stress. The molecular weights of the 156 P450 genes were in the range of 5,949.6-245,576.3 Da, and the length of the encoded amino acids for all the identified P450 genes ranged from 51 to 2,144. P450 proteins are divided into four different subfamilies based on phylogenetic relationship, gene structure, and chromosomal localization of gene replication. The length of P450 genes in upland cotton differs greatly, ranging from 1,500 to 13,000 bp. The number of exons in the P450 family genes ranged from 1 to 9, while the number of introns ranged from 0 to 8, and there were similar trends within clusters. A total of 31 cis-acting elements were identified by analyzing 1,500 bp promoter sequences. Differences were found in cis-acting elements among genes. The consistency between qRT-PCR and previous transcriptome analysis of salt tolerance DEGs indicated that they were likely to be involved in the salt tolerance of cotton seedlings. Our results provide valuable information on the evolutionary relationships of genes and functional characteristics of the gene family, which is beneficial for further study of the cotton P450 gene family.

Identification and characteristics of a novel gland-forming gene in cotton

The cotton (Gossypium hirsutum) pigment gland is a distinctive structure that functions as the main deposit organ of gossypol and its derivatives. It is also an ideal system in which to study cell differentiation and organogenesis. However, only a few genes that determine the process of gland formation have been reported, including GoPGF, CGP1, and CGFs; the molecular mechanisms underlying gland initiation are still largely unclear. Here, we report the discovery of the novel stem pigment gland-forming gene GoSPGF by map-based cloning; annotated as a GRAS transcription factor, this gene is responsible for the glandless trait specifically on the stem. In the stem glandless mutant T582, a point mutation (C to A) was found to create a premature stop codon and truncate the protein. Similarly, virus-induced gene silencing of GoSPGF resulted in glandless stems and dramatically reduced gossypol content. Comparative transcriptomic data showed that loss of GoSPGF significantly suppressed expression of many genes involved in gossypol biosynthesis and altered expression of genes involved in gibberellic acid signaling/biosynthesis. Overall, these findings provide more insight into the networks regulating glandular structure differentiation and formation in cotton, which will be helpful for understanding other plants bearing special gland structures such as tobacco (Nicotiana benthamiana), artemisia annua, mint (Mentha spp.), and rubber (Hevea brasiliensis).

Genetics and Breeding for Glandless Upland Cotton With Improved Yield Potential and Disease Resistance: A Review

Glandless cotton (devoid of toxic gossypol) can be grown as a triple-purpose crop for fiber, feeds, and food (as an oil and protein source). However, its sensitivity to insect pests and its low yield due to the lack of breeding activities has prevented the realization of its potential in commercial seed production and utilization. Since the mid-1990s, the commercialization of bollworm and budworm resistant Bt cotton and the eradication of boll weevils and pink bollworms have provided an opportunity to revitalize glandless cotton production in the United States. The objectives of this study were to review the current status of genetics and breeding for glandless cotton, with a focus on the progress in breeding for glandless Upland cotton in New Mexico, United States. Because there existed a 10-20% yield gap between the best existing glandless germplasm and commercial Upland cultivars, the breeding of glandless Upland cultivars with improved yield and disease resistance was initiated at the New Mexico State University more than a decade ago. As a result, three glandless Upland cultivars, i.e., long-staple Acala 1517-18 GLS, medium staple NuMex COT 15 GLS, and NuMex COT 17 GLS with Fusarium wilt race 4 resistance were released. However, to compete with the current commercial glanded cotton, more breeding efforts are urgently needed to introduce different glandless traits (natural mutations, transgenic or genome-editing) into elite cotton backgrounds with high yields and desirable fiber quality.