De Novo Transcriptome and Expression Profile Analysis to Reveal Genes and Pathways Potentially Involved in Cantharidin Biosynthesis in the Blister Beetle Mylabris cichorii

Genome-Wide Profiling of P450 Gene Expression Reveals Caste-Specific and Developmental Patterns in Solenopsis invicta

P450 enzymes are integral to insect physiology, metabolism, hormone regulation, and adaptation to environmental challenges. By leveraging transcriptomic and genomic data, this study characterized the expression of 68 unique P450 genes across developmental stages and castes in the red imported fire ant (Solenopsis invicta), uncovering stage- and caste-specific differential expression patterns. Genes from the CYP4, CYP6, and CYP9 families, known for metabolizing exogenous and endogenous compounds, were highly expressed in early larval stages and minim workers, underscoring their roles in supporting rapid growth, hormone metabolism, colony maintenance, and brood care. The overexpression of CYP4AA1-linked to pheromone production-in queens, female alates, and female alate pupae highlights its critical functions in reproductive dominance, social structure maintenance, and colony dynamics. Here, juvenile hormone biosynthesis genes, including CYP305A1 and CYP315A1, exhibited significant overexpression in later instar larvae and larger workers, emphasizing their roles in development and in fulfilling colony-wide physiological demands. The "Halloween genes" (CYP302A1, CYP306A1, CYP315A1, CYP307A1, and CYP314A1) and CYP18A1 demonstrated dynamic regulation across developmental stages and castes, reflecting their essential contributions to hormonal production and balance throughout S. invicta's lifecycle. These findings offer valuable insights into the molecular and biological mechanisms driving S. invicta's social organization, developmental transitions, physiological adaptations, and evolutionary success. They also provide a foundation for future research into the regulatory pathways governing P450 gene expression and function.

The first chromosome-level genome assembly and transcriptome sequencing provide insights into cantharidin production of the blister beetles

Blister beetles (Coleoptera: Meloidae) produce a natural defensive toxin cantharidin (CTD), which has been used for various cancer treatments and other diseases. Currently, the lack of chromosome-level reference genomes in Meloidae limits further understanding of the mechanism of CTD biosynthesis and environmental adaptation. In this study, the chromosome-level genome assembly of Mylabris phalerata was generated based on PacBio and Hi-C sequencing. This reference genome was about 136.68 Mb in size with contig N50 of 9.17 Mb and composed of 12 chromosomes. In comparison to six other Coleoptera insects, M. phalerata exhibited multiple expanded gene families enriched in juvenile hormone (JH) biosynthetic process pathway, farnesol dehydrogenase activity, and cytochrome P450, which may be related to CTD biosynthesis. Consistently, the transcriptomic analysis suggested the "terpenoid backbone biosynthesis" pathway and "the juvenile hormone" as putative core pathways of CTD biosynthesis and presented eight up-regulated differential expression genes in male adults as candidate genes. It is possible that the restricted feeding niche and lifestyle of M. phalerata were the cause of the gene family's contraction of odorant binding proteins. The ABC transporters (ABCs) related to exporting bound toxins out of the cell and the resistance to the self-secreted toxins (e.g. CTD) were also contracted, possibly due to other self-protection strategies in M. phalerata. A foundation of understanding CTD biosynthesis and environmental adaptation of blister beetles will be established by our reference genome and discoveries.

Natural Occurring Terpene Cyclic Anhydrides: Biosynthetic Origin and Biological Activities

Cyclic acid anhydride is a not very widespread structure in nature, but with a determining role in natural products possessing this functionality in their skeleton. To the best of our knowledge, no revision of terpenes containing cyclic anhydrides has been previously reported. The result was that more than 100 terpenic cyclic anhydrides and related compounds were found to be in need of being reported. This review has been systematically organized by terpene skeletons, from the smallest to largest, describing their sources and bioactivities. In addition, different biosynthetic pathways for their final oxidations, namely, routes A, B and C, leading to the formation of these heterocyclic natural products, have been proposed. We have also included the most plausible precursors of these natural products, which mostly happened to be present in the same natural source. Some molecules derived from terpene cyclic anhydrides, such as their natural imide derivatives, have also been described due to their significant biological activity. In this sense, special attention has been paid to cantharidin because of its historical relevance and its broad bioactivity. A plausible biosynthesis of cantharidin has been proposed for the first time. Finally, cyclic anhydride structures that were firstly assigned as anhydrides and later corrected have been also described.

