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

    Adaptive Evolution of the Ventral Scale Microornamentations among Three Snake Species

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

    Kexin WANG,Yutian ZHAO,Chaochao HU,Hong LI,Yanfu QU and Xiang JI

    Jiangsu Key Laboratory for Biodiversity and Biotechnology,College of Life Sciences,Nanjing Normal University,Nanjing 210046,Jiangsu,China

    Abstract Research on the ecological effect of microronamentations on the scale surface in reptiles has been carried out over the past few decades.It is found that the microornamentation pattern in reptiles is related to their habitats.This study examined the wettability on scale surface,as well as the differences in microornamentation on ventral scales from the mid-body region in three snake species,Hypsiscopus plumbea (aquatic),Oocatochus rufodorsata (semi-aquatic) and Elaphe carinata (terricolous).Moreover,the scale specimens were metallized and analyzed using scanning electron microscopy.Our results showed that there are microornamentations on the ventral scale surfaces of the tested species,which showed interspecific differences.To be specific,the aquatic snake shows the narrow,fine and regular denticulations which are connected to reduce friction and dirt shedding.By contrast,the terrestrial snake acquired the wider and shorter denticulation which would render more friction during locomotion but it shows greater water resistance to improve the capacity of dirt shedding and compensate for the disadvantage of short and wide denticulations.Additionally,the denticulation characteristics of the semi-aquatic snake fell in between those of aquatic and terrestrial snakes.Therefore,it is deduced in this study that the ventral scale microornamentations in snakes contribute to ecological adaptation to their preferential microhabitats.

    Keywords Elaphe carinata, Hypsiscopus plumbea, Oocatochus rufodorsata,scale microornamentation,snake,wettabiligy

    1.Introduction

    The scale is the most typical character distinguishing squamates from other amniotes,birds and mammals (Abo-Eleneen and Allam,2011).In squamates,the main function of scales is to assist animals in locomotion and to protect them from mechanical abrasion,solar radiation and dehydration (Abo-Eleneen and Allam,2011; Wegeneret al.,2014).In snakes,scales are composed of the outermost epidermal layer and inner corium layer;of them,the epidermal layer is hard,robust,and inflexible,which is composed of β-keratin,while the outermost structure of scales is called oberhautchen (Baden and Maderson,1970;Kleinet al.,2010; Krey and Farajallah,2013).All oberhautchen cells can directly contact the environment,and they are associated with various shapes and sizes in different taxa,such as the surface profile,imbricate structure,level of cell borders,presence or absence of pores/pits,and the occurrence and size of denticulation on posterior cell margins (Arnold,2002).Typically,the arrangement patterns of these cells and their surface microstructures are defined as the microornamentation(Gower,2003; Ruibal,1968).So far,a large number of studies on the microornamentation in squamates have been carried out using scanning electron microscopy.

    The earliest studies mainly focused on the relationship of microornamentation with taxa (Price,1982; Stewart and Daniel,1975).For instance,Price (1982) had found that the dorsal scale microornamentation patterns of snakes from aquatic,arboreal,burrowing and desert habitats were similar,and concluded that the microornamentation patterns had reflected the phylogenetic relationships instead of the ecological or habitat factors.Thereafter,numerous studies on scale microornamentations in squamates suggested that,the interspecific differences were increasingly correlated with ontogeny,body region,scale size and structure,as well as the taxon (Allam and Abo-Eleneen,2012; Gower,2003; Harvey,1993; Price and Kelly,1989).On the other hand,Renous and Gasc (1989) inferred that parts of the microornamentation patterns were related to the locomotor behavior.Over the past few decades,numerous studies have been carried out to examine the relationship of microornamentation with ecological adaptation (Allamet al.,2017; Arnold,2002; Gower,2003).For instance,lizards from dry habitats or vegetation away from the ground have possessed larger scales to suppress the shine generated by the smooth microornamentation(Arnold,2002).However,Arrigoet al.(2019) used 353 snake species to examine the relationship between dorsal scale microstructures between habitats and found the diversity and complexity of snake scale surface microstructures are related to the phylogeny,rather than their habitats.In addition,Arrigoet al.(2019) found that the scale surface microstructures could provide some level of hydrophobicity but it is still unclear which types of microornamentations are related to being higher hydrophobicity.Nevertheless,Gower (2003) predicts the narrow and regular microornamentation ridges,together with the possible pores/pits,can provide better hydrophobicity and dirt shedding capacity in burrowing snakes.

    The ecological functions of microornamentations are mainly to strengthen the dirt shedding capacity and to reduce the friction during locomotion and shine (Arnold,2002; Gower,2003; Rocha-Barbosa and Moraes-e-Silva,2009).Currently,two competing scientific hypotheses have been proposed to interpret the relationship of ecological adaptation with the microornamentation patterns.For one thing,the demands for friction reduction and dirt shedding may require that smooth microornamentations are favored by the natural selection (Gower,2003).Notably,a more regular surface can effectively reduce the friction produced when squamates are passing through the grassland or narrow burrows (Arnold,2002; Gower,2003; Rocha-Barbosa and Moraes-e-Silva,2009).Such a phenomenon is referred to as the “Lotus-Effect”,which suggests that the regular microstructure is becoming superior in dirt shedding on the surfaces of plant leaves and insect wings(Barthlott and Neinhuis,1997; Wagneret al.,1996).For another,the microornamentation should possess more complex threedimensional (3D) structures to reduce shine (Arnold,2002).For instance,the uneven microornamentation is favorable to reduce shine,thus strengthening the camouflage ability of animals(Arnold,2002).In brief,the microornamentation patterns are controlled by multiple conflicting selection pressures.

    Most previous studies have qualitatively analyzed microornamentations but most of them used a single sample and only a small proportion of them have involved quantitative analysis.For example,the data eigenvalues of the denticulation length are used to phylogenetically analyze the history of scale microornamentation in Lacertid lizards (Arnold,2002).Snake scales,especially for the ventral scales,are frequently,direct and continuous contact with their surrounding environment.They need to provide high friction for forward motion or low friction in sliding along the contact surface (Baumet al.,2014).Therefore,the microornamentation patterns,which display important functions,are related to specific adaptation to the ecological factors like habitat,substrate and locomotion (Klein and Gorb,2012; Krey and Farajallah,2013).

    Scale microornamentations can be regenerated repeatedly,which are similar to those present in the process of embryogenesis (Alibardi,2000).Meanwhile,they can also effectively reduce the wear resistance on natural substrates and help to adapt to the crawling environment (Klein and Gorb,2012).Taken together,these characteristics have made it possible to carry out quantitative study on microornamentation to explain the influence of habitats on snakes.Noteworthily,the arrangement,distribution and shape of microornamentation patterns may probably affect the macrostructures and function of scales in different species,even though they are similar.

    In this study,three most common snakes in China,includingHypsiscopus plumbea,Oocatochus rufodorsataandElaphe carinatawere selected.We use the microornamentation of ventral scales and wettability data for quantitative analysis to verify whether variations in microornamentation patterns on the ventral scales were related to their microhabitats.In these three species,H.plumbeaoccupies pond,river and paddy field,with the distribution scope ranging from south China to Indonesia(Jiang,1998).Meanwhile,O.rufodorsatais a semi-aquatic species,which can be extensively discovered in ditch,as well as the field along river banks from central China to Korea (Zhao,1998).E.carinatais a terrestrial snake species,which can be mainly found in most mountainous,hilly and plains in China and Vietnam(Huang,1998).

    2.Material and methods

    Seven male specimens for each species were collected in this study,among which,H.plumbeawas collected from Fujian Province of China,with the average snout-vent length of(42.90±0.99) cm (range:39.50-46.10 cm).Moreover,O.rufodorsataandE.carinatawere obtained from Zhejiang Province of China,with the average snout-vent lengths of (124.44±5.80) cm(range:108.60-149.30 cm) inO.rufodorsataand (125.49±7.40) cm(range:86.70-144.60 cm) inE.carinata,respectively.The ventral scales of these three snakes were shine and iridescence,but their abdomen colors were different from one to another.Specifically,the abdomen ofH.plumbeawas pale yellow,with a black vertical line in the center,while that ofO.rufodorsatawas yellow with black spots,and that ofElaphe carinatawas almost yellow,with dark spots on the posterior margin.All samples were then preserved in the refrigerator at-20 ℃ in the Laboratory of Nanjing Normal University.

    2.1.WettabiligyThe wettability on material surface is one of the important performance indicators,which can be determined based on the balance of surface energy in the interfaces between air,liquid,and solid materials (Mao,2016).In addition,it can also be evaluated experimentally through surveying the contact angle of one drop of water on the skin surface.Typically,the larger contact angles can be measured on the hydrophobic surfaces.In this study,skin from the center of venter was cut into 1cm×1cm pieces and rinsed with the ultrasonic bath (Bandelin Sonorex RK 100) for 5 min before drying for 24 h.Afterwards,a pipette was used to add 10 μl distilled water droplets onto the abdominal scales placed onto the microscope slides at room temperature.Subsequently,the camera was fixed in parallel with the slide,and the sample was photographed and labeled.Later,all photographs were measured for contact angle using the ImageJ version 1.43u.

    2.2.MicroornamentationAll ventral scales were selected from the midbody regions of the specimens,and cleaned before examination.Similarly,skin collected from the center of venter scales was cut into 5 mm×5 mm pieces and mounted on the microscope slides with string.Then,the samples were rinsed with distilled water for 20 min,fixed in 5% formalin for 24 h,rinsed for 1 h again,and cleaned in the ultrasonic bath for 5 min before they were immersed into a consecutive series of ethanol solutions (30%,50%,70%,80%,95% and 100%).Afterwards,the samples were maintained in each solution for 15 min at room temperature,followed by air-drying for 24 h and then affixing with the assistant double-sided adhesive tape to mount onto the aluminum scanning electron microscope(SEM) stubs.Subsequently,the mounted samples were coated with gold-palladium and observed under the SEM (JSM-5610LV) at the accelerating voltage of 8-10 kV,so as to select the appropriate magnifications for characterization and further analysis.Additionally,the dimensions of microornamentations were also measured in the SEM images using the ImageJ version 1.43u.Specifically,three pictures were obtained for each sample,and the average of every trait was calculated according to the 150 denticulations randomly selected in three pictures.Moreover,to quantitatively describe the microscopic morphology of the ventral scale,the length,width,spacing,area and row spacing of denticulations were measured,respectively.Additionally,the height of scales (from the base to the edge of scales) was also determined since the size of scale could affect that of denticulations.Figure 1 has depicted the measurement methods of these denticulation traits.

    Figure 1 Schematic diagram of the characteristic parameters of denticulation.The width of denticulation refers to the upper width of denticulation,as denoted by “W”.The length of denticulation represents the vertical distance from top to bottom of denticulation,as indicated by “L”.The center distance of two denticulations stands for the distance between two denticulations,as suggested by “D”.The area of denticulation is the area enclosed by “abc”.The row spacing suggests the distance between the two rows of denticulation.W:Width; L:Length; D:Distance.

    2.3.Data analysisThe STATISTICA 10.0 (StatSoft Inc.,Tulsa,Oklahoma,USA) was adopted for all statistical analyses,and the significance level was set at 0.05.Prior to parametric analyses,all data were tested for the homogeneity of variance using the Bartlett’s test,and the normally distributed data were analyzed using the Kolmogorov-Smirnov test.Meanwhile,the oneway analysis of variance (ANOVA) was utilized to examine the differences in the contact angle,the aspect ratio and the area of denticulations among these three species,as well as the difference in the length of denticulations (in which the height of denticulations was used as the covariate,while species as a factor),and the differences in the width,spacing and row spacing of denticulations (in which the length of denticulations was used as the covariate,whereas species as a factor).Besides,Tukey’s test was employed to test the post hoc comparisons for variables that were different among these three species.Throughout this article,all values are presented as the mean ±standard error (SE),and the significance level is set at a=0.05.

    3.Results

    3.1.WettabilityOur results suggested that contact angles were (91.20±2.03)° (80.88-96.62),(78.68±2.63)° (70.89-90.13) and(68.63±2.56)° (60.82-77.01) forE.carinata,O.rufodorsataandH.plumbea,respectively.Obviously,the terrestrial snake (E.carinata)had the largest contact angle,while the aquatic snake (H.plumbea) had the lowest one (F2,18=21.79,P<0.0001).

    3.2.Qualitative traits of scale microstructureNotably,the microornamentation on the scale surfaces of snakes seemed to be related to their living environments (Rocha-Barbosa and Moraes-e-Silva,2009).Therefore,we investigated whether the significant difference in scale microornamentations was found in three kinds of snakes.All ventral scales had no sensilla on the surface and the keel only existed on scales ofO.rufodorsataandH.plumbea.Under high magnification,the surface microstructures of three kinds of snakes displayed similar characteristics.As for the microscopic surface morphology in three kinds of snakes,a regularly arranged surface structure could be observed,which was composed of the small denticulations and micropores (Figure 2).Typically,the denticulations were caudally oriented and paralleled to the body axis,while the pores and pore-like structures were found on the scales in the three kinds of snakes,which were located at the junction of two adjacent denticulations(Figures 2B,2D and 2F).

    Additionally,the denticulation ofH.plumbeashowed a triangular structure,along with a long and thin tip (Figure 2B),and the distance between adjacent denticulations was close to each other,with almost no plate structure (Figure 2B).Similarly,the denticulation ofO.rufodorsataalso exhibited a triangular structure,but its length was shorter than the row spacing,which had resulted in a “plate” structure (Figure 2D).By contrast,the denticulation ofE.carinataappeared to be a finger-like structure,which was shorter and thicker,with a more obvious plate structure (Figure 2F).

    Figure 2 Ventral scales microstructure of ventral scales.A and B display the different magnifications of the ventral scales of Hypsiscopus plumbea.C and D exhibit the different magnifications of the ventral scales of Oocatochus rufodorsata.E and F present the different magnifications of the ventral scales of Elaphe carinata.Ke:Keel; De:Denticulation; Po:Pore; VS:Vertical stripe.

    Table 1 Values of surface microstructure characteristics for the ventral scale.

    3.3.Quantitative analysis of scale denticulationThe microscopic morphology of denticulations is depicted in Table 1,and the statistical analysis results were determined based on the measurement data.Using the height of scale as a covariate,the one-way ANCOVA revealed that there were significant differences in the length of denticulations among the three kinds of species (F2,17=5.17,P=0.02).Specifically,the denticulations inH.plumbea,andO.rufodorsatawere longer than that inE.carinata(Table 1).Besides,the width (F2,17=5.16,P=0.02) and row spacing (F2,17=10.90,P<0.001) of denticulations successively increased withH.plumbea,O.rufodorsata and E.carinata(Table 1).Obviously,the denticulation spacing ofH.plumbeawas shorter than those of the other two kinds of snakes (F2,17=4.79,P=0.02) but the difference of the denticulation spacing betweenO.rufodorsata and E.carinatawas not statistically significant (Table 1).In addition,there was not the difference in the area of denticulations among the three kinds of snakes(F2,18=3.53,P=0.05; Table 1).On the other hand,the width-tolength ratio of denticulations successively elevated withH.plumbea,O.rufodorsata and E.carinata(F2,18=46.13,P< 0.001;Table 1).

    4.Discussion

    Snakes live in a wide spectrum of habitats,and their limbless locomotor and/or behavioral specializations have reflected the differences in epidermal architecture,which are correlated with their predatory strategies and properties of the contact surface where they move (Klein and Gorb,2014; Rocha-Barbosa and Moraes-e-Silva,2009).The shape and arrangement of microornamentation may affect the friction and wear on the scale surface when the animal moves.Therefore,the length of denticulations,the number of spinule rows,and the presence or absence of pits/pores may be correlated with the preferred habitats of snakes (Gower,2003; Rocha-Barbosa and Moraese-Silva,2009).In addition,Klein and Gorb (2012) discovered that the epidermal thickness might be related to the specific adaptation to animal habitat and locomotion.On this account,scales and microornamentations of their surface not only represented the microstructure of animal,but also stood for its macrostructure (Rocha-Barbosa and Moraes-e-Silva,2009;Smithet al.,1982).

    Furthermore,the variations in microornamentation patterns of squamates have been extensively observed in intra-species and inter-species (Allam and Abo-Eleneen,2012;Allamet al.,2017; Gower,2003; Rocha-Barbosa and Moraes-e-Silva,2009).Such differences in microornamentations allow to explain the functional adaptation of microornamentations to their surrounding environment (Stewart and Daniel,1973).On the other hand,Williams and Peterson (1982) had proposed a novel evolution adaptation complex between adaptation and the morphological divergences of scales.In addition,Gower(2003) suggested that the variations in microornamentations could be ascribed to the selection pressure of various habitats.Besides,other studies also indicate that microornamentation variations on the scale surfaces are related to the natural selection pressures of their habitats (Allam and Abo-Eleneen,2012; Allamet al.,2017; Rocha-Barbosa and Moraes-e-Silva,2009). After analyzing surface nanostructures in 353 species spanning 19 of the 26 families of snakes and performing phylogenetic mapping on the snake phylogeny,Arrigoet al.(2019) did not find a correlation between nanomorphological characters and a simple life habit classification.Hence,these authors suggested that the diversity and complexity of snake skin surface nano-morphology are dominated by phylogenetic rather than habitat-specific functional constraints.

    It was found in this study that the ventral scales of the three investigated snakes had attractive iridescence and faint reflections but these colors were almost indistinguishable under the naked eye.Notably,the contact angle of one water droplet on the material surface is indicative of the wettability of the contact surface (Spinneret al.,2013).In this study,the contact angles of the ventral scales of three kinds of snakes were compared,which suggested that the terrestrial snake (E.carinata) had the greatest contact angle whereas the aquatic snake (H.plumbea) had the lowest one.Additionally,the microornamention patterns were similar among the three snake species,which were made up of digitations and pits(Figure 2).The presence of cell border digitations is the ancestral state for snake skin microstructures which were subsequently and independently lost in multiple lineages (Arrigoet al.,2019).The following phenomena could be found after comparing the microornamentation data of ventral scales among the three species.(1) terrestrial snakes (E.carinata) had the shortest denticulation length,while aquatic snakes (H.plumbea) had the longest one (Table 1 and Figure 2); (2) the denticulation width of terrestrial snakes was remarkably larger than those of aquatic and semi-aquatic snakes when the denticulation length was set at the same length (Table 1 and Figure 2); (3) the column space between the denticulations of terrestrial and aquatic snakes was greater than that of aquatic snakes (Table 1 and Figure 2);and (4) the row space between the denticulations of terrestrial snakes was the greatest,while that of aquatic snakes was the smallest (Table 1 and Figure 2).

    The big scale,especially for the ventral scale of snake,can generate shine (Arnold,2002),while the microornamentation contributes to reducing shine and generating iridescence,since it can form a reflective surface to reflect light rays in different ways or to scatter light,thus restricting or eliminating shine(Arnold,2002).Noteworthily,various mechanisms have been proposed to explain the iridescence in nature,such as the thinfilm interference,multilayer reflectors,photonic crystals,and diffraction gratings (Kinoshitaet al.,2008; Kroisset al.,2009;Seagoet al.,2009).On the other hand,the iridescent colors mainly allow the animals to communicate with conspecifics and avoid predators (Doucet and Meadows,2009).In squamates,a fine arrangement of microornamentation can reduce the shine on the scale surface and strengthen the dirt shedding capacity (Arnold,2002; Gower,2003).However,many researchers have supposed that iridescence is only regarded as a by-product of microornamentations,and there is no plausible explanation for the adaptive evolution.For instance,the fossorial snake spends most of their life time under the ground (Gower,2003).In addition,iridescent color has not been observed inT yphlops mirus(Gower,2003).Similarly,our results suggested no significant differences in shine and iridescence among these three kinds of snakes,although the microornamentation patterns among them showed significant differences in our study.

    The wettability of scale surface can be determined by the scale-liquid contact angle,which represents a quantitative measure of the wettability of a solid surface by a liquid.Differences in the wettability between animal skins can be ascribed to the chemical composition properties and the geometric structure on the skin surface (Spinneret al.,2013;Arrigoet al.,2019).In squamates,the outermost structure of scales is referred to as oberhautchen (Irishet al.,1988),and microornamentation includes the surface features of oberhautchen and the epidermal folding (Arnold,2002).Hence,the 3D structure and arrangement of microornamentation will affect the wettability of the scale surface.Specifically,the microornamentation,especially for the denticulations and ridges,can improve the wettability of scale surface,as previously reported in numerous publications (Autumn and Hansen,2006; Gower,2003; Spinneret al.,2014; Spinneret al.,2013; Arrigoet al.,2019).

    Generally,a contact angle of < 90° (low contact angle)has indicated the quite favorable wetting of the surface,and the fluid can spread over a large surface area.By contrast,a contact angle of > 90° (high contact angle) is indicative of the unfavorable wetting of the surface.Therefore,the fluid has minimal contact with the surface and will form a compact liquid droplet (Shiet al.,2016).Nowadays,the contact angle of scale surface has been measured in some species (Gower,2003;Spinneret al.,2013; Arrigoet al.,2019).According to our results,terrestrial snakes had the largest contact angle of > 90°,while aquatic snakes had the smallest one,and the contact angle of terrestrial and semi-aquatic snakes was < 90°.Therefore,it could be deduced that the scales of terrestrial snakes had higher water resistance.Besides,Gower (2003) suggested that the contact angles were > 90° forMelanophidium wynaudense,Uropeltis spp.,Rhinophis philippinus,andCylindrophis maculatus,while those were < 90° forTyphlops mirusandT.reticulatus.

    The microornamentation has possessed the narrow and regular ridges,which may be more advantageous in reducing friction and dirt shedding.Meanwhile,these microornamentation patterns may be correlated with the passively reduced adhesion of moist substrate onto the scales,since lower energy is required on the scale surface (Gower,2003).For example,the microornamentation with narrow,finely and regularly spaced ridges is correlated with the less wettable scales in Uropeltid snakes (Gower,2003); on this account,Uropeltid snakes had superior dirt shedding capacity under moist conditions (Gower,2003).In this study,aquatic snakes had the longest,while the finest and narrowest ridges,thus,they might perform better in dirt shedding (Table 1 and Figure 2).

    Scales of snakes,especially for terrestrial snakes,are involved in locomotion to promote movement through muscle activity(Wegeneret al.,2014).Specifically,the number and geometry structure of microornamentation will change the friction between scales and substrate,in the meantime of protecting their bodies from abrasion (Alibardi,2003; Oufieroet al.,2011).To enable the undulating locomotion,the snake skin has to offer high friction to support forward motion whereas low enough friction to enable the sliding along the substrate simultaneously (Baumet al.,2014).Many researchers have discovered that the greater frictional resistance of scales from the posterior to the anterior direction is expected when the long denticulation has overlapped the posterior borders of each spinule row (Gower,2003; Rocha-Barbosa and Moraes-e-Silva,2009).Different from previous results,we found that terrestrial snakes (E.carinata) owned the greater denticulation row space,along with the largest denticulation width.Therefore,it was supposed that increasing the denticulation width contributed to increasing the friction,which could in turn provide the snakes with the strong power to glide forward.Additionally,our results indicated that,terrestrial snakes could acquire friction during locomotor,however,they had to sacrifice their dirt shedding ability.It could be discovered based on the contact angle data that,it might be possible to change the scale properties and improve the water-resistance to compensate for the lack of dirt shedding.

    In conclusion,these three snake species investigated in this study have exhibited significant differences in the wettability and microornamentation pattern on the ventral scale surface,suggesting that these differences may actually be ascribed to the natural selection pressures and the tradeoff between wettability and friction,which may ultimately result in changes of the skin properties for terrestrial snakes.Besides,such differences are conducive to the survival of snakes,which help them to occupy more ecological niches.This study may potentially provide important quantitative data to explore a succession of questions about the relationship between the diversity and ecological adaptation of habitats among snakes.Despite this,we are still unable to draw definitive conclusions indicating that the microstructures of the three snakes are related to their habitats,the current results may also be related to the small number of species.Therefore,it is also necessary to increase the species and the exact physiological,mechanical and behavioral studies to confirm whether this result can be verified in a larger species range.

    AcknowledgementsThis study was carried out in compliance with the current laws of China,and was supported by grants from Natural Science Foundation of China(31200283 and 31770443),and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD)of Jiangsu Higher Education Institutions.The authors would like to thank Kun GUO and Li MA for help in sample collection during the research.

    亚洲精品日本国产第一区| 午夜激情福利司机影院| 国产精品蜜桃在线观看| www.av在线官网国产| av网站免费在线观看视频| 国产精品秋霞免费鲁丝片| 精品少妇久久久久久888优播| 国内精品宾馆在线| 成人综合一区亚洲| 永久免费av网站大全| 又大又黄又爽视频免费| 一区二区三区四区激情视频| 青春草国产在线视频| 在线观看免费视频网站a站| 黄色日韩在线| 99热网站在线观看| 一本大道久久a久久精品| 亚洲不卡免费看| 尾随美女入室| 国产极品天堂在线| 久久人人爽av亚洲精品天堂| 妹子高潮喷水视频| 亚洲国产精品999| 久久久亚洲精品成人影院| xxx大片免费视频| 亚洲国产最新在线播放| 免费黄频网站在线观看国产| 欧美日韩一区二区视频在线观看视频在线| 18禁裸乳无遮挡动漫免费视频| 91精品国产国语对白视频| 亚洲国产欧美日韩在线播放 | 国产精品成人在线| 欧美性感艳星| 91精品国产九色| 色5月婷婷丁香| 亚洲精品视频女| 日韩av免费高清视频| 国产免费一区二区三区四区乱码| 极品少妇高潮喷水抽搐| 夫妻性生交免费视频一级片| 99热网站在线观看| 在线 av 中文字幕| 日韩一区二区视频免费看| 国产精品一区二区性色av| 色吧在线观看| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 高清黄色对白视频在线免费看 | 99久久精品一区二区三区| 欧美+日韩+精品| 99re6热这里在线精品视频| 亚洲情色 制服丝袜| 春色校园在线视频观看| 晚上一个人看的免费电影| 国产精品人妻久久久久久| 色网站视频免费| 日本与韩国留学比较| 国产成人免费观看mmmm| 亚洲av.av天堂| 一本—道久久a久久精品蜜桃钙片| 啦啦啦在线观看免费高清www| 欧美性感艳星| 婷婷色综合www| 韩国av在线不卡| 内地一区二区视频在线| 亚洲经典国产精华液单| 爱豆传媒免费全集在线观看| 曰老女人黄片| 日韩电影二区| 亚洲真实伦在线观看| 久久久久人妻精品一区果冻| 极品教师在线视频| 肉色欧美久久久久久久蜜桃| 精品久久久久久久久av| 久久99热6这里只有精品| 国产毛片在线视频| 免费播放大片免费观看视频在线观看| 久久 成人 亚洲| 久久久久久久国产电影| 国产av码专区亚洲av| 亚洲一级一片aⅴ在线观看| 日韩免费高清中文字幕av| videos熟女内射| 国产精品久久久久久久电影| √禁漫天堂资源中文www| 国产黄色视频一区二区在线观看| 少妇的逼水好多| 日产精品乱码卡一卡2卡三| 青春草亚洲视频在线观看| 日本黄色片子视频| 国产黄片美女视频| 国产成人精品婷婷| 日韩一区二区三区影片| 国产精品一区二区在线观看99| 菩萨蛮人人尽说江南好唐韦庄| 春色校园在线视频观看| 国产精品不卡视频一区二区| 看非洲黑人一级黄片| 午夜福利视频精品| 中文字幕久久专区| av在线观看视频网站免费| 免费大片18禁| 18禁动态无遮挡网站| 成人亚洲精品一区在线观看| 国产一区亚洲一区在线观看| 国产精品久久久久久精品电影小说| 国产欧美日韩精品一区二区| 欧美日韩视频精品一区| 最近手机中文字幕大全| 熟女av电影| 如日韩欧美国产精品一区二区三区 | h日本视频在线播放| 99热国产这里只有精品6| 99久久综合免费| 国产日韩欧美亚洲二区| 日韩av不卡免费在线播放| 我要看日韩黄色一级片| 在线免费观看不下载黄p国产| 色5月婷婷丁香| 精品国产乱码久久久久久小说| 嫩草影院入口| 99久久精品一区二区三区| 国产亚洲5aaaaa淫片| 黄色怎么调成土黄色| 久久精品国产亚洲av天美| 伦理电影大哥的女人| 免费人妻精品一区二区三区视频| 看十八女毛片水多多多| 亚洲在久久综合| 大陆偷拍与自拍| 日韩精品有码人妻一区| 久久99热这里只频精品6学生| 午夜激情福利司机影院| 欧美精品高潮呻吟av久久| 成人免费观看视频高清| 日韩欧美 国产精品| 日韩中字成人| 国产精品免费大片| 嫩草影院入口| 久久精品久久久久久久性| 亚洲国产毛片av蜜桃av| 女性生殖器流出的白浆| 亚洲精品日韩在线中文字幕| 久久久久久伊人网av| 人人澡人人妻人| 亚洲国产色片| 高清毛片免费看| 国产伦精品一区二区三区四那| 寂寞人妻少妇视频99o| 日本av免费视频播放| 国产免费视频播放在线视频| 人妻 亚洲 视频| 国产精品无大码| 嫩草影院入口| 91精品一卡2卡3卡4卡| 最黄视频免费看| 性色avwww在线观看| 在线观看免费视频网站a站| 欧美日韩在线观看h| 国产成人a∨麻豆精品| 黄色欧美视频在线观看| 亚洲真实伦在线观看| 亚洲高清免费不卡视频| 欧美老熟妇乱子伦牲交| 国内揄拍国产精品人妻在线| 最近2019中文字幕mv第一页| av福利片在线观看| 国产日韩一区二区三区精品不卡 | 欧美日韩视频高清一区二区三区二| 另类精品久久| 高清视频免费观看一区二区| 中文天堂在线官网| 十分钟在线观看高清视频www | 建设人人有责人人尽责人人享有的| 新久久久久国产一级毛片| 欧美高清成人免费视频www| 高清午夜精品一区二区三区| 国产精品国产三级专区第一集| 一区二区av电影网| 人人澡人人妻人| 人妻少妇偷人精品九色| 国产淫语在线视频| 男女边摸边吃奶| 精品少妇久久久久久888优播| 亚洲av在线观看美女高潮| 国产女主播在线喷水免费视频网站| 国产精品一区二区在线不卡| av在线app专区| 日韩一本色道免费dvd| 观看美女的网站| 伊人久久国产一区二区| 日本与韩国留学比较| 精品国产国语对白av| 在线观看免费视频网站a站| 如日韩欧美国产精品一区二区三区 | 亚洲内射少妇av| 啦啦啦啦在线视频资源| 99精国产麻豆久久婷婷| 成人18禁高潮啪啪吃奶动态图 | 午夜视频国产福利| 极品人妻少妇av视频| a 毛片基地| 国产探花极品一区二区| 欧美精品国产亚洲| 九九爱精品视频在线观看| 国产精品一区二区性色av| 免费大片18禁| 国产av一区二区精品久久| 精华霜和精华液先用哪个| 亚洲四区av| 中文资源天堂在线| 一区二区av电影网| 亚洲av综合色区一区| 黑丝袜美女国产一区| 国产亚洲午夜精品一区二区久久| 亚洲av在线观看美女高潮| 久久免费观看电影| 在现免费观看毛片| 亚洲av二区三区四区| 97在线视频观看| 午夜精品国产一区二区电影| 免费人妻精品一区二区三区视频| 性色avwww在线观看| 亚洲国产精品999| 99久久中文字幕三级久久日本| 少妇人妻精品综合一区二区| 天堂中文最新版在线下载| 综合色丁香网| 欧美变态另类bdsm刘玥| 国产欧美另类精品又又久久亚洲欧美| 乱人伦中国视频| 男女啪啪激烈高潮av片| 国语对白做爰xxxⅹ性视频网站| 在线播放无遮挡| 精品一区二区免费观看| 亚洲av综合色区一区| 国产成人精品久久久久久| 亚洲色图综合在线观看| 男人添女人高潮全过程视频| 成人国产av品久久久| 夫妻午夜视频| 99久久人妻综合| 男女边摸边吃奶| av国产久精品久网站免费入址| 国产伦理片在线播放av一区| 午夜久久久在线观看| 亚洲精品一区蜜桃| 国产 一区精品| 男人爽女人下面视频在线观看| 国产一区二区三区综合在线观看 | 国产一区二区三区综合在线观看 | 欧美日韩综合久久久久久| 成人影院久久| 中文字幕久久专区| 99久久中文字幕三级久久日本| 久久精品久久久久久久性| 国产成人精品无人区| 国产成人a∨麻豆精品| a级片在线免费高清观看视频| 亚洲av男天堂| 国产精品蜜桃在线观看| 国产一区二区三区av在线| 中文乱码字字幕精品一区二区三区| 久久久国产精品麻豆| av一本久久久久| 少妇精品久久久久久久| 一区二区av电影网| 久久久久久久国产电影| 如日韩欧美国产精品一区二区三区 | 久久午夜福利片| 国产成人午夜福利电影在线观看| 2018国产大陆天天弄谢| 在线观看人妻少妇| 极品人妻少妇av视频| 久久精品久久久久久久性| 成人毛片a级毛片在线播放| 久久国产亚洲av麻豆专区| 少妇的逼水好多| 国产精品国产三级国产av玫瑰| 狠狠精品人妻久久久久久综合| 男女啪啪激烈高潮av片| 精品一区二区三区视频在线| 男女免费视频国产| 精品久久久精品久久久| 亚洲中文av在线| 国产亚洲最大av| 亚洲人成网站在线播| 午夜福利在线观看免费完整高清在| 日日撸夜夜添| 色网站视频免费| 中国三级夫妇交换| 国产精品久久久久成人av| 午夜影院在线不卡| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 婷婷色av中文字幕| 亚洲久久久国产精品| 成人国产麻豆网| 成人无遮挡网站| 高清毛片免费看| 夜夜看夜夜爽夜夜摸| 成人国产av品久久久| 久久这里有精品视频免费| 亚洲精品自拍成人| 青春草亚洲视频在线观看| 亚洲精品456在线播放app| 国产精品一区二区在线观看99| av免费在线看不卡| 中文乱码字字幕精品一区二区三区| 丝袜喷水一区| 亚洲综合色惰| 精品亚洲乱码少妇综合久久| 国产一区二区在线观看日韩| 纯流量卡能插随身wifi吗| 亚洲av国产av综合av卡| 久久久久国产精品人妻一区二区| 伊人亚洲综合成人网| 99久久精品国产国产毛片| 亚洲av二区三区四区| 亚洲情色 制服丝袜| 80岁老熟妇乱子伦牲交| 国产男女内射视频| 美女中出高潮动态图| 十分钟在线观看高清视频www | 国产一区二区三区av在线| 久久精品夜色国产| 精品人妻偷拍中文字幕| 国产老妇伦熟女老妇高清| 亚洲成人av在线免费| 亚洲精品久久午夜乱码| 亚洲国产精品999| 精品一区在线观看国产| h视频一区二区三区| 中国三级夫妇交换| 成人美女网站在线观看视频| 一级毛片 在线播放| 中文字幕制服av| 深夜a级毛片| 国产日韩欧美视频二区| 男女边摸边吃奶| 精品人妻熟女av久视频| 大话2 男鬼变身卡| 国产片特级美女逼逼视频| 人体艺术视频欧美日本| 嫩草影院入口| 国产有黄有色有爽视频| 美女cb高潮喷水在线观看| 女人久久www免费人成看片| 久久99热这里只频精品6学生| 亚洲国产精品999| 精品一区在线观看国产| 亚洲国产精品999| 精品一区在线观看国产| 尾随美女入室| 欧美少妇被猛烈插入视频| 日韩精品免费视频一区二区三区 | 在线看a的网站| 能在线免费看毛片的网站| 成人亚洲精品一区在线观看| 18禁在线无遮挡免费观看视频| 亚洲av二区三区四区| 啦啦啦视频在线资源免费观看| 国产成人一区二区在线| 黄色怎么调成土黄色| 国产探花极品一区二区| 亚洲,一卡二卡三卡| 十八禁网站网址无遮挡 | 亚洲经典国产精华液单| 最后的刺客免费高清国语| 99视频精品全部免费 在线| 国模一区二区三区四区视频| 大香蕉久久网| 激情五月婷婷亚洲| 日日啪夜夜爽| h日本视频在线播放| 如何舔出高潮| 美女大奶头黄色视频| 一区二区三区免费毛片| 日韩视频在线欧美| 国产黄频视频在线观看| 大片免费播放器 马上看| 色5月婷婷丁香| 亚洲内射少妇av| 久久久久久久久久久丰满| 国产黄色免费在线视频| 一区二区三区免费毛片| 免费人成在线观看视频色| 69精品国产乱码久久久| 男人爽女人下面视频在线观看| 日日摸夜夜添夜夜添av毛片| 99热网站在线观看| 在线 av 中文字幕| 麻豆成人午夜福利视频| 伊人久久精品亚洲午夜| 中文精品一卡2卡3卡4更新| 在线免费观看不下载黄p国产| 久久久久久久久久成人| 99热这里只有是精品50| 国产一区二区三区av在线| 美女cb高潮喷水在线观看| 欧美激情国产日韩精品一区| 国产成人91sexporn| 国产无遮挡羞羞视频在线观看| 国产免费一级a男人的天堂| 中国国产av一级| 高清在线视频一区二区三区| 日韩人妻高清精品专区| 亚洲精品乱码久久久久久按摩| 99热国产这里只有精品6| 99久久综合免费| 80岁老熟妇乱子伦牲交| 亚洲精品,欧美精品| 国产 精品1| 最近中文字幕2019免费版| 一级a做视频免费观看| 三上悠亚av全集在线观看 | 日韩中文字幕视频在线看片| 卡戴珊不雅视频在线播放| 欧美变态另类bdsm刘玥| 国产一区二区在线观看日韩| 日本欧美国产在线视频| 久久久欧美国产精品| 丰满人妻一区二区三区视频av| 免费看日本二区| 国产成人精品久久久久久| 一区二区三区乱码不卡18| 久久久久久久亚洲中文字幕| 中文欧美无线码| 国产精品国产三级国产专区5o| 亚洲欧美精品自产自拍| 国产亚洲最大av| 久热这里只有精品99| 久久久久视频综合| 晚上一个人看的免费电影| 男人舔奶头视频| 精品卡一卡二卡四卡免费| 欧美激情国产日韩精品一区| 亚洲精品国产成人久久av| 欧美97在线视频| 少妇人妻久久综合中文| av女优亚洲男人天堂| 久久久久国产网址| 日本av手机在线免费观看| 在线观看免费视频网站a站| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产亚洲av天美| av福利片在线| 99热6这里只有精品| 久久99一区二区三区| 一本久久精品| 亚洲欧美精品自产自拍| 美女cb高潮喷水在线观看| 成人免费观看视频高清| 2018国产大陆天天弄谢| 丁香六月天网| 亚洲av成人精品一二三区| 久久久久久久大尺度免费视频| 日本与韩国留学比较| 国产老妇伦熟女老妇高清| av福利片在线观看| 久久人人爽av亚洲精品天堂| 国产一区有黄有色的免费视频| 国产黄频视频在线观看| av不卡在线播放| 免费看光身美女| 色婷婷av一区二区三区视频| 啦啦啦啦在线视频资源| 我要看黄色一级片免费的| 你懂的网址亚洲精品在线观看| 夜夜骑夜夜射夜夜干| 寂寞人妻少妇视频99o| 中文乱码字字幕精品一区二区三区| 久久热精品热| 欧美3d第一页| 亚洲综合色惰| 日本91视频免费播放| xxx大片免费视频| 久热这里只有精品99| 亚洲精品国产色婷婷电影| 国产在线视频一区二区| 一本久久精品| 欧美精品一区二区大全| a级毛色黄片| 国产伦精品一区二区三区视频9| 亚洲av成人精品一区久久| 男女边吃奶边做爰视频| 欧美激情国产日韩精品一区| 亚洲欧美清纯卡通| 国产成人精品久久久久久| 久久精品夜色国产| 精华霜和精华液先用哪个| 免费黄频网站在线观看国产| 黑人猛操日本美女一级片| 高清视频免费观看一区二区| 亚洲成人手机| 国内少妇人妻偷人精品xxx网站| 久久精品久久久久久久性| 亚洲av成人精品一二三区| 精品久久久噜噜| 亚洲国产精品一区三区| 女人久久www免费人成看片| 久久国内精品自在自线图片| 欧美精品一区二区免费开放| av黄色大香蕉| 亚洲av电影在线观看一区二区三区| 亚洲中文av在线| 国产淫片久久久久久久久| 国产午夜精品久久久久久一区二区三区| 国产乱人偷精品视频| 亚洲欧洲日产国产| 国产白丝娇喘喷水9色精品| 老司机亚洲免费影院| 久久99蜜桃精品久久| 伊人久久国产一区二区| 欧美少妇被猛烈插入视频| 女人久久www免费人成看片| 男人和女人高潮做爰伦理| av播播在线观看一区| 十八禁高潮呻吟视频 | 日本黄色日本黄色录像| 老司机亚洲免费影院| 亚洲熟女精品中文字幕| 久久国产精品男人的天堂亚洲 | 国产成人freesex在线| 两个人免费观看高清视频 | 久久精品久久精品一区二区三区| a级一级毛片免费在线观看| av又黄又爽大尺度在线免费看| 国产欧美日韩一区二区三区在线 | 欧美区成人在线视频| 欧美日韩视频高清一区二区三区二| 久久这里有精品视频免费| 一本—道久久a久久精品蜜桃钙片| 中文字幕亚洲精品专区| 2022亚洲国产成人精品| 51国产日韩欧美| 中文天堂在线官网| 亚洲四区av| 日本免费在线观看一区| 99九九在线精品视频 | 一本色道久久久久久精品综合| 国产精品蜜桃在线观看| 亚洲真实伦在线观看| 日韩 亚洲 欧美在线| 成年美女黄网站色视频大全免费 | 久久毛片免费看一区二区三区| 亚洲精品色激情综合| 久久久久久人妻| 免费人妻精品一区二区三区视频| 国产成人aa在线观看| 亚洲精品aⅴ在线观看| 2021少妇久久久久久久久久久| 成人亚洲欧美一区二区av| 国产av国产精品国产| 18+在线观看网站| 激情五月婷婷亚洲| 欧美激情极品国产一区二区三区 | 亚洲av成人精品一二三区| 国产日韩欧美在线精品| 在线观看国产h片| 免费av不卡在线播放| 老司机亚洲免费影院| 婷婷色麻豆天堂久久| 如何舔出高潮| 亚洲av福利一区| 两个人免费观看高清视频 | 国产探花极品一区二区| 午夜91福利影院| 亚洲精品亚洲一区二区| 欧美 日韩 精品 国产| 亚洲综合色惰| 亚洲丝袜综合中文字幕| 久久国产亚洲av麻豆专区| 亚洲电影在线观看av| 国产 精品1| 欧美日韩在线观看h| 国产亚洲av片在线观看秒播厂| 免费av中文字幕在线| 免费看光身美女| 性色avwww在线观看| 街头女战士在线观看网站| 国产深夜福利视频在线观看| 精品久久久久久久久av| 欧美人与善性xxx| 一本久久精品| 最新的欧美精品一区二区| 自线自在国产av| 大码成人一级视频| 亚洲人成网站在线观看播放| av卡一久久| 亚洲av国产av综合av卡| 乱码一卡2卡4卡精品| 视频中文字幕在线观看| 欧美日韩一区二区视频在线观看视频在线| 久久99热6这里只有精品| 国产欧美另类精品又又久久亚洲欧美| 国产一区二区三区综合在线观看 | 好男人视频免费观看在线| 日产精品乱码卡一卡2卡三| 涩涩av久久男人的天堂| 噜噜噜噜噜久久久久久91| 免费观看av网站的网址| 男人狂女人下面高潮的视频| 777米奇影视久久| 大香蕉久久网| 秋霞在线观看毛片| 亚州av有码| 性高湖久久久久久久久免费观看| 国产精品久久久久久精品电影小说| 精品国产乱码久久久久久小说| 日韩在线高清观看一区二区三区| 高清欧美精品videossex| 国产一区二区在线观看日韩| 亚洲欧美成人精品一区二区| 最近手机中文字幕大全|