碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx

上传人:李司机 文档编号:7055702 上传时间:2024-06-05 格式:DOCX 页数:14 大小:161KB
返回 下载 相关 举报
碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx_第1页
第1页 / 共14页
碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx_第2页
第2页 / 共14页
碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx_第3页
第3页 / 共14页
碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx_第4页
第4页 / 共14页
碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx_第5页
第5页 / 共14页
点击查看更多>>
资源描述

《碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx》由会员分享,可在线阅读,更多相关《碳青霉烯类耐药性在人源和动物源产NDM酶大肠埃希菌间的传播.docx(14页珍藏版)》请在三一办公上搜索。

1、一、引言抗生素耐药性是21世纪全球最大的公共卫生挑战之一,碳青霉烯类药物是我们治疗危及生命的多重耐药菌感染的最有效抗生素之一。不幸的是,碳青霉烯耐药肠杆菌(CRE)的快速传播是目前困扰住院患者和医生的主要问题,因为CRE对碳青霉烯类药物、头抱菌素类药物和青霉素类药物均表现出高水平耐药性,严重限制了治疗方案。因此,CRE被称为“梦魇细菌,同时伴随着大量的经济负担,并被列为首要的公共健康威胁之一以及迫切需要优先研发新抗生素的病原体。目前,世界范围内医疗机构相关的CRE菌株主要产三种碳青霉烯酶:新德里金属-6内酰胺酶(NDM)、肺炎克雷伯菌碳青霉烯酶(KPC)和可水解碳青霉烯的氧青霉烯酶48型6-内

2、酰胺酶(OXA-48)产NDM菌株在50多个国家广泛流行,特别是在南亚(印度、巴基斯坦和孟加拉国)和中国,而产KPC菌株更普遍在北美洲和南美洲(美国、阿根廷和巴西)以及在中国和欧洲(希腊、意大利和以色列)流行。一项涉及36个欧洲国家的跨国研究显示,医院获得性碳青霉烯耐药大肠杆菌(CREC)中NDM和KPC阳性率分别为103%和7.2%。在中国,2014-2015年在25个省份收集的产NDM和KPC临床分离株分别占所有CREC的49.4%和39.6%。相反,在另一项中国19个省份健康人群的研究中,46.7%的CREC中blaNDM是唯一检测到的碳青霉烯酶基因。综上所述,这些数据表明,NDM是中国

3、大多数人源CREC菌株的碳青霉烯酶。在世界范围内来自食品、野生和伴侣动物的CREC菌株中,NDM仍然是流行最广的碳青霉烯酶。在中国,产NDM菌株正成为食品动物中最主要的产碳青霉烯酶菌株,并污染下游的肉制品,其中,boA一48和Hokpc极少出现。然而,尽管有一些证据表明CREC可在人与动物之间直接传播,但迄今为止,动物来源的CREC还没有被普遍认为是影响人类健康的问题之一。二、材料和方法(一)研究设计图1展示了本项研究的框架。文章由回顾性横断面研究开始。首先,从中国细菌耐药监测网(CHINET)中,获得在2016年1月1日至12月31日期间,中国30家医院的门诊患者和住院患者中收集的共计297