Functional studies of McSTE24, McCYP305a1, and McJHEH, three essential genes act in cantharidin biosynthesis in the blister beetle (Coleoptera: Meloidae)

Cantharidin is a toxic defensive substance secreted by most blister beetles when attacked. It has been used to treat many complex diseases since ancient times and has recently regained popularity as an anticancer agent. However, the detailed mechanism of the cantharidin biosynthesis has not been completely addressed. In this study, we cloned McSTE24 (encoding STE24 endopeptidase) from terpenoid backbone pathway, McCYP305a1 (encoding cytochrome P450, family 305) and McJHEH [encoding subfamily A, polypeptide 1 and juvenile hormone (JH) epoxide hydrolase] associated to JH synthesis/degradation in the blister beetle Mylabris cichorii (Linnaeus, 1758, Coleoptera: Meloidae). Expression pattern analyses across developmental stages in adult males revealed that the expressions of 3 transcripts were closely linked to cantharidin titer exclusively during the peak period of cantharidin synthesis (20-25 days old). In contrast, at other stages, these genes may primarily regulate different biological processes. When RNA interference with double-stranded RNA suppressed the expressions of the 3 genes individually, significant reductions in cantharidin production were observed in males and also in females following McJHEH knockdown, indicating that these 3 genes might primarily contribute to cantharidin biosynthesis in males, but not in females, while females could self-synthesis a small amount of cantharidin. These findings support the previously hypothesized sexual dimorphism in cantharidin biosynthesis during the adult phase. McCYP305a1 collaborates with its upstream gene McSTE24 in cantharidin biosynthesis, while McJHEH independently regulates cantharidin biosynthesis in males.

A minireview of the medicinal and edible insects from the traditional Chinese medicine (TCM)

Entomoceuticals define a subset of pharmaceuticals derived from insects. The therapeutic effect of insect-derived drugs has been empirically validated by the direct use of various folk medicines originating from three sources in particular: the glandular secretions of insects (e.g., silk, honey, venom), the body parts of the insect or the whole used live or by various processing (e.g., cooked, toasted, ground), and active ingredients extracted from insects or insect-microbe symbiosis. Insects have been widely exploited in traditional Chinese medicine (TCM) relative to other ethnomedicines, especially in the prospect of insect species for medicinal uses. It is noticeable that most of these entomoceuticals are also exploited as health food for improving immune function. In addition, some edible insects are rich in animal protein and have high nutritional value, which are used in the food field, such as insect wine, health supplements and so on. In this review, we focused on 12 insect species that have been widely used in traditional Chinese herbal formulae but have remained less investigated for their biological properties in previous studies. We also combined the entomoceutical knowledge with recent advances in insect omics. This review specifies the underexplored medicinal insects from ethnomedicine and shows their specific medicinal and nutritional roles in traditional medicine.

Unraveling the role of male reproductive tract and haemolymph in cantharidin-exuding Lydus trimaculatus and Mylabris variabilis (Coleoptera: Meloidae): a comparative transcriptomics approach

Background

Meloidae (blister beetles) are known to synthetize cantharidin (CA), a toxic and defensive terpene mainly stored in male accessory glands (MAG) and emitted outward through reflex-bleeding. Recent progresses in understanding CA biosynthesis and production organ(s) in Meloidae have been made, but the way in which self-protection is achieved from the hazardous accumulation and release of CA in blister beetles has been experimentally neglected. To provide hints on this pending question, a comparative de novo assembly transcriptomic approach was performed by targeting two tissues where CA is largely accumulated and regularly circulates in Meloidae: the male reproductive tract (MRT) and the haemolymph. Differential gene expression profiles in these tissues were examined in two blister beetle species, Lydus trimaculatus (Fabricius, 1775) (tribe Lyttini) and Mylabris variabilis (Pallas, 1781) (tribe Mylabrini). Upregulated transcripts were compared between the two species to identify conserved genes possibly involved in CA detoxification and transport.

