BPI基因生物信息學分析及其在不同豬群體中的遺傳變異情況

吳華莉 涂尾龍 曹建國 張鶯鶯 王洪洋 談永松

摘要：【目的】明確豬殺菌性/通透性增加蛋白（BPI）基因中可能存在的SNP位點，并探討BPI基因外顯子10在上海地區5個豬群體中的多態性及其與抗病相關的優勢基因型，為揭示BPI基因在仔豬腹瀉抗病機制中的作用及開展豬抗病育種提供理論依據。【方法】采用在線生物信息學軟件分別從堿基序列、SNP位點和蛋白氨基酸序列等角度比對分析豬BPI基因全長cDNA序列與電子克隆序列的差異，采用PCR-RFLP檢測分析豬BPI基因外顯子10的多態性，并檢測BPI基因外顯子10在杜洛克、長白、大白、申農和梅山豬群體中的基因型及基因頻率。【結果】豬BPI基因電子克隆編碼區（CDS）序列長度1452 bp，編碼483個氨基酸，存在44個SNP位點，其中22個為同義SNP位點、22個為非同義SNP位點。豬BPI基因全長cDNA序列（EF436278.1和FJ810853.1）與電子克隆序列相比存在14處堿基差異，而對應的翻譯蛋白序列存在4處氨基酸差異。豬BPI基因外顯子10存在AA、BB和AB 3種基因型，各基因型在杜洛克豬、長白豬、大白豬、申農豬和梅山豬群體中的分布情況各不相同。在杜洛克豬群體中AA基因型和AB基因型呈中度多態分布（0.404和0.416），BB基因型檢出比例較低（0.180）;在長白豬和申農豬群體中均未檢測出AA基因型，BB基因型檢出比例較高，而AB基因型檢出比例較低;在大白豬和梅山豬群體中均以BB基因型檢出比例較高，AB基因型次之，AA基因型檢出比例最低。總體來看，除了杜洛克豬群體的A等位基因頻率較高外，其他4個豬群體均以B等位基因頻率較高。【結論】豬BPI基因外顯子10 AA基因型頻率在杜洛克豬群體中最高，在大白豬和梅山群體中較低，在長白豬和申農豬群體中未檢出，即AA基因型在杜洛克豬群體中比例高與其仔豬腹瀉發病率低的情況一致。因此，在大白豬和梅山豬育種過程中需提高AA基因型比例，而在長白豬和申農豬育種過程中需引進AA基因型個體，降低BB和AB基因型比例，以增加仔豬的抗病能力。

關鍵詞： 豬;殺菌性/通透性增加蛋白（BPI）;SNP位點;基因型;多態性

中圖分類號： S828.2? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2019）03-0446-08

0 引言

