Mitochondrial DNA variation in the malaria vector Anopheles minimus across China, Thailand and Vietnam: evolutionary hypothesis, population structure and population history

Population genetic structure and wing geometric morphometrics of the filarial vector Armigeres subalbatus (Diptera: Culicidae) in Thailand

Armigeres subalbatus (Diptera: Culicidae) is a mosquito species of significant medical and veterinary importance. It is widely distributed across Southeast and East Asia and is commonly found throughout Thailand. This study assessed the genetic diversity and population structure of Ar. subalbatus in Thailand using the cytochrome c oxidase subunit I (COI) gene sequences. Additionally, wing shape variations among these populations were examined using geometric morphometrics (GM). Our results demonstrated that the overall haplotype diversity (Hd) was 0.634, and the nucleotide diversity (π) was 0.0019. Significant negative values in neutrality tests (p < 0.05) indicate that the Ar. subalbatus populations in Thailand are undergoing a phase of expansion following a bottleneck event. The mismatch distribution test suggests that the populations may have started expanding approximately 16,678 years ago. Pairwise genetic differentiation among the 12 populations based on Fst revealed significant differences in 32 pairs (p < 0.05), with the degree of differentiation ranging from 0.000 to 0.419. The GM analysis of wing shape also indicated significant differences in nearly all pairs (p < 0.05), except for between populations from Nakhon Pathom and Samut Songkhram, and between those from Chiang Mai and Mae Hong Son, suggesting no significant difference due to their similar environmental settings. These findings enhance our understanding of the population structure and phenotypic adaptations of mosquito vectors, providing vital insights for the formulation of more efficacious vector control strategies.

Insights into the mitochondrial cytochrome oxidase I (mt-COI) gene and wing morphometrics of Anopheles baimaii (Diptera: Culicidae) in malaria-endemic islands of Thailand

Anopheles baimaii (Diptera: Culicidae) significantly contributes to the transmission of parasites causing malaria in Southeast Asia and South Asia. This study examined the morphological (wing shape) and molecular (mitochondrial gene) variations of An. baimaii in four of Thailand's border islands, and also investigated the presence of Plasmodium parasites in these mosquitoes. No Plasmodium infections were detected in the samples. Significant differences in wing shape were observed in most island populations (p < 0.05). A single-linkage tree, constructed using Mahalanobis distances, clustered the populations into two groups based on geographical locations. Genetic variation in An. baimaii was also analyzed through cytochrome c oxidase subunit I (COI) gene sequences. This analysis identified 22 segregating sites and a low nucleotide diversity of 0.004. Furthermore, 18 distinct haplotypes were identified, indicating a high haplotype diversity of 0.825. Neutrality tests for the overall population revealed a significantly negative Fu's Fs value (-5.029), indicating a population expansion. In contrast, Tajima's D yielded a negative value (-1.028) that did not reach statistical significance. The mismatch distribution analysis exhibited a bimodal pattern, and the raggedness index was 0.068, showing no significant discrepancy (p = 0.485) between observed and expected distributions. Pairwise genetic differentiation assessments demonstrated significant differences between all populations (p < 0.05). These findings provide valuable insights into the COI gene and wing morphometric variations in An. baimaii across Thailand's islands, offering critical information for understanding the adaptations of this malaria vector and guiding future comprehensive research.

Revolutionizing Malaria Vector Control: The Importance of Accurate Species Identification through Enhanced Molecular Capacity

Many factors, such as the resistance to pesticides and a lack of knowledge of the morphology and molecular structure of malaria vectors, have made it more challenging to eradicate malaria in numerous malaria-endemic areas of the globe. The primary goal of this review is to discuss malaria vector control methods and the significance of identifying species in vector control initiatives. This was accomplished by reviewing methods of molecular identification of malaria vectors and genetic marker classification in relation to their use for species identification. Due to its specificity and consistency, molecular identification is preferred over morphological identification of malaria vectors. Enhanced molecular capacity for species identification will improve mosquito characterization, leading to accurate control strategies/treatment targeting specific mosquito species, and thus will contribute to malaria eradication. It is crucial for disease epidemiology and surveillance to accurately identify the Plasmodium spp. that are causing malaria in patients. The capacity for disease surveillance will be significantly increased by the development of more accurate, precise, automated, and high-throughput diagnostic techniques. In conclusion, although morphological identification is quick and achievable at a reduced cost, molecular identification is preferred for specificity and sensitivity. To achieve the targeted malaria elimination goal, proper identification of vectors using accurate techniques for effective control measures should be prioritized.

Cryptic Diversity and Demographic Expansion of Plasmodium knowlesi Malaria Vectors in Malaysia

