주메뉴 바로가기 본문 바로가기
  • KDCA
  • Contact us
  • E-Submission

PHRP : Osong Public Health and Research Perspectives

OPEN ACCESS. pISSN: 2210-9099. eISSN: 2233-6052
Search
Advanced Search
About
Browse articles
Editorial policies
For authors

Articles

Page Path

Original Article

Genetic diversity and spatiotemporal population structure of Anopheles sinensis in the Republic of Korea based on the mitochondrial cytochrome c oxidase subunit I (COI) marker

Haneul Jung1orcid, Bo Gyeong Han1orcid, Hyun-Il Shin2orcid, Myoung-Ro Lee1orcid, Jung-Won Ju1orcid, Hee-Il Lee1orcid

Published online: June 11, 2025

1Division of Vectors and Parasitic Diseases, Korea Disease Control and Prevention Agency, Cheongju, Republic of Korea

2Division of Zoonotic and Vector Borne Disease Control, Korea Disease Control and Prevention Agency, Cheongju, Republic of Korea

Corresponding author: Hee-Il Lee Division of Vectors and Parasitic Diseases, Korea Disease Control and Prevention Agency, 187 Osongsaengmyeong 2-ro, Osong-eup, Heungdeok-gu, Cheongju 28159, Republic of Korea E-mail: isak@korea.kr
• Received: February 18, 2025   • Revised: April 2, 2025   • Accepted: May 6, 2025

© 2025 Korea Disease Control and Prevention Agency.

This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

  • 354 Views
  • 18 Download

Full Article

Download PDF Download PDF
  • Objectives
    Anopheles sinensis is a predominant malaria vector found throughout the Republic of Korea (ROK). Population genetic analysis can provide insights into the origins and migration patterns of malaria vectors by assessing genetic variation and distribution among populations.
  • Methods
    In this study, we examined the population genetic structure of An. sinensis in the ROK using mitochondrial cytochrome c oxidase subunit I (COI) marker.
  • Results
    A total of 903 specimens of An. sinensis were collected from 9 sampling regions, including malaria-endemic and non-endemic areas. All populations displayed characteristics typical of migrant populations, with high haplotype diversity (Hd) and low nucleotide diversity (Pi). Network analysis identified 124 haplotypes grouped into 2 clusters. Both clusters included haplotypes from malaria-endemic and non-endemic areas. Cluster I shared its most recent common ancestry with Chinese reference sequences, while cluster II had related with Japanese reference sequences. Pairwise genetic distance (FST) analysis indicated generally low genetic differentiation among populations. Furthermore, FST values tended to increase proportionally with geographical distance between regions. Analysis of molecular variance confirmed that individual mosquitoes within the population had a large effect on the overall variation. Neutrality tests using 4 methods (Tajima’s D, Fu’s Fs, Fu and Li’s D, and Fu and Li’s F) yielded negative values, suggesting that An. sinensis populations are expanding in all studied regions.
  • Conclusion
    This study characterizes the genetic attributes of An. sinensis in the ROK, providing valuable insights into the biology of this important malaria vector and contributing useful data for malaria control strategies.
  • Keywords: Anopheles sinensis; Malaria; Population genetic structure
