孫大慶，李洪飛，楊 健，宋大巍,3

（1.黑龍江八一農(nóng)墾大學(xué) 國家雜糧工程技術(shù)研究中心，黑龍江 大慶 163319；2.黑龍江八一農(nóng)墾大學(xué) 牡丹江食品與生物技術(shù)創(chuàng)新研究院，黑龍江 牡丹江 157000；3.黑龍江八一農(nóng)墾大學(xué)食品學(xué)院，黑龍江 大慶 163319）

短乳桿菌是乳桿菌屬重要的菌種之一，它們廣泛分布于植物、動物及人體多種復(fù)雜環(huán)境中，是乳桿菌中環(huán)境適應(yīng)能力很強(qiáng)的菌種之一。短乳桿菌具有發(fā)酵能力強(qiáng)、耐酸、產(chǎn)細(xì)菌素等多種優(yōu)良表型[1-3]，至今已廣泛應(yīng)用于食品、醫(yī)療和養(yǎng)殖等行業(yè)，與人們飲食、衛(wèi)生和健康息息相關(guān)[4-6]，因此開展短乳桿菌相關(guān)研究具有重要的經(jīng)濟(jì)和社會應(yīng)用價(jià)值。

目前，NCBI的Genome數(shù)據(jù)庫已收錄60 個短乳桿菌天然質(zhì)粒基因組序列，從質(zhì)粒大小、數(shù)量和復(fù)制類型而言，短乳桿菌是乳桿菌屬中攜帶天然質(zhì)粒最豐富的菌種之一[7]，這一現(xiàn)象表明，質(zhì)粒在短乳桿菌生長、代謝、遺傳和進(jìn)化等生物學(xué)過程中可能扮演著特殊且重要的角色，此外，眾多天然質(zhì)粒的發(fā)現(xiàn)和測定為短乳桿菌基因工程改造和生物技術(shù)應(yīng)用提供了必要的研究和材料基礎(chǔ)。然而，目前對短乳桿菌天然質(zhì)粒的了解和研究比較有限，現(xiàn)有研究主要集中于載體構(gòu)建[8-10]和特殊功能質(zhì)粒[11-12]。本研究以短乳桿菌天然質(zhì)粒為研究對象，在原有質(zhì)粒復(fù)制起始蛋白（replication initiation protein，Rep）進(jìn)化樹分類方法[13-14]基礎(chǔ)上，進(jìn)一步結(jié)合質(zhì)?；蚪M共線性分析方法，探索建立一種更加準(zhǔn)確、有效的天然質(zhì)粒的分類方法，從而為今后短乳桿菌天然質(zhì)粒的科學(xué)分類和理性應(yīng)用提供有益的探索和研究基礎(chǔ)。

1 材料與方法

1.1 質(zhì)?；蚪M數(shù)據(jù)來源

短乳桿菌質(zhì)粒基因組由NCBI Genome數(shù)據(jù)庫下載獲得，質(zhì)粒詳細(xì)信息見表1。質(zhì)粒pC194、pMV158、pUCL287基因組登錄號分別為NC_002013.1、NC_010096.1、X75607.1。質(zhì)粒pAD1和pIP501基因組序列來自文獻(xiàn)[15-16]。

表1 短乳桿菌質(zhì)粒信息Table 1 Information about plasmids in Lactobacillus brevis

1.2 質(zhì)?；蚪M序列的收集、整理和初步分析

篩查Genome數(shù)據(jù)庫中所有短乳桿菌質(zhì)?；蚪M序列完整性和每個質(zhì)粒Rep編碼情況，對注釋或疑似的Rep進(jìn)行BLAST比對分析，鑒定其保守結(jié)構(gòu)域，確定其所屬蛋白家族。

1.3 方法

1.3.1 質(zhì)粒Rep進(jìn)化樹構(gòu)建及分類

利用MEGA 6.0軟件[17]中Muscle法進(jìn)行Rep氨基酸多序列比對，之后利用Neighbor-Joining法進(jìn)行質(zhì)粒Rep進(jìn)化樹構(gòu)建。進(jìn)化樹檢驗(yàn)采用Bootstrap法進(jìn)行，自展值設(shè)定為1 000。

