Molecular phylogeny of the Anopheles hyrcanus group (Diptera: Culicidae) based on rDNA-ITS2 and mtDNA-COII

Analysis of geometric morphometrics and molecular phylogeny for Anopheles species in the Republic of Korea

Human malaria, transmitted by Anopheles mosquitoes, is the most predominant mosquito-borne disease that is responsible for hundreds of thousands of deaths worldwide each year. In the Republic of Korea (ROK), there are currently several hundred malaria cases annually, mostly near the demilitarized zone (DMZ). Eight species of Anopheles mosquitoes are currently known to be present in the ROK. Similar to other major malaria vectors in Africa and India, it is very challenging to morphologically differentiate Anopheles mosquitoes in the ROK due to their extremely similar morphology. In this study, wing geometric morphometrics (WGM) were used to differentiate the eight Anopheles species collected at six locations near the DMZ, Seoul and Pyeongtaek from April-October 2021. Phylogenetic analysis was also performed using cytochrome c oxidase subunit 1 (COI), internal transcribed spacer 2 (ITS2), and tyrosine hydroxylase (TH) genes for comparison with WGM analysis and to infer evolutionary relationships. The results of cross-validation (overall accuracy = 74.8%) demonstrated that species identification using WGM alone was not possible with a high accuracy for all eight species. While phylogenetic analyses based on the COI region could not clearly distinguish some species, the analysis based on ITS2 and TH was more useful for resolving the phylogenetic correlation of the eight species. Our results may improve Anopheles species identification strategies for effective identification and control of malaria vectors in the ROK.

Eimeria zuernii (Eimeriidae: Coccidia): mitochondrial genome and genetic diversity in the Chinese yak

BACKGROUND

Coccidiosis caused by Eimeria zuernii (Eimeriidae: Coccidia) represents a significant economic threat to the bovine industry. Understanding the evolutionary and genetic biology of E. zuernii can assist in new interaction developments for the prevention and control of this protozoosis.

METHODS

We defined the evolutionary and genetic characteristics of E. zuernii by sequencing the complete mitogenome and analyzing the genetic diversity and population structure of 51 isolates collected from eight yak breeding parks in China.

RESULTS

The 6176-bp mitogenome of E. zuernii was linear and encoded typical mitochondrial contents of apicomplexan parasites, including three protein-coding genes [PCGs; cytochrome c oxidase subunits I and III (cox1 and cox3), and cytochrome b (cytb)], seven fragmented small subunit (SSU) and 12 fragmented large subunit (LSU) rRNAs. Genome-wide comparative and evolutionary analyses showed cytb and cox3 to be the most and least conserved Eimeria PCGs, respectively, and placed E. zuernii more closely related to Eimeria mephitidis than other Eimeria species. Furthermore, cox1-based genetic structure defined 24 haplotypes of E. zuernii with high haplotype diversities and low nucleotide diversities across eight geographic populations, supporting a low genetic structure and rapid evolutionary rate as well as a previous expansion event among E. zuernii populations.

CONCLUSIONS

To our knowledge, this is the first study presenting the phylogeny, genetic diversity, and population structure of the yak E. zuernii, and such information, together with its mitogenomic data, should contribute to a better understanding of the genetic and evolutionary biological studies of apicomplexan parasites in bovines.

Phylogenetic Analysis of Some Species of the Anopheles hyrcanus Group (Diptera: Culicidae) in China Based on Complete Mitochondrial Genomes

Some species of the Hyrcanus group are vectors of malaria in China. However, the member species are difficult to identify accurately by morphology. The development of sequencing technologies offers the possibility of further studies based on the complete mitochondrial genome. In this study, samples of mosquitoes of the Hyrcanus group were collected in China between 1997 and 2015. The mitochondrial genomes of ten species of the Hyrcanus group were analyzed, including the structure and base composition, codon usage, secondary structure of tRNA, and base difference sites in protein coding regions. Phylogenetic analyses using maximum-likelihood and Bayesian inference were performed based on mitochondrial genes and complete mitochondrial genomes The mitochondrial genome of 10 Hyrcanus group members ranged from 15,403 bp to 15,475 bp, with an average 78.23% (A + T) content, comprising of 13 PCGs (protein coding genes), 22 tRNAs, and 2 rRNAs. Site differences between some closely related species in the PCGs were small. There were only 36 variable sites between Anopheles sinensis and Anopheles belenrae for a variation ratio of 0.32% in all PCGs. The pairwise interspecies distance based on 13 PCGs was low, with an average of 0.04. A phylogenetic tree constructed with the 13 PCGs was consistent with the known evolutionary relationships. Some phylogenetic trees constructed by single coding regions (such as COI or ND4) or combined coding regions (COI + ND2 + ND4 + ND5 or ND2 + ND4) were consistent with the phylogenetic tree constructed using the 13 PCGs. The phylogenetic trees constructed using some coding genes (COII, ND5, tRNAs, 12S rRNA, and 16S rRNA) differed from the phylogenetic tree constructed using PCGs. The difference in mitochondrial genome sequences between An. sinensis and An. belenrae was very small, corresponding to intraspecies difference, suggesting that the species was in the process of differentiation. The combination of all 13 PCG sequences was demonstrated to be optimal for phylogenetic analysis in closely related species.