Malaria causes approximately half a million deaths worldwide each year and contributes significantly to morbidity and mortality in humans [1]. The Plasmodium parasite that causes malaria is transmitted by infected female Anopheles mosquitoes, and 5 species of malaria parasite can infect humans: Plasmodium falciparum, Plasmodium vivax, Plasmodium malariae, Plasmodium ovale, and Plasmodium knowlesi [2]. P. falciparum is predominantly found in malaria-endemic countries in Africa, whereas P. vivax is found in Southeast Asia. These 2 species occasionally cause co-infections in South America, Southeast Asia, and the Western Pacific regions [3]. P. malariae and P. ovale are broadly distributed geographically but often remain undiagnosed due to diagnostic difficulties, low parasitemia levels, and relatively mild virulence [4]. P. knowlesi occurs infrequently in Malaysia and Indonesia [3]. Notably, only P. vivax is prevalent in the Republic of Korea (ROK), with endemic areas primarily situated near the demilitarized zone bordering the Democratic People’s Republic of Korea (DPRK) [5].
Anopheline mosquitoes form genetic structures at different spatial scales, particularly in large, homogeneous geographic regions, due to recent expansions or extensive short-range gene flow [6]. Approximately 70 species of anopheline mosquitoes have been reported to have the ability to transmit malaria parasites to humans [7]. Among these, 8 species (Anopheles belenrae, Anopheles kleini, Anopheles koreicus, Anopheles lesteri, Anopheles lindesayi japonicus, Anopheles pullus, Anopheles sinensis, and Anopheles sineroides) of anopheline mosquitoes inhabit in the ROK [8]. Of these, An. sinensis has been confirmed to exhibit a broad geographical distribution and large population sizes [8,9]. An. sinensis is considered the primary malaria transmission vector in Asia [7], and its vector competence has been demonstrated through artificial sporozoite infection experiments [10]. Thus, An. sinensis is regarded as an important malaria vector in the ROK.
Mitochondrial DNA (mtDNA) serves as a valuable marker for analyzing population genetic structures in various species, including fish, arthropods, and amphibians [11–15]. A Chinese study examining the population genetics of An. sinensis with the mtDNA cytochrome c oxidase subunit I (COI) marker identified 2 regional clusters [16]. Similarly, in the ROK, genetic structuring of An. sinensis was investigated in 2015 using the mtDNA control region, revealing 2 clusters separated by the Sobaek and Taebaek mountain ranges, which acted as a genetic barrier [17,18].
Recent scientific information on the genetic structure of An. sinensis in the ROK remains insufficient. Approximately 10 years after the previous study, the present study aimed to confirm the genetic characteristics of An. sinensis using the COI marker and provide comprehensive molecular and biological data for this mosquito species, recognized as the main vector of malaria in the ROK.
Mosquito Collection and Identification
Mosquitoes were collected from 7 malaria-endemic and 2 malaria-non-endemic areas from June to October 2022 (Figure 1). In malaria-endemic regions, mosquitoes were collected from outdoor of civilian houses using black light traps (Shin-Young Comm. System), while in non-endemic areas, they were collected in livestock sheds using light-emitting diode traps (Biotrap). Collected mosquitoes were visually identified in the laboratory based on morphological taxonomic characteristics and subsequently confirmed through molecular identification by multiplex polymerase chain reaction (PCR) using primers targeting the internal transcribed spacer 2 region [19,20].
DNA Extraction and COI Gene Amplification
DNA was individually extracted from the whole body of mosquitoes using DNAzol (Thermo Fisher Scientific) or the automated nucleic acid extraction system QIAamp 96 DNA QIAcube HT Kit (QIAGEN). The mitochondrial COI gene (710 bp) was amplified using AccuPower PCR PreMix (Bioneer) with forward primer LCO1490 (5’-GGTCAACAAATCATAAAGATATTGG-3’) and reverse primer HCO2198 (5’-TAAACTTCAGGGTGACCAAAAAATCA-3’) [21]. PCR conditions were as follows: initial denaturation at 95 °C for 5 minutes; followed by 40 cycles consisting of denaturation at 95°C for 30 seconds, annealing at 53 °C for 30 seconds, elongation at 72 °C for 30 seconds; and a final extension at 72 °C for 10 minutes. The sizes of the amplified PCR products were confirmed using the QIAxcel capillary electrophoresis system (Qiagen). PCR amplicons were subjected to bidirectional sequencing using Sanger sequencing technology via a commercial sequencing service (BIOFACT).
Data Analyses
Chromatograms of the amplified COI sequences were assembled using DNASTAR Lasergene SeqMan Pro, and consensus sequences were aligned using MEGA 11 (MEGA software). Sequence variation metrics—including the number of haplotypes (H), number of segregating sites (S), average number of nucleotide differences (K), average number of mutations per sequence (θ), haplotype diversity (Hd), and nucleotide diversity (Pi)—were calculated using DnaSP (version 5.10.01; Universitat de Barcelona) [22]. The haplotype network of An. sinensis was generated using the median-joining method in Network 10.2 software (Fluxus Technology Ltd.) [23]. Pairwise genetic distances (FST) among regions were calculated using DnaSP to examine genetic differentiation, and gene flow (Nm) was inferred from the pairwise FST values using the formula Nm=(1–FST)/4FST [24]. Analysis of molecular variance (AMOVA) was performed using Arlequin (ver. 3.5.2.2; University of Geneva) [25] to quantify genetic variation within and among populations. Finally, neutrality tests were conducted using 4 different methods—Tajima’s D, Fu’s Fs, Fu and Li’s D, and Fu and Li’s F—using DnaSP software [22].
Characteristics of Sequences
All 903 mosquito samples were successfully amplified and sequenced with assembled chromatograms (Table 1). After trimming the consensus sequences, 625 bp of the initial 710 bp was used for analysis. Sequence alignment identified 79 variable and 546 conserved sites. S (number of segregating sites) varied from 18 to 49 sites based on the mosquito collection regions. K (average number of nucleotide difference) was 4.35360±0.63288 and θ (average number of mutations per sequence) was 7.38100±1.29023. Hd (haplotype diversity) was high (0.82159±0.08791), whereas Pi (nucleotide diversity) was low (0.00697±0.00101) (Table 2).
COI Haplotype Network
A total of 124 haplotypes were identified from the 903 analyzed An. sinensis sequences. The proportion of haplotype confirmed for the population by region was higher in the endemic regions of Goseong (14/22, 63.6%) and Josan-ri (13/21, 61.9%) compared to other regions (Table 2). Across all 9 collection regions, haplotype 10 (38.6%, 349/903) was most frequently observed, followed by haplotype 20 (6.6%, 60/903), being dominant and subdominant haplotypes, respectively. In total, 24 haplotypes were identified across multiple regions, while 75 haplotypes were represented by only a single individual each. The most dominant haplotype (haplotype 10), representing approximately one-third of the total population, was first detected in June in southern and central regions and became increasingly prevalent over time in northern regions. Similarly, haplotype 20 was initially identified in the southern and central regions in June and subsequently found in northern regions. Haplotype patterns found across multiple regions typically emerged first in the southern and central regions before July (summer) and were later confirmed in northern areas from summer onward (Figure 2). Haplotypes identified from single individuals occurred sporadically, predominantly after summer.
The COI haplotype network of An. sinensis was separated into 2 distinct clusters. Cluster I included mosquitoes collected from all 9 sampling locations and showed relationships to 7 Chinese reference sequences. In particular, haplotypes 6, 20, 34, and 38 shared identical COI sequences with mosquitoes from Jiangsu, Anhui, and Chongqing, China. Conversely, cluster II included mosquitoes from 8 locations (excluding Josan-ri) and was closely related to 1 Japanese reference sequence (Figure 3).
Population Genetic Structure
Pairwise FST was used to calculate the genetic distance between populations in the regions where mosquitoes were collected. According to the classification of He et al. [26], FST values between regions were categorized into small genetic differentiation (0–0.05) and moderate genetic differentiation (0.05–0.15). Small genetic differentiation predominated, indicating that the genetic distance between the regions was not large. Gene flow was calculated based on FST; however, it could not be computed when FST values were negative (Table 3). When comparing FST values to straight-line geographic distances between regions, a general trend of increasing FST values with increasing geographic distance was observed, although exceptions occurred (Table 4).
AMOVA was conducted by dividing the populations into 3 groups based on geographic barriers: east (Gijang and Goseong) and west (Beopheung-ri, Baekyeon-ri, Cheorwon, Gyeyang, Josan-ri, and Yesan) range of the mountain, and island (Baengnyeong). The AMOVA results indicated that genetic variation within populations was substantially greater than the variation among groups or among populations within groups (Table 5). This highlights that individual mosquito within populations significantly influenced the total genetic variation, more than among groups and among populations within groups.
Neutrality Test
Four neutrality tests—Tajima’s D, Fu’s Fs, Fu and Li’s D, and Fu and Li’s F—were performed using DnaSP, and produced all negative values (Table 6). These results indicate that the An. sinensis populations are undergoing population with expansion, characterized by an excess of low-frequency mutations. Notably, the Beopheung-ri region showed strong negative values (p<0.05) for Fu and Li’s D and Fu and Li’s F tests, suggesting heightened sensitivity to neutrality deviations compared to other regions.
Understanding the population genetic structure, migration patterns, and gene flow of vectors provides crucial information for effectively controlling vector-borne diseases [27]. The present study assessed the genetic structure of An. sinensis populations in the ROK using the COI gene, a representative marker in population genetics studies [28].
An. sinensis is widely distributed across most regions of the ROK [18,29] and recognized as the primary malaria vector in the region, as evidenced by its confirmed sporozoite infection capability in artificial infection experiments [10]. In this study, 903 specimens of An. sinensis collected from 9 different regions in the ROK were analyzed for population genetic structure using the COI marker. Our results indicate that the An. sinensis populations exhibit characteristics of migrant populations, characterized by high Hd and low Pi. In particular, Josan-ri, located near the DPRK border, likely represents a strong migrant population due to its notably high Hd and low Pi values (Table 2).
Haplotypes were initially identified in southern and central regions and subsequently appeared in northern regions over time (Figure 2). Although the number of collected mosquitoes is different by month, all collected individuals were analyzed without selective bias. An. sinensis is known to be cold-intolerant in the ROK [30], and adult mosquitoes in the field start emerging at a minimum temperature of approximately 15.3 °C [31]. In June, collection site temperatures were 18.6 °C in Gijang and 19.86 °C in Yesan, approximately 0.5–3.4 °C warmer than in northern regions. While various factors, including habitat and nutritional availability, influence adult mosquito emergence, temperature likely has a significant impact. Hence, earlier detection of haplotypes in southern regions is presumably due to warmer conditions. Additionally, by tracking the months in which haplotypes first appeared, we observed increases in mosquito numbers and haplotype diversity during the warmer summer months (July and August; Table 1). Given that developmental cycles of An. sinensis stimulate with rising temperatures [32], mutations may also occur more frequently during the summer due to accelerated mosquito growth cycles.
Network analysis of An. sinensis revealed 2 distinct clusters (Figure 3). Cluster I contained specimens from all 9 sampling regions and showed close relationships with 7 reference sequences from China. Specifically, haplotypes 20, 34, and 38 had identical COI sequences to those from Anhui and Chongqing provinces in China. Previous studies have documented insect pests, such as Nilaparvata lugens, Sogatella furcifera, and Laodelphax striatellus, introduced from China to the ROK via jet streams [33]; mosquitoes may have similarly migrated. Cluster II, lacking haplotypes from Josan-ri (the northernmost site), was related to a Japanese reference sequence and comprised diverse haplotypes from Gijang. This strong relationship likely arises from Gijang’s geographical proximity to Japan. However, a limitation of this study is the inclusion of only 1 southern sampling site (Gijang); thus, further studies from additional southern locations are necessary to confirm this connection. Both clusters exhibited associations with international sequences, likely reflecting the geographic position of the ROK between China and Japan.
FST was used to assess regional population differentiation. According to the classification by He et al. [26] in 2022, FST values indicate small (0–0.05), moderate (0.05–0.15), large (0.15–0.25), or very large (≥0.25) genetic differentiation. In our analysis, small genetic differentiation was predominant among the An. sinensis populations in the ROK (Table 3). Given the broad distribution of An. sinensis across the ROK, we analyzed FST values relative to straight-line geographic distances between sampling regions (Table 4). Among the sampled regions, Gijang exhibited moderate genetic differentiation, likely due to its significant distance from other regions. Interestingly, Josan-ri showed genetic differentiation from all other regions, including Baekyeon-ri, located less than 1 km away. This distinct differentiation suggests unique regional characteristics, identifying Josan-ri mosquitoes as a strong migrant group. Baengnyeong, an island in the northwest ROK, was anticipated to exhibit high genetic differentiation due to its isolation. However, the results differed from expectations, possibly due to inland mosquito movement influenced by prevailing westerly winds
AMOVA results revealed that genetic variation among regions was predominantly driven by individual mosquitoes within populations. Despite group divisions based on geographic barriers, these barriers did not significantly influence genetic differentiation. Neutrality tests (Tajima’s D, Fu’s Fs, Fu and Li’s D, and Fu and Li’s F) consistently yielded negative values, indicating that An. sinensis populations in the ROK are undergoing population expansion following a recent bottleneck or selection event. Such demographic expansions typically arise through haplotype mutations, but sufficient time may not have elapsed for extensive sequence differences to accumulate [28].
Previous studies examined the population genetic structure of An. sinensis in the ROK using the mitochondrial control region marker, reporting southern and northern genetic groups separated by the Sobaek and Taebaek mountain ranges [17,18]. Due to differences in genetic markers, direct comparisons are difficult; however, the COI marker employed in our study likely provided greater resolution, successfully identifying genetic differentiation even among regions previously classified into the same northern group. One limitation of our study was the occurrence of negative FST values. To overcome this, future studies may consider incorporating additional markers, such as alternative mitochondrial regions or microsatellite loci.
Mosquito movement and migration are challenging to predict, as they are influenced by geographical barriers, environmental factors, and human activities [34]. The present study utilized COI markers to elucidate the genetic population structure of An. sinensis in the ROK, confirming significant regional genetic differences. Such molecular biological analyses provide foundational data on malaria vectors, enhancing our understanding of malaria epidemiology and control strategies.
• The Republic of Korea is endemic for malaria, with Anopheles sinensis as a major vector.
• Genetic analysis of the Anopheles sinensis populations collected in this study showed small regional genetic differentiation.
• The cytochrome c oxidase subunit I haplotype network for An. sinensis comprised 2 distinct clusters, and certain haplotypes demonstrated homology with reference sequences from China.

