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

    Disordered Translocation is Hastening Local Extinction of the Chinese Giant Salamander

    2021-09-27 11:25:30GuochengSHUPingLIUTianZHAOChengLIYinmengHOUChunlinZHAOJieWANGXiaoxiaoSHUJiangCHANGJianpingJIANGandFengXIE
    Asian Herpetological Research 2021年3期

    Guocheng SHU ,Ping LIU ,Tian ZHAO ,Cheng LI ,Yinmeng HOU,4 ,Chunlin ZHAO,4 ,Jie WANG,Xiaoxiao SHU,4,Jiang CHANG,Jianping JIANG* and Feng XIE*

    1CAS Key Laboratory of Mountain Ecological Restoration and Bioresource Utilization &Ecological Restoration Biodiversity Conservation Key Laboratory of Sichuan Province,Chengdu Institute of Biology,Chinese Academy of Sciences,Chengdu 610041,Sichuan,China

    2Yibin University,Yibin 644000,Sichuan,China

    3University of Chinese Academy of Sciences,Beijing 100049,China

    4Sichuan University,Chengdu 610065,Sichuan,China

    5State Key Laboratory of Environmental Criteria and Risk Assessment,Chinese Research Academy of Environmental Sciences,Beijing 100012,China

    Abstract Biodiversity is declining globally by an unprecedented extinction rate.This is especially true for amphibians,accounting for 24.3% of all threatened vertebrates.As the largest extant amphibian species in the world,wild populations of the Chinese giant salamander (Genus Andrias) (CGS) have decreased dramatically because of overexploitation and habitat degradation.Translocation has become an important strategy for restoring threatened wild populations worldwide.However,disordered translocation usually has negative effects on the native populations.We provide an overview of CGS translocation and show that disordered translocation can increase local population extinction.Nearly four times the estimated number of wild individuals have been released across China.There are three types of translocation used for CGS,namely,reinforcement,reintroduction and ecological replacement,the last of which accounts for over one-third of translocations.Our genetic screening revealed that most released individuals were not from local populations,with one to four lineages detected in every release site (n =6).This disordered translocation can potentially reduce the genetic integrity of original populations.Hence,we suggest suspending current CGS translocation activities immediately,until more robust measures can be developed and implemented to improve the current translocation program,especially with respect to lineage identification and the identification of appropriate release sites.

    Keywords Chinese giant salamander,genetic test,translocation,conservation,wild population

    1.Introduction

    Global biodiversity is in decline and we are currently experiencing unprecedented rates of extinction,for example,current rates of extinction are about 1000 times the likely background rate of extinction and future rates are likely to be 10 000 times higher (Pimm

    et al.

    ,2014).Amphibians are disproportionately the most threatened vertebrate group,with 41% of species considered threatened (Hoffmann

    et al

    .,2010).In the face of rapid biodiversity loss,effective strategies for the conservation of threatened species are urgently required(Bainbridge

    et al

    .,2014;Taylor

    et al

    .,2017).Translocation has become an important tool to manage wildlife over the past few decades,which involves the humanmediated,intentional movements of living organisms from one area to another (Armstrong and Seddon,2007;IUCN/SSC,2013;Batson

    et al

    .,2015;Swan

    et al

    .,2019).If the primary objective is a conservation benefit,the translocation is called a conservation translocation.A conservation translocation is usually undertaken with the goal of improving the conservation status of the target species locally or globally,and/or restoring natural ecosystem functions or processes (IUCN/SSC,2013).Whilst translocations are being used with increased frequency,there is an obvious bias towards mammalian and avian taxa(Seddon

    et al

    .,2005;Bajomi

    et al

    .,2010).The Chinese giant salamander (CGS) have attracted widespread attention because they represent an evolutionary ancient and distinct amphibian lineage and have few close relatives.Also,they are the largest extant amphibian species in the world,with a maximum length and weight of 2 m and 50 kg,respectively (Zhao,1998).Historically,this species was widely distributed in the middle and lower tributaries of the Yangtze River,Yellow River,and Pearl River,including 17 provinces(or municipalities) in central,southern and southwestern China (Wang

    et al

    .,2004).However,the CGS populations have declined dramatically since the 1950s due to habitat loss and overharvesting for food (Wang

    et al

    .,2004).The CGS is listed as a National Special Protected Animal (category II) and as Critically Endangered on the IUCN Red List of Threatened Species (Liang

    et al

    .,2004),Biodiversity Red List of China and in included in CITES Appendix I (Jiang

    et al

    .,2016;CITES,2020).Therefore,the Chinese government has heavily invested to protect and conserve the CGS.So far,there are more than fifty nature reserves associated with protection of CGS in China(Liang

    et al.

    ,2013).Hunting wild CGS is prohibited and breeding the species in captivity requires permission.Recent field surveys have shown that the CGS has been critically depleted or extirpated in most of their original habitats (Xu

    et al

    .,2018;Turvey

    et al

    .,2018).With the development of CGS farming industry,the CGS farms can supply sufficient source populations for translocation programs of CGS,which provides hope and opportunity to save this threatened species.By the end of 2014,there were 2622 CGS farms in China and up to 1 249 000 individuals across the country (China Aquatic Wildlife Conservation,2015).As a result,government-promoted translocations for the protection of the CGS have been conducted widely in China.In 2014,the Ministry of Agriculture and Rural Affairs of the People’s Republic of China (MARAPRC) enacted the technical specification for the stock enhancement of Chinese giant salamander (SC/T 9414-2014).A lack of detailed knowledge on CGS population genetics meant that previous translocations did not consider CGS from different river basins represented as different lineages and the species was treated as a single population.Current studies revealed that the CGS is comprised of at least seven lineages (Yan

    et al

    .,2018;Liang

    et al

    .,2019).Turvey

    et al

    .(2019) argued that the wild populations from Pearl River/Nanling regions should be restored as a valid species

    Andrias sligoi

    (Boulenger,1924).Mix-source reintroductions can lead to outbreeding depression in hybrid offspring or reduce the genetic integrity of the original populations (Huff

    et al

    .,2011).Therefore,a detailed understanding of CGS translocations is important as it will guide future conservation management decisions.

    In this study,we aimed to (a) understand the current CGS translocations including the number of individuals involved and the locations where CGS were released,and the translocation type,(b) to conduct genetic screening on released individuals and assigned them to CGS lineages or species in order to assess the effectiveness and risk of current translocations,and (c) to make recommendations for future CGS translocations.

    2.Materials and Methods

    2.1.Release data collection

    The CGS often grabs the attention of the general public when released because it is regarded as a flagship species in China.Therefore,all CGS translocation news are generally released to the public.Data on CGS translocation were obtained from the internet (mainly media outlets) using the search terms Chinese giant salamander,giant salamander,and baby fish,along with the keywords reintroduction,translocation,release,and stock enhancement.The time frame of the search was from the last century to the end of October 2019.The language of the search was Chinese (Mandarin).We also obtained supplementary data from two fishery bureaus,including the Fishery Bureau of Tongjiang County,Sichuan and the Fishery Bureau of Cili County,Hunan.Following the definitions of IUCN translocations,we classified CGS translocations into 3 different types.Reinforcement,as the intentional movement and release of CGS into these existing populations of conspecifics;reintroduction,as the intentional movement and release of CGS inside their indigenous ranges from which they have disappeared;ecological replacement,as the intentional movement and release of CGS outside their indigenous ranges to perform a specific ecological function.

    2.2.Survey area and methodology

    We conducted sampling surveys at six release sites:Kaihua (4532 released individuals),Zhejiang;Yuexi (1890 released individuals),Anhui;Liannan(400 released individuals),Guangdong;Heyuan (1480 released individuals),Guangdong;Tongjiang (4600 released individuals),Sichuan,and Wanyuan (110 released individuals),Sichuan.These sites represent three distinct distribution regions of Chinese giant salamander in Turvey

    et al.

    (2019):the southwestern(Pearl/Nanling),northern (Yangtze/Sichuan) and southeastern(Huangshan and Zhejiang).Released individuals were collected in natural habitat through visual encounter methods and wire-mesh baited traps in September 2018 and 2019.The visual encounter methods consisted of wading,nocturnal spotlighting,and netting.Wire-mesh baited traps are reasonably effective for surveying

    Andrias japonicus

    and

    Cryptobranchus alleganiensis

    alleganiensis

    (Foster

    et al

    .,2008;Brigger

    et al

    .,2013;Turvey

    et al

    .,2018).We selected survey areas according to habitat conditions during the day and searched for CGSsin the evening (from 18:30 to 23:00).The traps consisted of 110 cm × 60 cm × 30 cm rectangular boxes made of 1.5 cm plastic-coated hardware cloth and a funnel trap.Fresh fish and chicken were used as bait for trapping (Browne

    et al

    .,2011).The traps were set in the upper,middle,and lower reaches of the surveyed river and checked them the next day (between 5:00 and 7:00 AM).

    2.3.Habitat suitability model

    To infer the potential distribution range of the CGS under current climatic conditions,we produced a habitat suitability model according to CGS locality records and environmental data (Chen

    et al

    .,2018).We combined GPS information on the current distribution of the species and environmental variable modeling to project its historical distribution.For distribution information,we obtained sampling records of wild individuals from the literatures (Murphy

    et al

    .,2000;Fei

    et al

    .,2006;Zheng,2006;Li

    et al

    .,2009;Luo

    et al

    .,2009;Guo

    et al

    .,2011;Yang

    et al

    .,2011;Chen

    et al

    .,2017;Wu

    et al

    .,2017;Wang

    et al

    .,2018;Liang

    et al

    .,2019) and the Global Biodiversity Information Facility (GBIF;https://www.gbif.org) and excluded duplicate GPS sites within each pixel (30 arc-seconds,1 km ×1 km).Climatic predictor variables were downloaded from WorldClim (www.worldclim.org/current).We downloaded maps of Chinese administrative areas,elevation,and land use from DIVA-GIS (available at http://www.diva-gis.org/Data).Vegetation cover was extracted from the land use data.Human footprint data were downloaded from http://sedac.ciesin.columbia.edu/wildareas/.To ensure consistency in climate layer resolution,we increased the resolution of last glacial maximum(LGM) layers from 2.5 min to 30 arc-seconds by the ‘resample’method in ArcGIS v.10.3.According to the sampling records,all of China was selected as the range of interest.

    Before constructing the habitat suitability model,bioclimatic raster layers were clipped into the study area and converted to ASC II format for use in MaxEnt.To prevent overfitting,we used the band collection statistic of the spatial analysis tool to calculate the correlation coefficients between bioclimatic layers in ArcGIS v.10.3 and removed the variables with strong correlations (correlation coefficient ≥ 0.8) and low contributions to the habitat suitability model with 19 bioclimatic factors used for preliminary analysis.Twelve bioclimatic predictor variables were ultimately included in our habitat suitability model:elevation,vegetation cover,human footprint data,annual mean temperature,mean diurnal range,isothermality,temperature seasonality,maximum temperature of the warmest month,minimum temperature of the coldest month,annual precipitation,precipitation seasonality,and precipitation of the coldest quarter.

    The habitat suitability model was constructed by MaxEnt v.3.4.1 (Phillips

    et al

    .,2018) with 8 uncorrelated bioclimatic data layers using the maximum entropy method based on the optimal model,with 25% of occurrence records as the test set and 75% as the training set.Other conditions included 10 000 background points,a logistic output format and 10 crossvalidation replicates for each run,with the average of those replicate models representing the final prediction results.The area under the receiver operating characteristic curve (AUC)was used to evaluate the overall fit of the model,and AUC values > 0.9 indicated that the model fit well (Araujo and Guisan,2006);the higher the AUC score is,the stronger the predictive power of the model.Jackknife analysis was employed to identify which bioclimatic factors were most important and contributed the most to the model.

    2.4.Molecular sampling and sequencing

    We swabbed oral mucosa cells or collected exfoliated skin cuticles from caught and trapped individuals.All animals were released after tissue collection at the original site of capture.The collected tissues were immediately stored in 95% ethanol at -20°C for DNA extraction.DNA was extracted using a TIANamp Genomic DNA Kit (Qiagen).Three mitochondrial DNA (mtDNA) fragments from the complete mitochondrial cytochrome b (

    Cytb

    ) gene,the partial control region (D-loop),and the partial cytochrome C oxidase 1 (

    COI

    ) gene were amplified by polymerase chain reaction (PCR).The primers and cycling parameters were as described by literatures (Tao

    et al

    .,2005,2006;Che

    et al

    .,2012).PCR products were purified using a Qiagen PCR purification kit and sequenced on an ABI 3100 automated sequencer.We also downloaded nucleotide sequences available in recent studies (Yan

    et al

    .,2018;Liang

    et al

    .,2019) and other references from GenBank (Table S1) and aligned them with ClustalW built into BioEdit version 7.1.9 (Thompson

    et al

    .,1997) using the default parameters.Alignments were also visually checked for ambiguous alignment in MEGA version 7.0.26 (Tamura

    et al

    .,2013).All new sequences were deposited in GenBank with the accession numbers listed in Table S1.

    2.5.Phylogenetic inference

    Bayesian inference (BI) was used to reconstruct phylogenetic relationships.The optimal partitioning scheme and best-fitting nucleotide substitution model for each partition of the mtDNA genomes were estimated using PartitionFinder v.2.1.1 (Lanfear

    et al

    .,2016).Bayesian analyses were conducted in MrBayes v.3.2.2 (Ronquist

    et al

    .,2012) executed in the CIPRES Science Gateway (http://www.phylo.org/index.php).All searches were performed with two independent runs,each initiated with a random tree.Each run consisted of four Markov chain Monte Carlo (MCMC) chains (three heated chains and one cold chain) estimated for 5 million generations and sampled every 1000 generations.Stationarity was assessed in Tracer 1.6,and the first 25% of samples were discarded as burn-in (Rambaut

    et al

    .,2014).The remaining trees were used to create a 50% majority-rule consensus tree and estimate Bayesian posterior probabilities (BPPs).

    3.Results

    We collected a total of 268 records of CGS translocations from official reports (internet searches,233 records),and fishery stations (35 records),and found that the first releases were undertaken in the Zhangjiajie National Nature Reserve,in Hunan Province,in 2002.By the end of October 2019,a total of 287 840 farmed individuals were released into the wild across 98 counties and cities (Figure 1A,1C) in central and southern China.All records of CGS being released were distrbuted in 16 provinces or municipalities (Figure 1B).In fact,the numbers of released animals may likely be an underestimate as not all releases were reported.Even so,this number is almost four times the estimated wild population size of

    Andrias

    species in China (75 000) (Liang

    et al

    .,2013).The number of released animals increased gradually between 2008 and 2012.This was followed by a sharp increase in the number of released animals between 2013 and 2015.The number of released animals peaked in 2016 and has subsequently declined (Figure 1A).The greatest number of individuals were released in Shaanxi Province since 2002,reaching up to 94 464 in total;by contrast,the minimum number of released individuals occurred in Beijing,where 200 individuals were released (Figure 1B).

    There are three types of CGS translocation according to definitions of IUCN:reinforcement (40.8%),reintroduction(24.5%) and ecological replacement (34.7%) (Figure 1C).This meant more than one-third of releases took place at sites outside the indigenous range of CGS.Our habitat suitability model showed that the suitable habitat of CGS was concentrated in central and southern China and was scattered in the east.However,more than 10% of the release sites were not inside predicted suitable habitat (Figure 1C).

    Figure 1 Translocation overview of the Chinese giant salamander.(A) Trends in the translocations of Chinese giant salamanders from 2002-2019.(B)The number of translocated giant salamanders across different Provinces in China.(C) Translocation sites (triangles) and types of the Chinese giant salamanders.The green areas present the predicted suitable habitat for Chinese giant salamander.The pie chart shows the proportion of three translocations.

    In total,48 CGS individuals were recaptured from five of the six sites we sampled for released animals.The details of sites and number of individuals captured at each site were as followed:Yuexi,Anhui (5 adults,2 larvae);Liannan,Guangdong (3 adults,2 larvae);Kaihua,Zhejiang (28 adults,5 larvae);Heyuan,Guangdong (1 adult,1 larva) and Wanyuan,Sichuan (1 adult).Our molecular data showed that all the animals we captured could be assigned to non-native lineages defined by recent studies (Yan

    et al

    .,2018;Liang

    et al

    .,2019).Released individuals from Kaihua belonged to four non-native clades,namely,B,C,D and U1,described in Yan

    et al

    .(2018) (Figure 2).A nonindigenous lineage from Shaanxi was released in Liannan and Heyuan (B clade;Figure 2).Also,Wanyuan’s animal clustered into a clade with northern Guangxi (U1 clade;Figure 2),and Yuexi’s individuals represent three different non-local clades(U1,U2 and B;Figure 2).

    Figure 2 Bayesian tree based on mtDNA and six sampling release sites of Chinese giant salamander.(A) A simplified Bayesian inference tree based on concatenated mtDNA (Cytb,D-loop and COI).Numbers near branches indicate Bayesian posterior probabilities.The capital letters outside and inside the brackets show corresponding clades for Chinese giant salamander in Yanet al. (2018)and Liang et al.(2019),respectively.Clade D corresponds to A.sligoi in Turvey et al. (2019).The distribution areas for different clades are marked after subulate branches.(B) The sampling release sites of Chinese giant salamander.The numbers represent sampling sites,including Wanyuan,Tongjiang,Yuexi,Kaihua,Liannan,and Heyuan in numerical order.The colored circles correspond to the different clades in the left.

    4.Discussion

    4.1.Factors that may influence CGS translocation

    The success of amphibian translocations are related to multiple factors,such as habitat condition,age,food limitation,disease,the number of animals released,etc.(Griffiths and Pavajeau,2008;Germano and Bishop,2009).Generally,translocations are most successful when more than 1000 individuals were released (Germano and Bishop,2009).On the contrary,most reports indicated that causes of amphibian translocation failure were homing and migration of introduced individuals out of release sites and poor habitat (Germano and Bishop,2009).Current studies show that some CGS translocations have not succeeded in restoring wild CGS populations (Luo

    et al

    .,2009;Lin

    et al

    .,2017;Lu

    et al

    .,2020).Turvey

    et al.

    (2018) reported that dead individuals of CGS were found at one translocation site.Researchers suggested that most translocations of CGS failed to restore the wild populations due to the limitation of factors such as water quality,food supply,human disturbance,release season,individual age,poaching,and subsequent management measures (Luo

    et al

    .,2009;Lin

    et al

    .,2017;Zhang

    et al

    .,2017;Turvey

    et al

    .,2018).Thus,a variety of factors influence the success rate of translocation of CGS.So far,the number of released individuals of CGS (287 840) has far exceeded the number of wild individuals (75 000) (Liang

    et al

    .,2013),but release sites are outside the indigenous range (ecological replacement) and outside the predicted suitable habitat which accounted for more than 33% and 10% of sites,respectively.In general,our habitat suitability model is similar to the previously published model (Chen

    et al

    .,2018),both models predicted that suitable CGS habitat was primarily located in central and southern China.However,there were some differences between our model and the Chen

    et al

    .(2018) model.Our suitable habitats were mainly distributed in Shaanxi,Guizhou,Hunan,Sichuan,Chongqing,Hubei,Henan,and Zhejiang,while the suitable habitats from Chen

    et al

    .(2018)were concentrated in Guizhou,Hunan,Zhejiang,Hubei,Anhui,and Sichuan.This difference could be the result of the additional predictive variables (i.e.human footprint) and the inclusion of more historical records in our model.In fact,the conservation biologists’ original intention was to restore wild CGS population by reinforcement or reintroduction and not to conduct ecological replacement.Currently,there are no good management measures (e.g.monitoring and habitat protection)after releasing CGS back into the wild.Therefore,whether these translocations will succeed or not is completely uncertain and requires mark-recapture studies.In recent years,the scale of CGS translocation has declined since its peak in 2016.This decline might be associated with fluctuations of CGS market value.From 2008 to 2010,the CGS price gradually increased and reached peak at the end of 2010,when the highest price was 3000 RMB/kg (Cheng,2015).During this period the CGS industry expanded and was able to provide a greater number of both larvae and adult animals for translocation.In 2011,driven by huge profits,many farmers started to raise CGS (Cheng,2015),and since 2012,the price of CGS drastically declined (Li

    et al

    .,2018).However,because the sexual maturity cycle of captive takes about 4 years (Luo

    et al

    .,2005).The supply of salamander exceeded the demand in 2016,which was when the CGS translocations reached the largest scale.Meanwhile,the price of CGS has almost bottomed out in 2017 (140 RMB/kg),therefore the CGS farming industry has gone into decline (Li

    et al

    .,2018).Subsequently,the number of CGS translocations also declined.

    4.2.The risk of current CGS translocation

    All the populations of CGS were treated as a single species during the previous translocations,but recent studies revealed that there are at least seven distinct lineages (Yan

    et al

    .,2018;Liang

    et al

    .,2019).Furthermore,a new study revealed that the wild populations from Pearl River/Nanling regions were a valid species

    Andrias sligoi

    (Turvey

    et al

    .,2019).In this study,even though we only recaptured 48 individuals from six release sites(Kaihua,Yuexi,Liannan,Heyuan,Tongjiang,and Wanyuan),all individuals belonged to non-native genetic lineages.For example,the native wild populations from Zhejiang belong to clade E according to the recent results (Yan

    et al

    .,2018;Liang

    et al

    .,2019),however,all sampled released individuals were divided into four non-native genetic lineages (B,C,D and U1) in Kaihua,Zhejiang,where the clade D is the representative of

    Andrias sligoi

    .This is quite contrary the genetic purity of karst cave populations (Wang

    et al

    .,2017).We don’t know if our captures were all the original released individuals or their offspring that were born in the wild.But our results suggest that some released individuals have survived at least four years (e.g.952 individuals were released in Kaihua,Zhejiang in October 2015.)and could compromise the genetic integrity of the native population.Meanwhile,some surveys suggest that the released individuals have established new populations in the wild,such as in Guangdong,but no evidence suggested that interbreeding has occurred between the released and wild populations (Su,2018).However,interbreeding has already occurred in

    Andrias japonicus

    in Japan owing to the introduction of

    Andrias davidianus

    (Wang

    et al

    .,2015),which resulted in severe genetic contamination in the wild.With so many non-native genetic lineages identified at our sampling sites,there is a strong indication that the extensive commercial farming has led to the movement of animals between farms throughout the range of CGS in China (Yan

    et al

    .,2018),and their progeny are being released into the wild without genetic testing.Due to their overwhelming advantage of population size,there is a great possibility that non-native strains or species may eventually eradicate the unique evolutionary characteristics of original populations.Although some hybrids are more fit than their parents and the pursuit of hybrid vigor or “heterosis” has been important in agriculture and animal breeding (Edmands,1999),empirical work indicates that multi-generational outbreeding depression can be sufficiently severe in some cases (e.g.,Houde

    et al

    .,2011;Huff

    et al

    .,2011).In addition,mixing genetically,divergent sources are often proposed to attenuate genetic diversity in reintroduced populations that may result from small effective population sizes (Huff

    et al

    .,2011).However,outbreeding depression is a possible negative tradeoff in hybrid offspring for mixing sources (Huff

    et al

    .,2011).Outbreeding becomes unhelpful with regard to the restoration of the wild population of CGS due to adaptability decline in hybrid offspring.Therefore,the current mixed-source translocation could accelerate the extinction of wild local CGS.

    Besides,the release of non-native lineages may lead to the risk of ecological invasion.In our investigation,more than onethird of release sites were not within the historical distribution of CGS according to the historical records.So,these released individuals should be treated as an alien invasive species at these sites.There may be competitive or predator-prey relationships between released individuals and local aquatic animals and this could destroy the local ecological balance.

    4.3.Translocation suggestions

    Based on our findings,we suggest that the administrative department suspend CGS translocation and collect genetic information for all released populations across the country in order to take immediate remedial measures.Meanwhile,for standardizing the procedures of CGS translocation,we recommend that the following points should be considered:

    Establish dedicated conservation breeding facilities.

    Although CGS farms can be used for translocation programs,commercial breeding activities have resulted in extensive movements of CGS between different farms (Cunningham

    et al

    .,2016).Some evidence suggests that the introgression has occurred in some sampled individuals from different localities(Yan

    et al

    .,2018).This study showed that different non-native lineages were present at the same translocation site,where released animals were all from CGS farms.Furthermore,the stocks held on the private farms may be unsuitable for translocation due to the potential risk of pathogens (e.g.

    Andrias

    davidianus ranavirus

    ,

    Aeromonas hydrophila

    ) and the unknown provenance of the salamanders held on the farms.Therefore,establishing dedicated conservation breeding facilities for CGS of known provenance should be conducted urgently (Turvey

    et al

    .,2018,2019).Current release programs must improve to identify the origin of captive animals and prevent extramural introductions of non-native lineages or different species (Turvey

    et al

    .,2019).

    Pre-preparation for releasing.

    Best practices for CGS translocation needs to involve selection of appropriate release sites (Luo

    et al

    .,2009;Lin

    et al

    .,2017),assessment of habitat (Lin

    et al

    .,2017),selection of animals to be released (Lin

    et al

    .,2017),genetic screening of animals to be released (Lin

    et al

    .,2017),prerelease health screening (SC/T 9414-2014),determination of release scale and time (Luo

    et al

    .,2009),adaptability training for releasing individuals (Gao

    et al

    .,2017;Kenison &Williams,2018),measurement of body parameters,marking with passive integrated transponder (PIT) tags (Dodd,2009;Liu

    et al

    .,2020)and preparation of appropriate food (e.g.small fish,shrimp)(Luo

    et al

    .,2009).Since the 1950s,and especially in the 1980s or human activities have caused damage,loss and reduction of

    Andrias davidianus

    habitat (Zhang

    et al

    .,2002).Our results indicate that some historical sites where CGS were present are no longer suitable for CGS due to the expansion of the human footprint and pollution.Release sites should be chosen according to the current distribution ranges and habitat suitability model,rather than historical distribution records of CGS.After determining the appropriate translocation sites,habitat assessement,including water quality,food supply,human disturbance and predators are required,because these environmental factors can affect the survival,growth and development (Luo

    et al

    .,2009;Zhang

    et al

    .,2017).Released individuals need be to strictly selected,as older individuals have higher survival (Luo

    et al

    .,2009;Zhang

    et al

    .,2016).As mentioned above,all animals destined for release should be subject to genetic screening to avoid genetic contamination.To increase survival,we suggest that the suitable release time in autumn and spring,because the food is plentiful and the water temperature is suitable during this period (Luo

    et al

    .,2009;Lin

    et al

    .,2017).Sex ratios of the released individuals should be between 1:1 and 1:1.5,which is consistent with the sex ratio of natural reproduction of giant salamanders (Luo

    et al

    .,2009).The number of translocated individuals should be determined by the carrying capacity of release site for CGS,because the translocation may fail if the scale of release is too large.The healthy and experienced individuals will likely adapt to new environments quickly.Quarantine and adaptability training are also critical for increasing the survival (Gao

    et al

    .,2017).Adaptability training mainly includes acclimatization to foods commonly found in their proposed release sites,the acclimatization to an array of environmental parameters that individuals will be exposed to at release sites (e.g.running water,temperature and photoperiod,etc.Gao

    et al

    .,2017),and behavioral adaptation (Crane and Mathis,2011;Gao

    et al

    .,2017).Body condition should be assessed for translocated animals,therefore the development and validation of a body scoring system is an important area of future research.Lastly,released individuals could be traced if they are marked with passive integrated transponders (PIT).

    Releasing and post release

    .Releasing cohorts should comprise of individuals of the same size or age class.Continuously track and monitor the released animals.We suggest that the animals should be monitored twice a week after releasing,and the monitoring frequency can be reduced after a month(based on our preliminary translocation project of CGS).The follow-up monitoring should continue for 5 years,since CGSs take at least three or four years to mature (Luo

    et al

    .,2005).Meanwhile,protecting the release sites from habitat destruction and poaching is essential.Any CGS deaths at release sites should be investigated and if possible,mitigated.Finally,all release programs should be systematically evaluated in order to continuously improve the translocations of CGS.Currently there is little evidence that CGS translocation results in the establishment of viable populations of native CGS;we suggest that translocations can only be deemed successful once there is evidence of breeding for several generations and the establishment of self-sustaining populations that have originated from translocated stock.

    Acknowledgements

    We would like to thank an anonymous reviewer and Nicholas Wu for their useful comments and suggestions.We thank all of the field assistants for participating in surveys and sampling.We appreciate the help in collecting data from the Sichuan Provincial Department of Agriculture and Rural Affairs,Gutian Mountain National Nature Reserve,Liannan Chinese Giant Salamander Provincial Nature Reserve,Heyuan Chinese Giant Salamander Provincial Nature Reserve,Heyuan Guishan Chinese Giant Salamander Municipal Nature Reserve.This work was supported by the National Key Programme of Research and Development,Ministry of Science and Technology (2016YFC0503200);the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA23080101);the Biodiversity Survey and Assessment Project of the Ministry of Ecology and Environment,China (2019HJ2096001006);the Open Project of CAS Key Laboratory of Mountain Ecological Restoration and Bioresource Utilization &Ecological Restoration Biodiversity Conservation Key Laboratory of Sichuan Province (KXYSWS2006);the Construction of Basic Conditions Platform of Sichuan Science and Technology Department(2019JDPT0020);and China Biodiversity Observation Networks (Sino BON).

    Appendix

    Table S1 Detailed information of samples and referenced sequences for genetic analysis in this study.Letters in front of brackets are corresponding clades in Yan et al. (2018).The capital letters in brackets present corresponding clades in Liang et al.(2019).

    (Continued Table S1)

    (Continued Table S1)

    91精品国产国语对白视频| 久久久久人妻精品一区果冻| 免费观看av网站的网址| 中文字幕久久专区| 国产精品一二三区在线看| 在线观看美女被高潮喷水网站| 国产精品国产三级专区第一集| 国产精品99久久久久久久久| 直男gayav资源| av卡一久久| 交换朋友夫妻互换小说| 高清视频免费观看一区二区| 日日啪夜夜爽| 麻豆乱淫一区二区| 91精品伊人久久大香线蕉| 亚洲欧美日韩无卡精品| 精品人妻偷拍中文字幕| 成人高潮视频无遮挡免费网站| 观看免费一级毛片| 欧美丝袜亚洲另类| 久久久色成人| 99久久综合免费| 黄色怎么调成土黄色| 亚洲国产欧美人成| 91aial.com中文字幕在线观看| 九九在线视频观看精品| 人妻制服诱惑在线中文字幕| 日韩电影二区| 国产精品久久久久久久久免| 丝袜脚勾引网站| 一二三四中文在线观看免费高清| 久久综合国产亚洲精品| av一本久久久久| 国产男女内射视频| 777米奇影视久久| 亚洲欧美成人综合另类久久久| 一区在线观看完整版| 这个男人来自地球电影免费观看 | 啦啦啦视频在线资源免费观看| 亚洲国产欧美人成| 亚洲综合色惰| 精品国产三级普通话版| 久久久久性生活片| 一区二区三区四区激情视频| 一区二区三区四区激情视频| 天堂8中文在线网| 久久久精品94久久精品| 国产精品国产三级国产av玫瑰| 亚洲综合色惰| www.色视频.com| 国产成人免费观看mmmm| 日韩精品有码人妻一区| 精品久久久久久久久亚洲| 一区二区三区精品91| 日本与韩国留学比较| 欧美日韩在线观看h| 亚洲真实伦在线观看| 国产黄频视频在线观看| av在线观看视频网站免费| 一区二区三区精品91| 欧美bdsm另类| 成人黄色视频免费在线看| 中文欧美无线码| 婷婷色综合大香蕉| 91aial.com中文字幕在线观看| 免费观看a级毛片全部| 亚洲欧美精品专区久久| 欧美亚洲 丝袜 人妻 在线| 久久久午夜欧美精品| 午夜日本视频在线| 我的老师免费观看完整版| 久久久久久久久久久免费av| 欧美激情国产日韩精品一区| 久久这里有精品视频免费| 国产一区有黄有色的免费视频| 日日啪夜夜爽| av在线老鸭窝| 国产精品麻豆人妻色哟哟久久| 大又大粗又爽又黄少妇毛片口| 中文在线观看免费www的网站| 国产91av在线免费观看| 国产69精品久久久久777片| 热re99久久精品国产66热6| 国产乱来视频区| 青春草国产在线视频| 国产美女午夜福利| 国产老妇伦熟女老妇高清| 久久久久精品性色| 亚洲精品日本国产第一区| av卡一久久| 王馨瑶露胸无遮挡在线观看| 日本爱情动作片www.在线观看| 色视频在线一区二区三区| 最近中文字幕高清免费大全6| 啦啦啦视频在线资源免费观看| videossex国产| 亚洲欧美清纯卡通| 精品久久久精品久久久| 国产精品熟女久久久久浪| 国产精品一二三区在线看| 久久精品国产鲁丝片午夜精品| 国产 一区 欧美 日韩| 亚洲不卡免费看| 亚洲精品国产av蜜桃| 女的被弄到高潮叫床怎么办| 女的被弄到高潮叫床怎么办| 国产爽快片一区二区三区| 尤物成人国产欧美一区二区三区| 免费黄频网站在线观看国产| 老师上课跳d突然被开到最大视频| 精品亚洲成国产av| 尤物成人国产欧美一区二区三区| 亚洲人成网站高清观看| 在线观看人妻少妇| 精品亚洲成a人片在线观看 | 亚洲国产高清在线一区二区三| 人妻制服诱惑在线中文字幕| 2021少妇久久久久久久久久久| 亚洲不卡免费看| 十分钟在线观看高清视频www | 亚洲va在线va天堂va国产| 成人18禁高潮啪啪吃奶动态图 | 一级av片app| 内地一区二区视频在线| 精品一区二区三卡| 亚洲性久久影院| 黄色日韩在线| 高清午夜精品一区二区三区| 亚洲av中文av极速乱| 自拍偷自拍亚洲精品老妇| 五月伊人婷婷丁香| 成人国产麻豆网| 一边亲一边摸免费视频| 国产精品国产三级专区第一集| av播播在线观看一区| 国产精品成人在线| av免费观看日本| 国产伦精品一区二区三区视频9| 一级毛片 在线播放| 亚洲一级一片aⅴ在线观看| 欧美精品国产亚洲| 国产精品精品国产色婷婷| 日本爱情动作片www.在线观看| 最近手机中文字幕大全| 国产精品一区www在线观看| 亚洲精品456在线播放app| 99热这里只有精品一区| 夜夜爽夜夜爽视频| 午夜视频国产福利| 国产亚洲av片在线观看秒播厂| 一级爰片在线观看| 欧美高清成人免费视频www| 欧美区成人在线视频| 婷婷色麻豆天堂久久| 精品国产乱码久久久久久小说| 亚洲精品乱久久久久久| 国产亚洲最大av| 91久久精品国产一区二区成人| 一本久久精品| 免费在线观看成人毛片| 特大巨黑吊av在线直播| 夜夜爽夜夜爽视频| 一个人看视频在线观看www免费| 成人免费观看视频高清| 国产高潮美女av| 欧美日韩视频精品一区| 日韩精品有码人妻一区| 亚洲内射少妇av| 女人十人毛片免费观看3o分钟| 国产精品熟女久久久久浪| 秋霞在线观看毛片| av线在线观看网站| 久久久精品免费免费高清| 亚洲中文av在线| 大话2 男鬼变身卡| av播播在线观看一区| 大码成人一级视频| 亚洲精品乱久久久久久| 中文欧美无线码| 丰满人妻一区二区三区视频av| 观看美女的网站| 亚洲,一卡二卡三卡| 丰满少妇做爰视频| 久久99热6这里只有精品| 伊人久久国产一区二区| 久久久久久久久久成人| 女性生殖器流出的白浆| 久久久久久九九精品二区国产| 日韩中字成人| 大片电影免费在线观看免费| 久久6这里有精品| 久久久色成人| 欧美日韩亚洲高清精品| 午夜视频国产福利| 精品酒店卫生间| 国内少妇人妻偷人精品xxx网站| 欧美bdsm另类| 国产亚洲午夜精品一区二区久久| 久久 成人 亚洲| 下体分泌物呈黄色| 三级国产精品片| 联通29元200g的流量卡| 在线播放无遮挡| h日本视频在线播放| 亚洲精品视频女| 最近最新中文字幕免费大全7| 国产精品爽爽va在线观看网站| 精品久久久精品久久久| 妹子高潮喷水视频| 男人添女人高潮全过程视频| 内射极品少妇av片p| 大话2 男鬼变身卡| 视频区图区小说| 久久久精品94久久精品| 成人无遮挡网站| 久久人人爽人人片av| 久久久a久久爽久久v久久| 十八禁网站网址无遮挡 | 国产日韩欧美在线精品| 伊人久久精品亚洲午夜| 国产永久视频网站| 亚洲综合色惰| 一级毛片久久久久久久久女| 国内精品宾馆在线| 欧美丝袜亚洲另类| 在线观看一区二区三区激情| 人妻一区二区av| 极品教师在线视频| 老师上课跳d突然被开到最大视频| 麻豆成人av视频| 纵有疾风起免费观看全集完整版| 99热国产这里只有精品6| 国产男女内射视频| 国语对白做爰xxxⅹ性视频网站| 久久久久性生活片| 精品一区二区三卡| 三级经典国产精品| 久久99蜜桃精品久久| 欧美精品国产亚洲| 亚洲va在线va天堂va国产| 日本欧美国产在线视频| 国产v大片淫在线免费观看| 久久精品熟女亚洲av麻豆精品| 肉色欧美久久久久久久蜜桃| 欧美日韩综合久久久久久| 久久99蜜桃精品久久| 九色成人免费人妻av| 一二三四中文在线观看免费高清| 熟女人妻精品中文字幕| 高清在线视频一区二区三区| 欧美+日韩+精品| 午夜福利网站1000一区二区三区| 亚洲美女黄色视频免费看| 亚洲内射少妇av| 亚洲精品国产av成人精品| 蜜桃亚洲精品一区二区三区| 亚洲精华国产精华液的使用体验| 欧美三级亚洲精品| av免费观看日本| 美女xxoo啪啪120秒动态图| 97精品久久久久久久久久精品| av在线老鸭窝| 免费av不卡在线播放| 男女下面进入的视频免费午夜| 丰满迷人的少妇在线观看| 三级国产精品片| kizo精华| 国内少妇人妻偷人精品xxx网站| 爱豆传媒免费全集在线观看| 亚洲国产最新在线播放| 纯流量卡能插随身wifi吗| 亚洲性久久影院| av在线观看视频网站免费| 夫妻午夜视频| 插阴视频在线观看视频| 午夜精品国产一区二区电影| 十八禁网站网址无遮挡 | 亚洲成人av在线免费| 美女主播在线视频| 久久影院123| 免费黄色在线免费观看| 久久精品熟女亚洲av麻豆精品| 大话2 男鬼变身卡| 我的女老师完整版在线观看| 精品久久久久久电影网| 久久久精品94久久精品| 亚洲欧美中文字幕日韩二区| 观看免费一级毛片| 亚洲成人av在线免费| 久久精品久久久久久噜噜老黄| 女人十人毛片免费观看3o分钟| 晚上一个人看的免费电影| 欧美老熟妇乱子伦牲交| 大话2 男鬼变身卡| 国产欧美另类精品又又久久亚洲欧美| 热99国产精品久久久久久7| 亚洲精品乱码久久久v下载方式| 亚洲在久久综合| 免费人妻精品一区二区三区视频| 一级片'在线观看视频| 国产免费一区二区三区四区乱码| 国产精品伦人一区二区| 久久婷婷青草| 国产乱人视频| 大片电影免费在线观看免费| 久久久久网色| 国产淫语在线视频| 这个男人来自地球电影免费观看 | 大片电影免费在线观看免费| 一本久久精品| 老师上课跳d突然被开到最大视频| 永久网站在线| 夫妻午夜视频| 亚洲人与动物交配视频| 久久久欧美国产精品| 蜜臀久久99精品久久宅男| 青春草亚洲视频在线观看| 免费在线观看成人毛片| 亚洲精品一二三| 亚洲欧美精品专区久久| 一级毛片久久久久久久久女| 久久精品久久久久久久性| 国产欧美日韩精品一区二区| 欧美精品一区二区大全| 日本猛色少妇xxxxx猛交久久| 乱系列少妇在线播放| 国产淫语在线视频| 在线观看av片永久免费下载| 国产精品欧美亚洲77777| 久久精品久久久久久久性| 久久鲁丝午夜福利片| 成人国产麻豆网| 18禁裸乳无遮挡免费网站照片| 欧美一级a爱片免费观看看| 永久免费av网站大全| 国产成人精品婷婷| 国产淫片久久久久久久久| 成人午夜精彩视频在线观看| 国产av精品麻豆| 精品久久久久久久久亚洲| av国产久精品久网站免费入址| 最近最新中文字幕大全电影3| 永久网站在线| 亚洲精品国产成人久久av| 妹子高潮喷水视频| 亚洲色图av天堂| 爱豆传媒免费全集在线观看| 精品一区二区三区视频在线| 国产精品国产av在线观看| 一区二区三区精品91| 亚洲欧美日韩卡通动漫| 最近最新中文字幕大全电影3| 在线看a的网站| 国产精品爽爽va在线观看网站| 欧美3d第一页| 99re6热这里在线精品视频| 亚洲天堂av无毛| 国产免费视频播放在线视频| 极品少妇高潮喷水抽搐| 九九在线视频观看精品| 黑人高潮一二区| 毛片一级片免费看久久久久| 精品亚洲成a人片在线观看 | 直男gayav资源| 两个人的视频大全免费| 一级二级三级毛片免费看| 美女视频免费永久观看网站| 免费大片黄手机在线观看| 视频区图区小说| 18禁在线无遮挡免费观看视频| 狂野欧美白嫩少妇大欣赏| 在线播放无遮挡| 久久精品国产自在天天线| 最后的刺客免费高清国语| 又黄又爽又刺激的免费视频.| 大陆偷拍与自拍| 国产精品国产三级国产av玫瑰| 老司机影院成人| 少妇的逼水好多| 婷婷色综合大香蕉| 能在线免费看毛片的网站| 国产淫语在线视频| 免费人妻精品一区二区三区视频| 女人十人毛片免费观看3o分钟| 久久久久久久亚洲中文字幕| 午夜福利高清视频| 舔av片在线| 欧美成人精品欧美一级黄| 岛国毛片在线播放| 日日撸夜夜添| 水蜜桃什么品种好| 亚洲无线观看免费| 日韩伦理黄色片| 大香蕉97超碰在线| 国产午夜精品一二区理论片| 最黄视频免费看| 亚洲av免费高清在线观看| 人人妻人人看人人澡| 最黄视频免费看| 我要看黄色一级片免费的| 亚洲国产精品999| 国产成人免费观看mmmm| 在线免费十八禁| 麻豆成人av视频| 国产免费一级a男人的天堂| 18禁裸乳无遮挡动漫免费视频| 大香蕉97超碰在线| 国产成人午夜福利电影在线观看| 欧美精品一区二区大全| 高清日韩中文字幕在线| 亚洲欧美一区二区三区国产| 在线亚洲精品国产二区图片欧美 | av在线蜜桃| 欧美一区二区亚洲| 老熟女久久久| 亚洲真实伦在线观看| 国产成人a∨麻豆精品| 又黄又爽又刺激的免费视频.| 日产精品乱码卡一卡2卡三| 国产成人免费无遮挡视频| a级一级毛片免费在线观看| 国内揄拍国产精品人妻在线| 日韩亚洲欧美综合| 午夜视频国产福利| 26uuu在线亚洲综合色| 亚洲四区av| 人人妻人人爽人人添夜夜欢视频 | 看免费成人av毛片| 大码成人一级视频| 色视频www国产| 身体一侧抽搐| 精品少妇久久久久久888优播| 亚洲国产毛片av蜜桃av| 日产精品乱码卡一卡2卡三| 精品久久久精品久久久| 欧美97在线视频| 欧美日韩在线观看h| 一本久久精品| 亚洲欧美成人精品一区二区| 久久久久精品久久久久真实原创| 成年人午夜在线观看视频| 亚洲国产精品999| 日日撸夜夜添| 国产精品一二三区在线看| 人妻夜夜爽99麻豆av| av卡一久久| 国产在线男女| 少妇丰满av| 乱系列少妇在线播放| 国产欧美亚洲国产| 精品久久久久久久久av| 男人狂女人下面高潮的视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品国产色婷婷电影| 中文字幕免费在线视频6| 边亲边吃奶的免费视频| 国产免费视频播放在线视频| 国产又色又爽无遮挡免| 一本色道久久久久久精品综合| 美女高潮的动态| 永久网站在线| 久久国产亚洲av麻豆专区| a级一级毛片免费在线观看| 亚洲国产精品专区欧美| 97在线视频观看| 国产视频首页在线观看| 新久久久久国产一级毛片| 免费人妻精品一区二区三区视频| 中文字幕精品免费在线观看视频 | 成人免费观看视频高清| 国产精品欧美亚洲77777| 97在线人人人人妻| 欧美亚洲 丝袜 人妻 在线| 亚洲国产最新在线播放| 亚洲,欧美,日韩| 欧美精品国产亚洲| 亚洲国产精品国产精品| 精品人妻视频免费看| 乱系列少妇在线播放| 人体艺术视频欧美日本| 青青草视频在线视频观看| a级一级毛片免费在线观看| 少妇丰满av| 国产av精品麻豆| 一级二级三级毛片免费看| 精品一区二区三卡| 九九久久精品国产亚洲av麻豆| 亚洲av综合色区一区| 亚洲欧美精品自产自拍| 久久人人爽av亚洲精品天堂 | 国产精品久久久久久久电影| 亚洲内射少妇av| 精品亚洲乱码少妇综合久久| 国产乱人视频| 日本爱情动作片www.在线观看| 青春草视频在线免费观看| 男女下面进入的视频免费午夜| 狂野欧美白嫩少妇大欣赏| av天堂中文字幕网| 欧美一区二区亚洲| 夜夜爽夜夜爽视频| 一区二区av电影网| 国产乱来视频区| 永久网站在线| 午夜福利视频精品| 欧美激情国产日韩精品一区| 最近最新中文字幕免费大全7| 午夜老司机福利剧场| 岛国毛片在线播放| 国产精品爽爽va在线观看网站| 男人添女人高潮全过程视频| 日日撸夜夜添| 免费黄网站久久成人精品| 黄片无遮挡物在线观看| 一区二区三区免费毛片| 老师上课跳d突然被开到最大视频| 精品国产一区二区三区久久久樱花 | 国产大屁股一区二区在线视频| 在线观看一区二区三区| 久久热精品热| 在线观看三级黄色| 亚洲精品乱码久久久v下载方式| 男人舔奶头视频| 久久久久国产精品人妻一区二区| 大陆偷拍与自拍| av国产久精品久网站免费入址| 国产毛片在线视频| 欧美少妇被猛烈插入视频| 久久国产精品男人的天堂亚洲 | 一级片'在线观看视频| 观看免费一级毛片| av免费在线看不卡| 国产乱人偷精品视频| 3wmmmm亚洲av在线观看| 久热久热在线精品观看| 自拍偷自拍亚洲精品老妇| 国产美女午夜福利| 国产高潮美女av| 日本一二三区视频观看| 国产精品久久久久久久电影| 日本一二三区视频观看| 亚洲av日韩在线播放| 国内精品宾馆在线| 国内少妇人妻偷人精品xxx网站| 丝袜脚勾引网站| 欧美精品一区二区大全| 国产成人aa在线观看| 男女边摸边吃奶| 最近的中文字幕免费完整| 久久99精品国语久久久| 美女视频免费永久观看网站| 精品国产三级普通话版| 亚洲图色成人| 亚洲精品乱久久久久久| 日日啪夜夜爽| 在线观看一区二区三区| 蜜桃亚洲精品一区二区三区| 在线观看av片永久免费下载| 免费观看的影片在线观看| 亚洲国产高清在线一区二区三| 高清欧美精品videossex| 精品国产露脸久久av麻豆| 国产免费一区二区三区四区乱码| 国产成人精品福利久久| 天天躁日日操中文字幕| 色综合色国产| 久久青草综合色| 少妇高潮的动态图| 国产av国产精品国产| 黄色一级大片看看| 在线观看免费高清a一片| 国产男女超爽视频在线观看| 欧美成人一区二区免费高清观看| 色网站视频免费| 性色avwww在线观看| 日本-黄色视频高清免费观看| 亚洲四区av| freevideosex欧美| 18禁动态无遮挡网站| 国产成人免费观看mmmm| 青春草视频在线免费观看| 国国产精品蜜臀av免费| 国产成人一区二区在线| 久久 成人 亚洲| 九九在线视频观看精品| 国产片特级美女逼逼视频| 日韩一本色道免费dvd| 久久人人爽人人片av| 街头女战士在线观看网站| 亚洲人成网站在线播| 一本—道久久a久久精品蜜桃钙片| 久久精品熟女亚洲av麻豆精品| 91久久精品国产一区二区三区| 日日摸夜夜添夜夜添av毛片| 国产黄片视频在线免费观看| 亚洲av二区三区四区| 亚洲国产最新在线播放| 国产熟女欧美一区二区| 夜夜骑夜夜射夜夜干| 男的添女的下面高潮视频| 国产视频内射| 亚洲av成人精品一二三区| 免费少妇av软件| 男女边吃奶边做爰视频| 伊人久久精品亚洲午夜| 午夜老司机福利剧场| 日韩欧美 国产精品| 久久人人爽人人片av| 视频区图区小说| 成人综合一区亚洲| 国产精品爽爽va在线观看网站| 亚洲性久久影院| 在线观看免费高清a一片| 日日摸夜夜添夜夜添av毛片| 日韩大片免费观看网站|