【研究意義】斷奶仔豬腹瀉一直是困擾養豬業健康發展的難題，主要引起仔豬消瘦，影響其生長指標，嚴重者甚至引起死亡。引起仔豬腹瀉原因有很多，其中大腸桿菌和沙門氏菌是最主要的病原菌（鄭龍龍等，2015）。殺菌性/通透性增加蛋白（Bactericidal/permeability-increasing protein，BPI）是脂多糖結合蛋白家族成員之一，定位于胞漿顆粒細胞中，為多形核中性粒細胞（PMNs）的主要抗菌成分（Canny and Levy，2008）。BPI在機體識別和調控革蘭氏陰性菌致病因子方面發揮重要作用，可中和內毒素而降低炎癥反應，其功能性N端還可殺死細菌，促進中性粒細胞對細菌的吞噬作用（Akin et al.，2011）。Tuggle等（2003）研究表明，BPI基因外顯子4和外顯子10的多態性與豬沙門氏菌易感性有關;朱璟等（2011）研究認為，BPI基因外顯子10 AA基因型可能具有較高的免疫應答能力，是提高對大腸桿菌抗性的有利基因型。因此，明確BPI基因外顯子在不同豬群體中的多態性對控制仔豬腹瀉具有重要意義。【前人研究進展】自Weiss等（1978）最先從PMNs中分離獲得人類BPI以來，目前已從人類及豬、牛、小鼠、獼猴等哺乳動物中克隆獲得BPI基因（陳薇等，1997;高恒，2008），其中，豬BPI基因定位于第17號染色體上，含有15個外顯子和14個內含子（魏麟，2012）。吳正常等（2013a）通過比較13個物種的BPI基因編碼區（CDS）序列，結果發現除獼猴外，人類及豬、鼠等動物均具有BPI1和BPI2兩種保守的結構域。人類與豬的BPI已被證實具有中和內毒素及抗革蘭氏陰性菌的作用（周紅等，1999;Levy et al.，2003）。此外，Tuggle等（2003）檢測到BPI基因外顯子4和外顯子10分別存在AvaⅡ和HPaⅡ酶切多態性，其基因型與豬沙門氏菌的易感性有關，即BPI基因為抗病育種候選基因。袁樹楷（2007）研究發現，榮昌豬BPI基因外顯子3第397位存在的堿基突變（G→A），引起氨基酸由精氨酸（Arg）替換為谷氨酸（Gln），且該突變位點可能對榮昌豬BPI蛋白功能及機體天然免疫力有重要影響;而榮昌豬BPI基因外顯子4存在4個SNP位點（T512C、G551T、C563T和G599A），說明該片段具有較大的遺傳選擇潛力。陳平（2015）研究證實，豬源BPI N端不僅能殺滅大腸桿菌、沙門氏菌等革蘭氏陰性菌，還對L-型化膿鏈球菌等革蘭氏陽性菌具有抗菌活性。吳正常等（2015）通過篩選靶向干擾豬BPI基因表達的siRNA，以期從細胞水平研究BPI基因對豬腸道革蘭氏陰性菌抗感染的作用機制。Bülow等（2018）研究發現，BPI通過結合革蘭氏陽性菌的脂肽、脂蛋白及脂磷壁酸而產生殺菌作用。【本研究切入點】目前，在豬場生產實踐中發現長白豬及其培育品種申農豬的仔豬腹瀉發病率較高，但是否與BPI基因外顯子多態性差異有關尚有待進一步探究。【擬解決的關鍵問題】通過分析豬BPI基因序列信息尋找可能存在變異位點（SNP位點），并檢測分析BPI基因外顯子10在上海地區5個豬群體中的多態性及其與抗病相關的優勢基因型，為揭示BPI基因在仔豬腹瀉抗病機制中的作用及開展豬抗病育種提供理論依據。

1 材料與方法

1. 1 試驗材料

供試的5個豬群體樣品中，杜洛克豬（89頭）、大白豬（154頭）和長白豬（79頭）由上海祥欣畜禽有限公司提供，梅山豬（130頭）由上海豬狀元畜牧綜合場提供，申農豬（75頭）由上海警備區富民農場提供，共計527頭，所取樣品均為豬耳組織。所采樣品置于預先添加75%酒精的1.5 mL離心管中，以冰盒保存帶回實驗室，-20 ℃保存備用。同時從NCBI網站（https：//www.ncbi.nlm.nih.gov/gene/100286800）下載豬BPI基因的4條電子克隆序列和2條全長cDNA克隆序列，具體信息如表1所示。

1. 2 試驗方法

1. 2. 1 生物軟件分析 采用BioEdit、Lasergene和DNAMAN對下載的6條BPI基因序列進行比對分析，選取共有的CDS序列，其長度為1452 bp;同時對其編碼蛋白氨基酸序列進行比對分析。

1. 2. 2 DNA提取 采用AXYGEN組織DNA提取試劑盒提取豬耳組織基因組DNA，并稀釋至30 ng/mL， -20 ℃保存備用。

1. 2. 3 引物設計合成及PCR擴增 參照朱璟等（2011）的方法設計擴增引物（上游引物：5'-CCCAAC ATGGAGATGCAGTTC-3';下游引物：5'-CAATGAA TCAATGAGCACACC-3'），引物由生工生物工程（上海）股份有限公司合成。PCR反應體系50.0 μL：Master Mix 25.0 μL，BPI-F和BPI-R各1.0 μL，DNA模板1.0 μL，ddH2O 22.0 μL。擴增程序：94 ℃預變性5 min;94 ℃ 30 s，58 ℃ 30 s，72 ℃ 30 s，進行30個循環;72 ℃延伸5 min。PCR產物采用1.0%瓊脂糖凝膠電泳進行檢測。

