Aromatization of natural products by a specialized detoxification enzyme.

In plants, lineage-specific metabolites can be created by activities derived from the catalytic promiscuity of ancestral proteins, although examples of recruiting detoxification systems to biosynthetic pathways are scarce. The ubiquitous glyoxalase (GLX) system scavenges the cytotoxic methylglyoxal, in which GLXI isomerizes the α-hydroxy carbonyl in the methylglyoxal-glutathione adduct for subsequent hydrolysis. We show that GLXIs across kingdoms are more promiscuous than recognized previously and can act as aromatases without cofactors. In cotton, a specialized GLXI variant, SPG, has lost its GSH-binding sites and organelle-targeting signal, and evolved to aromatize cyclic sesquiterpenes bearing α-hydroxyketones to synthesize defense compounds in the cytosol. Notably, SPG is able to transform acetylated deoxynivalenol, the prevalent mycotoxin contaminating cereals and foods. We propose that detoxification enzymes are a valuable source of new catalytic functions and SPG, a standalone enzyme catalyzing complex reactions, has potential for toxin degradation, crop engineering and design of novel aromatics.

Mediation of JA signalling in glandular trichomes by the woolly/SlMYC1 regulatory module improves pest resistance in tomato.

Almost all plants form trichomes, which protect them against insect herbivores by forming a physical barrier and releasing chemical repellents. Glandular trichomes produce a variety of specialized defensive metabolites, including volatile terpenes. Previous studies have shown that the defence hormone jasmonic acid (JA) affects trichome development and induces terpene synthases (TPSs) but the underlying molecular mechanisms remain unclear. Here, we characterized a loss-of-function allele of the HD-ZIP IV transcription factor woolly (wo) and analysed its role in mediating JA signalling in tomato. We showed that knockout of wo led to extensive trichome defects, including structural and functional changes in type VI glandular trichomes, and a dramatic reduction in terpene levels. We further found that wo directly binds to TPS gene promoters to recruit SlMYC1, a JA signalling modulator, and that together these transcription factors promote terpene biosynthesis in tomato trichomes. The wo/SlMYC1 regulatory module is inhibited by SlJAZ2 through a competitive binding mechanism, resulting in a fine-tuned JA response in tomato trichomes. Enhanced expression of SlMYC1 substantially increased terpene levels and improved tomato resistance to spider mites. Interestingly, we also found that SlMYC1 plays an additional role in glandular cell division and expansion in type VI trichomes, independent of JA. Together, our results reveal a novel, JA-mediated regulatory mechanism that promotes insect resistance in tomato.

1,10/1,11-Cyclization catalyzed by diverged plant sesquiterpene synthases is dependent on a single residue.

The exquisite chemodiversity of terpenoids is the product of the large diverse terpene synthase (TPS) superfamily. Here, by using structural and phylogenetic analyses and site-directed mutagenesis, we identified a residue (Cys440 in Nicotiana tabacum 5-epi-aristolochene synthase) proximal to an ion-binding motif common to all TPSs and named the preNSE/DTE residue, which determines the product specificity of sesquiterpene synthases from different plant species. In sesquiterpene synthases catalyzing 1,10-cyclization (1,10-cyclases) of farnesyl diphosphate, mutation of the residue in both specific and promiscuous 1,10-cyclases from different lineages leads to the accumulation of monocyclic germacrene A-11-ol, which is "short-circuited" from complex cyclization cascades, suggesting a key role of this residue in generating the first common intermediate of 1,10-cyclization. Altering this residue in a specific 1,11-cyclase results in alternative 1,10-cyclization products. Moreover, the preNSE/DTE residue can be harnessed to engineer highly specific sesquiterpene synthases for an improved proportion of high-value terpenoids, such as patchoulol, a main constituent of several traditional Chinese medicines that could treat SARS-CoV-2.

Characterization of gossypol biosynthetic pathway.

Gossypol and related sesquiterpene aldehydes in cotton function as defense compounds but are antinutritional in cottonseed products. By transcriptome comparison and coexpression analyses, we identified 146 candidates linked to gossypol biosynthesis. Analysis of metabolites accumulated in plants subjected to virus-induced gene silencing (VIGS) led to the identification of four enzymes and their supposed substrates. In vitro enzymatic assay and reconstitution in tobacco leaves elucidated a series of oxidative reactions of the gossypol biosynthesis pathway. The four functionally characterized enzymes, together with (+)-δ-cadinene synthase and the P450 involved in 7-hydroxy-(+)-δ-cadinene formation, convert farnesyl diphosphate (FPP) to hemigossypol, with two gaps left that each involves aromatization. Of six intermediates identified from the VIGS-treated leaves, 8-hydroxy-7-keto-δ-cadinene exerted a deleterious effect in dampening plant disease resistance if accumulated. Notably, CYP71BE79, the enzyme responsible for converting this phytotoxic intermediate, exhibited the highest catalytic activity among the five enzymes of the pathway assayed. In addition, despite their dispersed distribution in the cotton genome, all of the enzyme genes identified show a tight correlation of expression. Our data suggest that the enzymatic steps in the gossypol pathway are highly coordinated to ensure efficient substrate conversion.

