Metabolic engineering of gossypol in cotton

Probiotic Pediococcus pentosaceus restored gossypol-induced intestinal barrier injury by increasing propionate content in Nile tilapia

BACKGROUND

Intestinal barrier is a dynamic interface between the body and the ingested food components, however, dietary components or xenobiotics could compromise intestinal integrity, causing health risks to the host. Gossypol, a toxic component in cottonseed meal (CSM), caused intestinal injury in fish or other monogastric animals. It has been demonstrated that probiotics administration benefits the intestinal barrier integrity, but the efficacy of probiotics in maintaining intestinal health when the host is exposed to gossypol remains unclear. Here, a strain (YC) affiliated to Pediococcus pentosaceus was isolated from the gut of Nile tilapia (Oreochromis niloticus) and its potential to repair gossypol-induced intestinal damage was evaluated.

RESULTS

A total of 270 Nile tilapia (2.20 ± 0.02 g) were allotted in 3 groups with 3 tanks each and fed with 3 diets including CON (control diet), GOS (control diet containing 300 mg/kg gossypol) and GP (control diet containing 300 mg/kg gossypol and 108 colony-forming unit (CFU)/g P. pentosaceus YC), respectively. After 10 weeks, addition of P. pentosaceus YC restored growth retardation and intestinal injury induced by gossypol in Nile tilapia. Transcriptome analysis and siRNA interference experiments demonstrated that NOD-like receptors (NLR) family caspase recruitment domain (CARD) domain containing 3 (Nlrc3) inhibition might promote intestinal stem cell (ISC) proliferation, as well as maintaining gut barrier integrity. 16S rRNA sequencing and gas chromatography-mass spectrometry (GC-MS) revealed that addition of P. pentosaceus YC altered the composition of gut microbiota and increased the content of propionate in fish gut. In vitro studies on propionate's function demonstrated that it suppressed nlrc3 expression and promoted wound healing in Caco-2 cell model.

CONCLUSIONS

The present study reveals that P. pentosaceus YC has the capacity to ameliorate intestinal barrier injury by modulating gut microbiota composition and elevating propionate level. This finding offers a promising strategy for the feed industry to incorporate cottonseed meal into fish feed formulations.

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.

Biodegradation of Free Gossypol by Helicoverpa armigera Carboxylesterase Expressed in Pichia pastoris

Gossypol is a polyphenolic toxic secondary metabolite derived from cotton. Free gossypol in cotton meal is remarkably harmful to animals. Furthermore, microbial degradation of gossypol produces metabolites that reduce feed quality. We adopted an enzymatic method to degrade free gossypol safely and effectively. We cloned the gene cce001a encoding carboxylesterase (CarE) into pPICZαA and transformed it into Pichia pastoris GS115. The target protein was successfully obtained, and CarE CCE001a could effectively degrade free gossypol with a degradation rate of 89%. When esterase was added, the exposed toxic groups of gossypol reacted with different amino acids and amines to form bound gossypol, generating substances with (M + H) m/z ratios of 560.15, 600.25, and 713.46. The molecular formula was C₂₇H₂₈O₁₃, C₃₄H₃₆N₂O₆, and C₄₇H₅₉N₃O₃. The observed instability of the hydroxyl groups caused the substitution and shedding of the group, forming a substance with m/z of 488.26 and molecular formula C₃₁H₃₆O₅. These properties render the CarE CCE001a a valid candidate for the detoxification of cotton meal. Furthermore, the findings help elucidate the degradation process of gossypol in vitro.

Genome-wide association study reveals that GhTRL1 and GhPIN8 affect cotton root development

