• <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 the genus Theloderma Tschudi,1838(Amphibia:Anura:Rhacophoridae)from Tay Nguyen Plateau,central Vietnam

    2018-06-28 01:56:00NikolayPoyarkovJrIvanKropachevSvetlanaGogolevaNikolaiOrlov
    Zoological Research 2018年3期

    Nikolay A.Poyarkov,Jr.,Ivan I.Kropachev,Svetlana S.Gogoleva,Nikolai L.Orlov

    1Department of Vertebrate Zoology,Biological Faculty,Lomonosov Moscow State University,Moscow 119234,Russia

    2Joint Russian-Vietnamese Tropical Research and Technological Center,Nghia Do,Cau Giay,Hanoi,Vietnam

    3Tula Exotarium,Tula 300004,Russia

    4A.N.Severtsov Institute of Ecology and Evolution,Russian Academy of Sciences,Moscow 119071,Russia

    5Joint Russian-Vietnamese Tropical Research and Technological Center,South Branch,Ho Chi Minh City,Vietnam

    6Zoological Museum of the Lomonosov Moscow State University,Moscow 125009,Russia

    7Zoological Institute,Russian Academy of Sciences,St.Petersburg 199034,Russia

    INTRODUCTION

    Rhacophoridae is a diverse family of largely arboreal tree frogs consisting of 414 currently recognized species in 18 genera(AmphibiaWeb,2018;Che&Wang,2016;Frost,2018).Rhacophorids are distributed throughout Sub-Saharan Africa,China,Southern and Southeast Asia,Japan,Taiwan,the Philippines,and the Greater Sunda Islands(Frost,2018). Species diagnostics and generic allocation within Rhacophoridae can be difficult,and molecular phylogenetic analyses are essential for correct taxonomic decisions(Jiang et al.,2016;Li et al.,2008,2009;Orlov et al.,2012;Poyarkov et al.,2015;Rowley et al.,2011;Yu et al.,2009).

    Frogs of the genusThelodermaTschudi,1838 are highly arboreal small-to large-sized rhacophorid frogs,which were traditionally grouped in one genus based on the presence of dorsal asperities or tuberculate dorsal skin(Liem,1970;Rowley et al.,2011).However,recent studies have reported on several smooth species of rhacophorids that almost completely or totally lack dorsal asperities but are still assigned toThelodermabased on molecular data and reproductive biology(Nguyen et al.,2014;Orlov et al.,2012;Poyarkov et al.,2015).Thus,the taxonomy ofThelodermais in a current state of flux.Several studies have indicated a sister-clade relationship betweenNyctixalusBoulenger,1882 andTheloderma,with monophyly of the latter being poorly supported or not recovered(Dever et al.,2015;Dever,2017;Rowley et al.,2011).Thus,Poyarkov et al. (2015)recognizedNyctixalusas a subgenus ofThelodermaand assigned theTh.stellatumTaylor,1962-Th.horridum(Boulenger,1903)species group to the subgenusStelladermaPoyarkov,Orlov,Moiseeva,Pawangkhanant,Ruangsuwan,Vassilieva,Galoyan,Nguyen&Gogoleva,2015. Subsequent studies with larger taxon and gene sampling(Nguyen et al.,2015;Sivongxay et al.,2016)supportedNyctixalusas a sister clade with respect toThelodermaas well as monophyly ofThelodermaand its subdivision into two highly divergent clades(corresponding to subgeneraThelodermaandStelladerma sensuPoyarkov et al.,2015).Theloderma moloch(Annandale,1912),an enigmatic taxon whose phylogenetic position and assignment toThelodermabecame uncertain due to sequences published by Li et al.(2009),was simultaneously rediscovered in Chinese and Indian regions of the Himalayas and shown to be a member ofThelodermas. str. (Hou et al.,2017;Lalronunga&Lalrinchhana,2017;Li et al.,2016). In addition,Jiang et al.(2016)and Biju et al.(2016)concurrently demonstrated that sequences previously reported as “Th.moloch”by Li et al.(2009)actually corresponded to a distinct lineage of Rhacophoridae,representing new genusNasutixalusJiang,Yan,Wang&Che,2016(Jiang et al.,2016).The taxonomic position ofTh.andersoni(Ahl,1927)is also relatively unclear,with Hou et al.(2017)suggesting that this species might be a member of the genusRaorchestesBiju,Shouche,Dubois,Dutta&Bossuyt,2010.Taking these reassignments and new species descriptions into consideration,a total of 25 species ofThelodermaare recognized and distributed throughout Southeast Asia,from Assam in northeastern India to Myanmar,southern China and Indochina,the Malay Peninsula,and Sumatra and Borneo of the Greater Sunda Islands. New species in the genus continue to be discovered,with 13 described in the last 15 years alone(Frost,2018).

    During recent fieldwork in the montane forests of Tay Nguyen Plateau in central Vietnam,we encountered a small-sized rhacophorid species with a smooth dorsum.This tree frog was identified morphologically as a probable member of the genusThelodermabased on the presence of an intercalary cartilage between the terminal and penultimate phalanges of digits,a distinct tympanum,tips of the digits expanded into large disks bearing circummarginal grooves,absence of vomerine teeth,horizontal pupils,rounded canthus rostralis,terminal phalanx with a Y-shaped distal end,and the skin of the head not co-ossified to the skull(Liem,1970;McLeod&Ahmad,2007;Rowley et al.,2011).This species of frog was earlier reported by Orlov&Ananjeva(2007)and Orlov et al.(2012)asTheloderma laeve(Smith,1924),a taxon described from the Langbian Plateau of southern Vietnam and originally assigned to the genusPhilautusdue to the presence of smooth skin on the dorsum and the lack of warts,tubercles,and asperities(Smith,1924).Following examination of the holotype ofPhilautus laevisSmith,1924,Poyarkov et al.(2015)demonstrated it to be conspecific toTheloderma bambusicolumOrlov,Poyarkov,Vassilieva,Ananjeva,Nguyen,Nguyen&Geissler,2012 and confirmed its assignment to the genusThelodermaasTheloderma laeve(Smith,1924).Consequently,the taxon previously identified asTh.laeveby Orlov et al.(2012)is not assigned to any currently recognized species ofTheloderma.Close examination showed that this taxon can be easily distinguished from other known members of the genusThelodermaby a combination of several adult morphological characteristics.Molecular phylogenetic analysis of a 2 518-bp mtDNA fragment further indicated that this taxon is nested within the genusThelodermawith high support values and represents a sister species toTh.palliatumRowley,Le,Hoang,Dau&Cao,2011.In the present paper we discuss the phylogenetic position and taxonomic affiliation of this frog and describe it as a new species.

    MATERIALS AND METHODS

    Sample collection

    Fieldwork was carried out in central Vietnam in May 2016 and May 2017.Specimens ofThelodermasp.were collected by hand in the montane evergreen tropical forests of Kon Ka Kinh National Park and Kon Chu Rang Nature Reserve,both in Gia Lai Province;and in Thac Nham Forest and Kon Plinh Forest,Mang Canh District,Kon Tum Province;Tay Nguyen Plateau,central Vietnam(Figure 1).Geographic coordinates and altitude were obtained using a Garmin GPSMAP 60CSx GPS receiver(USA)and recorded in datum WGS 84.Specimens were euthanized by 20%benzocaine and tissue samples(femoral muscles)for genetic analysis were collected and stored in 96%ethanol prior to specimen preservation.Specimens were subsequently preserved in 70%ethanol and deposited in the herpetological collection of the Zoological Museum of Moscow State University(ZMMU)in Moscow,Russia.Comparative materials examined are stored in the herpetological collections of ZMMU and the Zoological Institute R.A.S.(ZISP)in St.Petersburg,Russia.A list of examined specimens is given in Appendix I.

    Morphological description

    Specimens ofThelodermasp.were photographed in life and after preservation.Sex was determined by direct observation of calls and by examination of the nuptial pads in collected males.Measurements were taken using a digital caliper under a light dissecting microscope to the nearest 0.01 mm,subsequently rounded to 0.1 mm.Statistical analyses were performed with Statistica 8.0(StatSoft,Inc.,2007).

    Figure 1 Known distribution of Theloderma auratum sp.nov.and Th.palliatum in Vietnam

    Adult morphology

    Morphometrics follow Orlov et al.(2012),Milto et al.(2013),and Poyarkov et al. (2015)and include the following 30 measurements:(1)SVL(snout-vent length);(2)A-G(axilla to groin,distance from posterior base of forelimb at its emergence from body to anterior base of hind limb at its emergence from body);(3)HW(head width at greatest cranial width);(4)HL(head length from rear of lower jaw to tip of the snout);(5)HD(head depth,greatest transverse depth of head,taken beyond interorbital region);(6)UEW(upper eyelid width,greatest width of upper eyelids);(7)IOD(interorbital distance);(8)ED(horizontal diameter of eye);(9)TD(horizontal diameter of tympanum);(10)ESL(tip of snout-eye distance);(11)IND(internarial distance between nostrils);(12)END(eye to nostril distance from anterior corner of eye to nostril);(13)TED(tympanum-eye distance from anterior edge of tympanum to posterior corner of eye);(14)NS(distance from nostril to tip of snout);(15)FLL(length of forelimb from tip of disk of finger III to axilla);(16)LAL(forearm length,from elbow to base of outer palmar tubercle);(17)ML(hand length from tip of third digit to base of outer palmar tubercle);(18)FFL( first finger length);(19)TFL(third finger length);(20)FTD(maximal diameter of disk of finger III);(21)NPL(nuptial pad length,measured for males only);(22)MKTe(length of external metacarpal tubercle);(23)HLL(length of hindlimb from tip of disk of toe IV to groin);(24)FL(femur length);(25)TL(tibia length);(26)FOT(length of hindlimb from tip of disk of toe IV to posterior edge of tibia);(27)FTL( first toe length);(28)FFTL(fourth toe length);(29)HTD(diameter of fourth toe tip,greatest diameter of disk on fourth toe);(30)MTTi(length of internal metatarsal tubercle).Skin texture,dorsal coloration,ventral coloration,and presence of supratympanic folds,circummarginal grooves,dorsolateral folds,vomerine teeth,hind limb and forelimb webbing,and dorsal,lateral,and ventral coloration were recorded.Webbing formula follows Savage&Heyer(1997).

    Data and adult morphological characters chosen for comparison with otherThelodermaspecies and allied taxa were taken from the following sources:Ahl(1927,1931);Anderson(1879);Annandale(1912);Bain et al. (2009);Boulenger(1882,1886,1893,1903);Bourret(1937,1942);Chan&Ahmad(2009);Chan-ard(2003);Das et al.(2013);Das(2007);Dever et al.(2015);Dever(2017);Fei et al.(1990,2009,2010);Hou et al.(2017);Inger(1954,1966);Lalronunga&Lalrinchhana(2017);Li et al.(2011,2016);Liem(1970);Liu&Hu(1962);Manthey&Grossmann(1997);Mathew&Sen.(2010);McLeod&Ahmad(2007);Nguyen et al.(2014);Orlov&Ananjeva(2007);Orlov&Ho(2005);Orlov et al.(2006,2010,2012);Orlov(2007);Peters(1871);Poyarkov et al.(2015);Rowley et al.(2011);Sivongxay et al.(2016);Smith(1924,1926,1930,1931);Stuart&Heatwole(2004);Taylor(1920,1962);Von Tschudi(1838);Wolf(1936);Zhao&Adler(1993).

    Larval morphology

    Morphometric data and abbreviations for larval characters follow Altig&McDiarmid(1999),Altig(2007),Grosjean(2001),and Hendrix et al.(2007).Labial tooth row formulae(LTRF)follow Altig&McDiarmid(1999). Tadpoles were staged according to the tables of Gosner(1960).All measurements were taken on freshly preserved specimens.Tadpoles were euthanized with ethyl acetate and preserved in 70%ethanol.The following 17 measurements were taken:(1)BH(maximum body height);(2)BL(body length);(3)BW(maximum body width);(4)ED(maximum diameter of eye);(5)NP(naro-pupilar distance);(6)IND(internarial distance,measured between centers of narial apertures);(7)IOD(interorbital distance);(8)ODW(oral disk width);(9)RN(rostro-narial distance);(10)SS(distance from tip of snout to spiracle opening);(11)LF(maximum height of lower tail fin);(12)MTH(maximum tail height);(13)TAL(tail length);(14)TL(total length);(15)TMH(height of tail musculature at base);(16)TMW(width of tail musculature at base);(17)UF(maximum height of upper tail fin).

    DNA isolation,PCR,and sequencing

    Total genomic DNA was extracted from the ethanol-preserved muscle tissues using standard phenol-chloroform extraction(Hillis et al.,1996).Total DNA concentration was estimated in 1μL using a NanoDrop 2000(Thermo Scientific,USA),and consequently adjusted to 100 ng DNA/μL.

    We amplified the mtDNA fragment consisting of partial sequences of 12S rRNA,tRNAval,and 16S rRNA mtDNA genes.These markers were chosen due to their successful use in studies of Rhacophorid taxonomic diversity(Bain et al.,2009;Dever et al.,2015;Li et al.,2008,2009,2013;Meegaskumbura et al.,2015;Nguyen et al.,2014,2015;Poyarkov et al.,2015;Rowley et al.,2011;Wilkinson&Drewes,2000;Wilkinson et al.,2002;Yu et al.,2007,2008;and references therein).PCR was performed in 20-μL reactions using 50 ng of genomic DNA,10 nmol of each primer,15 nmol of each dNTP,50 nmol additional MgCl2,Taq PCR buffer(10 mmol/L Tris-HCl,pH 8.3,50 mmol/L KCl,1.1 mmol/L MgCl2,and 0.01%gelatin),and 1 U of Taq DNA polymerase.Primers used in PCR and sequencing were obtained from previous studies(Hedges,1994;Li et al.,2008,2009;Wilkinson et al.,2002)and are summarized in Table 1.The PCR conditions for amplification of the 12S rRNA and 16S rRNA gene fragments followed those described in Poyarkov et al.(2015).

    PCR products were visualized by agarose electrophoresis in the presence of ethidium bromide and consequently purified using 2μL from a 1:4 dilution of ExoSapIt(Amersham,UK)per 5μL of PCR product prior to cycle sequencing.Sequence data collection and visualization were performed on an ABI 3730xl automated sequencer(Applied Biosystems,USA)in Evrogen Inc.,Moscow.The obtained sequences were aligned and deposited in GenBank under the accession numbers MG917762-MG917772(Table 2).

    Phylogenetic analyses

    Sequences of the 12S rRNA,tRNAval,and 16S rRNA mtDNA fragments from 135 Rhacophoridae specimens,including 122 representatives ofTheloderma(ca. 27 species),and 13 sequences of outgroup members of Rhacophoridae(generaRhacophorus,Nasutixalus,andNyctixalus)were included in the fi nal alignment with a total length of 2 518 bp.Details on voucher specimens and GenBank accession Nos.used in phylogenetic analyses are summarized in Table 2.Nucleotide sequences were initially aligned using ClustalX 1.81(Thompson et al.,1997)with default parameters,and then checked by eye in BioEdit 7.0.5.2(Hall,1999)and MEGA 7.0(Kumar et al.,2016)and slightly adjusted.

    Table 1 Primers used in this study

    Table 2 Sequences and voucher specimens of Theloderma and outgroup taxa used in this study

    Continued

    Continued

    Continued

    The dataset was divided into three partitions,12S rRNA,tRNAval,and 16S rRNA,with the optimal evolutionary models then estimated using MODELTEST v.3.06(Posada&Crandall,1998).For the 12S and 16S rRNA partitions,the best- fi tting model according to the Akaike information criterion(AIC)was the HKY+I+G model;whereas for the tRNAvalpartition,the Kimura 2-parameter model(+G+I)was selected as the one of best fi t.Mean uncorrected genetic distances(P-distances)between sequences were determined with MEGA 7.0(Kumar et al.,2016).

    Matrilineal genealogy was inferred using Bayesian inference(BI)and maximum likelihood(ML)algorithms.The BI analyses were conducted in MrBayes 3.1.2(Huelsenbeck&Ronquist,2001;Ronquist&Huelsenbeck,2003).Metropolis-coupled Markov chain Monte Carlo(MCMCMC)analyses were run with one cold chain and three heated chains for ten million generations and sampled every 1 000 generations. Five independent MCMCMC runs were performed and 1 000 trees were discarded as burn-in.Confidence in tree topology was tested by posterior probability(PP)for Bayesian inference(BI)trees(Huelsenbeck&Ronquist,2001). Nodes with posterior probability values over 0.95 werea prioriregarded as sufficiently resolved,whereas those between 0.95 and 0.90 were regarded as tendencies,and those below 0.90 were considered as unsupported.

    The ML analyses were conducted using Tree finder(Jobb et al.,2004)and confidence in node topology was tested by non-parametric bootstrapping with 1 000 replicates(ML BS,see Felsenstein,1985).Wea prioriregarded tree nodes with bootstrap(ML BS)values of 70%or greater and Bayesian posterior probabilities(BI PP)values over 0.95 as sufficiently resolved;ML BS values between 70%and 50%(BI PP between 0.95 and 0.90)were treated as tendencies and nodes with ML BS values below 50%(BI PP below 0.90)were regarded as unresolved(Felsenstein,2004;Huelsenbeck&Hillis,1993).

    Acoustic analyses

    Advertisement calls of theThelodermasp. were recorded in habitat in Kon Chu Rang Nature Reserve,Gia Lai Province,Tay Nguyen Plateau,Vietnam (N14°30′19.7′′,E108°32′30.0′′;950 m a.s.l.) on 26 May 2016 from 0224 h to 0300 h at an ambient temperature of 21.5°C(temperature was measured at the calling site immediately after recording with a digital thermometer KTJ TA218A Digital LCD Thermometer-Hydrometer)using a portable digital audio recorder Zoom h5(ZOOM Corporation,Tokyo,Japan)in stereo mode with 48 kHz sampling frequency and 16-bit precision.In total,two recordings from two males(ZMMU A-5828 and A-5829)were made.

    Calls were analyzed using Avisoft SASLab Pro software v.5.2.05(Avisoft Bioacoustics,Germany). Before analysis,we reduced background noise using low-pass(up to 500 Hz)and high-pass filters(down to 7 kHz).All parameters were measured in the spectrogram window of Avisoft.Spectrograms were created under Hamming window,FFT-length 1 024 points,frame 75%,and overlap 87.5%.For figure spectrograms,we lowered the sampling rate to 22.05 kHz.Figure spectrograms were created under Hamming window,FFT-length 512 points,frame 100%,and overlap 50%.In total,we measured 214Thelodermasp.calls.

    We measured four temporal parameters(duration of each call and series and interval between successive calls and series)and three frequency parameters(initial and final fundamental frequency and frequency of maximum amplitude).We also calculated the call repetition rate(calls/s)by counting the number of calls within each series,minus one,and dividing that number by the series duration.All numeral parameters are given as means±SE,and the minimum and maximum values are given in parentheses(min-max).

    RESULTS

    Phylogenetic analyses Sequences and statistics

    The final alignment contained 2 518 aligned characters,with 1 090 conserved sites and 1 395 variable sites,of which 1 170 were found to be parsimony-informative.The transition-transversion bias(R)was 2.372(data given for ingroup only). Nucleotide frequencies were 30.10%(A),20.18%(T/U),27.40%(C),and 22.31%(G).

    Position of Theloderma sp.in matrilineal genealogy

    The phylogenetic analyses results are presented in Figure 2.The ML and BI phylogenetic analyses resulted in essentially similar topologies.In general,the BI cladogram topology was consistent with results reported in previous work(Nguyen et al.,2015;Sivongxay et al.,2016),suggesting monophyly of the clade joiningNyctixalusandTheloderma(node support values 1.0/100,hereafter given for BI PP/ML BS,respectively)and monophyly ofTheloderma,though only with moderate support(0.90/85).

    Figure 2 Bayesian inference dendrogram of Theloderma and its relatives derived from analysis of 2 518-bp length 12SrRNA-16S rRNA mtDNA gene fragments

    The genusThelodermawas divided into two strongly-supported major clades:Clade I(1.0/100),joiningTh.stellatum,Th.Horridum,andTh.vietnamensePoyarkov,Orlov,Moiseeva,Pawangkhanant,Ruangsuwan,Vassilieva,Galoyan,Nguyen&Gogoleva,2015(corresponding to the subgenusStelladerma)and Clade II(1.0/90),joining all remaining species of the genus(subgenusThelodermas.str.)(Figure 2).Clade II was further subdivided into subclade II-A(Sundaland,Indochina,eastern Himalayas,southern China;1.0/96)and subclade II-B(Indochinese species;1.0/85).

    Subclade II-A joined two species groups:theTh.molochgroup consisting of large-sized warty species from Indochina and the eastern Himalayas(i.e.,Th.moloch,Th.phrynoderma(Ahl,1927),andTh.ryaboviOrlov,Dutta,Ghate&Kent,2006;1.0/100))and theTh.asperumgroup consisting of medium to small-sized species with contrasting black and white dorsal coloration(i.e.,Th.asperum(Boulenger,1886),Th.pyaukkyaDever,2017,Th.baibungense(Jiang,Fei&Huang,2009),Th.licinMcLeod&Ahmad,2007,Th.petilum(Stuart&Heatwole,2004),andTh.albopunctatum(Liu&Hu,1962);1.0/100)).Our data strongly suggested paraphyly ofTh.pyaukkyawith respect toTh.baibungense,with the former subdivided into two divergent lineages from northern and central Myanmar(A and B,see Figure 2)with poorly resolved genealogical relationships.Furthermore,Th.albopunctatumrepresents a species complex with at least four divergent lineages with unclear geographic structuring(see Figure 2,Table 2).

    Phylogenetic relationships within subclade II-B were essentially unresolved (Figure 2),though the following species groups were supported:Th.leporosumgroup(joiningTh.leporosumTschudi,1838 andTh.gordoniTaylor,1962;1.0/100),Th.lateriticumgroup(includingTh.lateriticumBain,Nguyen&Doan,2009 andTh.lacustrinumSivongxay,Davankham,Phimmachak,Phoumixay&Stuart,2016;1.0/100),Th.laevegroup(includingTh.laeve,Th.truongsonense(Orlov&Ho,2005),Th.nebulosumRowley,Le,Hoang,Dau&Cao,2011,andTh.annaeNguyen,Pham,Nguyen,Ngo&Ziegler,2016;1.0/100),andTh.corticalegroup(joiningTh.corticale(Boulenger,1903),Th.bicolor(Bourret,1937),Th.rhododiscus(Liu&Hu,1962),Th.palliatum,andThelodermasp.from the Tay Nguyen Plateau;0.99/94).A single specimen(VNMN 4403)from Gia Lai Province identified asTh.laeveby Nguyen et al.(2015)was grouped with samples ofTh.laeves.str.from southern Vietnam with high support(1.0/95)but represented a clearly distinct genealogical lineage and was identified asThelodermasp.in the present paper.

    TheThelodermasp. from Tay Nguyen Plateau was recovered as a sister species ofTh.palliatum,which inhabits the montane forests of Langbian Plateau,with strong node support(1.0/100)(Figure 2);however,no clear structuring was observed within theThelodermasp.clade.

    Sequence divergence

    The uncorrectedP-distances among and within the studied 2 518-bp mtDNA fragments for the examinedThelodermaspecies are shown in Table 3(data given for the ingroup only).The interspecific uncorrected geneticP-distances between the

    Thelodermasp.from Tay Nguyen Plateau and other congeners varied from 8.9%(betweenThelodermasp.and sister speciesTh.palliatum)to 17.9%(betweenThelodermasp.andTh.laeve)(Table 3). This degree of pairwise divergence was high,notably greater than the genetic divergence thresholds representing species level differentiation in frogs(Vences et al.,2005a,2005b;Vieites et al.,2009).

    Taxonomy

    Based on the phylogenetic analyses of the 2 518-bp length 12S rRNA and 16S rRNA mtDNA fragment sequences,the examined specimens ofThelodermasp.from Tay Nguyen Plateau in central Vietnam represented a highly divergent mtDNA lineage,clearly distinct from all otherThelodermaspecies for which comparable mtDNA sequences are available(Figure 2),and a sister species ofTh.palliatumfrom Langbian Plateau of southern Vietnam. The observed differences in mtDNA sequences were congruent with evidence from the diagnostic morphological characters(see “Comparisons”).These results support our hypothesis that the small-sized“smooth”Thelodermasp.from Tay Nguyen Plateau represents a previously unknown species,which we describe herein.

    Theloderma auratum sp.nov.

    Tables 4,5;Figures 3-8.

    Holotype:ZMMU A-5828( field number NAP-06351),adult male from montane evergreen tropical forest in Kon Chu Rang Nature Reserve,Gia Lai Province,Tay Nguyen Plateau,central Vietnam(N14°30′19.7′′,E108°32′30.0′′;elevation 950 m a.s.l.),collected when calling from bush leaves ca.40 cm above the ground on 26 May 2016 at 0235 h by Nikolay A.Poyarkov.

    Paratypes:ZMMU A-5829( field number NAP-06352)and ZMMU A-5830( field number NAP-06353),two adult males from the same locality and with the same collection information as the holotype;ZMMU A-5831( field number NAP-06354),adult male from montane evergreen tropical forest in Kon Ka Kinh National Park,Gia Lai Province,Tay Nguyen Plateau,central Vietnam(N14°13′02.8′′,E108°19′53.3′′;elevation 1 420 m a.s.l.),collected ca. 1 m above the ground from vegetation on 12 May 2016 at 2300 h by Nikolay A.Poyarkov;ZMMU A-5832( field number NAP-06402),adult male from montane evergreen tropical forest in Thac Nham Forest,Kon Plong District,Mang Canh Commune,Kon Tum Province,Tay Nguyen Plateau,Vietnam(N14°43′32.4′′,E108°18′04.5′′;elevation 1 230 m a.s.l.),collected ca. 30 cm above the ground from vegetation on 5 June 2016 at 2000 h by Nikolay A.Poyarkov.

    Referred specimens:ZISP 13422-13426, five tadpoles(Gosner stages 29-36)born and reared in captive laboratory conditions from adult specimens collected from montane evergreen tropical forest in Kon Plinh Village,Kon Plong District,Kon Tum Province,Tay Nguyen Plateau,Vietnam;elevation ca.1 000 m a.s.l.)in June 2014 by Ivan I.Kropachev and Nikolai L.Orlov.

    Table 3 Uncorrected P-distance(percentage)between 12S rRNA-16S rRNA 2 524 bp fragmentsequences of Theloderma species included in phylogeneticanalyses(below the diagonal)and standard error estimates(above the diagonal)

    Diagnosis:The new species was assigned to the genusThelodermaby its(1)distinct tympanum,(2)terminal phalanx with Y-shaped distal end,(3)intercalary cartilage between terminal and penultimate phalanges of digits,(4)tips of digits expanded into large disks bearing circummarginal grooves,(5)head skin not co-ossified to skull(Liem,1970;Nguyen et al.,2015;Poyarkov et al.,2015;Rowley et al.,2011),and molecular data(Figure 2).Theloderma auratumsp.nov.is distinguished from all otherThelodermaby a combination of the following morphological attributes:(1)absence of bony ridges from canthus rostralis to occiput;(2)completely smooth skin on dorsum lacking calcified warts or asperities;(3)pointed elongated tapering snout with distinct rounded canthus rostralis,nostrils dorsolateral,not protuberant;(4)vocal opening in males absent;(5)vomerine teeth absent;(6)small body size in males(SVL 21.8-26.4 mm);(7)head longer than wide,head width to head length ratio 78%-84%;eye diameter to SVL ratio 13%-15%;snout length to SVL ratio 16%-20%;(8)comparatively small tympanum:tympanum diameter half of eye diameter(TD/EL ratio 50%-60%),tympanum with few tiny tubercles;(9)supratympanic fold absent;(10)ventral surfaces completely smooth;(11)webbing between fingers absent;(12)external and internal metacarpal tubercles present,supernumerary metacarpal tubercle single,medial,and oval in shape;(13)toes half-webbed,toe-webbing formula: I 2-2? II 1?-2? III 2-3? IV 3-1? V;(14)inner metatarsal tubercle present,oval;outer metatarsal tubercle absent;(15)iris bicolored,golden-orange dorsally,black with copper-red flecks ventrally with a wide black longitudinal stripe running medially,incorporating black pupil;(16)dorsal surfaces golden-yellow with sparse golden-orange speckling or reticulations and few small dark-brown spots(usually located at mid-dorsum,on upper eyelids,snout tip,and sacral area);(17)lateral sides of head and body with wide dark reddish-brown to black lateral stripes,clearly separated from lighter dorsal coloration by straight contrasting edge;(18)ventral surfaces of body,throat,and chestgreyish-blue with indistinct brown confluent blotches;(19)upper eyelids with few(3-5)very small flat reddish superciliary tubercles;(20)limbs dorsally reddish-brown,ventrally brown with small bluish-white speckles.

    The new species is also markedly distinct from all congeners for which comparable sequences of the 2 518-bp length 12S rRNA to 16S rRNA mitochondrial DNA fragments are available(uncorrected genetic distanceP>8.9%).

    Description of holotype:Small-sized rhacophorid frog specimen in a good state of preservation;body slender,dorsoventrally compressed(Figure 3).Skin on ventral surface of left femur of holotype medially dissected 5 mm,with femoral muscles(partial)removed for molecular genetic analyses.Measurements of holotype are given in Table 4(in mm).

    Head:Head much longer than wide(HW/HL 83%),quite deep(HD/HL 42%),moderately flattened;dorsally smooth with skin not co-ossified to skull,lacking calcified warts,asperities,or sculpturing;snout long(ESL/HL 45%)and tapering,snout tip pointed in dorsal view(Figure 3A),rounded in pro file(Figure 3D),snout notably projecting beyond margin of lower jaw(Figure 3D);nostril ovoid,not protuberant(Figure 3D),oriented dorsolaterally,located much closer to tip of snout than to eye(END/ESL 69%)(Figure 3D);canthus rostralis distinct,rounded;loreal region slightly concave;eyes large(ED/HL 36%),eye diameter less than snout length(ED/ESL 78%),notably protuberant in dorsal view(Figure 3A)and pro file(Figure 3D),pupil horizontal,ovoid(Figure 3D;Figure 4);tympanum distinct and rounded,with vertical tympanum diameter equal to horizontal(TD)diameter(ratio 1.0);tympanic rim distinct,slightly elevated above skin of temporal region(Figure 3D),tympanum comparatively small comprising half of eye diameter(TD/ED 53%),located close to eye(TED/ED 22%);pineal ocellus absent;vomerine teeth absent;choanae ovoid,almost hidden under margins of mouth roof;vocal sac openings not discernable;tongue wide,spatulate,attached anteriorly with free posterior end,with no distinct notch at posterior end;supratympanic fold indistinct,smooth(Figure 3D,Figure 4).

    Table 4 Measurements of the type series of Theloderma auratum sp.nov.(all in mm)

    Figure 3 Male holotype of Theloderma auratum sp.nov.(ZMMU A-5828)in life(ex situ)

    Forelimbs:Forelimbs thin,slender;relative finger lengths:I<II<IV<III;tips of all fingers with well-developed disks with distinct circummarginal grooves(Figure 3C),disks two times wider than finger widths(third finger disk width(FTD)205%of third finger width at penultimate phalanx),disks rounded to slightly trapezoid,slightly expanded transversally,third finger disk width 125%of third finger disk length;third finger disk width(FTD)72%of tympanum diameter(TD);dermal fringing on fingers not developed(Figure 3C); finger webbing absent;subarticular tubercles rather large,rounded,slightly protruding,distinct on all fingers, finger subarticular formula:I(1),II(1),III(2),IV(2);supernumerary metacarpal(palmar)tubercle single,medial,and ovoid;nuptial pad present,ovoid,elongated,covering prepollex area,comprising 46%of finger I length;two metacarpal tubercles present,inner metacarpal small and rounded,outer metacarpal tubercle flattened(Figure 3C).

    Figure 4 Dorsolateral view of male holotype of Theloderma auratum sp.nov.(ZMMU A-5828)in life(in situ)(Photo by Nikolay A.Poyarkov)

    Hindlimbs:Hindlimbs slender,relatively long,heels overlap when legs are at right angles to body(Figure 3A),tibiotarsal articulations reach well beyond tip of snout;tibia slightly over half of snout-vent length(TL/SVL ratio 55%);dermal ridge along outer side of tibia or tarsal fold absent;toes half-webbed,toe-webbing formula:I 2-2? II 1?-2? III 2-3? IV 3-1? V;weak dermal fringes reaching to disks of all toes.Tips of toes bearing small disks with distinct circummarginal and transverse grooves;disks rounded,notably smaller than those of fingers(HTD/FTD ratio 91%);fourth toe disk width(HTD)105%of fourth toe disk length;relative toe lengths:I<II<V<III<IV;round,clearly distinct protuberant subarticular tubercles on all toes,toe subarticular formula:I(1),II(1),III(2),IV(3),V(2);inner metatarsal tubercle well pronounced and notably protuberant,oval-shaped,2.1 times longer than wide,outer metatarsal tubercle or supernumerary tubercles absent(Figure 3E).

    Skin texture and skin glands:Dorsal skin smooth,with numerous small flat tubercles irregularly scattered on dorsal surfaces of head and body,forming weak reticulate pattern in sacral area(Figure 3A),lateral surfaces of head and body completely smooth;tympanum with few tiny evenly scattered tubercles;upper eyelids with few(3-5)very small flat reddish superciliary tubercles(Figure 3D);calcified asperities and warts on dorsum absent;area above insertion of forelimbs,and lateral sides of belly smooth;supratympanic and dorsolateral folds absent;dorsal surface of limbs weakly shagreened,ventral sides of limbs smooth;throat and chest smooth,belly and ventral surface of thigh smooth;dermal appendage at vent,dermal fringes and dermal tibiotarsal projections absent.

    Color of holotype in life:Background of dorsal surface of head,body,forearms,thighs,and shanks uniform golden-yellow with numerous tiny flat golden-orange tubercles on dorsal surfaces and supraciliary tubercles,getting somewhat denser posteriorly and forming weak reticulations in sacral area(Figure 3A;Figure 4).Small dark-brown spots of irregular shape located on dorsal surface of snout between nostrils,upper eyelids(two small spots on left and four spots on right upper eyelid),two spots on orbit margins between eyes,few small brown flecks in scapular area,two larger brown spots on mid-dorsal line(one in middle of dorsum,one in sacral area)and a single small spot above cloaca(Figure 3A).Elbows dorsally dark orange-brown.Four dark transverse dark-brown lines on dorsal surface of each thigh,orange-brown spot on each knee;on shanks a single short brown line at proximal end,larger light-brown blotch medially and smaller brown spot at distal end near tibiotarsal joint(Figure 3A).

    Lateral sides of head and body with wide dark reddish-brown to black lateral stripes,clearly separated from lighter dorsal coloration by straight contrast edge,running from snout tip along canthus rostralis to anterior eye corner,then from posterior eye corner above tympanum and posteriorly to groin;dorsal edge of wide dark lateral stripe sharp and straight until groin,which shows very shallow dark inguinal loop with rounded edges. Tympanum uniformly dark-brown with no markings;tympanal area ventrally reddish-brown,dorsally black,sharply edged from light golden-yellow supratympanic area(Figure 3D).

    Background of dorsal surfaces of arms,hands,and feet reddish-brown,with dense bluish-white to turquoise speckles,which spread to fingers and toes(Figure 3A);disks on digits dorsally brownish with dense whitish marbling,ventrally grey.Ventral surfaces of body,including belly and chest greyish-blue with indistinct brown confluent blotches,getting smaller and denser anteriorly,forming dense brownish pattern on throat(Figure 3B).Two oblique brownish spots in chest area.Ventral surfaces of forelimbs grey with whitish speckling(Figure 3B,C).Ventral surfaces of hindlimbs reddish-brown with large bluish to turquoise blotches on groin,femur,and medial part of shank(Figure 3B).Throat,chest,and ventral surfaces of thighs with small bluish-white speckles.

    Pupil horizontal,oval-shaped;iris bicolored,golden-orange dorsally,copper-red with thick black intervening ventrally;with wide black longitudinal stripe running medially,incorporating black pupil(Figure 3D).

    Color of holotype in preservative:After preservation in ethanol for two years,coloration pattern of holotype resembles that observed in life;however,yellowish and reddish tints faded completely,turning beige-grey,brownish-black patterns on lateral surfaces appeared brownish-grey;ventral coloration lost purple and bluish tints and looked beige-grey with brown spots.

    Variation:All individuals in type series were very similar in morphology,body proportions,and body coloration;measurements of type series are shown in Table 4 and representative photographs showing variation in dorsal and ventral coloration of four male paratypes in life are given in Figure 5 and Figure 6.All specimens show certain variation in number and position of small dark-brown spots on dorsum(Figure 5A,Figure 6)and ventral pattern,including size of two oblique brownish spots in chest area(Figure 5B).Coloration ofTheloderma auratumsp. nov. showed slight variation in response to diel period and microhabitat conditions. In life,coloration of dorsum was somewhat lighter nocturnally than during daytime,with dorsal surfaces looking light-beige to cream.Variation in color related to diel period,stress,and other conditions has been reported for other species ofTheloderma(McLeod&Ahmad,2007;Rowley et al.,2011).

    Tadpole description:Description of larval morphology was based on five tadpoles (Gosnerstages 29-36)(ZISP 13422-13426)(see Referred specimens for details).Identification of tadpoles was confirmed by 16S rRNA partial sequencing(see Table 2).The main morphometric parameters of the tadpoles are given in Table 5.Ontogenetic stages ofTheloderma auratumsp.nov.are shown in Figure 7.Details of tadpole morphology are presented in Figure 8.

    Figure 5 Male paratypes of Theloderma auratum sp.nov.in life in dorsal(A)and ventral(B)views(ex situ)

    Figure 6 Dorsolateral views of male paratypes of Theloderma auratum sp.nov.in life(in situ)

    Table 5 Main morphometric parameters of Theloderma auratum sp.nov.tadpoles(n=5;all in mm);developmental stages based on tables of Gosner(1960)

    Figure 7 Ontogenetic stages of Theloderma auratum sp.nov.in life(ex situ)

    External morphology:Body oval,wider than high,body height 71%-78%of body width.Body longer than wide:body width 66%-74%of body length.Snout blunt,rounded.Tail more than two times longer than body,body length 55%-65%of tail length(Figure 7A,B);23-28 myotomes discernable in lateral view.Nostrils rounded,oriented anterodorsally,located almost at same distance to snout as to eye(RN/NP ratio 71%-105%).Internarial distance(IND)56%-77%of interorbital distance(IOD).Eyes with dorsal orientation. Eye diameter(ED)9%-18%of body length(BL).Eye-nostril distance(END)larger than eye diameter(ED).Spiracle single,sinistral,located closer to posterior portion of body,directed laterally.Distance from tip of snout to spiracle opening 67%-75%times body length.Vent tube in medial position,with aperture located on same line as margin of ventral tail fin.Height of tail musculature at highest portion 58%-74%of tail height and 50%-57%of maximum body height.Maximum height of dorsal tail fin 34%-42%of maximum tail height.Ventral tail fin same height as dorsal tail fin;maximum height of lower tail fin 35%-46%of maximum tail height.Ventral and dorsal tail fins start roughly at level of vent.Dorsal and ventral tail fins slightly higher in posterior one third of tail length;tail tip gently rounded(Figure 7C;Figure 8A).

    Oral disk:Mouth with ventral orientation(Figure 8C),oral disk wide,elliptical,width 38%-44%of body width.Lower labium fringed with double row of finger-shaped papillae(Figure 8D).Submarginal papillae located laterally of upper and lower labium. Number of marginal papillae in one row 80-85.LTRF:5(3-5)/3. Number of keratodonts 24-28 on 1 mm.Lower sheath “V”-shaped;upper sheath “U”-shaped,both with distinctly serrated edges(Figure 8D).

    Coloration in life:Body strongly pigmented:dorsally and ventrally uniformly dark-grey(Figure 7C,D).In dorsal view,eyes not totally visible on strongly pigmentated background(Figure 7C).Iris black.Spiracle slightly pigmented.Dorsal and ventral tail fins translucent at edges with greyish coloration,much lighter than on body.

    Figure 8 Tadpole morphology of Theloderma auratum sp.nov.(Gosner stage 31)

    Coloration in preservative:In ethanol,all parts of tadpoles were less intensively pigmented and turned greyish. In preservation,myotomes of tail muscles were more discernable than in life.

    Advertisement call:Advertisement calls ofTheloderma auratumsp.nov.were uttered in a series of 14.27±1.31(7-21,n=15)tonal calls(Figure 9),resembling an orthopteran call to the human ear.Series duration varied from 3.04 s to 11.34 s(7.88±0.62 s,n=15)and the interval between successive series comprised 21.61±2.57 s(12.16-44.4 s,n=13).Call duration varied from 30 ms to 75 ms(61±0.5 ms,n=214)and inter-call interval varied from 411 ms to 842 ms(529±6.76 ms,n=199)usually decreasing gradually from the beginning to end of the series(Figure 9). Call repetition rate within the series was 1.68±0.03 calls/s(1.53-1.98 calls/s,n=15).Initial fundamental frequency was 2 777±3 Hz(2 570-2 900 Hz,n=214)and final fundamental frequency was 2 856±4 Hz(2 710-3 090 Hz,n=214).The initial and final fundamental frequencies also represented the minimum and maximum fundamental frequencies,respectively.Frequency modulation was usually expressed in the weak lift of fundamental frequency during the whole call. The presence of harmonics varied between/within recordings and mostly depended on recording quality(e.g.,sensitivity of recording equipment,distance from vocalizing animal,signal volume,and background noise).Calls from the highest quality recording contained poorly visible second harmonics(Figure 9). The maximum amplitude frequency varied from 2 760 to 2 950 Hz(2 829±2 Hz,n=214)and coincided with the fundamental frequency.

    Figure 9 Oscillograms(top)and sonograms(bottom)of male advertisement calls of holotype of Theloderma auratum sp.nov.(ZMMU A-5828),recorded at 21.5°C(Kon Chu Rang Nature Reserve,Gia Lai Province,Tay Nguyen Plateau,central Vietnam)

    Position in mtDNA genealogy and sequence divergence:According to our mtDNA data,Theloderma auratumsp.nov.belongs to the subgenusThelodermas. str. (Poyarkov et al.,2015),and is reconstructed as a sister species ofTh.palliatum(Figure 2,clade II-B).Uncorrected geneticP-distances betweenTheloderma auratumsp.nov.12S rRNA and 16S rRNA sequences and all homologous sequences of congeners included in our analyses varied from 8.9%(with sister speciesTh.palliatum)to 17.9%(withTh.laeve)(see Table 3).This value is higher than that observed between several currently recognized species ofTheloderma(Table 3).

    Distribution and biogeography:The known distributions ofTheloderma auratumsp. nov. and its sister speciesTh.palliatumare shown in Figure 1.To date,the new species is known from montane evergreen tropical forests of Tay Nguyen Plateau in the central Annamite(Truong Son)Mountains,and has been recorded in Gia Lai,Kon Tum,and Thua Thien-Hue provinces.It is anticipated thatTheloderma auratumsp.nov.also occurs in the adjacent montane forests of Tay Nguyen Plateau;in particular,records from Quang Nam Province of Vietnam and Xekong Province of Laos are anticipated.

    Natural history notes:Our knowledge on the biology ofTheloderma auratumsp.nov.is scarce.The new species was recorded in primary polydominant tropical montane evergreen forests of Tay Nguyen Plateau at elevations ranging from 800 to 1 400 m a.s.l..Animals were recorded only in patches of primary undisturbed forest with complete multi-layered canopy and heavy undergrowth,suggesting the new species is a strict forest-dwelling specialist. At the type locality in Kon Chu Rang Nature Reserve(Gia Lai Province),the forest where the new species was recorded is dominated by large trees of the families Podocarpaceae(Dacrydium elatum,Dacrycarpus imbricatus),Magnoliaceae,Burseraceae (Canariumsp.),Myrtaceae(Syzygiumsp.),Hamamelidaceae(Simingtoniasp.),Lauraceae(Litsiasp.),Rhodoliaceae(Rhodoliasp.),Fagaceae,Sterculiaceae(Scaphiumsp.)(Figure 10).

    Theloderma auratumsp.nov.is a very secretive species with apparently nocturnal activity. Animals were usually encountered at night(between 1900 h and 0300 h)after or during heavy rain;males were recorded when calling from shrubs or small trees ca.20-90 cm above the ground.Diet of the new species remains unknown.In Kon Plong District,Kon Tum Province,both males and females were observed at an elevation of 1 000 m a.s.l.when sitting on leaves of shrubs at 40-50 cm from the ground;animals were active at an air temperature of 20-22°C and 100%humidity.

    Other species of anurans recorded syntopically with the new species at the type locality includedIngerophrynus galeatus(Günther,1864),Kurixalusbanaensis(Bourret,1939),Rhacophorus annamensisSmith,1924,Rhacophorus rhodopusLiu&Hu,1960,Rh.robertingeriOrlov,Poyarkov,Vassilieva,Ananjeva,Nguyen,Nguyen&Geissler,2012,Polypedates mutus(Smith,1940),Rana johnsiSmith,1921,Sylvirana nigrovittata(Blyth,1856),Microhyla pulverataBain&Nguyen,2004,andOccidozyga vittata(Andersson,1942).In Kon Ka Kinh National Park(Gia Lai Province),the new species was recorded in sympatry with three otherThelodermaspecies,includingTh.albopunctatum,Th.vietnamense,andTh.truongsonense,with the latter species sharing the same biotopes asTheloderma auratumsp.nov.In Thac Nham Forest(Kon Plong Province),the new species was recorded in sympatry withTh.albopunctatum,Th.gordoni,andTh.ryabovi.

    Reproductive biology:Reproduction of the new species was observed under laboratory conditions. In terrarium,clutches ofTheloderma auratumsp. nov. were found on the under-surface of snags(Figure 7)and were usually deposited 2-4 cm above the surface of the water.A clutch consisted of one or two,rarely three eggs,enveloped in slimy transparent external egg-capsules(see Figure 7A,B).Duration of embryonic development comprised 10-12 d at ambient air temperatures of 22-23°C.Tadpoles were fed on dry fi sh food and started active feeding 2-3 d after hatching.Duration of larval development from hatching to the end of metamorphosis was about 2.5 months.At metamorphosis,the main aspects of coloration and dorsal pattern corresponded to that of the adults,except the dorsal coloration pattern:metamorphs had a more intensive orange dorsum with orange reticulate pattern(Figure 7E),which disappeared with age.

    Comparisons:Theloderma auratumsp. nov.is easily distinguished from most species of the genusThelodermaby a combination of the following morphological attributes:(1)small to medium body size(SVL 21.8-26.4 mm);(2)absence of vomerine teeth;(3)smooth skin on dorsum with dorsal warts or asperities absent,(4)long obtuse snout with pointed snout tip,ESL/SVL ratio 16%-20%;(5)uniform golden-yellow dorsal coloration with weak golden-orange speckling and reticulations;(6)dark reddish-brown to black lateral stripes,separated from dorsal light coloration by sharp edge,forming very weak and shallow inguinal loop;and(7)bicolored iris,golden-orange dorsally,copper-red with thick black intervening ventrally, with a wide black longitudinal medial stripe.

    The small body size,absence of vomerine teeth,and smooth dorsal skin lacking asperities or warts distinguishesTheloderma auratumsp. nov. from the large-sizedThelodermaspecies:Th.bicolor(Bourret,1937),Th.corticale(Boulenger,1903),Th.gordoniTaylor,1962,Th.leporosumTschudi,1838,Th.moloch(Annandale,1912),andTh.nagalandenseOrlov,Dutta,Ghate et Kent,2006(vs.large body size,large dorsal asperities,and vomerine teeth present).

    Absence of hand webbing,reduced foot webbing,and completely smooth dorsum clearly distinguishes the new species fromTh.horridum(Boulenger,1903),Th.stellatumTaylor,1962,Th.vietnamensePoyarkov,Orlov,Moiseeva,Pawangkhanant,Ruangsuwan,Vassilieva,Galoyan,Nguyen&Gogoleva,2015,Th.ryaboviOrlov,Dutta,Ghate&Kent,2006,Th.phrynoderma(Ahl,1927),Th.asperum(Boulenger,1886),andTh.albopunctatum(Liu&Hu,1962)(vs.hand webbing well-developed inTh.horridum,rudimentary in remaining taxa;vs.foot webbing complete,dorsum with more or less developed warts and asperities in above-mentioned taxa).

    The Malayan speciesTh.licinandTh.asperum,Indochinese speciesTh.albopunctatum,and poorly knownTh.baibungense(Jiang,Fei&Huang,2009)andTh.pyaukkyaDever,2017 can be further distinguished fromTheloderma auratumsp. nov. by whitish coloration of whole dorsum or large white blotches on dark dorsal background color(vs.golden-yellow dorsum with golden-orange reticulations and few dark spots inTheloderma auratumsp. nov.),greyish or whitish belly with brown reticulations(vs.whitish to bluish belly with small brown spots and blotches inTheloderma auratumsp.nov.),and uniform red coloration of iris(vs.distinctly bicolored golden/black iris inTheloderma auratumsp.nov.).Furthermore,Th.licinMcLeod&Ahmad,2007 from Sundaland and the Malayan Peninsula can be distinguished from the new species by the uniform dark-red coloration of iris(vs.distinctly bicolored golden/black iris inTheloderma auratumsp.nov.)and short rounded snout with ESL/HL ratio below 40%(vs.long pointed snout with ESL/HL ratio over 43%inTheloderma auratumsp.nov.).

    Theloderma annaeNguyen,Pham,Nguyen,Ngo&Ziegler,2016 can be easily differentiated from the new species by a golden-green iris with black reticulations and greenish-grey dorsal pattern (vs. distinctly bicolored golden/black iris and golden-yellow dorsum with golden-orange reticulations inTheloderma auratumsp. nov.) and presence of supernumerary palmar tubercles(vs.absent inTheloderma auratumsp.nov.).Theloderma auratumsp.nov.can be distinguished from two other small-sizedThelodermaspecies found in northern Vietnam and southern China-Th.rhododiscus(Liu&Hu,1962)andTh.lateriticumBain,Nguyen&Doan,2009-by golden-yellow dorsum with golden-orange reticulations and few dark spots(vs.tea brown to dark grey dorsum inTh.rhododiscusand dorsum with deep brick-red background color with distinct black middorsal spot or blotches inTh.lateriticum),greyish belly or whitish belly with brown reticulations(vs.belly brownish-black scattered with grey-white network in both taxa),ventral surfaces of disks grey(vs.orange or red inTh.rhododiscus),completely smooth dorsum without asperities or warts(vs.distinctly white dorsal asperities present in both taxa),and distinctly bicolored golden/black iris(vs.uniformly dark reddish brown iris in both taxa).Theloderma lacustrinumSivongxay,Davankham,Phimmachak,Phoumixay&Stuart,2016 from Laos can be differentiated from the new species by its uniform red iris(vs.distinctly bicolored golden/black iris),beige-brown upper jaw(vs.uniformly dark reddish-brown upper jaw),and presence of large black blotches in inguinal region(vs.dark inguinal blotches absent).

    Thelodermapetilum(Stuart& Heatwole,2004)can be easily differentiated from the new species by two chocolate-brown bands running from lateral surfaces of head toward groin,distinctly separated from brownish-beige dorsum and edged with white(vs. dark lateral bands black to reddish-brown,sharply contrasting with golden-orange dorsum with golden-orange reticulations inTheloderma auratumsp.nov.).

    Superficially,Theloderma auratumsp.nov.morphologically resembles other small-bodiedThelodermaspecies found in central and southern Vietnam,namelyTh.laeve(Smith,1924),Th.nebulosumRowley,Le,Hoang,Dau&Cao,2011,Th.truongsonense(Orlov&Ho,2005),andTh.palliatumRowley,Le,Hoang,Dau&Cao,2011.Thus,comparisons with these species appear to be the most pertinent:Th.laeve(Smith,1924),Th.nebulosumRowley,Le,Hoang,Dau&Cao,2011,Th.truongsonense(Orlov&Ho,2005),andTh.palliatumRowley,Le,Hoang,Dau&Cao,2011(Figure 11).Table 6 summarizes the morphometric data on the small-bodiedThelodermaspecies of central and southern Vietnam.

    Figure 1 1 Small-bodied Theloderma species of central and southern Vietnam

    Theloderma nebulosum(Figure 11D)can be distinguished from the new species by small sparse dorsal asperities present(vs. absent),brown dorsum with darker patterning(vs.golden-yellow dorsum with golden-orange reticulations and few dark spots),and dark brown to black ventral surfaces with bluish marbling(vs.whitish to bluish belly with small brown spots and blotches).

    Theloderma truongsonense(Figure 11A)from the central and southern Annamite(Truong Son)Mountains can be distinguished fromTheloderma auratumsp. nov. by the following morphological attributes:dorsal skin shagreened with numerous small asperities scattered on dorsum and dorsal surfaces of head including upper eyelids(vs. completely smooth dorsum with only few(3-5) flat superciliary tubercles),comparatively shorter snout with rounded snout tip,ESL/SVL ratio below 16% (vs. tapering snout with pointed tip,comparatively longer,ESL/SVL ratio above 17%),dorsum brown with darker patterning(vs.golden-yellow dorsum with golden-orange reticulations and few dark spots),dorsal edge of dark lateral stripe irregular,deep inguinal loop or large black inguinal blotch surrounded by bluish to turquoise edging always present(vs. dark lateral stripe with smooth sharp edge,inguinal loop absent or shallow),comparatively wider and shorter head,HW/HL ratio above 85%(vs.comparatively narrower and longer head,HW/HL ratio below 84%),and comparatively shorter tibia,TL/SVL ratio below 52%(vs.comparatively longer tibia,TL/SVL ratio above 53%).

    Theloderma laeve(Figure 11C)from low to mid-elevations of Langbian Plateau,southern Vietnam,can be diagnosed fromTheloderma auratumsp.nov.by the following combination of morphological characters:dorsum beige with brownish patterning and thin light middorsal stripe(vs.golden-yellow dorsum with golden-orange reticulations and few dark spots),lateral stripe brownish-violet(vs.black to dark brown lateral stripe),ventral surfaces uniform violet-grey with no dark patterning(vs. whitish to bluish belly with small brown spots and blotches),head short,HL/SVL ratio below 35%(vs. head comparatively longer,HL/SVL ratio over 37%),head width slightly less than head length,HW/HL ratio 86%-93%(vs. head notably longer than wide,HW/HL ratio 78%-84%),comparatively shorter snout with rounded snout tip,ESL/SVL ratio 14%-16%(vs.tapering snout with pointed tip,comparatively longer,ESL/SVL ratio 17%-20%),and comparatively smaller interorbital distance,IOD/SVL ratio 9%-10%(vs.IOD/SVL ratio 11%-14%).

    From its sister species,Th.palliatum(Figure 11B),which inhabits high elevations of the Langbian Plateau,Theloderma auratumsp.nov.can be distinguished by the following morphological features:completely smooth dorsum(vs.notably tuberculated dorsum,with dorsal surfaces of body,head,and limbs covered with many tubercles and asperities),golden-yellow dorsum with golden-orange reticulations and few dark spots(vs.beige dorsum with large light-brown blotches on upper eyelids,scapular,and sacral area),dark lateral stripe with smooth sharp edge,inguinal loop shallow or absent,large black inguinal blotch absent(vs.dorsal edge of dark lateral stripe irregular,large black inguinal blotch or dark inguinal loop present),whitish to bluish belly with small brown spots and blotches(vs.almost black belly with thin bluish reticulations),snout dorsally pointed(vs.snout dorsally obtusely rounded),comparatively smaller tympanum,TD/SVL ratio 7%-8%(vs.comparatively larger tympanum,TD/SVL ratio 9%-10%),and comparatively larger eye to nostril distance,END/SVL ratio 11%-13%(END/SVL ratio 9%-10%).

    Tadpoles ofTheloderma auratumsp.nov.superficially resemble the tadpoles ofTh.nebulosumdescribed by Rowley et al.(2011);however,they can be distinguished from the latter by LTRF 5(3-5)/3(vs.LTRF 4(2-4)/3 inTh.nebulosum).Tadpoles ofTh.laeveandTh.palliatumalso exhibit LTRF 4(2-4)/3(see Orlov et al.,2012).

    Etymology:The specific name “auratum”is a Latin adjective in the nominative singular(neutral gender),derived from Latin“aurum”for “gold”,referring to the golden-yellowish dorsal coloration of the new species.

    Recommended vernacular name:We recommend the following common name in English:Golden Bug-Eyed Frog.Recommended vernacular name in Vietnamese:′êch Cay S`an Vàng.

    Conservation status:Theloderma auratumsp.nov.is,to date,known from five localities in the Tay Nguyen Plateau of central Vietnam.Further research is required to estimate its actual distribution,population trends,and possible threats.It appears that the new species is associated with primary undisturbed montane forests and may be affected by growing anthropogenic pressure and forest destruction,as observed in different areas of central Vietnam. Given the available information,we suggestTheloderma auratumsp.nov.be tentatively considered as a Data Deficient species following IUCN’s Red List categories(IUCN,2001).

    DISCUSSION

    Our phylogenetic data largely confirmed the phylogeny ofThelodermaas given by Poyarkov et al.(2015),Nguyen et al.(2015),and Sivongxay et al.(2016),but included several taxa not analyzed in these studies.Our data supported subdivision ofThelodermainto two major clades,corresponding to the subgeneraStelladerma(Clade I)andThelodermas. str.(Clade II).Furthermore,Th.annaewas reconstructed as a sister species ofTh.nebulosumwith high node support,andTh.moloch,Th.phrynoderma,andTh.ryaboviwere supported as a monophyletic sister group to theTh.asperumcomplex.

    Our research clearly demonstrated that current understanding of the diversity ofThelodermais far from complete.WithinTh.asperum,our data strongly suggested paraphyly ofTh.pyaukkyafrom Myanmar with respect toTh.baibungensefrom the eastern Himalayas.In addition,Th.pyaukkyacould be subdivided into two highly divergent non-monophyletic lineages from northern and central Myanmar(P-distance 3.9%).Dever(2017)describedTh.pyaukkyawithout includingTh.baibungensein phylogenetic analysis and without examination ofTh.baibungensespecimens;comparisons ofTh.pyaukkyain this paper do not allow the discrimination ofTh.pyaukkyafromTh.baibungense.Our data indicated that the taxonomic status ofTh.pyaukkyalineages requires careful reconsideration:synonymy of this species withTh.baibungensecan be assumed;however,high genetic differentiation between its lineages suggests that taxonomic reassessment of this group is necessary.Other examples of high intraspecific genetic distances included Indochinese speciesTh.lateriticum(P=3.6%),Th.truongsonense(P=3.8%),andTh.albopunctatum(P=2.5%);these results suggest the possible presence of cryptic species in these complexes and incomplete taxonomy of the group.Finally,Thelodermasp. from Gia Lai Province,previously identified asTh.laeveby Nguyen et al. (2015),clearly represented a deeply divergent mtDNA lineage sister toTh.laeves. str. (P=9.0%)and likely corresponds to a yet undescribed species.

    Table 6 Morphometric differentiation ofsmall-bodiedThelodermaspecies of central and southern Vietnam,resemblingThelodermaauratum sp.nov.

    Our study confirmed significant species richness underestimation inThelodermaand reported on a new species ofThelodermafrom the Tay Nguyen Plateau of the Annamite(Truong Son)Mountains in central Vietnam.This mountainous area harbours the highest diversity of amphibian species in Indochina(Geissler et al.,2015)and is recognized as a local center of herpetofaunal diversity with many new species of amphibians and reptiles described from the plateau in the last decade(Poyarkov et al.,2014,2017;Rowley et al.,2010,2011,2014,2016;Nguyen et al.,2018,Nguyen et al.,in press). Our data indicate that knowledge on amphibian diversity of this area is still far from complete and further diversity is likely to be revealed with additional survey efforts.Habitat loss and modification is a continued threat in the region(Meyfroidt&Lambin,2008),and is widely recognized as the greatest threat to amphibians in Southeast Asia.Range-restricted species are at greatest risk,and includeTheloderma auratumsp.nov.,which appears to be a forest specialist restricted to relatively undisturbed broadleaf evergreen montane forests.Further survey efforts and additional taxonomic studies to understand the true diversity of amphibians in the region are urgently required for effective conservation management.

    COMPETING INTERESTS

    The authors declare that they have no competing interests.

    AUTHORS’CONTRIBUTIONS

    N.A.P.and N.L.O.designed the study. N.A.P.and I.I.K.collected data.N.A.P.performed molecular experiments. N.A.P.and I.I.K.examined morphology.N.A.P.conducted phylogenetic analyses.S.S.G.conducted acoustic analysis.N.A.P.,I.I.K.,and S.S.G.wrote the manuscript,N.L.O.revised the manuscript.All authors read and approved the final manuscript.

    ACKNOWLEDGEMENTS

    Fieldwork was funded by the Joint Russian-Vietnamese Tropical and Technological Center and was conducted underpermission of the Department of Forestry,Ministry of Agriculture,and Rural Development of Vietnam(permit numbers 547/TCLN-BTTN,issued 21.04.2016,and 432/TCLN-BTTN,issued 30.03.2017).Forest Protection Departments of the Peoples’Committee of Gia Lai Province provided permits for fieldwork and sample collection(permit numbers 1951/UBND-NV,issued 04.05.2016,and 142/SNgV-VP,issued 11.04.2017).The authors are grateful to Andrey N.Kuznetsov and Leonid P.Korzoun for support and organization of the fieldwork. We are grateful to Valentina Kretova(Moscow)and Alexey A.Evsyunin(Tula Exotarium)for help with preparation of figures for this paper.Many thanks to Evgeniy Popov,Alexander Ostorshabov,Igor V.Palko,Evgeniya N.Solovyeva,and Duong Van Tang for their support during fieldwork and assistance in the laboratory. We are grateful to Evgeny Rybaltovsky and Marina Rybaltovskaya(Zoocom)for their help with animal care.We sincerely thank our Vietnamese colleagues Nguyen Dang Hoi,Hoang Minh Duc,Pham Thi Ha Giang,and Le Xuan Son for help and continued support.For permission to study specimens under their care and permanent support,we thank Valentina F.Orlova(ZMMU)and Natalia B.Ananjeva(ZISP).We are sincerely grateful to Samuel Shonleben and Egill Scallagrimsson for their kind help and useful comments,which helped us improve this manuscript.

    REFERENCES

    Ahl E.1927.Zur Systematik der asiatischen arten der froschgattung Rhacophorus.Sitzungsberichte der Gesellschaft Naturforschender Freunde zu Berlin,1927:35-47.

    Ahl E.1931.Amphibia,anura III,polypedatidae.Das Tierreich,55:477.

    Altig R,McDiarmid RW.1999.Body plan:development and morphology.In:McDiarmid RW,Altig R.Tadpoles:the Biology of Anuran Larvae.Chicago,IL:University of Chicago Press,24-51.

    Altig R.2007.A primer for the morphology of anuran tadpoles.Herpetological Conservation and Biology,2(1):71-74.

    AmphibiaWeb.2018.Amphibia Web provides Information on amphibian declines,natural history,conservation,and taxonomy.California:Berkeley,http://amphibiaweb.org/.(2018-01-15)

    Annandale N.1912.Zoological results of the Abor Expedition,1911-1912.I.Batrachia.Records of the Indian Museum,8:7-36.

    Bain RH,Nguyen TQ,Doan KV.2009.A new species of the genusThelodermaTschudi,1838(Anura:Rhacophoridae)from northwestern Vietnam.Zootaxa,2191:58-68.

    Biju SD,Senevirathne G,Garg S,Mahony S,Kamei RG,Thomas A,Shouche YS,Raxworthy CJ,Meegaskumbura M,Van Bocxlaer I.2016.Frankixalus,a new rhacophorid genus of tree hole breeding frogs with oophagous tadpoles.PLoS One,11(1):e0145727.

    Boulenger GA.1882.V.—Description of a new genus and species of frogs of the family Ranidae.Annals and Magazine of Natural History,10(55):35.

    Boulenger GA.1886.First report on additions to the batrachian collection in the Natural History Museum.Proceedings of the Zoological Society,1886:411-416.

    Boulenger GA.1893.Concluding report on the reptiles and batrachians obtained in Burma by Signor L.Fea dealing with the collection made in Pegu and the Karin Hills in 1887-88.Annali del Museo Civico di Storia Naturale di Genova,2(13):304-347.

    Boulenger GA.1903.IX.—Descriptions of three new batrachians from Tonkin.Annals and Magazine of Natural History,12(67):186-188.

    Bourret R.1937.Notes herpétologiques sur l’Indochine fran?aise.XIV.les batraciens de la collection du laboratoire des Sciences naturelles de l’université.descriptions de quinze especes ou variétés nouvelles.Annexe au Bulletin Général de l’Instruction Publique,Hanoi,1937:5-56.

    BourretR.1942.Les Batraciens de l’Indochine.Hanoi: Institut Océanographique de l’Indochine.

    Chan-ard T.2003.A Photographic Guide to Amphibians in Thailand.Bangkok,Thailand:Darnsutha Press Co.,Ltd.(In Thai)

    Che J,Wang K.2016.Amphibia China:an online database of Chinese Amphibians.Zoological Research,37(1):57-59.

    Das I.2007.Amphibians and Reptiles of Brunei:A Pocket Guide.Kota Kinabalu,Sabah,Malaysia:Natural History Publications(Borneo).

    Das I,Grinang J,Min PY.2013.Theloderma licinMcLeod and Ahmad,2007(Lissamphibia:Anura:Rhacophoridae),a new frog record for Borneo.Asian Herpetological Research,4(3):221-223.

    Dever JA,Nguyen H,Wilkinson JA.2015.Rediscovery and redescription ofTheloderma phrynoderma(Ahl,1927)(Anura:Rhacophoridae)from Myanmar.Copeia,103(2):402-415.

    Dever JA.2017.A new cryptic species of theTheloderma asperumcomplex(Anura:Rhacophoridae)from Myanmar.Journal of Herpetology,51(3):425-436.

    Fei L,Ye CY,Huang YZ.1990.Key to Chinese Amphibia.Chongqing,China:Publishing House for Scientific and Technological Literature.(in Chinese)

    Fei L,Hu SQ,Ye CY,Huang YZ.2009.Fauna Sinica.Amphibia.Volume 2,Anura Ranidae.Beijing:Science Press.(in Chinese)

    Fei L,Ye CY,Jiang JP.2010.Colored Atlas of Chinese Amphibians.Chengdu:Sichuan Publishing House of Science and Technology.(in Chinese)

    Felsenstein J.1985.Confidence limits on phylogenies:an approach using the bootstrap.Evolution,39(4):783-791.

    Felsenstein J.2004.Inferring phylogenies.Sunderland,Massachusetts:Sinauer Associates,465.

    Frost DR.2018.Amphibian species of the world 6.0,An Online Reference.New York:American Museum of Natural History,http://research.amnh.org/herpetology/amphibia/index.php/.(2018-01-15)

    Geissler P,Hartmann T,Ihlow F,R?dder D,Poyarkov NA Jr,Nguyen TQ,Ziegler T,B?hme W.2015.The Lower Mekong:an insurmountable barrier to amphibians in southern Indochina?Biological Journal of the Linnean Society,114(4):905-914.

    Gosner KL.1960.A simplified table for staging anuran embryos and larvae with notes on identification.Herpetologica,16(3):183-190.

    Grosjean S.2001.The tadpole ofLeptobrachium(Vibrissaphora)echinatum(Amphibia,Anura,Megophryidae).Zoosystema,23(1):143-156.

    Hall TA.1999.BioEdit:a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.Nucleic Acids Symposium Series,41(2):95-98.

    Hedges SB.1994.Molecular evidence for the origin of birds.Proceedings of the National Academy of Sciences of the United States of America,91(7):2621-2624.

    Hendrix R,Grosjean S,Le KQ,Vences M,Vu NT,Ziegler T.2007.Molecular identification and description of the tadpole of the Annam Flying Frog,Rhacophorus annamensisSMITH,1924(Anura:Rhacophoridae).Salamandra,43(1):11-19.

    Hillis DM,Moritz C,Mable BK.1996.Molecular Systematics:Molecular Systematics.2nded.Sunderland,Massachusetts:Sinauer Associates,655.Hou M,Yu GH,Chen HM,Liao CL,Ziang L,Chen J,Li PP,Orlov NL.2017.The taxonomic status and distribution range of sixThelodermaspecies(Anura:Rhacophoridae)with a new record in China.Russian Journal of Herpetology,24(2):99-127.

    Huelsenbeck JP,Hillis DM.1993.Success of phylogenetic methods in the four-taxon case.Systematic Biology,42(3):247-264.

    Huelsenbeck JP,Ronquist F.2001.MRBAYES:Bayesian inference of phylogenetic trees.Bioinformatics,17(8):754-755.

    Inger RF.1954.Systematics and zoogeography of Philippine Amphibia.Fieldiana:Zoology,33(4):183-531.

    Inger RF.1966.The systematics and zoogeography of the Amphibia of Borneo.Fieldiana:Zoology,52:1-402.

    IUCN.2001.IUCN Red List Categories and Criteria:Version 3.1.Gland,Switzerland and Cambridge,UK:IUCN Species Survival Commission,30.

    Jiang K,Yan F,Wang K,Zhou DH,Li C,Che J.2016.A new genus and species of treefrog from Medog,southeastern Tibet,China(Anura,Rhacophoridae).Zoological Research,37(1):15-20.

    JobbG,Von HaeselerA,StrimmerK.2004.TREEFINDER:a powerful graphical analysis environment for molecular phylogenetics.BMC Evolutionary Biology,4:18.

    Anderson J.1879.Anatomical and Zoological Researches:Comprising an Account of the Zoological Results of the Two Expeditions to Western Yunnan in 1868 and 1875;and a Monograph of the Two Cetacean Genera Platanista and Orcella.Volume V.2.London:Bernard Quaritch.

    Kumar S,Stecher G,Tamura K.2016.MEGA7:molecular evolutionary genetics analysis version 7.0 for bigger datasets.Molecular Biology and Evolution,33(7):1870-1874.

    Lalronunga S,Lalrinchhana C.2017.Morphologicaland molecular characterization ofTheloderma moloch(Anura: Rhacophoridae)from Indo-Burma biodiversity hotspot of northeast India.Science Vision,17(3):148-159.

    Li C,Hou M,Yu GH,Yan F,Li BZ,Jiang K,Li PP,Orlov NL.2016.Rediscovery ofTheloderma moloch(Amphibia:Anura:Rhacophoridae)from Southeast Tibet,China after 102 years.Russian Journal of Herpetology,23(1):41-54.Li JT,Che J,Bain RH,Zhao EM,Zhang YP.2008.Molecular phylogeny of Rhacophoridae(Anura):a framework of taxonomic reassignment of species within the generaAquixalus,Chiromantis,Rhacophorus,andPhilautus.Molecular Phylogenetics and Evolution,48(1):302-312.

    Li JT,Che J,Murphy RW,Zhao H,Zhao EM,Rao DQ,Zhang YP.2009.New insights to the molecular phylogenetics and generic assessment in the Rhacophoridae(Amphibia:Anura)based on five nuclear and three mitochondrial genes,with comments on the evolution of reproduction.Molecular Phylogenetics and Evolution,53(2):509-522.

    Li JT,Chen YY,Li SQ,Lv K,Wang YZ.2011.Catalogue of the type specimens of amphibians and reptiles in the herpetological museum of Chengdu Institute of Biology,Chinese Academy of Sciences: I.Rhacophoridae(Anura,Amphibia).Asian Herpetological Research,2(2):129-141.

    Li JT,Li Y,Klaus S,Rao DQ,Hillis DM,Zhang YP.2013.Diversification of rhacophorid frogs provides evidence for accelerated faunal exchange between India and Eurasia during the Oligocene.Proceedings of the National Academy of Sciences of the United States of America,110(9):3441-3446.

    Liem SS.1970.The morphology,systematics,and evolution of the old world treefrogs(Rhacophoridae and Hyperoliidae).Fieldiana:Zoology,57:1-145.Liu CC,Hu SQ.1962.A herpetological report of Kwangsi.Acta Zoologica Sinica,14(S1):73-104.

    Manthey U,Grossmann W.1997.Amphibien und Reptilien Südostasiens.Münster:Natur und Tier-Verlag,512.

    Mathew R,Sen N.2010.Pictorial Guide to Amphibians of North East India.Kolkata,India:Zoological Survey of India.

    McLeod DS,Ahmad N.2007.A new species ofTheloderma(Anura:Rhacophoridae)from southern Thailand and peninsular Malaysia.Russian Journal of Herpetology,14(1):65-72.

    Meegaskumbura M,Senevirathne G,Biju SD,Garg S,Meegaskumbura S,Pethiyagoda R,Hanken J,SchneiderCJ.2015.Patterns of reproductive-mode evolution in old world tree frogs(Anura,Rhacophoridae).Zoologica Scripta,44(5):509-522.

    Meyfroidt P,Lambin EF.2008.The causes of the reforestation in Vietnam.Land Use Policy,25(2):182-197.

    Milto KD,Poyarkov NA Jr,Orlov NL,Nguyen TT.2013.Two new Rhacophorid frogs from Cat Ba Island,Gulf of Tonkin,Vietnam.Russian Journal of Herpetology,20(4):287-300.

    Nguyen SN,Yang JX,Le TNT,Nguyen LT,Orlov NL,Hoang CV,Nguyen TQ,Jin JQ,Rao DQ,Hoang TN,Che J,Murphy RW,Zhang YP.2014.DNA barcoding of Vietnamese bent-toed geckos(Squamata:Gekkonidae:Cyrtodactylus)and the description of a new species.Zootaxa,3784(1):48-66.

    Nguyen LT,Poyarkov NA Jr,Le DT,Vo BD,Phan HT,Van Duong T,Murphy RW,Nguyen SN.2018.A new species ofLeptolalax(Anura:Megophryidae)from Son Tra Peninsula,central Vietnam.Zootaxa,4388(1):1-21.

    Nguyen TL,Poyarkov NA Jr,Nguyen TT,Nguyen TA,Nguyen VH,Gorin VA,Murphy RW,Nguyen SN.2018.A new species of the genusMicrohylaTschudi,1838(Amphibia:Anura:Microhylidae)from Tay Nguyen Plateau,Central Vietnam.Zootaxa,in press.

    Nguyen TT,Matsui M,Eto K.2015.Mitochondrial phylogeny of an Asian tree frog genusTheloderma(Anura:Rhacophoridae).Molecular Phylogenetics and Evolution,85:59-67.

    Chan K,Ahmad N.2009.Distribution and natural history notes on some poorly known frogs and snakes from Peninsular Malaysia.Herpetological Review,40(3):294-301.

    Orlov NL,Ho TC.2005.A new species ofPhilautusfrom Vietnam(Anura:Rhacophoridae).Russian Journal of Herpetology,12(2):135-142.

    Orlov NL,Dutta SK,Ghate HV,Kent Y.2006.New species ofThelodermafrom Kon Tum Province(Vietnam)and Nagaland State(India)[Anura:Rhacophoridae].Russian Journal of Herpetology,13(2):135-154.

    Orlov NL.2007.Redescription ofendemic species from Central Annam,Theloderma ryaboviOrlov,Dutta,Ghate et Kent,2006(Anura:Rhacophoridae)and some aspects of its natural history.Sovremennaya Gerpetologiya,7(1-2):98-110(in Russian).

    Orlov NL,Ananjeva NB.2007.Amphibians of South-East Asia.St.Petersburg,Russia:St.Petersburg University Press,270.(in Russian)

    Orlov NL,Ryabov SA,Ananjeva NB,Evsyunin AA.2010.Asian tree frogs of genusThelodermaTschudi,1838(Amphibia:Anura:Rhacophoridae:Rhacophorinae).St.Petersburg:Russian Academy of Sciences,52.

    Orlov NL,Poyarkov NA Jr,Vassilieva AB,Ananjeva NB,Nguyen TT,Nguyen VS,Geissler P.2012.Taxonomic notes on rhacophorid frogs(Rhacophorinae:Rhacophoridae:Anura)of southern part of Annamite Mountains(Truong Son,Vietnam),with description of three new species.Russian Journal of Herpetology,19(1):23-64.

    Peters WCH.1871.über neue reptilien Aus ostafrica und Sarawak(Borneo),vorzüglich Aus der Sammlung des Hrn.Marquis J.Doria zu Genua.Monatsberichte der K?niglichen Preussische Akademie des Wissenschaften zu Berlin,1871:566-581.

    Posada D,Crandall KA.1998.MODELTEST:testing the model of DNA substitution.Bioinformatics,14(9):817-818.

    Poyarkov NA Jr,Vassilieva AB,Orlov NL,Galoyan EA,Tran TAD,Le DTT,Kretova VD,Geissler P.2014.Taxonomy and distribution of narrow-mouth frogs of the genusMicrohylaTschudi,1838 (Anura: Microhylidae)from Vietnam with descriptions of fi ve new species.Russian Journal of Herpetology,21(2):89-148.

    Poyarkov NA Jr,Orlov NL,Moiseeva AV,Pawangkhanant P,Ruangsuwan T,Vassilieva AB,Galoyan EA,Nguyen TT,Gogoleva SS.2015.Sorting out moss frogs:mtDNA data on taxonomic diversity and phylogenetic relationships of the Indochinese species of the genusTheloderma(Anura,Rhacophoridae).Russian Journal of Herpetology,22(4):241-281.

    Poyarkov NA Jr,Duong TV,Orlov NL,Gogoleva SS,Vassilieva AB,Nguyen LT,Nguyen VHD,Nguyen SN,Che J,Mahony S.2017.Molecular,morphological and acoustic assessment of the genusOphryophryne(Anura,Megophryidae)from Langbian Plateau,southern Vietnam,with description of a new species.ZooKeys,672:49-120.

    Ronquist F,Huelsenbeck JP.2003.MRBAYES 3:Bayesian phylogenetic inference under mixed models.Bioinformatics,19(12):1572-1574.

    Rowley JJL,Hoang HD,Le DTT,Dau VQ,Cao TT.2010.A new species ofLeptolalax(Anura:Megophryidae)from Vietnam and further information on Leptolalax tuberosus.Zootaxa,2660:33-45.

    Rowley JJL,Le DTT,Hoang HD,Dau VQ,Cao TT.2011.Two new species ofTheloderma(Anura:Rhacophoridae)from Vietnam.Zootaxa,3098:1-20.

    Rowley JJL,Le DTT,Dau VQ,Hoang HD,Cao TT.2014.A striking new species of phytotelm-breeding tree frog(Anura:Rhacophoridae)from central Vietnam.Zootaxa,3785(1):25-37.

    Rowley JJL,Tran DTA,Le DTT,Dau VQ,Peloso PLV,Nguyen TQ,Hoang HD,Nguyen TT,Ziegler T.2016.Five new,microendemic Asian Leaf-litter Frogs(Leptolalax)from the southern Annamite mountains,Vietnam.Zootaxa,4085(1):63-102.

    Savage JM,Heyer WR.1997.Digital webbing formulae for anurans:a refinement.Herpetological Review,28(3):131.

    Sivongxay N,Davankham M,Phimmachak S,Phoumixay K,Stuart BL.2016.A new small-sizedTheloderma(Anura:Rhacophoridae)from Laos.Zootaxa,4147(4):433-442.

    Smith MA.1924.New tree-frogs from Indo-China and the Malay Peninsula.Proceedings of the Zoological Society of London,1924:225-234.

    Smith MA.1926.Spolia Mentawia:reptiles and amphibians.Annals and Magazine of Natural History,18(9):76-81.

    Smith MA.1930.The reptilia and amphibia of the Malay peninsula.a supplement to G.A.boulenger’s reptilia and batrachia,1912.Bulletin of the Raffles Museum,3:1-149.

    Smith MA.1931.Amphibia and Reptilia(1930).Zoological Record,67:1-40.StatSoft,Inc.2007.STATISTICA(data analysis software system),version 8.0.www.statsoft.com.

    Stuart BL,Heatwole HF.2004.A newPhilautus(Amphibia:Rhacophoridae)from northern Laos.Asiatic Herpetological Research,10:17-21.

    Taylor EH.1920.Philippine Amphibia.The Philippine Journal of Science,16:213-359.

    Taylor EH.1962.The amphibian fauna of Thailand.The University of Kansas Science Bulletin,43:265-599.

    Thompson JD,Gibson TJ,Plewniak F,Jeanmougin F,Higgins DG.1997.The Clustal_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.Nucleic Acids Research,25(24):4876-4882.

    Vences M,Thomas M,Bonett RM,Vieites DR.2005a.Deciphering amphibian diversity through DNA barcoding:chances and challenges.Philosophical Transactions of the Royal Society B,360(1462):1859-1868.Vences M,Thomas M,Van Der Meijden A,Chiari Y,Vieites DR.2005b.Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians.Frontiers in Zoology,2:5.

    Vieites DR,Wollenberg KC,Andreone F,Kohler J,Glaw F,Vences M.2009.Vast underestimation of Madagascar’s biodiversity evidenced by an integrative amphibian inventory.Proceedings of the National Academy of Sciences of the United States of America,106(20):8267-8272.

    Tschudi JJ von.1838.Classification der Batrachier:Mit Berücksichtigung der fossilen Thiere dieser Abtheilung der Reptilien.Neuchatel:Petitpierre.

    Wilkinson JA,Drewes RC.2000.Character assessment,genus level boundaries,and phylogenetic analyses of the family Rhacophoridae:a review and present day status.Francisco,California: Contemporary Herpetology,1-24.

    Wilkinson JA,Drewes RC,Tatum OL.2002.A molecular phylogenetic analysis of the family Rhacophoridae with an emphasis on the Asian and African genera.Molecular Phylogenetics and Evolution,24(2):265-273.

    Wolf S.1936.Revision der untergattungRhacophorus(ausschliesslich der Madagaskar-Formen).Bulletin of the Raffles Museum,12:137-217.

    Yu GH,Rao DQ,Yang JX,Zhang MW.2007.Non-monophyly ofRhacophorus rhodopus,ThelodermaandPhilautus albopunctatusinferred from mitochondrial 16S rRNA gene sequences.Zoological Research,28(4):437-442.

    Yu GH,Rao DQ,Yang JX,Zhang MW.2008.Phylogenetic relationships among Rhacophorinae (Rhacophoridae,Anura,Amphibia),withan emphasis on the Chinese species.Zoological Journal of the Linnean Society,153(4):733-749.

    Yu GH,Rao DQ,Zhang MW,Yang JX.2009.Re-examination of the phylogeny of Rhacophoridae(Anura)based on mitochondrial and nuclear DNA.Molecular Phylogenetics and Evolution,50(3):571-579.

    Zhao EM,Adler K.1993.Herpetology of China.Oxford,Ohio:Society for the Study of Amphibians&Reptiles.

    Appendix I

    List of examined specimens.

    Theloderma laeve:ZMMU A-4569(Bu Gia Map N.P.,Binh Phuoc Province,southern Vietnam;1 male,holotype ofTh.bambusicolum);ZMMU A-4570(Bu Gia Map N.P.,Binh Phuoc Province,southern Vietnam;2 males,paratypes ofTh.bambusicolum);ZMMU A-4573(Bu Gia Map N.P.,Binh Phuoc Province,southern Vietnam;3 females,paratypes ofTh.bambusicolum);ZMMU A-4571(Cat Loc area,Cat Tien N.P.,Lam Dong Province,southern Vietnam;4 males,1 female,paratypes ofTh.bambusicolum);

    Theloderma palliatum:ZMMU A-5835(Chu Yang Sin N.P.,Dak Lak Province,southern Vietnam;4 males);

    Theloderma truongsonense:ZMMU A-5827(SaoLa Reserve,ARoang area,ALuoi district,Thua-Thien-Hue Province,central Vietnam;5 males);

    Theloderma auratumsp.nov.:ZMMU A-5828(Kon Chu Rang N.R.,Gia Lai Province,central Vietnam;1 male,holotype);ZMMU A-5829(Kon Chu Rang N.R.,Gia Lai Province,central Vietnam;1 male,paratype);ZMMU A-5830(Kon Chu Rang N.R.,Gia Lai Province,central Vietnam;1 male,paratype);ZMMU A-5831(Kon Ka Kinh N.P.,Gia Lai Province,central Vietnam;1 male,paratype);ZMMU A-5832(Thac Nham Forest,Kon Plong area,Kon Tum Province,central Vietnam;1 male,paratype);ZMMU A-5833-A-5834(SaoLa Reserve,ARoang area,ALuoi district,Thua-Thien-Hue Province,central Vietnam;2 males).

    最好的美女福利视频网| 热99在线观看视频| 欧美日韩国产亚洲二区| 夜夜夜夜夜久久久久| 免费看美女性在线毛片视频| 国产三级在线视频| 欧美不卡视频在线免费观看| 亚洲不卡免费看| 嫁个100分男人电影在线观看| 十八禁网站免费在线| 欧美一区二区精品小视频在线| a级一级毛片免费在线观看| 国产亚洲欧美98| 免费搜索国产男女视频| 我的女老师完整版在线观看| av在线蜜桃| 在线观看舔阴道视频| 久久草成人影院| 午夜福利在线观看吧| 国产精品国产三级国产av玫瑰| 99视频精品全部免费 在线| 久久人妻av系列| 麻豆国产av国片精品| 中文字幕熟女人妻在线| 精品人妻视频免费看| 亚洲精品一卡2卡三卡4卡5卡| 精品一区二区三区av网在线观看| 中文字幕精品亚洲无线码一区| 一夜夜www| 国产私拍福利视频在线观看| 亚洲av中文av极速乱 | 成人特级黄色片久久久久久久| 国产探花在线观看一区二区| 国产亚洲欧美98| 午夜福利在线观看吧| 成人特级av手机在线观看| 国产视频一区二区在线看| av黄色大香蕉| 日本欧美国产在线视频| 国产国拍精品亚洲av在线观看| 97碰自拍视频| 不卡视频在线观看欧美| 别揉我奶头 嗯啊视频| 禁无遮挡网站| 久久久久久久午夜电影| 国产午夜精品久久久久久一区二区三区 | 又黄又爽又刺激的免费视频.| 99精品久久久久人妻精品| 露出奶头的视频| 一级a爱片免费观看的视频| 欧美不卡视频在线免费观看| 91麻豆av在线| 免费人成视频x8x8入口观看| 精品乱码久久久久久99久播| 欧美一区二区精品小视频在线| 亚洲成av人片在线播放无| 老司机福利观看| 国产亚洲av嫩草精品影院| 国产国拍精品亚洲av在线观看| 亚洲真实伦在线观看| 一个人看的www免费观看视频| 观看免费一级毛片| 国产成人av教育| 一个人免费在线观看电影| 我的女老师完整版在线观看| 国产男靠女视频免费网站| 一级黄色大片毛片| 真人做人爱边吃奶动态| 别揉我奶头 嗯啊视频| 少妇猛男粗大的猛烈进出视频 | 丰满的人妻完整版| 欧美又色又爽又黄视频| 欧美最新免费一区二区三区| 免费av观看视频| 国产精品久久电影中文字幕| 国内精品宾馆在线| 狠狠狠狠99中文字幕| 精品午夜福利在线看| 一级a爱片免费观看的视频| 悠悠久久av| 亚洲色图av天堂| 日韩欧美在线乱码| 久久人妻av系列| 国产精品久久久久久久久免| 麻豆国产av国片精品| 校园人妻丝袜中文字幕| 日韩欧美三级三区| 欧美日韩精品成人综合77777| 成人毛片a级毛片在线播放| 网址你懂的国产日韩在线| 国产精品人妻久久久久久| 国产在线精品亚洲第一网站| 亚洲精品456在线播放app | 国产精品国产高清国产av| 午夜视频国产福利| 乱码一卡2卡4卡精品| 一级av片app| 午夜爱爱视频在线播放| videossex国产| 搡老妇女老女人老熟妇| 有码 亚洲区| 91麻豆精品激情在线观看国产| 亚洲欧美日韩无卡精品| 三级男女做爰猛烈吃奶摸视频| 天堂av国产一区二区熟女人妻| 老师上课跳d突然被开到最大视频| 成人鲁丝片一二三区免费| 99精品久久久久人妻精品| 欧美成人一区二区免费高清观看| 国产精品自产拍在线观看55亚洲| 美女高潮喷水抽搐中文字幕| 最新中文字幕久久久久| 欧美精品啪啪一区二区三区| 联通29元200g的流量卡| 亚洲精品乱码久久久v下载方式| 国产视频内射| 精品人妻熟女av久视频| 亚洲欧美清纯卡通| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久中文| 精品99又大又爽又粗少妇毛片 | 一个人看的www免费观看视频| 欧美日本视频| 日本一本二区三区精品| 内地一区二区视频在线| 日本黄色片子视频| 欧美xxxx性猛交bbbb| 久久九九热精品免费| 欧美成人一区二区免费高清观看| 最近中文字幕高清免费大全6 | 精品欧美国产一区二区三| 色精品久久人妻99蜜桃| 亚洲欧美日韩卡通动漫| 听说在线观看完整版免费高清| 国产av不卡久久| 免费av毛片视频| 国产不卡一卡二| 久久国产乱子免费精品| 日韩精品青青久久久久久| 国产伦一二天堂av在线观看| av专区在线播放| 欧美一级a爱片免费观看看| 久久久久久久久久久丰满 | 欧美日本视频| 赤兔流量卡办理| 99精品在免费线老司机午夜| 淫妇啪啪啪对白视频| 久久精品国产鲁丝片午夜精品 | 人人妻,人人澡人人爽秒播| 欧美三级亚洲精品| 国产视频一区二区在线看| 欧美最新免费一区二区三区| 欧美不卡视频在线免费观看| 舔av片在线| 日本a在线网址| 少妇丰满av| 又粗又爽又猛毛片免费看| 三级国产精品欧美在线观看| 欧美+亚洲+日韩+国产| 又粗又爽又猛毛片免费看| 午夜免费成人在线视频| av在线观看视频网站免费| 午夜福利高清视频| 日韩精品青青久久久久久| 成年女人看的毛片在线观看| 欧美+亚洲+日韩+国产| 精品人妻偷拍中文字幕| 大型黄色视频在线免费观看| 人人妻,人人澡人人爽秒播| 日本五十路高清| 久久这里只有精品中国| 免费av毛片视频| 色哟哟·www| 成年女人永久免费观看视频| 噜噜噜噜噜久久久久久91| 国产精品自产拍在线观看55亚洲| 日韩欧美免费精品| 真人做人爱边吃奶动态| 99热6这里只有精品| 国产亚洲精品久久久com| 高清毛片免费观看视频网站| av在线天堂中文字幕| 日韩一本色道免费dvd| 久久久久国产精品人妻aⅴ院| 天堂网av新在线| 精品99又大又爽又粗少妇毛片 | 免费在线观看成人毛片| 午夜老司机福利剧场| 免费看美女性在线毛片视频| 亚洲狠狠婷婷综合久久图片| 欧美性猛交黑人性爽| 亚洲无线在线观看| 听说在线观看完整版免费高清| 国国产精品蜜臀av免费| 欧美性感艳星| 超碰av人人做人人爽久久| 淫秽高清视频在线观看| 男女边吃奶边做爰视频| 久久久久久久午夜电影| 少妇裸体淫交视频免费看高清| 国产精品精品国产色婷婷| 国产精品爽爽va在线观看网站| 免费大片18禁| а√天堂www在线а√下载| 最近最新免费中文字幕在线| 午夜福利在线观看吧| 国产精品精品国产色婷婷| 在线观看美女被高潮喷水网站| 啪啪无遮挡十八禁网站| 亚洲欧美清纯卡通| 精品人妻熟女av久视频| 桃色一区二区三区在线观看| 欧美xxxx性猛交bbbb| 国产69精品久久久久777片| 99热这里只有精品一区| 久久久久精品国产欧美久久久| 在线播放国产精品三级| 精品久久久久久久久久免费视频| 中文在线观看免费www的网站| 亚洲第一电影网av| 久久人妻av系列| 国产探花在线观看一区二区| 特级一级黄色大片| 日本黄色视频三级网站网址| 国产爱豆传媒在线观看| 久久人人精品亚洲av| 久久久久久九九精品二区国产| 91午夜精品亚洲一区二区三区 | 99riav亚洲国产免费| 久久九九热精品免费| 一进一出抽搐动态| 亚洲男人的天堂狠狠| 欧美一级a爱片免费观看看| 99久久精品一区二区三区| 午夜日韩欧美国产| 国产欧美日韩精品一区二区| 国产精品一区二区免费欧美| 国产av麻豆久久久久久久| 国产精品99久久久久久久久| 在线观看免费视频日本深夜| 悠悠久久av| 日本色播在线视频| 香蕉av资源在线| 久久精品综合一区二区三区| 免费人成视频x8x8入口观看| 午夜免费激情av| 亚洲欧美日韩高清专用| 亚洲成人久久爱视频| 久久精品国产99精品国产亚洲性色| 91精品国产九色| 免费观看精品视频网站| 伦精品一区二区三区| 熟女电影av网| 婷婷精品国产亚洲av在线| 99热网站在线观看| 舔av片在线| 嫩草影院新地址| 一本一本综合久久| 中出人妻视频一区二区| 久久国产精品人妻蜜桃| 亚洲午夜理论影院| 国产精品亚洲一级av第二区| 成人国产一区最新在线观看| 欧美bdsm另类| 精品午夜福利在线看| 我要搜黄色片| 久久久久九九精品影院| 国产成人aa在线观看| 亚洲欧美精品综合久久99| 国内精品久久久久精免费| 网址你懂的国产日韩在线| 亚洲精品一区av在线观看| 波多野结衣高清作品| 99久久无色码亚洲精品果冻| 成人鲁丝片一二三区免费| 黄色欧美视频在线观看| 婷婷精品国产亚洲av| 亚洲精品一区av在线观看| 亚洲,欧美,日韩| 神马国产精品三级电影在线观看| 亚洲av美国av| 国模一区二区三区四区视频| 韩国av在线不卡| 黄色丝袜av网址大全| 国产伦一二天堂av在线观看| 国产精品一区二区性色av| 亚洲综合色惰| 天美传媒精品一区二区| 午夜激情福利司机影院| 淫妇啪啪啪对白视频| 悠悠久久av| 久久人人爽人人爽人人片va| 日本欧美国产在线视频| 午夜精品一区二区三区免费看| 免费观看在线日韩| 日韩欧美一区二区三区在线观看| 午夜福利在线观看吧| 久久国产精品人妻蜜桃| 网址你懂的国产日韩在线| 他把我摸到了高潮在线观看| 欧美xxxx黑人xx丫x性爽| 国产69精品久久久久777片| 在线免费观看的www视频| 亚洲欧美日韩无卡精品| 超碰av人人做人人爽久久| 在线看三级毛片| 亚洲成人免费电影在线观看| 日韩av在线大香蕉| 欧美潮喷喷水| 婷婷六月久久综合丁香| 97超级碰碰碰精品色视频在线观看| 在线观看免费视频日本深夜| 亚洲av熟女| 欧美日本视频| a级毛片免费高清观看在线播放| 九色国产91popny在线| 在线观看一区二区三区| 99热6这里只有精品| 我的老师免费观看完整版| 久久天躁狠狠躁夜夜2o2o| av在线亚洲专区| 免费电影在线观看免费观看| 色尼玛亚洲综合影院| 少妇人妻一区二区三区视频| 桃色一区二区三区在线观看| 欧美又色又爽又黄视频| x7x7x7水蜜桃| 欧美又色又爽又黄视频| 国内精品久久久久久久电影| 欧美+亚洲+日韩+国产| 成年人黄色毛片网站| 一级黄色大片毛片| 中国美女看黄片| 欧美精品啪啪一区二区三区| 成人午夜高清在线视频| 国产免费男女视频| 一级毛片久久久久久久久女| 国产精品亚洲美女久久久| 成人特级av手机在线观看| 人妻丰满熟妇av一区二区三区| 久久久久久伊人网av| 国产v大片淫在线免费观看| 日日啪夜夜撸| 国产中年淑女户外野战色| 1000部很黄的大片| 在线免费十八禁| 女同久久另类99精品国产91| 日韩欧美精品免费久久| 九九在线视频观看精品| 嫁个100分男人电影在线观看| 国产精品1区2区在线观看.| 成人午夜高清在线视频| 看十八女毛片水多多多| 亚洲精华国产精华精| 天天一区二区日本电影三级| 69人妻影院| 久久久久久久久久成人| 搡老妇女老女人老熟妇| 人妻制服诱惑在线中文字幕| 狠狠狠狠99中文字幕| 欧美+日韩+精品| 99九九线精品视频在线观看视频| 亚洲成av人片在线播放无| 精品久久国产蜜桃| 小蜜桃在线观看免费完整版高清| 男女边吃奶边做爰视频| 免费一级毛片在线播放高清视频| 久久精品影院6| 淫妇啪啪啪对白视频| 久久中文看片网| av黄色大香蕉| 久久久久久久久久黄片| 男女那种视频在线观看| 波多野结衣高清作品| 亚洲人成网站高清观看| 天堂动漫精品| 窝窝影院91人妻| 有码 亚洲区| 99热精品在线国产| 色综合站精品国产| 搡老熟女国产l中国老女人| 极品教师在线免费播放| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 国内精品一区二区在线观看| 欧美在线一区亚洲| 小说图片视频综合网站| 午夜精品在线福利| 国产精品一区二区性色av| 亚洲av免费高清在线观看| 精品一区二区三区视频在线观看免费| 国产淫片久久久久久久久| 1000部很黄的大片| 99热6这里只有精品| av天堂在线播放| 日本撒尿小便嘘嘘汇集6| 99国产极品粉嫩在线观看| 日本撒尿小便嘘嘘汇集6| 欧美日本亚洲视频在线播放| 99久国产av精品| www日本黄色视频网| 婷婷丁香在线五月| 久久久成人免费电影| 99久久成人亚洲精品观看| 亚洲综合色惰| 亚洲人成网站高清观看| 女生性感内裤真人,穿戴方法视频| 免费一级毛片在线播放高清视频| 亚洲性夜色夜夜综合| 中国美白少妇内射xxxbb| 欧美日韩乱码在线| 免费无遮挡裸体视频| 午夜福利在线观看吧| 九九在线视频观看精品| 久9热在线精品视频| 成人毛片a级毛片在线播放| 精品人妻熟女av久视频| 国产av在哪里看| 欧美丝袜亚洲另类 | 午夜精品久久久久久毛片777| 精品一区二区三区视频在线观看免费| 乱系列少妇在线播放| 亚洲av中文av极速乱 | 日韩精品有码人妻一区| 九色国产91popny在线| 免费高清视频大片| 色综合婷婷激情| 一级毛片久久久久久久久女| 欧美激情久久久久久爽电影| 国产伦人伦偷精品视频| 日韩大尺度精品在线看网址| 一区福利在线观看| 桃色一区二区三区在线观看| 国内精品美女久久久久久| 丰满乱子伦码专区| 免费大片18禁| 亚洲精品456在线播放app | 免费av毛片视频| 亚洲久久久久久中文字幕| 日韩国内少妇激情av| 国产黄a三级三级三级人| 天堂影院成人在线观看| 在线天堂最新版资源| 成人午夜高清在线视频| 赤兔流量卡办理| 五月玫瑰六月丁香| 蜜桃久久精品国产亚洲av| 永久网站在线| 超碰av人人做人人爽久久| 精品久久久噜噜| 男女下面进入的视频免费午夜| 色视频www国产| 久久久久久九九精品二区国产| 色播亚洲综合网| 欧美+日韩+精品| 国产精品自产拍在线观看55亚洲| 国产高清三级在线| 草草在线视频免费看| 久99久视频精品免费| 不卡一级毛片| 精品久久久久久久久久久久久| 他把我摸到了高潮在线观看| av专区在线播放| 琪琪午夜伦伦电影理论片6080| 国产伦一二天堂av在线观看| 欧美成人一区二区免费高清观看| 国内精品美女久久久久久| aaaaa片日本免费| 午夜影院日韩av| 久久精品国产亚洲av天美| 亚洲熟妇中文字幕五十中出| 中国美白少妇内射xxxbb| 亚洲国产精品成人综合色| 夜夜看夜夜爽夜夜摸| 熟女人妻精品中文字幕| 亚洲成人中文字幕在线播放| 老女人水多毛片| 国产黄a三级三级三级人| 久久久久久久久久成人| 久久久久九九精品影院| 波多野结衣高清无吗| 老司机福利观看| 国产精品,欧美在线| 午夜视频国产福利| 女人十人毛片免费观看3o分钟| 全区人妻精品视频| 麻豆精品久久久久久蜜桃| 成人国产麻豆网| 国产在线精品亚洲第一网站| 99九九线精品视频在线观看视频| av天堂中文字幕网| 成人性生交大片免费视频hd| 一级a爱片免费观看的视频| 午夜久久久久精精品| 国产一区二区亚洲精品在线观看| 亚洲va日本ⅴa欧美va伊人久久| 99riav亚洲国产免费| 91午夜精品亚洲一区二区三区 | 久久久久久伊人网av| 真人做人爱边吃奶动态| 男女之事视频高清在线观看| 免费观看在线日韩| 国产69精品久久久久777片| 久久久久久伊人网av| 成人精品一区二区免费| 色播亚洲综合网| 12—13女人毛片做爰片一| 两人在一起打扑克的视频| 婷婷六月久久综合丁香| 男女那种视频在线观看| 十八禁网站免费在线| 国内精品久久久久精免费| 亚洲乱码一区二区免费版| 91麻豆精品激情在线观看国产| 欧美国产日韩亚洲一区| 国产精品一及| 久久草成人影院| 一a级毛片在线观看| 精品久久久久久成人av| 神马国产精品三级电影在线观看| 国内精品久久久久久久电影| 97碰自拍视频| 无人区码免费观看不卡| 91久久精品电影网| 成人av在线播放网站| 精品一区二区三区视频在线观看免费| 色播亚洲综合网| 人妻制服诱惑在线中文字幕| 男女之事视频高清在线观看| 搞女人的毛片| 制服丝袜大香蕉在线| 日日摸夜夜添夜夜添av毛片 | 中文字幕免费在线视频6| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 人人妻人人澡欧美一区二区| 国内精品久久久久久久电影| 日韩欧美国产一区二区入口| 国产中年淑女户外野战色| 嫩草影视91久久| 国产成人av教育| 波多野结衣高清无吗| 97超级碰碰碰精品色视频在线观看| 在线a可以看的网站| 日韩欧美精品免费久久| 国产午夜福利久久久久久| 国产淫片久久久久久久久| 嫩草影视91久久| 日日干狠狠操夜夜爽| 蜜桃久久精品国产亚洲av| 男人舔女人下体高潮全视频| 国产精品永久免费网站| 美女高潮喷水抽搐中文字幕| 日韩大尺度精品在线看网址| 国产av不卡久久| 男女那种视频在线观看| 精品午夜福利在线看| 中文字幕久久专区| 美女黄网站色视频| 国产av一区在线观看免费| 一区福利在线观看| 老司机福利观看| 热99在线观看视频| 欧美日韩亚洲国产一区二区在线观看| av国产免费在线观看| 蜜桃亚洲精品一区二区三区| 国产不卡一卡二| 天堂影院成人在线观看| 亚洲av成人av| 小蜜桃在线观看免费完整版高清| 日韩欧美免费精品| 少妇熟女aⅴ在线视频| 亚洲精华国产精华精| 亚洲第一区二区三区不卡| 国产麻豆成人av免费视频| 午夜影院日韩av| 一卡2卡三卡四卡精品乱码亚洲| 久久这里只有精品中国| 国产在线精品亚洲第一网站| 亚洲国产欧美人成| 欧美日本视频| 免费人成视频x8x8入口观看| 熟女电影av网| 如何舔出高潮| 九九热线精品视视频播放| av天堂中文字幕网| 免费看av在线观看网站| 变态另类丝袜制服| h日本视频在线播放| 国产日本99.免费观看| 天堂动漫精品| 亚洲经典国产精华液单| 国产视频内射| 日本精品一区二区三区蜜桃| 欧美不卡视频在线免费观看| 无人区码免费观看不卡| 亚洲av美国av| 看十八女毛片水多多多| 成年免费大片在线观看| 国产av不卡久久| 在线a可以看的网站| 少妇人妻精品综合一区二区 | 真人做人爱边吃奶动态| 美女黄网站色视频| 国内精品久久久久精免费| 国产精品日韩av在线免费观看| 中文字幕高清在线视频| 亚洲性夜色夜夜综合| 亚洲欧美日韩东京热| 最后的刺客免费高清国语| 男女视频在线观看网站免费| 麻豆久久精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 18禁黄网站禁片免费观看直播| 黄色丝袜av网址大全|