Next Article in Journal
Mapping QTLs Controlling Soybean Rust Disease Resistance in Chiang Mai 5, an Induced Mutant Cultivar
Next Article in Special Issue
Spatio-Temporal Evolutionary Patterns of the Pieridae Butterflies (Lepidoptera: Papilionoidea) Inferred from Mitogenomic Data
Previous Article in Journal
Comparison of Marker Effects and Breeding Values at Two Levels at THI for Milk Yield and Quality Traits in Brazilian Holstein Cows
Previous Article in Special Issue
Complete Mitogenome Analysis of Five Leafhopper Species of Idiocerini (Hemiptera: Cicadellidae)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phylogenetic Analysis and Genetic Structure of Schlegel’s Japanese Gecko (Gekko japonicus) from China Based on Mitochondrial DNA Sequences

1
Jiangsu Key Laboratory for Biodiversity and Biotechnology, College of Life Sciences, Nanjing Normal University, Nanjing 210023, China
2
School of Food Science, Nanjing Xiaozhuang University, Nanjing 211171, China
*
Author to whom correspondence should be addressed.
Genes 2023, 14(1), 18; https://doi.org/10.3390/genes14010018
Submission received: 10 October 2022 / Revised: 7 December 2022 / Accepted: 17 December 2022 / Published: 21 December 2022

Abstract

:
Gekko japonicus, i.e., Schlegel’s Japanese Gecko, is an important species which is widely distributed in East Asia. However, the information about population genetics of this species from China remains unclear. To address this issue, we used sequences from a fragment of the mitochondrial protein-coding gene cytochrome c oxidase I to estimate genetic diversity, genetic structure, and historical demography of G. japonicus populations from China. Phylogenetic analysis indicated that G. japonicus had a close relationship with Gekko wenxianensis. A total of 14 haplotypes were obtained, of which haplotype 1 was the most common and widely distributed. The genetic diversity of G. japonicus was comparatively low across different geographic populations. The populations of G. japonicus were divided into four groups which exhibited low levels of genetic differentiation, and expressed an unclear pattern of population structuring. In addition, potential population expansion of G. japonicus has occurred as well. Overall, these results demonstrate that the populations of G. japonicus reveal low genetic diversity in China, which is attributed to the founder and bottleneck events among populations. Our results will provide meaningful information on the population genetics of G. japonicus and will provide some insights into the study of origin of populations.

Graphical Abstract

1. Introduction

G. japonicus, known as Schlegel’s Japanese Gecko, is a nocturnal gecko which was reported to be distributed in China, Japan, and Korea [1]. Currently, studies on a variety of topics including population distribution, morphological characteristics, daily activity pattern, selection of oviposition site, and mitochondrial genome, have been performed on the G. japonicus [2]. Undoubtedly, the specialized body mechanism and ecological role make G. japonicus an excellent flagship species for many issues in ecological balance and biological evolution. Predictions of present and future distribution studies which were based on MaxEnt modeling, suggest that suitable habitats are currently located in coastal cities of Japan, China, and Korea, as well as in isolated patches of inland China. Due to climate change, suitable habitats are expected to shrink along coastlines, particularly at the coastal edge of climate change zones [3]. Recent concerns regarding the area of origin and genetic relationships among populations of G. japonicus are growing, as are the range expansions of G. japonicus in Korea and Japan [3]. Genetic diversity and inferred dispersal history demonstrated that the populations of G. japonicus from China is inferred to be the source population, which had higher genetic diversity and more private alleles compared to the populations of G. japonicus in Japan and Korea [4]. In brief, the genetic structure and distribution pattern of G. japonicus in China will provide important information about the origin and evolutionary of this species. Nevertheless, to date, information of its populations, particularly from China, is unavailable. The genetic diversity, genetic structure, and historical demography of G. japonicus from China remains unclear.
The current distribution of G. japonicus in China includes the eastern coast extending westward to eastern Sichuan Province, northward to southern Shaanxi, and Gansu Provinces, which overlaps with Gekko scabridus, Gekko hokouensis, Gekko subpalmatus, and Gekko chinensis [1,3]. G. japonicus and other species of gecko could be separated by morphological and molecular manners [1]. In some cases, however, species complex may display few or no phenotypic differences, making them impossible to be distinguished based on morphology [5,6]. By definition, species that are morphologically indistinguishable, but belong to distinct evolutionary lineages are referred to as cryptic species [7,8,9,10]. Molecular data, especially mitochondrial DNA (mtDNA), has shown to be an efficient tool for identification of species, owing to its relatively rapid evolution rate and lack of recombination [11,12]. Meanwhile, the mtDNA has been broadly applied in population genetics and phylogenetic studies [12,13]. Among these mitochondrial genes, the cytochrome c oxidase subunit 1 (COI) has been vastly employed in identification of species, study of genetic diversity and genetic structure on different lizard species [4,14]. This approach has aided the rapid identification of species and analysis of diversity patterns in diverse groups [15,16]. Specifically, COI genes could provide valuable information for understanding the population structure and variations.
To detect the genetic structure and diversity of G. japonicus, 325 specimens were collected from 37 sampling sites in China. The sequences of COI genes were obtained from these 325 specimens in China. Subsequently, phylogenetic analyses of G. japonicus with other species of Gekko were conducted to clarify the genetic relationship among the genus Gekko, which is necessary for the further understanding of the evolutionary history of Gekko. Genetic structure, genetic diversity, and historical demography of G. japonicus populations were determined. To our knowledge, this study is the first to understand the genetic structure, genetic diversity, and historical demography of G. japonicus populations in China based on the mitochondrial gene COI.

2. Materials and Methods

2.1. Sampling

A total of 325 individuals were collected from 37 sampling sites that covered most of the range of G. japonicus in China from 2008 to 2017 (Figure 1). They were identified as G. japonicus and preserved in 95% ethanol for genomic DNA extraction. Geographic coordinates of sampling sites were recorded with a Geographic Positioning System unit (Garmin GPSmap 60CSX, Shanghai, GARMIN, USA) using map datum WGS84. The sample information that was used is shown in Table S1.

2.2. DNA Extraction and Sequencing