Two regions located at chromosome A05 and D04 were found to be significantly associated with 0-0.5 mm and 0.5-2 mm diameter roots, respectively, and two candidate genes related to root development were identified. Roots absorb water and nutrients, and play an important role in plant growth. However, there are few genetic developmental studies on cotton root structural traits. In this study, we used 200 upland cotton (Gossypium hirsutum L.) varieties to analyze the phenotypic variation of 43 traits. A total of 2001 related single-nucleotide polymorphism (SNP) sites located within or near 1046 genes were detected through a genome-wide association study (GWAS). The 32 root traits were linked to SNPs that corresponded to 317 nonrepetitive genes. For SNPs associated with root length and 0-0.5 mm diameter root traits, a significant peak appeared on chromosome A05 (between 21.91 and 22.24 Mb). For SNPs associated with root surface area, root volume and 0.5-2 mm diameter root traits, a significant peak appeared on chromosome D04 (between 7.35 and 7.70 Mb). Within these two key regions, SNPs were detected in the promoter and coding regions of two candidate genes, GhTRL1-A05 and GhPIN8-D04. The expression levels of these two genes also changed significantly according to transcriptome sequencing and quantitative real-time PCR (qRT-PCR). After silencing the GhTRL1 and GhPIN8 genes via virus-induced gene silencing (VIGS), we found that the plants expressing TRV2::GhTRL1 and TRV2::GhPIN8 had a reduced root length, surface area. Moreover, the contents of cis-12-oxo-phytodienoic acid (cis-OPDA), isopentenyl adenosine (iPR) and cis-zeatin (cZ) in the roots of the plants expressing TRV2::GhTRL1 decreased. This study contributes to the cultivation and improvement of cotton varieties.

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.

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.

Root Illumination Promotes Seedling Growth and Inhibits Gossypol Biosynthesis in Upland Cotton

Gossypol, a terpenoid compound mainly synthesized in the cotton roots, acts as a phytoalexin in protecting the plants from biotic stress. Roots are critical for both the secondary metabolism and the growth of the plant. Light plays an important role in plant growth and material metabolism, however, the effect of root illumination (RI) on the cotton seedling growth and gossypol metabolism remains unclear. In the present study, the cotton genetic standard line TM-1 and four pairs of near-isogenic lines (NILs) were used as materials to study the impact of RI on cotton seedlings. Results showed that, compared with the cotton seedlings cultivated without RI, the photosynthetic rate, leaf area, and dry weight of roots and leaves were significantly increased, while the gossypol content in leaves and roots was significantly reduced in seedlings cultivated with RI. GO and KEGG enrichment analysis of the differentially expressed genes (DEGs) with and without RI both indicated that photosynthesis and terpenoid biosynthesis-related GO terms and pathways were significantly enriched, the expression profile confirmed that RI positively regulated the photosynthesis system and negatively affected the gossypol biosynthesis pathway in roots. This study revealed the effects of RI on seedlings' growth and gossypol biosynthesis in upland cotton, and provided important insights for the engineering of cotton with low gossypol accumulation.

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.

Cotton roots are the major source of gossypol biosynthesis and accumulation

BACKGROUND

Gossypol is a specific secondary metabolite in Gossypium species. It not only plays a critical role in development and self-protection of cotton plants, but also can be used as important anti-cancer and male contraceptive compound. However, due to the toxicity of gossypol for human beings and monogastric animals, the consumption of cottonseeds was limited. To date, little is known about the gossypol metabolism in cotton plants.

RESULTS

In this study, we found that cotyledon was the primary source of gossypol at the seed germination stage. But thereafter, it was mainly originated from developing roots. Grafting between glanded and glandless cotton as well as sunflower rootstocks and cotton scion revealed that gossypol was mainly synthesized in the root systems of cotton plants. And both glanded and glandless cotton roots had the ability of gossypol biosynthesis. But the pigment glands, the main storage of gossypol, had indirect effects on gossypol biosynthesis. In vitro culture of root and rootless seedling confirmed the strong gossypol biosynthesis ability in root system and the relatively weak gossypol biosynthesis ability in other organs of the seedling. Expression profiling of the key genes involved in the gossypol biosynthetic pathway also supported the root as the major organ of gossypol biosynthesis.

CONCLUSIONS

Our study provide evidence that the cotton root system is the major source of gossypol in both glanded and glandless cottons, while other organs have a relatively weak ability to synthesize gossypol. Gossypol biosynthesis is not directed related to the expression of pigment glands, but the presence of pigment glands is essential for gossypol accumulation. These findings can not only clarify the complex regulation network of gossypol metabolism, but it could also accelerate the crop breeding process with enhanced commercial values.