Ethics Approval

Not applicable.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Funding

This research was funded by the Korea Disease Control and Prevention Agency (Grant No. 6300-6332-305).

Availability of Data

All data generated or analyzed during this study are included in this published article. For other data, there may be requested through the corresponding author.

Authors’ Contributions

Conceptualization: HIS, HIL; Data curation: HJ; Formal analysis: HJ; Investigation: BGH; Methodology: HIS; Project administration: HIS, MRL, JWJ, HIL; Visualization: HJ; Writing–original draft: HJ; Writing–review & editing: all authors. All authors read and approved the final manuscript.

Additional Contributions

We are thankful to the Gangwon State Institute of Health and Environment, Gyeonggi Province Institute of Health and Environment, Incheon Metropolitan City Institute of Health and Environment, Chungcheongnam-do Institute of Health and Environment, and Busan Metropolitan City Institute of Health and Environment for their assistance in collecting mosquito samples.

Figure 1.
Anopheles sinensis sampling sites on the map. The site names are abbreviated.
figure
Figure 2.
Monthly occurrence of haplotypes simultaneously identified in multiple regions. Pink indicates southern and central regions (Gijang and Yesan), blue indicates northern regions (Beopheung-ri, Baengnyeong Island, Baekyeon-ri, Cheorwon, Goseong, Gyeyang, and Josan-ri), and red indicates both northern, central, and southern regions.
figure
Figure 3.
Haplotype network based on the cytochrome c oxidase subunit Ⅰ (COⅠ) gene as calculated by Network 10.2. Each color circle represents individuals from Beopheung-ri (BH), Baengnyeong island (BN), Baekyeon-ri (BY), Cheorwon (CW), Gijang (GJ), Gyeyang (GY), Josan-ri (JS), and Yesan (YS). The black circle symbolizes reference sequence: Anhui (AH1: MG816536, China; AH2: MG816544, China), Chongqing (CQ1: MG816556, China; CQ2: MG816550, China), Ibaraki (IB: LC054430, Japan), Jiangsu (Ji: MF322628, China), and Yunnan (YN1: MG816568, China; YN2: MG816562, China). White dots are a median vector (undetected hypothetical haplotypes). The size of the circle is proportional to the number of individuals.
figure
Table 1.
Information on the collection sites and collected sample number by month and sites
Table 1.
Location ID Coordinates No. of collected samples
 Total