Although Malaysia is considered free of human malaria, there has been a growing number of Plasmodium knowlesi cases. This alarming trend highlighted the need for our understanding of this parasite and its associated vectors, especially considering the role of genetic diversity in the adaptation and evolution among vectors in endemic areas, which is currently a significant knowledge gap in their fundamental biology. Thus, this study aimed to investigate the genetic diversity of Anopheles balabacensis, Anopheles cracens, Anopheles introlatus, and Anopheles latens-the vectors for P. knowlesi malaria in Malaysia. Based on cytochrome c oxidase 1 (CO1) and internal transcribed spacer 2 (ITS2) markers, the genealogic networks of An. latens showed a separation of the haplotypes between Peninsular Malaysia and Malaysia Borneo, forming two distinct clusters. Additionally, the genetic distances between these clusters were high (2.3-5.2% for CO1) and (2.3-4.7% for ITS2), indicating the likely presence of two distinct species or cryptic species within An. latens. In contrast, no distinct clusters were observed in An. cracens, An. balabacensis, or An. introlatus, implying a lack of pronounced genetic differentiation among their populations. It is worth noting that there were varying levels of polymorphism observed across the different subpopulations, highlighting some levels of genetic variation within these mosquito species. Nevertheless, further analyses revealed that all four species have undergone demographic expansion, suggesting population growth and potential range expansion for these vectors in this region.

Genetic Diversity, Haplotype Relationships, and kdr Mutation of Malaria Anopheles Vectors in the Most Plasmodium knowlesi-Endemic Area of Thailand

Plasmodium knowlesi, a malaria parasite that occurs naturally in long-tailed macaques, pig-tailed macaques, and banded leaf monkeys, is currently regarded as the fifth of the human malaria parasites. We aimed to investigate genetic diversity based on the cytochrome c oxidase subunit I (COI) gene, detect Plasmodium parasites, and screen for the voltage-gated sodium channel (VGSC)-mutation-mediated knockdown resistance (kdr) of Anopheles mosquitoes in Ranong province, which is the most P. knowlesi-endemic area in Thailand. One hundred and fourteen Anopheles females belonging to eight species, including An. baimaii (21.05%), An. minimus s.s. (20.17%), An. epiroticus (19.30%), An. jamesii (19.30%), An. maculatus s.s. (13.16%), An. barbirostris A3 (5.26%), An. sawadwongporni (0.88%), and An. aconitus (0.88%), were caught in three geographical regions of Ranong province. None of the Anopheles mosquitoes sampled in this study were infected with Plasmodium parasites. Based on the sequence analysis of COI sequences, An. epiroticus had the highest level of nucleotide diversity (0.012), followed by An. minimus (0.011). In contrast, An. maculatus (0.002) had the lowest level of nucleotide diversity. The Fu's Fs and Tajima's D values of the Anopheles species in Ranong were all negative, except the Tajima's D values of An. minimus (0.077). Screening of VGSC sequences showed no presence of the kdr mutation of Anopheles mosquitoes. Our results could be used to further select effective techniques for controlling Anopheles populations in Thailand's most P. knowlesi-endemic area.

Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border

BACKGROUND

Few studies have been conducted to investigate the distribution of mosquito vectors and the population structure of secondary vectors in the border region of Cambodia-Laos. The aim of this work was to study the mosquito diversity and molecular phylogeny of secondary vectors, i.e., Anopheles nivipes in this area.

METHODS

1440 adult mosquitoes were trapped in the Cambodia-Laos border. mtDNA-COII were amplified and sequenced from 53 An. nivipes DNA samples. Together with COII sequences deposited in GenBank, a total of 86 COII sequences were used for examining population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

The adult mosquitoes were classified into 5 genera and 27 species in this border region. The predominant genera were Culex (60.07%, 865/1440) and Anopheles (31.25%, 450/1440), and the major Anopheles species were An. nivipes (73.56%, 331/450) and Anopheles maculatus (14.22%, 64/450). Based on sequences analysis of COII, a high level of genetic differentiation was reported in two Northwest India (Cheema and Bathinda, Punjab) and Cambodia-Laos (Siem Pang, Stung treng) populations (F_ST = 0.97824, 0.97343, P < 0.05), as well as lower gene flow (Nm = 0.01112, 0.01365) in the An. nivipes populations. Phylogenetic analysis and SAMOVA revealed a gene barrier restricting gene flow among three An. nivipes populations. Mantel test suggested a significant correlation between geography and gene distance in all An. nivipes populations (Z = 44,983.1865, r = 0.5575, P = 0.0070). Neutrality test and Mismatch distribution revealed a recent population expansion of An. nivipes in the Cambodia-Laos population.

CONCLUSIONS

Anopheles nivipes was one of the major Anopheles species in the Cambodia-Laos border. Based on sequences analysis of COII, a genetic barrier between Cambodia-Laos and two Indian populations was found, and a recent population expanding or selecting of An. nivipes occurred in the Cambodia-Laos population, suggesting that COII might be an effective marker for describing the molecular phylogeny of An. nivipes. Further investigation and continuous surveillance of An. nivipes are warranted in this region.