Genome-wide analysis of genetic variations between dominant and recessive NILs of glanded and glandless cottons

Cotton is an important economic crop in worldwide. It produces fiber for the textile industry and provides cottonseeds with high-quality protein and oil. However, the presence of gossypol limits the utilization of cottonseed. Two pairs of cotton near isogenic lines (NILs) with different pigment glands, i.e., Coker 312 vs Coker 312 W and CCRI12 vs CCRI12W, exhibit different gossypol contents. The glandless traits of Coker 312 W and CCRI12W are controlled by recessive and dominant genes, respectively. However, knowledge regarding the genomic variations in the NILs is limited. Therefore, the NILs genomes were resequenced and the sequencing depths were greater than 34×. Compared with the TM-1 genome, numerous SNPs, Indels, SVs, and CNVs were discovered. KEGG pathway analysis revealed that genes with SNPs and Indels from the recessive NILs and genes with Indels from the dominant NILs shared only one enriched pathway, i.e., the sesquiterpenoid and triterpenoid biosynthesis pathway, which is relevant to gossypol biosynthesis. Expression analysis revealed that key genes with variations that participate in the gossypol biosynthesis and pigment gland formation pathways had different expression patterns among the dominant, recessive glandless and glanded plants. The expression levels in the glanded organs were higher than those in their NILs. Altogether, our results provide deeper insight into cotton NILs with different pigment glands.

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

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

Expression in Pichia pastoris and characterization of two novel dirigent proteins for atropselective formation of gossypol

We established an efficient fed-batch fermentation process for two novel dirigent proteins from cotton plants, GbDIR2 from Gossypium barbadense and GhDIR3 from G. hirsutum, using the engineered Pichia pastoris GlycoSwitch® SuperMan₅ strain to prevent hyperglycosylation. The two (His)₆-tagged proteins were purified by metal-chelate affinity chromatography and obtained in quantities of 12 and 15 mg L^-1 of culture volume, respectively. Glycosylation sites were identified for the native and for the enzymatically deglycosylated proteins by mass spectrometry, confirming five to six of the seven predicted glycosylation sites in the NxS/T sequence context. The predominant glycan structure was Man₅GlcNAc₂ with, however, a significant contribution of Man_4-10GlcNAc₂. Both dirigent proteins (DIRs) mediated the formation of (+)-gossypol by atropselective coupling of hemigossypol radicals. Similar to previously characterized DIRs, GbDIR2 and GhDIR3 lacked oxidizing activity and depended on an oxidizing system (laccase/O₂) for the generation of substrate radicals. In contrast to DIRs involved in the biosynthesis of lignans, glycosylation was not essential for function. Quantitative enzymatic deglycosylation yielded active GbDIR2 and GhDIR3 in excellent purity. The described fermentation process in combination with enzymatic deglycosylation will pave the way for mechanistic and structural studies and, eventually, the application of cotton DIRs in a biomimetic approach towards atropselective biaryl synthesis.

Rhizobacteria activates (+)-δ-cadinene synthase genes and induces systemic resistance in cotton against beet armyworm (Spodoptera exigua)

Gossypol is an important allelochemical produced by the subepidermal glands of some cotton varieties and important for their ability to respond to changing biotic stress by exhibiting antibiosis against some cotton pests. Plant growth-promoting rhizobacteria (PGPR) are root-colonizing bacteria that increase plant growth and often elicit defence against plant pathogens and insect pests. Little is known about the effect of PGPR on cotton plant-insect interactions and the potential biochemical and molecular mechanisms by which PGPR enhance cotton plant defence. Here, we report that PGPR (Bacillus spp.) treated cotton plants showed significantly higher levels of gossypol compared with untreated plants. Similarly, the transcript levels of the genes (i.e. (+)-δ-cadinene synthase gene family) involved in the biosynthesis of gossypol were higher in PGPR-treated plants than in untreated plants. Furthermore, the levels of jasmonic acid, an octadecanoid-derived defence-related phytohormone and the transcript level of jasmonic acid responsive genes were higher in PGPR-treated plants than in untreated plants. Most intriguingly, Spodoptera exigua showed reduced larval feeding and development on PGPR-treated plants. These findings demonstrate that treatment of plants with rhizobacteria may induce significant biochemical and molecular changes with potential ramifications for plant-insect interactions.