1. 2. 4 PCR-RFLP檢測 PCR擴增產物經1.0%瓊脂糖凝膠電泳檢測后進行酶切，酶切體系16.0 μL： PCR擴增產物4.0 μL，10×Buffer 2.0 μL，HpaⅡ酶0.5 μL，ddH2O 9.5 μL;酶切時間約30 min。

1. 2. 5 基因型分析 PCR擴增獲得的目的基因片段存在HpaⅡ酶切位點，采用內切酶HpaⅡ進行酶切后可產生3種基因型：AA型（445 bp）、BB型（304和142 bp）和AB型（445、304和142 bp）。根據酶切結果，不同基因型分別選取3個樣品進行測序。

1. 2. 6 基因型頻率和基因頻率計算 等位基因是指位于一對同源染色體相同位置上控制著相對性狀的一對基因，以A和B表示。依據群體遺傳學計算公式，分別計算等位基因頻率及其基因型頻率（王旭等，2017;王利平等，2018）。

2 結果與分析

2. 1 豬BPI基因CDS序列分析結果

豬BPI基因CDS序列長度1452 bp，編碼483個氨基酸，其信號肽區域長度為27個氨基酸（第1～27位氨基酸），功能區BPI1長度為226個氨基酸（第35～260位氨基酸），功能區BPI2長度為204個氨基酸（第275～478位氨基酸）。根據NCBI提供的SNP GeneView信息可知，豬BPI基因CDS序列存在44個SNPs位點，其中22個為同義SNP位點、22個為非同義SNP位點。

2. 2 豬BPI基因電子克隆序列與全長cDNA序列的SNP位點比對分析結果

DANMAN分析結果顯示，豬BPI基因全長cDNA序列（EF436278.1和FJ810853.1）與電子克隆序列相比存在14處堿基差異，具體信息表2。其中，4條電子克隆序列完全一致。從NCBI數據庫（https：//www.ncbi.nlm.nih.gov/gene/100286800）下載豬BPI基因DNA序列（NC_010459），然后根據其序列特征劃分BPI基因外顯子區域，15個外顯子長度分別為130、115、129、162、64、64、92、177、60、168、68、43、64、77和39 bp。電子克隆序列與全長cDNA序列的堿基差異分別位于：外顯子2（C401T）、外顯子3（A505G）、外顯子4（C582T、G607C、C615T、T654G、C666T和G702A）、外顯子6（C855T）、外顯子7（C925T）、外顯子8（G1005T）、外顯子10（A1258G）、外顯子13（T1533C）和外顯子15（C1647T）。

2. 3 豬BPI基因電子克隆蛋白序列與全長cDNA翻譯蛋白序列比對分析結果

豬BPI基因全長cDNA翻譯蛋白序列與電子克隆蛋白序列存在4處氨基酸差異，分別是第118位精氨酸突變為谷氨酰胺（Arg→Gln）、第154位精氨酸突變為甘氨酸（Arg→Gly）、第260位蘇氨酸突變為絲氨酸（Thr→Ser）和第371位丙氨酸突變為蘇氨酸（Ala→Thr）。其中，BPI1區域有3處氨基酸差異，BPI2區域僅有1處氨基酸差異。

2. 4 豬BPI基因PCR擴增結果

采用1.0%瓊脂糖凝膠電泳檢測PCR擴增產物，結果發現目的基因條帶大小為445 bp，且目的條帶單一明亮、無雜帶，符合后續的酶切要求。

2. 5 豬BPI基因測序結果

豬BPI基因PCR擴增產物測序結果，其序列中存在Hap II酶切位點（C/CGG），在測序結果中以加粗下劃線標注酶切位點位置。

2. 6 各基因型測序分析結果

