Genomes of the cosmopolitan fruit pest Bactrocera dorsalis (Diptera: Tephritidae) reveal its global invasion history and thermal adaptation
Yue Zhang,Shanlin Liu,Marc De Meyer,Zuxing Liao,Yan Zhao,Massimiliano Virgilio,Shiqian Feng,Yujia Qin,Sandeep Singh,Suk Ling Wee,Fan Jiang,Shaokun Guo,Hu Li,Pablo Deschepper,Sam Vanbergen,Hélène Delatte,Alies van Sauers-Muller,Tati Suryati Syamsudin,Anastasia Priscilla Kawi,Muo Kasina,Kemo Badji,Fazal Said,Lijun Liu,Zihua Zhao,Zhihong Li
DOI: https://doi.org/10.1016/j.jare.2022.12.012
IF: 12.822
2022-12-25
Journal of Advanced Research
Abstract:Highlights • The present study achieves a large-scale SNP data set, a total of 512 individual genomes comprising 487 B. dorsalis and 25 B. carambolae. • This study shows B. dorsalis originates from the Southern India region with three independent invasion and spread routes worldwide, which are mainly facilitated by human activities including trade and immigration. • CYP6a9 is identified as candidate gene that could has enhanced the thermal adaptation of B. dorsalis and thus has boosted its invasion to temperate regions. • The gene function is further verified using the RNAi technology. Introduction The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural pests worldwide, with highly debated species delimitation, origin, and global spread routes. Objectives Our study intended to (i) resolve the taxonomic uncertainties between B. dorsalis and B. carambolae , (ii) reveal the population structure and global invasion routes of B. dorsalis across Asia, Africa, and Oceania, and (iii) identify genomic regions that are responsible for the thermal adaptation of B. dorsalis. Methods Based on a high-quality chromosome-level reference genome assembly, we explored the population relationship using a genome-scale single nucleotide polymorphism dataset generated from the resequencing data of 203 B. dorsalis genomes and 25 B. carambolae genomes Genome-wide association studies and silencing using RNA interference were used to identify and verify the candidate genes associated with extreme thermal stress. Results We showed that B. dorsalis originates from the Southern India region with three independent invasion and spread routes worldwide: (i) from Northern India to Northern Southeast Asia, then to Southern Southeast Asia; (ii) from Northern India to Northern Southeast Asian, then to China and Hawaii; and (iii) from Southern India toward the African mainland, then to Madagascar, which is mainly facilitated by human activities including trade and immigration. Twenty-seven genes were identified by a genome-wide association study to be associated with 11 temperature bioclimatic variables. The Cyp6a9 gene may enhance the thermal adaptation of B. dorsalis and thus boost its invasion, which tended to be upregulated at a hardening temperature of 38°C. Functional verification using RNA interference silencing against Cyp6a9 , led to the specific decrease in Cyp6a9 expression, reducing the survival rate of dsRNA-feeding larvae exposed to extreme thermal stress of 45°C after heat hardening treatments in B. dorsalis . Conclusion This study provides insights into the evolutionary history and genetic basis of temperature adaptation in B. dorsalis. Graphical abstract Download : Download high-res image (124KB) Download : Download full-size image
multidisciplinary sciences