Population genetic structure of the malaria vector Anopheles minimus in Thailand based on mitochondrial DNA markers

BACKGROUND

The malaria vector Anopheles minimus has been influenced by external stresses affecting the survival rate and vectorial capacity of the population. Since An. minimus habitats have continuously undergone ecological changes, this study aimed to determine the population genetic structure and the potential gene flow among the An. minimus populations in Thailand.

METHODS

Anopheles minimus was collected from five malaria transmission areas in Thailand using Centers for Disease Control and Prevention (CDC) light traps. Seventy-nine females from those populations were used as representative samples. The partial mitochondrial cytochrome c oxidase subunit I (COI), cytochrome c oxidase subunit II (COII) and cytochrome b (Cytb) gene sequences were amplified and analyzed to identify species and determine the current population genetic structure. For the past population, we determined the population genetic structure from the 60 deposited COII sequences in GenBank of An. minimus collected from Thailand 20 years ago.

RESULTS

The current populations of An. minimus were genetically divided into two lineages, A and B. Lineage A has high haplotype diversity under gene flow similar to the population in the past. Neutrality tests suggested population expansion of An. minimus, with the detection of abundant rare mutations in all populations, which tend to arise from negative selection.

CONCLUSIONS

This study revealed that the population genetic structure of An. minimus lineage A was similar between the past and present populations, indicating high adaptability of the species. There was substantial gene flow between the eastern and western An. minimus populations without detection of significant gene flow barriers.

Isolation by geographical distance after release from Pleistocene refugia explains genetic and phenotypic variation in Xylotrupes siamensis (Coleoptera: Scarabaeidae)

Consistent and objective species delimitation is crucial to biodiversity studies, but challenges remain when conflicting taxonomic decisions have been made because different data sets and analytical methods were used to delineate species. In the rhinoceros beetle, Xylotrupes siamensis, the use of different morphological characters has resulted in taxonomic disagreement between studies. We used three molecular loci (mitochondrial CO1 and nuclear ITS2 and H3) to investigate the genetic divergence between populations exhibiting different male horn phenotypes. We also applied an approximate Bayesian computation approach to test alternative historical hypotheses that might explain the present genetic diversity among geographical populations. Furthermore, we used species distribution models to estimate the temporal variation in the geographical distribution of suitable habitats. The results show that the two phenotypic taxa within X. siamensis are not genetically structured and that their genetic structure can be explained using isolation by geographical distance. The emergence of the two phenotypic taxa might have been associated with historical isolation in separate refugia. However, spatial expansion and genetic interchange between populations might have gradually eroded the spatial genetic structure. We demonstrate that understanding the historical processes responsible for phenotypic divergence and genetic diversity among current populations could help with making evolutionarily coherent taxonomic decisions.

Population structure and genetic diversity of Aedes aegypti and Aedes albopictus in Penang as revealed by mitochondrial DNA cytochrome oxidase I

The population genetics study is crucial as it helps in understanding the epidemiological aspects of dengue and help improving a vector control measures. This research aims to investigate the population genetics structure of two common species of Aedes mosquitoes in Penang; Aedes aegypti and Aedes albopictus using Cytochrome Oxidase I (COI) mitochondrial DNA (mtDNA) marker. Molecular investigations were derived from 440 bp and 418 bp mtDNA COI on 125 and 334 larvae of Aedes aegypti and Aedes albopictus respectively, from 32 locations in Penang. All samples were employed in the BLASTn for species identification. The haplotype diversity, nucleotide diversity, neutrality test and mismatch distribution analysis were conducted in DnaSP version 5.10.1. AMOVA analysis was conducted in ARLEQUIN version 3.5 and the phylogenetic reconstructions based on maximum likelihood (ML) and neighbor-joining (NJ) methods were implemented in MEGA X. The relationships among haplotypes were further tested by creating a minimum spanning tree using Network version 4.6.1. All samples were genetically identified and clustered into six distinct species. Among the species, Ae. albopictus was the most abundant (67.2%), followed by Ae. aegypti (25.2%) and the rest were counted for Culex sp. and Toxorhynchites sp. Both Ae. aegypti and Ae. albopictus show low nucleotide diversity (π) and high haplotype diversity (h), while the neutrality test shows a negative value in most of the population for both species. There are a total of 39 and 64 haplotypes recorded for Ae. aegypti and Ae. albopictus respectively. AMOVA analysis revealed that most of the variation occurred within population for both species. Mismatch distribution analysis showed bimodal characteristic of population differentiation for Ae. aegypti but Ae. albopictus showed unimodal characteristics of population differentiation. Genetic distance based on Tamura-Nei parameter showed low genetic divergent within population and high genetic divergent among population for both species. The maximum likelihood tree showed no obvious pattern of population genetic structure for both Ae. aegypti and Ae. albopictus from Penang and a moderate to high bootstrap values has supported this conclusion. The minimum spanning network for Ae. aegypti and Ae. albopictus showed five and three dominant haplotypes respectively, which indicates a mixture of haplotypes from the regions analysed. This study revealed that there is no population genetic structure exhibited by both Ae. aegypti and Ae. albopictus in Penang. Mutation has occurred rapidly in both species and this will be challenging in controlling the populations. However, further analysis needed to confirm this statement.