4、99株大肠杆菌分离株。这些医院多数为大型诊疗机构,共覆盖22个省(直辖市),所在省(直辖市)总人口约9亿。为了避免重复,根据患者识别代码,临床标本中每个患者仅包含一个分离株。物种鉴定采用自动化系统,如Vitek2Compact(bioMerieux,法国)和PhoeniX-100(BectonzDickinsonandCompany,美国)系统或基质辅助激光解吸/电离飞行时间质谱(MALDI-TOfMS;BrukerDaltonics,德国)。使用Vitek2Compact自动系统(bioM4rieux)测定分离株对19种临床抗菌药物(头抱他噬、头抱嚷肠、头抱毗后、头胞哌酮/舒巴坦、哌拉西林/

5、他嚏巴坦、氨曲南、厄他培南、亚胺培南、美罗培南、阿米卡星、庆大霉素、环丙沙星、左氧氟沙星、甲氧茉咤/磺胺甲嗖、磷霉素、黏菌素、吠喃妥因和替加环素)的敏感性。药敏结果根据CLSl(CliniCalandLabOratOryStandardSInstitute,美国临床和实验室标准协会)标准进行判读,至少对一种碳青霉烯类药物具有抗性的菌株被定义为CREC(图1)。Fnrtpert29799S CO*,29 IMCSEC 1QUkatx31CC*dorecatonNCBl E M 2 31SNM-trmn18M) 043) ( from CIWM) Hmn(rt15) OcMn(n123 S曲 y(

6、4S)1 sqnc type e)s (osnga 1 n 6mM warns 2. Sxxo trg rMrfyws 31:1)Xbcepctate from 22 p Dec 316 2016)Anthropogenicdato图1.数据来源及研究流程图CMIC:最小抑菌浓度;CSEC:碳青霉烯敏感大肠杆菌;POS:阳性:neg:阴性。其次,于2017年4月1日至5月31日间在辽宁、湖南和陕西三个省各收集22个猪养殖场和16个鸡养殖场的粪便样本,以评估Hgndm和Rag的总体丰度;关于样本采集和处理过程的细节己在之前的研究中描述,如果提取的DNA质量不足以进行下一步分析,则排除该养殖场数据

7、(图1)。如前所述,使用qPCR方法以16SrRNA基因为对照,评估RaNDM和Raw的相对丰度。再次,对全球的动物源和人源CREC菌株进行了全面的基因组研究,以确定它们在基因、质粒和菌株水平上的基因组特征和关联(图1)。详细的过程将在下文中描述.本研究于2015年11月24日获得复旦大学附属华山医院伦理委员会的道德许可(项目号:KY2015-285)o个人同意书采用普通话版本,并从所有住院患者那里通过面对面或电话交流获得同意。研究排除了不愿意参与本研究的个体,并且所有参与者都被告知他们有权在任何阶段退出本研究。(二)数据来源中国22个省(直辖市)的人口统计数据来自中国国家统计局和中国农业年鉴

8、。人为因素共分为6类,包括:国内生产总值(GDP)、人口(数量和密度)、抗菌药物使用、动物生产(食品动物和淡水养殖,主要是鱼和虾)、动物源性食品消费和每日动物源性食品摄入量。大部分数据收集于2015年,抗菌药物使用数据来自2013年。在基因组分析中,我们从NCBI数据库中下载了所有可用的中国NDM阳性大肠杆菌全基因组序列(=463)用于综合分析。我们还从NCBl数据库中获得了来自4个宿主(人、猪、鸡、蝇)的1852株NDM阴性大肠杆菌的基因组序列,相当于NDM阳性序列总数的4倍(图1)。所有菌株均给予有序自然数,并使用EXCel中的RAND函数随机选择进行相应分析。(三)基因组分析基因组序列分

9、析包括基因组组装、多位点序列分型、质粒检测和系统发育种群预测,具体描述见已发表的文章。使用Bandage(0.8.1版)对不含质粒免制子基因的序列片段进行分析,以确定可能的质粒类型。桑基图由R3.6.1中的networkD3软件包(0.4版本)生成。所有组装后的NDM阳性基因组序列用于核心基因组比对,通过RedDogpipeline生成系统发育树。该树使用TreeofLifev5进行可视化处理。(四)溯源分析我们使用DAPC(discriminantanalysisofprincipalcomponents)模型来预测所有NDM阳性大肠杆菌的潜在来源。所有菌株(n=2315)的核心基因组单核甘

10、酸多态性(SNP)图谱,包括NDM阴性的大肠杆菌,均使用snippypipeline生成。用于DAPC模型的1389株NDM阴性大肠杆菌的SNP矩阵通过R3.6.1中的Adegenet软件包实现。其余463株阴性菌株用于评估上述模型的拟合度。最后,利用构建的DAPC模型预测了所有中国NDM阳性菌株(n=463)的遗传来源。预测组:训练组:测试组=1:3:1。(五)统计分析将人口统计学和非临床数据输入EXCeI2016,将变量“年龄分为7个亚组(婴儿、学龄前儿童、学龄儿童、青少年、青年、成人和老年人),同时将变量感染部位分为5个亚组(血液、尿液、痰液、脑脊液和其他)。此外,基于之前对鸡生产和国家

