Comparative genomics provides insights into molecular adaptation to hypermetamorphosis and cantharidin metabolism in blister beetles (Coleoptera: Meloidae)

Abstract

Blister beetles (Coleoptera: Meloidae) are currently subdivided into three subfamilies: Eleticinae (a basal group), Nemognathinae, and Meloinae. These are all characterized by the endogenous production of the defensive terpene cantharidin (CA), whereas the two most derived subfamilies show a hypermetamorphic larval development. Here, we provide novel draft genome assemblies of five species sampled across the three blister beetle subfamilies (Iselma pallidipennis, Stenodera caucasica, Zonitis immaculata, Lydus trimaculatus, and Mylabris variabilis) and performed a comparative analysis with other available Meloidae genomes and the closely-related canthariphilous species (Pyrochroa serraticornis) to disclose adaptations at a molecular level. Our results highlighted the expansion and selection of genes potentially responsible for CA production and metabolism, as well as its mobilization and vesicular compartmentalization. Furthermore, we observed adaptive selection patterns and gain of genes devoted to epigenetic regulation, development, and morphogenesis, possibly related to hypermetamorphosis. We hypothesize that most genetic adaptations occurred to support both CA biosynthesis and hypermetamorphosis, two crucial aspects of Meloidae biology that likely contributed to their evolutionary success.

Keywords

genome evolution / juvenile hormone / larval development / terpene production

References

[1] AhnSJ,VogelH, HeckelDG (2012). Comparative analysis of the UDP-glycosyltransferase multigene family in insects. Insect Biochemistry and Molecular Biology 42,133–147.

[2] ArreseEL,Soulages JL (2010). Insect fat body: energy, metabolism, and regulation. Annual Review of Entomology 55,207–225.

[3] BarskiOA,Tipparaju SM,BhatnagarA (2008). The aldo-keto reductase superfamily and its role in drug metabolism and detoxification. Drug Metabolism Reviews 40,553–624.

[4] BolgerAM,LohseM, UsadelB (2014). Trimmomatic: A flexible trimmer for Illumina Sequence Data. Bioinformatics 30,2114–2120.

[5] BolognaMA,Di Giulio A (2011). Biological and morphological adaptations in the pre-imaginal phases of the beetle family Meloidae. Atti Accademia Nazionale Italiana di Entomologia 59,141–152.

[6] BolognaMA,Di Giulio A,PintoJD (2002). Review of the genus Stenodera with a description of its first instar larva (Coleoptera: Meloidae). European Journal of Entomology 99,299–313.

[7] BolognaMA,Fattorini S,PintoJD (2001). Review of the primitive blister beetle genus Iselma with description of the first instar larva (Coleoptera: Tenebrionoidea: Meloidae). African Entomology 9,105–129.

[8] BolognaMA,TurcoF, PintoJD (2010). 11.19 Meloidae Gyllenhal, 1810. In: Leschen RAB,Beutel RG,Lawrence JF,Kristensen NP, Beutel RG, eds. Handbook of Zoology, Arthropoda: Insecta, vol. 2: Coleoptera, Beetles: Morphology and Systematics (Elateroidea, Bostrichiformia, Cucujiformia partim). The Gruyter, Berlin/New York, pp. 681–693.

[9] BuchfinkB,ReuterK, DrostHG (2021). Sensitive protein alignments at tree-of-life scale using DIAMOND. Nature Methods 18,366–368.

[10] CaiC,TihelkaE, GiacomelliM et al. (2022). Integrated phylogenomics and fossil data illuminate the evolution of beetles. Royal Society Open Science 9,211771.

[11] CarrelJE,EisnerT (1974). Cantharidin: Potent feeding deterrent to insects. Science 183,755–757.

[12] ChenH,ZhuYC, WhitworthRJ,Reese JC et al. (2013). Serine and cysteine protease-like genes in the genome of a gall midge and their interactions with host plant genotypes. Insect Biochemistry and Molecular Biology 43,701–711.

[13] CruseC,MouralTW, ZhuF (2023). Dynamic roles of insect carboxyl/cholinesterases in chemical adaptation. Insects 14,194.