Genomic DNA was extracted from liver tissue using DNA extraction kit (Qingke, Nanjing, China) according to the manufacturer’s instructions, and then the quality and integrity of the DNA samples were checked by agarose gel electrophoresis. The extracted DNA was stored at −20 °C. A partial fragment of the COI was amplified by the polymerase chain reaction (PCR) using primers: VF1-d (5′-TTCTCAACCAACCACAARGAYATYGG-3′) and VR1-d (5′-TAGACTTCTGGGTGGCCRAARAAYCA-3′) [17]. Amplifications were carried out in 20 µL reaction volumes containing 10 µL of 2 × Taq PCR Master Mix Ⅱ (TIANGEN, Beijing, China), 0.7 µL of genomic DNA, 0.4 µL of each primer, and 8.5 µL of deionized water. The cycle profiles for the PCR were as follows: initial denaturation of 3 min at 95 °C, followed by 5 cycles of 30 s at 95 °C, annealing for 1 min at 45 °C, an initial extension for 2 min at 72 °C; followed by 35 cycles of 30 s at 95 °C, annealing for 1 min at 51 °C, an initial extension for 2 min at 72 °C; and final extension for 5 min at 72 °C. PCR products were visualized on a 1% agarose gel stained with 4S Red Plus, Bi-directional sequencing was performed using the primers RepCOI-F (5′-TNTTMTCAACANACCACAAAGA-3′) and RepCOI-R (5′-ACTTCTGGRTGKCCAAARAATCA-3′) [18] by Nanjing Qingke Biology Co., Ltd. (Nanjing, China).

2.3. Sequence Alignment and Phylogenetic Analyses

We searched GenBank (National Center for Biotechnology Information, https://www.ncbi.nlm.nih.gov (accessed on 1 February 2022)) for published COI of some species of the genus Gekko. Then forty-four COI sequences were retrieved and downloaded for further analyses (Table S2). All 325 individuals were successfully amplified for COI. The chromatograms of each sequence were proofread and then assembled using SeqMan 12.3 [19]. Multiple sequences were aligned and trimmed using ClustalW implemented in MEGA 7.0 [20], aligned with the COI sequences of part species of G. japonicus group (G. auriverrucosus, G. chinensis, G. hokouensis, G. japonicus, G. scabridus, G. subpalmatus, G. swinhonis, G. wenxianensis), Gekko gecko group (Gekko gecko) and Gekko monarchus group (Gekko kikuchii), and Hemidactylus dracaenacolus and H. granti (Table S2), in which H. dracaenacolus and H. granti were used as outgroups for the phylogenetic analyses of the mtDNA.
The K2P model was used to calculate genetic distances using MEGA 7.0 [21]. For phylogenetic reconstruction, we performed maximum likelihood (ML) and Bayesian inference (BI) by using IQ TREE 1.6.12 [22] and MrBayes 3.2.6 [23]. PartitionFinder2 [24] was run to determine the appropriate model of molecular evolution in a likelihood ratio test framework based on the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) for the ML and BI method, respectively. For the maximum likelihood (ML), bootstrap analyses were performed on 1000 full heuristic replicates [25]. For the Bayesian inference, MarkovChain Monte Carlo chains were run for 10,000,000 generations (sampled every 1000 generations) to allow adequate time for convergence. The first 25% of the sampled trees were considered burn in. The final tree was visualized and edited by ITOL (https://itol.embl.de/ (accessed on 1 February 2022)).

2.4. Population Genetic Analysis

To identify similar groups of populations and to evaluate the amount of genetic variation among the partitions, a spatial analysis of molecular variance was conducted in SAMOVA v1.092 [26]. Several runs were performed using increasing numbers of groups (K = 1–20) and 100 annealing simulations for each K. In each run, populations were clustered into genetically and geographically homogenous groups. The number of groups was chosen so as to maximize genetic differentiation among the groups (FCT) [26]. The numbers of haplotypes, haplotype diversity (h) [27], and nucleotide diversity (π) for each population were estimated using DnaSP 6.0 program [28]. Spatial genetic variation was also quantified by estimating analog of FST using analysis of molecular variance (AMOVA) [29]. In addition, to identify the relationship among haplotypes, the mtDNA haplotypes of G. japonicus were estimated from a TCS network using PopART v1.7 [30]. The relationships of haplotypes were further elucidated using ML analysis implemented in IQTREE. The confidence levels at nodes after 1000 repetitions employed by the bootstrap method. The historical population demography was evaluated using Tajima’s D statistic [31], Fu’s Fs [32] and mismatch distribution test [33]. Mismatch distributions for the species were calculated with the expected frequency based on a population growth-decline model. The sum of squared deviations (SSDs) between observed and expected mismatch distribution and the raggedness index [34] were calculated to test the null hypothesis of spatial expansion using Arlequin 3.5 [35].

3. Results

3.1. Sequence Information

A total of 325 mitochondrial COI gene sequences of G. japonicus were obtained and further analyzed. The sequences were deposited in GenBank under project SUB12126169. The average frequencies of T, C, A, and G were 28.6, 30.7, 21.1, and 19.6%, respectively. The sequenced region contained 564 conserved sites, 75 variable sites, and 16 parsimony informative sites. The ratio of transform to transition (Ts/Tv) based on the K2P base substitution model is 2.2.

3.2. Phylogenetic Analysis

In the study, the ML tree and BI tree showed that all species formed a monophyletic group with high bootstrap and posterior probability values (Figure 2). Species of the G. japonicus group cluster well together and form sister groups with Gekko gecko group and Gekko monarchus group. ML tree and BI tree revealed that the G. swinhonis samples were separated into A1, A2, and B. Both in BI and ML tree, G. japonicus cluster formed a sister group with G. wenxianensis, which suggested that G. japonicus has a close relationship with G. wenxianensis. Then, G. swinhonis was clustered with G. auriverrucosus and the two formed a sister group with G. subpalmatus. In addition, a total of 12 individuals were collected from Anqing, Anhui Province, of which 11 individuals were clustered into G. japonicus, undoubtedly. AHAQ06 draws our attention by forming a sister group with all the G. japonicus. Subsequently, according to the COI genes of the AHAQ06 and 10 gecko species of China, interspecific and intraspecific genetic distance was calculated to further illustrate the reliability of the phylogenetic results (Table 1). Except for AHAQ06, the results indicated that the genetic distance between the G. japonicus and G. wenxianensis was the smallest (21.7–23.5), the genetic distance between the G. subpalmatus and G. gecko was the largest (42.0–44.9). Meanwhile, the genetic distance within the G. auriverrucosus was the smallest (0.3 ± 0.2) and the genetic distance within the G. scabridus was the largest (4.3 ± 0.6). Moreover, the genetic distance between AHAQ06 and G. japonicus had reached 12.0–12.9. The genetic distances largely coincide with the phylogenetic results. Therefore, there are large differences between AHAQ06 and G. japonicas at the molecular level, which indicate that AHAQ06 is a special individual to be further studied. Hence, AHAQ06 is excluded from the population-level data set. A total of 324 mitochondrial COI gene sequences was depicted in further population genetic analysis.

3.3. Genetic Diversity, Genetic Structure and Historical Demography of G. japonicus Population

3.3.1. Genetic Diversity

Fourteen haplotypes were identified among the 324 sequences (AHAQ06 are excluded) in G. japonicus (Table S3). The most frequent haplotype detected was Hap 1, shared by 233 individuals, followed by haplotypes Hap 5 shared by 55 individuals. The least frequent haplotypes were Hap 2, Hap 7, Hap 9, Hap 10, Hap 11, Hap 13, and Hap 14, each of which were recorded and designated as “private haplotypes” (Table 2). Geographic populations with a sample size of more than 5 were selected during the genetic diversity analysis to avoid errors, where possible. The population genetic analysis revealed that the nucleotide diversity (π) and average number of nucleotide differences (k) were highest in the SXYX population, based on COI genes at 0.00241, 1.600. Haplotype diversity between all sequences obtained ranged from 0 to 0.591. The highest h was the HNYZ population (h = 0.591), whereas the lowest was the HNYY and HNCB populations (h = 0). In general, total genetic diversity values were low between populations (π = 0.00176, h = 0.453), compared to Chilabothrus inornatus [36].

3.3.2. Genetic Structure

According to the COI data, SAMOVA revealed that the proportion of the total genetic differences between groups (FCT) was the highest (0.39527) when K = 4 (Figure S1). Four groups were separated as follows: A: Yongan, Fujian, Yangshuo, Guangxi, Huaihua, Hunan; B: Yangxian, Shaanxi. C: Anhui (Anqing, Lu’an, Wuhu), Fujian Nanping, Guangzhou Guiding, Guizhou (Huaxi, Longli, Libo), Guangxi (Guilin, Longsheng), Hubei (Jingmen, Wuhan), Hunan (Chengbu, Changde, Daoxian, Huayuan, Huayuan, Shuangpai, Shaoyang, Tong, Xinhuang, Yueyang, Yongzhou, Xinning, Yongzhou Lengshuitan, Zhuzhou, Yongzhou Qingtang Yueyyan Forest Farm), Jiangsu Rugao, Jiangxi Longnan, Zhejiang (Yongzhou Qingtang Yueyyan Forest Farm). Zhoushan, Lishui). D: Yongfu, Guangxi. The total variation was FST = 0.40362 and the intragroup variation was FSC = 0.01382 when K = 4. The AMOVA results revealed that the source of genetic differences emerged primarily from within the populations (Table 3). The groups differed significantly from each other (FCT = 0.39527) and were responsible for 39.527% of the variance. The non-significant differentiation among populations within groups (Fsc = 0.01382) were responsible for 0.84% of the variance. Additionally, the FST = 0.40362 (p < 0.05) indicated low genetic diferentiation among four groups (Table 3).