1.3.2 質(zhì)?；蚪M共線性分析及分類

利用Mauve 2.4.0軟件[18]進(jìn)行質(zhì)?；蚪M共線性分析，所有參數(shù)為默認(rèn)設(shè)置。DNA多序列對比分析利用DNAMAN 5.2.2軟件進(jìn)行。DNA重復(fù)序列分析利用DNASTAR 7.1.0 GeneQuest軟件包進(jìn)行。

2 結(jié)果與分析

2.1 質(zhì)?；蚪M基本特征

表2 短乳桿菌質(zhì)?；蚪M和Rep基本特征Table 2 Basic characteristics of plasmid genome and Rep in Lactobacillus brevis

經(jīng)檢索、比對和統(tǒng)計(jì)分析，Genome數(shù)據(jù)庫收錄的60 個質(zhì)粒，有9 個質(zhì)?；蚪M序列不完整或存在拼接錯誤，其余51 個質(zhì)粒為完整基因組序列，完整質(zhì)?；蚪M特征及編碼Rep情況見表2。由表2可知，短乳桿菌51 個完整質(zhì)?；蚪M大小為1.81～107.03 kb，GC含量為34.40%～45.63%，其中39 個質(zhì)粒編碼了Rep，35 個質(zhì)粒含有1 個Rep（質(zhì)粒pRH45II的Rep為本研究推定），4 個質(zhì)粒含有2 個Rep，12 個質(zhì)粒不含有已知Rep。根據(jù)Rep含有的保守結(jié)構(gòu)域，所有Rep可劃分為5 個Rep蛋白家族（表1）。

2.2 質(zhì)粒Rep進(jìn)化樹分類

圖1 短乳桿菌質(zhì)粒Rep氨基酸序列進(jìn)化樹Fig. 1 Phylogenetic tree based on amino acid sequences of plasmid Rep in Lactobacillus brevis

基于Rep氨基酸序列同源性，構(gòu)建了短乳桿菌質(zhì)粒和5 個蛋白家族參考質(zhì)粒Rep進(jìn)化樹（圖1）。進(jìn)化樹沒有包括假基因編碼和推定的Rep，這些質(zhì)粒將在后文論述。由圖1可以看出，短乳桿菌34 個質(zhì)粒編碼的38 個Rep和5 個參考質(zhì)粒編碼的Rep，根據(jù)各自含有的保守結(jié)構(gòu)域不同，能夠明顯聚類為5 個不同的質(zhì)粒家族，因此短乳桿菌30 個含有1 個Rep質(zhì)粒可以有效地劃分為5 個類型，即5 個質(zhì)粒家族。

家族3質(zhì)粒pl12108-8和pl12111-4的Rep（WP085769753.1和ARN96593.1）雖然都含有Rep_3家族蛋白共有的保守結(jié)構(gòu)域（pfam01051），但它們與該家族其他質(zhì)粒Rep存在很遠(yuǎn)的進(jìn)化距離和較低的自展值（40），進(jìn)一步序列比對發(fā)現(xiàn)，這2 個Rep不僅序列較短，而且與該家族其他質(zhì)粒Rep序列一致性較低，因此將家族3質(zhì)粒劃分為2 個分類亞家族。

由圖1可以看出，短乳桿菌4 個復(fù)合質(zhì)粒（2 個及以上復(fù)制子融合形成的質(zhì)粒）的8 個Rep與各自同源Rep明顯聚類。復(fù)合質(zhì)粒pl12108-2和pl12111-2含有2 個家族5 Rep，且聚類在一個分支的2 個Rep氨基酸序列一致性均為100%，這表明兩個復(fù)合質(zhì)粒由2 個不同的家族5質(zhì)粒融合形成，它們具有非常近的進(jìn)化關(guān)系，由于融合前兩個質(zhì)粒都屬于家族5，因此這兩個復(fù)合質(zhì)粒也屬于家族5。由于復(fù)合質(zhì)粒pL747-2和pl12111-4分別含有1 個家族3 Rep和1 個家族4 Rep，因此它們不能簡單歸屬于一個家族，只能暫時歸屬于兩個家族（家族3和4）。