Genome-wide characterization of the GRF family and their roles in response to salt stress in Gossypium.

BACKGROUND: Cotton (Gossypium spp.) is the most important world-wide fiber crop but salt stress limits cotton production in coastal and other areas. Growth regulation factors (GRFs) play regulatory roles in response to salt stress, but their roles have not been studied in cotton under salt stress. RESULTS: We identified 19 GRF genes in G. raimondii, 18 in G. arboreum, 34 in G. hirsutum and 45 in G. barbadense, respectively. These GRF genes were phylogenetically analyzed leading to the recognition of seven GRF clades. GRF genes from diploid cottons (G. raimondii and G. arboreum) were largely retained in allopolyploid cotton, with subsequent gene expansion in G. barbadense relative to G. hirsutum. Most G. hirsutum GRF (GhGRF) genes are preferentially expressed in young and growing tissues. To explore their possible role in salt stress, we used qRT-PCR to study expression responses to NaCl treatment, showing that five GhGRF genes were down-regulated in leaves. RNA-seq experiments showed that seven GhGRF genes exhibited decreased expression in leaves under NaCl treatment, three of which (GhGRF3, GhGRF4, and GhGRF16) were identified by both RNA-seq and qRT-PCR. We also identified six and three GRF genes that exhibit decreased expression under salt stress in G. arboreum and G. barbadense, respectively. Consistent with its lack of leaf withering or yellowing under the salt treatment conditions, G. arboreum had better salt tolerance than G. hirsutum and G. barbadense. Our results suggest that GRF genes are involved in salt stress responses in Gossypium. CONCLUSION: In summary, we identified candidate GRF genes that were involved in salt stress responses in cotton.

Core cis-element variation confers subgenome-biased expression of a transcription factor that functions in cotton fiber elongation.

Cotton cultivars have evolved to produce extensive, long, seed-born fibers important for the textile industry, but we know little about the molecular mechanism underlying spinnable fiber formation. Here, we report how PACLOBUTRAZOL RESISTANCE 1 (PRE1) in cotton, which encodes a basic helix-loop-helix (bHLH) transcription factor, is a target gene of spinnable fiber evolution. Differential expression of homoeologous genes in polyploids is thought to be important to plant adaptation and novel phenotypes. PRE1 expression is specific to cotton fiber cells, upregulated during their rapid elongation stage and A-homoeologous biased in allotetraploid cultivars. Transgenic studies demonstrated that PRE1 is a positive regulator of fiber elongation. We determined that the natural variation of the canonical TATA-box, a regulatory element commonly found in many eukaryotic core promoters, is necessary for subgenome-biased PRE1 expression, representing a mechanism underlying the selection of homoeologous genes. Thus, variations in the promoter of the cell elongation regulator gene PRE1 have contributed to spinnable fiber formation in cotton. Overexpression of GhPRE1 in transgenic cotton yields longer fibers with improved quality parameters, indicating that this bHLH gene is useful for improving cotton fiber quality.

Interaction between two timing microRNAs controls trichome distribution in Arabidopsis.

The miR156-targeted squamosa promoter binding protein like (SPL) transcription factors function as an endogenous age cue in regulating plant phase transition and phase-dependent morphogenesis, but the control of SPL output remains poorly understood. In Arabidopsis thaliana the spatial pattern of trichome is a hallmark of phase transition and governed by SPLs. Here, by dissecting the regulatory network controlling trichome formation on stem, we show that the miR171-targeted lost meristems 1 (LOM1), LOM2 and LOM3, encoding GRAS family members previously known to maintain meristem cell polarity, are involved in regulating the SPL activity. Reduced LOM abundance by overexpression of miR171 led to decreased trichome density on stems and floral organs, and conversely, constitutive expression of the miR171-resistant LOM (rLOM) genes promoted trichome production, indicating that LOMs enhance trichome initiation at reproductive stage. Genetic analysis demonstrated LOMs shaping trichome distribution is dependent on SPLs, which positively regulate trichome repressor genes TRICHOMELESS 1 (TCL1) and TRIPTYCHON (TRY). Physical interaction between the N-terminus of LOMs and SPLs underpins the repression of SPL activity. Importantly, other growth and developmental events, such as flowering, are also modulated by LOM-SPL interaction, indicating a broad effect of the LOM-SPL interplay. Furthermore, we provide evidence that MIR171 gene expression is regulated by its targeted LOMs, forming a homeostatic feedback loop. Our data uncover an antagonistic interplay between the two timing miRNAs in controlling plant growth, phase transition and morphogenesis through direct interaction of their targets.