3.3.3. Historical Demography

A total of fourteen haplotypes forms a network exhibiting a star like topology with the haplotype Hap 1 and Hap 5 at the center (Figure 3). Hap1 was found to be the most widespread, with all 37 populations investigated harboring haplotype Hap 1. A total of 14 populations investigated harbored haplotype Hap 5. Other haplotypes were formed by mutations of these two haplotypes. Hap 1 contains 233 individuals, which might be the ancestor haplotype which evolved into others. In addition, the phylogenetic relationship of all haplotypes based on COI gene was determined (Figure S2). The results showed that all haplotypes from G. japonicus formed one cluster and then cluster with G. swinhonis. There is no clear genetic lineage formed in populations of G. japonicus. Both Tajima’s D and Fu’s Fs were negative values (Tajima’s D = −2.77466, p < 0.05, and Fu’s Fs = −2.19112, p > 0.05, respectively), indicating a distinct departure from a null hypothesis with selective neutrality and population equilibrium. The results of the mismatch analysis of G. japonicus did not exclude the sudden expansion model due to non-convergence. The neutrality tests results implied the possibility of population expansion [37,38]. The nonsignificant SSD statistic (SSD = 0.0771, p > 0.05) and the HRag value (HRag = 0.259, p > 0.05) under the spatial expansion models failed to reject the spatial expansion model.

4. Discussion