根據Hap II酶切電泳結果，3種基因型各選擇3個樣品進行測序。AA基因型突變位點處堿基為T時未形成Hap II（C/CGG）酶切位點，酶切后PCR片段長度不變，仍為445 bp;BB基因型突變位點處堿基為G時形成Hap II（C/CGG）酶切位點，酶切后PCR片段被分為304和142 bp 2個片段;AB基因型為雜合子，酶切后PCR片段被分為445、304和142 bp 3個片段。

2. 7 PCR-RFLP檢測結果

采用1.5%瓊脂糖凝膠電泳檢測Hap II酶切結果，結果顯示獲得的目的條帶與預期結果一致。其中，AA基因型為1條條帶，BB基因型為2條條帶，AB基因型為3條條帶。說明待檢測樣品存在BPI基因外顯子10的3種基因型。當突變位點為T時未產生酶切位點，酶切產物長度445 bp;當突變位點為G時產生酶切位點，酶切產物長度為304和142 bp;當突變位點為T/G雜合子時，酶切產物長度分別為445、304和142 bp。

2. 8 BPI基因外顯子10在5個豬群體中的基因型及基因頻率

由表3可知，BPI基因外顯子10在杜洛克豬、長白豬、大白豬、申農豬和梅山豬群體中各基因型的分布情況各不相同。在杜洛克豬群體中AA和AB基因型呈中度多態分布，BB基因型檢出比例較低;在長白豬和申農豬群體中均未檢測出AA基因型，BB基因型檢出比例較高，而AB基因型檢出比例較低;在大白豬和梅山豬群體中均以BB基因型檢出比例較高，AB基因型次之，AA基因型檢出比例最低。總體來看，除了杜洛克豬群體的A等位基因頻率較高外，其他4個豬群體均以B等位基因頻率較高。

3 討論

BPI存在于多形核中性粒細胞（PMNs）中，能中和細菌釋放的脂多糖（LPS），并可激活磷脂，增強細胞膜通透性，導致細菌死亡，尤其對革蘭氏陰性菌具有廣泛的殺滅作用，在機體體內被稱為超級抗生素（朱璟，2012）。BPI的N端與病原菌LPS結合，而C端結合于吞噬細胞表面，增強中性粒細胞對細菌的吞噬作用（Nishimura et al.，2001;高恒等，2009）。本研究通過分析NCBI上已公布的豬BPI基因序列，發現其全長cDNA克隆序列與電子克隆序列相比存在14處堿基差異，分別位于：外顯子2（C401T）、外顯子3（A505G）、外顯子4（C582T、G607C、C615T、T654G、C666T和G702A）、外顯子6（C855T）、外顯子7（C925T）、外顯子8（G1005T）、外顯子10（A1258G）、外顯子13（T1533C）和外顯子15（C1647T）;而對應的翻譯蛋白序列存在4處氨基酸差異，分別是第118位Arg→Gln、第154位Arg→Gly、第260位Thr→Ser和第371位Ala→Thr。由于電子克隆序列是利用生物信息學技術組裝延伸ESTs序列所獲得，因此需采用全長cDNA克隆加以驗證。二者間存在的差異可能是突變位點所引起，其中在外顯子4上存在6處堿基差異和4個SNPs位點，與袁樹楷（2007）的研究結果一致，其余SNP位點則需在不同豬種群體進行差異位點檢測分析，以確定其遺傳變異規律。