Genetic homogeneity of Anopheles maculatus in Indonesia and origin of a novel species present in Central Java

Background

Anopheles maculatus (s.s.) is an important vector of malaria in Indonesia. Previously it was considered the only member of the Maculatus Group present in Indonesia. A novel species was recently identified in the Kulon Progo District in Central Java. Until recently, few investigations have been conducted looking at An. maculatus genetic diversity in Indonesia, including allopatric island populations.

Methods

Indonesian An. maculatus (s.l.) samples were collected in several locations in Java, Lesser Sunda Island group, Sumatra and in Kulon Progo (Yogyakarta, central Java) where a novel species has been identified. Samples from a 30-year-old colony of the Kulon Progo population were also included in the analysis. Maximum-likelihood analysis established the phylogenies of the ITS2 (nuclear) and cox1 (mitochondrial) markers. Putative times of separation were based on cox1 genetic distances.

Results

Two species of the Maculatus Group are present in Indonesia. The novel sibling species is more closely related to Anopheles dispar than to An. maculatus (s.s.). Anopheles maculatus (s.s.) samples are homogeneous based on the ITS2 sequences. Indonesian samples and An. dispar belong to the same cox1 maternal lineage and differ from all other known members of the Maculatus Group. Divergence time between the different populations found in Java was estimated using an established cox1 mutation rate.

Conclusions

A novel species within the Maculatus Group, most closely related to An. dispar, is confirmed present in the Kulon Progo area of Central Java. The divergence of this species from An. maculatus (s.s.) is explained by the stable refugia in the Kulon Progo area during the quaternary period of intense volcanic activity throughout most of Java. This novel species awaits detailed morphological description before applying a formal species name. For the interim, it is proposed that the Kulon Progo population be designated An. maculatus var. menoreh to distinguish it from An. maculatus (s.s.).

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3598-1) contains supplementary material, which is available to authorized users.

Intraspecific phylogeny of Anopheles (Kerteszia) neivai Howard, Dyar & Knab 1913, based on mitochondrial and nuclear ribosomal genes

Three mitochondrial regions and a fragment of a large nuclear ribosomal subunit was used to study the evolutionary patterns of An. neivai, a mosquito inhabiting mangroves and tropical forest in the lowland and coastal areas of the Yucatan Peninsula through the Pacific Ecuadorian coast. This species exhibits epidemiological importance regarding Malaria transmission in natural ecosystems, particularly in rural areas of the Pacific Colombian coast. The results based on phylogenetic networks and Bayesian inference showed no robust evidence supporting the existence of previously suggested cryptic species. Diversification patterns in geographically widespread species such as this one, are complex and therefore could impact malaria control strategies. Further studies focused on behavior, morphology, and phylogenomics will improve the understanding of the evolutionary patterns within An. neivai and its role as a disease vector.

Mitochondrial variation in subpopulations of Anopheles balabacensis Baisas in Sabah, Malaysia (Diptera: Culicidae)

Anopheles balabacensis, the primary vector of Plasmodium knowlesi in Sabah, Malaysia, is both zoophilic and anthropophilic, feeding on macaques as well as humans. It is the dominant Anopheles species found in Kudat Division where it is responsible for all the cases of P. knowlesi. However there is a paucity of basic biological and ecological information on this vector. We investigated the genetic variation of this species using the sequences of cox1 (1,383 bp) and cox2 (685 bp) to gain an insight into the population genetics and inter-population gene flow in Sabah. A total of 71 An. balabacensis were collected from seven districts constituting 14 subpopulations. A total of 17, 10 and 25 haplotypes were detected in the subpopulations respectively using the cox1, cox2 and the combined sequence. Some of the haplotypes were common among the subpopulations due to gene flow occurring between them. AMOVA showed that the genetic variation was high within subpopulations as compared to between subpopulations. Mantel test results showed that the variation between subpopulations was not due to the geographical distance between them. Furthermore, Tajima's D and Fu's Fs tests showed that An. balabacensis in Sabah is experiencing population expansion and growth. High gene flow between the subpopulations was indicated by the low genetic distance and high gene diversity in the cox1, cox2 and the combined sequence. However the population at Lipasu Lama appeared to be isolated possibly due to its higher altitude at 873 m above sea level.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

Anopheles Vectors in Mainland China While Approaching Malaria Elimination