In this study, the mitochondrial COI gene was explored for the identification and analysis of the genetic diversity and genetic structure of G. japonicus in China. Phylogenetic analyses of G. japonicus with 9 other species of geckos were also conducted to clarify the genetic relationship among the genus Gekko. Additionally, the genetic diversity and genetic structure of G. japonicus in China was evaluated via mitochondrial COI. To the best of our knowledge, this study is the first to examine the mitochondrial genetic diversity of G. japonicus in China using the COI gene. Our study contributes a taxonomic dataset and also provides significant insights into the population genetics of an ecologically vital reptile species in China.
The phylogenetic relationship of G. japonicus group was well supported, and the G. swinhonis were divided into A1, A2, and B. This result was consistent with a previous work using mitochondrial cytochrome b gene and 10 microsatellite loci, which revealed two genetic lineages in G. swinhonis [39]. Since AHAQ06 clustered with all the G. japonicus, we employed genetic distance to demonstrate the reliability of the phylogenetic analysis and the specificity of the AHAQ06. Phylogenetic analyses showed that G. japonicus has a close genetic relationship to G. wenxianensis in the genus Gekko (Figure 2), this result is the same as the one obtained in the previous paper [40]. The current distribution of G. japonicus is southern Shaanxi and Gansu provinces in China where it borders G. wenxianensis. Hence, their close genetic relationship might be explained by geographic isolation of an independent branch of G. japonicus in Wenxian county, which led to speciation, which will give us a new insight into the speciation and dispersal of G. japonicus.
The single individual AHAQ06, clustered with all other G. japonicus samples, but is genetically distinct, may represent a new species. However, this study is limited by the lack of samples. More samples are needed for further studies of molecular and morphological characteristics.
An overall low level of genetic diversity was detected for G. japonicus from China, considering the low π and h, despite a significant variability of environmental conditions between the geographic regions examined. Genetic diversity in populations is the result of the interplay between mutation, genetic drift, selection, and gene flow [41,42]. Genetic diversity can be lost from populations through various mechanisms. There may be a vital factor responsible for low genetic diversity. Low genetic diversity could be associated with founder and bottleneck events followed by a potential recent population expansion. Demographic bottleneck will increase the action of genetic drift and cause a general loss of genetic diversity, in a magnitude determined by its severity and duration [42,43], as has already been observed in African elephants [44], black-footed ferrets [45], and Arctic foxes [46]. Thus, low genetic diversity showed the possibility of a recent bottleneck for the populations of G. japonicus in China, resulting in a smaller original population before an expansion or bottleneck. Time lapse is not sufficient to accumulate sufficient genetic diversity as well as counter the effects of initial demographic bottleneck. In addition, more genetic data, especially the nuclear genome, need to be collected for further studies of the evolutionary history and genetic structure of G. japonicus. To better comprehend the genetic diversity of G. japonicus at the species level, collecting more samples of G. japonicus is therefore highly required.
SAMOVA provided the highest significant FCT (0.39527) when k = 4 was the partition (Figure S1). Subsequently, the AMOVA results demonstrate the distinct genetic structure when populations were grouped. The source of genetic differences emerged primarily from within the populations and rarely from among populations within groups. The absence of variation from among populations within groups may be due to human activities that have promoted gene flow. Previous studies have revealed that contemporary genetic flow is especially high between east coast China and west and central coast Japan, and between west and central coast Japan and Korea. Considering that geckos are frequently invasive species [47,48], it is highly probable that these contemporary mixing processes are a result of anthropogenic transportation [4]. Considering that high commercial trade exists within countries, human-assisted dispersal may have played a vital role in the diversification of G. japonicus in China. The genetic differences among groups were at a low level. Among them, Group B and Group D are located at the western and southwestern edges of the Chinese range of G. japonicus, respectively. Considering these, we infer that the low level of genetic differentiation between groups is possibly associated with the combination of geographic blockage and gene flow. Despite the human-assisted dispersal factor, the effects of climatic and geographical factors were not sufficient to be countered. Overall, the genetic differentiation of G. japonicus was significant, which was likely due to environmental heterogeneous selection. Founder and bottleneck events are further corroborated by the negative values in the neutrality tests (except Fu’s Fs). Long-distance dispersal events tend to result in reduced genetic diversity due to sequential founding events and genetic drift [49]. Moreover, the potential population expansion was supported by the “star-like” cluster originating from the haplotypes H1 and H5. Hap 1 might be the ancestor haplotype. Star-like patterns in sequence networks are generally considered to be characteristic of potential recent range expansion [50]. In the present study, the pattern of haplotype distribution is consistent with potential recent range expansion events. Therefore, our study shows that the G. japonicus populations have historically experienced founder and bottleneck events causing low genetic diversity. Additionally, phylogenetic relationships between haplotypes show that no clear genetic lineage formed in populations of G. japonicus. Seven private haplotypes are found, which may be partially explained by limited migration and gene flow among them.
It is worth noting that the sample sizes of G. japonicus are not equal among the groupings in our study and this may have an effect on the diversity measures. Hence, more elaborate studies involving a relatively uniform sample size may uncover the mitochondrial genetic diversity of G. japonicus from of China. Admittedly, as our results are based solely on mtDNA sequence data, they should be viewed with some caution, since our conclusions might be biased due to the intrinsic mutation rate, mode of inheritance, and effective population size of mtDNA markers. Thus, further studies integrating ecological, morphological, and molecular (using mtDNA and nuclear markers) data are necessary to understand the evolutionary history of G. japonicus from China. Such studies will also give deep insight into the evolution of G. japonicus from China. It will further set the base for further genetic studies that would help design effective Ecological Resource Management policies. Also, the genetic approach used in our study can be extended in studying the population genetics of other reptile species in China.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14010018/s1, Figure S1: SAMOVA results of genetic structure of G. japonicus. The grouping result is best when K = 4 with the maximum value of FCT. Figure S2: Phylogenetic relationships of 14 haplotypes resolved in G. japonicus based on COI gene. Same color indicates haplotypes of the same radiation set. The shade of color represents the number of haplotypes. Table S1: 325 Sample information of G. japonicus used in experiments in this study. Table S2: Species sources and accession numbers of COI gene of Gekko downloaded from GenBank. Table S3: Sample information included in each haplotype of the G. japonicus COI gene.

Author Contributions

P.L., J.Y. and K.Z. conceived and designed the experiments. F.C. performed the experiments. F.C., L.X., S.C. and Y.C. analyzed data. L.X., F.C., P.L. and J.Y. prepared the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Grant No. 31672269 to JY), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Institutional Review Board Statement

Animal welfare and experimental procedures were approved by the Institutional Animal Care and Use Committee of Nanjing Normal University [SYXK (Jiangsu) 2020-0047, IACUC-20210236]. The authors confirmed that animals did not suffer unnecessarily at any stage of experiments in this study.

Informed Consent Statement

Not applicable.

Data Availability Statement

All raw sequencing reads have been deposited in the NCBI Sequence Read Archive under project SUB12126169. This article contains supplementary figures and tables, which are available to authorized users.

Acknowledgments