June July August September October
BH 37°47'01.2"N 126°42'27.3"E 0 11 71 85 2 169
BN 37°57'28.9"N 124°39'52.5"E 0 96 96 0 0 192
BY 37°55'08.0"N 126°44'03.3"E 0 13 53 23 1 90
CW 38°15'49.0"N 127°09'51.7"E 1 8 81 51 0 141
GJ 35°11'55.7"N 129°12'09.1"E 20 21 23 24 0 88
GS 38°32'43.7"N 128°23'59.4"E 0 1 20 1 0 22
GY 37°34'49.0"N 126°44'51.6"E 0 11 17 52 0 80
JS 37°54'37.3"N 126°43'53.0"E 0 3 6 2 10 21
YS 36°40'27.5"N 126°41'39.4"E 44 52 4 0 0 100

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 2.
Summary data for sequence characteristics of Anopheles sinensis
Table 2.
Location ID n H S K θ Hd Pi
BH 169 39 43 4.18618 8.06430 0.81298 0.00670
BN 192 47 49 4.97044 9.08763 0.86971 0.00795
BY 90 24 34 3.60175 6.70418 0.72659 0.00576
CW 141 30 34 4.32006 6.33780 0.77285 0.00691
GJ 88 38 41 4.83490 8.71484 0.92659 0.00774
GS 22 14 28 5.36364 7.68100 0.90043 0.00858
GY 80 24 37 4.51867 8.07595 0.68006 0.00723
JS 21 13 18 3.49524 5.00314 0.91905 0.00559
YS 100 29 33 3.89152 6.76018 0.78606 0.00623