綜上可知，根據(jù)Rep進(jìn)化樹，32 個短乳桿菌質(zhì)?？梢苑诸悶? 個不同的質(zhì)粒家族，2 個復(fù)合質(zhì)粒pL747-2和pl12111-4同時歸屬于家族3和家族4，由參考質(zhì)?？芍易?和2質(zhì)粒屬于滾環(huán)型復(fù)制質(zhì)粒，家族3～5屬于theta型復(fù)制質(zhì)粒。

2.3 質(zhì)粒基因組共線性分類

復(fù)制子是質(zhì)粒復(fù)制和生存必不可少的保守區(qū)域，因此本實(shí)驗(yàn)利用基因組共線性分析結(jié)果進(jìn)行質(zhì)粒分類時，將是否具有同源性復(fù)制子作為分類的優(yōu)先原則，其次根據(jù)質(zhì)?；蚪M共線性程度大小進(jìn)行分類。根據(jù)質(zhì)?；蚪M的同源性和差異性，對短乳桿菌所有51 個質(zhì)粒進(jìn)行了基因組共線性分析，并按照上述分類原則進(jìn)行質(zhì)粒分類，結(jié)果見圖2。

由圖2可知，短乳桿菌質(zhì)粒家族1含有3 個成員，它們都含有RCR質(zhì)粒家族pC194同源復(fù)制子，該復(fù)制子由單鏈起點(diǎn)、rep和雙鏈起點(diǎn)組成[19]，是家族1質(zhì)粒共有序列，因此可以作為該家族質(zhì)粒分類的遺傳標(biāo)記。進(jìn)一步序列分析顯示，質(zhì)粒pSD11和pLB925A02均編碼一個高度同源的轉(zhuǎn)移蛋白，這使質(zhì)粒在一定條件下可以發(fā)生接合轉(zhuǎn)移，而質(zhì)粒pL747-7沒有編碼任何功能基因，因此家族1質(zhì)?？梢赃M(jìn)一步劃分為2 個亞家族類群。

由圖2可知，短乳桿菌質(zhì)粒家族2只含有1 個成員，它與RCR質(zhì)粒家族pMV158具有同源的復(fù)制子序列，該復(fù)制子由拷貝數(shù)調(diào)控基因、單鏈起點(diǎn)、rep和雙鏈起點(diǎn)組成[20]，是它們的共有序列，同樣可以作為該家族質(zhì)粒分類的遺傳標(biāo)記。

圖2 短乳桿菌質(zhì)?；蚪M共線性Fig. 2 Plasmid genome colinearity in Lactobacillus brevis

由圖2可知，短乳桿菌質(zhì)粒家族3含有4 個成員，均未編碼Rep，含有大于4 kb共有序列，共有序列注釋了2 個共有基因，分別編碼推定蛋白和移動蛋白。這些分析結(jié)果表明，家族3質(zhì)粒均為可轉(zhuǎn)移質(zhì)粒，很可能屬于完全依賴宿主酶系進(jìn)行theta復(fù)制的質(zhì)粒，如著名的ColE1家族質(zhì)粒[21-23]。

由圖2可知，短乳桿菌質(zhì)粒家族4含有12 個成員，均含有theta復(fù)制質(zhì)粒pUCL287[24-25]同源的Rep（屬于Rep_3家族），需要說明的是，質(zhì)粒pl12113-4 Rep注釋為假基因，經(jīng)序列比對，該基因提前出現(xiàn)終止密碼子，因此編碼的Rep不完整（290 aa）。家族4質(zhì)?；蚪M整體共線性較差，除了Rep編碼序列，12 個質(zhì)粒沒有發(fā)現(xiàn)共有序列，因此Rep編碼序列可以作為家族4質(zhì)粒分類的遺傳標(biāo)記。由參考質(zhì)粒pUCL287可知，rep基因上游重復(fù)序列（即復(fù)制起點(diǎn)）具有該家族復(fù)制子典型特征[26]，因此對所有14 個Rep_3家族Rep編碼序列進(jìn)行了比對分析，分析結(jié)果見圖3。由圖3可知，家族4質(zhì)粒復(fù)制起點(diǎn)可以分為2 種類型，其中11 個質(zhì)粒與pUCL287具有相同或相近的復(fù)制起點(diǎn)（4×11 bp＋4.5×22 bp）（圖3A），而家族4質(zhì)粒pl12108-8和復(fù)合質(zhì)粒pl12111-4具有另外一種類型的復(fù)制起點(diǎn)（3.5×21 bp）（圖3B），這一分析結(jié)果與Rep進(jìn)化樹（圖1）對該家族質(zhì)粒的分類結(jié)果一致，這進(jìn)一步證明家族4質(zhì)?？梢詣澐譃? 個亞家族類型。

