Integrating RNA-seq with functional expression to analyze the regulation and characterization of genes involved in monoterpenoid biosynthesis in Nepeta tenuifolia Briq

Growth inhibitory, immunosuppressive, cytotoxic, and genotoxic effects of γ-terpinene on Zeugodacus cucurbitae (Coquillett) (Diptera: Tephritidae)

γ-Terpinene, a monoterpene widely present in essential oils of many medicinal and aromatic plants with numerous biological properties, was evaluated for its insecticidal activity against melon fruit fly, Zeugodacus cucurbitae (Coquillett). Different concentrations (5, 25, 125, 625, and 3125 ppm) of γ-terpinene along with control were fed to larvae of melon fly. The number of pupae formed and adults emerged declined significantly after treatment. Morphologically deformed adults and pupae were also observed. The developmental duration too prolonged in treated larvae. Food assimilated, mean relative growth rate, larval weight gain, and pupal weight also declined. In the larvae treated with LC30 and LC50 concentrations, there was a decline in the titers of phenoloxidase and total hemocyte count, and variations were observed in the differential hemocyte count, suggesting an immunosuppressive effect of γ-terpinene on melon fly. Both concentrations also led to an increase in the apoptotic and necrotic cells as well as decrease in the viable hemocytes in the circulating hemolymph of treated larvae. Comet parameters (tail length, % tail DNA, tail moment, and olive tail moment) of γ-terpinene fed larvae increased significantly. Given the observed effects of γ-terpinene on normal developmental and nutritional physiology, its immunosuppressive properties, and its potential for genome damage, it can be considered for incorporation into integrated pest management strategies for controlling Z. cucurbitae.

A chromosome-level genome assembly reveals that a bipartite gene cluster formed via an inverted duplication controls monoterpenoid biosynthesis in Schizonepeta tenuifolia

Biosynthetic gene clusters (BGCs) are regions of a genome where genes involved in a biosynthetic pathway are in proximity. The origin and evolution of plant BGCs as well as their role in specialized metabolism remain largely unclear. In this study, we have assembled a chromosome-scale genome of Japanese catnip (Schizonepeta tenuifolia) and discovered a BGC that contains multiple copies of genes involved in four adjacent steps in the biosynthesis of p-menthane monoterpenoids. This BGC has an unprecedented bipartite structure, with mirrored biosynthetic regions separated by 260 kilobases. This bipartite BGC includes identical copies of a gene encoding an old yellow enzyme, a type of flavin-dependent reductase. In vitro assays and virus-induced gene silencing revealed that this gene encodes the missing isopiperitenone reductase. This enzyme evolved from a completely different enzyme family to isopiperitenone reductase from closely related Mentha spp., indicating convergent evolution of this pathway step. Phylogenomic analysis revealed that this bipartite BGC has emerged uniquely in the S. tenuifolia lineage and through insertion of pathway genes into a region rich in monoterpene synthases. The cluster gained its bipartite structure via an inverted duplication. The discovered bipartite BGC for p-menthane biosynthesis in S. tenuifolia has similarities to the recently described duplicated p-menthane biosynthesis gene pairs in the Mentha longifolia genome, providing an example of the convergent evolution of gene order. This work expands our understanding of plant BGCs with respect to both form and evolution, and highlights the power of BGCs for gene discovery in plant biosynthetic pathways.

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.

Genome-Wide Identification and Characterization of CPR5 Genes in Gossypium Reveals Their Potential Role in Trichome Development

