Molecular evidence for new sympatric cryptic species of Aedes albopictus (Diptera: Culicidae) in China: A new threat from Aedes albopictus subgroup?

Aedes (Stegomyia) albopictus (Skuse) is an indigenous species and the predominant vector of dengue fever in China. Understanding of genetic diversity and structure of the mosquito would facilitate dengue prevention and vector control. Sympatric cryptic species have been identified in the Ae. albopictus subgroup in Southeast Asia; however, little is known about the presence and distribution of cryptic species in China. This study aimed to examine the genetic diversity, evaluate potential new cryptic sibling species, and assess the prevalence of Wolbachia infections in field populations. Aedes adult female specimens were collected from five provinces in southern and central China during 2015–2016. Morphological identification was performed under dissection microscope. The mitochondrial DNA cytochrome c oxidase subunit 1 (cox1, DNA barcoding) locus and the ribosomal DNA internal transcribed spacer region 2 (ITS2) marker were used to examine the genetic variation, evaluate cryptic sibling species, and population structure in the field populations. Screening for the presence of Wolbachia was performed using multiplex PCR. A total of 140 individual specimens with morphological characteristics similar to Ae. albopictus were sequenced for DNA barcoding. Among these, 129 specimens (92.1%) were confirmed and identified as Ae. albopictus. The remaining 11 specimens, from 2 provinces, were identified as 2 distinct sequence groups, which were confirmed by ITS2 marker sequencing, suggesting the existence of potential cryptic species of Ae. albopictus. In Ae. albopictus, we found significant genetic differentiation and population structure between populations collected from different climate zones. Medium to high frequencies of Wolbachia infections were observed in natural Ae. albopictus populations, whereas Wolbachia was infrequent or absent in cryptic species populations. Our findings highlight the population differentiation by climate zone and the presence of novel, cryptic Aedes species in China. The low prevalence of Wolbachia infections in cryptic species populations could reflect either a recent invasion of Wolbachia in Ae. albopictus or different host immune responses to this symbiont in the cryptic species. The study provides useful information for vector control and host-symbiont coevolution. Further study is needed to investigate the potential for arbovirus infection and disease transmission in the emerged cryptic species.