目前，已有研究報道豬BPI基因外顯子1、外顯子3、外顯子4和外顯子10均存在SNP位點，且證實外顯子4和外顯子10的多態性與沙門氏菌易感性相關（袁樹楷，2007;曹曉華，2008;吳正常等，2013b;Wu et al.，2015）。朱璟等（2011）研究發現，豬BPI基因外顯子10 AA基因型是有利基因型，對大腸桿菌F18菌株的抗性顯著或極顯著高于AB和BB基因型個體，且證實這種抗性與豬BPI基因在空腸和十二指腸上調表達有關。本研究結果顯示，BPI基因外顯子10 AA基因型在外來豬種（杜洛克）群體中的比例較高，在地方品種（梅山豬）和大白豬群體中的比例較低，而在長白豬和申農豬群體中均未檢測出AA基因型，與朱璟（2012）的研究結果一致。申農豬是由長白×（大白×二花臉）經橫交固定、世代選育形成的母本新品系，其特點是耐粗飼、母豬體型較小、耗料較少、飼養管理要求不高。以其為母本、大白豬為父本雜交獲得的子代在生長指標、肉色、肌內脂肪等方面均優于親本（胡志剛等，2010;談永松等，2011），但該雜交品系在育種過程中常遇仔豬腹瀉頻發的問題，是否與BPI基因外顯子10 AA基因型在檢測群體中缺失有關還有待進一步探究。

本研究通過檢測豬BPI基因在5個豬群體的遺傳變異情況，一方面可尋找疾病易感基因型的遺傳標記應用于育種實踐，對提高豬的經濟價值具有重要意義;另一方面，基因型分布可反映豬群體遺傳結構穩定。因此，在申農豬育種過程中應重點關注BPI基因外顯子10 AA基因型比例的增加情況。此外，需進一步結合5個豬群體的仔豬腹瀉情況進行相關分析，從而判定仔豬腹瀉與豬BPI基因外顯子10基因型是否相關，并確定是否還存在控制疾病發生發展的其他協同基因。

4 結論

豬BPI基因外顯子10存在AA、BB和AB 3種基因型，AA基因型頻率在杜洛克豬群體中最高，在大白豬和梅山群體中較低，在長白豬和申農豬群體中未檢出，即AA基因型在杜洛克豬群體中比例高與其仔豬腹瀉發病率低的情況一致。因此，在大白豬和梅山豬育種過程中需提高AA基因型比例，而在長白豬和申農豬育種過程中需引進AA基因型個體，降低BB和AB基因型比例，以增加仔豬的抗病能力。

參考文獻：

曹曉華. 2008. 四個豬種間BPI蛋白基因外顯子3和4的SNP分析[D]. 重慶：西南大學. [Cao X H. 2008. SNP analysis of exon 3 and exon 4 of bactericidal/permeability-increasing protein gene in four swine breed[D]. Chongqing：Southwest University.]

陳平. 2015. 豬源殺菌通透性增加蛋白N端基因的克隆原核表達及抗菌活性測定[D]. 雅安：四川農業大學. [Chen P. 2015. A porcine bactericidal permeability increase protein N-terminal gene cloning and prokaryotic expression of antimicrobial activity[D]. Ya’an：Sichuan Agricultural University.]

陳薇，黃策，王海濤，傅玲，俞曉峰，徐靜，杜桂鑫. 1997. 中國人殺菌蛋白氨基端基因的克隆和序列分析[J]. 微生物學通報，24（6）：350-353. [Chen W，Huang C，Wang H T，Fu L，Yu X F，Xu J，Du G X. 1997. cDNA cloning and sequencing of the N-terminal gene of BPI from Chinese[J]. Microbiology，24（6）：350-353.]

高恒. 2008. 豬牛殺菌/通透性增加蛋白N端cDNA的克隆和原核表達[D]. 合肥：安徽農業大學. [Gao H. 2008. Cloning and prokaryotic expession of bactericidal/permeability-increasing proteins N-terminal cDNA in pig and cow[D]. Hefei：Anhui Agricultural University.]

高恒，溫納相，祁克宗. 2009. 殺菌通透性增加蛋白研究概況[J]. 中國獸醫雜志，45（11）：50-52. [Gao H，Wen N X，Qi K Z. 2009. Study on bactericidal/permeability-increasing protein[J]. Chinese Journal of Veterinary Medicine，45（11）：50-52.]

胡志剛，談永松，王林云. 2010. 申農1號雜交豬生長及胴體性能實驗[J]. 國外畜牧學（豬與禽），30（4）：67-68. [Hu Z G，Tan Y S，Wang L Y. 2010. Shennong hybridization pig growth and carcass performance experiment[J]. Animal Science Abroad（Pigs and Poultry），30（4）：67-68.]