11、猪生产计划的研究,我们根据年出栏量将中国大陆分别划分为2个鸡生产区和4个猪生产区。分类数据采用X2分析检验,连续数据采用Mann-WhitneyU检验。在SPSS23.0版本(IBM,美国)中进行单变量分析,随后使用Cramer系数(0)对P值0.05的变量进行共线性评估。如果一对变量高度相关(04.O),我们选择更具生物逻辑合理性的变量进行多变量逻辑回归分析。单变量分析中具有显著性意义的变量将被用于多变量分析。采用向前逐步(似然比)方法对多变量逻辑回归分析中的混杂因素进行控制。最终模型中保留P值小于005的变量。逻辑回归模型的拟合优度采用Hosmer-Lemeshow检验。森林图使用R语言中

12、的forestplot包生成。使用R3.6l对bNDM的丰度进行分析和可视化。当丰度值低于检测限时,添加一个伪数值IXIO-7(低于最低丰度一个数量级)后再进行IoglO变换。不同组间的比较采用Mann-WhitneyU检验或Kruskal-WallisH检验,取0.05为显著性水平。三、结果(一)CREC与临床和非临床因素的相关性本研究设计示意图如图1所示。在来自中国22个省(直辖市)的不同感染样本中共分离到29799株大肠杆菌,其中,631株(2.1%)被鉴定为CREG四川省的CREC流行率最高(4.7%),福建省无CREC分离株。对CREC感染患者的分析显示,男女比例为L4:1,但根据单