Background
Over the past decades, more than 11 sibling species or cryptic species have been identified and characterized in the Aedes albopictus subgroup of the Scutellaris group in the subgenus Stegomyia of Aedes [1][2][3]. Of these species, Ae. albopictus, originating from Asia, is the most widely distributed and has invaded on every continent except Antarctica [4,5]. Aedes albopictus is considered a medically important species and is a major vector of several human arboviruses, including dengue, Zika, chikungunya, yellow fever and West Nile viruses [5][6][7][8][9][10]. Dengue fever has experienced a 30-fold increase in incidence worldwide over the past 50 years and shows no signs of slowing down [11]. Since the 1970s there have been several major outbreaks of dengue fever in southern China, including in Hainan, Guangxi, Fujian, Zhejiang, Yunnan and Guangdong provinces [12][13][14]. Due to climate change, the transmission of the dengue virus has spread gradually from southern tropical or subtropical regions to the surrounding northern and western regions, and even to the central China Henan Province with a generally warm temperate continental climate [15]. The most recent outbreak of dengue fever occurred in 2014 in Guangdong Province, with a total of 45,224 dengue fever cases and 6 deaths [12,16,17]. Aedes albopictus mosquitoes are regarded as the sole vector for dengue transmission in nearly all these epidemics [14,18].
In China, Ae. albopictus is an indigenous species, closely associated with human migration, transportation, commerce and urbanization. It is the most important dengue vector species and has different susceptibilities to dengue virus in different geographical areas [19,20]. Due to the lack of effective treatments or vaccines for dengue fever, vector control through chemical or biological measures targeting mosquitoes or their breeding sites is essential for dengue prevention. With the progressive spread of insecticide resistance, the threat of Ae. albopictus is growing, and development of efficient surveillance methods is more urgent than ever before [21,22]. Population genetic studies of arthropod disease vectors can provide information about the transmission dynamics of specific pathogens, which aids in the design of strategies for controlling vector-borne disease epidemics [23,24]. The recent waves of dengue outbreak in China highlight the need to improve our knowledge of Ae. albopictus population distribution and dynamics. Although scientists have studied the diversity of the cox1 gene in Ae. albopictus in several localities in China [9,[25][26][27], there have been no systematic studies of the genetic diversity of Ae. albopictus field populations and its cryptic species.
Cryptic species are defined as sibling species of two or more morphologically indistinguishable biological groups that are closely related and live in the same habitat [28]. Cryptic species may be medically important in vectorborne disease transmission, vector ecology and evolutionary biology. A number of new cryptic species have been identified in mosquito genera (Diptera: Culicidae), including Culex [29][30][31], Anopheles [32][33][34][35][36][37][38][39] and Aedes [40,41]. In Ae. albopictus, a novel cryptic species has been reported in Vietnam [40], and the divergence between the cryptic species and Ae. albopictus was confirmed by analysis of nuclear ribosomal genes and mitochondrial genes. However, there are no reports of cryptic species of Ae. albopictus in other Asian countries, including China.
Natural infections of Wolbachia microbiota are common in Ae. albopictus, and the two Wolbachia biotypes, wAlbA and wAlbB, co-occur at a rate near 100% in many areas [42][43][44][45][46]. Maternally inherited Wolbachia bacteria can cause cytoplasmic incompatibility (CI) in many insect species, including Ae. albopictus mosquitoes [47,48]. Wolbachia mediates antiviral protection of Aedes mosquitoes against a broad range of RNA viruses, including dengue, yellow fever, chikungunya and Rift Valley fever virus [49]. Understanding the distribution and prevalence of Wolbachia in Ae. albopictus and its cryptic species will provide useful information for vector control and host-symbiont coevolution.
In this study, we investigated the genetic diversity and population structure of Ae. albopictus from different climate regions in China, uncovered and molecularly identified the cryptic Aedes species and its polymorphism, and detected Wolbachia infection in the natural Aedes populations. The prevalence of Wolbachia endosymbiont was evaluated by multiplex PCR genotyping and DNA sequencing of individuals in the natural Aedes populations.

Sample collection
Adult Aedes mosquito specimens were collected from April 2015 to October 2016 using BG-sentinel traps (Bioquip Products, Inc. California, USA) or electric aspirator mosquito catches [50] at 14 collection sites in five provinces in China: Henan, Guangdong, Guangxi, Yunnan and Hainan (Fig. 1). These sites are highly diverse in environmental conditions and most of them have experienced dengue epidemics in the past. The sampling site in Henan Province, located in central China, has a temperate climate with a distinct seasonality characterized by hot, humid summers and generally cold, windy and dry winters. The sampling sites in Guangdong and Guangxi provinces have a subtropical monsoon climate with long summers and year-round abundant precipitation. The sampling sites in Hainan and Yunnan provinces are tropical areas with a wet climate. Hainan Province had dengue epidemics in the late 1970s and early 1980s, however, no dengue epidemic has been reported from there since 1990. Guangdong Province has experienced multiple major dengue and chikungunya epidemics since 1980 and dengue has remained every year since 1994. In 2013, Guangdong, Yunnan and Henan provinces had dengue outbreaks. Dengue transmission occurred in southern China from July to November and the peak season is usually September and October whereas dengue epidemic in Henan Province is limited in summer from July to September, and showed earlier peaks and shorter epidemic periods [15]. Aedes albopictus mosquitoes are the primary vector of dengue virus across China, especially in urban areas. Aedes aegypti mosquitoes are only found in a small portion of southern China, including Hainan Province and small portions of Yunnan Province and southern tip of Guangdong Province [51]. Ten mosquito specimens from each collection site were used in the study. All mosquito specimens were morphologically identified under a stereomicroscope (Nikon) using morphological keys as described by Lu et al. [52]. All mosquito samples were stored at -20°C prior to DNA extraction.
PCR amplification and sequencing of mitochondrial DNA (mtDNA) Total DNA was extracted from individual adult mosquitoes using the Insect DNA Kit (OMEGA Bio-Tek, D0926-01, Guangzhou, China) according to the manufacturer's standard protocol. Extracted DNA was preserved at -20°C until molecular analysis. The mitochondrial gene cytochrome c oxidase subunit 1 (cox1) was used to examine sequence polymorphism among mosquito samples. PCR was performed to amplify a 651 bp fragment of the 5' cox1 region of mtDNA using the DNA primer pairs LCOI490 (5'-GGT CAA CAA ATC ATA AAG ATA TTG G-3') and HCO2198 (5'-TAA ACT TCA GGG TGA CCA AAA AAT CA-3') [53,54]. PCR amplification was performed in a 25 μl reaction volume with 12.5 μl GoTaq Green Master Mix (Promega, Guangzhou, China), 1 μl each of the forward and reverse primers at 10 μmol/l, 2 μl of template DNA and sufficient nuclease-free water to make 25 μl. PCR conditions were as follows: an initial denaturation at 94°C for 1 min followed by five cycles of 94°C for 40 s (denaturation), 45°C for 40 s (annealing), and 72°C for 1 min (extension); 30 cycles of 94°C for 40 s (denaturation), 53°C for 40 s (annealing), and 72°C for 1 min (extension); and a final extension at 72°C for 5 min. The amplified fragments were run on a 1% agarose gel to check integrity, stained with ethidium bromide and analyzed under UV light. PCR products were purified using a gel extraction kit (OMEGA Bio-Tek, D2500-02) and sequenced with PCR primers in both directions The internal transcribed spacer 2 (ITS2) region of ribosomal DNA was amplified from the DNA samples using the universal primers ITS2A (5'-ATC ACT CGG CTC GTG GAT CG-3') and ITS2B (5'-ATG CTT AAA TTT AGG GGG TAG TC-3'), which anneal to highly conserved sequences in the 5.8S and 28S rDNA genes flanking the entire ITS2 region [21,55]. PCR amplification was performed in a 25 μl reaction volume with 12.5 μl GoTaq Green Master Mix (Promega, Guangzhou, China), 1 μl each of the forward and reverse primers at 10 μmol/l, 2 μl of template DNA (1~2 ng/μl), and sufficient nuclease-free water to make 25 μl. PCR conditions were as follows: an initial denaturation at 94°C for 3 min followed by 30 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 1 min; and a final extension at 72°C for 5 min. Purification and sequencing of the PCR products were the same as described above for the cox1 gene. The sequences of ITS2 unique haplotypes were submitted to the GenBank database under the accession numbers MF623839-MF623851.

PCR detection of Wolbachia infection in mosquitoes
The Wolbachia infection status of individual mosquitoes was determined by PCR amplification of Wolbachia ribosomal DNA using primers specific for Wolbachia 16S rDNA (WF: 5'-CAT ACC TAT TCG AAG GGA TAG-3' and WR: 5'-AGC TTC GAG TGA AAC CAA TTC-3') [56]. To further classify infected mosquitoes by Wolbachia group, we amplified the Wolbachia surface protein gene (wsp) using wAlbA primers (328F: 5'-CCA GCA GAT ACT ATT GCG-3' and 691R: 5'-AAA AAT TAA ACG CTA CTC CA-3') for A group and wAlbB primers (183F: 5'-AAG GAA CCG AAG TTC ATG-3' and 691R: 5'-AAA AAT TAA ACG CTA CTC CA-3') for B group [57]. PCR amplification was performed in a 25 μl reaction volume with 12.5 μl GoTaq Green Master Mix (Promega, Guangzhou, China), 1 μl each of the forward and reverse primers at 10 μmol/l, 2 μl of template DNA, and sufficient nuclease-free water to make 25 μl. PCR conditions were as follows: an initial denaturation at 94°C for 3 min followed by 30 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 1 min; and a final extension at 72°C for 5 min. PCR-amplified fragments of 408 bp, 364 bp, and 509 bp for 16S rDNA, wAlbA and wAlbB, respectively, were revealed under UV light after electrophoresis on 1% agarose gel. Negative and positive controls for the PCR assay were included in each run.
To obtain the positive control, we sequenced PCR fragments from the 16S rDNA and wsp genes and confirmed that the amplified PCR product was Wolbachia by using BLAST search to compare it with existing sequences in the NCBI database.

Data analysis
The cox1 gene sequences from 140 mosquitoes were aligned using Clustal W multiple alignment in BioEdit (version 7.2.6.1) [58]. The number of segregating sites, haplotype diversity (H d ), and nucleotide diversity (π) within each population were determined using DnaSP version 5 [59]. Pairwise sequence divergences were calculated using a Kimura 2-parameter (K2P) distance model in MEGA 7.0.20 [60]. The K2P model was used to make our results comparable with most other studies on mosquito DNA barcoding. To examine population expansion, we also performed neutrality tests for each population. Deviations from selective neutrality were tested using Fu's F s statistic [61] and Tajima's D [62]. To determine the genealogical relationships among haplotypes, we constructed a haplotype network using a statistical parsimony algorithm implemented in TCS version 1.21 [63]. The minimum number of mutational steps between sequences was calculated with > 95% confidence. A haplotype network shows the haplotype frequencies in each population and their relatedness, which is useful in inferring the plausible geographical origin of a population [64]. Genetic differentiation among populations was estimated using Arlequin 3.5 [65]. Analysis of molecular variance (AMOVA) was conducted to determine the distribution of genetic variation within and among populations and among groups (tropical, subtropical, and temperate zone).

Results
Genetic polymorphism of Ae. albopictus and its cryptic species PCR amplification and sequencing of the mitochondrial cox1 gene resulted in a 651 bp fragment for each individual study subject, with no insertions or deletions. We compared the cox1 sequences with the existing sequences in the NCBI database by BLAST search. Of the 140 individuals, 129 sequences (92.1%) were identical or possessed > 98% similarity with Ae. albopictus (GenBank: KR068634) (Additional file 1: Table  S1). The remaining 11 individuals (7.9%) were approximately 10% K2P divergent from Ae. albopictus, indicating the existence of cryptic species (namely, Aedes sp.) of Ae. albopictus in China (Additional file 2: Table S2). Of these 11 individuals, 9 were from samples collected in Wuzhou, Guangxi Province, with one each from Baisha and Baoting, Hainan Province. Since there is only one individual of Ae. albopictus in GX-WZ population and one each of Aedes sp. from HN-BS, and HN-BT populations, which is insufficient for population genetic analysis, these individuals were excluded from the analysis of the population genetic structure and genetic diversity. Thus, only 9 individuals from each of these three populations (GX-WZ, HN-BS and HN-BT) were included in the genetic polymorphism analysis (Table 1) Varied genetic diversity was also found in the 8 populations from Hainan Province, with nucleotide diversity (π) ranging from 0.123 in specimens from Qiongzhong to 0. 666 in specimens from Haikou. Tajima's D tests for all study populations were not statistically significant ( Table  1), indicating that the populations are in genetic A total of 57 haplotypes of mtDNA cox1 were detected in the 140 specimens, including 47 haplotypes derived from 129 Ae. albopictus mosquitoes and 10 haplotypes derived from cryptic species (Additional file 3: Table S3, GenBank: KY765450-KY765506). Three predominant haplotypes were identified in Ae. albopictus populations: H01 (21.5%) from Guangdong Province, H19 (22.1%) from Yunnan and Hainan provinces, and H30 (10%) from Hainan Province. Other haplotypes were either unique to a specific population or had a limited geographical distribution (Additional file 3: Table S3). To determine the relationships among the samples, we constructed a median-joining network using haplotypes based on sequence variation. Haplotypes were connected when the probability of parsimony was at least 0.95. Two networks were constructed based on all the 57 haplotypes, one based on haplotypes from Ae. albopictus (Fig. 2a-c) and the other based on haplotypes from the cryptic species (Fig. 2d). Three haplotypes could not be connected to the networks at a 95% confidence level: H28 for Ae. albopictus, and H16 and H44 for the cryptic species. The Ae. albopictus haplotypes can be classified into three clusters corresponding to three climate zones (tropical, subtropical and temperate) as shown in Fig. 2a-c. Haplotype H01 from Guangdong Province (subtropical) was connected to haplotypes from Hainan Province (tropical) through haplotype H30 by one mutation step at nucleotide position 342 (T-342-C), and to haplotypes from Henan Province (temperate) through H26 by one mutation step at nucleotide position 624 (G-624-A). Similar patterns were also observed for haplotype H19 from Yunnan and Hainan provinces (tropical zones). These results may imply a multiple origin for Ae. albopictus populations in China. In the cryptic species, haplotype H17 had at least 4 connections with other haplotypes, suggesting it as a potential ancestral haplotype (Fig. 2d).
A phylogenetic tree based on cox1 sequence variation indicated three clades, which were assigned to three different species (Fig. 3). The first clade (solid blue circles in Fig. 3) corresponds to Ae. albopictus from China and Vietnam, the second clade (red squares) corresponds to cryptic Aedes species previously identified in Vietnam, and the third clade (purple diamond) corresponds to cryptic Aedes species previously identified in Pakistan. The haplotype Aedes sp. CH-H44 (BS07) (purple diamond) is a different cryptic species compared to the haplotypes from Wuzhou, Guangxi Province (red squares). A similar pattern was observed in the phylogenetic tree based on ITS2 sequence variation (Fig. 4). The first clade (solid blue circles in Fig. 4) corresponds Fig. 2 Cox1 haplotype networks showing the genealogical relationships. a-c Aedes albopictus. d Cryptic Aedes species. Each haplotype is represented by a pie chart with size proportional to its frequency in each population. A black dotted line indicates that a mutation step could not be determined between haplotypes at probability of parsimony above the 0.95 limit to Ae. albopictus from China and Vietnam and the second clade (red squares) corresponds to cryptic Aedes species previously identified in Vietnam. The ITS2 haplotype of Aedes sp. CH-H44 (BS07) (purple diamond) is also in different clade compared to the haplotypes from Wuzhou, Guangxi Province (red squares). The two taxa were distinguished using PCR length polymorphism at the ITS2 locus. Aedes albopictus had an amplicon size of~580 bp, whereas the cryptic species had an amplicon size of~415 bp, allowing for easy and accurate identification of the cryptic species through 1% agarose gel electrophoresis.  (Table 2). No genetic differentiation was observed between populations within Hainan Province or between populations from Yunnan and Hainan provinces, indicating strong gene flow between these populations. All pairwise F ST values of differentiation between cryptic population (GX-WZ) and the other 13 Ae. albopictus populations were highly significant and generally very high (F ST values > 0.9) ( Table 2). To further examine population structure and the extent of genetic variation between Ae. albopictus populations from tropical zones (Yunnan and Hainan provinces) and those from subtropical and temperate zones (Guangdong and Henan provinces), we conducted analysis of molecular variation (AMOVA) between the two groups of populations. Our results   indicated a significant overall population structure in Ae. albopictus (F ST = 0.38, P < 0.001). The majority of genetic variation (61.77%) was within populations, whereas approximately 34.21% was between the two groups, and only 4.01% of variation was among populations within groups (Table 3).

PCR detection of Wolbachia infection in natural mosquito populations
Wolbachia infections were detected in all 14 of the Ae. albopictus populations. Infection rates ranged from 50% (HN-DZ) to 100% (in 8 populations) with an average of 90% ( Table 4), suggesting that Wolbachia is highly prevalent in Ae. albopictus in China. In the cryptic Aedes sp. species, however, Wolbachia infection was absent (HN-BT and HN-BS) or occurred at low frequency (11%, GX-WZ) (Fig. 5a-d). Most infected individuals were infected with both the wAlbA and wAlbB strains of Wolbachia; the average superinfection rate in the 14 Ae. albopictus populations was 70%, with a range from 10% (HN-DZ) to 90% (YN-JH, HN-CJ and HN-QZ). Singlestrain Wolbachia infections were found in five populations for wAlbA and eight populations for wAlbB, with low prevalence (< 20%) (Fig. 5e, f ) except in HN-CM (70% wAlbA) ( Table 4).

Discussion
Aedes albopictus, one of the cryptic species of the Ae. albopictus subgroup, is an important vector for public health. It is highly invasive and is the most widely distributed mosquito species in the world. Aedes albopictus originated at the edges of forests and bred in natural habitats, but it has adapted to being a domestic mosquito [13]. Today, the species can be found throughout tropical, subtropical and temperate zones in China. Cryptic species often occur in sympatry and are so similar they cannot be distinguished via traditional species identification using morphological keys [68,69]. With the advent of molecular diagnosis technology and relatively inexpensive DNA sequencing technology, the discovery of cryptic species has become common in many insect groups [70]. However, little is known about the cryptic species within the Ae. albopictus subgroup in China. To the best of our knowledge, this is the first report describing the discovery of phylogenetically divergent cryptic species living in sympatry with Ae. albopictus in China. In this study, among 14 study populations collected across the tropical, subtropical and temperate climate zones of China, we found three populations in southern China in which cryptic species of Ae. albopictus coexisted in sympatry. Furthermore, multiple haplotypes of cryptic species were present in the Wuzhou population in Guangxi Province within the subtropical zone, and probably a novel cryptic species (KY765493) existed in the Baisha population in the tropical area of Hainan Province. Further investigation is needed to confirm the reproductive isolation of these species. These results may have important implications for vector control and for understanding the evolutionary processes of the species. For example, if the cryptic species is a novel disease vector with different biting or resting behaviors, current vector control interventions that target on species-specific vector behavior could lead directly to programme failure and thus these vector control strategies need to be adjusted.
The distribution of Aedes mosquito species is influenced by climatic, environmental and geographical factors, as well as by human behavior [71][72][73]. In this study, we found significantly higher genetic diversity in cryptic species populations than in Ae. albopictus populations. For example, the cryptic species population in Wuzhou (GX-WZ) has 33 segregation sites, with much higher nucleotide diversity and the highest average number of nucleotide differences compared with the Ae. albopictus populations. The high genetic diversity of cryptic species, as well as their coexistence with Ae. albopictus, may be explained by environmental heterogeneity in these densely forested, mountainous areas with a low level of gene flow and random genetic drift. The different population structure between tropical (Hainan and Yunnan) and subtropical (Guangdong)/ temperate (Henan) populations of Ae. albopictus might be due to selective pressures exerted by specific climate, environment and human activities. The mountainous areas of Yunnan and central Hainan provinces have lower human population density and more complex tropical environments than the densely populated Guangdong (coastal) and Henan (plain) provinces. The association between genetic population structure and   climate/environment has also been observed in other Ae. albopictus populations [19,74,75]. Gene flow between the two genetic clusters (tropical and subtropical/temperate) appeared to be restricted due to geographical isolation, which was also evidenced in the cox1 haplotype network analysis. Interestingly, both of the major tropical-zone haplotypes, H30 and H19, were just one mutation step (C-342-T) from subtropical (Guangdong Province) and temperate (Henan Province) haplotypes, and the same mutation step (G-624-A) was found between populations within tropical zones as well as between populations from subtropical and temperate zones. These results suggest that climate, geography, environment and human activity all play important roles in Ae. albopictus population structure. The maternally transmitted endosymbiotic bacterium Wolbachia is known to have an important impact on host reproduction in many insects, including mosquitoes (Diptera: Culicidae) [47,76,77]. Wolbachia can inhibit human pathogens transmitted by mosquitoes, including dengue virus [78][79][80], yellow fever [81], filarial nematodes [82,83], malaria parasites [84][85][86] and Zika virus [87]. Our results indicate significant variations in the frequency of Wolbachia infection in Ae. albopictus populations. Populations from Guangdong, Yunnan, and many sites in Hainan Province had a 100% infection rate, followed by a 90% infection rate in Henan Province in central China. Such high prevalence of Wolbachia has also been observed in other Asian countries, including Malaysia [43,88], Thailand [89], India [90] and Sri Lanka [91]. In our study, most populations were naturally infected with two Wolbachia strains (wAlbA and wAlbB) at high frequencies, suggesting that superinfection is common in Ae. albopictus, as observed in other studies [44,88,89,92]. Interestingly, a relatively low prevalence of Wolbachia was also observed in three Ae. albopictus populations (HN-DZ, HN-LS, and HN-BT) from Hainan Province, in which no more than 50% of the individuals were not detected with Wolbachia infection, indicating that there was variation in natural Wolbachia infection in these Ae. albopictus populations. Furthermore, no Wolbachia infection was detected in the Aedes cryptic species from Hainan Province, and Wolbachia infection was detected in just one individual of the cryptic species from Guangxi Province. A similar pattern was observed in the cryptic species in Vietnam [40], suggesting that the cryptic species may be resistant to Wolbachia infection or that Wolbachia are present at low cell density in the cryptic species that cannot be detected by PCR [93]. In addition, cryptic species may prevent Wolbachia introgression by reproductive isolation and maintaining ancestral levels of mitochondrial diversity. Wolbachia induced cytoplasmic incompatibility and mitochondrial selective sweep have been observed in the Ae. albopictus and other mosquito species [29,48]. Further studies are needed to confirm the reproductive isolation between Ae. albopictus and its cryptic species.

Conclusions
Our results indicated that the genetic diversity and population structure of Ae. albopictus between tropical, subtropical and temperate zones in China appeared to be separated by a single mutation step at the mitochondrial DNA barcoding cox1 gene. Sympatric, cryptic sibling species might be common in the Ae. albopictus subgroup in China. The prevalence of high-level Wolbachia infection in most of the Ae. albopictus populations, and the absence or low prevalence of Wolbachia in the sympatric cryptic species, possibly due to Wolbachia-induced genetic hitchhiking or selective sweep that has created a barrier to gene flow among the species. Elucidating the mechanisms of the observed absence or low prevalence of Wolbachia in sympatric cryptic species may provide insight toward the development of new vector control strategies. Finally, this study will have important implications for disease vector-based control programs, Wolbachia-based disease control strategies, and host evolutionary biology. Further study is needed to investigate the potential for arbovirus infection and disease transmission in the emerged cryptic species.