n, sample size; H, no. of haplotypes; S, no. of segregating sites; K, average no. of nucleotide difference; θ, average no. of mutations per sequence; Hd, haplotype diversity; Pi, nucleotide diversity; BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 3.
FST (pairwise genetic distance) and Nm (gene flow) values among populations
Table 3.
FST Nm
BH BN BY CW GJ GS GY JS YS
BH 63.5255 117.6745 –5,000.2500 4.0685 –26.1032 43.5328 14.5429 –115.4574
BN 0.0039 26.4309 51.8333 7.2598 –48.5125 16.8968 13.8028 92.0009
BY 0.0021 0.0094 –833.5833 2.7675 23.5369 –111.3611 5.5090 –76.9371
CW –0.0001 0.0048 –0.0003 3.4039 –52.5513 –80.6359 5.4526 567.9318
GJ 0.0579 0.0333 0.0829 0.0684 10.3793 2.6712 7.3465 3.5581
GS –0.0097 –0.0052 0.0105 –0.0048 0.0235 44.7950 861.8190 57.3537
GY 0.0057 0.0146 –0.0023 –0.0031 0.0856 0.0056 3.6963 36.8420
JS 0.0169 0.0178 0.0434 0.0438 0.0329 0.0003 0.0634 11.7577
YS –0.0022 0.0027 –0.0033 0.0004 0.0657 0.0043 0.0067 0.0208

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 4.
FST (pairwise genetic distance) and straight-line distance (km) between collection sites
Table 4.
FST Distance
BH BN BY CW GJ GS GY JS YS
BH 181.33 15.23 66.79 364.04 170.93 23.07 14.29 123.29
BN 0.0039 181.91 221.23 507.23 332.03 187.74 181.71 228.91
BY 0.0021 0.0094 53.71 374.31 161.43 37.63 0.98 138.04
CW –0.0001 0.0048 –0.0003 386.00 112.01 84.04 54.61 180.96
GJ 0.0579 0.0333 0.0829 0.0684 378.80 344.28 373.03 279.04
GS –0.0097 –0.0052 0.0105 –0.0048 0.0235 179.95 162.02 257.72
GY 0.0057 0.0146 –0.0023 –0.0031 0.0856 0.0056 36.73 100.80
JS 0.0169 0.0178 0.0434 0.0438 0.0329 0.0003 0.0634 137.32
YS –0.0022 0.0027 –0.0033 0.0004 0.0657 0.0043 0.0067 0.0208

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 5.
Analysis of molecular variance results for Anopheles sinensis
Table 5.
Source of variation Degrees of freedom Sum of squares Variance components Percentage of variance (%) F-index p
Among groups 2 27.930 0.04714 2.10 0.02094 0.11241±0.00945
Among populations within groups 6 17.926 0.00926 0.41 0.00420 0.06940±0.00785
Within populations 894 1,961.704 2.19430 97.49 0.02506 0.0000±0.0000
Total 902 2,007.560 2.25070 100
Table 6.
Neutrality test for Anopheles sinensis
Table 6.
Location ID Neutrality tests
Tajima’s D Fu’s Fs Fu and Li’s D Fu and Li’s F
BH –1.44182 –19.13100 –2.46241* –2.43457*
BN –1.35769 –24.04200 –1.25227 –1.56128
BY –1.44273 –8.64200 –1.52977 –1.79123
CW –0.94790 –9.51100 –0.52387 –0.83794
GJ –1.42160 –23.13900 –2.17915 –2.25224
GS –1.15273 –3.90100 –1.43444 –1.57528
GY –1.42160 –23.13900 –2.17915 –2.25224
JS –1.12075 –5.40000 –0.81795 –1.05626
YS –1.31062 –12.68300 –1.82358 –1.94258

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

