兒童復發性腮腺炎患者與正常兒童腮腺液微生物群落多樣性的差異

李 瑩 祁森榮 華 琳 尚佳健*

(1.首都醫科大學附屬北京口腔醫院兒童口腔科， 北京 100050；2.首都醫科大學附屬北京口腔醫院放射科， 北京 100050；3.首都醫科大學生物醫學工程學院數學與生物信息教研室，北京100069)

兒童復發性腮腺炎(juvenile recurrent parotitis，JRP，以下簡稱兒復腮)是發生于兒童腮腺的感染性疾病，發病率僅低于流行性腮腺炎。目前臨床的治療方法，特別是在急性期，主要采取應用廣譜抗生素治療，尚無針對性的抗炎治療方法[1]。兒復腮病因復雜，發病機制目前尚不十分清晰，可能與以下因素有關：細菌或病毒逆行性感染、導管系統發育異常、遺傳因素、自身免疫性疾病、變態反應等，其中細菌感染是最主要的原因之一[2]。目前對兒復腮的細菌學研究主要使用培養法，而其中一部分微生物，很難或者不能通過傳統的純培養分離方法得到。近年來，隨著分子生物學技術的發展，從基因學角度探索口腔微環境的組成和結構，使越來越多的未培養微生物得以檢測。16SrRNA克隆測序法，為各種微環境中許多未知菌種的鑒定與分析提供了非培養依賴的方法。近年來該技術已作為一項新技術逐漸引入到口腔微生物系統分析中，如牙周病、感染根管的細菌分析。但國內尚未見該方法用于兒復腮的細菌學檢測中的相關報道。本研究對兒復腮患兒的腮腺液進行微生物學檢測，并分析與正常兒童的差異，有助于了解其病因及發病機制，為尋找針對性治療方法提供參考。

1 對象與方法

1.1 研究對象

隨機選取2015年11月20日至2017年3月10日在首都醫科大學附屬北京口腔醫院頜面外科涎腺中心就診的兒復腮患者35例為疾病組(JRP組)，以及在兒童口腔科就診的無腮腺腫脹史的牙體疾病患者12例為對照組(normal children，NC組)。47名受試者中，男性30例，女性17例，年齡3～12歲，平均6.1歲。其中一部分患者取雙側腮腺液樣本(2個樣本)，一部分患者取單側腮腺液樣本(由于患兒配合程度和發病程度有差異)，共獲得60個臨床患者樣本，2例為空白對照樣本，操作方法相同但不收集唾液。

病例組納入標準：1)3～12歲兒童，腮腺反復腫脹(同一部位腮腺腫脹超過一次)或導管口有膿性分泌物(單側或雙側)者。2)腮腺造影表現為導管末梢點球狀擴張或腮腺超聲表現為腺體內點狀或竇狀低回聲區。3)血常規檢查結果顯示白細胞數增高，淋巴細胞正常者。

正常對照組納入標準：患兒年齡3～12歲。腮腺無腫脹史。排除標準：1)治療依從性差或伴有嚴重的全身性疾病者。2)確診為流行性腮腺炎患兒。3)腮腺只腫脹過一次，無法確診的。

符合要求的病例經患兒及家長知情同意，并由家長簽署《知情同意書》后納入本研究。

1.2 臨床檢查

收集每一位患兒的以下臨床資料：姓名、性別、出生日期、聯系方式、全身疾病情況、抗生素服用情況、主訴、現病史、?？茩z查、影像學檢查結果、診斷及治療史、取材情況、采集人等。

1.3 樣本采集

各組兒童腮腺導管口周圍直徑2 cm范圍內安爾碘消毒，在導管口處放置無菌消毒小棉粒，直徑約0.5 cm，于口外用手按摩腮腺區，并從耳后向前擠壓至腮腺導管口處，使擠出的腮腺液完全浸潤小棉粒。取出小棉粒放入無菌試管內，-80 ℃冰箱保存。

1.4 提取樣本中DNA，擴增樣本16SrRNA V3-V4區