Results

Based on our results, we hypothesize that, to avoid auto-intoxication, ABC, MFS or other solute transporters might sequester purported glycosylated CA precursors into MAG, and lipocalins could bind CA and mitigate its reactivity when released into the haemolymph during the autohaemorrhaging response. We also found an over-representation in haemolymph of protein-domains related to coagulation and integument repairing mechanisms that likely reflects the need to limit fluid loss during reflex-bleeding.

Conclusions

The de novo assembled transcriptomes of L. trimaculatus and M. variabilis here provided represent valuable genetic resources to further explore the mechanisms employed to cope with toxicity of CA in blister beetle tissues. These, if revealed, might help conceiving safe and effective drug-delivery approaches to enhance the use of CA in medicine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08118-8.

De novo transcriptome assemblies of Epicauta tibialis provide insights into the sexual dimorphism in the production of cantharidin

Blister beetles have medicinal uses for their defensive secretion cantharidin, which has curative effects on many cancers and other diseases. It was demonstrated that sexual dimorphism exists in the production of cantharidin between male and female adults. This study performed a de novo assembly of Epicauta tibialis transcriptomes and analyzed the differentially expressed genes (DEGs) between male and female adults to help to find genes and pathways associated with cantharidin biosynthesis. A total of 99,295,624 paired reads were generated, and more than 7 Gb transcriptome data for each sample were obtained after trimming. The clean data were used to de novo assemble and then cluster into 27,355 unigenes, with a mean length of 1442 bp and an N50 of 2725 bp. Of these, 14,314 (52.33%) unigenes were annotated by protein databases. Differential expression analysis identified 284 differentially expressed genes (DEGs) between male and female adults. Nearly 239 DEGs were up-regulated in male adults than in female adults, while 45 DEGs were down-regulated. The Kyoto Encyclopedia of Gene and Genomes pathway enrichment manifested that seven up-regulated DEGs in male adults were assigned to the terpenoid biosynthesis pathway, to which 19 unigenes were annotated. The DEGs in the terpenoid biosynthesis pathway between male and female adults may be responsible for the sexual dimorphism in cantharidin production. The up-regulated genes assigned to the pathway in male adults may play a significant role in cantharidin biosynthesis, and its biosynthesis process is probably via the mevalonate pathway. The results would be helpful to better understand and reveal the complicated mechanism of the cantharidin biosynthesis.

Transcriptome Sequencing of the Striped Cucumber Beetle, Acalymma vittatum (F.), Reveals Numerous Sex-Specific Transcripts and Xenobiotic Detoxification Genes

Acalymma vittatum (F.), the striped cucumber beetle, is an important pest of cucurbit crops in the contintental United States, damaging plants through both direct feeding and vectoring of a bacterial wilt pathogen. Besides providing basic biological knowledge, biosequence data for A. vittatum would be useful towards the development of molecular biopesticides to complement existing population control methods. However, no such datasets currently exist. In this study, three biological replicates apiece of male and female adult insects were sequenced and assembled into a set of 630,139 transcripts (of which 232,899 exhibited hits to one or more sequences in NCBI NR). Quantitative analyses identified 2898 genes differentially expressed across the male–female divide, and qualitative analyses characterized the insect’s resistome, comprising the glutathione S-transferase, carboxylesterase, and cytochrome P450 monooxygenase families of xenobiotic detoxification genes. In summary, these data provide useful insights into genes associated with sex differentiation and this beetle’s innate genetic capacity to develop resistance to synthetic pesticides; furthermore, these genes may serve as useful targets for potential use in molecular-based biocontrol technologies.

