• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Description of two new species of Hemiphyllodactylus(Reptilia: Gekkonidae) from karst landscapes in Yunnan, China, highlights complex conservation needs

    2022-10-17 03:27:58AdePrasetyoAgungAdaChorneliaLeeGrismerJesseGrismerEvanQuahJianMeiLuKyleTomlinsonAliceHughes
    Zoological Research 2022年5期
    關(guān)鍵詞:思茅

    Ade Prasetyo Agung, Ada Chornelia, L. Lee Grismer, Jesse L. Grismer, Evan S. H. Quah, Jian-Mei Lu,Kyle W. Tomlinson, Alice C. Hughes

    1 Landscape Ecology Group, Center for Integrative Conservation, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Mengla, Yunnan 666303, China

    2 University of Chinese Academy of Sciences, Beijing 101408, China

    3 Herpetology Laboratory, Department of Biology, La Sierra University, Riverside, California 92515, USA

    4 Institute for Tropical Biology and Conservation, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu, Sabah 88400, Malaysia

    5 Lee Kong Chian Natural History Museum, National University of Singapore, Singapore 117377, Singapore

    6 Community Ecology and Conservation Group, Center for Integrative Conservation, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Mengla, Yunnan 666303, China

    7 School of Biological Sciences, University of Hong Kong, Pokfulam, Hong Kong SAR, China

    ABSTRACT Karst habitats are hotspots of diversity and endemism. Their naturally fragmented distributions across broad geographic landscapes have led to the complex array of smaller evolutionary ecosystems that present unique challenges from a conservation perspective. Comprehensive biodiversity assessments of karst habitats have revealed that these ecosystems contain an almost unparalleled level of endemism, and many site-restricted species remain undescribed, thus posing considerable challenges for effective conservation management.Small rock-dwelling species, such as geckos, may be particularly prone to such isolation. In this paper, we discuss one such genus, i.e.,Hemiphyllodactylus, and explore its diversity across karst landforms in Yunnan Province, southwestern China. Based on morphological and genetic data, we describe two new species of Hemiphyllodactylus from karst habitats in Simao District and Yanshan County. A phylogenetic tree for Hemiphyllodactylus was constructed using 1 039 base pairs (bp) of the mitochondrial NADH dehydrogenase subunit 2 gene(ND2). The Simao and Yanshan specimens can be distinguished from all other congeners within their respective subclades based on uncorrected genetic pairwise distances greater than 6.3% and 4.3%respectively, as well as significant morphological differences. The discovery and description of these two new species brings the total number of described Hemiphyllodactylus species in China to 14 and indicates many more undescribed species from unsurveyed karst regions await discovery. Our findings suggest that karst ecosystems in Yunnan support a higher diversity of Hemiphyllodactylus than previously known. This study also highlights the importance of karst ecosystems as refugia for sitespecific endemic species and the need for heightened conservation efforts.

    Keywords: Discovery; Endemism; Geckos;Reptiles; South China

    lNTRODUCTlON

    Karst landscapes are characterized by high endemism due to their distinct ecological niches, which allow for the diversification of a wide variety of species (Clements et al.,2006; Grismer et al., 2021). However, these naturally fragmented ecosystems are challenging, as each lone karst formation may host species found nowhere else. Almost every research expedition of karst landforms (both caves and outer rocky surfaces) has uncovered new species with localized distributions (Agung et al., 2021; Dittmar et al., 2005; Huang et al., 2019; Quah et al., 2021; Tian & Huang, 2015), thereby identifying karsts as hotspots of endemism and biodiversity and priorities for conservation. Karst species are also often highly specialized, with poor dispersal capabilities due to their adaptations to the unique abiotic environments of karst ecosystems, such as microclimate, light intensity, and topography (e.g., fissured cliffs) (Whitten, 2009). Thus,understanding the distributions and patterns of endemism and diversity is critical for developing effective conservation plans for the region as many cryptic herptiles may have been overlooked (Vieites et al., 2009).

    As the largest family of geckos, Gekkonidae shows high levels of endemism in karst systems. For example, Grismer et al. (2014, 2018b, 2018c, 2020a, 2021) has reported over 100 gecko species endemic to karsts in Southeast Asia,underscoring the need for further work to gain a more complete understanding of the diversity and range of karstdependent taxa. Given their occurrence in fragmented karst hills and their limited dispersal capabilities, it is likely that more species are waiting to be described. In support of this, Grismer et al. (2018a) identified 12 new gecko species within two weeks in a single study of karsts in Myanmar, with similar patterns likely to exist in karsts across Southeast Asia.

    HemiphyllodactylusBleeker, 1860 (commonly known as half leaf-fingered geckos, dwarf geckos, or slender geckos),belongs to the family Gekkonidae. Recently, many new species of this genus have been discovered in vegetated karst ecosystems (Do et al., 2020; Grismer et al., 2018b; Nguyen et al., 2020; Zhang et al., 2020), many of which are endemic to single karst hills. However, like other small organisms (e.g.,snails, millipedes, and other invertebrates), this genus is often overlooked, and given the highly endemic nature of the group,more targeted protection is clearly needed.

    The genus is widely distributed across South Asia,Southeast Asia, South China, and the western Pacific islands(Agarwal et al., 2019; Grismer et al., 2013; Zug, 2010).Generally,Hemiphyllodactylusspecies are small in body size(snout-vent length<63 mm), nocturnal, scansorial, forestdwelling, and well camouflaged in their environments (Grismer et al., 2013; Zug, 2010), and thus easily overlooked unless specifically targeted. Most members of the genus are confined to tropical and subtropical montane regions in mainland Indochina, although some are also restricted to islands(Agarwal et al., 2019; Eliades et al., 2019; Grismer et al.,2013; Zug, 2010). To date,Hemiphyllodactylusconsists of two main groups, i.e.,hartertiandtypus. Thehartertigroup is composed of upland species from Peninsular Malaysia, while thetypusgroup is comprised of all other species from the entire range of the genus (Grismer et al., 2013). Many of these species are limited to a single site or a limited number of sites,and some “species” are suspected to represent complexes that require further work for accurate species delineation and description. This is especially important given the high rates of karst loss across Southeast Asia and South China, estimated to be 5.7% per year (Hughes, 2017). Thus, without information to ensure the identification of key sites, there is significant potential for species loss (Hughes, 2017).

    In recent years,Hemiphyllodactylusresearch and discoveries have experienced a renaissance, with species descriptions increasing every year from 2013. Since then, the number of species has jumped from 14 to 52, mostly from karst regions in Myanmar, Laos, and Vietnam (Uetz et al.,2021). In contrast, gekkonid research in karst regions of South China continues to lag, highlighting the need for in-depth exploration of these areas. This deficiency in field research has resulted in an underestimation ofHemiphyllodactylusdiversity in China. For example, in Yunnan Province, all populations ofHemiphyllodactyluswere previously considered to be a single widespread species (H. yunnanensis), until integrative taxonomic study revealed multiple species under the nomenH. yunnanensis(Grismer et al., 2013).

    At present, 12 described species ofHemiphyllodactylusare found in South China: i.e.,H. changningensisGuo, Zhou, Yan& Li, 2015;H. dupanglingensisZhang, Qian & Yang, 2020;H.dushanensisZhou & Liu, 1981;H. hongkongensisSung, Lee,Ng, Zhang & Yang, 2018;H. huishuiensisYan, Lin, Guo, Li &Zhou, 2016;H. jinpingensisZhou & Liu, 1981;H. longlingensisZhou & Liu, 1981;H. typusBleeker, 1860;H. yunnanensisBoulenger, 1903;H. zayuensisJiang, Wang, & Che, 2020,H.zhutangxiangensisAgung, Grismer, Grismer, Quah,Chornelia, Lu & Hughes, 2021; andH. zugiNguyen,Lehmann, Le Duc, Duong, Bonkowski & Ziegler, 2013.

    Based on fieldwork in the karst areas of Yunnan, we identified several potential new species, including the recently described and published speciesH. zhutangxiangensis(Agung et al., 2021). Here, we describe two new species, with our preliminary work also suggesting several other species awaiting description. In the current study, we phylogenetically delimited new evolutionary lineages based on molecular evidence and diagnosed those lineages based on morphological evidence, with descriptions of the new species.We also discuss the relevance of mapping distributions such as these as a basis for conservation of range-limited, siteendemic species that may be particularly vulnerable to extinction without better inclusion in conservation and environmental impact assessments.

    MATERlALS AND METHODS

    Field sites and specimen collection

    Surveys were conducted in August to September 2018 and June to July 2019 in karst landscapes across Yunnan.Permission to conduct surveys and collect samples was granted by relevant local protected area authorities and ethics approval for research was granted by the Ethics Committee of Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, China. In total, 139Hemiphyllodactylusspecimens were collected by hand from 13 different karst areas (Table 1). Specimens were collected at night (between 2000h-2200h) as geckos are nocturnal and actively forage during this time. AllHemiphyllodactylusspecimens were photographed in the morning after capture to record coloration and patterns. The specimens were euthanized using Tricaine MS-222 solution injected into the intracelomic cavity (Conroy et al., 2009). Tissue samples were obtained from the liver of each individual and stored in 95% ethanol separately for further genetic analysis. The specimens were then preserved in 10% formalin and transferred to 70% ethanol for storage prior to morphological investigation. After investigation, the specimens were deposited in the Kunming Natural History Museum of Zoology, Kunming Institute of Zoology (KIZ),Chinese Academy of Sciences, China. All animal procedures were performed in accordance with the ethical standards of the institution at which the study was conducted(Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences).

    DNA isolation, sequencing, and phylogenetic analyses

    Genomic DNA of the 139 newly collected specimens was isolated from liver tissue following proteinase K DNA extraction protocols using a QIAGEN Genomic-tip 2500(www.qiagen.com) DNA extraction kit. We amplified the complete mitochondrial NADH dehydrogenase subunit 2 gene(mtDNA-ND2), totaling 1 039 bp, using primers L4437b and H5934 following Macey et al. (1997) (5'-AAGCAGT TGGGCCCATACC-3' and 5'-AGRGTGCCAATGTCTTTG TGRTT-3', respectively). The protocol for polymerase chain reaction (PCR) amplifications followed Agung et al. (2021).The PCR processes and sequencing were executed at the South China DNA Barcoding Center.

    We constructed a dataset for phylogenetic analyses. We downloaded a total of 209ND2sequences from GenBank,containing 205 sequences of extantHemiphyllodactylusspecies and fourND2sequences of other outgroup taxa(Gehyra fehlmanni,G. mutilata,Hemidactylus frenatusandLepidodactylus lugubris), then added the 139 new sequences to the dataset. All downloaded sequences used in the analyses followed Agung et al. (2021), while the newly published sequences in this study are presented in Supplementary Table S1.

    Phylogenetic relationships were analyzed using maximumlikelihood (ML) and Bayesian inference (BI) in IQ-TREE(Nguyen et al., 2015; Trifinopoulos et al., 2016) and MrBayes 3.2.7a (Ronquist et al., 2012) on XSEDE using the CIPRES Science Gateway (Cyberinfrastructure for Phylogenetic Research; Miller et al., 2010), respectively. Prior to ML analysis, the best substitution model (TIM+F+R5) was selected for the non-partitioned dataset based on Bayesian information criterion (BIC) in ModelFinder (Kalyaanamoorthy et al., 2017). The ultrafast bootstrap approximation algorithm(UFBoot) was used with 1 000 bootstrap pseudoreplicates(Hoang et al., 2018), where nodes bearing values ≥95 were considered strongly supported (Minh et al., 2013). For BI analysis, default priors were selected, and two independent Markov Chain Monte Carlo (MCMC) algorithms were applied,with four chains in each (three hot and one cold), 50 million generations sampled every 1 000 generations, and the first 25% of samples discarded. All parameters from the two runs were checked in Tracer v1.7.1 (Rambaut et al., 2018),confirming convergences and effective sample sizes (ESS)were >200. Post burn-in sampled trees from both runs were combined and a 50% majority-rule consensus tree was produced. Nodes bearing Bayesian posterior probabilities(BPP) ≥0.95 were considered strongly supported(Huelsenbeck et al., 2001; Wilcox et al., 2002). Uncorrectedpairwise distances among and within species were computed using MEGA 7 (Kumar et al., 2016).

    Table 1 Location of the 13 selected karst field sites in Yunnan, China

    Morphological measurements and analyses

    We measured mensural and meristic traits of the 139 collected specimens following Zug (2010) and Grismer et al. (2013),with slight modifications following Agung et al. (2021).Mitutoyo Absolute Series-500 digital calipers (accuracy 0.01 mm) were used to measure mensural traits under a dissecting microscope (Nikon SMZ 445), on the left side of the body when possible. Recorded traits included: snout-vent length(SVL), trunk length (TL), head length (HL), head width (HW),eye diameter (ED), snout-eye length (SnEye), naris-eye length(NarEye), and snout width (SnW). Recorded meristic traits included: circumnasal scales (CN), internasal scales (IS),supralabial scales (SL), infralabial scales (IL), chin scales(Chin), ventral scales (VS), dorsal scales (DS), subdigital lamellae wider than long on first finger (SL1F) and toe (SL1T),subdigital lamellae formula determined as number of Ushaped digital pads on digits II-V of hands and feet, number of femoroprecloacal pores, and number of cloacal spurs(CloacS) on each side of hemipenial swelling. We also noted coloration and pattern on the dorsum, presence or absence of dark postorbital stripes extending at least to neck, presence or absence of dorsolateral and ventrolateral stripes, and presence or absence of anteriorly projecting arms of postsacral markings.

    We compared the morphology of each new lineage recovered from phylogenetic analysis against published morphological data for selected closely related species to establish significant differences between any of the measured traits. Prior to analysis, we corrected for the effects of body size on mensural traits in each new lineage combined with data from closely related species using the following equation:

    whereXadj is the adjusted value;Xis the measured value; β is the unstandardized regression coefficient for each operational taxonomic unit (OTU); SVL is the measured SVL;and SVLmean is the overall average SVL of all OTUs(Lleonart et al., 2000; Thorpe, 1975, 1983; Turan, 1999).

    All morphological analyses were computed in R v.4.0.1 (R Core Team, 2020). Each adjusted mensural trait was checked for equal variances across groups using Levene’s test. Traits with equal variances (P≥0.05) were analyzed using ordinary linear models, while traits with unequal variances (P≤0.05)were analyzed using generalized linear models with weighted least squares (gls) in the "nlme" package (Pinheiro et al.,2020). Error degrees of freedom (df) were calculated with the Satterthwaite approximation, applying theemmeansfunction in the package “emmeans” (Lenth, 2021), which compares appropriate estimates for uneven group variances. For meristic traits (count data), glm models were used. Quasi-Poisson errors were implemented in the glm models as all traits were under-dispersed when checked with thedispersiontestfunction in the “AER” package (Kleiber &Zeileis, 2008). For both mensural and meristic traits,significant differences were first evaluated with a variance test(analysis of variance (ANOVA) or Chi-square tests), and, if significant, subjected toposthocTukey’s HSD tests for mean comparisons involving three or more groups. We also used non-parametric permutation ANOVA (PERMANOVA) to determine whether the posited species differed from closely related species in multi-trait space using theadonisandpairwise.adonisfunctions in the “vegan” package (Oksanen et al., 2020) with 50 000 permutations. Lastly, principal component analysis (PCA) was run to test for group separation along the first two principal components using theprcompfunction.

    Designation of species-level lineages

    We used an integrated approach to delineate species-level lineages by consolidating phylogenetic position, genetic divergence, and morphological differences. We used three sequential criteria to designate species-level lineages: First,new lineages not clustered within the named species lineages in phylogenetic topology were marked as potential new species. Second, uncorrected pairwise genetic distances among the new lineages with either known species or other putative species lineages were measured, with a 3.0%difference in mtDNAND2considered minimal to define a potential new species. This cut-off was based on Zhang et al.(2020), who reported a 3.6% sequence divergence between

    H. linnwayensisandH. ywanganensisafter uncorrected pairwise comparison of the 670 bp (partial)ND2gene across numerous species, with these two species known to be morphologically distinct (Grismer et al., 2018c). Third, those lineages showing >3.0% genetic difference from their nearest relatives were examined for morphological distinctiveness from closely related species. If this third criterion was also met, the lineage was considered a confirmed new species.Lineages satisfying all three criteria were described.Individuals not satisfying all criteria may also be distinct species, but further data and specimens are needed for this to be established.

    RESULTS

    Phylogeny

    The ML and BI analyses recovered similar tree topologies with strong nodal support, but with a slight difference in the position ofH. khlonglanensis, i.e., sister species to clades 3 and 4 in the ML tree, but sister species to clade 4 only in the BI tree.Figure 1 shows the topology of the ML tree, including both nodal support values (UFBoot/BPP). Our inferred ML topology was highly consistent with the topology reported by Grismer et al. (2020b). We named all clades following Grismer et al.(2020b) and renamed the South Myanmar lineage, which consisted ofH. zwegabinensis,H. kyaiktiyoensis, andH.pinlaungensis, as clade 8.

    Figure 1 Maximum-likelihood topology illustrating phylogenetic relationships among Hemiphyllodactylus species

    All 139 Yunnan karst specimens belonged to thetypusgroup and were placed into four clades: i.e., clades 3, 4, 6,and 7 (in bold in Figure 1). All specimens were nested and formed OTUs as follows: 25 specimens were from Yongde County, recovered asH. changningensis(clade 3); 47 specimens were from Jian Shui County, Yimen County, and Gejiu City, recovered asH. jinpingensisand sister to the newly recovered OTU comprised of 22 specimens from Simao District and Ning’er County (clade 4); 13 specimens were from Yanshan County, recovered as a new OTU (clade 6) and sister toH. huishuiensis; and 32 specimens were from multiple locations, including Mengla County, Simao District, Luxi County, Jianshui County, Gejiu City, and Tonghai County, and nested within theyunnanensisclade complex (clade 7). Theyunnanensisclade (clade 7) consisted of multiple distinct evolutionary lineages, and further analysis is needed to resolve their relationships and taxonomy.

    The ML and BI trees recovered two distinct evolutionary lineages in our Yunnan samples, i.e., Simao and Ning’er County population (hereafter referred to as Simao population)in clade 4 and Yanshan County population in clade 6, as described below.

    Genetic distance

    Uncorrected geneticP-distances among and within theND2gene of the newly recovered OTUs (within clades 4 and 6) are presented in Tables 2, 3. The interspecific genetic distances within clade 4 ranged from 2.7% (betweenH. ywanganensisandH. linnwayensis) to 14.3% (betweenH. chiangmaiensisandH. jinpingensis), and the new Simao population showed at least 6.3% genetic distance to its sister speciesH. jinpingensis(Table 2).

    In clade 6, interspecific genetic distances ranged from 4.3%(between new Yanshan population and sister speciesH.huishuiensis) to 23.0% (betweenH. kizirianiandH. banaensis;Table 3). The genetic distances of the two new populations(Simao and Yanshan) within their clades (clades 4 and 6,respectively) were above 3%, the threshold used to delineate a new species in this study (see Methods).

    Morphological analysis of Simao population

    We compared the Simao population to members within clade 4, i.e.,H. jinpingensis,H. montawaensis,H. ngwelwini, andH.tonywhitteni, as they were the closest relatives to the Simao population containing more than two individuals in each species (P-distance=6.3%, 10.2%, 10.5%, and 10.8%,respectively; Table 2). The rawH. jinpingensisdata used for statistical analysis were obtained from our specimens (45 individuals out of 47 collected from various locations in Yunnan) as no appropriate raw data were available from previous studies. The rawH. jinpingensisdata are provided in Supplementary Table S2. The raw data forH. montawaensisandH. tonywhitteniwere obtained from Grismer et al. (2018b),and forH. ngwelwiniwere obtained from Grismer et al.(2020a). Linear and general linear model analyses of each morphological trait in the Simao population and its congeners showed that seven mensural traits and eight meristic traits differed significantly among the five groups (Table 4).

    Post hocmultiple comparison tests showed significant differences in two mensural and five meristic traits between theH. jinpingensisand Simao specimens, although it should be noted that there was overlap in the recorded range of values among individuals across species (see Supplementary Table S3). Of the two mensural traits,H. jinpingensishad greater ED, whereas the Simao specimens had greater HW.Furthermore,H. jinpingensishad higher values in five meristic traits (VS, SL, IL, SL1F, and SL1T) compared to the Simaospecimens. Two mensural and five meristic traits differed significantly between theH. montawaensisand Simao specimens. Of the two mensural traits,H. montawaensishad greater HL, whereas Simao specimens had greater HW. Of the five meristic traits,H. montawaensishad more VS, but fewer SL, IL, SL1T, and Chin. One mensural and three meristic traits differed significantly between the Simao specimens andH. ngwelwini, notably the Simao specimens had shorter HL and fewer VS and Chin, but more IL. Two mensural and seven meristic traits differed significantly between the Simao specimens andH. tonywhitteni.Of the two mensural traits, the Simao specimens had a shorter HL, but greater HW. For the meristic traits, the Simao specimens had more CN, SL, IL, SL1F, and Chin, but fewer DS and VS(Figure 2).

    Table 2 Uncorrected genetic P-distances (%) in ND2 (1 039 bp) gene of Hemiphyllodactylus in clade 4

    Table 3 Uncorrected genetic P-distances (%) in ND2 (1 039 bp) gene of Hemiphyllodactylus in clade 6

    Table 4 Summary statistics for each trait of compared species in clade 4

    Multi-trait spatial analysis using PERMANOVA pairwise comparisons indicated that the five groups (Simao specimens,

    H. jinpingensis,H. montawaensis,H. ngwelwini, andH.tonywhitteni) differed significantly from each other (Table 5). In addition, clustering based on PCA showed that the Simao specimens were morpho-spatially distinct and clustered separately along PC1 with respect toH. montawaensisandH.tonywhitteni, but partially overlapped with respect toH.jinpingensisandH. ngwelwini, with 26.7% of the variation loaded most heavily for HL and IL (see Supplementary Table S4 for PCA scores on each trait). For PC2, all groups overlapped with each other, except forH. montawaensisandH. tonywhitteni(Figure 3). Details on statistical results for group separation along the two axes are presented in Table 6.

    Morphological analyses of Yanshan population

    We compared the Yanshan population toH. huishuiensis(raw data obtained from Yan et al., 2016) andH. nahangensis(raw data obtained from Do et al., 2020), as the latter two species were the closest relatives of the Yanshan population based on phylogeny and uncorrected genetic distance (Pdistance=4.3% and 8.1%, respectively; Table 3). Linear and general linear model analyses of each morphological trait in the Yanshan population and its two closely related species identified significant differences in seven mensural and seven meristic traits (Table 7).

    Post hocmultiple comparison tests showed that six mensural and five meristic traits differed significantly between theH. huishuiensisand Yanshan specimens, although it should be noted that there was overlap in the recorded range of values among individuals across species (see Supplementary Table S5). Of the six mensural traits,H.huishuiensishad greater TrunkL and HL, whereas Yanshan specimens had greater HW, NarEye, SnEye, and SnW. Of the five meristic traits,H. huishuiensishad more IS, whereas Yanshan specimens had more CN, CloacS, SL1F, and SL1T.Seven mensural and five meristic traits differed significantly between the Yanshan specimens andH. nahangensis. Of the seven mensural traits, Yanshan specimens had greater HW,NarEye, SnEye, and SnW, whereasH. nahangensishad greater TrunkL, ED, and HL. Of the five meristic traits,Yanshan specimens had more CN, SL1F, and SL1T, whereasH. nahangensishad more DS and IS (Figure 4).

    Figure 2 Differences in adjusted mensural (left) and meristic (right) traits between Simao population and closely related species

    Table 5 PERMANOVA pairwise comparisons of multi-trait space of compared species within clade 4

    Multi-trait spatial analysis using PERMANOVA pairwise comparison indicated that the three groups (Yanshan population,H. huishuiensis, andH. nahangensis) differed significantly from each other (Table 8). In addition, clustering based on PCA showed that specimens from Yanshan,H.huishuiensis, andH. nahangensiswere morpho-spatially distinct and clustered separately along PC1, which accounted for 45.7% of variation in the dataset and was loaded most heavily for SL1T, CN, and SnW (see Supplementary Table S6 for PCA scores on each trait). The Yanshan specimens overlapped withH. huishuiensisandH. nahangensisalong PC2, which accounted for 15.5% of variation in the dataset and was loaded most heavily for DS and VS (Figure 5).Details on statistical results for group separation on the two axes are presented in Table 9.

    Systematics

    Based on phylogenetic evidence, uncorrected genetic distance, and lineage distinction from significantly different morphological characters, the Simao and Yanshan populations likely constitute new evolutionary lineages, and thus represent two new species, as described below.

    Hemiphyllodactylus simaoensis sp. nov. (Figure 6)

    Suggested English name: Simao slender gecko Suggested Chinese name: 思茅半葉趾虎

    Figure 3 Principal component analysis (left) and biplot (right) of compared species within clade 4

    Table 6 Test results for group separation of compared species within clade 4 along first two PCA axes (PC1 and PC2)

    Holotype: Adult male (KIZ 062064) collected on 6 July 2019 by Ade P. Agung, Jian-Mei Lu, and Zong-Bao Yang from forested karst hills in Simao District, Pu’er City, Yunnan Province, China (N22.735°, E100.803°; 1 129 m a.s.l.).

    Paratypes:Four adult males (KIZ 062063, KIZ 062066-KIZ 062068) and three adult females (KIZ 062065, KIZ 062069-KIZ 062070), same data as holotype. Five adult males (KIZ 062072, KIZ 062073, KIZ 062077, KIZ 062079,KIZ 062087) and nine adult females (KIZ 062074-KIZ 062076,KIZ 062080-KIZ 062084, KIZ 062088) were collected on 16-17 August 2018 by Ade P. Agung, Ada Chornelia, Jian-Mei Lu, L. Lee Grismer, Jesse L. Grismer, Evan S.H. Quah, Brian Folt, and Myin Kyaw Thura from forested karst in Ning’er County, Pu’er City, Yunnan Province, China (N23.083°,E101.019°; 1 436 m a.s.l.).

    Diagnosis:Hemiphyllodactylus simaoensissp. nov.can be distinguished from all congeners by a unique combination of the following characters: maximum SVL 40.87 mm; chin scales 7-10; enlarged postmentals; circumnasal scales 5-6;internasal scales 1-4; supralabial scales 8-12; infralabial scales 8-11; subdigital lamellae on fingers II-V (3 or 4)-(3-5)-(3-5)-(3 or 4); subdigital lamellae on toes II-V (3 or 4)-(3-5)-(3-5)-(3 or 4); dorsal scales 11-15; ventral scales 5-7; palegray base color on body, two lines of dark blotches running from neck to sacrum on dorsal side; dark postorbital stripe extending at least to base of neck; dorsolateral stripe indistinct or completely absent; ventrolateral stripe on trunk absent; dark postsacral markings bearing anteriorly projecting arms.

    Table 7 Summary statistics for analysis of each trait of Yanshan population and selected congeners

    Figure 4 Differences in adjusted mensural (left) and meristic (right) traits between specimens from Yanshan population and closely related species (H. huishuiensis and H. nahangensis)

    Table 8 PERMANOVA pairwise comparisons of multi-trait space of compared species within clade 6

    Figure 5 PCA of Yanshan population, H. huishuiensis, and H.nahangensis showing morphospatial relationships along first two components

    Description of holotype: Adult male, SVL 34.06 mm; head triangular in dorsal profile, depressed, distinct from neck (HL 6.03 mm; HW 6.77 mm); lores flat; snout long (SnEye 3.48 mm; SnEye/HL 58%) and narrow (SnW 1.26 mm; SnW/HW 19%); eyes large (ED 1.96 mm; ED/HL 33%); rostral scale wider than high, bordered posteriorly by two large supranasals and three internasals (IS); nares bordered anteriorly by rostral scale, ventrally by first supralabial scale, dorsally by supranasal scale, posteriorly by three postnasals; supralabials square, 10/9 (left/right), tapering from rostral to point in line with posterior margin of orbit (SL); infralabials square, 10/9(left/right), tapering from mental to point in line with posterior margin of orbit (IL); scales on head small, rounded, largest on rostrum; mental triangular, bordered by first infralabials and posteriorly by two enlarged postmentals; each postmental bordered anterolaterally by first infralabial; eight chin scales touching internal edges of infralabials from juncture of 2ndand 3rdinfralabial scales on left of mental scale to same juncture on right (Chin); scales in gular region rounded, nonoverlapping, becoming larger and more ovoid on venter. Body type robust and small, (TrunkL/SVL 49%), dorsoventrally compressed; dorsal body scales small, granular, 12 dorsal scales at midbody contained within one eye diameter; ventral body scales smooth and flat, much larger than dorsal scales,subimbricate, seven ventral scales at midbody contained within one eye diameter; forelimbs relatively short, covered dorsally with granular, subimbricate scales, smaller smooth scales ventrally; palmar scales flat, unevenly shaped, nonoverlapping; finger I vestigial, clawless, fingers II-V well developed; proximal subdigital lamellae undivided,rectangular; distal subdigital lamellae divided and undivided,angular, U-shaped, except terminal lamellae rounded,undivided; lamellar formula on fingers II-V (3)-(4)-(4)-(3) on both hands; claws on fingers II-V well developed, unsheathed,strongly curved; hind limbs short, covered dorsally with granular, subimbricate scales, smaller smooth scales ventrally; plantar scales flat, unevenly shaped, nonoverlapping; toe I vestigial, clawless, toes II-V well developed;proximal subdigital lamellae undivided, rectangular; distal subdigital lamellae divided and undivided, angular, U-shaped,except terminal lamellae rounded, undivided; lamellar formula on toes II-V (3)-(4)-(4)-(4) on both feet; claws on toes II-V well developed, unsheathed, strongly curved; one cloacal spur(CloacS) on each side; tail long, original (TL 32.32 mm;TL/SVL 94%), round in cross-section, dorsal scales on taillarger than on body and head, smaller than subcaudals, no plate-like subcaudal scales.

    Table 9 Test results for group separation of compared species within clade 6 along first two PCA axes (PC1 and PC2)

    Figure 6 Hemiphyllodactylus simaoensis sp. nov.

    Coloration in life(Figure 6): AllHemiphyllodactylusspecies can change the intensity and boldness of their coloration and patterns. The description below was taken when the holotype was photographed in the morning after capture. Base color of dorsal side of head, body, and limbs pale-gray, with two lines of dark blotches running from neck to sacrum on back; no dark marking on top of head; thin and indistinct dark pre- and postorbital stripes extending from external nares, through eye to just anterior of forelimb insertion on body; dark postsacral marking bearing anteriorly projecting arms; limbs with irregularly shaped, diffuse, and indistinct dark markings; dorsal side of tail with brown to black banded pattern; abdomen unicolor gray.

    Variation: Coloration of this species varies considerably(Figure 6). Variations in mensural and meristic data are presented in Supplementary Table S7.

    Distribution:Hemiphyllodactylus simaoensissp. nov.is currently only known from Simao District and Ning’er County,Pu’er City, Yunnan Province, China.

    Natural history:The holotype and seven paratypes (KIZ 062063, KIZ 062065-KIZ 062070) were collected on the evening of 6 July 2019 from a wall of a vacant building rarely used by humans in a forested karst area. The left and right sides of the building were densely covered by shrubs and trees. The back of the building was bordered by forested karst hills. The front of the building faced an old basketball court with several artificial lights, with the main road located further in front (Figure 7). The other 14 paratypes from Ning’er County (KIZ 062072-KIZ 062077, KIZ 062079-KIZ 062084,KIZ 062087-KIZ 062088) were collected one year earlier, on the evenings of 16 and 17 August 2018, from a cement wall in a forested karst area.

    Figure 7 Habitat of Hemiphyllodactylus simaoensis sp. nov.

    Etymology:The specific epithetsimaoensisrefers to the name of the district where the holotype originates, Simao District, Pu’er City, Yunnan Province, China.

    Morphological comparisons:A full list of trait comparisons is provided in Supplementary Table S3. Here, we describe morphological variations inHemiphyllodactylus simaoensis

    sp. nov.and differences from its congeners for traits that differed between species. However, in terms of the range of trait values,Hemiphyllodactylus simaoensissp. nov.was indistinguishable from theH. jinpingensisspecimen examined but distinct fromH. jinpingensisreported in Guo et al. (2015)(Supplementary Table S3). We consider this to be an artifact of researcher bias in the way data were taken. In terms of body ratios,Hemiphyllodactylus simaoensissp. nov.differs fromH. jinpingensis(values obtained from Guo et al. (2015)),H.montawaensis,H. tonywhitteni, andH. ngwelwiniby shorter head (HL/SVL), wider head (HW/HL), greater SnEye distance(SnEye/HL), greater NarEye distance (NarEye/HL), larger eyes (ED/HL), and wider snout (SnW/HL) (Table 10).

    In terms of scalation, the new species differs fromH.montawaensis,H. tonywhitteni, andH. ngwelwiniby fewer VS(5-7 vs. 7 or 8, 7-9, 7 or 8, respectively), and differs fromH.tonywhitteniby more CN (5 or 6 vs. 3-5). For body coloration and pattern, the new species differs fromH. montawaensisby presence of dark transverse blotches on dorsum (vs. reticulate pattern). The new species differs fromH. tonywhitteniby presence of dark transverse blotches on dorsum (vs. absent)and absence of light-colored dorsolateral spots on trunk (vs.present). The new species differs fromH. ngwelwiniby presence of dark transverse blotches on dorsum (vs. absent)(see Supplementary Table S3 for all comparative values).

    Table 10 Key trait differences between Hemiphyllodactylus simaoensis sp. nov. and its congeners

    Figure 8 Variations in color and patterns of Hemiphyllodactylus yanshanensis sp. nov.

    Hemiphyllodactylus yanshanensis sp. nov. (Figure 8)

    Suggested English name: Yanshan slender gecko

    Suggested Chinese name: 硯山半葉趾虎

    Holotype: Adult male (KIZ 062090) collected on 1 July 2019 by Ade P. Agung, Jian-Mei Lu, and Zong-Bao Yang from forested karst hills in Yanshan County, Yunnan Province,China (N23.61680°, E104.41669°; 1 536 m a.s.l.).

    Paratypes: Nine adult females (KIZ 062091-KIZ 062096, KIZ 062100-KIZ 062102) and three adult males (KIZ 062097-KIZ 062099), same data as holotype.

    Diagnosis:Hemiphyllodactylus yanshanensissp. nov. can be distinguished from all congeners by a unique combination of the following characters: maximum SVL 46.28 mm; chin scales 8-11; enlarged postmentals; circumnasal scales 5-6;internasal scales 1-3; supralabial scales 9-12; infralabial scales 9-12; ventral scales 7-13; dorsal scales 11-22;subdigital lamellae on fingers II-V (4 or 5)-(5-7)-(5-7)-(4 or 5);subdigital lamellae on toes II-V (4 or 5)-(5 or 6)-(5-7)-(5 or 6);subdigital lamellae wider than long on first finger (4 or 5);subdigital lamellae wider than long on first toe (5 or 6); pale brown color base on body with various transverse blotched patterns on dorsum; dark postorbital stripe extending at least to base of neck; dorsolateral stripe on trunk present;ventrolateral stripe on trunk absent; postsacral marking bearing anteriorly projecting arms.

    Description of holotype: Adult male, SVL 40.03 mm; head triangular in dorsal profile, depressed, distinct from neck (HL 7.27 mm; HW 7.68 mm); lores flat; snout long (SnEye 4.07 mm; SnEye/HL 56%) and narrow (SnW 1.52 mm; SnW/HW 20%); eyes large (ED 2.27 mm; ED/HL 31%); rostral scale wider than high, bordered posteriorly by two large supranasals and two internasals (IS); nares bordered anteriorly by rostral scale, ventrally by first supralabial scale, dorsally by supranasal scale, posteriorly by three postnasals; supralabials square, 9/10 (left/right), tapering from rostral to point in line with posterior margin of orbit (SL); infralabials square, 10/10(left/right), tapering from mental to point in line with posterior margin of orbit (IL); scales on head small, rounded, largest on rostrum; mental triangular, bordered by first infralabials and posteriorly by two enlarged postmentals; each postmental bordered anterolaterally by first infralabial; nine chin scales touching internal edges of infralabials from juncture of 2ndand 3rdinfralabial scales on left of mental scale to same juncture on right (Chin); scales in gular region rounded, nonoverlapping, becoming larger and more ovoid on venter. Body type robust and small (TrunkL/SVL 47%), dorsoventrally compressed; dorsal body scales small, granular, 14 dorsal scales at midbody contained within one eye diameter; ventral body scales smooth and flat, much larger than dorsal scales,subimbricate, eight ventral scales at midbody contained within one eye diameter; forelimbs relatively short, covered dorsally with granular, subimbricate scales, smaller smooth scales ventrally; palmar scales flat, unevenly shaped, nonoverlapping; finger I vestigial, clawless, fingers II-V well developed; proximal subdigital lamellae undivided,rectangular; distal subdigital lamellae divided and undivided,angular, U-shaped, except terminal lamellae rounded,undivided; lamellar formula on fingers II-V (4)-(5)-(5)-(4) on both hands; claws on fingers II-V well developed, unsheathed,strongly curved; hind limbs short, covered dorsally with granular, subimbricate scales, smaller smooth scales ventrally; plantar scales flat, unevenly shaped, nonoverlapping; toe I vestigial, clawless, toes II-V well developed;proximal subdigital lamellae undivided, rectangular; distal subdigital lamellae divided and undivided, angular, U-shaped,except terminal lamellae rounded, undivided; lamellar formula on toes II-V (4)-(5)-(5)-(5) on left foot and (4)-(6)-(6)-(5) on right foot; claws on toes II-V well developed, unsheathed,strongly curved; one cloacal spur (CloacS) on each side; tail long, original (TL 34.01 mm; TL/SVL 85%), round in crosssection, dorsal scales on tail larger than on body and head,smaller than subcaudals, no plate-like subcaudal scales.

    Coloration in life(Figure 8): AllHemiphyllodactylusspecies can change the intensity and boldness of their coloration and patterns. The description below was taken when the holotype was photographed in the morning after capture. Base color of dorsal side of head, body, and limbs pale-brown and densely mottled with darker markings; top of head overlain with indistinct blotches; indistinct pre- and postorbital stripes extending from external nares, through eye to just anterior of forelimb insertion on body; postsacral marking bearing anteriorly projecting arms; limbs with irregularly shaped,diffuse, dark markings; tail pale-gray, with several irregularly shaped dark markings on dorsal side, diffused on lateral sides;abdomen unicolor gray.

    Variation: Coloration of this species varies considerably(Figure 8). Variations in mensural and meristic data are presented in Supplementary Table S8.

    Distribution:Hemiphyllodactylus yanshanensissp. nov.is only known from the type locality in Yanshan County, Yunnan Province, China (Figure 9).

    Natural history:All specimens were collected during the evening of 1 July 2019 on the surface of a forested karst hill in Yanshan County, Yunnan Province, China (Figure 10). There were several gravid females among the specimens collected,each containing two eggs. We also found two eggs placed together loosely in the crevices of the same hill, which we assumed were laid byHemiphyllodactylus yanshanensissp.nov.based on their size (Figure 10), indicating that its reproductive season extends into July (see reports on otherHemiphyllodactylusspecies in Cobos et al. (2016) and Grismer et al. (2015)). The karst hills were surrounded by paddy fields and several huts and houses were nearby.

    Etymology:The specific epithetyanshanensisrefers to the name of Yanshan County where the specimens were found.

    Morphological comparisons:A full list of trait comparisons is provided in Supplementary Table S5. Here, we describe morphological variations inHemiphyllodactylus yanshanensis

    Figure 9 Map of Yunnan with type localities of Hemiphyllodactylus simaoensis sp. nov. (white circles) and Hemiphyllodactylus yanshanensis sp. nov. (black star)

    sp. nov.and differences fromH. huishuiensisandH.nahangensisfor traits that differed between species. In terms of body ratios,Hemiphyllodactylus yanshanensissp. nov.differs fromH.huishuiensisandH. nahangensisby shorter head (HL/SVL), wider head (HW/HL), greater SnEye distance(SnEye/HL), greater NarEye distance (NarEye/HL), larger eyes (ED/HL), and wider snout (SnW/HL) (Table 11).

    In terms of scalation, the new species differs fromH.huishuiensisandH. nahangensisby more circumnasal scales(5 or 6 vs. 3 and 2 or 3, respectively) and more subdigital lamellae wider than long on first finger and toe (SL1F=4 or 5 vs. 3; SL1T=5 or 6 vs. 3). Furthermore, the new species also differs fromH. nahangensisby fewer internasal scales(IS=1-3 vs. 4-6) and fewer femoroprecloacal pores or pitted scales in females (10-16 vs. 22-24). In coloration, the new species differs fromH. huishuiensisby presence of lightcolored dorsolateral spots on trunk (vs. absent). It also differs fromH. nahangensisby presence of dark postorbital stripe and dark transverse dorsal blotches (vs. absent) (see Supplementary Table S5 for all comparative values).

    DlSCUSSlON

    Diversity of Hemiphyllodactylus in southern Chinese karsts

    Most Chinese specimens ofHemiphyllodactyluswere previously identified asH. yunnanensisBoulenger 1903. In 2013, however, Grismer et al. (2013) conducted an integrative taxonomic study and revealed that this taxon is much more diverse than previously reported. A number of species have since been described in China. Most recently,H.zhutangxiangensiswas described from Lancang Lahu Autonomous County (Agung et al., 2021) andH.zayuensiswas described from Zayu County, Tibet (Che et al., 2020),with both species considered likely members of clade 3, given that all members are distributed in western Yunnan and northern Myanmar. The population in Dao County, Hunan Province, was also recently identified asH. dupanglingensis(Zhang et al., 2020) and as sister toH. zugi, a species that occurs on the border of Vietnam and Guangxi Province, China(Nguyen et al., 2013). Both are members of clade 6, together withH. hongkongensis(Sung et al., 2018),H. dushanensis(Zhou et al., 1981), andH. huishuiensis(Yan et al., 2016),which occur in southeastern China. Thus, together withH.yunnanensis,H. longlingensis,H. jinpingensis,H.changningensis(Guo et al., 2015), andH. typusBleeker (Uetz et al., 2021), a total of 12 species are currently known in China. However, given the extensive and unsurveyed karst habitat across southern China, many more species are likely to be discovered.

    Table 11 Key trait differences between Hemiphyllodactylus yanshanensis sp. nov. and its congeners

    MitochondrialND2analysis indicated thatHemiphyllodactylussimaoensissp.nov.andHemiphyllodactylus yanshanensissp. nov.are members of clades 4 and 6, respectively, and are embedded within thetypusgroup.Hemiphyllodactylus simaoensissp. nov.appears to be most closely related toH. jinpingensis(6.3%; Table 2) in clade 4, with most members of this clade distributed on Shan Plateau, eastern Myanmar.Hemiphyllodactylus yanshanensis

    sp. nov.appears to be most closely related toH. huishuiensisfrom southeastern China (4.3%; Table 3) in clade 6. These distributions and relationships highlight the complex biogeography of Yunnan. The discovery of these two new species increases the total number ofHemiphyllodactylusspecies in China to 14, four of which have been described in the last five years. Thus, further research is required as many species likely remain to be discovered in China and other parts of Asia.

    Our Yunnan samples were nested in four different clades(clades 3, 4, 6 and 7; Figure 1), suggesting that Yunnan may have been colonized on at least four separate occasions, as the sister species to clades 3 and 4 (H. khlonglanensis) and clades 6 and 7 (H. pardalis) are only known from Thailand.Nevertheless, more population data are needed to determine if the splits between clades 3 and 4 and clades 6 and 7 occurred before or after colonization of Yunnan.

    The description ofHemiphyllodactylus yanshanensissp.nov.in this paper is the first evidence of colonization in Yunnan by a member of clade 6, with most clade members occurring in eastern Indochina (Vietnam and eastern Laos)and southeastern China (Guizhou, Hunan, and Hongkong SAR). These species discoveries emphasize thatHemiphyllodactylusdiversity within China and neighboring regions is underestimated.

    Trait reliability

    The newly described species (Hemiphyllodactylus simaoensis

    sp. nov.andHemiphyllodactylus yanshanensissp. nov.)differed significantly in several morphological traits from their sister species (see Results). Nevertheless, these morphological traits are challenging to use in the field without specimen comparisons, and thus, obtaining tissue samples for genetic analysis is crucial to corroborate what is noted in the field. The use of integrated taxonomic approaches with genetic and morphological data, coupled with robust statistical analyses, will increase the accuracy of species delimitation.Therefore, given that many species inHemiphyllodactylusare morphologically convergent and cryptic (Grismer et al., 2013),single morphological evidence may not be sufficient for species identification inHemiphyllodactylus.

    Conservation implications

    Our study suggests that karst regions in Yunnan, which cover nearly 44% of the province (Nester, 2021), are likely to harbor additional undescribed species ofHemiphyllodactylus.However, despite the recent delineation of geoparks (e.g.,Shilin, Dali-Cangshan), these hotspots of endemism are disproportionately under-protected. Furthermore,H.yunnanensismay be a species complex (clade 7) and require further taxonomic revision. Based on this preliminary analysis,we hypothesize thatHemiphyllodactylusin the karst regions of Yunnan has had a complex colonization history, given that we potentially found multiple species in one site as well as one species distributed in multiple karst hills (Agung et al.,unpublished data).

    The continuation of species discoveries with increasing field surveys may clarify the status of isolated populations ofHemiphyllodactylusspecies groups. Karst outcrops provide biodiversity arks and focal points for speciation (Clements et al., 2006; Grismer et al., 2021), with a high degree of isolation forHemiphyllodactyluspopulations. Due to their poor dispersal abilities, different populations in different karst areas may have evolved into separate lineages, and thus could be recognized as distinct species (Grismer, 1999). Considering the allopatric populations ofHemiphyllodactylus, with presumably no or low rates of gene flow among them, they may have become phenetically and genetically distinguishable over time (as reviewed in de Queiroz, 2007). In the last 10 years, approximately 31 new species have been identified throughout Asia, with each species usually described from a single karst area or complex (i.e., Grismer et al., 2020a),highlighting the high rate of localHemiphyllodactylusendemism in each karst region.Hemiphyllodactylus simaoensissp. nov.andHemiphyllodactylus yanshanensis

    sp. nov.exhibited clear genetic and phenotypic differences from their sister lineages,H. jinpingensis(located more than 200 km away) andH. huishuiensis(located more than 350 km away), respectively, thus supporting the assumption of high local endemism for each karst region, even though some areas may have multiple sympatric species. High endemism in karst geckos is not only reported forHemiphyllodactylus, but also forCyrtodactylus(Davis et al., 2019; Grismer et al., 2021;Luu et al., 2016; Nazarov et al., 2018; Nguyen et al., 2017;Pauwels et al., 2016) andCnemaspis(Grismer et al., 2014;Wood et al., 2017). In addition, high numbers of endemic flora and invertebrates have also been reported from limestone forest habitats (Clements et al., 2006; Marzuki et al., 2021;Nguyen et al., 2021). However, the genusHemiphyllodactylusremains understudied and increased survey efforts with broader geographic coverage will likely increase the number of new species discovered. Notably, many karsts still lack inventories and given the high rate of endemism, much work is needed to map taxa across the region and develop appropriate plans for protection and representative coverage of species across this naturally fragmented system.

    Neglected habitats, such as karst ecosystems, require indepth investigation and survey. Quantifying the biodiversity of a region also requires accurate taxonomy, as species are the fundamental units of conservation planning (Mace, 2004).Conserving karst landscapes and the communities that depend upon them cannot be achieved without adequate knowledge of existing species. Major threats to gecko species,including the two described here, are habitat destruction from quarrying and deforestation. Furthermore, with the rapid rate of karst loss across much of Southeast Asia, failure to identify and map the ranges of species may ultimately mean species become extinct without formal description. Therefore,appropriate conservation management plans are needed to reduce the loss of key habitats and risk of potential invasion by non-native species. Without identifying species,understanding their ranges, and increasing systematic surveys of these diverse and complex regions, the extinction of species may occur before they are described, or we may only observe specimens whose habitats have already been destroyed. Mapping the distributions of range-limited herptiles and understanding their vulnerabilities are crucial factors for targeted management and conservation. For example,research on the Chinese salamander (Andrias jiangxiensis,Chai et al., 2022) has highlighted how few sites retain the species, underscoring the need for the identification and protection of these sites. Our study and that of Chai et al.(2022) demonstrate the need for an integrated approach,including genetics and morphology, to identify species and form the basis for management and conservation. Therefore,further work is urgently needed to understand cryptic distributions and enable targeted management.

    Synthesis and future directions

    Prior to 2021, only 12 species ofHemiphyllodactyluswere known in southern China. Recent studies have transformed our knowledge regarding the group and added new species based on a variety of different traits and molecular data.However, these recent studies also highlight the need to better understand regional diversity patterns and biogeography and provide a basis for future management and conservation given the high rates of habitat loss across the range. The description of multiple species from Yunnan in the last 10 years also emphasizes the need to redouble our efforts to ensure sufficient representation in sampling. Furthermore, given the lack of any single trait to differentiate between species (and the need to analyze multiple traits simultaneously), integration of morphological and genetic analyses is required to better identify species, although field identification may continue to be a challenge. Key regions where we expect high complexity in species populations (due to high environmental heterogeneity forming barriers to gene flow) include the Yuanjiang valley and Ailao mountain range in Yunnan. Studies in countries such as Malaysia and Myanmar have demonstrated high levels of regional endemism within the group, indicating that more work is needed to sample and describeHemiphyllodactylusin Yunnan. Without knowledge of where species are found, developing effective and efficient management will remain difficult. Given the rate of species description in this neglected taxon, as well as the high rates of habitat loss, further research is needed to describe species and map their ranges before they can be protected, especially for genera with high levels of site-specific endemism in inventory-poor areas.

    NOMENCLATURAL ACTS REGlSTRATlON

    The electronic version of this article in portable document format will represent a published work according to the International Commission on Zoological Nomenclature (ICZN),and hence the new names contained in the electronic version are effectively published under that Code from the electronic edition alone (see Articles 8.5-8.6 of the Code). This published work and the nomenclatural acts it contains have been registered in ZooBank, the online registration system for the ICZN. The ZooBank LSIDs (Life Science Identifiers) can be resolved and the associated information can be viewed through any standard web browser by appending the LSID to the prefix http://zoobank.org/.

    Publication LSID: urn:lsid:zoobank.org:pub:763DE16C-BA54-41ED-910A-F0423AA62A98.

    Nomenclatural act LSID forHemiphyllodactylus simaoensissp. nov.: urn:lsid:zoobank.org:act:A863C059-F959-4450-BC D9-6DF3B92D5E68.Nomenclatural act LSID forHemiphyllodactylus yanshanensissp. nov.: urn:lsid:zoobank.org:act:02D7C54D-EFEC-4225-86 49-BC1D7FED3617.

    SClENTlFlC FlELD SURVEY PERMlSSlON lNFORMATlON

    The Ethics Committee of Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, approved the study and provided ethics permission. All local regulations were followed in the collection of samples and specimens for the study.

    SUPPLEMENTARY DATA

    Supplementary data to this article can be found online.

    COMPETlNG lNTERESTS

    The authors declare that they have no competing interests

    AUTHORS’ CONTRlBUTlONS

    A.P.A., L.L.G., J.L.G., A.C.H. designed the study. A.P.A., A.C.,L.L.G., J.L.G., E.S.H.Q., and J.M.L. collected data. A.P.A.,A.C., and L.L.G. conducted phylogenetic analyses. A.P.A.examined morphology. A.P.A and K.W.T conducted statistical analyses. A.P.A. and E.S.H.Q. wrote the manuscript. A.P.A.,L.L.G., and A.C.H. revised the manuscript. All authors read and approved the final version of the manuscript.

    ACKNOWLEDGEMENTS

    We thank Brian Folt, Myint Kyaw Thura, Zhong-Bao Yang, and Shi-Fu Pu for help during fieldwork, Tuanjit Sritongchuay for loaning equipment for specimen measurements, and Yan-Hua Chen and Zi-Nan Ding for help with administrative work.

    猜你喜歡
    思茅
    思茅關(guān)貿(mào)易研究
    2016—2022年云南思茅酸雨變化特征分析
    思茅區(qū)總工會(huì):舉辦游園活動(dòng)
    空港
    青年生活(2019年35期)2019-09-10 00:25:43
    思茅山橙根中生物堿類成分及其抗腫瘤活性研究
    不同波長(zhǎng)誘蟲(chóng)燈對(duì)思茅松毛蟲(chóng)監(jiān)測(cè)研究
    區(qū)長(zhǎng)陳奇調(diào)研思茅區(qū)綜合檔案館建設(shè)
    云南檔案(2018年3期)2018-01-24 15:43:18
    思茅趕街天 端午
    普洱(2016年12期)2017-05-26 16:16:53
    云南櫻桃樹(shù)羅塘銅礦床地質(zhì)特征及找礦標(biāo)志
    西部資源(2016年1期)2016-07-30 07:00:30
    思茅區(qū)生豬產(chǎn)業(yè)發(fā)展的現(xiàn)狀分析
    国产黄频视频在线观看| 中文字幕av电影在线播放| 免费黄色在线免费观看| 亚洲国产看品久久| 天堂中文最新版在线下载| 婷婷色综合大香蕉| 午夜免费鲁丝| 日韩中文字幕欧美一区二区 | 三级国产精品片| 男女午夜视频在线观看| 国产熟女欧美一区二区| 日韩伦理黄色片| 丰满迷人的少妇在线观看| 人妻人人澡人人爽人人| 一个人免费看片子| 成年美女黄网站色视频大全免费| 中文精品一卡2卡3卡4更新| 欧美日韩国产mv在线观看视频| 在线亚洲精品国产二区图片欧美| 成人二区视频| 免费看不卡的av| 久久99一区二区三区| 亚洲国产av新网站| 午夜老司机福利剧场| 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人爽人人夜夜| 这个男人来自地球电影免费观看 | 国产精品熟女久久久久浪| 久久精品国产a三级三级三级| 久久鲁丝午夜福利片| 亚洲精品,欧美精品| 免费av中文字幕在线| 亚洲欧美中文字幕日韩二区| 国产色婷婷99| 国产成人精品久久久久久| av片东京热男人的天堂| 色婷婷久久久亚洲欧美| h视频一区二区三区| 午夜福利在线观看免费完整高清在| 少妇人妻久久综合中文| 国产精品久久久久成人av| 亚洲精品一区蜜桃| 狠狠婷婷综合久久久久久88av| 26uuu在线亚洲综合色| xxxhd国产人妻xxx| 伊人久久大香线蕉亚洲五| 日韩,欧美,国产一区二区三区| 国产精品成人在线| av在线播放精品| 国产一级毛片在线| 黄色一级大片看看| 亚洲美女视频黄频| www.熟女人妻精品国产| 亚洲欧美日韩另类电影网站| 伦精品一区二区三区| 一级,二级,三级黄色视频| 国产亚洲午夜精品一区二区久久| 日韩一区二区三区影片| 国产伦理片在线播放av一区| 纵有疾风起免费观看全集完整版| 久久精品国产自在天天线| 午夜久久久在线观看| 亚洲国产av影院在线观看| 18禁国产床啪视频网站| 人妻一区二区av| 99精国产麻豆久久婷婷| 国产成人精品久久二区二区91 | 青春草国产在线视频| 曰老女人黄片| 午夜激情av网站| 婷婷成人精品国产| 亚洲少妇的诱惑av| 亚洲人成网站在线观看播放| videosex国产| 日韩av不卡免费在线播放| 亚洲av成人精品一二三区| 国产激情久久老熟女| 毛片一级片免费看久久久久| 97精品久久久久久久久久精品| av线在线观看网站| 亚洲视频免费观看视频| 老司机影院毛片| 午夜日本视频在线| 久久这里只有精品19| 精品一区在线观看国产| 精品卡一卡二卡四卡免费| 亚洲av综合色区一区| 亚洲人成网站在线观看播放| 国产日韩一区二区三区精品不卡| 9热在线视频观看99| 80岁老熟妇乱子伦牲交| 五月伊人婷婷丁香| 午夜福利影视在线免费观看| 久久久精品免费免费高清| 久久久亚洲精品成人影院| 激情五月婷婷亚洲| 亚洲精品国产av蜜桃| av有码第一页| 成年女人毛片免费观看观看9 | √禁漫天堂资源中文www| 水蜜桃什么品种好| 十八禁网站网址无遮挡| 最新中文字幕久久久久| 欧美国产精品va在线观看不卡| 男人添女人高潮全过程视频| 如日韩欧美国产精品一区二区三区| 看免费av毛片| 国产又爽黄色视频| 国产 一区精品| 中文字幕av电影在线播放| 日韩制服丝袜自拍偷拍| 男女下面插进去视频免费观看| www.熟女人妻精品国产| 老汉色∧v一级毛片| 国产在线免费精品| 啦啦啦在线观看免费高清www| 亚洲欧美中文字幕日韩二区| 不卡视频在线观看欧美| 七月丁香在线播放| 亚洲国产色片| 99国产精品免费福利视频| 欧美bdsm另类| 日本欧美国产在线视频| 欧美+日韩+精品| 日韩伦理黄色片| 日韩免费高清中文字幕av| 久久婷婷青草| 中文字幕人妻熟女乱码| 午夜91福利影院| 天天影视国产精品| 这个男人来自地球电影免费观看 | 三上悠亚av全集在线观看| 亚洲人成77777在线视频| 亚洲成国产人片在线观看| a 毛片基地| 国产精品不卡视频一区二区| 日韩制服骚丝袜av| 九色亚洲精品在线播放| 成人免费观看视频高清| 黑人猛操日本美女一级片| 热re99久久精品国产66热6| 女人久久www免费人成看片| 少妇的丰满在线观看| 久久久久精品性色| www.av在线官网国产| 亚洲综合色惰| 人妻少妇偷人精品九色| 成人免费观看视频高清| 国产精品 欧美亚洲| 在线观看www视频免费| av有码第一页| av一本久久久久| 婷婷色av中文字幕| 寂寞人妻少妇视频99o| 欧美精品av麻豆av| 极品少妇高潮喷水抽搐| 最新的欧美精品一区二区| 亚洲 欧美一区二区三区| 男女午夜视频在线观看| 青草久久国产| 国产精品久久久久成人av| 国产成人午夜福利电影在线观看| 亚洲综合色惰| 超碰97精品在线观看| 老司机影院毛片| 亚洲精品国产一区二区精华液| 国产日韩欧美视频二区| 丝袜美腿诱惑在线| 久久综合国产亚洲精品| 精品一区二区三卡| 久久久久精品久久久久真实原创| 日韩 亚洲 欧美在线| 免费大片黄手机在线观看| av免费在线看不卡| 侵犯人妻中文字幕一二三四区| 精品人妻熟女毛片av久久网站| 性少妇av在线| 建设人人有责人人尽责人人享有的| 亚洲三级黄色毛片| 亚洲av中文av极速乱| 制服人妻中文乱码| 热re99久久精品国产66热6| www.av在线官网国产| 亚洲综合色网址| 一区二区三区激情视频| 国产成人精品久久久久久| 黑人欧美特级aaaaaa片| 国产又爽黄色视频| 日韩不卡一区二区三区视频在线| 国产黄色视频一区二区在线观看| 久久国产精品大桥未久av| 熟女少妇亚洲综合色aaa.| 免费看av在线观看网站| 男女下面插进去视频免费观看| 日韩成人av中文字幕在线观看| 国产乱人偷精品视频| 亚洲精品视频女| 国产1区2区3区精品| 国产成人精品一,二区| 成人手机av| 日韩欧美精品免费久久| 十八禁网站网址无遮挡| 激情五月婷婷亚洲| 天天影视国产精品| 亚洲国产av影院在线观看| 少妇 在线观看| 午夜福利在线免费观看网站| 人人澡人人妻人| 极品少妇高潮喷水抽搐| 最黄视频免费看| 国产精品久久久久久精品电影小说| 99re6热这里在线精品视频| 久久久久视频综合| 久久久久久人妻| 边亲边吃奶的免费视频| 国精品久久久久久国模美| 日韩精品免费视频一区二区三区| 日日撸夜夜添| 女人被躁到高潮嗷嗷叫费观| 日韩人妻精品一区2区三区| 中文字幕av电影在线播放| 午夜免费观看性视频| 亚洲美女视频黄频| 国产亚洲一区二区精品| 国产免费福利视频在线观看| 色播在线永久视频| 宅男免费午夜| av免费观看日本| 久久久久国产一级毛片高清牌| 久久久久久久久久久久大奶| 日产精品乱码卡一卡2卡三| 久久精品久久精品一区二区三区| 丝瓜视频免费看黄片| 久久人妻熟女aⅴ| 亚洲av国产av综合av卡| 日韩av免费高清视频| 久久人人爽人人片av| 18禁国产床啪视频网站| 精品福利永久在线观看| 久久久久国产一级毛片高清牌| 亚洲一区中文字幕在线| 人妻 亚洲 视频| 日韩av在线免费看完整版不卡| 91久久精品国产一区二区三区| 青春草国产在线视频| 国产探花极品一区二区| 欧美人与性动交α欧美精品济南到 | 日本wwww免费看| 一级毛片我不卡| 丝袜喷水一区| 亚洲国产精品一区二区三区在线| 在线观看三级黄色| 国产老妇伦熟女老妇高清| 日本欧美视频一区| 国产高清不卡午夜福利| 美女国产视频在线观看| av电影中文网址| 亚洲三级黄色毛片| 亚洲人成电影观看| 一二三四中文在线观看免费高清| 成年av动漫网址| 最近中文字幕2019免费版| 一本—道久久a久久精品蜜桃钙片| www日本在线高清视频| 考比视频在线观看| 街头女战士在线观看网站| 看免费av毛片| 免费黄网站久久成人精品| 你懂的网址亚洲精品在线观看| √禁漫天堂资源中文www| 啦啦啦啦在线视频资源| 女性生殖器流出的白浆| 久久99热这里只频精品6学生| h视频一区二区三区| 日韩不卡一区二区三区视频在线| 国产无遮挡羞羞视频在线观看| av网站免费在线观看视频| 午夜老司机福利剧场| 国产欧美亚洲国产| 色视频在线一区二区三区| 两个人免费观看高清视频| 一本色道久久久久久精品综合| 女的被弄到高潮叫床怎么办| 免费观看a级毛片全部| 久久久久国产精品人妻一区二区| 在线看a的网站| 好男人视频免费观看在线| 在线精品无人区一区二区三| 久久影院123| 一区福利在线观看| 最近最新中文字幕免费大全7| 午夜精品国产一区二区电影| 哪个播放器可以免费观看大片| 97精品久久久久久久久久精品| 国产免费福利视频在线观看| 不卡av一区二区三区| 伦理电影大哥的女人| 青青草视频在线视频观看| h视频一区二区三区| 18禁裸乳无遮挡动漫免费视频| 伊人久久大香线蕉亚洲五| 国产精品无大码| 亚洲五月色婷婷综合| 国产成人精品久久久久久| 丰满少妇做爰视频| 老司机影院毛片| 日日摸夜夜添夜夜爱| 宅男免费午夜| 精品一区在线观看国产| 男人添女人高潮全过程视频| 国产亚洲欧美精品永久| 日韩电影二区| 久久久国产一区二区| 亚洲av在线观看美女高潮| 欧美人与性动交α欧美精品济南到 | 欧美成人精品欧美一级黄| 免费人妻精品一区二区三区视频| 精品少妇一区二区三区视频日本电影 | 亚洲在久久综合| www.熟女人妻精品国产| 久久毛片免费看一区二区三区| 2018国产大陆天天弄谢| 亚洲精品久久成人aⅴ小说| 一本—道久久a久久精品蜜桃钙片| 亚洲伊人色综图| 久久精品久久久久久噜噜老黄| 另类精品久久| 亚洲色图 男人天堂 中文字幕| 在线免费观看不下载黄p国产| 高清黄色对白视频在线免费看| 一级毛片我不卡| 成人亚洲精品一区在线观看| 王馨瑶露胸无遮挡在线观看| 精品国产一区二区三区久久久樱花| 91精品伊人久久大香线蕉| 精品少妇黑人巨大在线播放| 亚洲色图综合在线观看| 五月开心婷婷网| 校园人妻丝袜中文字幕| 国产欧美亚洲国产| 91精品伊人久久大香线蕉| 看免费成人av毛片| 日日爽夜夜爽网站| 久久久久久久久久久免费av| 毛片一级片免费看久久久久| 国产av精品麻豆| 在线观看免费高清a一片| 久久久久精品人妻al黑| 国产成人免费无遮挡视频| 精品少妇黑人巨大在线播放| 啦啦啦在线免费观看视频4| 亚洲,欧美,日韩| 黄片播放在线免费| 91精品三级在线观看| 国产一区二区 视频在线| 国产精品不卡视频一区二区| 我的亚洲天堂| 免费大片黄手机在线观看| 亚洲 欧美一区二区三区| 日韩,欧美,国产一区二区三区| 天堂8中文在线网| 熟女av电影| 青春草亚洲视频在线观看| 亚洲精华国产精华液的使用体验| 国产片内射在线| 国产av一区二区精品久久| 亚洲内射少妇av| 精品一品国产午夜福利视频| 国产人伦9x9x在线观看 | 一区二区三区四区激情视频| 黄色一级大片看看| videossex国产| 国产精品一区二区在线观看99| 国产免费视频播放在线视频| 最近中文字幕2019免费版| 国语对白做爰xxxⅹ性视频网站| 99精国产麻豆久久婷婷| 一级毛片我不卡| 国产精品av久久久久免费| 亚洲综合精品二区| 久久99一区二区三区| 亚洲国产色片| 亚洲精品久久久久久婷婷小说| 国产乱来视频区| 街头女战士在线观看网站| 免费看av在线观看网站| 久久久久精品久久久久真实原创| 精品人妻一区二区三区麻豆| 中文字幕色久视频| 美女中出高潮动态图| 久久久久视频综合| 午夜免费观看性视频| 亚洲欧洲精品一区二区精品久久久 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲综合色网址| 秋霞伦理黄片| 久久综合国产亚洲精品| 国产精品.久久久| 丝袜美足系列| 一级毛片黄色毛片免费观看视频| 久久久久精品性色| 亚洲精品久久久久久婷婷小说| 在线观看www视频免费| 一区在线观看完整版| 狠狠精品人妻久久久久久综合| 免费人妻精品一区二区三区视频| 亚洲精品在线美女| 99热国产这里只有精品6| 99久久人妻综合| 桃花免费在线播放| 精品一品国产午夜福利视频| 精品第一国产精品| 女人久久www免费人成看片| av国产久精品久网站免费入址| 精品国产露脸久久av麻豆| 一本大道久久a久久精品| 成年美女黄网站色视频大全免费| 亚洲欧美精品综合一区二区三区 | 黄色视频在线播放观看不卡| 久久精品亚洲av国产电影网| 中文字幕制服av| 搡女人真爽免费视频火全软件| 丝瓜视频免费看黄片| 十分钟在线观看高清视频www| tube8黄色片| 日本猛色少妇xxxxx猛交久久| 热99久久久久精品小说推荐| 91在线精品国自产拍蜜月| 性高湖久久久久久久久免费观看| 性色avwww在线观看| 久久精品国产亚洲av高清一级| 亚洲欧美一区二区三区国产| 亚洲天堂av无毛| 日韩制服丝袜自拍偷拍| xxx大片免费视频| 热99久久久久精品小说推荐| 国产精品人妻久久久影院| tube8黄色片| 国产男女内射视频| 亚洲第一青青草原| 卡戴珊不雅视频在线播放| 夫妻性生交免费视频一级片| 免费日韩欧美在线观看| 亚洲精品视频女| 老汉色∧v一级毛片| 人人妻人人澡人人看| 男女午夜视频在线观看| 国产精品 欧美亚洲| 亚洲一区二区三区欧美精品| 99九九在线精品视频| 国产福利在线免费观看视频| 性色av一级| 国产国语露脸激情在线看| 波多野结衣av一区二区av| 高清av免费在线| 亚洲欧美成人精品一区二区| 久久精品久久精品一区二区三区| 精品国产一区二区三区久久久樱花| 亚洲精品日本国产第一区| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美日韩在线播放| av片东京热男人的天堂| 一区福利在线观看| 黄片播放在线免费| 90打野战视频偷拍视频| 国产精品熟女久久久久浪| freevideosex欧美| 中文字幕制服av| 免费不卡的大黄色大毛片视频在线观看| 精品人妻在线不人妻| 人人妻人人添人人爽欧美一区卜| 国产精品一区二区在线观看99| 老司机影院成人| 国产精品偷伦视频观看了| 香蕉精品网在线| 性少妇av在线| 亚洲人成网站在线观看播放| 国产av国产精品国产| 中文字幕人妻丝袜制服| 亚洲欧美日韩另类电影网站| 亚洲欧洲国产日韩| 午夜影院在线不卡| 亚洲精品国产色婷婷电影| 国产欧美亚洲国产| 欧美 亚洲 国产 日韩一| 婷婷色麻豆天堂久久| 久久久久久人妻| 国产成人精品婷婷| 欧美+日韩+精品| 咕卡用的链子| 婷婷色麻豆天堂久久| 熟妇人妻不卡中文字幕| 一级片'在线观看视频| 欧美人与善性xxx| 色播在线永久视频| 电影成人av| 在线亚洲精品国产二区图片欧美| 中文精品一卡2卡3卡4更新| 18禁国产床啪视频网站| 色播在线永久视频| 电影成人av| 国产欧美日韩一区二区三区在线| 可以免费在线观看a视频的电影网站 | 女人久久www免费人成看片| 老司机亚洲免费影院| 超色免费av| 十八禁高潮呻吟视频| 欧美另类一区| 国产成人精品在线电影| 欧美精品亚洲一区二区| 99热国产这里只有精品6| 青春草视频在线免费观看| 久久久久久久精品精品| av.在线天堂| 亚洲三区欧美一区| av片东京热男人的天堂| 国产又色又爽无遮挡免| 卡戴珊不雅视频在线播放| 人成视频在线观看免费观看| 纯流量卡能插随身wifi吗| 大码成人一级视频| 深夜精品福利| 亚洲四区av| 欧美日韩视频精品一区| 日日啪夜夜爽| 国产日韩一区二区三区精品不卡| 丝袜美腿诱惑在线| 一区二区三区激情视频| 成人国语在线视频| 亚洲男人天堂网一区| 亚洲精品国产av成人精品| 免费人妻精品一区二区三区视频| 国产亚洲av片在线观看秒播厂| 欧美精品av麻豆av| 最近中文字幕2019免费版| 三级国产精品片| 国产97色在线日韩免费| 人妻一区二区av| 国产精品久久久久成人av| 国产精品一国产av| 一级爰片在线观看| 最近中文字幕高清免费大全6| 高清黄色对白视频在线免费看| 香蕉精品网在线| 熟女少妇亚洲综合色aaa.| 女性生殖器流出的白浆| 久久精品夜色国产| 日韩不卡一区二区三区视频在线| av视频免费观看在线观看| 好男人视频免费观看在线| 最新的欧美精品一区二区| 涩涩av久久男人的天堂| 亚洲国产精品一区三区| 免费女性裸体啪啪无遮挡网站| 999精品在线视频| 国产亚洲欧美精品永久| 国产白丝娇喘喷水9色精品| 狂野欧美激情性bbbbbb| 国产精品秋霞免费鲁丝片| 国产成人午夜福利电影在线观看| 国语对白做爰xxxⅹ性视频网站| 91国产中文字幕| 日本欧美视频一区| 啦啦啦视频在线资源免费观看| 亚洲av欧美aⅴ国产| www.自偷自拍.com| 青春草亚洲视频在线观看| 免费女性裸体啪啪无遮挡网站| 黄色 视频免费看| av网站在线播放免费| 国产白丝娇喘喷水9色精品| 少妇熟女欧美另类| 免费不卡的大黄色大毛片视频在线观看| 精品少妇黑人巨大在线播放| 9热在线视频观看99| 极品人妻少妇av视频| 国产精品不卡视频一区二区| 亚洲四区av| 99热国产这里只有精品6| 国产精品国产av在线观看| 最近手机中文字幕大全| 国产精品一国产av| 精品一品国产午夜福利视频| 亚洲婷婷狠狠爱综合网| 亚洲美女黄色视频免费看| 久久国产精品男人的天堂亚洲| 在线观看免费高清a一片| 国产极品天堂在线| 巨乳人妻的诱惑在线观看| 国产成人a∨麻豆精品| 日产精品乱码卡一卡2卡三| 免费少妇av软件| 国产成人a∨麻豆精品| 免费人妻精品一区二区三区视频| 国产色婷婷99| 午夜日韩欧美国产| 香蕉国产在线看| 美女中出高潮动态图| 亚洲国产欧美网| 制服诱惑二区| 国产av一区二区精品久久| 夜夜骑夜夜射夜夜干| 成人黄色视频免费在线看| 国产av一区二区精品久久| 又黄又粗又硬又大视频| 观看av在线不卡| 蜜桃在线观看..| 亚洲av免费高清在线观看| 最新中文字幕久久久久| 免费在线观看完整版高清| 亚洲国产精品成人久久小说| 99热网站在线观看| 成年人免费黄色播放视频| 99热全是精品| 精品午夜福利在线看|