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

    Thermal Biology of Cold-climate Distributed Heilongjiang Grass Lizard,Takydromus amurensis

    2020-12-30 06:59:26XinHAOShiangTAOYuMENGJingyangLIULuoxinCUIWanliLIUBaojunSUNPengLIUandWengeZHAO
    Asian Herpetological Research 2020年4期

    Xin HAO,Shiang TAO,Yu MENG,Jingyang LIU,Luoxin CUI,Wanli LIU,Baojun SUN,Peng LIU*,# and Wenge ZHAO*,#

    1 College of Life Science and Technology,Harbin Normal University,Harbin 150025,Heilongjiang,China

    2 Key Laboratory of Animal Ecology and Conservation Biology,Institute of Zoology,Chinese Academy of Sciences,Beijing 100101,China

    3 College of Chemistry and Life sciences,Zhejiang Normal University,Jinhua 321004, Zhejiang,China

    4 School of Chemical Engineering,University of Science and Technology Liaoning,Anshan 114051,Liaoning,China

    Abstract Thermal biology traits reflect thermal adaptations to an environment and can be used to infer responses to climate warming in animal species.Within a widespread genus or species,assessing the latitudinal or altitudinal gradient of thermal physiological traits is essential to reveal thermal adaptations and determine future vulnerability to climate warming geographically.We determined the thermal biology traits of a cold-climate distributed liza rd,Takydromus amurensis,a nd integra ted published thermal biology traits within the genus Takydromus to reveal a preliminary geographical pattern in thermal adaptation.The mean selected body temperature (cloaca temperature; Tsel),critical thermal maximum (CTmax),critical thermal minimum(CTmin),and optimal temperature for locomotion (i.e.,sprint speed; Topt) of T.amurensis were 32.6,45.1,3.1,and 33.4 °C,respectively.The resting metabolic rates of T.amurensis were positively related to temperature from 18 °C to 38 °C.We compared the traits of tropical T.sexlineatus,subtropical T.septentrionalis,and T.wolteri with T.amurensis and found that the CTmax and thermal tolerance range (the difference between CTmax and CTmin; TTR) increased toward high latitudes,whereas CTmin increased toward low latitudes in these four Takydromus lizards.According to this preliminary pattern,we speculate the species at medium and low latitudes would be more vulnerable to extreme heat events caused by ongoing climate warming.We highlight the importance of integrating thermal biology traits along geographical clues,and its potential contribution to evaluate the vulnerabilities of species in the context of climate warming.

    Keywords counter gradient,CTmax,CTmin,thermal biological trait,thermal tolerance range,Tsel,Takydromus

    1.Introduction

    Ongoing climate warming has imposed novel stressors on animals (e.g.,Thomaset al.,2004),which has affected their populations in a number of ways.For example,many species have altered their distributions,and some populations face risk of collapse (e.g.,Poundset al.,1999;Wilsonet al.,2005).Recent studies have shown that the thermal safety margin (difference between environmental temperature and critical maximum temperature,CTmax) of low-latitude animals is narrower than that of animals living in mid and high latitudes,and the thermal biological characteristics of low-latitude animals are more susceptible to temperature variation,indicating that low-latitude animals are more vulnerable to climate warming (e.g.,Dillonet al.,2010; Güsewellet al.,2017;Sundayet al.,2014).Thus,a detailed understanding of the thermal biology traits and the thermal responses of animals is not only useful for studying adaptations to thermal environments,but also necessary for determining the vulnerability to thermal variation caused by climate warming.

    Thermal biology traits have been shown to be important proxies of the thermal adaptation of animals to their environments,as measured by selected body temperatures (Tsel),CTmax,CTmin,and the thermal dependency of performance such as locomotion (e.g.,Huey and Hertz,1984; Huey and Kingsolver,1989; Angillettaet al.,2002a).Unlike endotherms,which can maintain a stable body temperature in a wide temperature range at great energetic costs (Bennett and Ruben,1979),the body temperatures of reptiles and other ectotherms are regulated by thermal interactions with the environment and thus are more susceptible to thermal variations and their effects upon behavior and physiological processes(Huey and Stevenson,1979; Angillettaet al.,2002a,2002b).Accordingly,reptiles have a high risk of extinction owing to climate warming (e.g.,Sinervoet al.,2010).In addition,they are affected by internal and external factors (Wanget al.,2017).Through the analysis of the thermal biology traits of widespread genera or species,inferences can be made regarding the geographical patterns and interspecific variation of traits.This can aid in determining the vulnerabilities of widespread genera to climate change along geographical gradients (Hueyet al.,2009; Sinervoet al.,2010).

    Reptiles can maintain their body temperature within a moderate range by thermoregulation,with a range of optimal body temperatures for different processes(Van Dammeet al.,1991; Bauwenset al.,1995).The selected body temperatures of a reptile represent the body temperatures without biotic or abiotic constraits(Hueyet al.,1979; Hertzet al.,1993),which can be determined by the temperature gradient in the laboratory,although extrapolation of laboratory data to field populations should be interpreted cautiously.Because Tselrepresents the “ideal” range of body temperatures ectothermic animals aim to maintain through behavioral thermoregulation,it is normally different from active body temperatures (Lichtet al.,1966; Avery,1978).CTmaxand CTminare the maximum and minimum temperatures at which vitality can be sustained,respectively (Hertzet al.,1983; Gilbert and Miles,2017).If the body temperature is lower than CTminor higher than CTmax,reptiles would lose motor function and suffer cold narcosis (temperature lower than CTmin) or muscle spasms (temperature higher than CTmax) (Lutterschmidt and Hutchison,1997;Angilletta,2009; Camacho and Rusch,2017).Therefore,CTmaxand CTminare used to determine the viability of reptiles in various thermal environments (Hertzet al.,1983; Deutschet al.,2008).

    Metabolism is the primary physiological process that determines an organism’s demand from the environment and the allocation of energy among body functions(Brownet al.,2004).In ectotherms,the resting metabolic rate is often used as a proxy for metabolism and is defined as the energy expenditure of animals while in the resting,fasting,and nonproductive states (Andrews and Pough,1985; Christianet al.,1999; Berget al.,2017).In reptiles,resting metabolic rates are highly sensitive to body temperature,requiring the organism to respond quickly to variations in the thermal environment (Maet al.,2018a,2018b; Sunet al.,2018).As another critical function,locomotion is also significantly sensitive to body temperature and is related to fitness because of its importance for foraging,avoiding predators,and choosing mating partners (Bauwenset al.,1995; Chenet al.,2003; Sunet al.,2018).In summary,evaluating the thermal sensitivity of critical indices (i.e.,metabolic rates and locomotion) is essential to understand thermal adaptation in ectotherms (Angilletta,2009).

    Takydromusliza rds have a la rge geographical distribution,with a significant latitudinal span,ranging from cold-temperate to tropical regions.As a species that inhabits cold climates,the Heilongjiang grass lizard,Takydromus amurensis,is a small lacertid lizard (up to 70 mm snout-vent length [SVL]) which is distributed in northeast China and the boundary of China,Russia,and the Korean Peninsula (Xuet al.,2017).Like otherTakydromuslizards in mainland China,T.amurensisinhabits areas at the boundary of forests and grasslands or scrub habitats.Females are able to lay one or two clutches from May to July,with an average clutch size of six eggs (Zhaoet al.,1999).T.amurensisis distributed in high latitudes,which typically exhibit low average temperatures and drastic thermal fluctuations (Figure 1).Thermal biology traits should be investigated inT.amurensisto obtain a comprehensive understanding of thermal adaptation in cold-climate ectotherms.The thermal biology of several other congener species has been previously studied,includingT.wolteri,T.septentrionalis,andT.sexlineatus,which inhabit subtropical and tropical areas,respectively (Jiet al.,1996a; Chenet al.,2003; Zhang and Ji,2004).By assessing the same traits inT.amurensis,i.e.,those investigated in these previous studies,we can improve our understanding of the specific variations and preliminary patterns of thermal biology traits inTakydromusspecies along latitudinal and altitudinal gradients.

    Figure 1 Annual temperature environments for Anshan population of T.amurensis,where we collected the lizards.The red line indicates the average air temperature along the year,and grey area indicates the daily fluctuations of air temperature.

    We determined selected body temperatures,thermal tolerance (i.e.,CTmaxand CTmin),and thermally sensitive traits including resting metabolic rate and sprint speed as focal study traits for the Heilongjiang grass lizard,T.amurensis.We aimed to (1) contribute with new data on a cold-climate distributedTakydromusspecies (i.e.,T.amurensis) for the first time; (2) compare thermal biology traits amongTakydromuslizards with published data concerningT.septentrionalis,T.wolteri,andT.sexlineatus;and (3)summarize preliminary patterns of variation in Tsel,CTmax,CTmin,and optimal body temperature for performance inTakydromuslizards along a latitudinal or altitudinal gradient.

    2.Materials and Methods

    2.1.Study species collectionWe collected adultT.amurensis(n=13; 4 males and 9 females) in late April 2016,using either a noose or by hand,from Anshan (see details in Figure 1 for annual thermal environments),Liaoning,China (41°01’ N,120°7’ E;~210 m).After collection,we transported the lizards back to the laboratory in Beijing.The lizards were individually housed in plastic terraria(350 mm × 250 mm × 220 mm,length × width × height).The terraria were set up with a mixture of moist soil and sand,with randomly placed patches of grass.The lizards were housed in a temperature-controlled room at 18°C under a natural photoperiod,with a supplementary heating lamp suspended above one end of the terrarium from 06:00 to 20:00.During the hea ting period,temperatures in the terraria ranged from 18 °C to 40 °C.Food (crickets and larvalTenebrio molitor) and water were providedad libitum.After three days of captivity,thermal biology traits were started to be determined.

    2.2.Selected body temperatures (Tsel)Selected body temperatures (Tsel) were determined according to the established protocol,with minor modifications (Shuet al.,2010).In brief,Tselwas measured in a custom-made terrarium (1000 mm × 500 mm × 300 mm,length × width× height) placed in a temperature-controlled room at 18 °C.The terrarium was fitted with grass substrate and‘natural’ retreats to mimic the natural environment.One heating light bulb (275 W) was suspended above one end of the terrarium,creating a thermal gradient from 18 °C to 60 °C during the heating period.The heating period lasted 14 h daily from 06:00 to 20:00.On the first day,a group of lizards (four to five) were introduced into the terrarium from the cool end at 16:00 for acclimation to the surroundings.Then,on the second day,the body temperatures of the lizards were measured twice at 09:00 and 15:00,with an UNT-325 electronic thermometers(UNT T-325,Shanghai,China) by inserting the tip of the probe into cloacas.During the measurements,each lizard was captured by hand and its body temperature was measured immediately (within 30 s),without disturbing the other lizards in the terrarium.Then,the average of the two body temperature records of each lizard was estimated and used as the lizard’s selected body temperature.

    2.3.Thermal toleranceThe critical thermal minimum(CTmin) a nd critical thermal maximum (CTmax)temperatures of the lizards were determined in a programmed incubator (KB 240,Binder,Germany).We cooled (for CTmindetermination) or heated (for CTmaxdetermination) lizards from a starting temperature of 28 °C at a rate of 1 °C per 10 min.During the cooling and heating,we continuously monitored the behavior of the lizards.When the lizards lost the capacity to respond to intense stimulation and could not right themselves after being turned over,the body temperatures were recorded as CTminor CTmax,respectively (Zhang and Ji,2004; Xu and Ji,2006).The lizards were then moved to a 28 °C chamber to recover.The protocol stands that if a lizard cannot recover from the cold or heat shock,its record is eliminated for further analysis (Liet al.,2017).However,in this study all the lizards were recovered from the shock

    2.4.Locomotor performanceLocomotor performance was estimated using sprint speed.The sprint speed was measured at five test temperatures ranging from 18 °C to 38 °C (18 °C,23 °C,28 °C,33 °C,and 38 °C),based on the reported range of field body temperatures of the species(Xuet al.,2017),in a randomized sequence.Locomotor performance was measured once per day under one temperature treatment (18 °C,33 °C,28 °C,23 °C and 38 °C),in a sequence order.Before the test,the lizards were placed in an incubator at the test temperature for approximately 2 h for acclimation.To ensure that body temperatures of lizards matched those of the respective test treatment,we measured the body temperatures of a subset of lizards before starting stimulating them to run(Sunet al.,2014).Locomotor performance was tested by stimulating the lizard to run through a racetrack,which was recorded by an HD video camera (Sony,DCRSR220E,Japan).The racetrack was 1500 mm × 100 mm× 150 mm,with intervals marked every 250 mm.Each lizard was stimulated by a paintbrush to run twice at each temperature with an interval of 1 h for rest.The videos were analyzed by Windows Movie Maker.For each lizard,the fastest speed through 250 mm for each time was recorded,and the average of the two fastest records was used as the sprint speed (Sunet al.,2014).

    2.5.Resting metabolic rate (RMR)The resting metabolic rates were estimated by respiratory gas exchange rates measured at five random test temperatures (18 °C,23 °C,28°C,33 °C,and 38 °C).Before the test,the lizards were fasted for at least 12 h.All lizards were acclimated at the correspondent test temperature for 2 h in an incubator (KB 240,Binder,Germany).Then,the lizard was enclosed within the respirometry chamber placed in the incubator (KB 240,Binder,Germany).The respiratory gas exchange was determined using a closedflow respirometry system with a volume of 281.4 mL(Stable System International Inc.Las Vegas,NV,USA),and the resting metabolic rate was estimated via the CO2production rate using a previously established method(Sunet al.,2018).For stabilizing the gas composition of the system (i.e.,baseline),we opened the system to the air for approximately 5 min before changing the system to closed-circuit respirometry.Then,the carbon dioxide production rates (i.e.,VCO2) in the closed circuit were continuously recorded for approximately 10 min.To minimize the effects of circadian rhythms,measurements were conducted from 10:00 to 18:00.The metabolic rates were calculated as the CO2production per gram of body mass per hour (mL/g/ h),following the equation metabolic rate=VCO2× volume/body mass,whereVCO2is the CO2production rate in percentage (%/h) in the closed circuit with a volume of 281.4 mL.

    2.6.Statistical analysisThe normality of distributions and homogeneity of variance were tested using the Kolmogorov-Smirnov test and Bartlett’s test prior to analysis.First,we tested for sex differences in Tsel,CTmax,and CTminusing one-way ANOVA,with sex included as a factor.As there were no significant differences (allP> 0.05),we pooled the data of both sexes for all these traits.For the analysis of locomotion,we performed the Gaussian (R2=0.473),modified Gaussian (R2=0.523),Lorentzian (R2=0.513),Weibull (R2=0.515),and Pseudo-Voigt regressions (R2=0.512) to analyze the sprint speed to body temperature,respectively.With the best fitting,the modified Gaussian regression was employed to analyze the sprint speed,with the following equation:Sprint speed=y0+ae ^ {-0.5 [(Tb-b)/c]2}; whereTbis body temperature; a,b,c andy0are parameters; and e is the natural constant.This was used to calculate the optimal temperature for sprint speed.Repeated measures ANOVAs were used to analyze the effect of test temperature and sex on resting metabolic rate and sprint speed.

    Data of previously illustrated traits ofT.septentrionalisfrom Zhoushan,Zhejiang (29°32’-31°04’N,121°30’-123°25’E;~200 m),T.wolterifrom Chuzhou,Anhui(32°15’-32°21’ N,118°07’-118°18’ E;~260 m),andT.sexlineatusfrom Shaoguan,Guangdong (24°50’ N,113°30’ E;~290 m) (Jiet al.,1996a; Chenet al.,2003;Zhang and Ji,2004; unpublished data) were collected from published figures with Plot Digitizer (http://plotdigitizer.sourceforge.net/).The thermal traits ofT.septentrionalis(collected throughout May and July),T.wolteri(collected in May),andT.sexlineatus(collected in mid-April) were obtained from the literature.Only data collected during the reproductive season were selected to ensure comparability with the results of the current study(Jiet al.,1996a; Chenet al.,2003; Zhang and Ji,2004).General linear models were fitted to assess the latitudinal tendency ofTsel,CTmax,CTmin,and Toptfor sprint speed ofTakydromuslizards.As the four species inhabit similar altitudes (i.e.,from~200 m to~290 m),we did not include this variable in the model.The fact that elevation is fixed in the design of the study provides clarity to be able to detect the effect of latitude.

    3.Results

    The Tsel,CTmin,and CTmaxforT.amurensiswere 32.55 ±0.46 °C,3.06 ± 0.13 ℃,and 45.10 ± 0.11℃,respectively(see details in Table 1).The sprint speed was significantly affected by body temperature,increasing with body temperature from 18 °C to around 33 °C and then decreasing at higher temperatures (F4,48=5.495,P<0.0001; Figure 2a).The optimal body temperature forT.amurensisthat maximized sprint speed was 33.20 °C according to the modified Gaussian regression.Resting metabolic rate was significantly positively related to temperature from 18 °C to 38 °C (F4,48=40.708,P<0.0001; Figure 2b).The Sex had no effect on the resting metabolic rate or sprint speed (bothP> 0.05).

    In the fourTakydromusspecies,Tselwere similar,as indicated by the non-significant tendency along latitudes(y=32.99-0.55x,R2=0.37,P=0.38; Figure 3a).When we integrated the published thermal tolerance data from otherTakydromuslizards,we observed that the CTmaxvalues (y=45.43-0.93x,R2=0.78,P=0.11) remained similar; however,CTmindecreased as latitude increased(y=2.137+1.04x,R2=0.95,P=0.02) (Figure 3b),indicating that the thermal tolerance range increased with latitude.

    Figure 2 Locomotor performance (sprint speed) (a) and resting metabolic rate (RMR) (b) of T.amurensis at different temperatures.RMR was expressed as CO2 production per gram body mass per hour (mL/g/h).Data are expressed as mean ± SEM.

    Table 1 Snout-vent length,body mass,selected body temperatures,the critical thermal maximum,critical thermal minimum and thermal tolerance range (TTR) of Takydromus amurensis.The sample size for all traits are equally 13.

    4.Discussion

    In the context of climate change,understanding the thermal biology traits of animals is fundamental to determining organismal vulnerabilities (e.g.,Sundayet al.,2011,2012,2014).In the current study,the thermal biology traits ofT.amurensiswere quantified at the beginning of the reproductive season,to facilitate a comparison with available data of congeners collected using similar methods.The Tsel,CTmin,CTmax,and sprint speed among different species of the genusTakydromuswere shown to be compa ra ble a nd ecologically meaningful (Jiet al.,1995,1996a; Duet al.,2000; Zhang and Ji,2004).Since the four populations live at similar latitudes (i.e.,ranged from~200 m to~290 m),we only tested the effects of latitude.

    Figure 3 Tsel (a) and CTmax and CTmin (b) of four Takydromus lizards.Data are expressed as mean ± SEM.Species are listed from high to low latitude in X axis.Data are collected by Plot Digitizer(http://plotdigitizer.sourceforge.net/),from the figures and available data of references (Ji et al.,1996a; Xu and Ji,2006; Zhang and Ji,2004).

    Tselis the body temperature at which biological processes function at an optimal/suboptimal level (Van Dammeet al.,1991; Hertzet al.,1993; Blouin-Demerset al.,2000; Angillettaet al.,2002a; Ortega and Martín-Vallejo,2019).InT.amurensis,the optimal temperature for sprint speed was similar to Tsel(33.20 °Cvs.32.55°C),indicating thatT.amurensiscan accurately regulate its body temperature to perform optimal/suboptimal functions.However,as we did not measure operative tempera ture (Te) a nd active body tempera tures(Ta),we were unable to determine the efficiency of thermoregulation.Numerous reptiles have been reported to have temporal,spatial,and individual variation in Tsel(Gatten Jr,1974; Andrews,1998; Stellatelliet al.,2018;Refsnideret al.,2019).There is also increasing evidence demonstrating that thermal biology traits vary among seasons (e.g.,Yanget al.,2008; Sunet al.,2014; see details in below).However,Tsel is normally only determined during the reproductive season and assumed representing the Tsel of the given species or population (Jiet al.,1995,1996a; Shine and Madsen,1996; ?orovi? and Crnobrnja-Isailovi?,2018; Refsnideret al.,2019).It should be noted that althoughT.amurensisis a cold-climate species,the TselinT.amurensisis very similar to published data forT.wolteri,T.septentrionalis,andT.sexlineatus(Figure 3a).Additionally,the lizardsZootoca vivipara,Lacerta agilis,andIberolacerta bonnaliwhich also inhabit cold climates,use thermoregulation to maintain their body temperature warm in responding to cold climates (Herczeget al.,2003;Yanget al.,2015;Ortega et al.,2016).Influence of the thermal environment and activity patterns could induce variations in Tsel(Jiet al.,1996a,1996b; Duet al.,2000;Zhang and Ji,2004).For example,lizards in forested habitats are likely to have lower Tselthan those using more open habitats.The TselofSphenomorphus indicus,whose habitats are forested,was as low as 25.7℃ (Jiet al.,1996b).On the contrary,the lizards using more open habitats and employing basking behaviors would exhibit higher selected body temperatures (e.g.,Liet al.,2017;Wanget al.,2019).T.amurensisis usually active in open microhabitat,with shuttling between open and shaded patches for thermoregulation (Zhao,2002),which may plausibly explain high TselofT.amurensis.Admittedly,we conducted selected body temperatures measurements ofT.amurensiswith established method according to those used in otherTakydromusspecies (e.g.,Jiet al.,1995;Zhang and Ji,2004),increasing literatures indicate that the selected body temperatures of ectotherms should be determined with multiple repeated measurenets (i.e.,eight to ten measurements per individual) (e.g.,Hertzet al.,1993; Liet al.,2017).We also encourage further researches should employ multiple measurements to make a better understanding of thermal preferences in lizards,includingTakydromulizards.

    The thermal tolerance range increased toward high latitudes inTakydromuslizards,with relative conservative CTmaxand decreased CTmintowards high latitudes.(Figure 3b).Microhabitat usage can significantly influence CTmaxand CTmin.For example,CTmaxwas lower inS.indicus,which uses shaded habitats likely to be cooler (37.6 ℃;Jiet al.,1996b) than species occupying open habitats,such asEumeces elegans(41.9 ℃; Duet al.,2000) andT.septentrionalis(42.3 ℃; Jiet al.,1996a).Among sympatric lizards,Eremias argus(45.8 ℃) andE.multiocellta(45.1 °C)that occupying filtered microhabitats,had lower CTmaxthanPhrynocephalus przewalskii(47.1 °C),which uses open habitats (Liet al.,2017; Wanget al.,2019).Analogously,CTminwas noticeably greater in lizards that use open habitats.The CTminofE.chinensis(6.3 °C; Jiet al.,1995) andE.elegans(9.3 ℃; Duet al.,2000),species occupying open habitats,were greater than species using more spatially complex habitats,such asE.argus(1.0 ℃; Luoet al.,2006).Although all the four species compared in the current study use open habitats (Zhaoet al.,1999),T.amurensismay spend more time on open habitats in response to lower temperature conditions,as happens in other highlatitude distributed lizards (Sundayet al.,2011; Sundayet al.,2012).This activity pattern may result in increased time invested in thermoregulation (i.e.,basking) and thus result in higher CTmax.Based on this assumption,global meta-analyses have demonstrated that CTmaxincreases with latitude in lizards (Sundayet al.,2011; Sundayet al.,2012).In contrast,CTminis lower at high latitudes,likely because species inhabiting lower latitudinal or warmer climatic regions do not face the challenges of cold and seasonally varying temperatures (Janzen,1967;Ghalamboret al.,2006).Although we studied the thermal biology traits ofTakydromuslizards across a wide latitude range,our results are preliminary and can be improved upon by incorporating more species and populations.

    The resting metabolic rates ofT.amurensiswere thermally dependent,with allometric enhancement as the temperature increased from 18 °C to 38 °C (Figure 2b).As the basic biological rate,the resting metabolic rates of lizards vary among species,populations,and individuals.The resting metabolic rates ofT.amurensiswere lower than those of the desert-dwelling lizardsE.argus,E.multiocellata,andP.przewalskiiat even temperatures,although they are distributed at similar latitudes (Liet al.,2017; Wanget al.,2019).Unfortunately,the resting metabolic rate ofT.septentrionalis,T.wolteri,andT.sexlineatushas not been studied extensively; this means we were unable to make comparisons amongTakydromuslizards.The sprint speeds ofT.amurensisand otherTakydromusspecies were enhanced as the temperatures increased until the optimal temperature was reached,and then decreased as the temperature increased beyond the optimal temperature (Jiet al.,1996a; Zhang and Ji,2004).However,in all four species we compared,T.amurensishad the lowest sprint speed (Figure 4a).Notably,the optimal body temperature for the sprint speed ofT.amurensiswas 33.20 °C.In the other three species,the optimal temperatures for sprint speed increased as the latitude decreased (Figure 4b).For species from warm and stable thermal environments,sprint speed is more thermally dependent,and the optimal temperatures are higher,as shown inT.sexlineatus(Figure 4; Zhang and Ji,2004).WhyT.amurensis,a cold-climate species,has such a high optimal temperature for sprint speed is still unclear.However,the optimal temperature for sprint speed (i.e.,33.20 ℃),which is roughly equivalent to Tsel(i.e.,32.55℃),imposed an optimal/suboptimal sprint speed at TselinT.amurensis.Future studies could investigate the effect of temperature at a finer scale (e.g.,2 °C between each test temperature),to further clarify the thermal dependence of sprint speed inT.amurensisand thus determine the optimal temperature more precisely.Although the mechanism is unclear,Tseland Toptmight be subject to both thermal environment and confounding evolutionary factors (Huey and Kingsolver,1989).Alternatively,maintaining Tbat a level that appears suboptimal may allow the animal to remain responsive to environmental variation and,if required,increase activity such as those associated with foraging and escaping,following the‘suboptimal is optimal’ hypothesis (Martin and Huey,2008).For example,the range of Teand Tbfor some Liolaemidae lizards is below both their preferred body temperature and the optimal temperature at which they can reach maximum locomotor performance (Boninoet al.,2011).Addtionally,Tbat suboptimal temperatures is also beneficial to avoid over-heating effects in performance,since the performance would drop sharply if the optimal temperatures are surpassed (e.g.,Martin and Huey,2008).

    Based on our current preliminary pattern of Tsel,CTmin,CTmax,and Toptfor sprint speed inTakydromuslizards,it is plausible thatT.amurensiswould be less vulnerable to the warming temperatures imposed by climate change than otherTakydromusspecies in habitat in low latitudes.First,the CTmaxforT.amurensisis high; thus,this may indicate that the thermal safety margin is large.However,we were unable to verify this as we did not calculate the thermal safety margin owing to the lack of precise data on environmental temperatures (Sundayet al.,2014).Second,asT.amurensisinhabit areas with cooler temperatures than their physiological optima,it is possible that warming temperatures may enhance their fitness,as in high-latitude species (e.g.,Cabezas Carteset al.,2019;Gómez Aléset al.,2019).

    Figure 4 Sprint speed of Takydromus lizards (a) and their Topt (b).(a) Arrows indicate the Topt for sprint speed,respectively.Data are expressed as mean ± SEM.Species are listed from high to low latitude in X axis.Data are collected by Plot Digitizer (http://plotdigitizer.sourceforge.net/),from the figures of references (Ji et al.,1996a; Xu and Ji.,2006; Zhang and Ji,2004).

    The thermal biology traits ofT.amurensiswere collected during the reproductive season (i.e.,late April).However,it is possible that seasonal plasticity or acclimatization/acclimation can influence thermal biology traits and their associated performance curve in reptiles and thus cause variations of thermal biology traits inTakydromuslizards.For example,Takydromus septentrionalisfrom different locations did not differ in their selected body temperatures and thermal tolerance after identical thermal acclimation(Yanget al.,2008).Furthermore,the thermal performance curve of the sprint speed ofPlestiodon chinensiswas shown to vary between seasons (Sunet al.,2014).The metabolic enzyme activity in the muscles ofAlligator mississippiensisalso responds to temperatures differently between seasonal acclimation (Seebacheret al.,2003).Determining the seasonal variation of biological traits in cold-climate species such asT.amurensisis important because thermal fluctuations are predicted to be more significant at high latitudes (IPCC 2013).Another potential cause of the observed variation in the thermal biology traits ofT.amurensisin the current study may be the limited sample size (i.e.,n=13).Larger sample sizes of individuals will be required to investigate seasonal acclimation and confirm our preliminary results on latitudinal trends.However,as a cold-climate species,the natural population density ofT.amurensisis low and is threatened by habitat destruction (Zhao,2002; Portniaginaet al.,2019).In future studies,meta-analysis methods may be more suitable for revealing the geographical patterns of thermal biology traits within congeners.It is possible that the inter-specific variation in thermal biology traits along latitudinal gradients is attributed to other factors,such as precipitation,vegetation cover,or food availability,which also varies across altitudes.However,as theTakydromusspecies occupy similar elevations,it is plausible that the observed pattern among species is induced by latitudinal differences.Future studies,investigating more species distributed over a large geographical span,should be conducted to test altitudinal pattern.

    Our study determined the thermal biology traits of cold-climate distributedT.amurensis,summarized the interspecies variation,and provides preliminary findings suggesting that Tsel,CTmax,and CTminvary along latitudinal gradients inTakydromuslizards.These findings will contribute to the understanding of the thermal adaptation of reptiles and establish basic criteria for determining the climatic vulnerability ofTakydromusspecies.Future studies should investigate more details of the biological response to thermal environments,including acclimation and life history variations.Furthermore,building on the data collected in this study with additional measurements of operative temperature and active body temperature in field,we will be able to calculate the thermal safety margin and the effectiveness of thermoregulation,which would be helpful in evaluating the degree of safety or beharioral flexibility to thermal variation,in the context of climate warming.Additional research into thermal biological responses will facilitate projection of the vulnerabilities ofTakydromusspecies to climate warming.

    AcknowledgementsWe thank Xingzhi HAN,Teng LI and Tingting WANG for their assistance in the field and Lab.Ethics approval and protocol (IOZ14001) for the collection,handling,and husbandry of the study animals was given by Animal Ethics Committees at Institute of Zoology,Chinese Academy of Sciences.This work was supported by National Natural Science Foundation of China (31870391 and 31500324).Baojun SUN is supported by Youth Innovation Promotion Association CAS (No.2019085).

    中出人妻视频一区二区| 一进一出好大好爽视频| 久久99一区二区三区| 一区在线观看完整版| 欧美成人免费av一区二区三区| 我的亚洲天堂| 亚洲精品国产区一区二| 久久精品国产清高在天天线| 十八禁人妻一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 人人澡人人妻人| 亚洲国产精品合色在线| 在线观看一区二区三区| 国产精品 国内视频| 一级作爱视频免费观看| 精品国内亚洲2022精品成人| 超碰97精品在线观看| 精品国产美女av久久久久小说| 欧美激情高清一区二区三区| 亚洲第一欧美日韩一区二区三区| 精品国产国语对白av| 亚洲中文字幕日韩| 十八禁网站免费在线| 青草久久国产| 免费在线观看日本一区| 三上悠亚av全集在线观看| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人| 亚洲男人的天堂狠狠| 老司机亚洲免费影院| 欧美日韩亚洲高清精品| 黄色视频不卡| 97超级碰碰碰精品色视频在线观看| 午夜免费鲁丝| 亚洲精品中文字幕在线视频| 国产91精品成人一区二区三区| 性色av乱码一区二区三区2| 老司机午夜十八禁免费视频| 高清欧美精品videossex| 伊人久久大香线蕉亚洲五| 久久精品影院6| 午夜精品在线福利| 国产1区2区3区精品| 黄色 视频免费看| 男女床上黄色一级片免费看| 制服诱惑二区| 夜夜躁狠狠躁天天躁| 日韩中文字幕欧美一区二区| 亚洲熟妇熟女久久| 热99re8久久精品国产| 黄色毛片三级朝国网站| 男人舔女人的私密视频| 丝袜美腿诱惑在线| 一级片'在线观看视频| 757午夜福利合集在线观看| 丁香欧美五月| 一夜夜www| 日本vs欧美在线观看视频| 看免费av毛片| 婷婷精品国产亚洲av在线| 视频区欧美日本亚洲| 女警被强在线播放| 欧美日韩国产mv在线观看视频| 老汉色∧v一级毛片| 人人妻人人爽人人添夜夜欢视频| 两个人看的免费小视频| 夜夜看夜夜爽夜夜摸 | 99riav亚洲国产免费| 亚洲国产中文字幕在线视频| 色婷婷久久久亚洲欧美| 久久欧美精品欧美久久欧美| 国产伦人伦偷精品视频| 欧美成人免费av一区二区三区| 欧美乱码精品一区二区三区| 欧美在线黄色| 又大又爽又粗| 50天的宝宝边吃奶边哭怎么回事| 久久婷婷成人综合色麻豆| 亚洲熟妇中文字幕五十中出 | 午夜视频精品福利| 色综合站精品国产| 日本撒尿小便嘘嘘汇集6| 美国免费a级毛片| 日韩成人在线观看一区二区三区| 欧美成人性av电影在线观看| 成人影院久久| 国产精品98久久久久久宅男小说| 亚洲视频免费观看视频| 午夜精品久久久久久毛片777| 最近最新中文字幕大全电影3 | 精品欧美一区二区三区在线| 亚洲av五月六月丁香网| 国产一卡二卡三卡精品| 欧美色视频一区免费| 看免费av毛片| 亚洲美女黄片视频| a级毛片在线看网站| 精品久久蜜臀av无| 老司机深夜福利视频在线观看| 欧美日本亚洲视频在线播放| 老鸭窝网址在线观看| 国内久久婷婷六月综合欲色啪| 亚洲一区中文字幕在线| 他把我摸到了高潮在线观看| 香蕉国产在线看| 国产99久久九九免费精品| 三上悠亚av全集在线观看| 99在线人妻在线中文字幕| 欧美日韩福利视频一区二区| 9191精品国产免费久久| 国产不卡一卡二| av超薄肉色丝袜交足视频| 国产高清国产精品国产三级| 麻豆成人av在线观看| 午夜免费观看网址| 又大又爽又粗| 欧美日韩亚洲国产一区二区在线观看| www日本在线高清视频| 欧美日韩黄片免| av视频免费观看在线观看| 一个人观看的视频www高清免费观看 | 真人一进一出gif抽搐免费| 日韩大尺度精品在线看网址 | 亚洲精品一二三| 午夜免费激情av| 怎么达到女性高潮| 久久久国产欧美日韩av| 夜夜爽天天搞| 日韩大尺度精品在线看网址 | 国产精品免费一区二区三区在线| 久久精品成人免费网站| 一级毛片精品| 午夜a级毛片| 超色免费av| 黑人欧美特级aaaaaa片| 12—13女人毛片做爰片一| 亚洲av日韩精品久久久久久密| 操出白浆在线播放| 男男h啪啪无遮挡| 精品电影一区二区在线| 免费女性裸体啪啪无遮挡网站| 亚洲第一欧美日韩一区二区三区| 手机成人av网站| 少妇被粗大的猛进出69影院| 国产精品美女特级片免费视频播放器 | 午夜福利在线观看吧| 成人手机av| 欧美日韩福利视频一区二区| www.精华液| 99国产精品99久久久久| 亚洲一区二区三区不卡视频| 999精品在线视频| 久久国产乱子伦精品免费另类| 免费女性裸体啪啪无遮挡网站| 亚洲欧美精品综合久久99| 亚洲欧美日韩另类电影网站| 日本精品一区二区三区蜜桃| 可以免费在线观看a视频的电影网站| 女性被躁到高潮视频| 精品一品国产午夜福利视频| 一边摸一边抽搐一进一小说| 亚洲av片天天在线观看| 亚洲精品一二三| 视频区图区小说| 咕卡用的链子| 少妇被粗大的猛进出69影院| 国产国语露脸激情在线看| 热re99久久国产66热| 久久人人精品亚洲av| 男女下面插进去视频免费观看| 精品一区二区三卡| 久9热在线精品视频| 久久精品国产综合久久久| 99re在线观看精品视频| bbb黄色大片| 757午夜福利合集在线观看| 在线免费观看的www视频| 女警被强在线播放| 在线观看免费视频日本深夜| 亚洲九九香蕉| 我的亚洲天堂| 久久亚洲精品不卡| av片东京热男人的天堂| 亚洲五月色婷婷综合| 国产成年人精品一区二区 | 日韩精品青青久久久久久| 午夜亚洲福利在线播放| 成人18禁高潮啪啪吃奶动态图| 一区二区三区激情视频| 精品免费久久久久久久清纯| 可以在线观看毛片的网站| 午夜影院日韩av| 国产黄a三级三级三级人| 亚洲免费av在线视频| 变态另类成人亚洲欧美熟女 | 在线观看午夜福利视频| 美国免费a级毛片| 国产又爽黄色视频| 波多野结衣一区麻豆| 欧美性长视频在线观看| 国产一区二区三区综合在线观看| 精品日产1卡2卡| 激情视频va一区二区三区| 俄罗斯特黄特色一大片| 看免费av毛片| 久久天堂一区二区三区四区| 亚洲欧美精品综合久久99| 999久久久精品免费观看国产| 亚洲精品国产精品久久久不卡| 亚洲成a人片在线一区二区| 一级a爱视频在线免费观看| 美女午夜性视频免费| 久久精品亚洲av国产电影网| 欧美国产精品va在线观看不卡| 在线永久观看黄色视频| 久久久久国产一级毛片高清牌| 精品福利观看| 国产精品av久久久久免费| 国产高清videossex| 国产精品影院久久| 一级片免费观看大全| 亚洲中文日韩欧美视频| 少妇被粗大的猛进出69影院| 老司机午夜十八禁免费视频| 久久精品91蜜桃| 黄频高清免费视频| 久久国产亚洲av麻豆专区| 午夜影院日韩av| 大香蕉久久成人网| 午夜日韩欧美国产| 91精品三级在线观看| 97超级碰碰碰精品色视频在线观看| 岛国在线观看网站| 国产一区二区三区在线臀色熟女 | 亚洲熟妇中文字幕五十中出 | 亚洲视频免费观看视频| 亚洲第一欧美日韩一区二区三区| 丝袜在线中文字幕| 交换朋友夫妻互换小说| 我的亚洲天堂| 久久精品国产清高在天天线| 亚洲avbb在线观看| av中文乱码字幕在线| 黄片大片在线免费观看| 国产精品1区2区在线观看.| 亚洲av电影在线进入| 国产精品一区二区在线不卡| 国产成人系列免费观看| x7x7x7水蜜桃| 午夜视频精品福利| 亚洲一区二区三区欧美精品| 亚洲一区二区三区不卡视频| 欧美日韩中文字幕国产精品一区二区三区 | 一区在线观看完整版| 无限看片的www在线观看| 在线天堂中文资源库| 色尼玛亚洲综合影院| av天堂在线播放| 久久性视频一级片| 又紧又爽又黄一区二区| 欧美成人性av电影在线观看| 精品久久久精品久久久| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久| 天天躁夜夜躁狠狠躁躁| 久久久久国内视频| 国产伦一二天堂av在线观看| 久久精品人人爽人人爽视色| 两个人看的免费小视频| 国产欧美日韩一区二区三| 成年女人毛片免费观看观看9| 精品国产超薄肉色丝袜足j| 狠狠狠狠99中文字幕| 一级毛片高清免费大全| a级片在线免费高清观看视频| 日本vs欧美在线观看视频| 国产视频一区二区在线看| 午夜亚洲福利在线播放| 国产亚洲精品一区二区www| 国产欧美日韩精品亚洲av| 久久伊人香网站| 男女午夜视频在线观看| 久久久精品国产亚洲av高清涩受| 麻豆久久精品国产亚洲av | 女同久久另类99精品国产91| av网站免费在线观看视频| 麻豆一二三区av精品| 欧美老熟妇乱子伦牲交| 人妻久久中文字幕网| 99国产极品粉嫩在线观看| 欧美精品一区二区免费开放| 女性被躁到高潮视频| 国产一区二区在线av高清观看| 亚洲精品国产一区二区精华液| 成人国产一区最新在线观看| 国产熟女午夜一区二区三区| 亚洲精品中文字幕在线视频| 日韩大码丰满熟妇| 国产精品一区二区精品视频观看| 精品免费久久久久久久清纯| 老司机深夜福利视频在线观看| 日本黄色日本黄色录像| 亚洲国产中文字幕在线视频| 国产精品香港三级国产av潘金莲| 国产单亲对白刺激| 国产深夜福利视频在线观看| 欧美精品一区二区免费开放| 老熟妇乱子伦视频在线观看| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区不卡视频| 香蕉丝袜av| 99国产精品99久久久久| 国产亚洲欧美在线一区二区| 国产精品电影一区二区三区| 18禁国产床啪视频网站| 国产片内射在线| 亚洲国产毛片av蜜桃av| 日本wwww免费看| 老汉色∧v一级毛片| 99riav亚洲国产免费| 777久久人妻少妇嫩草av网站| 亚洲人成77777在线视频| 人人澡人人妻人| av欧美777| 狂野欧美激情性xxxx| 99国产极品粉嫩在线观看| 国产精品美女特级片免费视频播放器 | 日本免费一区二区三区高清不卡 | 久久人人爽av亚洲精品天堂| 午夜a级毛片| 精品久久蜜臀av无| 村上凉子中文字幕在线| 桃色一区二区三区在线观看| 嫩草影视91久久| 麻豆国产av国片精品| 国产亚洲精品久久久久久毛片| 国产伦人伦偷精品视频| 在线十欧美十亚洲十日本专区| 五月开心婷婷网| 国产高清videossex| 男女做爰动态图高潮gif福利片 | 久久青草综合色| 亚洲激情在线av| 欧美中文综合在线视频| 日本 av在线| 高清av免费在线| 亚洲国产精品999在线| 欧美丝袜亚洲另类 | videosex国产| 人妻丰满熟妇av一区二区三区| 久久久久久久午夜电影 | 国产亚洲欧美在线一区二区| 亚洲中文日韩欧美视频| 亚洲性夜色夜夜综合| 国产精品野战在线观看 | 可以在线观看毛片的网站| 99热只有精品国产| 久久久久久久久免费视频了| 天堂√8在线中文| 国产成人av激情在线播放| 一区二区日韩欧美中文字幕| tocl精华| a级毛片黄视频| 国产欧美日韩一区二区三| 性欧美人与动物交配| 最好的美女福利视频网| 97超级碰碰碰精品色视频在线观看| 欧美一区二区精品小视频在线| 国产精华一区二区三区| 亚洲国产欧美日韩在线播放| 亚洲色图 男人天堂 中文字幕| av免费在线观看网站| 最近最新免费中文字幕在线| 免费日韩欧美在线观看| а√天堂www在线а√下载| 国产成+人综合+亚洲专区| 亚洲九九香蕉| 又黄又爽又免费观看的视频| 久久精品91无色码中文字幕| 国产精品亚洲一级av第二区| 中文欧美无线码| 亚洲熟女毛片儿| 超色免费av| 亚洲精品成人av观看孕妇| 久久久久国内视频| 超碰97精品在线观看| 男人操女人黄网站| 黄色 视频免费看| 一区在线观看完整版| 热re99久久精品国产66热6| 久久这里只有精品19| 中文亚洲av片在线观看爽| www.www免费av| 亚洲人成网站在线播放欧美日韩| 在线观看午夜福利视频| 亚洲欧美日韩另类电影网站| 久久天躁狠狠躁夜夜2o2o| 国产精品98久久久久久宅男小说| 久久精品91无色码中文字幕| 亚洲精品在线美女| 一进一出抽搐动态| 午夜影院日韩av| 精品久久蜜臀av无| 久久狼人影院| 国产激情久久老熟女| 婷婷六月久久综合丁香| 最近最新中文字幕大全电影3 | 在线十欧美十亚洲十日本专区| 免费女性裸体啪啪无遮挡网站| 亚洲专区字幕在线| 免费人成视频x8x8入口观看| 黑人巨大精品欧美一区二区mp4| 母亲3免费完整高清在线观看| 久久精品亚洲熟妇少妇任你| 一级,二级,三级黄色视频| 老熟妇乱子伦视频在线观看| 一区二区三区精品91| 亚洲欧美激情在线| 久9热在线精品视频| 男人舔女人的私密视频| 男女之事视频高清在线观看| 岛国在线观看网站| 日韩高清综合在线| 亚洲国产欧美日韩在线播放| 国产国语露脸激情在线看| 亚洲 国产 在线| 性少妇av在线| 日韩精品免费视频一区二区三区| √禁漫天堂资源中文www| 国产不卡一卡二| 亚洲欧洲精品一区二区精品久久久| 精品国产美女av久久久久小说| 久久久久精品国产欧美久久久| 日韩欧美一区二区三区在线观看| 国产高清videossex| 美女高潮喷水抽搐中文字幕| 一级毛片高清免费大全| 18禁美女被吸乳视频| 一边摸一边抽搐一进一出视频| 久久99一区二区三区| 亚洲午夜精品一区,二区,三区| 欧美成人性av电影在线观看| 国产亚洲精品久久久久久毛片| 高清毛片免费观看视频网站 | 国产精品久久久久久人妻精品电影| 亚洲国产精品合色在线| 男人的好看免费观看在线视频 | 人人妻人人澡人人看| 嫩草影视91久久| 久久精品国产99精品国产亚洲性色 | 亚洲一区二区三区色噜噜 | 欧美精品亚洲一区二区| 18禁裸乳无遮挡免费网站照片 | 天堂影院成人在线观看| 成年人免费黄色播放视频| 国产成人免费无遮挡视频| 美女高潮到喷水免费观看| 欧美乱码精品一区二区三区| 一级片免费观看大全| 久久久久久免费高清国产稀缺| 制服诱惑二区| 波多野结衣高清无吗| 中出人妻视频一区二区| a级片在线免费高清观看视频| 国产在线精品亚洲第一网站| 午夜福利在线免费观看网站| 色老头精品视频在线观看| 超色免费av| 亚洲成人久久性| 日本撒尿小便嘘嘘汇集6| 国产精品香港三级国产av潘金莲| 91大片在线观看| 亚洲中文av在线| 在线国产一区二区在线| 亚洲黑人精品在线| 国产99白浆流出| 美女 人体艺术 gogo| 在线视频色国产色| 最近最新中文字幕大全电影3 | 悠悠久久av| 性少妇av在线| 少妇被粗大的猛进出69影院| 国产精品野战在线观看 | 99国产精品一区二区蜜桃av| av免费在线观看网站| 欧美乱妇无乱码| 人人妻,人人澡人人爽秒播| 婷婷丁香在线五月| 99热只有精品国产| 精品久久久精品久久久| 又大又爽又粗| 日韩欧美国产一区二区入口| 欧洲精品卡2卡3卡4卡5卡区| 免费高清视频大片| 99国产精品一区二区三区| 桃红色精品国产亚洲av| 国产男靠女视频免费网站| 多毛熟女@视频| 亚洲精品一二三| 亚洲中文日韩欧美视频| 一级片免费观看大全| 免费在线观看影片大全网站| 久久青草综合色| av中文乱码字幕在线| 国产激情久久老熟女| 亚洲熟妇中文字幕五十中出 | 精品福利永久在线观看| 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9| 男女做爰动态图高潮gif福利片 | 超碰成人久久| 欧美乱妇无乱码| 亚洲成人免费av在线播放| 亚洲欧洲精品一区二区精品久久久| 亚洲精品美女久久久久99蜜臀| 一级毛片精品| 人成视频在线观看免费观看| 欧美日韩黄片免| 久久久国产成人精品二区 | 看免费av毛片| 亚洲七黄色美女视频| 夜夜看夜夜爽夜夜摸 | 久久香蕉精品热| 精品久久久久久电影网| 少妇被粗大的猛进出69影院| 午夜日韩欧美国产| 亚洲一区二区三区欧美精品| 精品电影一区二区在线| 欧美中文综合在线视频| 视频区欧美日本亚洲| 国产精品1区2区在线观看.| 视频区图区小说| 欧美在线黄色| 国产激情欧美一区二区| 国产精品久久久久久人妻精品电影| 妹子高潮喷水视频| 精品无人区乱码1区二区| 午夜福利一区二区在线看| 成人精品一区二区免费| 麻豆av在线久日| 亚洲伊人色综图| 国产高清国产精品国产三级| 久久精品人人爽人人爽视色| 欧美不卡视频在线免费观看 | 涩涩av久久男人的天堂| 多毛熟女@视频| 19禁男女啪啪无遮挡网站| 亚洲一区中文字幕在线| 女警被强在线播放| 黄色女人牲交| 国产熟女xx| 18禁裸乳无遮挡免费网站照片 | 男男h啪啪无遮挡| 国产成人精品无人区| 亚洲精品在线美女| 91麻豆av在线| 美女高潮到喷水免费观看| 欧美精品亚洲一区二区| 国产极品粉嫩免费观看在线| 在线观看免费视频网站a站| 国产99久久九九免费精品| cao死你这个sao货| 18禁观看日本| 免费久久久久久久精品成人欧美视频| 亚洲av成人av| 又黄又爽又免费观看的视频| 国产精品美女特级片免费视频播放器 | 亚洲精品一二三| 亚洲专区字幕在线| 久9热在线精品视频| av欧美777| 亚洲av成人一区二区三| 国产精品偷伦视频观看了| 国产男靠女视频免费网站| 久久人人爽av亚洲精品天堂| 免费观看人在逋| 国产无遮挡羞羞视频在线观看| 久久久久久亚洲精品国产蜜桃av| 久久中文看片网| 欧美 亚洲 国产 日韩一| 欧美乱妇无乱码| 在线观看一区二区三区| 亚洲av片天天在线观看| 久久性视频一级片| 精品一区二区三卡| 国产成人精品久久二区二区免费| 久久天躁狠狠躁夜夜2o2o| 久久人人精品亚洲av| 怎么达到女性高潮| 一级片免费观看大全| 男女午夜视频在线观看| 亚洲精品久久午夜乱码| 国产精品 国内视频| 97碰自拍视频| 中文字幕最新亚洲高清| 99精国产麻豆久久婷婷| 成人特级黄色片久久久久久久| 亚洲一区高清亚洲精品| 日韩三级视频一区二区三区| 一级a爱视频在线免费观看| 十分钟在线观看高清视频www| 午夜福利,免费看| 老司机深夜福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 欧美在线黄色| 亚洲av片天天在线观看| 午夜精品在线福利| 两个人免费观看高清视频| 成人精品一区二区免费| 国产精品成人在线| 亚洲自拍偷在线| 黄色毛片三级朝国网站| 久久精品91无色码中文字幕| 精品久久久久久,| 免费一级毛片在线播放高清视频 |