The Chromosome-Level Genome Assembly of Bean Blossom Thrips (Megalurothrips usitatus) Reveals an Expansion of Protein Digestion-Related Genes in Adaption to High-Protein Host Plants

Abstract

Megalurothrips usitatus (Bagnall) is a destructive pest of legumes, such as cowpea. The biology, population dynamics and control strategies of this pest have been well studied. However, the lack of a high-quality reference genome for M. usitatus has hindered the understanding of key biological questions, such as the mechanism of adaptation to feed preferentially on high-protein host plants and the resistance to proteinase inhibitors (PIs). In this study, we generated a high-resolution chromosome-level reference genome assembly (247.82 Mb, 16 chromosomes) of M. usitatus by combining Oxford Nanopore Technologies (ONT) and Hi-C sequencing. The genome assembly showed higher proportions of GC and repeat content compared to other Thripinae species. Genome annotation revealed 18,624 protein-coding genes, including 4613 paralogs that were preferentially located in TE-rich regions. GO and KEGG enrichment analyses of the paralogs revealed significant enrichment in digestion-related genes. Genome-wide identification uncovered 506 putative digestion-related enzymes; of those, proteases, especially their subgroup serine proteases (SPs), are significantly enriched in paralogs. We hypothesized that the diversity and expansion of the digestion-related genes, especially SPs, could be driven by mobile elements (TEs), which promote the adaptive evolution of M. usitatus to high-protein host plants with high serine protease inhibitors (SPIs). The current study provides a valuable genomic resource for understanding the genetic variation among different pest species adapting to different plant hosts.

Chromosome-level genome assembly of the western flower thrips Frankliniella occidentalis

“…The functionally annotated terms were discrepant according to the reference databases, ranged from 15619 PCGs for Nr database to 370 domains for InterPro database (Table 2). The G + C content of the final genome assembly was 50.75% (Table 1), which is similar to that of the published WFT genome 18 , lower than that of Frankliniella intonsa 30,31 , Megalurothrips usitatus 32,33 , Stenchaetothrips biformis 34 and Thrips palmi 35 , while slightly higher than those of Frankliniella fusca 36 .…”