Targeting insect mitochondrial complex I for plant protection.

Plant engineered to express double-stranded RNA (dsRNA) can target the herbivorous insect gene for silencing. Although mounting evidence has emerged to support feasibility of this new pest control technology, field application is slow largely due to lack of potent targets. Here, we show that suppression of the gene encoding NDUFV2, a subunit of mitochondrial complex I that catalyses NADH dehydrogenation in respiratory chain, was highly lethal to insects. Feeding cotton bollworm (Helicoverpa armigera) larvae with transgenic cotton tissues expressing NDUFV2 dsRNA led to mortality up to 80% within 5 days, and almost no larvae survived after 7 days of feeding, due to the altered mitochondrial structure and activity. Transcriptome comparisons showed a drastic repression of dopa decarboxylase genes. Reciprocal assays with Asian corn borer (Ostrinia furnacalis), another lepidopteran species, revealed the sequence-specific effect of NDUFV2 suppression. Furthermore, the hemipteran lugus Apolygus lucorum was also liable to NDUFV2 repression. These data demonstrate that the mitochondrial complex I is a promising target with both sequence specificity and wide applicability for the development of new-generation insect-proof crops.

Functional characterization of a basic helix-loop-helix (bHLH) transcription factor GhDEL65 from cotton (Gossypium hirsutum).

Cotton fiber is proposed to share some similarity with the Arabidopsis thaliana leaf trichome, which is regulated by the MYB-bHLH-WD40 transcription complex. Although several MYB transcription factors and WD40 family proteins in cotton have been characterized, little is known about the role of bHLH family proteins in cotton. Here, we report that GhDEL65, a bHLH protein from cotton (Gossypium hirsutum), is a functional homologue of Arabidopsis GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) in regulating trichome development. Transcripts of GhDEL65 were detected in 0 ∼ 1 days post-anthesis (DPA) ovules and abundant in 3-DPA fibers, implying that GhDEL65 may act in early fiber development. Ectopic expression of GhDEL65 in Arabidopsis gl3 egl3 double mutant partly rescued the trichome development, and constitutive expression of GhDEL65 in wild-type plants led to increased trichome density on rosette leaves and stems, mainly by activating the transcription of two key positive regulators of trichome development, GLABRA1 (GL1) and GLABRA2 (GL2), and suppressed the expression of a R3 single-repeat MYB factor TRIPTYCHON (TRY). GhDEL65 could interact with cotton R2R3 MYB transcription factors GhMYB2 and GhMYB3, as well as the WD40 protein GhTTG3, suggesting that the MYB-bHLH-WD40 protein complex also exists in cotton fiber cell, though its function in cotton fiber development awaits further investigation.

Arabidopsis MYC2 interacts with DELLA proteins in regulating sesquiterpene synthase gene expression.

Arabidopsis thaliana flowers emit volatile terpenes, which may function in plant-insect interactions. Here, we report that Arabidopsis MYC2, a basic helix-loop-helix transcription factor, directly binds to promoters of the sesquiterpene synthase genes TPS21 and TPS11 and activates their expression. Expression of TPS21 and TPS11 can be induced by the phytohormones gibberellin (GA) and jasmonate (JA), and both inductions require MYC2. The induction of TPS21 and TPS11 results in increased emission of sesquiterpene, especially (E)-ß-caryophyllene. DELLAs, the GA signaling repressors, negatively affect sesquiterpene biosynthesis, as the sesquiterpene synthase genes were repressed in plants overaccumulating REPRESSOR OF GA1-3 (RGA), one of the Arabidopsis DELLAs, and upregulated in a penta DELLA-deficient mutant. Yeast two-hybrid and coimmunoprecipitation assays demonstrated that DELLAs, represented by RGA, directly interact with MYC2. In yeast cells, the N terminus of MYC2 was responsible for binding to RGA. MYC2 has been proposed as a major mediator of JA signaling and crosstalk with abscisic acid, ethylene, and light signaling pathways. Our results demonstrate that MYC2 is also connected to GA signaling in regulating a subset of genes. In Arabidopsis inflorescences, it integrates both GA and JA signals into transcriptional regulation of sesquiterpene synthase genes and promotes sesquiterpene production.

