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

    A New Species of Odorous Frog Genus Odorrana (Anura,Ranidae) from Southern Guizhou Province,China

    2021-12-30 08:24:20TaoLUOSiweiWANGNingXIAOYaliWANGandJiangZHOU
    Asian Herpetological Research 2021年4期
    關(guān)鍵詞:荔波

    Tao LUO ,Siwei WANG ,Ning XIAO ,Yali WANG and Jiang ZHOU

    1 School of Karst Science,Guizhou Normal University,Guiyang 550001,Guizhou,China

    2 School of Life Sciences,Guizhou Normal University,Guiyang 550001,Guizhou,China

    3 Guiyang Nursing Vocational College,Guiyang 550003,Guizhou,China

    Abstract We describe Odorrana liboensis sp.nov., a new species from the Maolan National Nature Reserve,Libo County,Guizhou Province,China.Phylogenetic analyses based on DNA sequences of the mitochondrial 12S rRNA,16S rRNA,and ND2 genes supported the new species as an independent lineage,closely related to O.lipuensis.The uncorrected genetic distances between the 12S rRNA and 16S rRNA in the new species and its closest congener, O.lipuensis, were 6.06% and 5.19%,respectively.The new species is distinguished from its congeners by a combination of the following morphological characters:(1) having medium body size,with the snout-vent length (SVL) of adult females approximately 1.2 times as long as males at 56.9 ±1.0 (55.8–58.2 mm,n=9) in females and 48.7 ± 1.2 (47.1–49.9 mm,n=5) in males;(2) head length greater than width in males and females;(3) tympanum distinctly visible,greater than one-half the diameter of the eye;(4) eyes big and prominent,width of upper eyelid (UEW) approximately 3/4 of eye diameter (ED);(5) dorsolateral folds absent;(6)two metacarpal tubercles;(7) relative finger lengths:II

    Keywords Taxonomy,morphology,Odorrana,Odorrana liboensis sp.nov., karst cave,Guizhou

    1.Introduction

    The genus

    Odorrana

    was first recognized by Fei

    et al

    .(1990)with the type species

    Odorrana margaretae

    (Liu,1950).Despite previous controversy over classification (Dubois,1992;Matsui,1994),molecular phylogenetic analyses suggested that

    Odorrana

    forms a well-supported monophyletic group (Matsui

    et al

    .,2005;Ngo

    et al

    .,2006;Cai

    et al.

    ,2007;Che

    et al

    .,2007;Stuart,2008;Chen

    et al

    .,2013;Li

    et al

    .,2018a).According to Frost (2021),

    Odorrana

    contains at least 61 recognized species and is widely distributed in montane streams and rivers in the subtropical and tropical regions of East and Southeast Asia (Fei

    et al

    .,

    2012;AmphibiaChina,2021;Frost,2021).The range of species in this genus covers the Ryukyu Archipelago,southern China,northeastern India,and the Thai-Malay Peninsula,and extends southwards to Sumatra and Borneo (Frost,2021).All species in the genus are associated with mountain streams except for

    O.wuchuanensis

    (Xu,1983) and

    O.lipuensis

    Mo,Chen,Wu,Zhang &Zhou,2015,which occur in dark caves.Among the

    Odorrana

    species,39 occur in China and 24 of these are endemic to China (Fei

    et al

    .,

    2012;AmphibiaChina,2021;Frost,2021;Chen

    et al

    .,

    2020;Shen

    et al

    .,

    2020;Liu

    et al

    .,2021).Fei

    et al

    .(1990) established

    Odorrana

    ,and later Ye and Fei (2001) suggested four species groups (

    O.andersonii

    ,

    O.kuangwuensis

    ,

    O.schmacheri

    ,

    and

    O.livida

    ) within

    Odorrana

    based on a morphological phylogenetic study.Fei

    et al

    .(2005)established two subgenera (

    Odorrana

    and

    Bamburana

    ) within

    Odorrana

    and recognized

    Odorrana versabilis

    as the type species of

    Bamburana

    .Molecular phylogenetic studies support the division of species groups within

    Odorrana

    but not the two subgenera (Che

    et al

    .,

    2007).Subsequently,Fei

    et al

    .(2009) divided the Chinese

    Odorrana

    species into six species groups (

    O.tormota

    ,

    O.taiwaniana

    ,

    O.graminea

    ,

    O.margaretae

    ,

    O.schmackeri

    ,and

    O.andersonii

    ;

    O.taiwaniana

    is a synonym for

    O.swinhoana

    ) based on morphological characters.These divisions have been accepted by other researchers and applied to the

    Odorrana

    species distributed outside China (Pham

    et al

    .,2016;Li

    et al

    .,2018a).However,recent studies have rejected the monophyly of the

    O.margaretae

    ,

    O.schmackeri

    ,and

    O.andersonii

    species groups (Chen

    et al

    .,2013).The phylogenetic relationships between these species groups,the species included in the species groups,and the species group affiliations of new species published in recent years are unclear.The montane river and stream habitats of most

    Odorrana

    may promote intraspecific divergence through geographical isolation,especially given the wide geographical range of the complex (Chen

    et al

    .,2013;Li

    et al

    .,2015).Indeed,many cryptic species of

    Odorrana

    have been discovered.For example,25 new species of

    Odorrana

    have been described since 2005 (See species list of the genus

    Odorrana

    in Frost(2021)and AmphibiaChina(2021)).This indicates that further discoveries or potential taxonomic changes within

    Odorrana

    are likely.From 2016 to 2018,14 specimens of an unknown species of

    Odorrana

    were collected inside a completely dark karst cave in the Maolan National Nature Reserve,Libo County,Guizhou Province,China.Morphologically,these specimens most closely resemble

    O.lipuensis

    and

    O.kweichowensis

    Li,Xu,Lv,Jiang,Wei &Wang,2018 (Mo

    et al

    .,2015;Li

    et al

    .,2018a),but differ from

    O.lipuensis

    ,

    O.kweichowensis

    ,and all other

    Odorrana

    from China and adjoining countries.They inhabit a dark karst cave environment similar to that of

    O.

    wuchuanensis

    and

    O.lipuensis

    but are morphologically very different from

    O.

    wuchuanensis

    and

    O.lipuensis

    .To distinguish these specimens,we conducted phylogenetic analyses based on mitochondrial DNA and morphological comparisons.All of the analyses consistently indicated that the specimens from Maolan National Nature Reserve are a new taxon.We describe this taxon here as a new species.

    2.Materials and Methods

    2.1.Sampling

    A total of 72 specimens were collected in this study.Fourteen specimens were of the undescribed species from Maolan National Nature Reserve in Libo County,Guizhou Province;35 were

    O.kweichowensis

    from Lengshuihe Nature Reserve in Jinsha County,Guizhou Province;four were

    O.yizhangensis

    Fei,Ye &Jiang,2007 from Kuankoushui National Nature Reserve in Suiyang County,Guizhou Province,and Yuntai Mountain,Shibing County;one was

    O.huanggangensis

    Chen,Zhou &Zheng,2010 from Yueliangshan Nature Reserve in Conjiang County,Guizhou Province;one was

    O.kweichowensis

    from Maolan National Nature Reserve in Libo County,Guizhou Province;and 17 were

    O.wuchuanensis

    from Maolan National Nature Reserve,Libo County,Guizhou Province (Figure 1).All of the specimens were fixed in 10%buffered formalin and later transferred to 75% ethanol for preservation.The muscles used for molecular analysis were preserved in 95% alcohol at ?20 °C.All of the specimens were deposited in Guizhou Normal University (GZNU),Guiyang City,Guizhou Province,China.

    2.2.DNA Extraction,PCR amplification,and sequencing

    Genomic DNA was extracted from muscular tissue using a DNA extraction kit from Tiangen Biotech Co.,Ltd.(Beijing).All samples were sequenced for three mitochondrial genes:partial 12S ribosomal RNA gene (12S rRNA),16S ribosomal RNA gene (16S rRNA),and NADH dehydrogenase subunit 2 (ND2).The primers used for 12S rRNA were 12SF (5’-GGTTTGRTCCTRGCCTTAC-3’)and 12SR (5’-CCATGTTACGACTTGCCTCT-3’) following Chen

    et al

    .(2013),the primers used for 16S rRNA were 16SF (5’-ACGAGCCTAGTGATAGCTGGTT-3’) and 16SR (5’-CGGTCTGAACTCAGATCACGT-3’) following Chen

    et al

    .(2013),and the primers used for ND2 were Gln-LND2 (5’-CCCTTTGCACTTCCTTTATGC-3’) and Ala-HND2 (5’-GGCCTGAGTTGCATTCATG-3’) following Li

    et al

    .(2015).PCR amplifications were performed in a 25 μL reaction volume with the following cycling conditions:an initial denaturing step at 95 °C for five min;36 cycles of denaturing at 95 °C for 40 s,annealing at 52 °C (for 12S rRNA)/49 °C (for 16S rRNA)/60 °C (for ND2) for 40 s and extending at 72 °C for 1 min,and a final extending step of 72 °C for 10 min.The PCR products were sequenced on an ABI Prism 3730 automated DNA sequencer in Chengdu TSING KE Biological Technology Co.Ltd.(Chengdu,China).All of the sequences have been deposited in GenBank (Table 1).Some homologous DNA sequences of voucher specimens of related species were downloaded from GenBank and incorporated into the phylogenetic analyses.

    2.3.Phylogenetic analyses

    We used a total of 187 sequences (including 72 12S rRNA,79 16S rRNA,and 36 ND2 gene sequences) for molecular analyses.Three mitochondrial genes from 11 muscle tissues were sequenced and 154 sequences downloaded from GenBank from 48 species of the genus

    Odorrana

    were used.These included the undescribed species from China that was named in this study (Figure 1).A total of 26 sequences were downloaded from GenBank as out-groups(Chen

    et al

    .,2013).Detailed information on these materials is shown in Table 1.All of the sequences were assembled and aligned using the MUSCLE (Edgar,2004) module in MEGA 7.0 (Kumar

    et al

    .,2016)with default settings.Alignments were checked by eye and revised manually if necessary.Trimming with the gaps partially deleted was performed in GBLOCKS 0.91b (Castresana,2000).Phylogenetic trees were constructed with both maximum likelihood (ML) and Bayesian inference (BI).The ML was conducted in IQ-TREE(Nguyen

    et al

    .,2015) with 2000 ultrafast bootstrapping (Hoang

    et al

    .,2018) and was performed until a correlation coefficient of at least 0.99 was reached.The BI was performed in MrBayes 3.2.1 (Ronquist

    et al

    .,2012),and the best-fit model was obtained by the Bayesian inference criteria (BIC) computed with PartitionFinder 2 (Lanfear

    et al

    .,2016).For this analysis,12S rRNA,16S rRNA,and ND2 genes were defined.The analysis suggested that the best partition scheme was 12S rRNA/16S rRNA/ND2 genes.We selected GTR+I+G as the best model for 12S rRNA and 16S rRNA and the TIM+I+G model as the best model for the ND2 gene.Two independent runs were conducted in BI analysis,each of which was performed for 2,000,000 generations and sampled every 1000 generations.The first 25% of the samples were discarded as burn-in,resulting in a potential scale reduction factor (PSRF) of <0.01.Nodes in the trees were considered well supported when Bayesian posterior probabilities (BPP) were ≥0.95 and the ML ultrafast bootstrap value(UFB) was ≥95%.Uncorrected

    p

    -distances (1000 replicates) based on 12S rRNA and 16S rRNA were calculated in MEGA 7.0 (Kumar

    et al

    .,2016).

    Figure 1 Sampling collection localities and distribution of the Odorrana liboensis sp.nov.,O.lipuensis,O.wuchuanensis,O.kweichowensis,O.huanggangensis,and O.yizhangensis in southwest China.1.Maolan National Nature Reserve,Libo County,Guizhou,China.2.Lipu County,Guangxi,China.3.Jiarong Town,Libo County,Guizhou,China.4.Maolan National Nature Reserve,Libo County,Guizhou,China.5.Lengshuihe Nature Reserve,Jinsha County,Guizhou,China.6.Yueliangshan Nature Reserve,Congjiang County,Guizhou,China.7.Yuntai Mountains,Shibing County,Guizhou,China.8.Kuankuoshui National Nature Reserve,Suiyang County,Guizhou,China.The base maps are from Standard Map Service website (http://bzdt.ch.mnr.gov.cn/index.html).

    Table 1 Localities,voucher information,and GenBank numbers for all samples used in this study.

    2.4.Morphological analysis

    Morphometric data were taken from 74 well-preserved adult specimens (voucher information in Table 2 and Table S1).Measurements were recorded to the nearest 0.1 mm with digital calipers by Tao Luo following Fei

    et al

    .(2009) and Li

    et al

    .(2018a).A total of 24 morphological characteristics were measured.These measurements were as follows:

    SVL

    =snout-vent length(distance from the tip of the snout to the posterior edge of the vent),

    HDL

    =head length (distance from the tip of the snout to the articulation of jaw),

    HDW

    =maximum head width (greatest width between the left and right articulations of jaw),

    SL

    =snout length(distance from the tip of the snout to the anterior corner of the eye),

    IND

    =internasal distance (minimum distance between the inner margins of the external nares),

    NED

    =nasal to eye distance(distance between the nasal and the anterior corner of the eye),

    NSD

    =nasal to snout distance (distance between the nasal the posterior edge of the vent),

    IOD

    =interorbital distance (minimum distance between the inner edges of the upper eyelids),

    IFE

    =distance between anterior corner of eye,

    IAE

    =distance between posterior corner of eyes,

    UEW

    =upper eyelid width (greatest width of the upper eyelid margins measured perpendicular to the anterior-posterior axis),

    ED

    =eye diameter (distance from the anterior corner to the posterior corner of the eye),

    TD

    =tympanum diameter,

    LAHL

    =length of lower arm and hand (distance from the elbow to the distal end of the finger IV),

    LAN

    =width of lower arm,

    HAL

    =hand length(distance from the posterior end of the inner metacarpal tubercle to the distal tip of Finger IV),

    HLL=

    hindlimb length (maximum length from the vent to the distal tip of the toe IV),

    TL=

    tibia length (distance from knee to tarsus),

    TW=

    maximal tibia width,

    TFL=

    length of foot and tarsus (distance from the tibiotarsal articulation to the distal end of the toe IV),

    FOL=

    foot length (from the base of inner metatarsal tubercle to the tip of fourth toe),

    TED

    =tympanum-eye distance (from anterior edge of tympanum to posterior corner of the eye),

    FDW

    =finger disk width (width at the widest part of the disk of finger III),and

    DPW

    =distal phalanx width (maximal width of the distal phalanx of finger III).To reduce the impact of allometry,a size-corrected value from the ratio of each character to SVL was calculated for the following morphometric analyses.Principal component analysis (PCAs) of size-corrected variables and simple bivariate scatterplots was used to explore and reflect the morphometric differences between the new species and

    O.lipuensis

    .Mann-Whitney

    U

    tests were used to test the significance of differences on morphometric characters between the new species and

    O.lipuensis

    ,

    O.kweichowensis

    and

    O.wuchuanensis

    .Mann-Whitney

    U

    tests were also conducted to test the morphometric differences between the males and the females of the new species.The statistical analyses were performed using SPSS 21.0 (SPSS,Inc.,Chicago.IL,USA),and differences were considered to be significant at

    P

    <0.05.

    Sex was determined by direct observation of calling behavior and the presence of internal vocal sac openings for males,as well as the presence of eggs on the abdomen for females.The presence or absence of nuptial pads/spines was examined by optical microscopy.

    We compared the morphological characters of the new taxon with other species of

    Odorrana

    .Comparative data were obtained from the literature for 61 species of

    Odorrana

    (Table 3).For comparison,we examined the type and/or topotype materials for

    O.lipuensis

    ,

    O.kweichowensis

    ,and

    O.

    wuchuanensis

    (Appendix I).

    3.Results

    3.1.Phylogenetic analyses and genetic divergence

    The maximum likelihood (ML) and Bayesian inference (BI)phylogenetic trees were constructed based on concatenated DNA sequences of the mitochondrial 12S rRNA (752 bp),16S rRNA (1006 bp),and ND2 (1030 bp) genes with a total length of 2788 bp.ML and BI analyses resulted in a largely identical topology (Figure 2).The new taxon is a sister taxon to

    O.lipuensis

    with high node support values (1.00 in BI and 100% in ML;Figure 2).The smallest

    p

    -distance divergences between the new lineage and other species of

    Odorrana

    were 6.06% in 12S rRNA (between new taxon and

    O.lipuensis

    ) and 3.90% in 16S rRNA (between new taxon and

    O.geminata

    ),which were at the same divergence level as those among recognized congeners (1.94% and 2.79% in 12S rRNA between

    O.

    wuchuanensis

    vs.O.mutschmanni

    and

    O.kweichowensis

    vs.O.schmackeri

    ,respectively;and 1.30% and 1.95%in 16S rRNA between

    O.lungshengensis

    vs.O.yizhanggensis

    and

    O.nanjiangensis

    vs.O.hejiangensis

    ,respectively),indicating that the new taxon represents an independent evolutionary lineage(Tables S2–S3).

    3.2.Morphological analyses

    The Mann-Whitney

    U

    tests indicated that males of the new species were significantly different from

    O.lipuensis

    and

    O.kweichowensis

    in many morphometric characters (all

    P

    values <0.05;Table 4).PCA extracted three and four principal component factors with Eigenvalues greater than two in males and females,respectively(Table S4).The first two principal components explained 61.34% and 63.28%,of the total variation in males and females,respectively.These differences were mainly influenced by limb and head characteristics (Table S4).The new taxon was distinctly separated from

    O.lipuensis

    and

    O.kweichowensis

    on the two-dimensional plots of PC1

    vs.

    PC2,in both males and females(Figure 3).

    3.3.Taxonomic account

    sp.nov.(Tables 2 and 4,Figures 4–5)

    Holotype.

    GZNU20180608004 (Figure 4),adult male,collected by Tao LUO on June 8,2018 from Maolan National Nature Reserve (25.481711° N,108.078003° E,ca.715 m a.s.l.),Libo County,Guizhou Province,China.

    Figure 2 Phylogenetic tree based on mitochondrial 12S rRNA+16S rRNA+ND2 genes.In the phylogenetic tree,ultrafast bootstrap supports (UFB) from ML analyses and Bayesian posterior probabilities (BPP) from BI analyses were noted beside nodes.The scale bar represents 0.05 nucleotide substitutions per site.Numbers at the tips of branches correspond to the ID numbers in Table 1.

    Figure 3 Plots of the first principal component (PC 1) versus the second (PC 2) for Odorrana liboensis sp.nov.and O.lipuensis from a principal component analysis.A,male.B,female.

    Paratypes.

    Thirteen individuals collected at the same locality as the holotype.GZNU20180608002,GZNU20180608003,GZNU20180608006,GZNU20180608007,GZNU20180608009,GZNU20180608010,GZNU20180608011 and GZNU20180608012,adult males collected by Tao Luo on June 8;GZNU20180815001(Figure 5),adult females collected by Xiang ZENG on August 15,2018;GZNU20160802001,GZNU20160802002,GZNU20160802003,and GZNU20160729009,adult females collected by Tao Luo on August 2,2016.

    Etymology.

    The specific epithet“

    liboensis

    ”is in reference to the type locality of the new species:Libo County,Guizhou Province,China.We propose the common English name“Libo Odorous Frog”and Chinese name“Li Bo Chou Wa (荔波臭蛙)”.

    Differential diagnosis.

    Odorrana liboensis

    sp.nov.

    is assigned to genus

    Odorrana

    based upon molecular phylogenetic analyses and the following morphological characters:(1) dorsum is mostly green;(2) tips of digits dilated,tapering,disks with circum-marginal grooves or lateroventral grooves,and vertical diameter longer than horizontal diameter in the disks;(3)supernumerary tubercle below the base of fingers III and IV;(4) feet fully webbed to disks,without tarsal fold;(5) the first finger thick and nuptial pad distinct and (6) dorsal skin smooth,dorsolateral folds absent or fine (Fei

    et al

    .,2009;Li

    et al

    .,2018a).

    Odorrana liboensis

    sp.nov.

    is distinguished from congeners by a combination of the following characters:(1) having medium body size,with the snout-vent length (SVL) of adult females approximately 1.2 times as long as that of males at 56.9 ± 1.0(55.8–58.2 mm,

    n=

    9) in females and 48.7 ± 1.2 (47.1–49.9 mm,

    n=

    5) in males;(2) head length greater than width in males and females;(3) tympanum distinctly visible,greater than one-half the diameter of the eye;(4) eyes big and prominent,with the width of the upper eyelid (UEW) approximately 3/4 of the eye diameter (ED);(5) dorsolateral folds absent;(6) two metacarpal tubercles;(7) relative finger lengths:II Description of holotype.

    GZNU 20180608004 (Figure 4),adult male.Medium body size,SVL 49.8 mm;habitus slender;head length slightly larger than head width (HDL=17.6 mm;HDW=15.6 mm;HDL/HDW 1.12);snout short,rounded,and projecting beyond the lower jaw in dorsal view,longer than the diameter of the eye (SL/ED 1.30);nostril rounded,distinct,closer to the snout tip than the eye (NEL/NSL 0.52);internasal distance less than interorbital distance (IND/IOD 0.91) and greater than upper eyelid width (IND/UEW 1.06);pineal body invisible;tympanum distinct,rounded,75% of eye diameter,depressed relative to the skin of the temporal region,tympanic rim slightly elevated relative to tympanum;vomerine teethmoderately developed,on two oblique ridges;tongue cordiform,deeply notched posteriorly;and eyes large,slightly protuberant in dorsal view,eye diameter 26.70% of head length,pupils transverse;supratympanic fold absent.

    Table 3 References for morphological characters for congeners of the genus Odorrana.

    Forelimbs are slender and comparatively short,the length of the lower arm and hand is 41.97% of SVL;fingers slender,relative finger lengths:I

    Hindlimbs are slender (HLL/SVL 1.86);heels overlapping when thighs are positioned at right angles to the body;tibiotarsal articulation reaching the between eye to nostril when leg stretched forward;foot length less than tibia length(FL/TL 0.87);relative toe lengths I

    Figure 4 Morphological features of the live adult male holotype GZNU20180608004 of Odorrana liboensis sp.nov. and O.wuchuanensis.(A) Dorsal view;(B) Ventral view;(C) Dorsolateral view;(D) Tadpole;(E) Type locality of Odorrana liboensis sp.nov. in June,2018 (dry season);(F) Dorsolateral view of O.wuchuanensis in life.

    Skin is smooth on upper surfaces;venter smooth;tiny spinules on lateral body,upper edge of lid,temporal region,and anterior and posterior edge of tympanum;weak supratympanic fold from the posterior edge of the eye to the posterior edge of the tympanum;and dorsolateral fold absent.

    Coloration of holotype in life

    (Figure 4).Dorsal parts of head and dorsum,flank,forelimb,thigh,tibia,and foot are grass-green with irregular brown mottling;the throat,chest,and abdomen lack black spots;dorsal surfaces of limbs have indistinct brownish-black bands,the thighs with four brownishblack bands and tibias with three;ventral surface of the limbs is pink;the iris is black,surrounded by a gold-green network;and the tympanum is dark brown.

    Preserved holotype coloration

    .After preservation in 75%ethanol,the dorsal surface of the body coloration changed to dark brown grey;the dorsal surface of the head coloration changed to dark grey;the transverse bands on limbs and digits were not distinct and the coloration changed to lighter colors;the throat was light grayish yellow;the chest was light ash black;the belly was light gray-white;the posterior of ventral surface of body,inner thigh,and upper part of the tibia were creamy yellow;the palms and metatarsal tubercles were ash black;and the ventral surface of the forelimbs and hindlimbs were creamy yellow with brown mottling.

    Table 4 Morphological comparison of Odorrana liboensis sp.nov.(OB),O.lipuensis (OP),O.kweichowensis (OK),and O.wuchuanensis (OW).All units are in mm.P-values are at 95% significance.Morphometric characters are explained in the methods section.BM and BF are the abbreviations for male and female from Odorrana liboensis sp.nov.

    (Continued Table 4)

    Variations.

    The basic statistics for measurements are presented in Table 2.In life,all paratypes matched the overall basic morphological characters of the holotype.Females (SVL 56.9 ±1.0 mm,

    n=

    5) have larger body size than males (SVL 48.7 ± 1.2 mm,

    n=

    9),with the SVL in females approximately 1.2 times that in males,but the relative sizes of ED and TD are obviously larger in males than in females (Table 4);the dorsal surfaces of limbs with indistinct brownish-black bands (Figure 5).

    Sexual dimorphism.

    Adult females have an SVL of 55.8–58.2 mm,larger than adult males,with an SVL of 47.1–49.9 mm.Adult males lack vocal sacs.In breeding,light white nuptial pads are present on finger I in males.

    Comparisons.

    Comparative data of

    Odorrana liboensis

    sp.nov.

    with 61 recognized congeners of

    Odorrana

    are given in Table S5.By having medium body size

    (minimum SVL >47.0 mm in males),

    Odorrana liboensis

    sp.nov.

    differs from

    O.absita

    ,

    O.anlungensis

    ,

    O.gigatympana

    ,

    O.huanggangensis

    ,

    O.khalam

    ,

    O.kweichowensis

    ,

    O.monjerai

    ,

    O.morafkai

    ,

    O.nasica

    ,

    O.orba

    ,and

    O.tormota

    (

    vs.

    maximum SVL <47.0 mm in males).By having medium body size (maximum SVL <50.0 mm in males),

    Odorrana liboensis

    sp.nov.

    differs from

    O.amamiensis

    ,

    O.andersonii

    ,

    O.aureola

    ,

    O.bacboensis

    ,

    O.cangyuanensis

    ,

    O.chapaensis

    ,

    O.geminata

    ,

    O.grahami

    ,

    O.hainanensis

    ,

    O.hosii

    ,

    O.ishikawae

    ,

    O.jingdongensis

    ,

    O.junlianensis

    ,

    O.kuangwuensis

    ,

    O.leporipes

    ,

    O.livida

    ,

    O.lungshengensis

    ,

    O.macrotympana

    ,

    O.margaretae

    ,

    O.mutschmanni

    ,

    O.nanjiangensis

    ,

    O.nasuta

    ,

    O.sinica

    ,

    O.splendida

    ,

    O.supranarina

    ,

    O.swinhoana

    ,

    O.tiannanensis

    ,

    O.trankieni

    ,

    O.versabilis

    ,

    O.wuchuanensis

    ,and

    O.zhaoi

    (

    vs.

    minimum SVL >50.0 mm).By having medium body size

    (

    vs.

    maximum SVL <59.0 mm in females),

    Odorrana liboensis

    sp.nov.

    differs from

    O.amamiensis

    ,

    O.andersonii

    ,

    O.anlungensis

    ,

    O.aureola

    ,

    O.bacboensis

    ,

    O.banaorum

    ,

    O.bolavensis

    ,

    O.chapaensis

    ,

    O.chloronota

    ,

    O.dulongensis

    ,

    O.fengkaiensis

    ,

    O.geminata

    ,

    O.grahami

    ,

    O.hainanensis

    ,

    O.hejiangensis

    ,

    O.hosii

    ,

    O.huanggangensis

    ,

    O.ishikawae

    ,

    O.jingdongensis

    ,

    O.junlianensis

    ,

    O.kuangwuensis

    ,

    O.kweichowensis

    ,

    O.livida

    ,

    O.lungshengensis

    ,

    O.macrotympana

    ,

    O.margaretae

    ,

    O.monjerai

    ,

    O.morafkai

    ,

    O.mutschmanni

    ,

    O.nanjiangensis

    ,

    O.narina

    ,

    O.nasuta

    ,

    O.orba

    ,

    O.rotodora

    ,

    O.schmackeri

    ,

    O.splendida

    ,

    O.supranarina

    ,

    O.swinhoana

    ,

    O.tianmuii

    ,

    O.tiannanensis

    ,

    O.tormota

    ,

    O.versabilis

    ,

    O.wuchuanensis

    ,

    O.yentuensi

    ,

    O.yizhangensis

    ,

    and

    O.ichangensis

    (

    vs.

    maximum SVL >59.0 mm).By lacking dorsolateral folds,

    Odorrana liboensis

    sp.nov.

    differs from

    O.absita

    ,

    O.amamiensis

    ,

    O.banaorum

    ,

    O.bolavensis

    ,

    O.exiliversabilis

    ,

    O.gigatympana

    ,

    O.graminea

    ,

    O.hosii

    ,

    O.indeprensa

    ,

    O.khalam

    ,

    O.leporipes

    ,

    O.livida

    ,

    O.monjerai

    ,

    O.narina

    ,

    O.nasica

    ,

    O.nasuta

    ,

    O.orba

    ,

    O.supranarina

    ,

    O.tormota

    ,

    O.trankieni

    ,

    O.utsunomiyaorum

    ,

    O.yentuensis

    ,and

    O.zhaoi

    (

    vs.

    present).

    Figure 5 Morphological features of the live adult female GZNU20180815001 of Odorrana liboensis sp.nov. (A) Dorsal view;(B) Dorsolateral view;(C)Ventral view;(D) Ventral view of right foot;(E) Ventral view of right hand;(F) The entrance habitat of the karst cave at the type locality of Odorrana liboensis sp.nov.;(G) Water outlet in the cave.

    By the tibiotarsal articulation reaching to between the eye and the nostril when the leg is stretched forward,

    Odorrana liboensis

    sp.nov.

    differs from

    O.bacboensis

    ,

    O.jingdongensis

    ,

    O.lungshengensis

    ,

    O.margaretae

    ,

    O.mutschmanni

    ,

    O.nanjiangensis

    ,

    O.narina

    ,

    O.orba

    ,

    O.sinica

    ,

    O.swinhoana

    ,

    O.tormota

    ,and

    O.yizhangensis

    (

    vs.

    reaching the tip of the snout),from

    O.nasica

    and

    O.nasuta

    (

    vs.

    reaching the tip of the snout or a little beyond),from

    O.hainanensis

    (

    vs.

    reaching the tip of the snout or the anterior corner of eye),from

    O.junlianensis

    (

    vs.

    reaching the tip of the snout or between the nostril and the snout),from

    O.cangyuanensis

    ,

    O.exiliversabilis

    ,

    O.fengkaiensis

    ,

    O.gigatympana

    ,

    O.grahami

    ,

    O.graminea

    ,

    O.tiannanensis

    ,

    O.versabilis

    ,

    O.yentuensis

    ,and

    O.zhaoi

    (

    vs.

    reaching to or beyond the tip of the snout),from

    O.amamiensis

    (

    vs.

    reaching far beyond the tip of the snout),from

    O.amamiensis

    ,

    O.anlungensis

    ,

    O.huanggangensis

    ,

    O.kuangwuensis

    ,

    O.macrotympana

    ,

    O.wuchuanensis

    ,

    and

    O.ichangensis

    (

    vs.

    reaching the nostril or beyond the tip of the snout),from

    O.lipuensis

    ,

    O.splendida

    ,and

    O.supranarina

    (

    vs.

    reaching the anterior corner of the eye),from

    O.rotodora

    and

    O.trankieni

    (

    vs.

    reaching beyond the eye),and from

    O.utsunomiyaorum

    (

    vs.

    reaching between the anterior corner of the eye and the nostril).By lacking black bars on the lips,

    Odorrana liboensis

    sp.nov.

    differs from

    O.andersonii

    ,

    O.anlungensis

    ,

    O.bacboensis

    ,

    O.bolavensis

    ,

    O.chapaensis

    ,

    O.dulongensis

    ,

    O.fengkaiensis

    ,

    O.geminata

    ,

    O.grahami

    ,

    O.hainanensis

    ,

    O.hejiangensis

    ,

    O.huanggangensis

    ,

    O.ishikawae

    ,

    O.junlianensis

    ,

    O.kuangwuensis

    ,

    O.kweichowensis

    ,

    O.lungshengensis

    ,

    O.margaretae

    ,

    O.mutschmanni

    ,

    O.nanjiangensis

    ,

    O.schmackeri

    ,

    O.splendida

    ,

    O.supranarina

    ,

    O.tianmuii

    ,

    O.tiannanensis

    ,

    O.wuchuanensis

    ,

    O.yizhangensis

    ,

    and

    O.ichangensis

    (

    vs.

    the presence of black bars).By the absence of vocal sacs in males,

    Odorrana liboensis

    sp.nov.

    differs from

    O.absita

    ,

    O.amamiensis

    ,

    O.andersonii

    ,

    O.anlungensis

    ,

    O.aureola

    ,

    O.bacboensis

    ,

    O.banaorum

    ,

    O.bolavensis

    ,

    O.cangyuanensis

    ,

    O.chapaensis

    ,

    O.chloronota

    ,

    O.dulongensis

    ,

    O.exiliversabilis

    ,

    O.fengkaiensis

    ,

    O.geminata

    ,

    O.gigatympana

    ,

    O.grahami

    ,

    O.graminea

    ,

    O.hainanensis

    ,

    O.hejiangensis

    ,

    O.huanggangensis

    ,

    O.indeprensa

    ,

    O.ishikawae

    ,

    O.jingdongensis

    ,

    O.junlianensis

    ,

    O.khalam

    ,

    O.kweichowensis

    ,

    O.lungshengensis

    ,

    O.macrotympana

    ,

    O.morafkai

    ,

    O.nanjiangensis

    ,

    O.nasica

    ,

    O.nasuta

    ,

    O.orba

    ,

    O.sinica

    ,

    O.swinhoana

    ,

    O.tianmuii

    ,

    O.tiannanensis

    ,

    O.tormota

    ,

    O.trankieni

    ,

    O.utsunomiyaorum

    ,

    O.versabilis

    ,

    O.yentuensis

    ,

    O.yizhangensis

    ,

    O.ichangensis

    ,and

    O.zhaoi

    (

    vs.

    the presence of vocal sacs).By the relative lengths of fingers IOdorrana liboensis

    sp.nov.

    differs from

    O.lipuensis

    and

    O.dulongensis

    (

    vs.

    I=II O.absita

    ,

    O.amamiensis

    ,and

    O.anlungensis

    (

    vs.

    II O.andersonii

    (

    vs.

    I O.aureola

    ,

    O.bacboensis

    ,

    O.banaorum

    ,

    O.bolavensis

    ,

    O.chloronota

    ,

    O.fengkaiensis

    ,

    O.gigatympana

    ,

    O.grahami

    ,

    O.jingdongensis

    ,

    O.junlianensis

    ,

    O.kweichowensis

    ,

    O.margaretae

    ,

    O.mutschmanni

    ,

    O.orba

    ,

    O.schmackeri

    ,

    O.swinhoana

    ,

    O.tiannanensis

    ,

    O.wuchuanensis

    ,and

    O.yentuensis

    (

    vs.

    II O.cangyuanensis

    (

    vs.

    I O.exiliversabilis

    ,

    O.hainanensis

    ,

    O.hejiangensis

    ,

    O.kuangwuensis

    ,

    O.nasica

    ,

    O.nasuta

    ,and

    O.versabilis

    (

    vs.

    IIO.huanggangensis

    (

    vs.

    I≤II=IV O.splendida

    ,

    O.supranarina

    ,and

    O.trankieni

    (

    vs.

    IVO.zhaoi

    (

    vs.

    I=IVOdorrana liboensis

    sp.nov.

    differs from

    O.andersonii

    ,

    O.fengkaiensis

    ,

    O.grahami

    ,

    O.graminea

    ,

    O.hainanensis

    ,

    O.huanggangensis

    ,

    O.jingdongensis

    ,

    O.junlianensis

    ,

    O.kweichowensis

    ,

    O.lungshengensis

    ,

    O.margaretae

    ,

    O.tianmuii

    ,

    O.yizhangensis

    ,and

    O.zhaoi

    (

    vs.

    presence of white pectoral spinules).By having two metacarpal tubercles on the base of the hand,

    Odorrana liboensis

    sp.nov.

    differs from

    O.absita

    ,

    O.aureola

    ,and

    O.bacboensis

    (

    vs.

    lacking metacarpal tubercles),from

    O.cangyuanensis

    ,

    O.dulongensis

    ,

    O.lungshengensis

    ,

    O.nanjiangensis

    ,

    O.orba

    ,

    O.schmackeri

    ,

    O.tiannanensis

    ,and

    O.yizhangensis

    (

    vs.

    one metacarpal tubercle),from

    O.andersonii

    ,

    O.exiliversabilis

    ,

    O.fengkaiensis

    ,

    O.hainanensis

    ,

    O.hejiangensis

    ,

    O.huanggangensis

    ,

    O.jingdongensis

    ,

    O.junlianensis

    ,

    O.nasica

    ,

    O.nasuta

    ,

    O.tianmuii

    ,

    O.tormota

    ,and

    O.ichangensis

    (

    vs.

    three metacarpal tubercles).By having one metacarpal tubercle on the base of the metatarsus,

    Odorrana liboensis

    sp.nov.

    differs from

    O.absita

    ,

    O.amamiensis

    ,

    O.exiliversabilis

    ,

    O.nasica

    ,and

    O.nasuta

    (

    vs.

    two metatarsal tubercles)

    .

    The congeners

    O.kweichowensis

    and

    O.wuchuanensis

    have a sympatric distribution with

    Odorrana liboensis

    sp.nov.

    The new species can be distinguished from these species by a series of morphological characters as follows.The new species differs from

    O.kweichowensis

    by having a larger body size (adult males with a length of 47.1–49.9 mm

    vs.

    adult

    O.kweichowensis

    males with a length of 36.2–43.3 mm) and having a lower ratio of IND/IOD at 0.89 in males and 0.81 in females (

    vs.

    1.14 in males and 1.19 in females for

    O.kweichowensis

    ),lacks black bars on the lips (

    vs.

    the presence of black bars),an absence of vocal sacs in males (

    vs.

    the presence of vocal sacks),lacks large black spots on the dorsum (

    vs.

    large black spots in the center of the dorsum),has thighs with four brown bands and tibias with three (

    vs.

    thighs with five brown bands and tibias with six),lacks white pectoral spinules in mature males (

    vs.

    the presence of white pectoral spinules),and lacks toes with entire webbing (

    vs.

    full webbing).The new species differs from

    O.wuchuanensis

    by having a smaller body size (SVL 47.1–49.9 mm in adult males and 55.8–58.2 mm in adult females

    vs.

    71.1–76.5 mm in adult males and 75.8–90.0 mm in adult females) and a lower ratio of IND/IOD 0.89 in males and 0.81 in females (

    vs.

    1.63 in males and 1.14 in females in

    O.wuchuanensis

    ).The new species has a tibiotarsal articulation reaching to between the eye and the nostril when the leg is stretched forward (

    vs.

    reaching the nostril),lacks black bars on the lips (

    vs.

    the presence of black bars),lacks large black spots on the ventral surface (

    vs.

    the presence of large black spots),has relative finger lengths of Ivs.

    IIvs.

    having large black spots in the dorsum and ventral surfaces),and features males without white spines on the dorsal surface of the arm (

    vs.

    the presence of large white spines).

    Odorrana liboensis

    sp.nov.

    is phylogenetically closest to

    O.lipuensis

    ,and this new species could be distinguished from the latter by tibiotarsal articulation reaching to between the eye and the nostril when the leg is stretched forward (

    vs.

    reaching the anterior corner of the eye);relative finger lengths Ivs.

    I=IIvs.O.lipuensis

    ,mean 1.15);males and females with a higher ratio of IOD/HDW (in males 0.31–0.47,mean 0.40,and in females 0.39–0.41,mean 0.40

    vs.

    in males 0.24–0.31,mean 0.27,and in females 0.23–0.31,mean 0.27);males with a higher ratio of TED/HDL (0.11–0.13,mean 0.12

    vs.

    0.08–0.10,mean 0.09);the absence of small white spiny grains from the anterior corner of the eye along the eye and via the tympanic membrane below to the anterior part of the cloacal foramen (

    vs.

    the presence of grains);indistinct brownish-black bands on limbs (

    vs.

    distinct);the throat,chest,and abdomen lacking black spots (

    vs.

    the presence of black spots);and the ventral surface of the limbs is pink (

    vs.

    light purple).

    Distribution and ecology.

    Odorrana liboensis

    sp.nov.

    is known only from the type locality,Maolan National Nature Reserve,Libo County,Guizhou Province,China,at elevations between 645 and 728 m.The new species has only been found in one cave in the area located far from the village.There is no light in the cave,and the new species was found in a small pool approximately 1.5 m wide and 79 cm deep,about 145 m from the entrance of the cave,where the water temperature is approximately 20 °C all year round.The tadpoles were collected on July 23,2016,but no pairs of male and female adults were found to hold them,and with four years of survey data,the adults were collected only from mid-July to mid-August.Therefore,we speculate that the breeding period begins in late June and continues until about mid-August.Inside this cave,

    Odorrana liboensis

    sp.nov.

    is sympatric with

    Chinapotamon maolanense

    ,

    Sinocyclocheilus longibarbatus

    ,

    Hipposideros armiger,Aselliscus stoliczkanus

    ,

    Murina liboensis

    ,and

    Leopoldamys edwardsi

    .Outside the cave,no adults and tadpoles of the new species were found during a herpetological survey in the vicinity.However,a number of amphibians and reptiles can be found outside the cave,including

    Tylototriton asperrimus

    ,

    Quasipaa boulengeri

    ,

    Bufo gargarizans

    ,

    Kurixalus eif fingeri

    ,

    Rhacophorus dennysi

    ,

    Goniurosaurus liboensis

    ,

    Sinomicrurus

    macclellandi

    ,

    Lycodon flavozonatum

    ,and

    Orthriophis moellendorf fi.

    4.Discussion

    Most species of

    Odorrana

    live in montane streams.Previously,only two species (

    O.wuchuanensis

    and

    O.lipuensis

    ) had been reported to live in cave environments

    (Fei

    et al

    .,2012;Mo

    et al

    .,2015).We discovered,and report on,a third species of

    Odorrana

    ,

    Odorrana liboensis

    sp.nov.

    ,that inhabits cave environments.Two of these species (including

    O.wuchuanensis

    ) are found in karst caves in Guizhou.Phylogenetic analyses based on three mitochondrial genes suggested that

    Odorrana liboensis

    sp.nov.

    belongs to

    Odorrana

    but is distinct from its congeners.The genetic distances of 12S rRNA and 16S rRNA between the

    Odorrana liboensis

    sp.nov.

    and the closely related

    O.lipuensis

    were 6.06% and 5.19%.The genetic distances of 16S rRNA was greater than 5%,and this is greater than the distance that typically represents differentiation,at the species level,in frogs (>3%)(Vences

    et al

    .,2005;Fouquet

    et al

    .,2007).The new species was different from its congeners on the basis of many morphological characters,and this supports its validity.

    Odorrana liboensis

    sp.nov.

    described here increases the number of

    Odorrana

    species to 62,with 39 recorded from China (Frost,2021;AmphibiaChina,2021).Based on 49 previously named species,one new species,and three mitochondrial genes,we conducted a phylogenetic study of the genus

    Odorrana.

    The number of species covered here and the amount of data analyzed exceeds previously reported data(Chen

    et al

    .,2013;Li

    et al

    .,2018a).Eight highly supported major clades were identified in

    Odorrana

    (Clades A–H).The topology of this tree differed significantly from previous studies,focusing mainly on the root evolutionary branches of

    Odorrana

    (Figure 2).In the phylogenetic tree,

    O.lipuensis

    +

    Odorrana liboensis

    sp.nov.

    was the first species to diverge from the genus

    Odorrana

    ,whereas in other studies by Chen

    et al

    .(2013) and Li

    et al

    .(2018a),the first species to diverge was

    O.chapaensis

    .In our phylogenetic tree,

    O.chapaensis

    +O.

    geminata

    (Clade D1) serves as a sister taxonomic unit to Clade D2 consisting of

    O.andersonii

    as well as

    O.wuchuanensis

    ,

    O.mutschmanni

    ,

    O.yizhangensis

    ,

    O.lungshengensis

    ,and

    O.anlungensis

    .In contrast,in previous studies,

    O.chapaensis

    appeared as a sister taxon in a clade that included all other

    Odorrana

    (Chen

    et al

    .,2013;Li

    et al

    .,2018a),or as a sister taxon to

    Odorrana

    other than

    O.lipuensis

    (He,2017).Ye and Fei(2001) suggested that the primitive taxa of

    Odorrana

    may have originated from the Hengduan Mountains and the plateau of western Yunnan.The Guizhou plateau may have been the center of differentiation of

    Odorrana

    (Ye and Fei,2001).

    O.lipuensis

    and

    Odorrana liboensis

    sp.nov.

    appear to have diverged from

    Odorrana

    and formed the ancestral evolutionary branch of the genus

    Odorrana

    .Therefore,we believe that the ancestral distribution of

    Odorrana

    may be the region south of Guizhou and northwest of Guangxi.However,this preliminary speculation needs to be supported by data from additional species and nuclear genes.Biodiversity conservation in southwestern China is a priority of the Chinese government (Ministry of Environmental Protection,2015).Biodiversity conservation programs in this region play an important role in maintaining the stability of mountain ecosystems as well as protecting biodiversity.In the past three years alone,16 new amphibian species have been described from Guizhou Province,China (Zhang

    et al

    .,2017;Li

    et al

    .,2018a,b;Li

    et al

    .,2019a,b;Li

    et al

    .,2020a,b;Lyu

    et al

    .,2019;Wang

    et al

    .,2019;Luo

    et al

    .,2020;Liu

    et al

    .,2020;Lyu

    et al

    .,2020;Su

    et al

    .,2020;Wei

    et al

    .,2020;Wang

    et al

    .,2020;Cheng

    et al

    .,2021).The discovery of these new species suggests that amphibian species diversity in this region is severely underestimated.In the context of global warming,there is an urgent need for a comprehensive,systematic,and in-depth survey of the impacts of climate change on terrestrial vertebrates to provide a basis for scientific decisions regarding amphibian conservation (IPCC,2014).

    Acknowledgements This study was supported by programs of the Biodiversity Survey and Assessment Project of the Ministry of Ecology and Environment,China(2019HJ2096001006),the Strategic Priority Research Program B of the Chinese Academy of Sciences (CAS) (No.XDB31000000),the National Natural Science Foundation of China (NSFC 31460091),the National Animal Collection Resource Center,China (No.2005DKA21402),the Application of Amphibian Natural Antioxidant Peptides as Cosmetic Raw Material Antioxidants (No.QKZYD [2020]4002).We thank Gang LI,Xiang ZENG,Jun ZHOU,Hongtao CUI,Tao LIU,and Jian CHEN for their help during specimen collection.We thank Mr.Yunming MO,Ms.Xiaowen LIAO,and Prof.Weicai CHEN for their help in reviewing the specimens.We thank LetPub (www.letpub.com) for its linguistic assistance during manuscript preparation.

    Appendix

    Specimens examined

    Odorrana lipuensis

    (

    n

    =8):China:Guangxi:Lipu County (type locality):four males:NHMG 1306001,NHMG 1306002,NHMG 20140702,and NHMG 20 140703;four females:NHMG 1303019,NHMG 1303018,NHMG 1306003,and NHMG 20140701.

    Odorrana kweichowensis

    (

    n

    =35):China:Guizhou Province:Jinsha County (type locality).18 males:GZNU 20170717008–010,GZNU 20170717014,GZNU 20170717016–019,GZNU 20170717021–023,GZNU 20170717025,GZNU 20170717026,GZNU 20170717028–7030,GZNU 20170718003,and GZNU 2017 071800;17 females:GZNU 20170717001–007,GZNU 20170718001–002,GZNU 20170725001–007,and GZNU 20170725013.

    Odorrana

    wuchuanensis

    (

    n

    =17):China:Guizhou Province:Libo County (topotype locality):12 males:GZNU 20160729001–006,GZNU 20160729008,GZNU 20160806001–005;five females:GZNU 20160729007,GZNU 20160809001,GZNU 20160809003,GZNU 20160809006,and GZNU 20160809011.About the front cover:

    Bronchocela cristatella

    (Kuhl,1820) is a moderately slender arboreal lizard existing throughout Southeast Asia (Malaysia,Myanmar,Thailand,Singapore,India,Philippines,Indonesia,Papua New Guinea,Andaman and Nicobar Islands of India).Eggshell structure and composition of this lizard were reported by Kun GUO

    et al.

    in this issue.This lizard was drawn by Timan YANG and Hanwen KANG.

    猜你喜歡
    荔波
    荔波縣委常委會(huì)研究部署離退休干部黨建工作
    晚晴(2021年12期)2021-01-10 19:03:04
    美麗的荔波
    荔波國際兒童動(dòng)漫節(jié)
    荔波——?jiǎng)e樣的美
    梅花香自荔波來
    西南最大的冰雪主題樂園 荔波冰雪水世界樂園
    秘境荔波
    荔波縣委離退局干部職工參加全縣植樹活動(dòng)
    晚晴(2016年3期)2016-05-30 20:40:40
    荔波縣委離退局傳達(dá)學(xué)習(xí)州紀(jì)委文件精神
    晚晴(2016年5期)2016-05-14 08:01:10
    荔波三荔水庫庫區(qū)右岸巖溶滲漏分析
    十八禁网站免费在线| 午夜福利在线观看吧| 美女黄网站色视频| 国产亚洲精品综合一区在线观看| 国产单亲对白刺激| 日韩 亚洲 欧美在线| 1024手机看黄色片| 欧美xxxx黑人xx丫x性爽| 欧美一区二区亚洲| 五月玫瑰六月丁香| 日日撸夜夜添| 最新中文字幕久久久久| 国产真实伦视频高清在线观看| 亚洲激情五月婷婷啪啪| 一进一出抽搐gif免费好疼| 成人特级黄色片久久久久久久| 又黄又爽又免费观看的视频| 国产精品乱码一区二三区的特点| 不卡视频在线观看欧美| 真人做人爱边吃奶动态| 成人美女网站在线观看视频| 日本三级黄在线观看| av黄色大香蕉| 又黄又爽又刺激的免费视频.| 精品人妻一区二区三区麻豆 | 97超级碰碰碰精品色视频在线观看| 成人av在线播放网站| 亚洲四区av| 国产精品嫩草影院av在线观看| 国产亚洲av嫩草精品影院| 久久国内精品自在自线图片| 久久国内精品自在自线图片| 欧美色视频一区免费| 精品一区二区三区视频在线| 国产黄色视频一区二区在线观看 | 午夜福利高清视频| 12—13女人毛片做爰片一| 国产熟女欧美一区二区| 寂寞人妻少妇视频99o| 在线天堂最新版资源| 一区二区三区高清视频在线| 两个人的视频大全免费| 男女啪啪激烈高潮av片| 日本-黄色视频高清免费观看| 国产成人精品久久久久久| 亚洲成人久久性| 97碰自拍视频| 欧美又色又爽又黄视频| 国产 一区精品| 热99在线观看视频| 色综合亚洲欧美另类图片| 免费不卡的大黄色大毛片视频在线观看 | 日韩av在线大香蕉| 日韩欧美在线乱码| 国产aⅴ精品一区二区三区波| 精品一区二区三区av网在线观看| 久久久久久久久久久丰满| 久久久久精品国产欧美久久久| 国产高清三级在线| 国内精品久久久久精免费| 国产伦精品一区二区三区四那| 欧美高清成人免费视频www| 国产免费男女视频| 国产精品一区二区三区四区免费观看 | 男人舔奶头视频| 日韩,欧美,国产一区二区三区 | av在线观看视频网站免费| 精品一区二区三区视频在线观看免费| 亚洲自拍偷在线| 国产精品亚洲一级av第二区| 久久久久国产网址| 亚洲人成网站在线播放欧美日韩| 国产免费一级a男人的天堂| 亚洲精品日韩在线中文字幕 | 91av网一区二区| 久久久久久久久大av| 亚洲av.av天堂| 变态另类丝袜制服| 欧美性猛交╳xxx乱大交人| 亚洲av熟女| 99久久精品国产国产毛片| 日韩 亚洲 欧美在线| 精品日产1卡2卡| 亚洲美女搞黄在线观看 | 高清日韩中文字幕在线| 长腿黑丝高跟| 黄色视频,在线免费观看| 最近视频中文字幕2019在线8| av视频在线观看入口| 一边摸一边抽搐一进一小说| 亚洲精品国产av成人精品 | 特级一级黄色大片| 熟女电影av网| 大型黄色视频在线免费观看| 国产一区二区三区av在线 | 亚州av有码| 丰满乱子伦码专区| 国产高清有码在线观看视频| 18禁裸乳无遮挡免费网站照片| 别揉我奶头~嗯~啊~动态视频| 少妇熟女欧美另类| 蜜桃久久精品国产亚洲av| 午夜精品一区二区三区免费看| 在线国产一区二区在线| 亚洲专区国产一区二区| 欧美一区二区精品小视频在线| avwww免费| 日韩欧美一区二区三区在线观看| 九色成人免费人妻av| 免费观看人在逋| 国产亚洲精品久久久com| 久久草成人影院| 精品久久久噜噜| 欧美日韩国产亚洲二区| 国产白丝娇喘喷水9色精品| 亚洲中文字幕一区二区三区有码在线看| 嫩草影视91久久| 欧美日韩精品成人综合77777| 国产成年人精品一区二区| 超碰av人人做人人爽久久| 男人的好看免费观看在线视频| 日韩 亚洲 欧美在线| 久久久久久久久久久丰满| 久久久久久久久中文| 国产毛片a区久久久久| 欧美一区二区精品小视频在线| 黄片wwwwww| 长腿黑丝高跟| 欧美在线一区亚洲| 亚洲成人久久爱视频| 国产精品久久电影中文字幕| 看片在线看免费视频| 国产精品野战在线观看| 欧美激情久久久久久爽电影| 亚洲欧美精品自产自拍| eeuss影院久久| 亚洲无线在线观看| 日韩欧美精品免费久久| 欧美3d第一页| 欧美一级a爱片免费观看看| 内射极品少妇av片p| 国产爱豆传媒在线观看| 精品人妻视频免费看| 变态另类丝袜制服| 嫩草影院新地址| 一本一本综合久久| 极品教师在线视频| 看黄色毛片网站| 99久久久亚洲精品蜜臀av| eeuss影院久久| 91久久精品国产一区二区三区| 男插女下体视频免费在线播放| 成人特级av手机在线观看| 国产 一区精品| 亚州av有码| 黄片wwwwww| 国产精品人妻久久久影院| 色噜噜av男人的天堂激情| 精品乱码久久久久久99久播| 国产精品精品国产色婷婷| 春色校园在线视频观看| 91精品国产九色| 人人妻人人看人人澡| 国产精品亚洲一级av第二区| 18禁黄网站禁片免费观看直播| 亚洲国产欧美人成| 午夜免费激情av| 啦啦啦啦在线视频资源| 人妻丰满熟妇av一区二区三区| 亚洲经典国产精华液单| 亚洲av中文av极速乱| 高清日韩中文字幕在线| а√天堂www在线а√下载| 老司机影院成人| 国内久久婷婷六月综合欲色啪| 午夜日韩欧美国产| 国产一级毛片七仙女欲春2| 热99在线观看视频| 精品免费久久久久久久清纯| 国产精品一区www在线观看| 高清毛片免费看| 欧美激情国产日韩精品一区| 1000部很黄的大片| www.色视频.com| 成人漫画全彩无遮挡| 国产亚洲精品久久久com| 精品99又大又爽又粗少妇毛片| 国产精品一区二区免费欧美| 精品福利观看| 亚洲乱码一区二区免费版| 床上黄色一级片| 日本爱情动作片www.在线观看 | 欧美三级亚洲精品| 午夜福利在线在线| 亚洲av免费在线观看| 搡老妇女老女人老熟妇| 最近2019中文字幕mv第一页| 亚洲在线自拍视频| 小说图片视频综合网站| 黄片wwwwww| 一夜夜www| 成人鲁丝片一二三区免费| 日韩欧美国产在线观看| 国产午夜精品论理片| 欧美一区二区精品小视频在线| 亚洲国产欧洲综合997久久,| 午夜老司机福利剧场| 亚洲国产精品久久男人天堂| 国产成人freesex在线 | 国产视频内射| 日韩在线高清观看一区二区三区| 精品熟女少妇av免费看| 99久久精品热视频| 夜夜爽天天搞| 国产单亲对白刺激| 日日摸夜夜添夜夜添av毛片| 日韩一本色道免费dvd| 我的女老师完整版在线观看| 久久6这里有精品| 99热全是精品| 尾随美女入室| 久久久久久久久中文| 看免费成人av毛片| 亚洲成人精品中文字幕电影| 亚洲av成人av| 日本一二三区视频观看| 精品久久久久久久久久久久久| 国产精品国产高清国产av| 在线观看一区二区三区| 夜夜爽天天搞| 老师上课跳d突然被开到最大视频| 黄色日韩在线| 菩萨蛮人人尽说江南好唐韦庄 | 精品久久久噜噜| 熟女电影av网| 丰满人妻一区二区三区视频av| 免费搜索国产男女视频| 色尼玛亚洲综合影院| 久久鲁丝午夜福利片| av国产免费在线观看| 久久韩国三级中文字幕| 色av中文字幕| 一个人看的www免费观看视频| 久久99热这里只有精品18| 一区二区三区免费毛片| av黄色大香蕉| 黄色日韩在线| 久久中文看片网| 久久国内精品自在自线图片| 99热这里只有是精品50| 高清午夜精品一区二区三区 | 国产成年人精品一区二区| 午夜免费激情av| 免费观看精品视频网站| 亚洲激情五月婷婷啪啪| 综合色av麻豆| 日本撒尿小便嘘嘘汇集6| 久久久精品大字幕| av在线蜜桃| 天堂动漫精品| 亚洲av免费高清在线观看| 久久人妻av系列| 亚洲中文字幕日韩| 波多野结衣高清作品| 国产麻豆成人av免费视频| 免费一级毛片在线播放高清视频| 中文字幕av成人在线电影| 亚洲最大成人av| 俺也久久电影网| 中国国产av一级| 麻豆av噜噜一区二区三区| 91狼人影院| 欧美又色又爽又黄视频| 欧美一区二区亚洲| 欧美中文日本在线观看视频| 一个人看视频在线观看www免费| 国产一区二区三区av在线 | 成人午夜高清在线视频| 国产视频内射| 婷婷亚洲欧美| www日本黄色视频网| 校园人妻丝袜中文字幕| 中出人妻视频一区二区| 久久久久国产精品人妻aⅴ院| 国产伦精品一区二区三区视频9| 日韩av在线大香蕉| 亚洲人成网站在线播| 亚洲中文字幕日韩| 国产一区二区激情短视频| 国产亚洲欧美98| 成年版毛片免费区| 久久这里只有精品中国| 一本久久中文字幕| 亚洲人成网站高清观看| 国内少妇人妻偷人精品xxx网站| 久久久久久国产a免费观看| 麻豆乱淫一区二区| 久久久久久久久久成人| 91av网一区二区| 日韩高清综合在线| 精品久久久久久久久久免费视频| av在线蜜桃| 九九热线精品视视频播放| 51国产日韩欧美| 国产伦一二天堂av在线观看| avwww免费| 给我免费播放毛片高清在线观看| 99热这里只有是精品50| 国内少妇人妻偷人精品xxx网站| 久久精品综合一区二区三区| 日本黄色片子视频| 黄色日韩在线| 久久人人爽人人爽人人片va| 韩国av在线不卡| 国语自产精品视频在线第100页| 久久中文看片网| 亚洲人与动物交配视频| 国产精品亚洲美女久久久| 色噜噜av男人的天堂激情| 亚洲精品久久国产高清桃花| 久久久久久久亚洲中文字幕| 国产精品久久久久久亚洲av鲁大| 综合色丁香网| 干丝袜人妻中文字幕| 人人妻人人澡人人爽人人夜夜 | 日本黄大片高清| 黑人高潮一二区| 国产单亲对白刺激| 在线天堂最新版资源| 国产一区二区激情短视频| 夜夜夜夜夜久久久久| 国产老妇女一区| 精品久久国产蜜桃| 麻豆精品久久久久久蜜桃| 日本成人三级电影网站| av在线播放精品| 91精品国产九色| 全区人妻精品视频| 成人漫画全彩无遮挡| 别揉我奶头 嗯啊视频| 午夜福利在线观看免费完整高清在 | 成人亚洲欧美一区二区av| 亚洲自偷自拍三级| 久久人人爽人人片av| 看免费成人av毛片| 欧美精品国产亚洲| 如何舔出高潮| 精品一区二区三区视频在线观看免费| 激情 狠狠 欧美| 精品午夜福利视频在线观看一区| av.在线天堂| 激情 狠狠 欧美| 亚洲图色成人| ponron亚洲| 岛国在线免费视频观看| 观看免费一级毛片| 欧美+亚洲+日韩+国产| 69av精品久久久久久| 亚洲国产精品久久男人天堂| av在线蜜桃| 天美传媒精品一区二区| 久久久a久久爽久久v久久| 国产av在哪里看| videossex国产| 99精品在免费线老司机午夜| 此物有八面人人有两片| 美女xxoo啪啪120秒动态图| 一级av片app| 免费观看在线日韩| 综合色丁香网| 乱码一卡2卡4卡精品| ponron亚洲| 亚洲婷婷狠狠爱综合网| 久久久欧美国产精品| 草草在线视频免费看| 午夜激情福利司机影院| 日韩精品中文字幕看吧| 长腿黑丝高跟| 欧美极品一区二区三区四区| 中文资源天堂在线| 91麻豆精品激情在线观看国产| 久久人人爽人人片av| 亚洲精品一卡2卡三卡4卡5卡| 国产高清视频在线观看网站| 欧美丝袜亚洲另类| 欧美日韩国产亚洲二区| 日日撸夜夜添| 国产精品一二三区在线看| 免费观看在线日韩| 欧美性猛交╳xxx乱大交人| 淫妇啪啪啪对白视频| 日日啪夜夜撸| 最近手机中文字幕大全| 国产在线男女| 色综合亚洲欧美另类图片| 免费在线观看成人毛片| 国产探花极品一区二区| 有码 亚洲区| 乱码一卡2卡4卡精品| 在线观看一区二区三区| 久久精品夜色国产| 你懂的网址亚洲精品在线观看 | 欧美区成人在线视频| 亚洲精品一区av在线观看| 夜夜爽天天搞| 精品久久久久久久人妻蜜臀av| 久久精品夜夜夜夜夜久久蜜豆| 国产片特级美女逼逼视频| 欧美激情久久久久久爽电影| 成年女人永久免费观看视频| 99久久成人亚洲精品观看| 亚洲乱码一区二区免费版| 在线免费十八禁| 久久精品国产亚洲av涩爱 | a级毛片a级免费在线| 最近最新中文字幕大全电影3| 美女大奶头视频| 精品一区二区三区视频在线| 免费av不卡在线播放| 日本免费a在线| 亚洲av中文字字幕乱码综合| 国产私拍福利视频在线观看| 在线免费观看不下载黄p国产| av天堂中文字幕网| 国产精品女同一区二区软件| 淫妇啪啪啪对白视频| 九九久久精品国产亚洲av麻豆| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 免费搜索国产男女视频| 欧美日本视频| 黑人高潮一二区| 亚洲不卡免费看| 日产精品乱码卡一卡2卡三| 久久久色成人| 国产精品1区2区在线观看.| 免费看美女性在线毛片视频| 日韩大尺度精品在线看网址| 国产伦精品一区二区三区视频9| 成年女人毛片免费观看观看9| 国产爱豆传媒在线观看| 最近手机中文字幕大全| 欧美不卡视频在线免费观看| or卡值多少钱| 欧美潮喷喷水| 久久精品国产清高在天天线| aaaaa片日本免费| 1000部很黄的大片| 亚洲最大成人手机在线| 亚洲,欧美,日韩| 国国产精品蜜臀av免费| 久久久精品大字幕| 99久久中文字幕三级久久日本| 欧美高清性xxxxhd video| 村上凉子中文字幕在线| 最近2019中文字幕mv第一页| 欧美日韩一区二区视频在线观看视频在线 | 精品久久久久久成人av| 2021天堂中文幕一二区在线观| 麻豆国产av国片精品| 亚洲电影在线观看av| 嫩草影院入口| 看十八女毛片水多多多| 最近中文字幕高清免费大全6| 尾随美女入室| 嫩草影院精品99| 美女被艹到高潮喷水动态| 国产高清视频在线观看网站| 天堂影院成人在线观看| 人人妻人人澡欧美一区二区| 国产精品一区www在线观看| 亚洲色图av天堂| 岛国在线免费视频观看| 熟妇人妻久久中文字幕3abv| 俄罗斯特黄特色一大片| 午夜免费男女啪啪视频观看 | 亚洲欧美日韩高清在线视频| 久久中文看片网| 真实男女啪啪啪动态图| 丰满人妻一区二区三区视频av| 国产精品久久久久久亚洲av鲁大| 91麻豆精品激情在线观看国产| 亚洲久久久久久中文字幕| 日日啪夜夜撸| 能在线免费观看的黄片| 国产极品精品免费视频能看的| 97热精品久久久久久| 久久鲁丝午夜福利片| 偷拍熟女少妇极品色| 18+在线观看网站| 久久久精品94久久精品| 欧美最黄视频在线播放免费| 精品乱码久久久久久99久播| 夜夜夜夜夜久久久久| 免费看a级黄色片| 最新在线观看一区二区三区| 久久久午夜欧美精品| 久久久久性生活片| 12—13女人毛片做爰片一| 18禁裸乳无遮挡免费网站照片| 亚洲精品一卡2卡三卡4卡5卡| 美女内射精品一级片tv| 看十八女毛片水多多多| 日韩欧美免费精品| 久久久欧美国产精品| 国产精品日韩av在线免费观看| av在线老鸭窝| 在线观看一区二区三区| 色av中文字幕| 在线观看一区二区三区| 亚洲av电影不卡..在线观看| 黄片wwwwww| 内地一区二区视频在线| 黄片wwwwww| 欧美成人免费av一区二区三区| 日产精品乱码卡一卡2卡三| 亚洲人与动物交配视频| 男女啪啪激烈高潮av片| 狂野欧美白嫩少妇大欣赏| 国内少妇人妻偷人精品xxx网站| 乱系列少妇在线播放| 两个人的视频大全免费| 亚洲av免费高清在线观看| 日本成人三级电影网站| 成人无遮挡网站| 成人亚洲欧美一区二区av| 亚洲高清免费不卡视频| av福利片在线观看| 18禁在线无遮挡免费观看视频 | 免费观看人在逋| 日本熟妇午夜| 亚洲av一区综合| 国产成人一区二区在线| 噜噜噜噜噜久久久久久91| 91av网一区二区| aaaaa片日本免费| 久久热精品热| 精品一区二区三区视频在线| av在线蜜桃| 成人性生交大片免费视频hd| 精品人妻偷拍中文字幕| 午夜激情福利司机影院| 国产高潮美女av| 国产精品美女特级片免费视频播放器| 男女啪啪激烈高潮av片| 久久人人爽人人片av| 午夜福利18| 极品教师在线视频| 黄片wwwwww| 尤物成人国产欧美一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 免费看a级黄色片| 最新在线观看一区二区三区| 最好的美女福利视频网| 一级av片app| 黄片wwwwww| 久久精品国产鲁丝片午夜精品| 国产蜜桃级精品一区二区三区| 一本精品99久久精品77| 国产精品一及| 亚洲在线观看片| 男女边吃奶边做爰视频| 午夜日韩欧美国产| 日韩欧美免费精品| 亚洲国产精品合色在线| 国模一区二区三区四区视频| 亚洲一区二区三区色噜噜| 人妻少妇偷人精品九色| 嫩草影视91久久| 麻豆av噜噜一区二区三区| 欧美日韩精品成人综合77777| 欧美另类亚洲清纯唯美| 国产精品久久久久久精品电影| 日本在线视频免费播放| 久久中文看片网| 免费观看精品视频网站| 亚洲精品国产成人久久av| 欧美日本视频| 亚洲av熟女| 在线a可以看的网站| 国内精品宾馆在线| 男女啪啪激烈高潮av片| 欧美3d第一页| 久久这里只有精品中国| 国产久久久一区二区三区| 在线观看一区二区三区| 国产乱人偷精品视频| 国产精品,欧美在线| 啦啦啦观看免费观看视频高清| 亚洲精品国产成人久久av| 国产成人freesex在线 | 亚洲经典国产精华液单| 给我免费播放毛片高清在线观看| 天天躁夜夜躁狠狠久久av| 美女高潮的动态| 久久亚洲精品不卡| 国产在线精品亚洲第一网站| 国产精品久久久久久精品电影| 欧美一级a爱片免费观看看| 亚洲无线在线观看| 亚洲国产精品合色在线| 高清毛片免费观看视频网站| 小说图片视频综合网站| 午夜爱爱视频在线播放| 日本与韩国留学比较| 老司机福利观看| 高清日韩中文字幕在线| 国产av麻豆久久久久久久| 日本免费a在线| 亚洲人成网站在线播放欧美日韩| 人妻少妇偷人精品九色| 91av网一区二区| 日韩欧美国产在线观看| 国产av不卡久久| 91在线观看av|