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.

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.