談永松，胡志剛，涂尾龍，曹建國，高勤學. 2011. 申農1號豬不同雜交組合試驗報告[J]. 國外畜牧學（豬與禽），31（3）：75-80. [Tan Y S，Hu Z G，Tu W L，Cao J G，Gao Q X. 2011. Experimental report on different hybrid combinations of Shennong pig[J]. Animal Science Abroad（Pigs and Poultry），31（3）：75-80.]

王利平，席冬梅，熊和麗，李國治，李鵬，鄧衛東. 2018. 中甸犏牛SLC11A1基因3'非翻譯區多態位點研究及其與抗病性關聯的評估[J]. 河南農業科學，47（9）：131-137. [Wang L P，Xi D M，Xiong H L，Li G Z，Li P，Deng W D. 2018. Study on polymorphisms of 3' untranslated region of SLC11A1 gene and evaluation of disease resistance in Zhongdian yakow[J]. Journal of Henan Agricultural Scien-ces，47（9）：131-137.]

王旭，白俊艷，楊又兵，雷雪芹，龐有志，李宏偉，王歡玲，湯昱琳，邵俊紅，陰懷源. 2017. 小尾寒羊MyoG外顯子2和MyoD 5'側翼區的多態性研究[J]. 河南農業科學，46（6）：138-141. [Wang X，Bai J Y，Yang Y B，Lei X Q，Pang Y Z，Li H W，Wang H L，Tang Y L，Shao J H，Yin H Y. 2017. Study on the polymorphism of MyoG exon 2 and MyoD 5' flanking region of small tail han sheep[J]. Journal of Henan Agricultural Sciences，46（6）：138-141.]

魏麟. 2012. 豬LBP和BPI基因多態性及其蛋白質N端功能研究[D]. 長沙：湖南農業大學. [Wei L. 2012. Study on the polymorphism of porcine LBP and BPI genes and their N-terminal proteins’ function[D]. Changsha：Hunan A-gricultural University.]

吳正常，蘇先敏，王瑾，鄭先瑞，吳圣龍，包文斌. 2013a. 13個物種BPI基因編碼區生物信息學分析[J]. 中國畜牧獸醫，40（3）：28-32. [Wu Z C，Su X M，Wang J，Zheng X R，Wu S L，Bao W B. 2013a. Bioinformatics analysis on the coding region of BPI gene among thirteen different species[J]. China Animal Husbandry & Veterinary Medicine，40（3）：28-32.]

吳正常，王靖，朱世平，趙喬輝，劉璐，訾臣，吳圣龍，包文斌. 2013b. 豬BPI基因部分cSNPs檢測及其對蛋白結構功能的影響[J]. 畜牧獸醫學報，44（7）：1000-1007. [Wu Z C，Wang J，Zhu S P，Zhao Q H，Liu L，Zi C，Wu S L，Bao W B. 2013b. The partial cSNPs detection of pig BPI gene and its influence on protein structure and function[J]. Acta Veterinaria et Zootechnia Scienca，44（7）：1000-1007.]

吳正常，殷學梅，夏日煒，孫壽永，朱國強，吳圣龍，包文斌. 2015. 豬殺菌/通透性增加蛋白基因siRNA載體構建及干擾效果評價[J]. 畜牧獸醫學報，46（3）：491-496. [Wu Z C，Yin X M，Xia R W，Sun S Y，Zhu G Q，Wu S L，Bao W B. 2015. Construction and evaluation of porcine bactericidal/permeabilit-increasing protein gene（BPI） siRNA expression vector[J]. Acta Veterinaria et Zootechnia Scien-ca，46（3）：491-496.]

袁樹楷. 2007. 榮昌豬BPI基因全長cDNA克隆及SNP分析[D]. 重慶：西南大學. [Yuan S K. 2007. Cloning and SNP analysis of bactericidal/permeability-increasing protein gene in Rongchang pig[D]. Chongqing：Southwest University.]

