Genome of tropical bed bug Cimex hemipterus (Cimicidae, Hemiptera) reveals tetraspanin expanded in bed bug ancestor

Abstract

Cimex species are ectoparasites that exclusively feed on warm-blooded animals such as birds and mammals. Three cimicid species are known to be persistent pests for humans, including the tropical bed bug Cimex hemipterus, common bed bug Cimex lectularius, and Eastern bat bug Leptocimex boueti. To date, genomic information is restricted to the common bed bug C. lectularius, which limits understanding their biology and to provide controls of bed bug infestations. Here, a chromosomal-level genome assembly of C. hemipterus (495 Mb [megabase pairs]) contained on 16 pseudochromosomes (scaffold N50 = 34 Mb), together with 9 messenger RNA and small RNA transcriptomes were obtained. In comparison between hemipteran genomes, we found that the tetraspanin superfamily was expanded in the Cimex ancestor. This study provides the first genome assembly for the tropical bed bug C. hemipterus, and offers an unprecedented opportunity to address questions relating to bed bug infestations, as well as genomic evolution to hemipterans more widely.

Keywords

bedbug / genome / microRNA / sesquiterpenoid / tetraspanin

References

[1] Andrews, S. (2010) FastQC: a quality control tool for high throughput sequence data. Available from: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Accessed 28 January 2023.

[2] Baker, F.C.,Mauchamp, B.,Tsai, L.W. and Schooley, D.A. (1983) Farnesol and farnesal dehydrogenase (s) in corpora allata of the tobacco hornworm moth,Manduca sexta. Journal of Lipid Research,24,1586–1594.

[3] Baril, T. and Hayward, A. (2022) Migrators within migrators: exploring transposable element dynamics in the monarch butterfly,Danaus plexippus. Mobile DNA,13,5.

[4] Baril, T.,Galbraith, J. and Hayward, A. (2024) Earl Grey: a fully automated user-friendly transposable element annotation and analysis pipeline. Molecular Biology and Evolution,41,msae068.

[5] Benoit, J.B.,Adelman, Z.N.,Reinhardt, K.,Dolan, A.,Poelchau, M.,Jennings, E.C. et al. (2016) Unique features of a global human ectoparasite identified through sequencing of the bed bug genome. Nature Communications,7,10165.