China is approaching malaria elimination; however, well-documented information on malaria vectors is still missing, which could hinder the development of appropriate surveillance strategies and WHO certification. This review summarizes the nationwide distribution of malaria vectors, their bionomic characteristics, control measures, and related studies. After several years of effort, the area of distribution of the principal malaria vectors was reduced, in particular for Anopheles lesteri (synonym: An. anthropophagus) and Anopheles dirus s.l., which nearly disappeared from their former endemic regions. Anopheles sinensis is becoming the predominant species in southwestern China. The bionomic characteristics of these species have changed, and resistance to insecticides was reported. There is a need to update surveillance tools and investigate the role of secondary vectors in malaria transmission.

An effective method for the identification and separation of Anopheles minimus, the primary malaria vector in Thailand, and its sister species Anopheles harrisoni, with a comparison of their mating behaviors

BACKGROUND

Species of the Anopheles minimus complex are considered to be the primary vectors of malaria in South and Southeast Asia. Two species of the complex, Anopheles minimus and Anopheles harrisoni, occur in Thailand. They are sympatric and difficult to accurately distinguish based on morphological characters. The aim of this study was to investigate the potential of antennal sensory organs to distinguish these two species. Additionally, we investigated their ability to mate in cages of different sizes, as well as the possible mechanism(s) that evokes stenogamous behavior.

METHODS

Large sensilla coeloconica present on the antennae of females of An. minimus and An. harrisoni were counted under a conventional light microscope and various types of antennal sensilla were examined under a scanning electron microscope (SEM). Determinations of mating ability were carried out in 20 and 30 cm³ cages with a density resting surface (DRS) of 7.2. The insemination rate, frequency of clasper (gonocoxopodite) movement of the male genitalia during induced copulation and duration of mating of the two species were compared.

RESULTS

The mean numbers of large sensilla coeloconica on antennal flagellomeres 1-8 and the mean number of large sensilla coeloconica on each flagellum in An. minimus (26.25) and An. harrisoni (31.98) were significantly different. Females of both species bear five types of antennal sensilla: chaetica, trichodea, basiconica, coeloconica and ampullacea. Marked differences in the structure of the large sensilla coeloconica were observed between the two species. Furthermore, only An. minimus could copulate naturally in the small cages. The frequency of clasper movement in the stenogamous An. minimus was significantly higher than in An. harrisoni, but there was no difference in the duration of mating.

CONCLUSIONS

To our knowledge, this study is the first to examine and discover the usefulness of large sensilla coeloconica on the antennae of females and the frequency of clasper movement in males for distinguishing the sibling species An. minimus and An. harrisoni. The discovery provides an effective and relatively inexpensive method for their identification. Additionally, the greater frequency of clasper movement of An. minimus might influence its ability to mate in small spaces.

The complete mitochondrial genome of Anopheles minimus (Diptera: Culicidae) and the phylogenetics of known Anopheles mitogenomes

Anopheles minimus is an important vector of human malaria in southern China and Southeast Asia. The phylogenetics of mosquitoes has not been well resolved, and the mitochondrial genome (mtgenome) has proven to be an important marker in the study of evolutionary biology. In this study, the complete mtgenome of An. minimus was sequenced for the first time. It is 15 395 bp long and encodes 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs) and a non-coding region. The gene organization is consistent with those of known Anopheles mtgenomes. The mtgenome performs a clear bias in nucleotide composition with a positive AT-skew and a negative GC-skew. All 13 PCGs prefer to use the codon UUA (Leu), ATN as initiation codon but cytochrome-oxidase subunit 1 (COI) and ND5, with TCG and GTG, and TAA as termination codon, but COI, COII, COIII and ND4, all with the incomplete T. tRNAs have the typical clover-leaf structure, but tRNA(Ser(AGN)) is consistent with known Anopheles mtgenomes. The control region includes a conserved T-stretch and a (TA)n stretch, and has the highest A+T content at 93.1%. The phylogenetics of An. minimus with 18 other Anopheles species was constructed by maximum likelihood and Bayesian inference, based on concatenated PCG sequences. The subgenera, Cellia and Anopheles, and Nyssorhynchus and Kerteszia have mutually close relationships, respectively. The Punctulatus group and Leucosphyrus group of Neomyzomyia Series, and the Albitarsis group of Albitarsis Series were suggested to be monophyletic. The monophyletic status of the subgenera, Cellia, Anopheles, Nyssorhynchus and Kerteszia need to be further elucidated.

Northern Bobwhite (Colinus virginianus) Mitochondrial Population Genomics Reveals Structure, Divergence, and Evidence for Heteroplasmy