We thank four anonymous reviewers, guest editors (Jiasheng Hao and Jiayong Zhang) and editors for helpful comments on the manuscript. We would like to thank Yonglei Zhang, Zhiming Lei, Qi Xiao, Sijia Yuan, Luwen Zhang, Shan Su, Hui Chen and Chaochao Hu for their help in the experiments and preparation of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rösler, H.; Bauer, A.M.; Heinicke, M.P.; Greenbaum, E.; Jackman, T.; Nguyen, T.Q.; Ziegler, T. Phylogeny, Taxonomy, and Zoogeography of the Genus Gekko Laurenti, 1768 with the Revalidation of G. reevesii Gray, 1831 (Sauria: Gekkonidae). Zootaxa 2011, 2989, 1. [Google Scholar] [CrossRef]
  2. Kim, I.-H.; Park, J.; Cheon, K.-S.; Lee, H.-J.; Kim, J.-K.; Park, D. Complete Mitochondrial Genome of Schlegel’s Japanese Gecko Gekko japonicus (Squamata: Gekkonidae). Mitochondrial DNA Part A 2016, 27, 3684–3686. [Google Scholar] [CrossRef] [PubMed]
  3. Kim, D.-I.; Park, I.-K.; Bae, S.-Y.; Fong, J.J.; Zhang, Y.-P.; Li, S.-R.; Ota, H.; Kim, J.-S.; Park, D. Prediction of Present and Future Distribution of the Schlegel’s Japanese Gecko (Gekko japonicus) Using MaxEnt Modeling. J. Ecol. Environ. 2020, 44, 5. [Google Scholar] [CrossRef]
  4. Kim, J.-S.; Park, J.; Fong, J.J.; Zhang, Y.-P.; Li, S.-R.; Ota, H.; Min, S.-H.; Min, M.-S.; Park, D. Genetic Diversity and Inferred Dispersal History of the Schlegel’s Japanese Gecko (Gekko japonicus) in Northeast Asia Based on Population Genetic Analyses and Paleo-Species Distribution Modelling. Mitochondrial DNA A DNA Mapp. Seq. Anal. 2020, 31, 120–130. [Google Scholar] [CrossRef]
  5. Nava, S.; Beati, L.; Labruna, M.B.; Cáceres, A.G.; Mangold, A.J.; Guglielmone, A.A. Reassessment of the Taxonomic Status of Amblyomma cajennense (Fabricius, 1787) with the Description of Three New Species, Amblyomma tonelliae n. sp, Amblyomma interandinum n. sp. and Amblyomma patinoi n. sp., and Reinstatement of Amblyomma mixtum, and Amblyomma sculptum (Ixodida: Ixodidae). Ticks Tick Borne Dis. 2014, 5, 252–276. [Google Scholar] [CrossRef]
  6. Elmer, K.R.; Kusche, H.; Lehtonen, T.K.; Meyer, A. Local Variation and Parallel Evolution: Morphological and Genetic Diversity across a Species Complex of Neotropical Crater Lake Cichlid Fishes. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2010, 365, 1763–1782. [Google Scholar] [CrossRef] [Green Version]
  7. Struck, T.H.; Feder, J.L.; Bendiksby, M.; Birkeland, S.; Cerca, J.; Gusarov, V.I.; Kistenich, S.; Larsson, K.-H.; Liow, L.H.; Nowak, M.D.; et al. Finding Evolutionary Processes Hidden in Cryptic Species. Trends Ecol. Evol. 2018, 33, 153–163. [Google Scholar] [CrossRef] [Green Version]
  8. De Sousa-Paula, L.C.; Pessoa, F.A.C.; Otranto, D.; Dantas-Torres, F. Beyond Taxonomy: Species Complexes in New World Phlebotomine Sand Flies. Med. Vet. Entomol. 2021, 35, 267–283. [Google Scholar] [CrossRef]
  9. Fišer, C.; Robinson, C.T.; Malard, F. Cryptic Species as a Window into the Paradigm Shift of the Species Concept. Mol. Ecol. 2018, 27, 613–635. [Google Scholar] [CrossRef] [Green Version]
  10. Gill, B.A.; Kondratieff, B.C.; Casner, K.L.; Encalada, A.C.; Flecker, A.S.; Gannon, D.G.; Ghalambor, C.K.; Guayasamin, J.M.; Poff, N.L.; Simmons, M.P.; et al. Cryptic Species Diversity Reveals Biogeographic Support for the “mountain Passes Are Higher in the Tropics” Hypothesis. Proc. Biol. Sci. 2016, 283, 20160553. [Google Scholar] [CrossRef]
  11. Piganeau, G.; Gardner, M.; Eyre-Walker, A. A Broad Survey of Recombination in Animal Mitochondria. Mol. Biol. Evol. 2004, 21, 2319–2325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Galtier, N.; Nabholz, B.; Glémin, S.; Hurst, G.D.D. Mitochondrial DNA as a Marker of Molecular Diversity: A Reappraisal. Mol. Ecol. 2009, 18, 4541–4550. [Google Scholar] [CrossRef] [PubMed]
  13. Singh, T.R. Mitochondrial Gene Rearrangements: New Paradigm in the Evolutionary Biology and Systematics. Bioinformation 2008, 3, 95–97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Solov’eva, E.N.; Poiarkov, N.A.; Dunaev, E.A.; Duĭsebaeva, T.N.; Bannikova, A.A. Molecular differentiation and taxonomy of the sunwatcher toadheaded agama species complex Phrynocephalus superspecies helioscopus (Pallas, 1771) (reptilia: Agamidae). Russ. J. Genet. 2011, 47, 842–856. [Google Scholar] [CrossRef]
  15. Laopichienpong, N.; Muangmai, N.; Supikamolseni, A.; Twilprawat, P.; Chanhome, L.; Suntrarachun, S.; Peyachoknagul, S.; Srikulnath, K. Assessment of Snake DNA Barcodes Based on Mitochondrial COI and Cytb Genes Revealed Multiple Putative Cryptic Species in Thailand. Gene 2016, 594, 238–247. [Google Scholar] [CrossRef]
  16. Hebert, P.D.N.; Cywinska, A.; Ball, S.L.; deWaard, J.R. Biological Identifications through DNA Barcodes. Proc. Biol. Sci. 2003, 270, 313–321. [Google Scholar] [CrossRef] [Green Version]
  17. Ivanova, N.V.; Dewaard, J.R.; Hebert, P.D.N. An Inexpensive, Automation-Friendly Protocol for Recovering High-Quality DNA. Mol. Ecol. Notes 2006, 6, 998–1002. [Google Scholar] [CrossRef]
  18. Nagy, Z.T.; Sonet, G.; Glaw, F.; Vences, M. First Large-Scale DNA Barcoding Assessment of Reptiles in the Biodiversity Hotspot of Madagascar, Based on Newly Designed COI Primers. PLoS ONE 2012, 7, e34506. [Google Scholar] [CrossRef]
  19. Swindell, S.R.; Plasterer, T.N. SEQMAN. Contig Assembly. Methods Mol. Biol. 1997, 70, 75–89. [Google Scholar]
  20. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  21. Gintis, H. The Hitchhiker’s Guide to Altruism: Gene-Culture Coevolution, and the Internalization of Norms. J. Theor. Biol. 2003, 220, 407–418. [Google Scholar] [CrossRef] [PubMed]
  22. Nguyen, L.-T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  23. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice across a Large Model Space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lanfear, R.; Frandsen, P.B.; Wright, A.M.; Senfeld, T.; Calcott, B. PartitionFinder 2: New Methods for Selecting Partitioned Models of Evolution for Molecular and Morphological Phylogenetic Analyses. Mol. Biol. Evol. 2017, 34, 772–773. [Google Scholar] [CrossRef] [Green Version]
  25. Adeniran, A.A.; Hernández-Triana, L.M.; Ortega-Morales, A.I.; Garza-Hernández, J.A.; de la Cruz-Ramos, J.; Chan-Chable, R.J.; Vázquez-Marroquín, R.; Huerta-Jiménez, H.; Nikolova, N.I.; Fooks, A.R.; et al. Identification of Mosquitoes (Diptera: Culicidae) from Mexico State, Mexico Using Morphology and COI DNA Barcoding. Acta Trop. 2021, 213, 105730. [Google Scholar] [CrossRef]
  26. Dupanloup, I.; Schneider, S.; Excoffier, L. A Simulated Annealing Approach to Define the Genetic Structure of Populations. Mol. Ecol. 2002, 11, 2571–2581. [Google Scholar] [CrossRef]
  27. Mohd Yusoff, N.I.S.; Mat Jaafar, T.N.A.; Vilasri, V.; Mohd Nor, S.A.; Seah, Y.G.; Habib, A.; Wong, L.L.; Danish-Daniel, M.; Sung, Y.Y.; Mazlan, A.G.; et al. Genetic Diversity, Population Structure and Historical Demography of the Two-Spined Yellowtail Stargazer (Uranoscopus cognatus). Sci. Rep. 2021, 11, 13357. [Google Scholar] [CrossRef]
  28. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  29. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of Molecular Variance Inferred from Metric Distances among DNA Haplotypes: Application to Human Mitochondrial DNA Restriction Data. Genetics 1992, 131, 479–491. [Google Scholar] [CrossRef]
  30. Leigh, J.; Bryant, D. PopART: Full-Feature Software for Haplotype Network Construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  31. Tajima, F. DNA Polymorphism in a Subdivided Population: The Expected Number of Segregating Sites in the Two-Subpopulation Model. Genetics 1989, 123, 229–240. [Google Scholar] [CrossRef] [PubMed]
  32. Fu, Y.X. Statistical Tests of Neutrality of Mutations against Population Growth, Hitchhiking and Background Selection. Genetics 1997, 147, 915–925. [Google Scholar] [CrossRef]
  33. Rogers, A.R.; Harpending, H. Population Growth Makes Waves in the Distribution of Pairwise Genetic Differences. Mol. Biol. Evol. 1992, 9, 552–569. [Google Scholar] [CrossRef] [PubMed]
  34. Harpending, H.C. Signature of Ancient Population Growth in a Low-Resolution Mitochondrial DNA Mismatch Distribution. Hum. Biol. 1994, 66, 591–600. [Google Scholar] [CrossRef] [PubMed]
  35. Excoffier, L.; Lischer, H.E.L. 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–567. [Google Scholar] [CrossRef] [PubMed]
  36. Aungst, E.R.; Puente-Rolón, A.R.; Reynolds, R.G. Genetic Diversity in the US Ex Situ Populations of the Endangered Puerto Rican Boa, Chilabothrus inornatus. Zoo Biol. 2020, 39, 205–213. [Google Scholar] [CrossRef] [PubMed]
  37. Shi, X.-J.; Zhang, M.-L. Phylogeographical Structure Inferred from CpDNA Sequence Variation of Zygophyllum xanthoxylon across North-West China. J. Plant Res. 2015, 128, 269–282. [Google Scholar] [CrossRef] [PubMed]
  38. Ajao, A.M.; Nneji, L.M.; Adeola, A.C.; Oladipo, S.O.; Ayoola, A.O.; Wang, Y.-Y.; Adeniyi, A.V.; Olademeji, Y.U. Genetic Diversity and Population Structure of the Native Western African Honeybee (Apis mellifera Adansonii Latreille, 1804) in Nigeria Based on Mitochondrial COI Sequences. Zool. Anz. 2021, 293, 17–25. [Google Scholar] [CrossRef]
  39. Yan, J.; Wang, Q.; Chang, Q.; Ji, X.; Zhou, K. The Divergence of Two Independent Lineages of an Endemic Chinese Gecko, Gekko Swinhonis, Launched by the Qinling Orogenic Belt. Mol. Ecol. Mol. Ecol. 2010, 19, 2490–2500. [Google Scholar] [CrossRef]
  40. Zhou, K.; Wang, Q. New Species Of Gekko (Squamata: Sauria: Gekkonidae) From China: Morphological And Molecular Evidence. Zootaxa 2008, 1778, 59–68. [Google Scholar] [CrossRef]
  41. Ballenghien, M.; Faivre, N.; Galtier, N. Patterns of Cross-Contamination in a Multispecies Population Genomic Project: Detection, Quantification, Impact, and Solutions. BMC Biol. 2017, 15, 25. [Google Scholar] [CrossRef] [PubMed]
  42. Ellegren, H.; Galtier, N. Determinants of Genetic Diversity. Nat. Rev. Genet. 2016, 17, 422–433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Garza, J.C.; Williamson, E.G. Detection of Reduction in Population Size Using Data from Microsatellite Loci. Mol. Ecol. 2001, 10, 305–318. [Google Scholar] [CrossRef] [PubMed]
  44. Whitehouse, A.M.; Harley, E.H. Post-Bottleneck Genetic Diversity of Elephant Populations in South Africa, Revealed Using Microsatellite Analysis. Mol. Ecol. 2002, 10, 2139–2149. [Google Scholar] [CrossRef] [PubMed]
  45. Wisely, S.M.; Buskirk, S.W.; Fleming, M.A.; McDonald, D.B.; Ostrander, E.A. Genetic Diversity and Fitness in Black-Footed Ferrets before and during a Bottleneck. J. Hered. 2002, 93, 231–237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Ploshnitsa, A.; Goltsman, M.; Happ, G.; Macdonald, D.; Kennedy, L.; Volff, J. Historical and Modern Neutral Genetic Variability in Mednyi Arctic Foxes Passed through a Severe Bottleneck. J. Zool. 2013, 289, 68–76. [Google Scholar] [CrossRef]
  47. Newbery, B.; Jones, D. Presence of Asian House Gecko Hemidactylus Frenatus across an Urban Gradient in Brisbane: Influence of Habitat and Potential for Impact on Native Gecko Species. In Pest or Guest: The Zoology of Overabundance; Lunney, D., Eby, P., Hutchings, P., Burgin, S., Eds.; Royal Zoological Society of New South Wales: Mosman, NSW, Australia, 2007; pp. 59–65. [Google Scholar] [CrossRef] [Green Version]
  48. Yang, D.; González-Bernal, E.; Greenlees, M.; Shine, R. Interactions between Native and Invasive Gecko Lizards in Tropical Australia. Austral Ecol. 2012, 37, 592–599. [Google Scholar] [CrossRef]
  49. Werneck, F.P.; Leite, R.N.; Geurgas, S.R.; Rodrigues, M.T. Biogeographic History and Cryptic Diversity of Saxicolous Tropiduridae Lizards Endemic to the Semiarid Caatinga. BMC Evol. Biol. 2015, 15, 94. [Google Scholar] [CrossRef] [Green Version]
  50. Excoffier, L.; Foll, M.; Petit, R.J. Genetic Consequences of Range Expansions. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 481–501. [Google Scholar] [CrossRef]