*p<0.05.

  • 1. Garrido-Cardenas JA, Gonzalez-Ceron L, Manzano-Agugliaro F, et al. Plasmodium genomics: an approach for learning about and ending human malaria. Parasitol Res 2019;118:1−27.
  • 2. Milner DA. Malaria pathogenesis. Cold Spring Harb Perspect Med 2018;8:a025569.
  • 3. World Health Organization (WHO) (CH). World malaria report 2023 [Internet]. WHO; 2023 [cited 2025 Feb 2]. Available from: https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2023.
  • 4. Poespoprodjo JR, Douglas NM, Ansong D, et al. Malaria. Lancet 2023;402:2328−45.
  • 5. Seong H, Suh J, Choi JY, et al. Development of a Plasmodium vivax malaria model for evaluating the effects of control strategies on the malaria burden in democratic people's Republic of Korea. Front Public Health 2024;12:1423004.
  • 6. Besansky NJ, Lehmann T, Fahey GT, et al. Patterns of mitochondrial variation within and between African malaria vectors, Anopheles gambiae and An. arabiensis, suggest extensive gene flow. Genetics 1997;147:1817−28.
  • 7. Sinka ME, Bangs MJ, Manguin S, et al. A global map of dominant malaria vectors. Parasit Vectors 2012;5:69.
  • 8. Rueda LM, Brown TL, Kim HC, et al. Species composition, larval habitats, seasonal occurrence and distribution of potential malaria vectors and associated species of Anopheles (Diptera: Culicidae) from the Republic of Korea. Malar J 2010;9:55.
  • 9. Foley DH, Klein TA, Kim HC, et al. Geographic distribution and ecology of potential malaria vectors in the Republic of Korea. J Med Entomol 2009;46:680−92.
  • 10. Lee HW, Cho SH, Shin EH, et al. Experimental infection of Anopheles sinensis with Korean isolates of Plasmodium vivax. Korean J Parasitol 2001;39:177−83.
  • 11. Bunmee K, Thaenkham U, Saralamba N, et al. Population genetic structure of the malaria vector Anopheles minimus in Thailand based on mitochondrial DNA markers. Parasit Vectors 2021;14:496.
  • 12. Joshi BD, Johnson JA, Negi T, et al. Understanding genetic diversity and population genetic structure of three Cyprinidae fishes occupying the same habitat from Uttarakhand, India. Mitochondrial DNA B Resour 2019;4:2956−61.
  • 13. King'ori EM, Obanda V, Nyamota R, et al. Population genetic structure of the elephant tick Amblyomma tholloni from different elephant populations in Kenya. Ticks Tick Borne Dis 2022;13:101935.
  • 14. Nakamura Y, Yamagishi J, Hayashida K, et al. Genetic diversity and population structure of Glossina morsitans morsitans in the active foci of human African trypanosomiasis in Zambia and Malawi. PLoS Negl Trop Dis 2019;13:e0007568.
  • 15. Zhang L, Sechi P, Xie J, et al. Genetic diversity and population structure of the endangered salamander hynobius yiwuensis inferred from mitochondrial DNA Sequences. Life (Basel) 2024;14:769.
  • 16. Feng X, Huang L, Lin L, et al. Genetic diversity and population structure of the primary malaria vector Anopheles sinensis (Diptera: Culicidae) in China inferred by cox1 gene. Parasit Vectors 2017;10:75.
  • 17. Jung J, Jung Y, Min GS, et al. Analysis of the population genetic structure of the malaria vector Anopheles sinensis in South Korea based on mitochondrial sequences. Am J Trop Med Hyg 2007;77:310−5.
  • 18. Kang S, Jung J, Kim W. Population genetic structure of the malaria vector anopheles sinensis (Diptera: Culicidae) sensu stricto and evidence for possible introgression in the Republic of Korea. J Med Entomol 2015;52:1270−81.
  • 19. Multiplex PCR assay for the identification of eight Anopheles species belonging to the Hyrcanus, Barbirostris and Lindesayi groups. Malar J 2021;20:287.
  • 20. Joshi D, Park MH, Saeung A, et al. Multiplex assay to identify Korean vectors of malaria. Mol Ecol Resour 2010;10:748−50.
  • 21. Folmer O, Black M, Hoeh W, et al. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol 1994;3:294−9.
  • 22. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009;25:1451−2.
  • 23. Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 1999;16:37−48.
  • 24. Slatkin M. Gene flow and the geographic structure of natural populations. Science 1987;236:787−92.
  • 25. Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 2010;10:564−7.
  • 26. He ZQ, Zhang QQ, Wang D, et al. Genetic polymorphism of circumsporozoite protein of Plasmodium falciparum among Chinese migrant workers returning from Africa to Henan Province. Malar J 2022;21:248.
  • 27. Zhang HD, Gao J, Li CX, et al. Genetic diversity and population genetic structure of Aedes albopictus in the Yangtze River Basin, China. Genes (Basel) 2022;13:1950.
  • 28. Chaiphongpachara T, Laojun S, Changbunjong T, et al. Genetic diversity, haplotype relationships, and kdr mutation of Malaria Anopheles Vectors in the most plasmodium knowlesi-Endemic Area of Thailand. Trop Med Infect Dis 2022;7:412.
  • 29. Lee JH, Kim HW, Mustafa B, et al. The relationships between microbiome diversity and epidemiology in domestic species of malaria-mediated mosquitoes of Korea. Sci Rep 2023;13:9081.
  • 30. Lee J, Lee IY, Yong TS, et al. Preliminary observation on the gradient distribution of hibernating females of Anopheles sinensis and Anopheles pullus in the Republic of Korea. J Am Mosq Control Assoc 2008;24:121−2.
  • 31. Mogi M, Okazawa T. Development of Anopheles sinensis immatures (Diptera: Culicidae) in the field: effects of temperature and nutrition. Med Entomol Zool 1996;47:355−62.
  • 32. Feng X, Zhang S, Huang F, et al. Biology, bionomics and molecular biology of Anopheles sinensis Wiedemann 1828 (Diptera: Culicidae), main Malaria vector in China. Front Microbiol 2017;8:1473.
  • 33. Park T, Choe H, Jeong H, et al. Spatial pattern analysis for distribution of migratory insect pests at paddy field in Jeolla-province. Korean J Appl Entomol 2018;57:361−72. Korean.
  • 34. Md Naim D, Kamal NZ, Mahboob S. Population structure and genetic diversity of Aedes aegypti and Aedes albopictus in Penang as revealed by mitochondrial DNA cytochrome oxidase I. Saudi J Biol Sci 2020;27:953−67.

