Genes regulating gland development in the cotton plant

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.

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.

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.

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.

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.

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.

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.

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.

Highly efficient Agrobacterium rhizogenes-mediated hairy root transformation for gene editing analysis in cotton

CRIPSR/Cas9 gene editing system is an effective tool for genome modification in plants. Multiple target sites are usually designed and the effective target sites are selected for editing. Upland cotton (Gossypium hirsutum L., hereafter cotton) is allotetraploid and is commonly considered as difficult and inefficient to transform, it is important to select the effective target sites that could result in the ideal transgenic plants with the CRISPR-induced mutations. In this study, Agrobacterium rhizogenes-mediated hairy root method was optimized to detect the feasibility of the target sites designed in cotton phytoene desaturase (GhPDS) gene. A. rhizogenes showed the highest hairy root induction (30%) when the bacteria were cultured until OD₆₀₀ reached to 0.8. This procedure was successfully applied to induce hairy roots in the other three cultivars (TM-1, Lumian-21, Zhongmian-49) and the mutations were detected in GhPDS induced by CRISPR/Cas9 system. Different degrees of base deletions at two sgRNAs (sgRNA5 and sgRNA10) designed in GhPDS were detected in R15 hairy roots. Furthermore, we obtained an albino transgenic cotton seeding containing CRISPR/Cas9-induced gene editing mutations in sgRNA10. The hairy root transformation system established in this study is sufficient for selecting sgRNAs in cotton, providing a technical basis for functional genomics research of cotton.

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.

USDA's revised biotechnology regulation's contribution to increasing agricultural sustainability and responding to climate change

Biotechnology can provide a valuable tool to meet UN Sustainable Development Goals and U.S. initiatives to find climate solutions and improve agricultural sustainability. The literature contains hundreds of examples of crops that may serve this purpose, yet most remain un-launched due to high regulatory barriers. Recently the USDA revised its biotechnology regulations to make them more risk-proportionate, science-based, and streamlined. Here, we review some of the promising leads that may enable agriculture to contribute to UN sustainability goals. We further describe and discuss how the revised biotechnology regulation would hypothetically apply to these cases.

Cottonseed Meal Protein Isolate as a New Source of Alternative Proteins: A Proteomics Perspective

Cottonseed meal (CSM) is a good source of dietary proteins but is unsuitable for human consumption due to its gossypol content. To unlock its potential, we developed a protein extraction process with a gossypol removal treatment to generate CSM protein isolate (CSMPI) with ultra-low gossypol content. This process successfully reduced the free and total gossypol content to 4.8 ppm and 147.2 ppm, respectively, far below the US FDA limit. In addition, the functional characterisation of CSMPI revealed a better oil absorption capacity and water solubility than pea protein isolate. Proteome profiling showed that the treatment improved protein identification, while SDS-PAGE analysis indicated that the treatment did not induce protein degradation. Amino acid analysis revealed that post-treated CSMPI was rich in branched-chain amino acids (BCAAs). Mass spectrometry analysis of various protein fractions obtained from an in vitro digestibility assay helped to establish the digestibility profile of CSM proteins. Several potential allergens in CSMPI were also found using allergenic prediction software, but further evaluation based on their digestibility profiles and literature reviews suggests that the likelihood of CSMPI allergenicity remains low. Overall, our results help to navigate and direct the application of CSMPIs as alternative proteins toward nutritive human food application.

Comparative Transcriptome Analysis Reveals Genes Associated with the Gossypol Synthesis and Gland Morphogenesis in Gossypium hirsutum

Gossypium hirsutum is an important source of natural textile fibers. Gossypol, which is a sesquiterpenoid compound mainly existing in the cotton pigment glands, can facilitate resistance to the stress from diseases and pests. The level of gossypol in the cotton is positively correlated to the quantity of pigment glands. However, the underlying regulatory mechanisms of gossypol synthesis and gland morphogenesis are still poorly understood, especially from a transcriptional perspective. The transcripts of young leaves and ovules at 30 DPA of the glanded plants and glandless plants were studied by RNA-Seq and 865 million clean reads were obtained. A total of 34,426 differentially expressed genes (DEGs) were identified through comparative transcriptome analysis. Genes related to gossypol synthesis or gland morphogenesis displayed significant differential expression between the two cultivars. Functional annotation revealed that the candidate genes related to catalytic activity, the biosynthesis of secondary metabolites, and biomolecular decomposition processes. Our work herein unveiled several potential candidate genes related to gossypol synthesis or gland morphogenesis and may provide useful clues for a breeding program of cotton cultivars with low cottonseed gossypol contents.

