引言
2021年12月4日，大連海洋大學(xué)海洋科學(xué)與環(huán)境學(xué)院在JCR Q1區期刊《Water Research》上發(fā)表了題為“Metagenomics combined with comprehensive validation as a public health risk assessment tool for urban and agricultural run-off”的研究論文，該論文提出了一個(gè)由病原體篩選、病原體基因分型和風(fēng)險評估三個(gè)階段組成的公共健康風(fēng)險評估流程：1）運用宏基因組測序、微流控芯片和qPCR對病原體進(jìn)行篩選，2）多位點(diǎn)序列分型(MLST)或常規PCR進(jìn)行病原菌分離和基因分型；3）對宏基因組數據的毒力基因進(jìn)行致病性評估。
譯名：宏基因組學(xué)結合綜合驗證作為城市和農業(yè)徑流公共健康風(fēng)險評估工具
期刊：Water Research
影響因子：11.236
發(fā)表時(shí)間：2021年12月4日
合作單位：大連海洋大學(xué)海洋科技與環(huán)境學(xué)院
研究方法：宏基因組測序；微流控芯片；qPCR；MLST分型
研究背景
水生病原體可引起嚴重疾病，如霍亂、鉤端螺旋體病、傷寒或腹瀉。鑒于這些新出現的疾病的不可預測性和高度多樣性，水傳播病原體監測的關(guān)鍵挑戰之一是優(yōu)先為局部暴發(fā)做好準備。為應對這些挑戰，該研究分三個(gè)步驟進(jìn)行公共衛生風(fēng)險評估：潛在病原體篩選、病原基因型確認、與當地臨床資料比較。
摘要
研究者選擇了兩個(gè)具有潛在公共衛生問(wèn)題的地區(東港和湛江)進(jìn)行了宏基因組測序，對毒力基因進(jìn)行了致病性評估。結果表明，共檢測到90多種細菌病原體，其中56種含有毒力基因。在東港，發(fā)現了一株致病性不動(dòng)桿菌屬(Acinetobacter sp.)，在2018年和2020年達到了較高的豐度，所有致病性弧菌都在2019年10月達到頂峰。在湛江，鮑曼不動(dòng)桿菌（A. baumanni）等腸桿菌科細菌種類(lèi)在2019年和2020年大量存在，而氣單胞菌屬（Aeromonas）和弧菌屬（Vibrio spp）在2017年11月濃度達到頂峰。通過(guò)MLST對40個(gè)菌種進(jìn)行了分離與分型，其中一半是臨床上的流行基因型。此外，在2017年采集的水樣中發(fā)現了非洲豬瘟病毒(ASFV)，這在中國首次報告非洲豬瘟（2018年）之前，還檢測到甲型肝炎病毒(HAV)和腸道病毒71型(EV71)等RNA病毒，濃度分別在2020年4月和2018年4月達到峰值。HAV和EV71的動(dòng)態(tài)與當地流行趨勢一致。最后，基于毒力基因圖譜，該研究確定了兩個(gè)城市污水中的風(fēng)險水平。這一工作流程說(shuō)明了對當地流行病早期預警的可能性，有助于提前為當地的特定病原體做好準備。
研究方法
采樣點(diǎn)：東港市大洋河(DY)和東溝(DG)，湛江城區(DA)和農業(yè)區(DB)
采樣時(shí)間：2017-2020年間，每年采集兩次水樣（半年一次）
檢測方法：MLST基因分型，微流控芯片篩選，病毒基因分型，宏基因組測序分析
實(shí)驗流程：
圖1.廢水中新出現病原體風(fēng)險評估的工作流程示意圖。
工作流程包括潛在病原體的篩選、病原體基因分型和風(fēng)險評估。
主要研究結果
1.東港水樣中潛在致病菌的篩選
在宏基因組數據中，89.3%的reads被分配到細菌中。參考毒力因子數據庫(VFDB)，使用MetaPhlAn2對每個(gè)樣品進(jìn)行分類(lèi)鑒定以進(jìn)一步研究潛在病原體。
在DY的宏基因組數據中大量檢測到40種病原體，包括水傳播致病菌和呼吸道致病菌。2019年后，鮑曼不動(dòng)桿菌、副溶血性弧菌、霍亂弧菌、溶藻弧菌、創(chuàng )傷弧菌、肺炎克雷伯菌、腸沙門(mén)氏菌、結核分枝桿菌和空腸彎曲桿菌等幾種檢測到的病原體的相對豐度有所增加。相比之下，其他物種在2017年至2020年期間有所減少，例如韋羅尼氣單胞菌、嗜肺氣單胞菌和邦格里沙門(mén)氏菌(圖2B)，隨后培養實(shí)驗，檢測出其中32種細菌。
在DY數據中，城市污水中的細菌組成在2017年至2018年期間保持穩定(圖4)。此后，變形桿菌的相對豐度在2019年4月至2019年下降到17.3%，但在2020年9月后成為優(yōu)勢類(lèi)群。宏基因組數據中共檢測到76種細菌病原體，其中32種用于可視化(圖2D)。致病性不動(dòng)桿菌屬在2018年和2020年達到了較高的相對豐度，所有致病性弧菌都在2019年10月達到頂峰，其他病原體的相對豐度沒(méi)有明顯的波動(dòng)。在32個(gè)選定的細菌病原菌中，在不同采樣時(shí)間的培養檢測中獲得17種(圖2D)。在研究期間，結核分枝桿菌、梅毒螺旋體和腸炎沙門(mén)氏菌的相對豐度保持不變，它們分別是結核病、梅毒和傷寒/副傷寒的臨床病原體。
圖2. 東港市農業(yè)(DY)和城市污水(DG)采樣點(diǎn)細菌群落2017-2020年動(dòng)態(tài)
A：DY細菌群落相對豐度；B：DG中細菌群落相對豐度；C：DY樣品中細菌病原體的動(dòng)態(tài)；D：DG樣品中細菌病原體的動(dòng)態(tài)
2.東港地區細菌病原菌基因分型
對單個(gè)菌種進(jìn)行MLST基因分型，以確定是否屬于致病基因型。共分離出16種細菌進(jìn)行MLST，東港地區霍亂弧菌和副溶血性弧菌分別鑒定出3個(gè)和5個(gè)序列類(lèi)型（STs），其中僅副溶血性弧菌ST3在當地臨床中被確認。7個(gè)肺炎克雷伯菌STs（包含ST11和ST23）被發(fā)現，并且在當地醫院相關(guān)感染中占主導地位。檢出5株嗜肺鏈球菌，但未檢出臨床病例，臨床和環(huán)境標本中均檢出北京分支(ST4)結核分枝桿菌，環(huán)境中的腸鏈球菌ST34和ST19在當地患者中被檢出。對于其他病原體，除陰溝腸桿菌外，每種都檢測到6種STs，且至少有一種STs出現在臨床中。雖然一些致病基因型在當地的臨床樣本中沒(méi)有被鑒定出來(lái)，但它們是普遍存在的，如金黃色葡萄球菌ST7，因此仍需要進(jìn)一步的研究。
3.東港水樣中致病病毒的篩選
選擇了60種與公共衛生風(fēng)險相關(guān)的致病病毒，其中，通過(guò)將宏基因組reads比對到參考基因組確認了20種DNA病毒，用微流控芯片對40種RNA病毒進(jìn)行了篩選，最后通過(guò)qPCR檢測確認。
在DY發(fā)現了非洲豬瘟病毒(ASFV)和豬痘病毒(Swinpox Virus)等幾種病毒。ASFV在2017年10月被檢測到，并在2018年達到頂峰(圖6A)，并通過(guò)qPCR證實(shí)。qPCR在2018年4月后連續檢測到ASFV，其中2018年4月的病毒載量為3.0×103拷貝/L。在DG地點(diǎn)未檢測到ASFV和Swinpox病毒，在兩個(gè)地點(diǎn)均未檢測到其它人畜共患病DNA病毒。通過(guò)高通量微流控芯片檢測檢測了40種致病性RNA病毒的存在。DY樣品中，經(jīng)qPCR確認，2019年僅檢測到諾如病毒、流行性腮腺炎病毒、鼻病毒和Human herpesvirus 6病毒。DG樣品中，2017年至2019年在樣本中檢測到human enterovirus 71 (EV71)，2020年4月檢測到Coxsackievirus 16 (CA16)。
4.東港水樣中病毒基因分型研究
對陽(yáng)性樣本進(jìn)行病毒基因分型，使用濃縮的DNA樣本對DY 2017和2018年水樣中ASFV B646L基因進(jìn)行PCR擴增。從DY和其他國家的32個(gè)ASFV B646L序列數據中檢測到的B646L的系統發(fā)育分析表明，前者與從遼寧和其他東歐國家的死豬中獲得的第一個(gè)ASFV分離物(CN201801/2018) 在基因上是相同的。對東港市手足口病流行時(shí)2019年4月達到峰值(5.8×102拷貝/L)的EV71濃度進(jìn)行了定量檢測(圖3B)，其中EV71的VP1基因PCR為陽(yáng)性，系統發(fā)育分析證實(shí)東港EV71屬于在遼寧省流行的C4a型，但PCR未擴增出CA16基因片段。在DG樣品中，通過(guò)qPCR檢測，從2018年10月連續檢出了甲型肝炎病毒(HAV)，峰值在2020年4月(4.7×103拷貝/L)(圖5C)，這與甲型肝炎病例大幅增加的時(shí)間吻合。進(jìn)一步通過(guò)PCR擴增得到甲型肝炎病毒VP3/P2A區，對100個(gè)公開(kāi)提供的甲型肝炎病毒基因組系統發(fā)育分析，發(fā)現DG中甲型肝炎病毒與大連分離株DL3(圖5C)有遺傳關(guān)系。其余陽(yáng)性檢測的病毒可能是由于病毒豐度低所致未被常規PCR擴增出。
圖3.代表病毒的動(dòng)態(tài)變化和系統發(fā)育關(guān)系。
5.湛江水樣中潛在致病菌的篩選
宏基因組數據中大約90%的reads與細菌相對應。DA樣品中，變形菌門(mén)數量最多，其次是疣微菌門(mén)和放線(xiàn)菌(分別為35%和58%)(圖4A)。DB樣品中，變形菌門(mén)從2017年的50.1%增加到2020年的84.7%(圖4A)。兩個(gè)地點(diǎn)的微生物組成自2018年11月后出現了差異。在每個(gè)采樣時(shí)間，豐度 top20屬代表了微生物群落的35%到57%(圖5B)。在2019年(5.8%)和2020年(4.8%)采集的樣品中發(fā)現了弧菌屬(Arcobacter)，而氣單胞菌屬(Aeromonas)和弧菌屬(Vibrio Spp)在2017年11月達到峰值(分別為0.43%和4.1%)，在DA中還檢測到其它幾種細菌病原體(圖5C)，所有時(shí)間均檢出假結核耶爾森菌（Yersinia pseudotuberculosis），2018年5月達到峰值(0.00013%)，其他菌種也達到峰值，如不動(dòng)桿菌屬、肺炎克雷伯菌屬等?；谂囵B的分析證實(shí)，DA和DB位點(diǎn)分別存在22種和18種致病物種(圖5C，D)。對兩個(gè)城市的宏基因組數據進(jìn)行直接分類(lèi)。每個(gè)數據集中的宏基因組組裝基因組(MAG)數量從7到23個(gè)不等，相對豐度在0.215%到1.325%之間。然而，對MAGs的注釋結果表明，它們都不含有毒力因子，說(shuō)明該方法只能檢測高豐度的病原體。
圖4.湛江地區細菌群落動(dòng)態(tài)研究
6.湛江水樣細菌病原菌基因分型
在DA和DB中分別鑒定出4種和13種副溶血性弧菌的序列類(lèi)型 (sequence types，STs)。在10份湛江樣品中有3份分離到ST1740和ST1799。湛江城區共分離到4株霍亂弧菌STs(ST69、ST1170、ST1171和1株新的STs)，其中ST69屬于O1大流行血清群。在城市和農業(yè)樣品中均檢出肺炎克雷伯菌ST23、ST35和ST86。城區樣品中檢出三株流行基因型ST11、ST34和ST13沙門(mén)氏菌，并出現在廣東省臨床樣品中。
7.湛江水樣中病毒的檢測和基因分型
湛江宏基因組數據中3.2%的reads分類(lèi)為DNA病毒序列，其中97.8%對應細菌噬菌體。2017年和2019年還檢測到四種人類(lèi)病原體，其中包括與兒童腹瀉有關(guān)的Tanapox病毒。2018年5月后，湛江也檢測到低豐度的猴痘病毒和Human herpesvirus 6，2019年5月，檢測到很低(0.001%)的Swinpox病毒，可能由于其低豐度然而，qPCR檢測未檢測到上述病毒。對于RNA病毒，使用微流控芯片和qPCR技術(shù)在2018年11月、2019年5月和11月湛江樣品中檢測出諾如病毒，載量分別為3.4×104、1.78×103和2.9×103拷貝/L，2018年11月檢測到豬札幌病毒(Sapovirus)，為3.2×103拷貝/L，但常規PCR沒(méi)有擴增出這些病毒的DNA片段。
8.機會(huì )性病原體豐度定量和WGS
兩個(gè)城市的樣品宏基因組數據中給定病原體在多個(gè)點(diǎn)的相對豐度顯著(zhù)增加。例如，2018年4月，DG樣本中A. baumannii, Acinetobacter lwoffi和 Acinetobacter junii增加了100倍以上，2018年10月，DY樣本中發(fā)現培養基Aeromonas media增加了10倍以上(圖2)。東港水樣的qPCR檢測表明，污水中病原菌豐度在7.1×104~1.33×109拷貝/L之間。2018年4月和2020年4月，A. junii, A. baumannii和A. lowffii三種細菌在DG中的豐度超過(guò)了半數致死量（LD50）水平。
從湛江樣品中選擇了4株無(wú)MLST方案的病原菌，A. junii , V.alginolyticus, S. algae 和A. hemolyticus進(jìn)行全基因組測序（WGS）。結果表明，V. alginolyticus與本地水產(chǎn)養殖場(chǎng)的菌株有親緣關(guān)系，表明它們是該地區特有的(圖5A)。S. algae ZJ61與我國山東的4株臨床分離株有親緣關(guān)系，而S. algae ZJ-SH在一個(gè)獨立的分支上(圖6B)。結合75個(gè)A. junii基因組進(jìn)行系統發(fā)育分析，2018年?yáng)|港水樣分離株JZ20與2020分離株JZ102(圖6C)聚為一類(lèi)。對57個(gè)A. hemolyticus基因組分析，湛江分離株ZJ54與從醫院分離的SZ3452有親緣關(guān)系(圖6D)，分離株HK30與醫院的菌株XH900聚在一起。
圖5.根據Vibrio alginolyticus (A)Shewanella algae (B) Acinetobacter junii (C) Acinetobacter haemolyticus (D)核心基因組序列，
用最大似然法構建系統發(fā)育樹(shù)本研究中獲得的分離株用星號表示
9.污水樣本的公共衛生風(fēng)險評估
細菌：低風(fēng)險，在宏基因組數據中可以檢測到病原體，但不能培養；中風(fēng)險，宏基因組數據中檢索到病原體及其關(guān)鍵毒力基因，但沒(méi)有分離到致病基因型；高風(fēng)險，在隨后的分離和MLST能確認致病基因型。病毒：中風(fēng)險，qPCR/微流控芯片/宏基因組數據檢測到；高風(fēng)險，用傳統的PCR可進(jìn)一步擴增。
對兩個(gè)城市采樣地點(diǎn)的公共衛生風(fēng)險進(jìn)行了評估，通過(guò)與VFDB數據庫比對，確定了95種參考病原體。主要毒力因子包括(1)毒素，(2)鞭毛和菌毛，(3)脂多糖和抗原，(4)致病島。
24個(gè)物種含有關(guān)鍵的毒力因子，這些因子在不同的水樣中情況不同，因此城市和農業(yè)水樣有不同的健康風(fēng)險特征。2020年4月收集的東港城市污水總體健康風(fēng)險最高(圖6A)，包含1種病毒病原體高風(fēng)險和6種細菌病原體高風(fēng)險。2019年至2020年，報告的甲型肝炎病例有所上升，同時(shí)甲型肝炎病毒的在樣品中濃度很高。從2017年-2020年病例數增加數量（183-622例）和病毒檢測量（增加10倍）不一致，可能是因為城市臨床監測還沒(méi)建立起來(lái)，無(wú)癥狀或未被識別的患者未被報告，這突顯了環(huán)境監測的重要性。此外，農業(yè)污水中ASFV、蠟狀芽孢桿菌和腸球菌有高風(fēng)險，需要干預以保證食品安全。
湛江水樣結果表明A. baumanni, B. cereus, C. jejuni, E. coli, K. pneumoniae, S. enterica, S. sonii, S. pneumoniae和V. cholerae在城市地區具有很高的健康風(fēng)險(圖6B)。其中，2018年5月采集的廢水中含有7種高風(fēng)險病原體，Aeromonas sp., B. cereus和V. parahaemolyticus是農業(yè)污水的主要關(guān)注對象，湛江的城市和農業(yè)廢水中沒(méi)有發(fā)現高風(fēng)險病毒病原體。
圖6.東港(A)和湛江(B)污水樣品中細菌和病毒病原體的公共衛生風(fēng)險評價(jià)
藍色：低健康風(fēng)險；灰色：中等健康風(fēng)險；紅色：高健康風(fēng)險
總結
該研究提出了一種同時(shí)檢測多種毒力基因的污水公共健康風(fēng)險評估框架，該框架具有很大的病原體環(huán)境監測應用潛力。廢水監測是可視化病原體動(dòng)態(tài)和提供病原體與疾病之間關(guān)系的重要洞察力的有力平臺。在農村，臨床監測或報告系統不夠充分，這一風(fēng)險評估工具還可以監測農村地區病原體的流行病學(xué)。
該篇論文運用了百邁客云分析平臺進(jìn)行宏基因數據分析，如您有任何疑問(wèn)，可以討論區留言或者發(fā)郵件給我們（郵箱地址：tech&@biomarker.com.cn)
正值年末，百邁客微生物系列產(chǎn)品促銷(xiāo)活動(dòng)火熱進(jìn)行中，無(wú)論是測序還是分析，價(jià)格優(yōu)惠，歡迎掃碼聯(lián)系我們~
百邁客微生物組一站式服務(wù)
百邁客擁有完備的實(shí)驗分析平臺，可以開(kāi)展宏基因組（Illumina/ONT），二代微生物多樣性，全長(cháng)微生物多樣性，細菌/真菌基因組測序及數據分析服務(wù),百邁客在微生物組領(lǐng)域深耕多年，致力于提供高質(zhì)量的科研服務(wù)。目前已搭建微生物多樣性分析平臺、宏基因組分析平臺、微生物基因組和宏基因組binning分析四大平臺，助力微生物組研究。
京公網(wǎng)安備 11011302003368號