Download Citation

Download a citation file in RIS format that can be imported by all major citation management software, including EndNote, ProCite, RefWorks, and Reference Manager.

Format:

Include:

Genetic diversity and spatiotemporal population structure of Anopheles sinensis in the Republic of Korea based on the mitochondrial cytochrome c oxidase subunit I (COI) marker
Download Citation

Download a citation file in RIS format that can be imported by all major citation management software, including EndNote, ProCite, RefWorks, and Reference Manager.
Format:
Include:
Genetic diversity and spatiotemporal population structure of Anopheles sinensis in the Republic of Korea based on the mitochondrial cytochrome c oxidase subunit I (COI) marker
Close

Figure

  • 0
  • 1
  • 2
Genetic diversity and spatiotemporal population structure of Anopheles sinensis in the Republic of Korea based on the mitochondrial cytochrome c oxidase subunit I (COI) marker
Image Image Image
Figure 1. Anopheles sinensis sampling sites on the map. The site names are abbreviated.
Figure 2. Monthly occurrence of haplotypes simultaneously identified in multiple regions. Pink indicates southern and central regions (Gijang and Yesan), blue indicates northern regions (Beopheung-ri, Baengnyeong Island, Baekyeon-ri, Cheorwon, Goseong, Gyeyang, and Josan-ri), and red indicates both northern, central, and southern regions.
Figure 3. Haplotype network based on the cytochrome c oxidase subunit Ⅰ (COⅠ) gene as calculated by Network 10.2. Each color circle represents individuals from Beopheung-ri (BH), Baengnyeong island (BN), Baekyeon-ri (BY), Cheorwon (CW), Gijang (GJ), Gyeyang (GY), Josan-ri (JS), and Yesan (YS). The black circle symbolizes reference sequence: Anhui (AH1: MG816536, China; AH2: MG816544, China), Chongqing (CQ1: MG816556, China; CQ2: MG816550, China), Ibaraki (IB: LC054430, Japan), Jiangsu (Ji: MF322628, China), and Yunnan (YN1: MG816568, China; YN2: MG816562, China). White dots are a median vector (undetected hypothetical haplotypes). The size of the circle is proportional to the number of individuals.

Figure 1.

Figure 2.

Figure 3.

Genetic diversity and spatiotemporal population structure of Anopheles sinensis in the Republic of Korea based on the mitochondrial cytochrome c oxidase subunit I (COI) marker
Location ID Coordinates No. of collected samples
 Total