Identification and characterization of a novel gene involved in glandular trichome development in Nepeta tenuifolia

Nepeta tenuifolia is a medicinal plant rich in terpenoids and flavonoids with antiviral, immunoregulatory, and anti-inflammatory activities. The peltate glandular trichome (PGT) is a multicellular structure considered to be the primary storage organ for monoterpenes; it may serve as an ideal model for studying cell differentiation and the development of glandular trichomes (GTs). The genes that regulate the development of GTs have not yet been well studied. In this study, we identified NtMIXTA1, a GT development-associated gene from the R2R3 MYB SBG9 family. NtMIXTA1 overexpression in tobacco resulted in the production of longer and denser GTs. Virus-induced gene silencing of NtMIXTA1 resulted in lower PGT density, a significant reduction in monoterpene concentration, and the decreased expression of genes related to monoterpene biosynthesis. Comparative transcriptome and widely targeted metabolic analyses revealed that silencing NtMIXTA1 significantly influenced the expression of genes, and the production of metabolites involved in the biosynthesis of terpenoids, flavonoids, and lipids. This study provides a solid foundation describing a mechanism underlying the regulation of GT development. In addition, this study further deepens our understanding of the regulatory networks involved in GT development and GT development-associated metabolite flux, as well as provides valuable reference data for studying plants with a high medicinal value without genetic transformation.

VIGS Goes Viral: How VIGS Transforms Our Understanding of Plant Science

Virus-induced gene silencing (VIGS) has developed into an indispensable approach to gene function analysis in a wide array of species, many of which are not amenable to stable genetic transformation. VIGS utilizes the posttranscriptional gene silencing (PTGS) machinery of plants to restrain viral infections systemically and is used to downregulate the plant's endogenous genes. Here, we review the molecular mechanisms of DNA- and RNA-virus-based VIGS, its inherent connection to PTGS, and what is known about the systemic spread of silencing. Recently, VIGS-based technologies have been expanded to enable not only gene silencing but also overexpression [virus-induced overexpression (VOX)], genome editing [virus-induced genome editing (VIGE)], and host-induced gene silencing (HIGS). These techniques expand the genetic toolbox for nonmodel organisms even more. Further, we illustrate the versatility of VIGS and the methods derived from it in elucidating molecular mechanisms, using tomato fruit ripening and programmed cell death as examples. Finally, we discuss challenges of and future perspectives on the use of VIGS to advance gene function analysis in nonmodel plants in the postgenomic era.

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.

Leveraging National Germplasm Collections to Determine Significantly Associated Categorical Traits in Crops: Upland and Pima Cotton as a Case Study

Observable qualitative traits are relatively stable across environments and are commonly used to evaluate crop genetic diversity. Recently, molecular markers have largely superseded describing phenotypes in diversity surveys. However, qualitative descriptors are useful in cataloging germplasm collections and for describing new germplasm in patents, publications, and/or the Plant Variety Protection (PVP) system. This research focused on the comparative analysis of standardized cotton traits as represented within the National Cotton Germplasm Collection (NCGC). The cotton traits are named by 'descriptors' that have non-numerical sub-categories (descriptor states) reflecting the details of how each trait manifests or is absent in the plant. We statistically assessed selected accessions from three major groups of Gossypium as defined by the NCGC curator: (1) "Stoneville accessions (SA)," containing mainly Upland cotton (Gossypium hirsutum) cultivars; (2) "Texas accessions (TEX)," containing mainly G. hirsutum landraces; and (3) Gossypium barbadense (Gb), containing cultivars or landraces of Pima cotton (Gossypium barbadense). For 33 cotton descriptors we: (a) revealed distributions of character states for each descriptor within each group; (b) analyzed bivariate associations between paired descriptors; and (c) clustered accessions based on their descriptors. The fewest significant associations between descriptors occurred in the SA dataset, likely reflecting extensive breeding for cultivar development. In contrast, the TEX and Gb datasets showed a higher number of significant associations between descriptors, likely correlating with less impact from breeding efforts. Three significant bivariate associations were identified for all three groups, bract nectaries:boll nectaries, leaf hair:stem hair, and lint color:seed fuzz color. Unsupervised clustering analysis recapitulated the species labels for about 97% of the accessions. Unexpected clustering results indicated accessions that may benefit from potential further investigation. In the future, the significant associations between standardized descriptors can be used by curators to determine whether new exotic/unusual accessions most closely resemble Upland or Pima cotton. In addition, the study shows how existing descriptors for large germplasm datasets can be useful to inform downstream goals in breeding and research, such as identifying rare individuals with specific trait combinations and targeting breakdown of remaining trait associations through breeding, thus demonstrating the utility of the analytical methods employed in categorizing germplasm diversity within the collection.