Rational engineering of plasticity residues of sesquiterpene synthases from Artemisia annua: product specificity and catalytic efficiency.

Most TPSs (terpene synthases) contain plasticity residues that are responsible for diversified terpene products and functional evolution, which provide a potential for improving catalytic efficiency. Artemisinin, a sesquiterpene lactone from Artemisia annua L., is widely used for malaria treatment and progress has been made in engineering the production of artemisinin or its precursors. In the present paper, we report a new sesquiterpene synthase from A. annua, AaBOS (A. annua α-bisabolol synthase), which has high sequence identity with AaADS (A. annua amorpha-4,11-diene synthase), a key enzyme in artemisinin biosynthesis. Comparative analysis of the two enzymes by domain-swapping and structure-based mutagenesis led to the identification of several plasticity residues, whose alteration changed the product profile of AaBOS to include γ-humulene as the major product. To elucidate the underlying mechanisms, we solved the crystal structures of AaBOS and a γ-humulene-producing AaBOS mutant (termed AaBOS-M2). Among the plasticity residues, position 399, located in the substrate-binding pocket, is crucial for both enzymes. In AaBOS, substitution of threonine for leucine (AaBOSL339T) is required for γ-humulene production; whereas in AaADS, replacing the threonine residue with serine (AaADST399S) resulted in a substantial increase in the activity of amorpha-4,11-diene production, probably as a result of accelerated product release. The present study demonstrates that substitution of plasticity residues has potential for improving catalytic efficiency of the enzyme.

miRNA control of vegetative phase change in trees.

After germination, plants enter juvenile vegetative phase and then transition to an adult vegetative phase before producing reproductive structures. The character and timing of the juvenile-to-adult transition vary widely between species. In annual plants, this transition occurs soon after germination and usually involves relatively minor morphological changes, whereas in trees and other perennial woody plants it occurs after months or years and can involve major changes in shoot architecture. Whether this transition is controlled by the same mechanism in annual and perennial plants is unknown. In the annual forb Arabidopsis thaliana and in maize (Zea mays), vegetative phase change is controlled by the sequential activity of microRNAs miR156 and miR172. miR156 is highly abundant in seedlings and decreases during the juvenile-to-adult transition, while miR172 has an opposite expression pattern. We observed similar changes in the expression of these genes in woody species with highly differentiated, well-characterized juvenile and adult phases (Acacia confusa, Acacia colei, Eucalyptus globulus, Hedera helix, Quercus acutissima), as well as in the tree Populus x canadensis, where vegetative phase change is marked by relatively minor changes in leaf morphology and internode length. Overexpression of miR156 in transgenic P. x canadensis reduced the expression of miR156-targeted SPL genes and miR172, and it drastically prolonged the juvenile phase. Our results indicate that miR156 is an evolutionarily conserved regulator of vegetative phase change in both annual herbaceous plants and perennial trees.

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.

Cysteine protease enhances plant-mediated bollworm RNA interference.

Oral ingestion of plant-expressed double stranded RNA (dsRNA) triggers target gene suppression in insect. An important step of this process is the transmission of dsRNA from plant to midgut cells. Insect peritrophic matrix (PM) presents a barrier that prevents large molecules from entering midgut cells. Here, we show that uptake of plant cysteine proteases, such as GhCP1 from cotton (Gossypium hirsutum) and AtCP2 from Arabidopsis, by cotton bollworm (Helicoverpa armigera) larvae resulted in attenuating the PM. When GhCP1 or AtCP2 pre-fed larvae were transferred to gossypol-containing diet, the bollworm accumulated higher content of gossypol in midgut. Larvae previously ingested GhCP1 or AtCP2 were more susceptible to infection by Dendrolimus punctatus cytoplasmic polyhedrosis virus (DpCPV), a dsRNA virus. Furthermore, the pre-fed larvae exhibited enhanced RNAi effects after ingestion of the dsRNA-expressing plant. The bollworm P450 gene CYP6AE14 is involved in the larval tolerance to gossypol; cotton plants producing dsRNA of CYP6AE14 (dsCYP6AE14) were more resistant to bollworm feeding (Mao et al. in Transgenic Res 20:665-673, 2011). We found that cotton plants harboring both 35S:dsCYP6AE14 and 35S:GhCP1 were better protected from bollworm than either of the single-transgene lines. Our results demonstrate that plant cysteine proteases, which have the activity of increasing PM permeability, can be used to improve the plant-mediated RNAi against herbivorous insects.