June July August September October
BH 37°47'01.2"N 126°42'27.3"E 0 11 71 85 2 169
BN 37°57'28.9"N 124°39'52.5"E 0 96 96 0 0 192
BY 37°55'08.0"N 126°44'03.3"E 0 13 53 23 1 90
CW 38°15'49.0"N 127°09'51.7"E 1 8 81 51 0 141
GJ 35°11'55.7"N 129°12'09.1"E 20 21 23 24 0 88
GS 38°32'43.7"N 128°23'59.4"E 0 1 20 1 0 22
GY 37°34'49.0"N 126°44'51.6"E 0 11 17 52 0 80
JS 37°54'37.3"N 126°43'53.0"E 0 3 6 2 10 21
YS 36°40'27.5"N 126°41'39.4"E 44 52 4 0 0 100
Location ID n H S K θ Hd Pi
BH 169 39 43 4.18618 8.06430 0.81298 0.00670
BN 192 47 49 4.97044 9.08763 0.86971 0.00795
BY 90 24 34 3.60175 6.70418 0.72659 0.00576
CW 141 30 34 4.32006 6.33780 0.77285 0.00691
GJ 88 38 41 4.83490 8.71484 0.92659 0.00774
GS 22 14 28 5.36364 7.68100 0.90043 0.00858
GY 80 24 37 4.51867 8.07595 0.68006 0.00723
JS 21 13 18 3.49524 5.00314 0.91905 0.00559
YS 100 29 33 3.89152 6.76018 0.78606 0.00623
FST Nm
BH BN BY CW GJ GS GY JS YS
BH 63.5255 117.6745 –5,000.2500 4.0685 –26.1032 43.5328 14.5429 –115.4574
BN 0.0039 26.4309 51.8333 7.2598 –48.5125 16.8968 13.8028 92.0009
BY 0.0021 0.0094 –833.5833 2.7675 23.5369 –111.3611 5.5090 –76.9371
CW –0.0001 0.0048 –0.0003 3.4039 –52.5513 –80.6359 5.4526 567.9318
GJ 0.0579 0.0333 0.0829 0.0684 10.3793 2.6712 7.3465 3.5581
GS –0.0097 –0.0052 0.0105 –0.0048 0.0235 44.7950 861.8190 57.3537
GY 0.0057 0.0146 –0.0023 –0.0031 0.0856 0.0056 3.6963 36.8420
JS 0.0169 0.0178 0.0434 0.0438 0.0329 0.0003 0.0634 11.7577
YS –0.0022 0.0027 –0.0033 0.0004 0.0657 0.0043 0.0067 0.0208
FST Distance
BH BN BY CW GJ GS GY JS YS
BH 181.33 15.23 66.79 364.04 170.93 23.07 14.29 123.29
BN 0.0039 181.91 221.23 507.23 332.03 187.74 181.71 228.91
BY 0.0021 0.0094 53.71 374.31 161.43 37.63 0.98 138.04
CW –0.0001 0.0048 –0.0003 386.00 112.01 84.04 54.61 180.96
GJ 0.0579 0.0333 0.0829 0.0684 378.80 344.28 373.03 279.04
GS –0.0097 –0.0052 0.0105 –0.0048 0.0235 179.95 162.02 257.72
GY 0.0057 0.0146 –0.0023 –0.0031 0.0856 0.0056 36.73 100.80
JS 0.0169 0.0178 0.0434 0.0438 0.0329 0.0003 0.0634 137.32
YS –0.0022 0.0027 –0.0033 0.0004 0.0657 0.0043 0.0067 0.0208
Source of variation Degrees of freedom Sum of squares Variance components Percentage of variance (%) F-index p
Among groups 2 27.930 0.04714 2.10 0.02094 0.11241±0.00945
Among populations within groups 6 17.926 0.00926 0.41 0.00420 0.06940±0.00785
Within populations 894 1,961.704 2.19430 97.49 0.02506 0.0000±0.0000
Total 902 2,007.560 2.25070 100
Location ID Neutrality tests
Tajima’s D Fu’s Fs Fu and Li’s D Fu and Li’s F
BH –1.44182 –19.13100 –2.46241* –2.43457*
BN –1.35769 –24.04200 –1.25227 –1.56128
BY –1.44273 –8.64200 –1.52977 –1.79123
CW –0.94790 –9.51100 –0.52387 –0.83794
GJ –1.42160 –23.13900 –2.17915 –2.25224
GS –1.15273 –3.90100 –1.43444 –1.57528
GY –1.42160 –23.13900 –2.17915 –2.25224
JS –1.12075 –5.40000 –0.81795 –1.05626
YS –1.31062 –12.68300 –1.82358 –1.94258
Table 1. Information on the collection sites and collected sample number by month and sites

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 2. Summary data for sequence characteristics of Anopheles sinensis

n, sample size; H, no. of haplotypes; S, no. of segregating sites; K, average no. of nucleotide difference; θ, average no. of mutations per sequence; Hd, haplotype diversity; Pi, nucleotide diversity; BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 3. FST (pairwise genetic distance) and Nm (gene flow) values among populations

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 4. FST (pairwise genetic distance) and straight-line distance (km) between collection sites

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

Table 5. Analysis of molecular variance results for Anopheles sinensis
Table 6. Neutrality test for Anopheles sinensis

BH, Beopheung-ri; BN, Baengnyeong Island; BY, Baekyeon-ri; CW, Cheorwon; GJ, Gijang; GS, Goseong; GY, Gyeyang; JS, Josan-ri; YS, Yesan.

p<0.05.

Table 1.

Table 2.

Table 3.

Table 4.

Table 5.

Table 6.

PHRP : Osong Public Health and Research Perspectives

About

Browse articles 

Editorial policies

For authors  

phrp 사이트 QRCODE
© Korea Disease Control and Prevention Agency.