Herein, we evaluated the concordance of population inferences and conclusions resulting from the analysis of short mitochondrial fragments (i.e., partial or complete D-Loop nucleotide sequences) versus complete mitogenome sequences for 53 bobwhites representing six ecoregions across TX and OK (USA). Median joining (MJ) haplotype networks demonstrated that analyses performed using small mitochondrial fragments were insufficient for estimating the true (i.e., complete) mitogenome haplotype structure, corresponding levels of divergence, and maternal population history of our samples. Notably, discordant demographic inferences were observed when mismatch distributions of partial (i.e., partial D-Loop) versus complete mitogenome sequences were compared, with the reduction in mitochondrial genomic information content observed to encourage spurious inferences in our samples. A probabilistic approach to variant prediction for the complete bobwhite mitogenomes revealed 344 segregating sites corresponding to 347 total mutations, including 49 putative nonsynonymous single nucleotide variants (SNVs) distributed across 12 protein coding genes. Evidence of gross heteroplasmy was observed for 13 bobwhites, with 10 of the 13 heteroplasmies involving one moderate to high frequency SNV. Haplotype network and phylogenetic analyses for the complete bobwhite mitogenome sequences revealed two divergent maternal lineages (d_XY = 0.00731; F_ST = 0.849; P < 0.05), thereby supporting the potential for two putative subspecies. However, the diverged lineage (n = 103 variants) almost exclusively involved bobwhites geographically classified as Colinus virginianus texanus, which is discordant with the expectations of previous geographic subspecies designations. Tests of adaptive evolution for functional divergence (MKT), frequency distribution tests (D, F_S) and phylogenetic analyses (RAxML) provide no evidence for positive selection or hybridization with the sympatric scaled quail (Callipepla squamata) as being explanatory factors for the two bobwhite maternal lineages observed. Instead, our analyses support the supposition that two diverged maternal lineages have survived from pre-expansion to post-expansion population(s), with the segregation of some slightly deleterious nonsynonymous mutations.

Population Genetic Structure of the Malaria Vector Anopheles sinensis (Diptera: Culicidae) Sensu Stricto and Evidence for Possible Introgression in the Republic of Korea

Anopheles sinensis Wiedemann sensu stricto (s.s.) is a dominant mosquito and considered a secondary malaria vector in the Republic of Korea (ROK). Despite the potential significance for malaria control, population genetics studies have been conducted using only mitochondrial DNA (mtDNA), and studies of the genetics of hybridization have never been attempted. In this study, 346 specimens from 23 localities were subject to experiments. Among them, 305 An. sinensis s.s. specimens from 20 localities were used for mtDNA analysis, and 346 specimens comprising 341 An. sinensis s.s. from 22 localities and five Anopheles kleini Rueda from one locality were examined in the microsatellite study. Neighbor-joining analysis of pairwise FST and RST based on microsatellite results showed that the populations are divided into two groups, as did the mtDNA results. However, the Bayesian analysis and factorial correspondence analysis plots showed three distinct clusters. Among the mtDNA and microsatellite results, only microsatellites represented small but positive and significant isolation-by-distance patterns. Both molecular markers show the Taebaek and Sobaek Mountain ranges as barriers between the northern and southern parts of the ROK. The newly recognized third group suggests possible introgressive hybridization of An. sinensis s.s. with closely related species. The slightly different composition of populations in each group based on different markers is probably because of different population dynamics in each group. These results imply that there is restricted gene flow of epidemiologically important malaria-related genes between the northern and southern parts of the ROK.

Multilocus nuclear DNA markers reveal population structure and demography of Anopheles minimus

Utilization of multiple putatively neutral DNA markers for inferring evolutionary history of species population is considered to be the most robust approach. Molecular population genetic studies have been conducted in many species of Anopheles genus, but studies based on single nucleotide polymorphism (SNP) data are still very scarce. Anopheles minimus is one of the principal malaria vectors of Southeast (SE) Asia including the Northeastern (NE) India. Although population genetic studies with mitochondrial genetic variation data have been utilized to infer phylogeography of the SE Asian populations of this species, limited information on the population structure and demography of Indian An. minimus is available. We herewith have developed multilocus nuclear genetic approach with SNP markers located in X chromosome of An. minimus in eight Indian and two SE Asian population samples (121 individual mosquitoes in total) to infer population history and test several hypotheses on the phylogeography of this species. While the Thai population sample of An. minimus presented the highest nucleotide diversity, majority of the Indian samples were also fairly diverse. In general, An. minimus populations were moderately substructured in the distribution range covering SE Asia and NE India, largely falling under three distinct genetic clusters. Moreover, demographic expansion events could be detected in the majority of the presently studied populations of An. minimus. Additional DNA sequencing of the mitochondrial COII region in a subset of the samples (40 individual mosquitoes) corroborated the existing hypothesis of Indian An. minimus falling under the earlier reported mitochondrial lineage B.

Molecular population genetics of the NADPH cytochrome P450 reductase (CPR) gene in Anopheles minimus