[6] Bolger, A.M.,Lohse, M. and Usadel, B. (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics,30,2114–2120.

[7] Chang, Z.X.,Tang, N.,Wang, L.,Zhang, L.Q.,Akinyemi, I.A. and Wu, Q.F. (2016) Identification and characterization of microRNAs in the white-backed planthopper,Sogatella furcifera. Insect Science,23,452–468.

[8] Cheong, S.P.,Huang, J.,Bendena, W.G.,Tobe, S.S. and Hui, J.H. (2015). Evolution of ecdysis and metamorphosis in arthropods: the rise of regulation of juvenile hormone. Integrative and Comparative Biology,55,878–890.

[9] Chen, C.,Chen, H.,Zhang, Y.,Thomas, H.R.,Frank, M.H.,He, Y. et al. (2020) Tbtools: an integrative toolkit developed for interactive analyses of big biological data. Molecular Plant,13,1194–1202.

[10] Eddy, S.R. (2011) Accelerated profile HMM searches. PLoS Computational Biology,7,e1002195.

[11] Emms, D.M. and Kelly, S. (2019) OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biology,20,238.

[12] Fromm, B.,Høye, E.,Domanska, D.,Zhong, X.,Aparicio-Puerta, E.,Ovchinnikov, V. et al. (2022) MirGeneDB 2.1: toward a complete sampling of all major animal phyla. Nucleic Acids Research,50,D204–D210.

[13] Goddard, J. and deShazo, R. (2009) Bed bugs (Cimex lectularius) and clinical consequences of their bites. The Journal of the American Medical Association,301,1358–1366.

[14] Hackl, T.,Duponchel, S.,Barenhoff, K.,Weinmann, A. and Fischer, M.G. (2021) Virophages and retrotransposons colonize the genomes of a heterotrophic flagellate. eLife,10,e72674.

[15] Hallgren, J.,Tsirigos, K.D.,Pedersen, M.D.,Almagro Armenteros, J.J.,Marcatili, P.,Nielsen, H. et al. (2022) DeepTMHMM predicts alpha and beta transmembrane proteins using deep neural networks. BioRxiv. Available from:

[16] Huang, S.,Yuan, S.,Dong, M.,Su, J.,Yu, C.,Shen, Y. et al. (2005) The phylogenetic analysis of tetraspanins projects the evolution of cell–cell interactions from unicellular to multicellular organisms. Genomics,86,674–684.

[17] Huang, Y.,Sun, M.,Zhuang, L. and He, J. (2021) Molecular phylogenetic analysis of the AIG family in vertebrates. Genes,12,1190.

[18] Hubley, R.,Finn, R.D.,Clements, J.,Eddy, S.R.,Jones, T.A.,Bao, W. et al. (2016) The Dfam database of repetitive DNA families. Nucleic Acids Research,44,D81–D89.

[19] Ikeda, K.T.,Hirose, Y.,Hiraoka, K.,Noro, E.,Fujishima, K.,Tomita, M. et al. (2015) Identification, expression, and molecular evolution of microRNAs in the “living fossil” Triops cancriformis (tadpole shrimp). RNA,21,230–242.

[20] Katoh, K. and Standley, D.M. (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Molecular Biology and Evolution,30,772–780.

[21] Kim, D.,Paggi, J.M.,Park, C.,Bennett, C. and Salzberg, S.L. (2019) Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nature Biotechnology,37,907–915.

[22] Kotaki, T.,Shinada, T.,Kaihara, K.,Ohfune, Y. and Numata, H. (2009) Structure determination of a new juvenile hormone from a heteropteran insect. Organic Letters,11,5234–5237.

[23] Kozomara, A.,Birgaoanu, M. and Griffiths-Jones, S. (2019) miRbase: from microRNA sequences to function. Nucleic Acids Research,47,D155–D162.

[24] Krinsky, W.L. (2019) True bugs (Hemiptera). In: Medical and Veterinary Entomology (eds. G.R. Mullen & L.A. Durden), pp. 107–127. Cambridge, MA: Academic Press.

[25] Kumar, S.,Stecher, G. and Tamura, K. (2016) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molecular Biology and Evolution,33,1870–1874.

[26] Letunic, I. and Bork, P. (2021) Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Research,49,W293–W296.

[27] Li, Y.,Nong, W.,Baril, T.,Yip, H.Y.,Swale, T.,Hayward, A. et al. (2020) Reconstruction of ancient homeobox gene linkages inferred from a new high-quality assembly of the Hong Kong oyster (Magallana hongkongensis) genome. BMC Genomics,21,713.

[28] Lieberman-Aiden, E.,Van Berkum, N.L.,Williams, L.,Imakaev, M.,Ragoczy, T.,Telling, A. et al. (2009) Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science,326,289–293.

[29] Mackowiak, S.D. (2011) Identification of novel and known miRNAs in deep-sequencing data with miRDeep2. Current Protocols in Bioinformatics,36,12.10.1–12.10.15.

[30] Marco, A.,Ninova, M.,Ronshaugen, M. and Griffiths-Jones, S. (2013) Clusters of microRNAs emerge by new hairpins in existing transcripts. Nucleic Acids Research,41,7745–7752.

[31] Mayoral, J.G.,Nouzova, M.,Navare, A. and Noriega, F.G. (2009) NADP⁺-dependent farnesol dehydrogenase, a corpora allata enzyme involved in juvenile hormone synthesis. Proceedings of the National Academy of Sciences USA,106,21091–21096.

[32] Mei, Y.,Jing, D.,Tang, S.,Chen, X.,Chen, H.,Duanmu, H. et al. (2022) InsectBase 2.0: a comprehensive gene resource for insects. Nucleic Acids Research,50,D1040–D1045.

[33] Mendes, F.K.,Vanderpool, D.,Fulton, B. and Hahn, M.W. (2020) CAFE 5 models variation in evolutionary rates among gene families. Bioinformatics,36,5516–5518.

[34] Namiki, T.,Niwa, R.,Sakudoh, T.,Shirai, K.I.,Takeuchi, H. and Kataoka, H. (2005) Cytochrome P450 CYP307A1/Spook: a regulator for ecdysone synthesis in insects. Biochemical and Biophysical Research Communications,337,367–374.

[35] Niwa, R.,Namiki, T.,Ito, K.,Shimada-Niwa, Y.,Kiuchi, M.,Kawaoka, S. et al. (2010) Non-molting glossy/shroud encodes a short-chain dehydrogenase/reductase that functions in the ’Black Box’ of the ecdysteroid biosynthesis pathway. Development,137,1991–1999.

[36] Niwa, R. and Niwa, Y.S. (2014) Enzymes for ecdysteroid biosynthesis: their biological functions in insects and beyond. Bioscience, Biotechnology, and Biochemistry,78,1283–1292.

[37] Nong, W.,Cao, J.,Li, Y.,Qu, Z.,Sun, J.,Swale, T. et al. (2020) Jellyfish genomes reveal distinct homeobox gene clusters and conservation of small RNA processing. Nature Communications,11,3051.

[38] Nyati, P.,Nouzova, M.,Rivera-Perez, C.,Clifton, M.E.,Mayoral, J.G. and Noriega, F.G. (2013) Farnesyl phosphatase, a corpora allata enzyme involved in juvenile hormone biosynthesis in Aedes aegypti. PLoS ONE,8,e71967.

[39] Ono, H.,Rewitz, K.F.,Shinoda, T.,Itoyama, K.,Petryk, A.,Rybczynski, R. et al. (2006) Spook and Spookier code for stage-specific components of the ecdysone biosynthetic pathway in Diptera. Developmental Biology,298,555–570.

[40] Ou, Q.,Magico, A. and King-Jones, K. (2011) Nuclear receptor DHR4 controls the timing of steroid hormone pulses during Drosophila development. PLoS Biology,9,e1001160.

[41] Pertea, M.,Pertea, G.M.,Antonescu, C.M.,Chang, T.C.,Mendell, J.T. and Salzberg, S.L. (2015) StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nature Biotechnology,33,290–295.

[42] Peters, B.A.,Liu, J. and Drmanac, R. (2015) Co-barcoded sequence reads from long DNA fragments: a cost-effective solution for “perfect genome” sequencing. Frontiers in Genetics,5,466.

[43] Platt, R.N.,Blanco-Berdugo, L. and Ray, D.A. (2016) Accurate transposable element annotation is vital when analyzing new genome assemblies. Genome Biology and Evolution,8,403–410.

[44] Price, M.N.,Dehal, P.S. and Arkin, A.P. (2010) FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS ONE,5,e9490.

[45] Putnam, N.H.,O’Connell, B.L.,Stites, J.C.,Rice, B.J.,Blanchette, M.,Calef, R. et al. (2016) Chromosome-scale shotgun assembly using an in vitro method for long-range linkage. Genome Research,26,342–350.

[46] Qu, Z.,Nong, W.,So, W.L.,Barton-Owen, T.,Li, Y.,Leung, T.C.N. et al. (2020) Millipede genomes reveal unique adaptations during myriapod evolution. PLoS Biology,18,e3000636.

[47] Qu, Z.,Bendena, W.G.,Tobe, S.S. and Hui, J.H. (2018) Juvenile hormone and sesquiterpenoids in arthropods: biosynthesis, signaling, and role of microRNA. The Journal of Steroid Biochemistry and Molecular Biology,184,69–76.

[48] Reinhardt, K. and Siva-Jothy, M.T. (2007) Biology of the bed bugs (Cimicidae). Annual Review of Entomology,52,351–374.

[49] Rivera-Perez, C.,Nouzova, M.,Clifton, M.E.,Garcia, E.M.,LeBlanc, E. and Noriega, F.G. (2013) Aldehyde dehydrogenase 3 converts farnesal into farnesoic acid in the corpora allata of mosquitoes. Insect Biochemistry and Molecular Biology,43,675–682.

[50] Romero, A.,Potter, M.F.,Potter, D.A. and Haynes, K.F. (2007) Insecticide resistance in the bed bug: a factor in the pest’s sudden resurgence? Journal of Medical Entomology,44,175–178.

[51] Rosenfeld, J.A.,Reeves, D.,Brugler, M.R.,Narechania, A.,Simon, S.,Durrett, R. et al. (2016) Genome assembly and geospatial phylogenomics of the bed bug Cimex lectularius. Nature Communications,7,10164.

[52] Sadílek, D.,Urfus, T. and Vilímová J. (2019) Genome size and sex chromosome variability of bed bugs feeding on animal hosts compared to Cimex lectularius parasitizing human (Heteroptera: Cimicidae). Cytometry Part A,95,1158–1166.

[53] Sanderson, M.J. (2003) R8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics,19,301–302.

[54] Satyaveanthan, M.V.,Suhaimi, S.A.,Ng, C.L.,Muhd-Noor, N.D.,Awang, A.,Lam, K.W. et al. (2021) Purification, biochemical characterisation and bioinformatic analysis of recombinant farnesol dehydrogenase from Theobroma cacao. Plant Physiology and Biochemistry,161,143–155.

[55] Schrempf, D.,Minh, B.Q.,De Maio, N.,von Haeseler, A. and Kosiol, C. (2016) Reversible polymorphism-aware phylogenetic models and their application to tree inference. Journal of Theoretical Biology,407,362–370.

[56] Sievers, F.,Wilm, A.,Dineen, D.,Gibson, T.J.,Karplus, K.,Li, W. et al. (2011) Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Molecular Systems Biology,7,539.

[57] Singh, N.,Wang, C. and Cooper, R. (2015) Role of vision and mechanoreception in bed bug,Cimex lectularius L. behavior. PLoS ONE,10,e0118855.

[58] So, W.L.,Kai, Z.,Qu, Z.,Bendena, W.G. and Hui, J.H. (2022) Rethinking sesquiterpenoids: a widespread hormone in animals. International Journal of Molecular Sciences,23,5998.

[59] Soderlund, C.,Bomhoff, M. and Nelson, W.M. (2011) SyMAP v3. 4: a turnkey synteny system with application to plant genomes. Nucleic Acids Research,39,e68.

[60] Štefka, J.,Votýpka, J.,Lukeš J. and Balvín, O. (2022) Cimex lectularius and Cimex hemipterus (bed bugs). Trends in Parasitology,38,919–920.

[61] Subramanian, B.,Gao, S.,Lercher, M.J.,Hu, S. and Chen, W.H. (2019) Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees. Nucleic Acids Research,47,W270–W275.

[62] Termini, C.M. and Gillette, J.M. (2017) Tetraspanins function as regulators of cellular signaling. Frontiers in Cell and Developmental Biology,5,34.

[63] Thomas, I.,Kihiczak, G.G. and Schwartz, R.A. (2004) bed bug bites: a review. International Journal of Dermatology,43,430–433.

[64] Tobe, S.S. and Bendena, W.G. (1999) The regulation of juvenile hormone production in arthropods: functional and evolutionary perspectives. Annals of the New York Academy of Sciences,897,300–310.

[65] Tsang, S.S.,Law, S.T.S.,Li, C.,Qu, Z.,Bendena, W.G.,Tobe, S.S. et al. (2020) Diversity of insect sesquiterpenoid regulation. Frontiers in Genetics,11,1027.

[66] Villalobos-Sambucaro, M.J.,Nouzova, M.,Ramirez, C.E.,Eugenia Alzugaray, M.,Fernandez-Lima, F.,Ronderos, J.R. et al. (2020) The juvenile hormone described in Rhodnius prolixus by Wigglesworth is juvenile hormone III skipped bisepoxide. Scientific Reports,10,3091.

[67] Vurture, G.W.,Sedlazeck, F.J.,Nattestad, M.,Underwood, C.J.,Fang, H.,Gurtowski, J. et al. (2017) GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics,33,2202–2204.

[68] Wang, O.,Chin, R.,Cheng, X.,Wu, M.K.Y.,Mao, Q.,Tang, J. et al. (2019) Efficient and unique cobarcoding of second-generation sequencing reads from long DNA molecules enabling cost-effective and accurate sequencing, haplotyping, and de novo assembly. Genome Research,29,798–808.

[69] Wang, Y.,Tang, H.,DeBarry, J.D.,Tan, X.,Li, J.,Wang, X. et al. (2012) MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Research,40,e49.

[70] Warren, J.T.,O’Connor, M.B. and Gilbert, L.I. (2009) Studies on the Black Box: incorporation of 3-oxo-7-dehydrocholesterol into ecdysteroids by Drosophila melanogaster and Manduca sexta. Insect Biochemistry and Molecular Biology,39,677–687.

[71] Weisenfeld, N.I.,Kumar, V.,Shah, P.,Church, D.M. and Jaffe, D.B. (2017) Direct determination of diploid genome sequences. Genome Research,27,757–767.

[72] Wood, D.E.,Lu, J. and Langmead, B. (2019) Improved metagenomic analysis with Kraken 2. Genome Biology,20,257.

[73] Xu, H.,Lee, S.J.,Suzuki, E.,Dugan, K.D.,Stoddard, A.,Li, H.S. et al. (2004) A lysosomal tetraspanin associated with retinal degeneration identified via a genome-wide screen. The EMBO Journal,23,811–822.

[74] Yang, M.,Derbyshire, M.K.,Yamashita, R.A. and Marchler-Bauer, A. (2020) NCBI’s conserved domain database and tools for protein domain analysis. Current Protocols in Bioinformatics,69,e90.

[75] Zhang, F.L. and Casey, P.J. (1996) Protein prenylation: molecular mechanisms and functional consequences. Annual Review of Biochemistry,65,241–269.

[76] Zhong, Y.F. and Holland, P.W. (2011) HomeoDB2: functional expansion of a comparative homeobox gene database for evolutionary developmental biology. Evolution & Development,13,567–568.

RIGHTS & PERMISSIONS

2024 The Authors. Insect Science published by John Wiley & Sons Australia, Ltd on behalf of Institute of Zoology, Chinese Academy of Sciences.