Figure 1. Sampling sites of G. japonicus. Different colors represent different SAMOVA grouping results The blue circles represent the sampling points of Group A, red circles represent the sampling points of Group B, yellow circles represent the sampling points of Group C, green circles represent the sample sampling points of Group D. Group A: Yongan, Fujian, Yangshuo, Guangxi, Huaihua, Hunan; Group B: Yangxian, Shaanxi. Group C: Anhui (Anqing, Lu’an, Wuhu), Fujian Nanping, Guangzhou Guiding, Guizhou (Huaxi, Longli, Libo), Guangxi (Guilin, Longsheng), Hubei (Jingmen, Wuhan), Hunan (Chengbu, Changde, Daoxian, Huayuan, Huayuan, Shuangpai, Shaoyang, Tong, Xinhuang, Yueyang, Yongzhou, Xinning, Yongzhou Lengshuitan, Zhuzhou, Yongzhou Qingtang Yueyyan Forest Farm), Jiangsu Rugao, Jiangxi Longnan, Zhejiang (Yongzhou Qingtang Yueyyan Forest Farm). Zhoushan, Lishui). Group D: Yongfu, Guangxi.
Figure 1. Sampling sites of G. japonicus. Different colors represent different SAMOVA grouping results The blue circles represent the sampling points of Group A, red circles represent the sampling points of Group B, yellow circles represent the sampling points of Group C, green circles represent the sample sampling points of Group D. Group A: Yongan, Fujian, Yangshuo, Guangxi, Huaihua, Hunan; Group B: Yangxian, Shaanxi. Group C: Anhui (Anqing, Lu’an, Wuhu), Fujian Nanping, Guangzhou Guiding, Guizhou (Huaxi, Longli, Libo), Guangxi (Guilin, Longsheng), Hubei (Jingmen, Wuhan), Hunan (Chengbu, Changde, Daoxian, Huayuan, Huayuan, Shuangpai, Shaoyang, Tong, Xinhuang, Yueyang, Yongzhou, Xinning, Yongzhou Lengshuitan, Zhuzhou, Yongzhou Qingtang Yueyyan Forest Farm), Jiangsu Rugao, Jiangxi Longnan, Zhejiang (Yongzhou Qingtang Yueyyan Forest Farm). Zhoushan, Lishui). Group D: Yongfu, Guangxi.
Genes 14 00018 g001
Figure 2. BI and ML tree constructed based on the COI gene sequence of 11 species and AHA06 of Gekko from China. Members of the G. japonicus group species are shown in orange, Gekko gecko group species are shown in green, Gekko monarchus group species are shown in blue and the purple part is the outgroup. The numbers above branches specify the posterior probability and bootstrap percentages (1000 replicates).
Figure 2. BI and ML tree constructed based on the COI gene sequence of 11 species and AHA06 of Gekko from China. Members of the G. japonicus group species are shown in orange, Gekko gecko group species are shown in green, Gekko monarchus group species are shown in blue and the purple part is the outgroup. The numbers above branches specify the posterior probability and bootstrap percentages (1000 replicates).
Genes 14 00018 g002
Figure 3. The relationship diagram of the haplotype network of the G. japonicus based on the COI gene. The circles indicate the haplotypes; circle sizes indicate the probability of haplotypes. The different colors in the circles indicate the distribution in different populations, and the oblique lines indicate mutations between haplotypes.
Figure 3. The relationship diagram of the haplotype network of the G. japonicus based on the COI gene. The circles indicate the haplotypes; circle sizes indicate the probability of haplotypes. The different colors in the circles indicate the distribution in different populations, and the oblique lines indicate mutations between haplotypes.
Genes 14 00018 g003
Table 1. Analysis of genetic distance within and between species of Chinese Gekko.
Table 1. Analysis of genetic distance within and between species of Chinese Gekko.
Species1234567891011
1. G. japonicus0.6 ± 0.2
2. AHAQ0612.0–12.9-
3. G.wenxianensis21.7–23.523.6–24.11.0 ± 0.3
4. G. chinensis31.7–33.328.3–31.025.4–28.73.4 ± 0.6
5. G. swinhonis30.0–34.427.2–31.629.2–32.825.1–30.09.1 ± 0.9
6. G. subpalmatus31.6–33.032.1–33.028.4–30.629.1–30.222.5–25.81.3 ± 0.4
7. G.auriverrucosus30.2–31.430.5–31.331.5–32.730.2–33.024.0–27.224.5–26.20.3 ± 0.2
8. G. hokouensis33.5–35.429.8–31.328.7–30.830.6–33.728.9–31.427.9–30.130.5–33.10.9 ± 0.3
9. G. scabridus31.4–35.233.8–36.226.8–30.724.1–29.323.2–29.525.3–29.027.3–31.729.2–33.04.3 ± 0.6
10. G.gecko36.2–38.635.1–35.835.4–37.336.4–38.534.2–37.842.0–44.939.1–40.737.0–40.835.3–38.52.1 ± 0.4
11. G.kikuchii33.0–35.235.537.3–38.128.7–30.428.2–33.628.1–28.236.1–37.336.2–37.828.5–29.732.6–33.4-
Note: The table below the diagonal line indicates the interspecific genetic distance, the bold numbers represent the average genetic distance within each species, and the genetic distance is expressed as a percentage.
Table 2. Genetic variation parameters of 37 various geographic populations.
Table 2. Genetic variation parameters of 37 various geographic populations.
PopulationNumber of SamplesSegregating Site (S)HaplotypesHaplotype Diversity (h)Nucleotide Diversity (π)Average Number of Difference (K)
Anqing, Anhui (AHAQ)112Hap1, 2, 30.4730.000930.618
Lin’an, Anhui (AHLA)10Hap1NANANA
Wuhu, Anhui (AHWH)122Hap1, 40.1670.000500.333
Yong’an, Fujian (FJYA)33Hap1, 5NANANA
Nanping, Fujian (FJNP)10Hap1NANANA
Guilin, Guangxi (GXGL)193Hap1, 50.1050.000470.316
Longsheng, Guangxi (GXLS)141Hap1, 60.2640.000400.264
Yongfu, Guangxi (GXYF)53Hap1, 50.4000.001801.2
Guiding, Guangzhou (GZGD)20Hap1NANANA
Libo, Guizhou (GZLB)30Hap1NANANA
Huaxi, Guizhou (GZHX)10Hap1NANANA
Longli, Guizhou (GZLL)20Hap1NANANA
Wuhan, Hubei (HBWH)133Hap1, 4, 50.5900.001971.308
Hubeijing (HBJM)10Hap1NANANA
Chengbu, Hunan (HNCB)130Hap1000
Changde, Hunan (HNCD)82Hap1, 70.2500.000750.5
Daoxian, Hunan (HNDX)194Hap1, 5, 60.5560.002371.579
Huaihua, Hunan (HNHH)153Hap1, 50.2480.001120.743
Huitong, Hunan (HNHT)162Hap1, 8, 90.5420.000910.583
Huayuan, Hunan (HNHY)163Hap1, 50.2330.001050.7
Shuangpai, Hunan (HNSP)44Hap1, 5, 10NANANA
Shaoyang, Hunan (HNSY)51Hap1, 110.4000.000600.4
Tongdao, Hunan (HNTD)62Hap1, 40.3330.001000.667
Xinkou, Hunan (HNXH)203Hap1, 50.5210.002351.563
Xinning, Hunan (HNXN)40Hap1NANANA
Yueyang, Hunan (HNYY)120Hap1000
Yongzhou, Hunan (HNYZ)123Hap1, 4, 50.5910.002321.545
Zhuzhou, Hunan (HNZZ)74Hap1, 5, 60.5240.001721.143
Rugao, Jiangsu (JSRG)20Hap1NANANA
Longnan, Jiangxi (JXLN)246Hap1, 120.1590.001440.957
Yueyan Forestry, Qingtang, Yongzhou, Hunan (QTYY)53Hap1, 50.4000.001881.200
Yangxian, Shaanxi (SXYX)153Hap1, 50.5330.002411.600
Lingshutan, Yongzhou, Hunan (YZLST)122Hap1, 40.3030.000910.606
Hangzhou, Zhejiang (ZJHZ)85Hap1, 5, 130.4640.001881.25
Zhoushan, Zhejiang (ZJZS)114Hap1, 4, 140.5640.002081.382
Lishui, Zhejiang (ZJLS)10Hap1NANANA
Yangshuo, Guangxi (GXYS)10Hap1NANANA
Overall32416140.4530.001761.171
Table 3. AMOVA analysis of the G. japonicus population.
Table 3. AMOVA analysis of the G. japonicus population.
Source of
Variations
Sum of
Squares
Variance
Components
Percentage of VariationFixation IndicesSignificance Tests
Among groups47.7620.5742039.53FCT = 0.395270.00000
Among populations within groups32.0260.012140.84FSC = 0.013820.02639
Within populations252.1080.8663559.64
Total333.8961.45269 FST = 0.403620.03030
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xia, L.; Cai, F.; Chen, S.; Cai, Y.; Zhou, K.; Yan, J.; Li, P. Phylogenetic Analysis and Genetic Structure of Schlegel’s Japanese Gecko (Gekko japonicus) from China Based on Mitochondrial DNA Sequences. Genes 2023, 14, 18. https://doi.org/10.3390/genes14010018

AMA Style

Xia L, Cai F, Chen S, Cai Y, Zhou K, Yan J, Li P. Phylogenetic Analysis and Genetic Structure of Schlegel’s Japanese Gecko (Gekko japonicus) from China Based on Mitochondrial DNA Sequences. Genes. 2023; 14(1):18. https://doi.org/10.3390/genes14010018

Chicago/Turabian Style

Xia, Longjie, Fengna Cai, Shasha Chen, Yao Cai, Kaiya Zhou, Jie Yan, and Peng Li. 2023. "Phylogenetic Analysis and Genetic Structure of Schlegel’s Japanese Gecko (Gekko japonicus) from China Based on Mitochondrial DNA Sequences" Genes 14, no. 1: 18. https://doi.org/10.3390/genes14010018

APA Style

Xia, L., Cai, F., Chen, S., Cai, Y., Zhou, K., Yan, J., & Li, P. (2023). Phylogenetic Analysis and Genetic Structure of Schlegel’s Japanese Gecko (Gekko japonicus) from China Based on Mitochondrial DNA Sequences. Genes, 14(1), 18. https://doi.org/10.3390/genes14010018

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop