METHODS: To verify the causative agent of this outbreak and characterise the viral genes, the genes encoding the structural proteins C/prM/E of viruses isolated from local residents were sequenced followed by mutation and phylogenetic analysis. Recombination, selection pressure, potential secondary structure and three-dimensional structure analyses were also performed.
RESULTS: Phylogenetic analysis revealed that all epidemic strains were of the cosmopolitan DENV-2 genotype and were most closely related to the Zhejiang strain (MH010629, 2017) and then the Malaysia strain (KJ806803, 2013). Compared with the sequence of DENV-2SS, 151 base substitutions were found in the sequences of 89 isolates; these substitutions resulted in 20 non-synonymous mutations, of which 17 mutations existed in all samples (two in the capsid protein, six in the prM/M proteins, and nine in the envelope proteins). Moreover, amino acid substitutions at the 602nd (E322:Q → H) and 670th (E390: N → S) amino acids may have enhanced the virulence of the epidemic strains. One new DNA binding site and five new protein binding sites were observed. Two polynucleotide binding sites and seven protein binding sites were lost in the epidemic strains compared with DENV-2SS. Meanwhile, five changes were found in helical regions. Minor changes were observed in helical transmembrane and disordered regions. The 429th amino acid of the E protein switched from a histamine (positively charged) to an asparagine (neutral) in all 89 isolated strains. No recombination events or positive selection pressure sites were observed. To our knowledge, this study is the first to analyse the genetic characteristics of epidemic strains in the first dengue outbreak in Hunan Province in inland China.
CONCLUSIONS: The causative agent is likely to come from Zhejiang Province, a neighbouring province where dengue fever broke out in 2017. This study may help clarify the intrinsic geographical relatedness of DENV-2 and contribute to further research on pathogenicity and vaccine development.
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.