由圖2可知，短乳桿菌質(zhì)粒家族5含有23 個成員，包括15 個編碼Rep（屬于RepA_N家族[27-28]）質(zhì)粒和8 個非編碼Rep質(zhì)粒。根據(jù)基因組共線性和序列比對，家族5質(zhì)?？梢赃M(jìn)一步分為3 個亞家族。亞家族1含有18 個質(zhì)粒，其中3 個質(zhì)粒沒有發(fā)現(xiàn)Rep，大多數(shù)質(zhì)粒編碼了接合轉(zhuǎn)移相關(guān)基因，因此推測這些質(zhì)粒起源于一個接合質(zhì)粒。亞家族2含有3 個質(zhì)粒，均為非編碼Rep質(zhì)粒，它們共有一個約20 kb噬菌體相關(guān)基因編碼區(qū)域，因此推測它們可能起源于一個編碼噬菌體質(zhì)粒。亞家族3含有2 個同源性很高的質(zhì)粒，它們沒有編碼Rep，與家族5其他質(zhì)粒具有少量的同源性區(qū)域，但都編碼的是一些未知蛋白。

由圖2可知，短乳桿菌質(zhì)粒家族6含有8 個成員，均含有theta復(fù)制質(zhì)粒pIP501同源的Rep（屬于RepR家族[29]），多數(shù)為大型質(zhì)粒。根據(jù)基因組共線性差異，8 個質(zhì)?？梢悦黠@劃分為2 個亞家族。此外，由基因組共線性可知，雖然家族6質(zhì)粒普遍較大，編碼了眾多功能基因，但與家族4、5質(zhì)粒相比，家族6質(zhì)粒基因組普遍顯示了良好共線性，這一結(jié)果表明，家族6質(zhì)粒之間可能具有親密的進(jìn)化關(guān)系，并且/或者它們可能較少發(fā)生基因水平轉(zhuǎn)移、轉(zhuǎn)座重組等進(jìn)化事件。

圖3 短乳桿菌家族3質(zhì)粒rep基因上游重復(fù)序列Fig. 3 Repeat sequences of rep gene flank of Family 3 plasmids in Lactobacillus brevis

3 討 論

乳桿菌屬菌種數(shù)量和菌種多樣性是乳酸菌中最豐富的菌屬[30]，同時乳桿菌屬攜帶的天然質(zhì)粒數(shù)量和多樣性也是乳酸菌中最豐富的菌屬[31]，乳桿菌屬這兩種特性是一種巧合，還是存在某種聯(lián)系，目前還是未知的。為了探究質(zhì)粒多樣性與乳桿菌生理、遺傳等生物多樣性之間是否存在聯(lián)系，本課題組對乳桿菌屬中質(zhì)粒數(shù)量最多的植物乳桿菌進(jìn)行了質(zhì)粒進(jìn)化、分類和起源研究，建立了一種準(zhǔn)確、有效的分析方法，揭示了植物乳桿菌天然質(zhì)粒系統(tǒng)進(jìn)化、分類和起源[13-14]。本研究以短乳桿菌（乳桿菌屬中質(zhì)粒多樣性僅次于植物乳桿菌）天然質(zhì)粒為研究對象，在以往質(zhì)粒Rep進(jìn)化樹分類方法研究基礎(chǔ)上，進(jìn)一步結(jié)合質(zhì)?；蚪M共線性分析方法，建立了一種更加準(zhǔn)確、有效的短乳桿菌天然質(zhì)粒分類方法，清晰、準(zhǔn)確的揭示了短乳桿菌天然質(zhì)粒的分類關(guān)系，為今后短乳桿菌天然質(zhì)?？茖W(xué)分類和理性應(yīng)用奠定了理論研究基礎(chǔ)。