利用340F (5′-CCCTAYGGGGYGCASCAG-3′)[3]和805R (5′-GACTACHVGGGTATCTAATCC-3′)[4]引物對對樣本微生物16SrRNA V3-V4區進行PCR擴增。

1.5 文庫構建、上機測序和Alpha多樣性指標

使用Illumina Hiseq 2500 PE250高通量測序平臺(北京諾賽基因組研究中心有限公司)對擴增成功樣本的16S-rRNA V3-V4區進行測序得到原始數據，對原始數據進行拼接過濾得到有效數據，根據97%的相似性水平對有效數據進行歸類，得到多個序列聚類操作分類單元(operational taxonomic units, OTUs)。利用QIIME[5]軟件包中的uclust[6]聚類方法將移除嵌合體后的序列歸為多個OTUs，得到OTU分組情況。然后基于SILVA參考數據庫，利用RDP Classifier[7](采用confidence cutoff=0.8)對每個樣品的OTUs進行分類學(Taxonomy)注釋，并對OTUs進行生物信息分析。

以隨機抽取序列數為40 000提取Alpha多樣性指標，包括樣品文庫的覆蓋率(goods_coverage)、豐富度指數(包括Chao1[8]與菌群豐富度指數(abundance coverage estimator,ACE)、多樣性指數(包括Shannon等[9]與Simpson[10])、observed_species等。

1.6 統計學方法

采用SPSS 19.0統計軟件進行統計分析。兒復腮患者與正常兒童腮腺液的Alpha多樣性指標的組間比較采用兩獨立樣本t檢驗，以P<0.05為差異有統計學意義。

2 結果

2.1 樣本擴增情況

對62個樣本進行DNA提取、擴增，其中50例樣本擴增成功(JRP組34例，NC組16例)，12個樣本(包含2個空白對照樣本)擴增未成功。擴增成功的50個樣本總共獲得 7 303 281條有效數據，平均讀長為459個堿基。共得到 19 399個OTUs。

2.2 稀釋性曲線(rarefaction curve)

各樣本的稀釋曲線(圖1)大多趨于平緩，表明大部分物種已經被檢測到，實驗得到的數據可用于后續分析。

2.3 Alpha多樣性指標計算

統計分析結果顯示，在豐富度指數(Chao1)、多樣性指數(Shannon與Simpson)中，除了Simpson兩組比較沒有統計學意義，其他指標組間差異均有統計學意義(P<0.001)，詳見圖2～4和表1。

圖1 各樣本稀疏曲線Fig.1 The rarefaction curve of samples

The rarefaction curves of most samples tended to be flatten, which indicates that most of species hade been detected and the obtained data could be used for further analysis.

圖2 各樣本Chao1曲線Fig.2 The Chao1 curve of samples

The difference of Chao1 between the disease group and normal group was statistically significant (P<0.001)

圖3 各樣本Shannon指數曲線Fig.3 The Shannon Index Curve of Samples

The curve of each sample tends to be smooth and reach saturation, which indicates that the sequencing quantity of this experiment can basically reflect the diversity of microflora. In this study, the difference of Shannon index between the disease group and normal group was statistically significant (P<0.01).

圖4 simpson指數曲線Fig.4 The Simpson index curve of samples

There was no significant difference between the disease group and normal group.

表1 兒童復發性腮腺炎患者與正常兒童腮腺液的Alpha多樣性指標比較Tab.1 Comparison of Alpha diversity in parotid fluid between JRP and the normal children index

2.4 微生物構成的豐度排序(屬水平)

JRP組按豐度高低排名前5位的微生物為：鏈球菌屬、嗜血桿菌屬、 奈瑟氏菌屬、孿生球菌屬、氏菌屬。NC組的排序為：鏈球菌屬、嗜血桿菌屬、孿生球菌屬、奈瑟氏菌屬、羅氏菌屬。

3 討 論

在本研究中，12個樣本擴增未成功。其中2例空白樣本均未擴增成功，可以說明獲取樣本方法可靠，沒有可檢測出來的生物污染。另外10個樣本(JRP組2例，NC組8例)未擴增成功的原因可能是腮腺液內微生物含量低或者是取樣本時所獲得的腮腺液較少。

口腔是消化道的起始端，口腔微生物群落是指定植于人體口腔的微生物集合[11]。2008年人類口腔微生物組數據庫(the human oral microbiome database, HOMD)正式啟動，該數據庫主要建立在16SrRNA檢測的基礎上，記錄了13門619種口腔微生物的分類[12]。徐欣等[13]通過收集不同年齡、不同牙列階段的臨床樣本進行細菌學檢測，結果顯示人群中口腔微生物的分布模式差異存在統計學意義；隨著年齡改變，口腔細菌在“門”水平呈現動態改變；嬰兒口腔內的微生物多樣性較成人低。腮腺結構和功能存在個體差異，腮腺導管開口于口腔內，口內微生物可進入腮腺導管系統。不同個體腮腺微生物是否具有多樣性差異，相關研究較少。本實驗通過研究口腔生物群落的多樣性，希望深入理解口腔微生物與疾病之間的相互作用關系。

人體口腔微生物群落的穩定是指不同種群間達到動態的生態平衡的結果[14]。當口腔微生物群落與宿主處于平衡狀態時，發揮屏障作用，阻止外源性致病菌的入侵；當平衡狀態被破壞時，可能會發生口腔疾病[15]?？谇晃⑸镞€可作為病灶，與全身系統性疾病關系密切。Becker等[16]對齲病微生物多樣性進行研究，發現多種細菌的組成和結構在齲病患者與無齲人群之間差異有統計學意義, 國內研究[17]也有相同結論。臨床研究[1]顯示兒復腮患兒的腮腺腫脹發作與上呼吸道感染等疾病相關，目前普遍認為微生物是從口腔逆行進入腮腺導管系統，導致腮腺炎性反應的發生。本研究顯示正常兒童腮腺內有復雜的微生物環境。兒復腮患兒腮腺液中豐度較高的細菌都可以在正常兒童腮腺內發現，說明兒復腮可能是由條件致病菌參與的感染性炎性反應，而且致病菌可能并非是在患病時由口腔內逆行進入腮腺內的。

本研究把兒復腮患者與正常兒童腮腺液兩組樣本的微生物群落分別作為一整體，研究兩組之間的Alpha多樣性指標的差異。Alpha多樣性指標包括：Goods-coverage、豐富度指數(Chao1與ACE)、多樣性指數(Shannon與Simpson)、Observed_species等。其中Goods_coverage表示各樣品文庫的覆蓋率，本實驗中各樣本的覆蓋率均達到97%以上；豐富度指數衡量樣本中物種種類個數；多樣性指數用來衡量微生物群落的異質性，其中ACE是用來估計群落中含有 OTU 數目的指數；Observed_species表示樣品中OTU的個數。本研究結果表明兒復腮患者組中微生物群落的豐富度、多樣性低于正常兒童組。提示兒復腮不僅僅是單一菌種感染致病，而是腮腺內原有條件致病菌占比增多造成生態關系失衡的結果[18]。

采用廣譜抗生素全身治療，可能會進一步干擾菌群關系的平衡。李琰[19]的研究顯示中藥長皮膏在治療慢性下肢靜脈性潰瘍中療效顯著，分析機制可能是治療增加了創面菌種的多樣性，通過競爭性生長，弱化致病菌的致病作用，促進潰瘍面的愈合。兒復腮臨床治療中碘化油導管灌注的臨床效果較好，原因可能是碘化油為廣譜抑菌劑，有利于抑制致病菌，恢復正常菌群環境[20]。因此兒復腮的治療策略需要將患兒與微生物群落作為一整體加以考慮，以恢復或重建適宜的腮腺微生態環境為目標。

微生物群落的失衡是導致口腔多種感染性疾病的重要因素，如齲病[21]、牙周疾病[22-23]等，本研究提示兒復腮也具有相同機制。治療策略及預防復發以恢復或重建的腮腺微生態環境為目標。兒復腮患者腮腺內菌群失衡的機制及主要致病菌的構成、豐度等仍需進一步研究。

[1] 謝李松.兒童復發性腮腺炎病因學研究及治療進展[J].中國實用口腔科雜志,2012,5(12):759-762.

[2] Brook I. The bacteriology of salivary gland infections[J]. Oral Maxillofac Surg Clin North Am, 2009, 21 (3):269-274.

[3] Gantner S, Andersson A F, Alonso-Sáez L, et al. Novel primers for 16S rRNA-based archaeal community analyses in environmental samples[J]. J Microbiol Methods, 2011 ,84(1):12-18.

[4] Herlemann D P, Labrenz M, Jürgens K, et al. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea[J].ISME J, 2011 ,5(10):1571-1579.

[5] Caporaso J G, Kuczynski J, Stombaugh J, et al. QIIME allows analysis of high-throughput community sequencing data[J]. Nat Methods, 2010, 7(5): 335-336.

[6] Edgar R C, Haas B J, Clemente J C, et al. UCHIME improves sensitivity and speed of chimera detection[J].Bioinformatics，2011,27(16):2194-2200.

[7] Wang Q, Garrity G M, Tiedje J M, et al. Naivebayesian classifier for rapid assignment of rrnasequences into the new bacterial taxonomy[J]. Applied and environmental microbiology, 2007,73(16):5261-5267.

[8] Chao A. Non-parametric estimation of the number of classes in a population[J]. Nat Methods,1984,11(4):265-270.

[9] Shannon C E, Weaver W. A mathematical theory of communication[J]. Bell Sys Technic J, 1948,27(4): 379-423.

[10] Simpson E H. The Measurement of diversity[J]. Nature, 1949, 163(4148):688.

[11] Dewhirst F E, Chen T, Izard J, et al. The human oral microbiome[J]. J Bacteriol, 2010,192(19): 5002-5017.

[12] 何金枝,徐欣,周學東.口腔微生物與全身健康研究進展[J].微生物與感染,2017,12(3):139-145.

[13] 徐欣,何金枝,周學東.口腔微生物群落在口腔與全身疾病預警中的作用[J].華西口腔醫學雜志,2015,33(6):555-560.

[14] Merritt J, Kreth J, Shi W, et al. LuxS controls bacteriocin production in Streptococcus mutans through a novel regulatory component [J]. Mol Microbiol,2005,57(4):960-969.

[15] 施文元，周學東.人體口腔微生物群的相互作用[J].華西口腔醫學雜,2010,28(1):1-4.

[16] Becker M R，Paster B J，Leys E J，et al. Molecular analysis of bacterial species associated with childhood caries[J].J Clin Microbiol，2002，40(3)：1001-1009.

[17] 李志艷,趙翀,劉靜,等.青海漢族和藏族高齲患者口腔微生物多樣性差異研究[J].牙體牙髓牙周病學雜志,2017,27(5):250-253,265.

[18] Papadopoulou-Alataki E, Chatziavramidis A, Vampertzi O, et al. Evaluation and management of juvenile recurrent parotitis in children from northern Greece[J]. Hippokratia, 2015, 19(4):356-359.

[19] 李琰.中藥長皮膏對慢性下肢靜脈性潰瘍的微生態調控作用的研究[D].上海： 上海中醫藥大學,2015:1-140.

[20] Roby B B, Mattingly J, Jensen E L, et al. Treatment of juvenile recurrent parotitis of childhood:an analysis of effectiveness[J]. JAMA Otolaryngol Head Neck Surg, 2015, 141(2):126-129.

[21] 王少果,游祥磊,車春曉,等.16SrRNA高通量測序研究有齲者唾液微生物群落結構及多樣性[J].口腔醫學研究,2015,31(5):461-465.

[22] 陳慧.口腔微生物組的結構分析及與疾病的關系研究[D].上海：上海交通大學,2016:1-189.

[23] 許志亮. 非手術根管再治療-根尖手術聯合治療難治性根尖周炎[J].首都醫科大學學報，2015，36(3)：478-482.