Temporal control of trichome distribution by microRNA156-targeted SPL genes in Arabidopsis thaliana.

The production and distribution of plant trichomes is temporally and spatially regulated. After entering into the flowering stage, Arabidopsis thaliana plants have progressively reduced numbers of trichomes on the inflorescence stem, and the floral organs are nearly glabrous. We show here that SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) genes, which define an endogenous flowering pathway and are targeted by microRNA 156 (miR156), temporally control the trichome distribution during flowering. Plants overexpressing miR156 developed ectopic trichomes on the stem and floral organs. By contrast, plants with elevated levels of SPLs produced fewer trichomes. During plant development, the increase in SPL transcript levels is coordinated with the gradual loss of trichome cells on the stem. The MYB transcription factor genes TRICHOMELESS1 (TCL1) and TRIPTYCHON (TRY) are negative regulators of trichome development. We show that SPL9 directly activates TCL1 and TRY expression through binding to their promoters and that this activation is independent of GLABROUS1 (GL1). The phytohormones cytokinin and gibberellin were reported to induce trichome formation on the stem and inflorescence via the C2H2 transcription factors GIS, GIS2, and ZFP8, which promote GL1 expression. We show that the GIS-dependent pathway does not affect the regulation of TCL1 and TRY by miR156-targeted SPLs, represented by SPL9. These results demonstrate that the miR156-regulated SPLs establish a direct link between developmental programming and trichome distribution.

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.

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.

Transcriptional regulation of plant secondary metabolism.

Plant secondary metabolites play critical roles in plant-environment interactions. They are synthesized in different organs or tissues at particular developmental stages, and in response to various environmental stimuli, both biotic and abiotic. Accordingly, corresponding genes are regulated at the transcriptional level by multiple transcription factors. Several families of transcription factors have been identified to participate in controlling the biosynthesis and accumulation of secondary metabolites. These regulators integrate internal (often developmental) and external signals, bind to corresponding cis-elements--which are often in the promoter regions--to activate or repress the expression of enzyme-coding genes, and some of them interact with other transcription factors to form a complex. In this review, we summarize recent research in these areas, with an emphasis on newly-identified transcription factors and their functions in metabolism regulation.

The sage genome provides insight into the evolutionary dynamics of diterpene biosynthesis gene cluster in plants.

The widely cultivated medicinal and ornamental plant sage (Salvia officinalis L.) is an evergreen shrub of the Lamiaceae family, native to the Mediterranean. We assembled a high-quality sage genome of 480 Mb on seven chromosomes, and identified a biosynthetic gene cluster (BGC) encoding two pairs of diterpene synthases (diTPSs) that, together with the cytochromes P450 (CYPs) genes located inside and outside the cluster, form two expression cascades responsible for the shoot and root diterpenoids, respectively, thus extending BGC functionality from co-regulation to orchestrating metabolite production in different organs. Phylogenomic analysis indicates that the Salvia clades diverged in the early Miocene. In East Asia, most Salvia species are herbaceous and accumulate diterpenoids in storage roots. Notably, in Chinese sage S. miltiorrhiza, the diterpene BGC has contracted and the shoot cascade has been lost. Our data provide genomic insights of micro-evolution of growth type-associated patterning of specialized metabolite production in plants.

Cotton plants expressing CYP6AE14 double-stranded RNA show enhanced resistance to bollworms.

RNA interference (RNAi) plays an important role in regulating gene expression in eukaryotes. Previously, we generated Arabidopsis and tobacco plants expressing double-stranded RNA (dsRNA) targeting a cotton bollworm (Helicoverpa armigera) P450 gene, CYP6AE14. Bollworms fed on transgenic dsCYP6AE14 plants showed suppressed CYP6AE14 expression and reduced growth on gossypol-containing diet (Mao et al., in Nat Biotechnol 25: 1307-1313, 2007). Here we report generation and analysis of dsRNA-expressing cotton (Gossypium hirsutum) plants. Bollworm larvae reared on T2 plants of the ds6-3 line exhibited drastically retarded growth, and the transgenic plants were less damaged by bollworms than the control. Quantitative reverse-transcription polymerase chain reaction (RT-PCR) showed that the CYP6AE14 expression level was reduced in the larvae as early as 4 h after feeding on the transgenic plants; accordingly, the CYP6AE14 protein level dropped. These results demonstrated that transgenic cotton plants expressing dsCYP6AE14 acquired enhanced resistance to cotton bollworms, and that RNAi technology can be used for engineering insect-proof cotton cultivar.