CRISPR/Cas9-Mediated Mutagenesis of the Granule-Bound Starch Synthase Gene in the Potato Variety Yukon Gold to Obtain Amylose-Free Starch in Tubers

Potato (Solanum tuberosum L.) is the third most important food crop after rice and wheat. Its tubers are a rich source of dietary carbohydrates in the form of starch, which has many industrial applications. Starch is composed of two polysaccharides, amylose and amylopectin, and their ratios determine different properties and functionalities. Potato varieties with higher amylopectin have many food processing and industrial applications. Using Agrobacterium-mediated transformation, we delivered Clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9) reagents to potato (variety Yukon Gold) cells to disrupt the granule-bound starch synthase (gbssI) gene with the aim of eliminating the amylose component of starch. Lugol-Iodine staining of the tubers showed a reduction or complete elimination of amylose in some of the edited events. These results were further confirmed by the perchloric acid and enzymatic methods. One event (T2-7) showed mutations in all four gbss alleles and total elimination of amylose from the tubers. Viscosity profiles of the tuber starch from six different knockout events were determined using a Rapid Visco Analyzer (RVA), and the values reflected the amylopectin/amylose ratio. Follow-up studies will focus on eliminating the CRISPR components from the events and on evaluating the potential of clones with various amylose/amylopectin ratios for food processing and other industrial applications.

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.

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.

Glandular trichomes: new focus on horticultural crops

Plant glandular trichomes (GTs) are epidermal outgrowths with the capacity to biosynthesize and secrete specialized metabolites, that are of great scientific and practical significance. Our understanding of the developmental process of GTs is limited, and no single plant species serves as a unique model. Here, we review the genetic mechanisms of GT initiation and development and provide a summary of the biosynthetic pathways of GT-specialized metabolites in nonmodel plant species, especially horticultural crops. We discuss the morphology and classification of GT types. Moreover, we highlight technological advancements in methods employed for investigating GTs. Understanding the molecular basis of GT development and specialized metabolites not only offers useful avenues for research in plant breeding that will lead to the improved production of desirable metabolites, but also provides insights for plant epidermal development research.

Recent Advances and Future Perspectives in Cotton Research

Cotton is not only the world's most important natural fiber crop, but it is also an ideal system in which to study genome evolution, polyploidization, and cell elongation. With the assembly of five different cotton genomes, a cotton-specific whole-genome duplication with an allopolyploidization process that combined the A- and D-genomes became evident. All existing A-genomes seemed to originate from the A₀-genome as a common ancestor, and several transposable element bursts contributed to A-genome size expansion and speciation. The ethylene production pathway is shown to regulate fiber elongation. A tip-biased diffuse growth mode and several regulatory mechanisms, including plant hormones, transcription factors, and epigenetic modifications, are involved in fiber development. Finally, we describe the involvement of the gossypol biosynthetic pathway in the manipulation of herbivorous insects, the role of GoPGF in gland formation, and host-induced gene silencing for pest and disease control. These new genes, modules, and pathways will accelerate the genetic improvement of 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.

Molecular cloning and characterization of GhERF105, a gene contributing to the regulation of gland formation in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Gossypium hirsutum L. (cotton) is one of the most economically important crops in the world due to its significant source of fiber, feed, foodstuff, oil and biofuel products. However, the utilization of cottonseed was limited due to the presence of small and darkly pigmented glands that contain large amounts of gossypol, which is toxic to human beings and non-ruminant animals. To date, some progress has been made in the pigment gland formation, but the underlying molecular mechanism of its formation was still unclear.