[14] EisnerT,SmedleySR, YoungDK et al. (1996). Chemical basis of courtship in a beetle (Neopyrochroa flabellata): Cantharidin as precopulatory “enticing” agent. PNAS 93,6494–6498.

[15] EmmsDM,KellyS (2019). OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biology 20,238.

[16] FratiniE,RossiMN, SpagoniL et al. (2022). Molecular characterization of Kunitz-type protease inhibitors from blister beetles (Coleoptera, Meloidae). Biomolecules 12,988.

[17] FratiniE,Salvemini M,LombardoF et al. (2021). Unraveling the role of male reproductive tract and haemolymph in cantharidin-exuding Lydus trimaculatus and Mylabris variabilis (Coleoptera: Meloidae): A comparative transcriptomics approach. BMC Genomics 22,808.

[18] García-FernándezR,PeigneurS,PonsT et al. (2016). The Kunitz-type protein ShPI-1 inhibits serine proteases and voltage-gated potassium channels. Toxins 8,110.

[19] HaasBJ,DelcherAL, MountSM et al. (2003). Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Research 31,5654–5666.

[20] HaasBJ,Salzberg SL,ZhuW et al. (2008). Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biology 9,R7.

[21] HarrisonMC,HammondRL, MallonEB (2015). Reproductive workers show queenlike gene expression in an intermediately eusocial insect, the buff-tailed bumble bee Bombus terrestris. Molecular Ecology 24,3043–3063.

[22] JiangM,LüSM, QiZY et al. (2019). Characterized cantharidin distribution and related gene expression patterns in tissues of blister beetles,Epicauta chinensis. Insect Science 26,240–250.

[23] KajitaniR,Toshimoto K,NoguchiH et al. (2014). Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Research 24,1384–1395.

[24] KatohK,Standley DM (2013). MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Molecular Biology and Evolution 30,772–780.

[25] KolbergL,Raudvere U,KuzminI,AdlerP,ViloJ, PetersonH (2023). g:Profiler—interoperable web service for functional enrichment analysis and gene identifier mapping. Nucleic Acids Research 51,W207–W212.

[26] KozlovAM,DarribaD, FlouriT et al. (2019). RAxML-NG: A fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35,4453–4455.

[27] LiQ,WangD, LvS et al. (2014). Comparative proteomics and expression analysis of five genes in Epicauta chinensis larvae from the first to fifth instar. PLoS ONE 9,e89607.

[28] LiuSP,PanZ, RenGD (2016). Identification of three morphologically indistinguishable Epicauta species (Coleoptera, Meloidae, Epicautini) through DNA barcodes and morphological comparisons. Zootaxa 4103,361–373.

[29] Marçais G, KingsfordC (2011). A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 27,764–770.

[30] MendesFK,Vanderpool D,FultonB,HahnMW (2020). CAFE 5 models variation in evolutionary rates among gene families. Bioinformatics 36,5516–5518.

[31] MolfiniM,MuzziM, ManciniE et al. (2023). The cranial apparatus glands of the canthariphilou Pyrochroa coccinea (Coleoptera: Pyrochroidae: Pyrochroinae), and their implications in sexual behaviour. Arthropod Structure & Development 77,101316.

[32] MuzziM,ManciniE, FratiniE et al. (2022). Male accessory glands of blister beetles and cantharidin release: A comparative ultrastructural analysis. Insects 13,132.

[33] PanZ,BolognaMA (2014). Taxonomy, bionomics and faunistics of the nominate subgenus of Mylabris Fabricius, 1775, with the description of five new species (Coleoptera: Meloidae: Mylabrini). Zootaxa 3806,1–78.

[34] PintoJD,BolognaMA, BousemanJK (1996). First-instar larvae, courtship and oviposition in Eletica: Amending the definition of the Meloidae (Coleoptera: Tenebrionoidea). Systematic Entomology 21,63–74.

[35] Ranallo-BenavidezTR,Jaron KS,SchatzMC (2020). GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nature Communications 11,1432.