13、变量分析,男性CREC的感染率明显高于女性(P0Q5表明所选的模型适合当前的数据。图2.单变量(a)和多变量(b)分析的森林图。(一)敏感性试验药敏试验表明,CREC(n=631)对头抱菌素、头抱哌酮/舒巴坦、哌拉西林/他嚏巴坦和氟噗诺酮类(58.8%95.7%)的耐药性明显高于碳青霉烯敏感大肠杆菌(Carbapenem-SusceptibleE.colitCSEC;n=29168,2.9%59%)(P0.0001)(表1)。CREC对阿米卡星和磷霉素的耐药率分别为16.8%和28.5%,而CSEC对阿米卡星和磷霉素的耐药率分别为2.7%和5.6%。CREC和CSEC对黏菌素和替加环素的敏感性

14、均较高(大于96.7%)。值得注意的是,59%的CSEC对头抱口塞后表现出耐药性(表1)。虽然无法获得CHINET分离株的耐药基因型,但可利用NCBI数据库中中国临床大肠杆菌分离株(NDM阳性菌株和阴性菌株分别为195株和494株)的序列检测两组间的耐药基因阳性率。总的来说,除和OqXA8基因外,NDM阳性菌中大多数抗菌药物的耐药基因,特别是6-内酰胺耐药基因(Hacrx-MbGoA和b0cM)的阳性率显著高于阴性菌株。最近上传至NCBl的不携带NDM的mcr阳性菌株数量的增加可能是导致NDM阴性菌株中mcr基因阳性率更高的原因;该结果与我们的药敏试验结果相一致一一NDM阴性菌株对有菌素的耐药

15、性高于阳性菌株。表1临床CREC和CSEC分离株的药敏数据AntimicrobialagentsCREC(n=631)CSEC(n=29168)PvalueCefoperazone/sulbactam413(65.5%)1313(4.5%)0.0001Piperacillin/tazobactam371(58.8%)846(2.9%)0.0001Ceftazidime536(84.9%)7380(25.3%)0.0001Cefotaxime604(95.7%)17209(59.0%)0.0001Cefepime483(76.5%)7671(26.3%)0.0001Cefoxitin413(65

16、.5%)3413(11.7%)0.0001Aztreonam481(76.2%)10150(34.8%)0.0001Ertapenem569(90.2%)00.0001Imipenem395(62.6%)00.0001Meropenem487(77.2%)00.0001Amikacin106(16.8%)788(2.7%)0.0001Gentamicin339(53.7%)12367(42.4%)0.0001Ciprofloxacin482(76.4%)16247(55.7%)0.0001Levofloxacin454(71.9%)15459(53.0%)0.0001Trimethoprims

17、ulfamethoazole414(65.6%)16392(56.2%)0.0001Fosfomycin180(28.5%)1633(5.6%)0.0001Colistin11(1.7%)700(2.4%)0.2849Nitrofurantoin107(17.0%)1108(3.8%)0.0001Tigecycline21(3.3%)117(0.4%)11(n=90)和10(n=96)种STS在人源NDM阳性大肠杆菌与鸡、猪和苍蝇源NDM阳性大肠杆菌中同时存在,而在人源NDM阴性菌中,分别有25(=93)和19(=77)种STS与鸡和猪的相关STS重叠。在所有4个来源中均检出了属于STI67、

18、ST206、STlO和ST48的NDM阳性菌株。其中,85.9%的(55/64)STI67和62.1%(18/29)的STlO分别来自人和猪。(五)动物源和人源NDM阳性菌株/质粒之间的基因组关联研究利用463株NDM阳性大肠杆菌间的362650个SNPs构建了核心基因组系统发育树,研究了菌株与质粒水平上的关联性(图3和图4).通过种群结构贝叶斯分析(BAPS),将4个来源的菌株划分为11个谱系。人源菌株分布在所有谱系中,其中,L4谱系(n=19)只有人源菌株(图3)。L8谱系中含有的NDM阳性大肠杆菌数量最多,共Ill株(24.0%),其中33(29.7%)、35(31.5%)、23(20.

19、7%)和20株(18.0%)分别来自人、猪、鸡和苍蝇(图3)。LIl谱系中的87株菌来自上述4个来源,其中,63株(72.4%)来自人(图3)。同一来源的菌株通常形成相似性极高的分支(大于99.95%,图4中以粉红色三角形表示),而在C1C5分支中存在来自两个或两个以上不同来源的亲缘关系紧密的菌株(图4)。图3.463株NDM阳性大肠杆菌的系统发育树。谱系由不同颜色的树枝和树右边的竖条表示。右侧桑基图展示了各谱系、分离来源、携带blaNDM的质粒类型和NDM蛋白亚型之间的分布和联系。所有的竖条都标记有相应的标签。图4.系统发育树显示了5个主要分支C1-C5中不同宿主间最接近的遗传关系。基于46

20、3株NDM阳性菌株的SNP差异利用最大似然法构建了中值有根系统发育树。平均分支长度小于00005(99.95%)的簇被折叠并用粉色三角形表示。三角形的大小与演化支中菌株的数量成正比。不同宿主之间亲缘关系近的菌株用红色三角形表示,同时5个遗传关系最近的分支标记为C1C5.属于这5个分支的菌株以粉红色圆圈背景显示,并标记了确切的分支长度。这些菌株的名称和宿主分别用字符和不同颜色的圆点表示。除105个不确定质粒类型的RaNDM短序列和一个位于染色体上的Rgndmb序列外,其余357个菌株中HaNDM的质粒属于19种不同的InC类型(图3);其中,b/。NDMs(332,71.7%)为最常见NDM亚型

21、,主要位于InCX3型质粒上(255,76.8%)。值得注意的是,lncB/0/K/Z型质粒(=40,8.6%)只携带bNDM9。总体而言,nc3(275,59.4%)是所有来源和谱系的主要质粒类型,在117(42.5%)、74(26.9%)、50(18.2%)和34株(12.4%)分别来自人、猪、鸡和苍蝇的菌株中发现。NDM阳性的InCB/O/K/Z质粒主要在鸡源菌株中发现(35,87.5%),而人源菌株中不存在该类型质粒(图3)。(六)NDM阳性大肠杆菌的来源预测采用DAPC方法预测463株NDM阳性分离株的来源,其中,测试组共463株NDM阴性株,训练组共1389株NDM阴性株(比例为1

22、:1:3)。虽然DAPC模型预测了来自人、鸡和猪的菌株之间存在密切的遗传关系,但仍可观察到一个可区分的边界图5(a),随后选择了代表结果预测效果最好的前600个主成分(共1200个)进行分析图5(b)0在测试组中该模型预测的菌株来自鸡、苍蝇、人和猪的准确率分别为87.4%、83.9%、88.3%和83.5%图5(c).通过该模型,19.5%S=24)、8.1%(/?=10)和1.6%S=2)的鸡源NDM阳性大肠杆菌(n=123)可能来自人、猪和苍蝇,而27.3%(27/99)的猪源NDM阳性大肠杆菌可能来自人图5(d。虽然30.8%(60/195)的人源分离株似乎没有改变宿主,但预计分别有53

23、.8%(n=105)和14.9%(n=29)的菌株分别来自鸡和猪图5(d)。所有苍蝇源菌株(n=46)分别来自人(n=5,10.9%)、鸡(”22,47.8%)和猪(”19,41.3%)图5(d),图5.使用DAPC模型对NDM阳性大肠杆菌溯源。(a)训练集的散点图。每个点代表一个菌株,不同的颜色表示不同的来源。椭圆表示每个原点的95%Q.(b)DAPC模型的交叉验证。(C)测试集的表格图。(d)预测集的表格图。行和列分别表示报告的来源和预测的来源,每个正方形的大小表示相应的菌株数量。PCA:主成分分析。四、讨论研究数据涵盖来自22个省(直辖市)的样本,代表了中国最大的临床CREC监测研究。来

24、自中国30家医院的感染相关大肠杆菌分离株中CREC的总流行率(2.1%f63129799zCHINEL2016)高于欧洲(0.1%,未提供/121582,EARS-Net,2016)和美国(0.3%,5/1916,30家医院,20152017年)。本研究中CREC感染与男性存在正相关,这与之前的报道一致。在我们目前的研究中发现,婴儿(1岁)被CREC感染的概率最高,近期一篇关于美国儿童中有较高CRE感染率的报道也印证了这一点。婴儿的微生物群落变化较大,容易受到抗生素的干扰,这可能为CRE定植提供了机会,这已被证明是老年重症监护室(ICU)患者内源性感染的风险因素。有趣的是,人类和动物中使用的6

25、.内酰胺类药物的排放与CREC感染没有显著相关性,这意味着6内酰胺类药物向环境中的排放可能不是CREC传播的主要驱动因素。研究表明,6内酰胺类药物在室温会迅速降解,这可能导致它们在环境中几乎不产生选择压力。研究首次确定了畜禽养殖量与人类CREC感染之间的正相关关系。微观基因组分析提供了几个支持这种联系的证据。首先,我们观察到在中国三个省的畜禽养殖场粪便样本中,Rondm的丰度很高,而不含有bpc.其次,从人、鸡、猪和苍蝇分离株中鉴定出具有相似的携带Rondm的InCX3质粒,然而在目前和以前的研究中bpc很少在动物中被发现。第三,鉴定了来自人和动物源NDM阳性大肠杆菌核心基因组序列之间的密切关

26、系。第四,溯源分析显示,人和动物源的NDM阳性大肠杆菌之间的分界模糊。综上所述,这些观察表明NDM阳性CREC菌株在动物和人类之间存在显著的传播风险。值得注意的是,动物源性食品的消费与CREC感染间的关系仅在单变量分析时呈现显著相关,在多变量逻辑回归分析中则无显著相关,表明其他可能的传播途径(如污染物传播)也是RaNDM质粒或CREC扩散的一个因素。例如,在包括医院污水、废水、饮用水、自然水道和其他来源的多个环境样本中均检测到Ronom以及伴侣动物、苍蝇和野生动物中也分离到携带NDM的菌株。此外,先前研究报道了CREC通过环境在动物和人类之间传播。以上证据均表明环境污染物是潜在的传播驱动因素。

27、RaN.在动物和人类之间的成功传递可能包括以下证据。首先,RQNDM-InCX3质粒的流行率为76.8%,出现在各主要谱系的CREC内(图3),提示该质粒可能是非宿主特异性质粒;此外,本研究中可能低估了b0NDMncX3质粒的流行率,因为大量携带NDM(n=105)的序列太短,无法确认其质粒类型,许多可能是InCX3型质粒。其次,携带HaNDM的质粒一一特别是RaNDM-InCX3一对宿主细菌几乎没有适应性代价,增加了其在不同生态系统中的稳定性和持续性。第三,除了兽医禁用的碳青霉烯类药物外,包括青霉素和头抱菌素在内的&内酰胺类药物,是中国(2018年31921)和其他116个世界动物健康组织成

28、员国(2015-2017年间占总抗菌使用量的16.1%)用于食品动物养殖的第二大类抗菌药物。这些常用的B-内酰胺药物主要被用于饮用水中,可能对畜禽中NDM阳性细菌的选择和持续存在发挥重要作用。第四,RaNDM和包括mcr;、/oR和tet等其他耐药基因在可转移质粒上共存,可分别通过使用黏菌素、氟苯尼考和四环素将blaNDM共同筛选出来。尽管我们的研究证据表明了CREC在动物和人类之间的传播,但仍存在一些局限性。第一,采样方案可能在时间跨度方面不完全一致;然而,该研究的数据集足够大,并通过流行病学和基因组学分析,提供了CREC在动物和人类之间传播的有力证据。第二,尽管以前的大规模监测表明NDM可

29、能是大肠杆菌中最普遍的碳青霉.烯酶类型,但研究中仍缺少临床CREC菌株中碳青霉烯酶类型的数据。第三,我们直接分析了动物和人类感染之间的关联,但没有分析患者的肠道定植情况和相应的动物源性食品。然而,先前的研究表明,血液和尿路感染的病原体经常来自肠道共生菌或食源性致病菌。第四,包括大肠杆菌、肺炎克雷伯菌和梭状芽也杆菌在内的肠杆菌类细菌是临床感染菌中RaNDM的主要宿主。本研究只关注了NDM阳性的大肠杆菌,没有关注其他宿主菌,因为它是动物中主要的碳青霉烯耐药细菌。第五,部分基因组数据来自其他研究,存在信息缺失的情况(例如,分离时间和地点、真实来源和报告来源之间不一致等),可能给DAPC模型引入了系统

30、误差。然而,在系统发育树的几个簇中发现了菌株间的强关联,并且在不同的宿主中发现了相同的质粒,表明了CREC在动物和人类之间的传播。从“全健康的理念出发,为对抗CREC在动物和人类之间的传播,应该考虑以下的行动:由于撤销抗菌药物已被证明可以减少细菌耐药性的流行率,因此可以通过更合理地使用B-内酰胺类药物和其他可能产生共筛选效应的抗菌药物,进而减少畜禽养殖场携带RaNDM细菌的流行。中国畜禽养殖业正向大规模集约化养殖转变,并逐渐减少农村后院养殖系统,这将有助于加强对卫生、抗菌药物消费和感染控制措施的监管。此外,还应减少仅为促生长或预防疾病使用的药物量。食物供应链、环境和医院在耐药病原体的传播中发挥

31、着重要作用;因此,建立完善的动物(农场和周围环境)和人类(定植、感染和生活环境)耐药菌/耐药基因的监测体系是限制CREC传播的必要前提。在动物和人源菌株中,有很大比例的RaNDM基因(70%)与InCX3质粒相关(图3);因此,需进一步研究宿主的偏好性、持久性以及消除b*DMncX3质粒的干预措施。农场、屠宰场、食品加工厂和零售店应通过有效的清洁和消毒程序切断CREC传播途径,以减少CREC对非生物表面的污染。此外,由于苍蝇和鸟类在CREC的传播中发挥着重要作用,需要充分了解它们对CREC传播的贡献。最后,中国农业农村部于2020年7月停止使用所有抗生素作为促生长剂之后,其他替抗方法如CRlSPR-CaS9、噬菌体和传统中药是有希望对抗这种无处不在的病原体的新武器。五、结论研究揭示了全国范围内临床CREC感染和地区畜禽养殖产量间的显著相关性,同时揭示了鸡、猪养殖场粪便样本中blaNDM丰度的普遍性。此外,在猪养殖场发现了长期存在的b/ONDM阳性大肠杆菌,并阐明了不同宿主来源菌株基因组间的界限模糊。综上所述,这些证据表明了CREC在动物和人类之间的关联和传播。综合所有证据,我们假设CREC首先出现在医学临床中,然后被引入至畜禽养殖业中,而猪和鸡正是CREC持续存在的有利宿主,这导致了CREC通过食物链或通过环境载体生物和(或)非生物在人类和动物之间的循环传播。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 生活休闲 > 在线阅读


备案号:宁ICP备20000045号-2

经营许可证:宁B2-20210002

宁公网安备 64010402000987号