本研究經(jīng)多序列比對和重復(fù)序列分析發(fā)現(xiàn)，短乳桿菌家族4質(zhì)粒復(fù)制子可以分為明顯不同的2 個亞家族類型（圖3），亞型1質(zhì)粒與參考質(zhì)粒pUCL287復(fù)制子的重復(fù)序列相似，均由兩組不同長度的串聯(lián)重復(fù)序列組成，該特征由Benachour等[26]1997年首次鑒定，而亞型2質(zhì)粒重復(fù)序列僅由一組串聯(lián)重復(fù)序列組成，這在以往有關(guān)乳桿菌質(zhì)粒的文獻(xiàn)中未發(fā)現(xiàn)報(bào)道，因此這個新的Rep_3家族復(fù)制子是否具有復(fù)制功能仍需后續(xù)實(shí)驗(yàn)驗(yàn)證。此外，亞型1質(zhì)粒復(fù)制子重復(fù)序列也并非完全一致，質(zhì)粒pL747-2、pLB925A03和pL747-6與其他質(zhì)粒重復(fù)序列也存在不同程度的差異性（圖3A）。這些分析結(jié)果表明，短乳桿菌家族4質(zhì)粒Rep及其結(jié)合位點(diǎn)（即重復(fù)序列）并非嚴(yán)格保守，顯示出較豐富的多樣性，因此短乳桿菌家族4質(zhì)粒的復(fù)制起始機(jī)制很可能因此存在一定的靈活性和魯棒性。

4 結(jié) 論

本研究通過質(zhì)粒Rep進(jìn)化樹分析，將短乳桿菌34 個編碼Rep質(zhì)粒劃分為5 個家族類型，之后通過質(zhì)粒基因組共線性分析，將短乳桿菌所有51 個質(zhì)粒劃分為6 個家族類型，以及更詳細(xì)的11 個亞家族類型，推測家族1、2質(zhì)粒屬于滾環(huán)型復(fù)制質(zhì)粒，家族3～6質(zhì)粒屬于theta型復(fù)制質(zhì)粒。基因組共線性質(zhì)粒分類方法克服了Rep進(jìn)化樹分類方法對質(zhì)粒編碼Rep的依賴和不足，Rep進(jìn)化樹和基因組共線性兩種分析方法對編碼Rep質(zhì)粒的分類結(jié)果一致，兩種方法互為補(bǔ)充和證明，因此本研究表明，Rep進(jìn)化樹和基因組共線性相結(jié)合的分析方法可以準(zhǔn)確、有效的分析和分類短乳桿菌天然質(zhì)粒，該方法是短乳桿菌天然質(zhì)粒系統(tǒng)分類的一種有效方法。

參考文獻(xiàn)：

[1]GUO T, ZHANG L, XIN Y, et al. Oxygen-inducible conversion of lactate to acetate in heterofermentative Lactobacillus brevis ATCC367[J]. Applied and Environmental Microbiology, 2017, 83(21)：1659. DOI：10.1128/AEM.01659-17.

[2]NODA M, MIYAUCHI R, DANSHIITSOODOL N, et al.Characterization and mutational analysis of a two-polypeptide bacteriocin produced by citrus lyo-derived Lactobacillus brevis 174A[J]. Biological and Pharmaceutical Bulletin, 2015, 38(12)： 1902-1909. DOI：10.1248/bpb.b15-00505.

[3]BANERJEE S P, DORA K C, CHOWDHURY S. Detection,partial purification and characterization of bacteriocin produced by Lactobacillus brevis FPTLB3 isolated from freshwater fish[J]. Journal of Food Science and Technology, 2013, 50(1)： 17-25. DOI：10.1007/s13197-011-0240-4.