Jasmonate-responsive transcription factors regulating plant secondary metabolism

Plants produce a large variety of secondary metabolites including alkaloids, glucosinolates, terpenoids and phenylpropanoids. These compounds play key roles in plant-environment interactions and many of them have pharmacological activity in humans. Jasmonates (JAs) are plant hormones which induce biosynthesis of many secondary metabolites. JAs-responsive transcription factors (TFs) that regulate the JAs-induced accumulation of secondary metabolites belong to different families including AP2/ERF, bHLH, MYB and WRKY. Here, we give an overview of the types and functions of TFs that have been identified in JAs-induced secondary metabolite biosynthesis, and highlight their similarities and differences in regulating various biosynthetic pathways. We review major recent developments regarding JAs-responsive TFs mediating secondary metabolite biosynthesis, and provide suggestions for further studies.

Production and transcriptional regulation of proanthocyanidin biosynthesis in forage legumes

Proanthocyanidins (PA), also known as condensed tannins, contribute to important forage legumes traits including disease resistance and forage quality. PA in forage plants has both positive and negative effects on feed digestibility and animal performance. The analytical methods and their applicability in measuring the contents of PA in forage plants are essential to studies on their nutritional effects. In spite of important breakthroughs in our understanding of the PA biosynthesis, important questions still remain to be answered such as the PA polymerization and transport. Recent advances in the understanding of transcription factor-mediated gene regulation mechanisms in anthocyanin and PA biosynthetic pathway in model plants suggest new approaches for the metabolic engineering of PA in forage plants. The present review will attempt to present the state-of-the-art of research in these areas and provide an update on the production and metabolic engineering of PA in forage plants. We hope that this will contribute to a better understanding of the ways in which PA production to manipulate the content of PA for beneficial effects in forage plants.

Cotton proteomics for deciphering the mechanism of environment stress response and fiber development

Cotton fiber is considered as the backbone of the textile industry. The productivity of cotton crop is severely hampered by the occurrence of pathogens, pests, and various environmental factors. Nevertheless, cotton plant has developed sophisticated mechanisms to respond to environment stresses to avoid detrimental effects on its growth and development. Therefore, understanding the mechanisms of cotton fiber development and environment stress response is of considerable interest for designing agriculture breeding strategies to ensure sustainable productivity. The application of proteomics technologies to advance our knowledge in cotton fiber development and abiotic/biotic stress tolerance has increased dramatically in the last 5years as evidenced by the large amount of publications in this area. This review summarizes the work which has been reported for cotton proteomics and evaluates the findings in context of the approaches that are widely employed with the aim to generate novel insight useful for cotton improvement.

BIOLOGICAL SIGNIFICANCE

Cotton (Gossypium spp.) is considered as the foremost commercially important fiber crop grown all over the world and is deemed as the backbone of the textile industry. Cotton is also an important source of edible oil seed and a nutrient-rich food crop as cottonseed contains high-quality protein and oil. The growth and productivity of cotton crop are often hampered by various biotic stress factors, such as insect pests and pathogens. In addition, cotton plants are frequently subjected to unavoidable environmental factors that cause abiotic stress, such as salt, heat and drought. Proteomic techniques provide one of the best options for understanding the gene function and phenotypic changes during cotton fiber development and stress response. This review first summarizes the work which has been reported for cotton proteomics about cotton fiber development and abiotic/biotic stress tolerance, and also evaluates the findings in context of the approaches that are widely employed during last 5years, with the aim to generate novel insight useful for cotton improvement. This article is part of a Special Issue entitled: Proteomics of non-model organisms.