[36] RiccieriA,CapognaE, PintoJD,Bologna MA (2023). Molecular phylogeny, systematics and biogeography of the subfamily Nemognathinae (Coleoptera, Meloidae). Invertebrate Systematics 37,101–116.

[37] RiccieriA,ManciniE, PitzalisM,Salvi D,BolognaMA (2022). Multigene phylogeny of blister beetles (Coleoptera, Meloidae) reveals extensive polyphyly of the tribe Lyttini and allows redefining its boundaries. Systematic Entomology 47,569–580.

[38] RiccieriA,ManciniE, SalviD et al. (2020). Phylogeny, biogeography and systematics of the hyper-diverse blister beetle genus Hycleus (Coleoptera: Meloidae). Molecular Phylogenetics and Evolution 144,106706.

[39] ScottJG,LiuN, WenZ (1998). Insect cytochromes P450: Diversity, insecticide resistance and tolerance to plant toxins. Comparative Biochemistry and Physiology Part C: Pharmacology, Toxicology and Endocrinology 121,147–155.

[40] ScottJG,WenZ (2001). Cytochromes P450 of insects: The tip of the iceberg. Pest Management Science 57,958–967.

[41] SimonsenM,MailundT, PedersenCN (2011). Inference of large phylogenies using neighbour-joining. In: Fred A,Felipe J,Gamboa H, eds. Biomedical Engineering Systems and Technologies. BIOSTEC 2010; 20–23 Jan 2010, Valencia, Spain. Communications in Computer and Information Science, Vol. 127. Springer,Berlin Heidelberg, pp. 334–344.

[42] SongY,GuF, LiuZ et al. (2022). The key role of fatty acid synthase in lipid metabolism and metamorphic development in a destructive insect pest. Spodoptera litura (Lepidoptera: Noctuidae). International Journal of Molecular Sciences 23,9064.

[43] SpellmonN,HolcombJ, TrescottL et al. (2015). Structure and function of SET and MYND domain-containing proteins. International Journal of Molecular Sciences 16,1406–1428.

[44] SuyamaM,Torrents D,BorkP (2006). PAL2NAL: Robust conversion of protein sequence alignments into the corresponding codon alignment. Nucleic Acids Research 34,609–612.

[45] TianX,SuX, LiC et al. (2021). Draft genome of the blister beetle,Epicauta chinensis. International Journal of Biological Macromolecules 193,1694–1706.

[46] UniProt Consortium (2019). UniProt: A worldwide hub of protein knowledge. Nucleic Acids Research 8,506–515.

[47] WoodDE,LuJ, LangmeadB (2019). Improved metagenomic analysis with Kraken 2. Genome Biology 20,257.

[48] WuYM,LiJ, ChenXS (2018). Draft genomes of two blister beetles Hycleus cichorii and Hycleus phaleratus. Gigascience 7,giy006.

[49] WuYM,LiJR, LiJ et al. (2023). Investigation of sex expression profiles and the cantharidin biosynthesis genes in two blister beetles. PLoS ONE 18,e0290245.

[50] WuYM,LiuYY, ChenXS (2020). Genomic content of chemosensory receptors in two sister blister beetles facilitates characterization of chemosensory evolution. BMC Genomics 21,589.

[51] YangZ (2007). PAML 4: Phylogenetic analysis by maximum likelihood. Molecular Biology and Evolution 24,1586–1591.

[52] YokoiT,NabeT, IshizukaC et al. (2020). Transcription-inducing activity of natural and synthetic juvenile hormone agonists through the Drosophila Methoprene-tolerant protein. Pest Management Science 76,2316–2323.

[53] ZhouC,ZhengX, WangL et al. (2023). The first chromosome-level genome assembly and transcriptome sequencing provide insights into cantharidin production of the blister beetles. Integrative Zoology,

RIGHTS & PERMISSIONS

2024 The Authors. Integrative Zoology published by International Society of Zoological Sciences, Institute of Zoology/Chinese Academy of Sciences and John Wiley & Sons Australia, Ltd.