[4]KIM S Y. Production of fermented kale juices with Lactobacillus strains and nutritional composition[J]. Preventive Nutrition and Food Science, 2017, 22(3)： 231-236. DOI：10.3746/pnf.2017.22.3.231.

[5]OGAWA M, SAIKI A, MATSUI Y, et al. Effects of oral intake of heat-killed Lactobacillus brevis SBC8803 (SBL88?) on dry skin conditions： a randomized, double-blind, placebo-controlled study[J].Experimental and Therapeutic Medicine, 2016, 12(6)： 3863-3872.DOI：10.3892/etm.2016.3862.

[6]ACOSTA M P, RUZAL S M, CORDO S M. S-layer proteins from Lactobacillus sp. inhibit bacterial infection by blockage of DC-SIGN cell receptor[J]. International Journal of Biological Macromolecules,2016, 92： 998-1005. DOI：10.1016/j.ijbiomac.2016.07.096.

[7]FUKAO M, OSHIMA K, MORITA H, et al. Genomic analysis by deep sequencing of the probiotic Lactobacillus brevis KB290 harboring nine plasmids reveals genomic stability[J]. PLoS One, 2013,8(3)： 60521. DOI：10.1371/journal.pone.0060521.

[8]PARK J Y, JEONG S J, SA H D, et al. Construction of a shuttle vector based on the small cryptic plasmid pJY33 from Weissella cibaria 33[J]. Plasmid, 2015, 79： 30-36. DOI：10.1016/j.plasmid.2015.03.008.

[9]WADA T, NODA M, KASHIWABARA F, et al. Characterization of four plasmids harboured in a Lactobacillus brevis strain encoding a novel bacteriocin, brevicin 925A, and construction of a shuttle vector for lactic acid bacteria and Escherichia coli[J]. Microbiology, 2009,155(5)： 1726-1737. DOI：10.1099/mic.0.022871-0.

[10]JEONG S J, PARK J Y, LEE H J, et al. Characterization of pFMBL1, a small cryptic plasmid isolated from Leuconostoc mesenteroides SY2[J]. Plasmid, 2007, 57(3)： 314-323. DOI：10.1016/j.plasmid.2006.09.003.

[11]BERGSVEINSON J, ZIOLA B. Comparative genomic and plasmid analysis of beer-spoiling and non-beer-spoiling Lactobacillus brevis isolates[J]. Canadian Journal of Microbiology, 2017, 63(12)： 970-983.DOI：10.1139/cjm-2017-0405.

[12]BERGSVEINSON J, BAECKER N, PITTET V, et al. Role of plasmids in Lactobacillus brevis BSO 464 hop tolerance and beer spoilage[J]. Applied and Environmental Microbiology, 2015, 81(4)：1234-1241. DOI：10.1128/AEM.02870-14.

[13]孫大慶, 李洪飛, 宋大巍, 等. 植物乳桿菌天然質(zhì)粒分類[J]. 食品科學(xué), 2017, 38(12)： 69-74. DOI：10.7506/spkx1002-6630-201712011.

[14]孫大慶, 李洪飛, 楊健. 植物乳桿菌編碼復(fù)制起始蛋白天然質(zhì)粒的系統(tǒng)進(jìn)化分析[J]. 微生物學(xué)通報(bào), 2017, 44(5)： 1047-1055.DOI：10.13344/j.microbiol.china.160940.

[15]FRANCIA M V, HAAS W, WIRTH R, et al. Completion of the nucleotide sequence of the Enterococcus faecalis conjugative virulence plasmid pAD1 and identification of a second transfer origin[J].Plasmid, 2001, 46(2)： 117-127. DOI：10.1006/plas.2001.1533.

[16]THOMPSON J K, COLLINS M A. Completed sequence of plasmid pIP501 and origin of spontaneous deletion derivatives[J]. Plasmid,2003, 50(1)： 28-35. DOI：10.1016/S0147-619X(03)00042-8.