RESULTS

In this study, we identified an AP2/ERF transcription factor named GhERF105 (GH_A12G2166), which was involved in the regulation of gland pigmentation by the comparative transcriptome analysis of the leaf of glanded and glandless plants. It encoded an ERF protein containing a converved AP2 domain which was localized in the nucleus with transcriptional activity, and showed the high expression in glanded cotton accessions that contained much gossypol. Virus-induced gene silencing (VIGS) against GhERF105 caused the dramatic reduction in the number of glands and significantly lowered levels of gossypol in cotton leaves. GhERF105 showed the patterns of spatiotemporal and inducible expression in the glanded plants.

CONCLUSIONS

These results suggest that GhERF105 contributes to the pigment gland formation and gossypol biosynthesis in partial organs of glanded plant. It also provides a potential molecular basis to generate 'glandless-seed' and 'glanded-plant' cotton cultivar.

The gland localized CGP1 controls gland pigmentation and gossypol accumulation in cotton

Pigment glands, also known as black glands or gossypol glands, are specific for Gossypium spp. These glands strictly confine large amounts of secondary metabolites to the lysigenous cavity, leading to the glands' intense colour and providing defence against pests and pathogens. This study performed a comparative transcriptome analysis of glanded versus glandless cotton cultivars. Twenty-two transcription factors showed expression patterns associated with pigment glands and were characterized. Phenotypic screening of the genes, via virus-induced gene silencing, showed an apparent disappearance of pigmented glands after the silencing of a pair of homologous MYB-encoding genes in the A and D genomes (designated as CGP1). Further study showed that CGP1a encodes an active transcription factor, which is specifically expressed in the gland structure, while CGP1d encodes a non-functional protein due to a fragment deletion, which causes premature termination. RNAi-mediated silencing and CRISPR knockout of CGP1 in glanded cotton cultivars generated a glandless-like phenotype, similar to the dominant glandless mutant Gl₂^e . Microscopic analysis showed that CGP1 knockout did not affect gland structure or density, but affected gland pigmentation. The levels of gossypol and related terpenoids were significantly decreased in cgp1 mutants, and a number of gossypol biosynthetic genes were strongly down-regulated. CGP1 is located in the nucleus where it interacts with GoPGF, a critical transcription factor for gland development and gossypol synthesis. Our data suggest that CGP1 and GoPGF form heterodimers to control the synthesis of gossypol and other secondary metabolites in cotton.

Genome sequencing of the Australian wild diploid species Gossypium australe highlights disease resistance and delayed gland morphogenesis

The diploid wild cotton species Gossypium australe possesses excellent traits including resistance to disease and delayed gland morphogenesis, and has been successfully used for distant breeding programmes to incorporate disease resistance traits into domesticated cotton. Here, we sequenced the G. australe genome by integrating PacBio, Illumina short read, BioNano (DLS) and Hi-C technologies, and acquired a high-quality reference genome with a contig N50 of 1.83 Mb and a scaffold N50 of 143.60 Mb. We found that 73.5% of the G. australe genome is composed of various repeat sequences, differing from those of G. arboreum (85.39%), G. hirsutum (69.86%) and G. barbadense (69.83%). The G. australe genome showed closer collinear relationships with the genome of G. arboreum than G. raimondii and has undergone less extensive genome reorganization than the G. arboreum genome. Selection signature and transcriptomics analyses implicated multiple genes in disease resistance responses, including GauCCD7 and GauCBP1, and experiments revealed induction of both genes by Verticillium dahliae and by the plant hormones strigolactone (GR24), salicylic acid (SA) and methyl jasmonate (MeJA). Experiments using a Verticillium-resistant domesticated G. barbadense cultivar confirmed that knockdown of the homologues of these genes caused a significant reduction in resistance against Verticillium dahliae. Moreover, knockdown of a newly identified gland-associated gene GauGRAS1 caused a glandless phenotype in partial tissues using G. australe. The G. australe genome represents a valuable resource for cotton research and distant relative breeding as well as for understanding the evolutionary history of crop genomes.