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

    Reversal of Tetracycline Resistance by Cepharanthine,Cinchonidine,Ellagic Acid and Propyl Gallate in a Multi-drug Resistant Escherichia coli

    2021-06-10 10:02:42DarkoJenicHelenWallerHelenCollinsClettErridge
    Natural Products and Bioprospecting 2021年3期

    Darko Jenic·Helen Waller·Helen Collins·Clett Erridge

    Abstract Bacterial resistance to antibiotics is an increasing threat to global healthcare systems.We therefore sought compounds with potential to reverse antibiotic resistance in a clinically relevant multi-drug resistant isolate of Escherichia coli (NCTC 13400).200 natural compounds with a history of either safe oral use in man,or as a component of a traditional herb or medicine,were screened.Four compounds; ellagic acid,propyl gallate,cinchonidine and cepharanthine,lowered the minimum inhibitory concentrations (MICs) of tetracycline,chloramphenicol and tobramycin by up to fourfold,and when combined up to eightfold.These compounds had no impact on the MICs of ampicillin,erythromycin or trimethoprim.Mechanistic studies revealed that while cepharanthine potently suppressed effl ux of the marker Nile red from bacterial cells,the other hit compounds slowed cellular accumulation of this marker,and/or slowed bacterial growth in the absence of antibiotic.Although cepharanthine showed some toxicity in a cultured HEK-293 mammalian cell-line model,the other hit compounds exhibited no toxicity at concentrations where they are active against E.coli NCTC 13400.The results suggest that phytochemicals with capacity to reverse antibiotic resistance may be more common in traditional medicines than previously appreciated,and may off er useful scaff olds for the development of antibiotic-sensitising drugs.

    Keywords Antibiotic resistance ·Natural products ·Phytochemical ·Screening ·Effl ux pump inhibitor

    1 Introduction

    Gram-negative bacteria displaying resistance to multiple classes of antibiotics are increasingly prevalent and difficult to treat in modern healthcare systems [1].Historically,this issue has been countered by the discovery and development of new classes of antibiotics.However,as only two new classes of antibiotic have been discovered since the 1960s,attention has focussed more recently on the identification of agents which target bacterial mechanisms of resistance [2,3].

    This approach has proven successful for the restoration of utility of some β-lactam antibiotics via lactamase inhibitors,such as clavulanic acid [3].Some success has also been reported in the identification of compounds with capacity to reverse resistance to various non-lactam antibiotics in Gram-positive bacteria,particularlyStaphylococcus aureus[4,5].However,resistance reversal in Gram negative bacteria has proven significantly more challenging,and comparatively little progress has been made in identifying novel agents with sensitising properties [6].

    We therefore sought to identify compounds with capacity to reverse resistance to six antibiotics of different classes using a multi-drug-resistantEscherichia coli.Our strategy was based on the hypothesis that since plants have successfully countered the threat of Gram-negative infection for millennia through their production of antibacterial secondary metabolites,it is likely that they have also faced the issue of resistance to these defences via the same resistance mechanisms currently posing a threat to man.This raises the possibility that plant secondary metabolites with potential to combat bacterial resistance mechanisms may have also co-evolved in parallel with compounds with direct bactericidal or bacteriostatic properties,and may therefore have potential as scaffolds for the development of novel classes of therapeutic resistance-reversing drugs.

    As small molecules intended for clinical use should also be oflow toxicity in mammalian systems,and ideally with potential for rapid translation to human trials,we chose to focus only on compounds isolated from plants with a history of use in man as traditional herbs or medicines.A natural product library comprising 200 such phytochemicals was therefore screened to identify compounds with capacity to reverse resistance using the model multi-drug resistantE.coliisolate NCTC 13400.Potential mechanisms for the sensitising properties of the hit compounds,and their impacts on the viability of mammalian cells cultured in vitro,were then explored.

    2 Results and Discussion

    2.1 Primary Screen of 200 Natural Compounds for Reversal of Resistance to Tetracycline

    At least four major mechanisms are currently thought to mediate bacterial resistance to antibiotics; namely drug destruction,efflux pumps,target modification and reduced permeability of the cell wall.Although some success has been achieved in the development of drugs with capacity to reverse resistance via the first of these mechanisms (i.e.β-lactamase inhibitors),comparatively little progress has been made in the discovery of molecules with the capacity to block other forms of resistance,and sensitise to other classes of antibiotic,particularly in Gram-negative bacteria [1,3].

    Plants have faced millennia of exposure to Gram-negative pathogens,many of which are likely to have evolved diverse resistance mechanisms to counter bacteriostatic or bactericidal plant compounds [7,8].As plants may have also faced selection in turn to develop phytochemical secondary metabolites targeting such mechanisms of resistance,we hypothesised that screening a natural product library may offer potential leads for resistance reversing agents.

    The library screened in the present study contains several compounds of natural origin that are used as antibiotics,and three of these (chloramphenicol,doxycycline and rifampicin) fully prevented growth ofE.coliNCTC 13400 at 100 μM,in the presence or absence of 0.25× minimum inhibitory concentration (MIC) tetracycline (Fig.1 a).None of the non-antibiotic compounds tested achieved the same level ofinhibition at the 21 h timepoint.However,a closer inspection of the data from earlier timepoints indicated that several compounds slowed growth to a greater extent in the presence of tetracycline than without,suggesting potential for resistance reversal (Fig.1 b).The four compounds that scored the highest according to this metric (after exclusion of existing antibiotics) were cepharanthine,cinchonidine,ellagic acid and propyl gallate.

    Fig.1 Primary screen of 200 phytochemicals for capacity to reverse tetracycline resistance in E.coli NCTC 13400.E.coli NCTC 13400 was cultured in microtitre plates in the presence or absence of 0.25× MIC (16 μg/mL) tetracycline with 100 μM of each phytochemical,or vehicle control (1% DMSO).Absorbance at 600 nm was measured at 21 h (a) and 5 h (b).Each individual point represents one of the tested phytochemicals.Hit compounds chosen for further study are highlighted

    Checkerboard assays then showed that each of these comp ounds reduced the MIC of tetracycline by up to twofold at both the 8 and 21 h timepoints (Fig.2).We hypothesised that the hit compounds may exhibit synergistic or additive effects on resistance reversal,so a mixture of the 4 phytochemicals together (at 64 μM each) was also tested,and this reduced the MIC of tetracycline by fourfold at both timepoints (Table 1).

    Fig.2 Checkerboard assays of the top four hits from the primary screen.E.coli NCTC 13400 was cultured in the presence of doubling dilutions of tetracycline,in combination with doubling dilutions of cepharanthine (a),cinchonidine (b),ellagic acid (c) propyl gallate (d),or vehicle control (DMSO).Absorbance was measured at 600 nm at 21 h.Means of three independent experiments ± SD are shown

    2.2 Eff ects of Hit Compounds on Sensitivity of E.coli NCTC 13400 to Other Antibiotics

    Table 1 MICs and IC50 of tetracycline in the presence of maximal inhibitory concentrations of each phytochemical

    The four hit compounds had no impact on the MICs of ampicillin,erythromycin or trimethoprim (Fig.3).However,the MIC of chloramphenicol was reduced twofold by cepharanthine and cinchonidine,and fourfold by ellagic acid.Propyl gallate also reduced the MIC of tobramycin by twofold,at 21 h (Table 2).A mixture of the 4 phytochemicals together (at 64 μM each) did not impact on the MICs of ampicillin,erythromycin or trimethoprim,but lowered the MIC of tetracycline by fourfold,chloramphenicol by eightfold and tobramycin by fourfold.Disc diffusion assays using each antibiotic or mixture revealed a similar,although less pronounced,pattern of sensitisation (Table 3).

    Fig.3 Eff ects of top hits on sensitivity of E.coli NCTC 13400 to other antibiotics.E.coli NCTC 13400 was cultured for 21 h in the presence of 64 μM of each phytochemical,vehicle control (0.64% DMSO),and indicated concentrations of ampicillin (a),chloramphenicol (b),erythromycin (c),tobramycin (d),or trimethoprim (e).Absorbance at 600 nm was measured at 21 h.Means of three independent experiments ± SD are shown.Abbreviations:Ellagic acid (Ell),cepharanthine (Cep),cinchonidine (Cin),propyl gallate (Pro),mix of all 4 phytochemicals (Mix4)

    Two of the top four hits from the present study,or their close derivatives,have been reported previously to exhibit antibiotic potentiating effects.Ellagic acid was reported to potentiate the antibiotics novobiocin,coumermycin,chlorobiocin,rifampicin and fusidic acid inAcinetobacter baumannii[9].Gallic acid,which is structurally similar to propyl gallate,was also shown to potentiate the activities of sulfamethoxazole and tetracycline inPseudomonas aeruginosa[10].However,to our knowledge,no sensitising properties of cepharanthine or cinchonidine have been reported previously.

    2.3 Potential Mechanisms of Action of Hit Phytochemicals

    Because the hit compounds lowered the MICs of antibiotics from several different classes,we hypothesised that they may act at least partly as inhibitors of multidrug effl ux systems.To test this,we fi rst measured accumulation of the lipophilic dye Nile red usingE.coliNCTC 13400 cells.Cepharanthine markedly accelerated cellular uptake of this dye,consistent with inhibition of effl ux,while the other hit compounds and mixtures consistently reduced dye accumulation compared to vehicle control (Fig.4 a).Effl ux assays then confi rmed that cepharanthine markedly inhibited effl ux of Nile red from loaded cells,while the other hit compounds and mixtures had little impact on this pathway (Fig.4 b).

    Fig.4 Eff ects of phytochemicals on uptake and effl ux of Nile red in E.coli NCTC 13400.For Nile red accumulation assays (a),E.coli NCTC 13400 was cultured with 64 μM of each phytochemical or vehicle control (0.64% DMSO).After 2 min,5 μM of the marker Nile red was added.For effl ux assays (b),CCCP de-energised cells were pre-loaded with Nile red,and re-energised after 2 min with 50 mM glucose.In both cases,fl uorescence at 620 nm,which is proportional to cellular content of Nile red,was measured every 5 s for indicated periods.The kinetics of bacterial growth in the absence of antibiot-ics were also assayed by measuring absorbance at 600 nm for E.coli NCTC 13400 (c) or E.coli DH5-α (d),in the presence of 64 μM of each phytochemical,or vehicle control (0.64% DMSO).The effects of cepharanthine,artesunate and berberine (64 μM) on efflux were compared with 50 μM CCCP (a positive control for efflux inhibition) in E.coli NCTC 13400 (e).Results are representative of at least three independent experiments.Abbreviations:Ellagic acid (Ell),cepharanthine (Cep),cinchonidine (Cin),propyl gallate (Pro),mix of all 4 phytochemicals (Mix4),mix of Ell,Cin and Pro only (Mix3)

    Table 2 Eff ects of hit phytochemicals on sensitivity of E.coli NCTC 13400 to antibiotics of other classes

    Table 3 Disk diffusion assays of the eff ects of hit phytochemicals on antibiotic sensitivity in E.coli NCTC 13400

    We next tested the hypothesis that the hit compounds may limit cellular metabolism or availability of energy,since these may contribute to resistance mechanisms.In the absence of antibiotics,ellagic acid and propyl gallate both markedly slowed the rate of growth of the test strain,while cepharanthine and cinchonidine had little or no impact on growth rate (Fig.4 c).To test whether the reversal of resistance may be specific toE.coliNCTC 13400,the effects of each phytochemical on antibiotic sensitivity were re-examined using theE.colistrain DH5-α,which lacks pEK499 and resistance to the tested antibiotics.While ellagic acid generally sensitised this strain to most antibiotics to a similar degree to that observed inE.coliNCTC 13400,the other compounds were largely ineffective at reversing resistance to the tested antibiotics in DH5-α (Table 4).In the absence of antibiotics,ellagic acid and propyl gallate had a greater impact on the growth rate ofE.coliDH5-α than cepharanthine and cinchonidine,but the degree ofinhibition was not as pronounced as that seen inE.coliNCTC 13400 (Fig.4 d).

    Taken together,these findings suggest that cepharanthine likely mediates its sensitising effects at least partly via its action as an efflux pump inhibitor (EPI),since this compound both accelerated accumulation of Nile red and slowed its efflux from bacterial cells.Notably,although numerous natural EPIs for Gram-positive bacteria have been described [4,11],very few have been reported for Gram-negative bacteria [2].Two of the best known examples of such molecules are artesunate,which inhibits AcrAB-TolC inE.coli[12],and berberine,which inhibits MexXY-OprM inP.aeruginosa[13].Interestingly,both of these compounds are present in the library,but they did not score well for the primary phenotypic endpoint of tetracycline resistance reversal,ranking in 156th,and 65th places,respectively.Further experiments revealed that neither of these compounds inhibits efflux of Nile red fromE.coliNCTC 13400 (Fig.4 e).Thus,EPI activity against RND type pumps,such as MexXYOprM and AcrAB-TolC,may not necessarily be a strong indicator of potential to reverse resistance to tetracycline.It is also notable that while artesunate is a proven inhibitor of the AcrAB-TolC system inE.coli[12],it did not sensitise to tetracycline in our assays.Taken together,these data suggest that cepharanthine,and the other hits from the present study,are not likely to reverse tetracycline resistance by acting as EPIs for AcrAB-TolC [14].

    Table 4 Eff ects of hit phytochemicals on antibiotic sensitivity of E.coli DH5-α

    It is also interesting to note that cepharanthine was reported previously to exhibit EPI activity in cultured mammalian cells,by inhibiting the human P-glycoprotein multidrug efflux pump [15].This raises the possibility that the target of cepharanthine activity inE.colimay exhibit structural similarity with P-glycoprotein,which is a member of the ATP-binding cassette (ABC)-family of efflux pumps.However,although there are many predicted examples of ABC type pumps in theE.coligenome [16],no such pumps of this class are represented on the pEK-499 plasmid [17].The TetA efflux system,which confers the tetracycline resistance provided by this plasmid,is a member of the structurally distinct major facillitator family of efflux pumps [2].We note thatE.coliDH5-α,which lacks the TetA efflux system,was not sensitised to tetracycline by cepharanthine.However,further studies will be required to ascertain more def niitively whether the TetA system is a target of cepharanthine EPI activity.

    As the other hit compounds did not inhibit Nile red efflux,it is likely that their sensitising properties are conferred via other mechanisms.One possibility is that,as some resistance mechanisms are relatively energy-intensive [18],agents which non-specifically reduce the availability of energy or nutrients could inhibit such mechanisms.Notably,both ellagic acid and propyl gallate reduced the growth rate ofE.coliNCTC 13400 andE.coliDH5-α in the absence of antibiotics.However,further studies will be required to explore this possibility and the mechanism of action of cinchonidine,since this compound had little or no impact on growth in either strain.

    2.4 Eff ects of Hit Compounds on the Viability of Cultured Mammalian Cells

    To gain preliminary insight into the potential impact of the hit compounds on viability of mammalian cells,MTT assays using the HEK-293 (human embryonic kidney) cell line were performed in vitro.Although there was no signifi cant impact of cinchonidine,ellagic acid or propyl gallate on cell viability up to 128 μM,cepharanthine exhibited toxicity with an IC50 of 24.6 μM (Fig.5).The mixture of 4 phytochemicals also showed some toxicity at the highest dose tested,presumably due to content of cepharanthine.A mixture of the 3 non-toxic phytochemicals,lacking cepharanthine,was therefore tested,and this mixture showed no eff ect on viability in this assay at up to 32 μM of each compound.Retesting the 3 compound mixture in theE.colisensitisation assays showed that it had the same eff ect as the 4 compound mixture on the MICs of the tested antibiotics,except for tobramycin,for which sensitivity was increased by only twofold,rather than fourfold for the 4 compound mixture (Table 2).

    Fig.5 Eff ects of hit phytochemicals on viability of human HEK-293 cells cultured in vitro.HEK-293 cells were cultured at 2 × 10 4 cells per well of 96-well plates.After 24 h,cells were treated with indicated concentrations of each phytochemical or DMSO control,and viability was quantifi ed by MTT assay after a further 24 h.Concentrations are given as total phytochemical content for the mixtures.Means of three independent experiments ± SD are shown.Abbreviations:Ellagic acid (Ell),cepharanthine (Cep),cinchonidine (Cin),propyl gallate (Pro),mix of all 4 phytochemicals (Mix4),mix of Ell,Cin and Pro only (Mix3)

    Low toxicity in mammalian systems is a key property of any compound intended for use as a scaff old for drug development.It is promising to note that in our HEK-293 system,ellagic acid,cinchonidine and propyl gallate exhibited no toxicity at concentrations up to 128 μM.This is consistent with the established use of ellagic acid and propyl gallate as approved food additives in various territories.Dietary studies in man have also confirmed that doses of up to 500 mg pure ellagic acid are well tolerated with no obvious side e ffects,and the molecule can reach micromolar concentrations in serum following oral supplementation [19].Cinchonidine is also likely to be of relatively low toxicity,as it is isolated from the bark of the medicinal plantCinchona officinalis,which has been used to treat fever and malaria for centuries.The low toxicity of these hits supports the hypothesised advantage of screening compounds isolated from traditional medicines,and lends further weight to earlier observations that such compounds are likely to exhibit lower toxicity,and greater ease of translation to trials,than purely synthetic compounds,both of which were key objectives of the present study [20].

    Nevertheless,we recognise that cepharanthine exhibited significant toxicity in our HEK-293 system at doses necessary to enhance antibiotic activity.It is possible that this toxicity may relate directly to its activity as an EPI,as numerous attempts have been made to use inhibitors of human multidrug efflux pumps to reverse chemotherapy resistance of tumours in clinical trials.To date,such trials have typically failed due to excessive toxicity of the EPI used,suggesting that inhibition of certain efflux pumps (such as P-glycoprotein,which is a proven target of cepharanthine [1 5]),may be inherently toxic to some human cells [21].

    In summary,the natural compounds ellagic acid,propyl gallate and cinchonidine are identified as potentially useful sensitisers ofE.colicarrying the pEK499 plasmid to antibiotics of several classes.Cepharanthine is also identified as a novel inhibitor of efflux inE.coli.In terms of potential for further discovery,the relatively fruitful ‘hit rate’ of the primary screen (~ 2%),also suggests that phytochemicals isolated from traditional medicines may be more rich in such compounds than previously appreciated.These molecules could have potential as scaffolds for the development of antibiotic resistance reversing drugs with favourable safety profiles.

    3 Experimental Section

    3.1 General Experimental Procedures

    EscherichiacoliNCTC 13400 (Public Health England) was selected for study as it displays resistance to eight different classes of antibiotic,via carriage of the pEK-499 plasmid,which contains ten resistance genes of well-established function [17].The pEK-499 plasmid is also a common cause of antibiotic resistance in urinary tract infections in the UK [17].The control strain,lacking resistance to the tested antibiotics,wasE.coliDH5-α.The phytochemical library used (Puretitre library,Caithness Biotechnologies,UK) comprises 200 compounds of natural origin,and was chosen for screening as it comprises only compounds with a history of either safe oral use in man,or as a component of a traditional herb or medicine.

    3.2 Primary Screen for Phytochemicals with Capacity to Reverse Resistance to Tetracycline

    Tetracycline was used as the model antibiotic in the primary screen,as the focus of this study was on non-lactamase based resistance mechanisms,and because pEK-499 confers a 32-fold increase in minimum inhibitory concentration (MIC) for tetracycline inE.colistrain J53,via expression of the tet(A) gene [17].An overnight culture ofE.coliNCTC 13400 was used to inoculate a test culture to a density of ~ 4 × 10 5 cells/mL in Luria broth (LB) supplemented with 16 μg/mL tetracycline,equivalent to one quarter of the MIC for this antibiotic in this strain.180 μL of this suspension was then plated in 96-well plates,and 20 μL of each natural compound (at 1 mM with 32 μg/mL tetracycline) was added to each well for a fi nal concentration of 100 μM.8 wells per plate were used as vehicle controls (1% dimethyl sulphoxide,DMSO).Plates were then incubated at 37 °C,and absorbance at 600 nm was measured at 0,5,8 and 21 h using a microplate reader (Tecan).All values were corrected for baseline absorbance by subtraction of the value measured at t = 0 to account for the potential impact of test compounds on absorbance at 600 nm.

    Two scores were used to rank the eff ectiveness of the test compounds in the primary screen.The fi rst score was calculated as the ratio of growth rate between 0 and 5 h with tetracycline to growth rate over the same period without tetracycline.The second score was calculated as the ratio of the growth rate between 5 and 8 h in the presence and absence of tetracycline.The sum of the two scores was then used for the fi nal ranking.Existing antibiotics were removed from the list,and the top four remaining hits were taken forward for replication and further studies.

    3.3 Checkerboard Assays and Measurement of Growth Kinetics

    The broth microdilution method was used,with doubling dilutions of compound and each antibiotic ofinterest arranged across a 96-well microplate in a total volume of 150 μL.Plates were incubated at 37 °C,and absorbance at 600 nm was measured at 0,5,8 and 21 h using a microplate reader (Tecan).For measurement of growth kinetics in the absence of antibiotics,bacteria were plated with 64 μM of each phytochemical (or 16 μM each for the mixtures),and cultured at 37 °C with measurements of absorbance at 600 nm every hour to 8 h post plating,and also at 10 and 24 h.Baseline (t = 0) absorbance was subtracted from each value to account for pigments associated with the phytochemicals or mixtures.As the hit phytochemicals alone were unable to prevent growth ofE.colieven at the maximum practicable concentration (1 mM),MICs could not be calculated for these compounds,and therefore fractional inhibitory concentration index (FICI) values could not be determined.

    3.4 Bacterial Accumulation and Effl ux Assays

    Nile red was used as a marker of cellular uptake and effl ux in these experiments,as it is a substrate for theE.coliAcrABTolC effl ux system,its polarity is similar to that of many antibiotics and its entry to and exit from the bacterial cell can be conveniently monitored by changes in fl uorescence.The method of Bohnert et al.was used,with the following modifi cations [14].1.4 mL of an overnight culture ofE.coliNCTC 13400 was pelleted (4000gfor 5 min),then resuspended in phosphate buff ered saline (PBS) with 1 mM MgCl2(PBS-M) to an absorbance at 600 nm of 1.0.For accumulation assays,200 μL of this suspension was plated directly into each well of a 96-well microtitre plate,and phytochemicals were added to a concentration of 64 μM (16 μM of each compound for the mixtures).Fluorescence values for each well were then measured every 5 s for 2 min using excitation of 544 nm and emission of 620 ± 10 nm with a fl uorescence microplate reader (Fluostar Omega,BMG Labtech).20 μL of Nile red stock (1 mM) was then injected into each well to a concentration of 5 μM.Fluorescence values were then measured every 5 s for 40 min.For effl ux assays,bacterial suspensions were prepared in PBS-M as above,but were fi rst supplemented with 10 μM carbonyl cyanide m-chlorophenylhydrazone (CCCP,to de-energise cells and halt effl ux) for 15 min,before addition of Nile red stock to a concentration of 5 μM.After 3 h incubation at room temperature,cells were pelleted (4000gfor 5 min),resuspended in 1.4 mL PBS-M,plated at 200 μL per well,and supplemented with phytochemicals as above.Fluorescence values were measured every 5 s for 2 min.Each well was then supplemented with glucose (to 50 mM) to re-energise the cells and restart effl ux,and fl uorescence measurements were immediately resumed for a further 40 min.

    3.5 Mammalian Cell Viability Assays

    The HEK-293 cell-line was chosen as a model cell-line to investigate the potential toxicity of the hit compounds in mammalian systems because it is human,not of tumour origin (it is virally transformed),and commonly used to test the eff ects of new compounds on cell viability [22].Cells were cultured in Dulbecco’s Modifi ed Eagles Medium (DMEM) with 10% serum and plated at 2 × 104cells per well of 96-well plates.After 24 h of culture,cells were challenged with phytochemicals at concentrations from 0.5 to 128 μM,or vehicle control (equivalent concentration of DMSO).Cell viability was then measured 24 h later by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay,as described previously [23].

    3.6 Statistical Analyses

    One-way ANOVA with Dunnett’s post-hoc test was used to compare the means of test conditions with those of the control condition.IC50values,and assessment of potential signifi cance of differences between dose curves,were calculated using GraphPad Prism.Differences were assumed to be signifi cant at P < 0.05.

    AcknowledgementsThis work was supported by an Erasmus+ Mobility Studentship (ref DK KOBENHA 57) awarded to DJ.The funder was not involved in the interpretation of results or writing of the article.

    Compliance with Ethical Standards

    Conflict ofinterestCE is director of Caithness Biotechnologies Ltd,which provided the library for the study.

    Research Involving Human Participants and/or AnimalsThis article does not contain any studies with human participants or animals performed by any of the authors.

    Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License,which permits use,sharing,adaptation,distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence,and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence,unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use,you will need to obtain permission directly from the copyright holder.To view a copy of this licence,visit http://creat iveco mmons.org/licen ses/by/4.0/.

    嫁个100分男人电影在线观看| 久9热在线精品视频| 他把我摸到了高潮在线观看| 99热网站在线观看| 日本与韩国留学比较| 一个人看的www免费观看视频| 身体一侧抽搐| 小说图片视频综合网站| 麻豆成人av在线观看| 中文资源天堂在线| 亚洲国产精品久久男人天堂| 国产成人aa在线观看| 久久天躁狠狠躁夜夜2o2o| 天天躁日日操中文字幕| 99国产精品一区二区蜜桃av| 91在线精品国自产拍蜜月| 在线看三级毛片| 91麻豆av在线| 国产成年人精品一区二区| 少妇熟女aⅴ在线视频| 日本三级黄在线观看| 久久天躁狠狠躁夜夜2o2o| 丰满乱子伦码专区| 国产高清激情床上av| 国产女主播在线喷水免费视频网站 | 人妻制服诱惑在线中文字幕| 大又大粗又爽又黄少妇毛片口| 内地一区二区视频在线| 非洲黑人性xxxx精品又粗又长| videossex国产| 18禁在线播放成人免费| 琪琪午夜伦伦电影理论片6080| 高清日韩中文字幕在线| 精品人妻视频免费看| 成人欧美大片| 真人一进一出gif抽搐免费| 国产精品久久久久久精品电影| 少妇熟女aⅴ在线视频| 联通29元200g的流量卡| 特大巨黑吊av在线直播| 在线观看舔阴道视频| 99久久久亚洲精品蜜臀av| 久久久成人免费电影| 日韩欧美免费精品| 男人狂女人下面高潮的视频| netflix在线观看网站| 12—13女人毛片做爰片一| 日韩欧美三级三区| 中文在线观看免费www的网站| 搡老熟女国产l中国老女人| av天堂中文字幕网| 国产免费av片在线观看野外av| 亚洲无线在线观看| 免费观看的影片在线观看| 亚洲av不卡在线观看| 色综合婷婷激情| 久久香蕉精品热| 精品国产三级普通话版| 一个人观看的视频www高清免费观看| 久99久视频精品免费| 国产在视频线在精品| 又爽又黄无遮挡网站| 亚洲av二区三区四区| 久久婷婷人人爽人人干人人爱| 在线观看舔阴道视频| 中文字幕精品亚洲无线码一区| 国产一区二区三区av在线 | 亚洲精品亚洲一区二区| 九九久久精品国产亚洲av麻豆| 日本a在线网址| 色视频www国产| 看黄色毛片网站| 赤兔流量卡办理| 欧美xxxx黑人xx丫x性爽| 久久久久久国产a免费观看| 亚洲专区中文字幕在线| 色综合色国产| 国产亚洲精品av在线| 婷婷丁香在线五月| av视频在线观看入口| 国产一区二区亚洲精品在线观看| 久久婷婷人人爽人人干人人爱| 无遮挡黄片免费观看| 综合色av麻豆| 成人二区视频| 国产精品98久久久久久宅男小说| 亚州av有码| 国产亚洲av嫩草精品影院| 别揉我奶头 嗯啊视频| 欧美不卡视频在线免费观看| 久久久久久久久大av| 国产高清有码在线观看视频| 亚洲图色成人| av.在线天堂| 久久亚洲精品不卡| 亚洲欧美清纯卡通| 窝窝影院91人妻| 12—13女人毛片做爰片一| 有码 亚洲区| 简卡轻食公司| 亚洲性久久影院| 免费看a级黄色片| 一级毛片久久久久久久久女| 亚洲男人的天堂狠狠| 国产精品一区二区三区四区免费观看 | 国内精品一区二区在线观看| 久久婷婷人人爽人人干人人爱| 欧美成人一区二区免费高清观看| 可以在线观看毛片的网站| 成人av一区二区三区在线看| or卡值多少钱| 午夜福利成人在线免费观看| 日韩中文字幕欧美一区二区| 国产一区二区亚洲精品在线观看| 精品午夜福利视频在线观看一区| 国产视频一区二区在线看| 美女黄网站色视频| 美女 人体艺术 gogo| 在线国产一区二区在线| 国产一区二区在线观看日韩| 精品免费久久久久久久清纯| 一进一出抽搐动态| 老司机深夜福利视频在线观看| av天堂中文字幕网| 亚洲精品色激情综合| 成年人黄色毛片网站| 啦啦啦韩国在线观看视频| 97碰自拍视频| 亚洲精品影视一区二区三区av| 欧美高清性xxxxhd video| 黄色日韩在线| 99在线人妻在线中文字幕| 可以在线观看的亚洲视频| 18+在线观看网站| 午夜福利18| 亚洲aⅴ乱码一区二区在线播放| 欧美国产日韩亚洲一区| 22中文网久久字幕| 亚洲国产精品成人综合色| 高清在线国产一区| 日韩一本色道免费dvd| 极品教师在线免费播放| 直男gayav资源| 欧美三级亚洲精品| 亚洲av一区综合| 亚洲人成网站在线播| .国产精品久久| 久久精品影院6| 日本撒尿小便嘘嘘汇集6| 伊人久久精品亚洲午夜| 亚洲18禁久久av| 日本一二三区视频观看| 在线免费观看不下载黄p国产 | 亚洲国产日韩欧美精品在线观看| 大又大粗又爽又黄少妇毛片口| 午夜日韩欧美国产| 欧美成人免费av一区二区三区| 九九爱精品视频在线观看| 亚洲人与动物交配视频| 国产视频内射| av天堂中文字幕网| 亚洲性夜色夜夜综合| 色av中文字幕| 12—13女人毛片做爰片一| 老司机午夜福利在线观看视频| 亚洲性久久影院| 男女下面进入的视频免费午夜| 午夜福利视频1000在线观看| 中文字幕高清在线视频| 久久99热这里只有精品18| 久久久色成人| 免费无遮挡裸体视频| 亚洲欧美日韩无卡精品| 桃红色精品国产亚洲av| 色综合婷婷激情| 日韩精品有码人妻一区| 亚洲四区av| 国产亚洲精品综合一区在线观看| 九九久久精品国产亚洲av麻豆| 亚洲最大成人手机在线| 国产精品人妻久久久久久| 亚洲乱码一区二区免费版| 两个人视频免费观看高清| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品亚洲一区二区| 精品久久久久久久末码| 亚洲中文日韩欧美视频| 热99在线观看视频| 在线免费十八禁| 少妇熟女aⅴ在线视频| 亚洲av免费在线观看| 99热6这里只有精品| 欧美激情久久久久久爽电影| 久久久久久国产a免费观看| 欧美极品一区二区三区四区| 日韩欧美 国产精品| 91在线精品国自产拍蜜月| 久久人妻av系列| 国产三级中文精品| 69av精品久久久久久| 成人一区二区视频在线观看| 热99在线观看视频| 亚洲欧美清纯卡通| 一级黄片播放器| 免费电影在线观看免费观看| 日本与韩国留学比较| 亚洲精品成人久久久久久| 婷婷六月久久综合丁香| 精品久久久久久久久av| 丝袜美腿在线中文| 日本黄色片子视频| 91在线观看av| 三级国产精品欧美在线观看| 波多野结衣高清无吗| 男女做爰动态图高潮gif福利片| 露出奶头的视频| 麻豆av噜噜一区二区三区| 免费看光身美女| 国产私拍福利视频在线观看| 男人的好看免费观看在线视频| 乱码一卡2卡4卡精品| 欧美三级亚洲精品| 日日啪夜夜撸| 搡老岳熟女国产| 久久中文看片网| 免费大片18禁| 免费观看的影片在线观看| 51国产日韩欧美| 亚洲美女搞黄在线观看 | 成人二区视频| 国产精品一区二区免费欧美| 成人国产综合亚洲| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站在线播| 中文字幕av成人在线电影| 色播亚洲综合网| 日日摸夜夜添夜夜添小说| 久久婷婷人人爽人人干人人爱| x7x7x7水蜜桃| ponron亚洲| 欧美xxxx黑人xx丫x性爽| 免费av观看视频| 香蕉av资源在线| 夜夜夜夜夜久久久久| 熟女电影av网| 中亚洲国语对白在线视频| 人人妻人人看人人澡| 日韩欧美一区二区三区在线观看| 夜夜夜夜夜久久久久| 国内精品美女久久久久久| 成人午夜高清在线视频| netflix在线观看网站| 亚洲av免费在线观看| 91麻豆av在线| 亚洲精品在线观看二区| 麻豆成人午夜福利视频| 在线天堂最新版资源| 最近中文字幕高清免费大全6 | 国内揄拍国产精品人妻在线| 色视频www国产| 国内精品美女久久久久久| 99久久无色码亚洲精品果冻| 欧美人与善性xxx| 国产男人的电影天堂91| 久久久久精品国产欧美久久久| 国产三级中文精品| 乱系列少妇在线播放| 国产精品98久久久久久宅男小说| 亚洲欧美日韩卡通动漫| 搞女人的毛片| 美女xxoo啪啪120秒动态图| 精品久久久久久久久久免费视频| 国产 一区精品| 别揉我奶头~嗯~啊~动态视频| 一个人看的www免费观看视频| 黄色视频,在线免费观看| 老女人水多毛片| 国产av在哪里看| 国内精品宾馆在线| av.在线天堂| 成人午夜高清在线视频| 特级一级黄色大片| 深爱激情五月婷婷| 国产伦精品一区二区三区四那| 一a级毛片在线观看| 春色校园在线视频观看| 欧美丝袜亚洲另类 | 亚洲国产高清在线一区二区三| 国产高清激情床上av| 一个人观看的视频www高清免费观看| 亚洲男人的天堂狠狠| 国产一区二区三区视频了| 免费电影在线观看免费观看| 成人午夜高清在线视频| 在线看三级毛片| 久久精品国产亚洲av涩爱 | a级一级毛片免费在线观看| 九九久久精品国产亚洲av麻豆| 亚洲精品456在线播放app | 一级av片app| 国产爱豆传媒在线观看| 成年人黄色毛片网站| av福利片在线观看| 国产精品久久视频播放| 黄色配什么色好看| 欧美性猛交╳xxx乱大交人| 可以在线观看毛片的网站| 小蜜桃在线观看免费完整版高清| 久久精品国产亚洲av香蕉五月| 中出人妻视频一区二区| 男人舔奶头视频| 成年人黄色毛片网站| 在线观看66精品国产| 日本爱情动作片www.在线观看 | 精品不卡国产一区二区三区| 97热精品久久久久久| 成人永久免费在线观看视频| 乱人视频在线观看| 成人一区二区视频在线观看| 久久久久久久久久黄片| 在线播放国产精品三级| 国产伦精品一区二区三区视频9| 免费看av在线观看网站| 国产精品久久久久久亚洲av鲁大| 在线免费十八禁| 亚洲成人免费电影在线观看| 99国产精品一区二区蜜桃av| 国产伦在线观看视频一区| 日韩欧美三级三区| 国产精品伦人一区二区| 国产亚洲精品av在线| 又爽又黄无遮挡网站| 欧美不卡视频在线免费观看| 久久精品久久久久久噜噜老黄 | 波多野结衣高清作品| 精品一区二区三区视频在线| 亚洲成av人片在线播放无| 精品乱码久久久久久99久播| 国产 一区 欧美 日韩| 日本 av在线| 蜜桃亚洲精品一区二区三区| 人人妻,人人澡人人爽秒播| 天天躁日日操中文字幕| 亚洲av二区三区四区| 亚洲国产精品成人综合色| 亚洲va日本ⅴa欧美va伊人久久| 最近视频中文字幕2019在线8| 午夜精品一区二区三区免费看| 高清日韩中文字幕在线| 日韩欧美国产在线观看| 高清日韩中文字幕在线| 精品一区二区三区av网在线观看| 国产成年人精品一区二区| av黄色大香蕉| 九九久久精品国产亚洲av麻豆| 九九爱精品视频在线观看| 久久精品久久久久久噜噜老黄 | 在线免费观看的www视频| 午夜久久久久精精品| 国产精品av视频在线免费观看| 性欧美人与动物交配| 国产一区二区三区av在线 | 欧美xxxx性猛交bbbb| 日本黄色视频三级网站网址| 真人一进一出gif抽搐免费| 亚洲一区二区三区色噜噜| 高清毛片免费观看视频网站| 欧美日韩黄片免| 嫁个100分男人电影在线观看| 国产视频一区二区在线看| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线播放欧美日韩| 日本欧美国产在线视频| 免费在线观看影片大全网站| 在线免费十八禁| ponron亚洲| 国产成人aa在线观看| 大型黄色视频在线免费观看| 两人在一起打扑克的视频| 亚洲精品在线观看二区| 亚洲精品456在线播放app | 淫妇啪啪啪对白视频| 我的女老师完整版在线观看| 色尼玛亚洲综合影院| 国产精品av视频在线免费观看| 久久精品国产亚洲av天美| 久久亚洲精品不卡| 国产 一区 欧美 日韩| 97超视频在线观看视频| 可以在线观看毛片的网站| 琪琪午夜伦伦电影理论片6080| 亚洲中文日韩欧美视频| 在线观看美女被高潮喷水网站| av在线亚洲专区| 欧美一区二区亚洲| 亚洲欧美日韩东京热| 欧美成人性av电影在线观看| 看十八女毛片水多多多| 日韩欧美精品免费久久| 国产午夜福利久久久久久| 毛片女人毛片| 国产高清视频在线播放一区| 午夜激情福利司机影院| 97人妻精品一区二区三区麻豆| 人妻久久中文字幕网| 亚洲美女视频黄频| 国产又黄又爽又无遮挡在线| 亚洲国产欧洲综合997久久,| 午夜日韩欧美国产| 日日干狠狠操夜夜爽| 国产黄片美女视频| 成年人黄色毛片网站| 午夜福利高清视频| aaaaa片日本免费| 成人国产一区最新在线观看| 国产亚洲欧美98| 精华霜和精华液先用哪个| 女的被弄到高潮叫床怎么办 | 日韩一区二区视频免费看| 琪琪午夜伦伦电影理论片6080| 亚洲成av人片在线播放无| 午夜亚洲福利在线播放| 亚洲专区国产一区二区| 嫁个100分男人电影在线观看| av在线蜜桃| 简卡轻食公司| 亚洲美女黄片视频| 一卡2卡三卡四卡精品乱码亚洲| 国产主播在线观看一区二区| 成人美女网站在线观看视频| 色综合站精品国产| 国产精品乱码一区二三区的特点| 在线观看午夜福利视频| 成人毛片a级毛片在线播放| 亚洲av免费高清在线观看| 国产精品一区二区性色av| 三级毛片av免费| av在线观看视频网站免费| 如何舔出高潮| 色吧在线观看| 窝窝影院91人妻| 欧美性猛交╳xxx乱大交人| 成人高潮视频无遮挡免费网站| 亚洲七黄色美女视频| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩卡通动漫| 日韩在线高清观看一区二区三区 | 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品色激情综合| 男女那种视频在线观看| www.色视频.com| 一卡2卡三卡四卡精品乱码亚洲| 国产三级在线视频| 禁无遮挡网站| 亚洲国产日韩欧美精品在线观看| 亚洲av免费在线观看| bbb黄色大片| 国语自产精品视频在线第100页| 国产精品三级大全| 日韩欧美国产在线观看| xxxwww97欧美| 偷拍熟女少妇极品色| 日日夜夜操网爽| 国产私拍福利视频在线观看| 99热这里只有精品一区| 日韩欧美精品v在线| 两人在一起打扑克的视频| 午夜a级毛片| 99热这里只有是精品50| 国产高清视频在线观看网站| 91狼人影院| 亚洲av一区综合| 成年版毛片免费区| 亚洲欧美日韩高清专用| 自拍偷自拍亚洲精品老妇| 亚洲精品乱码久久久v下载方式| 日本精品一区二区三区蜜桃| 色5月婷婷丁香| 亚洲aⅴ乱码一区二区在线播放| 成人特级黄色片久久久久久久| 欧美xxxx性猛交bbbb| av中文乱码字幕在线| 久久99热6这里只有精品| 亚洲 国产 在线| 色精品久久人妻99蜜桃| 69av精品久久久久久| 联通29元200g的流量卡| 欧美黑人欧美精品刺激| 蜜桃久久精品国产亚洲av| 久久人人精品亚洲av| 国产中年淑女户外野战色| 99久久久亚洲精品蜜臀av| 国产精品嫩草影院av在线观看 | 97碰自拍视频| 精品久久久久久成人av| 国产精品国产三级国产av玫瑰| 中文字幕高清在线视频| 一个人免费在线观看电影| 亚洲色图av天堂| 桃红色精品国产亚洲av| av在线观看视频网站免费| 精品99又大又爽又粗少妇毛片 | 国产 一区 欧美 日韩| 91久久精品国产一区二区三区| 老熟妇乱子伦视频在线观看| 成人av在线播放网站| 97人妻精品一区二区三区麻豆| 九色国产91popny在线| 丰满乱子伦码专区| 婷婷精品国产亚洲av| 欧美色视频一区免费| 亚洲国产精品sss在线观看| 国产大屁股一区二区在线视频| 成年女人永久免费观看视频| 欧美日韩瑟瑟在线播放| 日本欧美国产在线视频| 91av网一区二区| 亚洲在线自拍视频| 男女之事视频高清在线观看| 欧美日韩黄片免| 久久久国产成人精品二区| 国产精品国产三级国产av玫瑰| 中文字幕av在线有码专区| 性插视频无遮挡在线免费观看| 1000部很黄的大片| 啦啦啦啦在线视频资源| 亚州av有码| 国产伦一二天堂av在线观看| 日本一二三区视频观看| 亚洲欧美日韩无卡精品| av在线老鸭窝| 一个人免费在线观看电影| 国产午夜福利久久久久久| 国产视频内射| 成人国产麻豆网| 级片在线观看| 久久精品综合一区二区三区| 白带黄色成豆腐渣| 国产麻豆成人av免费视频| 搡女人真爽免费视频火全软件 | 成人永久免费在线观看视频| 色av中文字幕| 国产精品三级大全| 一个人看的www免费观看视频| 日韩 亚洲 欧美在线| 国产精品av视频在线免费观看| 赤兔流量卡办理| 欧美xxxx黑人xx丫x性爽| 人妻夜夜爽99麻豆av| 蜜桃亚洲精品一区二区三区| 久久精品影院6| 婷婷亚洲欧美| 久99久视频精品免费| 草草在线视频免费看| 国产伦精品一区二区三区四那| 亚洲aⅴ乱码一区二区在线播放| 成人毛片a级毛片在线播放| 女人被狂操c到高潮| 97超视频在线观看视频| 老司机深夜福利视频在线观看| 久久国内精品自在自线图片| 亚洲成人中文字幕在线播放| 可以在线观看毛片的网站| 一个人观看的视频www高清免费观看| 男女啪啪激烈高潮av片| av黄色大香蕉| 免费观看人在逋| 亚洲第一区二区三区不卡| 18禁黄网站禁片午夜丰满| 精品国产三级普通话版| 亚洲三级黄色毛片| 日韩欧美免费精品| 国产女主播在线喷水免费视频网站 | 成人二区视频| 国产国拍精品亚洲av在线观看| 亚洲国产精品合色在线| 热99re8久久精品国产| 我的女老师完整版在线观看| 国产探花极品一区二区| 熟女人妻精品中文字幕| 搡女人真爽免费视频火全软件 | 国语自产精品视频在线第100页| 午夜激情福利司机影院| 国产精品国产三级国产av玫瑰| 日本熟妇午夜| 99精品久久久久人妻精品| 国产免费av片在线观看野外av| 日日摸夜夜添夜夜添av毛片 | 午夜福利在线在线| 乱码一卡2卡4卡精品| 日本爱情动作片www.在线观看 | a级毛片a级免费在线| 小说图片视频综合网站| 亚洲精品456在线播放app | 国产精品国产高清国产av| 少妇人妻精品综合一区二区 | 久久精品91蜜桃| 人妻久久中文字幕网| 日韩强制内射视频| 中文亚洲av片在线观看爽| 免费看av在线观看网站| 亚洲自偷自拍三级| 亚洲性久久影院| 日韩高清综合在线| 1000部很黄的大片| 亚洲欧美激情综合另类| 日本 av在线| 免费人成在线观看视频色| 一区二区三区免费毛片| 国产 一区精品| 波多野结衣巨乳人妻| 国产精华一区二区三区| www日本黄色视频网| 欧美高清成人免费视频www| 一本一本综合久久|