[17]TAMURA K, STECHER G, PETERSON D, et al. MEGA6： molecular evolutionary genetics analysis version 6.0[J]. Molecular Biology and Evolution, 2013, 30(12)： 2725-2729. DOI：10.1093/molbev/mst197.

[18]DARLING A C E, MAU B, BLATTNER F R, et al. Mauve： multiple alignment of conserved genomic sequence with rearrangements[J].Genome Research, 2004, 14(7)： 1394-1403. DOI：10.1101/gr.2289704.

[19]SEERY L T, NOLAN N C, SHARP P M, et al. Comparative analysis of the pC194 group of rolling circle plasmids[J]. Plasmid, 1993, 30(3)：185-196. DOI：10.1006/plas.1993.1051.

[20]MOSCOSO M, ERITJA R, ESPINOSA M. Initiation of replication of plasmid pMV158： mechanisms of DNA strand-transfer reactions mediated by the initiator RepB protein[J]. Journal of Molecular Biology, 1997, 268(5)： 840-856. DOI：10.1006/jmbi.1997.1012.

[21]WANG Z, YUAN Z, HENGGE U R. Processing of plasmid DNA with ColE1-like replication origin[J]. Plasmid, 2004, 51(3)： 149-161.DOI：10.1016/j.plasmid.2003.12.002.

[22]CESARENI G, HELMER-CITTERICH M, CASTAGNOLI L. Control of ColE1 plasmid replication by antisense RNA[J]. Trends in Genetics,1991, 7(7)： 230-235. DOI：10.1016/0168-9525(91)90370-6.

[23]JUNG Y H, LEE Y. RNases in ColE1 DNA metabolism[J]. Molecular Biology Reports, 1995, 22(2/3)： 195-200. DOI：10.1007/BF00988728.

[24]BENACHOUR A, FRèRE J, NOVEL G. pUCL287 Plasmid from Tetragenococcus halophila (Pediococcus halophilus) ATCC 33315 represents a new theta-type replicon family of lactic acid bacteria[J].FEMS Microbiology Letters, 1995, 128(2)： 167-175. DOI：10.1111/j.1574-6968.1995.tb07518.x.

[25]BENACHOUR A, FRèRE J, BOUTIBONNES P, et al.Characterization and replication mode determination of the minimal replicon of Tetragenococcus halophila ATCC33315 plasmid pUCL287[J]. Biochimie, 1995, 77(11)： 868-874. DOI：10.1016/0300-9084(95)90005-5.

[26]BENACHOUR A, FRèRE J, FLAHAUT S, et al. Molecular analysis of the replication region of the theta-replicating plasmid pUCL287 from Tetragenococcus (Pediococcus) halophilus ATCC33315[J].Molecular and General Genetics MGG, 1997, 255(5)： 504-513.DOI：10.1007/s004380050523.

[27]WEAVER K E, KWONG S M, FIRTH N, et al. The RepA_N replicons of Gram-positive bacteria： a family of broadly distributed but narrow host range plasmids[J]. Plasmid, 2009, 61(2)： 94-109.DOI：10.1016/j.plasmid.2008.11.004.

[28]SCHUMACHER M A, TONTHAT N K, KWONG S M, et al.Mechanism of staphylococcal multiresistance plasmid replication origin assembly by the RepA protein[J]. Proceedings of the National Academy of Sciences, 2014, 111(25)： 9121-9126. DOI：10.1073/pnas.1406065111.

[29]BRANTL S. Antisense-RNA mediated control of plasmid replicationpIP501 revisited[J]. Plasmid, 2015, 78： 4-16. DOI：10.1016/j.plasmid.2014.07.004.

[30]SUN Z, HARRIS H M B, MCCANN A, et al. Expanding the biotechnology potential of Lactobacilli through comparative genomics of 213 strains and associated genera[J]. Nature Communications, 2015,6： 8322. DOI：10.1038/ncomms9322.

[31]孫大慶, 李洪飛, 宋大巍, 等. 乳桿菌屬天然質(zhì)粒研究進(jìn)展[J]. 食品科學(xué), 2015, 36(11)： 251-255. DOI：10.7506/spkx1002-6630-201511047.