Trichomes protect plants against insects, microbes, herbivores, and abiotic damages and assist seed dispersal. The function of CPR5 genes have been found to be involved in the trichome development but the research on the underlying genetic and molecular mechanisms are extremely limited. Herein, genome wide identification and characterization of CPR5 genes was performed. In total, 26 CPR5 family members were identified in Gossypium species. Phylogenetic analysis, structural characteristics, and synteny analysis of CPR5s showed the conserved evolution relationships of CPR5. The promoter analysis of CPR5 genes revealed hormone, stress, and development-related cis-elements. Gene ontology (GO) enrichment analysis showed that the CPR5 genes were largely related to biological regulation, developmental process, multicellular organismal process. Protein-protein interaction analysis predicted several trichome development related proteins (SIM, LGO, and GRL) directly interacting with CPR5 genes. Further, nine putative Gossypium-miRNAs were also identified, targeting Gossypium CPR5 genes. RNA-Seq data of G. arboreum (with trichomes) and G. herbaceum (with no trichomes) was used to perform the co-expression network analysis. GheCPR5.1 was identified as a hub gene in a co-expression network analysis. RT-qPCR of GheCPR5.1 gene in different tissues suggests that this gene has higher expressions in the petiole and might be a key candidate involved in the trichome development. Virus induced gene silencing of GheCPR5.1 (Ghe02G17590) confirms its role in trichome development and elongation. Current results provide proofs of the possible role of CPR5 genes and provide preliminary information for further studies of GheCPR5.1 functions in trichome development.

Comparison of Pulegone and Estragole Chemotypes Provides New Insight Into Volatile Oil Biosynthesis of Agastache rugosa

The aerial parts of Agastache rugosa are rich in essential oils containing monoterpenoids, phenylpropanoids, and aromatic compounds. These are used as herbs, perfume plants, and ornamental plants. Based on the difference in the constituents of the essential oil, A. rugosa is divided into pulegone and estragole chemotypes, but the mechanism of key metabolite biosynthesis in these two A. rugosa chemotypes remains unclear. In this study, we compared the morphological differences, metabolite constituents, and transcriptomic data between the two chemotypes of A. rugosa. Monoterpenoid was the main compound in the pulegone chemotype, and phenylpropanoid was the main compound in the estragole chemotype; however, limonene was detected in both chemotypes. Furthermore, 46 genes related to pulegone and estragole biosynthesis were identified. Limonene synthase, limonene-3-hydroxylase, and isopiperitenol dehydrogenase were upregulated in the pulegone chemotype, while phenylalanine ammonia-lyase, 4-coumarate: CoA ligase, CYP73A, coumaroyl-aldehyde dehydrogenase, and eugenol synthase were downregulated in the pulegone chemotype. We identified chavicol methyl transferase and limonene-3-hydroxylase in A. rugosa. This work not only provides the difference in morphology and metabolites in pulegone and estragole chemotypes, but also offers a molecular mechanism of volatile oil biosynthesis, which could be a basis for specialized metabolites in specialized chemotypes.

Selectivity of enzymes involved in the formation of opposite enantiomeric series of p-menthane monoterpenoids in peppermint and Japanese catnip

Peppermint (Mentha x piperita L.) and Japanese catnip (Schizonepeta tenuifolia (Benth.) Briq.) accumulate p-menthane monoterpenoids with identical functionalization patterns but opposite stereochemistry. In the present study, we investigate the enantioselectivity of multiple enzymes involved in monoterpenoid biosynthesis in these species. Based on kinetic assays, mint limonene synthase, limonene 3-hydroxylase, isopiperitenol dehydrogenase, isopiperitenone reductase, and menthone reductase exhibited significant enantioselectivity toward intermediates of the pathway that proceeds through (-)-4S-limonene. Limonene synthase, isopiperitenol dehydrogenase and isopiperitenone reductase of Japanese catnip preferred intermediates of the pathway that involves (+)-4R-limonene, whereas limonene 3-hydroxylase was not enantioselective, and the activities of pulegone reductase and menthone reductase were too low to acquire meaningful kinetic data. Molecular modeling studies with docked ligands generally supported the experimental data obtained with peppermint enzymes, indicating that the preferred enantiomer was aligned well with the requisite cofactor and amino acid residues implicated in catalysis. A striking example for enantioselectivity was peppermint (-)-menthone reductase, which binds (-)-menthone with exquisite affinity but was predicted to bind (+)-menthone in a non-productive orientation that positions its carbonyl functional group at considerable distance to the NADPH cofactor. The work presented here lays the groundwork for structure-function studies aimed at unraveling how enantioselectivity evolved in closely related species of the Lamiaceae and beyond.