Population genetic structure and evolutionary genetics of Anopheles sinensis based on knockdown resistance (kdr) mutations and mtDNA-COII gene in China-Laos, Thailand-Laos, and Cambodia-Laos borders

BACKGROUND

Vector control is still a pivotal method for preventing malaria, and its potency is weakened by the increasing resistance of vectors to chemical insecticides. As the most abundant and vital malaria vector in Southeast Asia, the chemical insecticide resistance status in Anopheles sinensis remains elusive in Laos, which makes it imperative to evaluate the true nature of chemical insecticide resistance-associated genetic mutations in An. sinensis in Laos.

METHODS

Adult An. sinensis were collected from three border regions in Laos. DNA was extracted from individual mosquitoes. PCR amplification and DNA sequencing of a fragment containing codon 1014 of the voltage-gated sodium channel (vgsc) gene were completed to study the kdr allele frequency distribution, kdr intron polymorphism, population genetic diversity, and the evolutionary status of the kdr codon. The mitochondrial cytochrome c oxidase subunit II gene (COII) was amplified and sequenced to examine population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

Nine wild kdr haplotypes of the vgsc gene were detected in this study, and eight of them, namely 1014L1, 1014L2, 1014L4, 1014L7, 1014L9, 1014L10, 1014L11, and 1014L21, were discovered in the China-Laos border (northern Laos), while 1014L3 was only detected in the Thailand-Laos border (northwestern Laos) and Cambodia-Laos border (southern Laos). The newly identified haplotype, 1014L21, was uniquely distributed in the China-Laos border and was not identified in other countries. Based on sequence analysis of the mitochondrial COII genes, significant genetic differentiation and limited gene flow were detected between the China-Laos and Cambodia-Laos An. sinensis populations, which suggested that those two regions were genetically isolated. The distinct distribution of the kdr haplotype frequencies is probably the result of geographical isolation in mosquito populations.

CONCLUSIONS

Lack of kdr mutations in the vgsc gene was probably due to genetic isolation and the absence of intense selection pressure in the three border regions of Laos. This study reveals that pyrethroid-based chemical insecticides are still appropriate for battling An. sinensis in parts of Laos, and routine monitoring of chemical insecticide resistance should be continuously implemented and focused on more restricted areas as part of chemical insecticide resistance management.

Morphological and molecular identification reveals a high diversity of Anopheles species in the forest region of the Cambodia-Laos border

Background

To develop an effective malaria vector intervention method in forested international border regions within the Greater Mekong Subregion (GMS), more in-depth studies should be conducted on local Anopheles species composition and bionomic features. There is a paucity of comprehensive surveys of biodiversity integrating morphological and molecular species identification conducted within the border of Laos and Cambodia.

Methods

A total of 2394 adult mosquitoes were trapped in the Cambodia–Laos border region. We first performed morphological identification of Anopheles mosquitoes and subsequently performed molecular identification using 412 recombinant DNA–internal transcribed spacer 2 (rDNA-ITS2) and 391 mitochondrial DNA–cytochrome c oxidase subunit 2 (mtDNA-COII) sequences. The molecular and morphological identification results were compared, and phylogenetic analysis of rDNA-ITS2 and mtDNA-COII was conducted for the sequence divergence among species.

Results

Thirteen distinct species of Anopheles were molecularly identified in a 26,415 km² border region in Siem Pang (Cambodia) and Pathoomphone (Laos). According to the comparisons of morphological and molecular identity, the interpretation of local species composition for dominant species in the Cambodia–Laos border (An. dirus, An. maculatus, An. philippinensis, An. kochi and An. sinensis) achieved the highest accuracy of morphological identification, from 98.37 to 100%. In contrast, the other species which were molecularly identified were less frequently identified correctly (0–58.3%) by morphological methods. The average rDNA-ITS2 and mtDNA-COII interspecific divergence was respectively 318 times and 15 times higher than their average intraspecific divergence. The barcoding gap ranged from 0.042 to 0.193 for rDNA-ITS2, and from 0.033 to 0.047 for mtDNA-COII.

Conclusions

The Cambodia–Laos border hosts a high diversity of Anopheles species. The morphological identification of Anopheles species provides higher accuracy for dominant species than for other species. Molecular methods combined with morphological analysis to determine species composition, population dynamics and bionomic characteristics can facilitate a better understanding of the factors driving malaria transmission and the effects of interventions, and can aid in achieving the goal of eliminating malaria.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05167-0.