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 487 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.
RESULTS: Mosquitoes were collected from a total of 15 sites using gravid traps and a backpack aspirator around Kampong Puruh Karu, Sarawak, Malaysian Borneo, where sylvatic DENV spillover has been documented. A total of 2447 mosquitoes comprising 10 genera and 4 species of Aedes, were collected over the three years, 2013, 2014 and 2016, in the three major land cover types in the area, homestead, agriculture and forest. Mosquitoes were identified morphologically, pooled by species and gender, homogenized, and subject to DNA barcoding of each Aedes species and to arbovirus screening. As predicted, Ae. niveus was found almost exclusively in forests whereas Ae. albopictus was collected in all land cover types. Aedes albopictus was significantly (P = 0.04) more abundant in agricultural fields than forests. Sylvatic DENV was not detected in any Aedes mosquito pools, however genomes of 14 viruses were detected using next generation sequencing.
CONCLUSIONS: Land cover type affects the abundance and distribution of the most likely bridge vectors of sylvatic DENV in Malaysia Borneo. Conversion of forests to agriculture will likely decrease the range and abundance of Ae. niveus but enhance the abundance of Ae. albopictus.