Development of insecticide resistance (IR) in mosquito vectors is a primary huddle to malaria control program. Since IR has genetic basis, and genes constantly evolve with response to environment for adaptation to organisms, it is important to know evolutionary pattern of genes conferring IR in malaria vectors. The mosquito Anopheles minimus is a major malaria vector of the Southeast (SE) Asia and India and is susceptible to all insecticides, and thus of interest to know if natural selection has shaped variations in the gene conferring IR. If not, the DNA fragment of such a gene could be used to infer population structure and demography of this species of malaria vector. We have therefore sequenced a ~569 bp DNA segment of the NADPH cytochrome P450 reductase (CPR) gene (widely known to confer IR) in 123 individuals of An. minimus collected in 10 different locations (eight Indian, one Thai and one Vietnamese). Two Indian population samples were completely mono-morphic in the CPR gene. In general, low genetic diversity was found with no evidence of natural selection in this gene. The data were therefore analyzed to infer population structure and demography of this species. The 10 populations could be genetically differentiated into four different groups; the samples from Thailand and Vietnam contained high nucleotide diversity. All the 10 populations conform to demographic equilibrium model with signature of past population expansion in four populations. The results in general indicate that the An. minimus mosquitoes sampled in the two SE Asian localities contain several genetic characteristics of being parts of the ancestral population.

Population explosion in the yellow-spined bamboo locust Ceracris kiangsu and inferences for the impact of human activity

Geographic distance and geographical barriers likely play a considerable role in structuring genetic variation in species, although some migratory species may have less phylogeographic structure on a smaller spatial scale. Here, genetic diversity and the phylogenetic structure among geographical populations of the yellow-spined bamboo locust, Ceracris kiangsu, were examined with 16S rDNA and amplified fragment length polymorphisms (AFLPs). In this study, no conspicuous phylogeographical structure was discovered from either Maximum parsimony (MP) and Neighbor-joining (NJ) phylogenetic analyses. The effect of geographical isolation was not conspicuous on a large spatial scale.At smaller spatial scales local diversity of some populations within mountainous areas were detected using Nei's genetic distance and AMOVA. There is a high level of genetic diversity and a low genetic differentiation among populations in the C. kiangsu of South and Southeast China. Our analyses indicate that C. kiangsu is a monophyletic group. Our results also support the hypothesis that the C. kiangsu population is in a primary differentiation stage. Given the mismatch distribution, it is likely that a population expansion in C. kiangsu occurred about 0.242 Ma during the Quaternary interglaciation. Based on historical reports, we conjecture that human activities had significant impacts on the C. kiangsu gene flow.

Population genetics of the malaria vector Anopheles aconitus in China and Southeast Asia

Anopheles aconitus is a well-known vector of malaria and is broadly distributed in the Oriental Region, yet there is no information on its population genetic characteristics. In this study, the genetic differentiation among populations was examined using 140 mtDNA COII sequences from 21 sites throughout Southern China, Myanmar, Vietnam, Thailand, Laos and Sri Lanka. The population in Sri Lanka has characteristic rDNA D3 and ITS2, mtDNA COII and ND5 haplotypes, and may be considered a distinct subspecies. Clear genetic structure was observed with highly significant genetic variation present among population groups in Southeast Asia. The greatest genetic diversity exists in Yunnan and Myanmar population groups. All population groups are significantly different from one another in pairwise Fst values, except Northern Thailand with Central Thailand. Mismatch distributions and extremely significant F(s) values suggest that the populations passed through a recent demographic expansion. These patterns are discussed in relation to the likely biogeographic history of the region and compared to other Anopheles species.

Glacial history of a modern invader: phylogeography and species distribution modelling of the Asian tiger mosquito Aedes albopictus

Background

The tiger mosquito, Aedes albopictus, is one of the 100 most invasive species in the world and a vector of human diseases. In the last 30 years, it has spread from its native range in East Asia to Africa, Europe, and the Americas. Although this modern invasion has been the focus of many studies, the history of the species’ native populations remains poorly understood. Here, we aimed to assess the role of Pleistocene climatic changes in shaping the current distribution of the species in its native range.

Methodology/Principal Findings

We investigated the phylogeography, historical demography, and species distribution of Ae. albopictus native populations at the Last Glacial Maximum (LGM). Individuals from 16 localities from East Asia were analyzed for sequence variation at two mitochondrial genes. No phylogeographic structure was observed across the study area. Demographic analyses showed a signature of population expansion that started roughly 70,000 years BP. The occurrence of a continuous and climatically suitable area comprising Southeast China, Indochinese Peninsula, and Sundaland during LGM was indicated by species distribution modelling.

Conclusions/Significance

Our results suggest an evolutionary scenario in which, during the last glacial phase, Ae. albopictus did not experience a fragmentation phase but rather persisted in interconnected populations and experienced demographic growth. The wide ecological flexibility of the species probably played a crucial role in its response to glacial-induced environmental changes. Currently, there is little information on the impact of Pleistocene climatic changes on animal species in East Asia. Most of the studies focused on forest-associated species and suggested cycles of glacial fragmentation and post-glacial expansion. The case of Ae. albopictus, which exhibits a pattern not previously observed in the study area, adds an important piece to our understanding of the Pleistocene history of East Asian biota.