Diversity of Insect Sesquiterpenoid Regulation

Insects are arguably the most successful group of animals in the world in terms of both species numbers and diverse habitats. The sesquiterpenoids juvenile hormone, methyl farnesoate, and farnesoic acid are well known to regulate metamorphosis, reproduction, sexual dimorphism, eusociality, and defense in insects. Nevertheless, different insects have evolved with different sesquiterpenoid biosynthetic pathway as well as products. On the other hand, non-coding RNAs such as microRNAs have been implicated in regulation of many important biological processes, and have recently been explored in the regulation of sesquiterpenoid production. In this review, we summarize the latest findings on the diversity of sesquiterpenoids reported in different groups of insects, as well as the recent advancements in the understanding of regulation of sesquiterpenoid production by microRNAs.

Draft Genome of a Blister Beetle Mylabris aulica

Mylabris aulica is a widely distributed blister beetle of the Meloidae family. It has the ability to synthesize a potent defensive secretion that includes cantharidin, a toxic compound used to treat many major illnesses. However, owing to the lack of genetic studies on cantharidin biosynthesis in M. aulica, the commercial use of this species is less extensive than that of other blister beetle species in China. This study reports a draft assembly and possible genes and pathways related to cantharidin biosynthesis for the M. aulica blister beetle using nanopore sequencing data. The draft genome assembly size was 288.5 Mb with a 467.8 Kb N50, and a repeat content of 50.62%. An integrated gene finding pipeline performed for assembly obtained 16,500 protein coding genes. Benchmarking universal single-copy orthologs assessment showed that this gene set included 94.4% complete Insecta universal single-copy orthologs. Over 99% of these genes were assigned functional annotations in the gene ontology, Kyoto Encyclopedia of Genes and Genomes, or Genbank non-redundant databases. Comparative genomic analysis showed that the completeness and continuity of our assembly was better than those of Hycleus cichorii and Hycleus phaleratus blister beetle genomes. The analysis of homologous orthologous genes and inference from evolutionary history imply that the Mylabris and Hycleus genera are genetically close, have a similar genetic background, and have differentiated within one million years. This M. aulica genome assembly provides a valuable resource for future blister beetle studies and will contribute to cantharidin biosynthesis.

Draft genomes of two blister beetles Hycleus cichorii and Hycleus phaleratus

Background

Commonly known as blister beetles or Spanish fly, there are more than 1500 species in the Meloidae family (Hexapoda: Coleoptera: Tenebrionoidea) that produce the potent defensive blistering agent cantharidin. Cantharidin and its derivatives have been used to treat cancers such as liver, stomach, lung, and esophageal cancers. Hycleus cichorii and Hycleus phaleratus are the most commercially important blister beetles in China due to their ability to biosynthesize this potent vesicant. However, there is a lack of genome reference, which has hindered development of studies on the biosynthesis of cantharidin and a better understanding of its biology and pharmacology.

Results

We report 2 draft genomes and quantified gene sets for the blister beetles H. cichorii and H. phaleratus, 2 complex genomes with >72% repeats and approximately 1% heterozygosity, using Illumina sequencing data. An integrated assembly pipeline was performed for assembly, and most of the coding regions were obtained. Benchmarking universal single-copy orthologs (BUSCO) assessment showed that our assembly obtained more than 98% of the Endopterygota universal single-copy orthologs. Comparison analysis showed that the completeness of coding genes in our assembly was comparable to other beetle genomes such as Dendroctonus ponderosae and Agrilus planipennis. Gene annotation yielded 13 813 and 13 725 protein-coding genes in H. cichorii and H. phaleratus, of which approximately 89% were functionally annotated. BUSCO assessment showed that approximately 86% and 84% of the Endopterygota universal single-copy orthologs were annotated completely in these 2 gene sets, whose completeness is comparable to that of D. ponderosae and A. planipennis.

Conclusions

Assembly of both blister beetle genomes provides a valuable resource for future biosynthesis of cantharidin and comparative genomic studies of blister beetles and other beetles.