鄭龍龍，朱玲云，陳關雄，劉萍丹，唐寧，趙德育，劉佳，楊貴樹，白衛兵. 2015. 斷奶仔豬腹瀉病因分析[J]. 中國畜牧獸醫，42（10）：2772-2778. [Zheng L L，Zhu L Y，Chen G X，Liu P D，Tang N，Zhao D Y，Liu J，Yang G S，Bai W B. 2015. Cause analysis of weaning piglets diarrhea[J]. China Animal Husbandry & Veterinary Medicine，42（10）：2772-2778.]

周紅，袁建成，周立新，鄭江，蕭光夏. 1999. 殺菌性/通透性增加蛋白中和內毒素的作用的體內外研究[J]. 中華醫學雜志，79（4）：304-305. [Zhou H，Yuan J C，Zhou L X，Zheng J，Xiao G X. 1999. In vitro and in vivo studies of bactericidal/permeabilityincreasing protein neutralize endotoxins[J]. National Medical Journal of China，79（4）：304-305.]

朱璟. 2012. BPI基因對斷奶仔豬F18大腸桿菌抗性的功能分析[D]. 揚州：揚州大學. [Zhu J. 2012. Effect of BPI gene involved in resistance to E. coli F18 infection of weaning piglets[D]. Yangzhou：Yangzhou University.]

朱璟，潘章源，訾臣，劉璐，杜子棟，黃小國，吳圣龍，包文斌. 2011. 11個豬群體BPI基因第10外顯子HpaⅡ遺傳變異分析[J]. 中國畜牧雜志，47（23）：14-17. [Zhu J，Pan Z Y，Zi C，Liu L，Du Z D，Huang X G，Wu S L，Bao W B. 2011. BPI gene exon 10 HpaⅡ genetic variation analysis in eleven pig groups[J]. Chinese Journal of Animal Scien-ce，47（23）：14-17.]

Akin H，Tahan G，Türe F，Eren F，Atu?a O，Tahand V，Hamzao?lu I，?meryüza N，T?züna N，Hamzaoglu H O. 2011. Associa-tion between bactericidal/permeability increasing protein （BPI） gene polymorphism（Lys216Glu） and in?ammatory bowel disease[J]. Journal of Crohn’s and Colitis，5（1）：14-18.

Bülow S，Zeller L，Werner M，Toelge M，Holzinger J，Entzian C，Schubert T，Waldow F，Gisch N，Hammerschmidt S，Gessner A. 2018. Bactericidal/permeability-increasing protein is an enhancer of bacterial lipoprotein recognition[J].Frontiers in Immunology，9：2768. doi：10.3389/fimmu. 2018.02768.

Canny G，Levy O. 2008. Bactericidal/permeability-increasing protein（BPI） and BPI homologs at mucosal sites[J]. Trends in Immunology，29（11）：541-547.

Levy O，Canny G，Serhan C N，Colgan S P. 2003. Expression of BPI（bactericidal/permeability-increasing protein） in human mucosal epithelia[J]. Biochemical Society Tran-sactions，31（4）：795-800.

Nishimura H，Gogami A，Miyagawa Y，Nanbo A，Murakami Y，Baba T，Nagasawa S. 2001. Bactericidal/permeability-increasing protein promotescomplement activation for neutrophil-mediated phagocytosison bacterial surface[J].Immunology，103（4）：519-525.

Tuggle C K，Stabel T J， Shi X W，Mellencamp M A. 2003. Genetic markers for improved disease resistance in animals（BPI）. United States，Kind Code：A1，Patent Application： 20030104424[P]. 2003-11-25.

Weiss J，Elsbach P，Olsson I，Odeberg H. 1978. Purification and characterization of a potent bactericidal and membrane active protein from the granules of human polymorphonuclear leukocytes[J]. The Journal Biological Chemi-stry，253（8）：2664-2672.

Wu Z C，Liu Y，Zhao Q H，Zhu S P，Huo Y J，Zhu G Q，Wu S L，Bao W B. 2015. Association between polymorphisms in exons4 and 10 of the BPI gene and immune indicesin Sutai pigs[J]. Genetics and Molecular Research，14（2）：6048-6058.

（責任編輯 蘭宗寶）