The polymorphism and the geographical distribution of the knockdown resistance (kdr) of Anopheles sinensis in the Republic of Korea

Background

In the Republic of Korea (ROK), six sibling species of the Anopheles sinensis complex are considered the vector species of malaria, but data on their susceptibilities to malaria and vector capacities have been controversial. The intensive use of insecticides has contributed to the rapid development and spread of insecticide resistance in the An. sinensis complex. Knockdown resistance (kdr) to pyrethroids and DDT in the An. sinensis complex is associated with a mutation in codon 1014 of the voltage-gated sodium channel (VGSC) gene. Because the degree of insecticide resistance varies among mosquito species and populations, the detection of kdr mutations among the six sibling species of the An. sinensis complex is a prerequisite for establishing effective long-term vector control strategies in the ROK

Methods

In order to investigate species-specific kdr mutations, An. sinensis complex specimens have been collected from 22 sites in the ROK. Because of the difficulties with species identifications that are based only on morphological characteristics, molecular identification methods have been conducted on every specimen. Part of the IIS6 domain of the VGSC was polymerase chain reaction-amplified and directly sequenced.

Results

The molecular analyses revealed that mutations existed at codon 1014 only in An. sinensis sensu stricto and no mutations were found in the other five Anopheles species. In An. sinensis s.s., one wild type (TTG L1014) and three mutant types (TTT L1014F, TTC L1014F, and TGT L1014C) of kdr alleles were detected. The TTC L1014F mutation was observed for the first time in this species.

Conclusions

The fact that the highly polymorphic kdr gene is only observed in An. sinensis s.s., out of the six Anopheles species and their geographical distribution suggest the need for future studies of insecticide resistance monitoring and investigations of species-specific resistance mechanisms in order to build successful malaria vector control programmes in the ROK.

Mitochondrial genetic differentiation across populations of the malaria vector Anopheles lesteri from China (Diptera: Culicidae)

Background

Anopheles lesteri is a primary vector of Plasmodium spp. in central China. A complete understanding of vector population structure and the processes responsible for the differentiation is important to the vector-based malaria control programmes and for identifying heterogeneity in disease transmission as a result of discrete vector populations. There is no adequate An. lesteri population genetic data available.

Methods

Polymorphism of sequence variations in mitochondrial COII and Cytb genes were assessed to explore the level of genetic variability and differentiation among six populations of An. lesteri from China.

Results

There were 30 (4.37%) and 21 (5.33%) polymorphic sites for mtDNA-COII and Cytb gene, respectively. Totally 31 COII and 30 Cytb haplotypes were obtained. The range of F_ST values was from 0.101 to 0.655 by mtDNA-COII, and 0.029 to 0.231 by Cytb gene. The analysis of molecular variance (AMOVA) showed that the percentage of variation within populations (65.83%, 88.48%) was greater than that among populations (34.17%, 11.52%) using both genes. The Tajima's D and Fu's Fs values were all negative, except Tajima's D values of YN and HNB populations, which suggest a large number of low-frequency mutations in populations and the populations were in expansion proceeding.

Conclusions

Levels of genetic variation within An. lesteri populations were higher than among them. While these results may suggest considerable levels of gene flow, other explanations, such as the effect of historical population perturbations can also be hypothesized.

Genetic and phenotypic variation of the malaria vector Anopheles atroparvus in southern Europe

BACKGROUND

There is a growing concern that global climate change will affect the potential for pathogen transmission by insect species that are vectors of human diseases. One of these species is the former European malaria vector, Anopheles atroparvus. Levels of population differentiation of An. atroparvus from southern Europe were characterized as a first attempt to elucidate patterns of population structure of this former malaria vector. Results are discussed in light of a hypothetical situation of re-establishment of malaria transmission.

METHODS

Genetic and phenotypic variation was analysed in nine mosquito samples collected from five European countries, using eight microsatellite loci and geometric morphometrics on 21 wing landmarks.

RESULTS

Levels of genetic diversity were comparable to those reported for tropical malaria vectors. Low levels of genetic (0.004 <FST <0.086) and phenotypic differentiation were detected among An. atroparvus populations spanning over 3,000 km distance. Genetic differentiation (0.202 <FST <0.299) was higher between the sibling species An. atroparvus and Anopheles maculipennis s.s. Differentiation between sibling species was not so evident at the phenotype level.

CONCLUSIONS

Levels of population differentiation within An. atroparvus were low and not correlated with geographic distance or with putative physical barriers to gene flow (Alps and Pyrenées). While these results may suggest considerable levels of gene flow, other explanations such as the effect of historical population perturbations can also be hypothesized.