李季 黃天帶 蔡海濱 華玉偉 黃華孫

摘 要 以橡膠樹(shù)轉(zhuǎn)基因材料為研究對(duì)象，首次基于載體左右邊界序列設(shè)計(jì)的特異引物進(jìn)行載體骨架序列的初步PCR檢測(cè)，結(jié)果表明，15個(gè)轉(zhuǎn)基因陽(yáng)性株系中有5個(gè)和2個(gè)株系分別含有左側(cè)和右側(cè)的邊界序列，占33.3%和13.3%，結(jié)合這些株系PCR產(chǎn)物的序列分析結(jié)果進(jìn)一步證實(shí)除T-DNA區(qū)以外的載體骨架序列確實(shí)整合到這些株系的植物基因組中。同時(shí)對(duì)載體骨架序列存在的弊端及可能的解決途徑做了相應(yīng)的分析。

關(guān)鍵詞 橡膠樹(shù)；轉(zhuǎn)基因；載體骨架序列；PCR

中圖分類號(hào) S794.1 文獻(xiàn)標(biāo)識(shí)碼 A

傳統(tǒng)觀點(diǎn)普遍認(rèn)為在遺傳轉(zhuǎn)化過(guò)程中，T-DNA邊界之外的載體序列或載體骨架序列（vector backbone sequences）不可能轉(zhuǎn)移至宿主植物細(xì)胞中[1]，然而，針對(duì)多種物種的轉(zhuǎn)基因整合模式的研究結(jié)果表明，T-DNA區(qū)域以外的載體骨架序列經(jīng)常插入到植物基因組中[2-8]。不過(guò)有關(guān)骨架序列的整合機(jī)制還有待于深入的研究。

研究結(jié)果表明，載體骨架序列整合進(jìn)植物染色體中可能會(huì)產(chǎn)生重組熱點(diǎn)或?qū)е庐惓Ｖ亟M，進(jìn)而影響目標(biāo)基因的表達(dá)[9]，轉(zhuǎn)基因后代中也易發(fā)生目標(biāo)位點(diǎn)的丟失與表達(dá)沉默等[10-11]。在轉(zhuǎn)基因植物批準(zhǔn)商業(yè)化過(guò)程中，決策部門要求提供插入到植物染色體上DNA的全部序列信息，且不能含有載體骨架序列[6]；而且非目的基因外源DNA片段的存在是否造成生物安全性尚無(wú)定論，這些都給轉(zhuǎn)基因過(guò)程中載體骨架序列的存在產(chǎn)生了負(fù)面影響，也成為了目前植物基因工程研究的熱點(diǎn)之一。

橡膠樹(shù)轉(zhuǎn)基因起步相對(duì)較晚，目前大多數(shù)的研究主要集中在遺傳轉(zhuǎn)化過(guò)程的優(yōu)化、目的基因的整合及表達(dá)情況等方面，對(duì)于載體骨架序列的發(fā)現(xiàn)及其整合現(xiàn)象還未見(jiàn)相關(guān)的研究報(bào)道。

本研究在對(duì)橡膠樹(shù)轉(zhuǎn)基因植株進(jìn)行陽(yáng)性鑒定的同時(shí)，首次根據(jù)載體左右邊界序列設(shè)計(jì)特異引物對(duì)載體骨架序列進(jìn)行了初步的PCR分析，結(jié)果表明在遺傳轉(zhuǎn)化過(guò)程中，載體骨架序列確實(shí)隨著T-DNA區(qū)一同轉(zhuǎn)移至橡膠樹(shù)基因組中，這也是目前在橡膠樹(shù)遺傳轉(zhuǎn)化研究中首次報(bào)道。由于載體骨架序列的存在可能會(huì)影響基因的表達(dá)及在批準(zhǔn)商業(yè)化時(shí)不能存在載體骨架序列，因此，有必要對(duì)轉(zhuǎn)基因植物的載體骨架序列進(jìn)行研究。筆者對(duì)橡膠樹(shù)轉(zhuǎn)基因植株中的載體骨架序列進(jìn)行初步的分析，并對(duì)其消除的意義及解決措施做了相應(yīng)的討論，以期在后續(xù)的實(shí)驗(yàn)過(guò)程中采取措施加以剔除，達(dá)到僅保留目標(biāo)基因片段與植物染色體融合，最終降低后續(xù)的相關(guān)風(fēng)險(xiǎn)。

1 材料與方法

1.1 材料

本研究使用的實(shí)驗(yàn)材料為本課題組于2007～2009年以橡膠樹(shù)（Hevea brasiliensis）花藥愈傷組織為外植體，通過(guò)農(nóng)桿菌介導(dǎo)的遺傳轉(zhuǎn)化方法獲得的轉(zhuǎn)基因植株穩(wěn)定期葉片。26個(gè)轉(zhuǎn)基因株系于2011年11月定植于中國(guó)熱帶農(nóng)業(yè)科學(xué)院橡膠研究所國(guó)家橡膠樹(shù)種質(zhì)資源圃，每個(gè)株系定植3株，共成活59株；同時(shí)以未轉(zhuǎn)化花藥胚為初生胚通過(guò)次生體胚發(fā)生增殖后再生的3株組培苗作為對(duì)照株系，總計(jì)62株。所有材料均經(jīng)過(guò)載體上的抗生素特異引物檢測(cè)以及GUS染色分析的方法確定其是否為陽(yáng)性單株，共檢測(cè)出15個(gè)轉(zhuǎn)基因陽(yáng)性株系（37份植株），具體檢測(cè)結(jié)果見(jiàn)文獻(xiàn)報(bào)道[12]。Tj200863、TP0828、T2009519、Tj2008171、T2009626和T2009401六個(gè)株系經(jīng)southern檢驗(yàn)確定目的基因已整合進(jìn)基因組中（結(jié)果未發(fā)表），這62株均作為載體骨架序列的分析材料。

PCR反應(yīng)所需要的Taq酶購(gòu)于Fermentas公司，100 bp DNA Ladder購(gòu)于北京全式金生物技術(shù)有限公司。瓊脂糖凝膠回收試劑盒E.Z.N.A.TM Gel Extraction Kit購(gòu)于OMEGA公司，pMDR-19T vector（TaKaRa）購(gòu)于大連寶生物公司。

1.2 方法

1.2.1 DNA制備 葉片（穩(wěn)定期）基因組DNA的制備方法參照文獻(xiàn)報(bào)道的方法[12]。（1）取兩個(gè)離心管蓋大小的葉片直接置于研缽中，加入65 ℃預(yù)熱的600 μL裂解液，研磨至勻漿狀，轉(zhuǎn)移勻漿至2.0 mL的離心管內(nèi)；（2）加入等體積的24 ∶ 1（氯仿 ∶ 異戊醇）進(jìn)行抽提，無(wú)需先用25 ∶ 24 ∶ 1（酚 ∶ 氯仿 ∶ 異戊醇）抽提；（3）加入100 μL的ddH2O和1 μL的RNaseA（10 mg/mL），37 ℃溫浴至DNA完全溶解，不需要進(jìn)行純化處理，直接置于-20 ℃?zhèn)溆?。DNA的提取質(zhì)量采用1%的瓊脂糖凝膠電泳檢測(cè)。

1.2.2 PCR分子檢測(cè) 根據(jù)表達(dá)載體左右邊界（LB和RB）附近區(qū)間的序列采用Primer5設(shè)計(jì)PCR擴(kuò)增引物（表1），其中1條引物設(shè)定在T-DNA邊界區(qū)附近（LB引物設(shè)計(jì)在T-DNA區(qū)內(nèi)，RB引物設(shè)計(jì)在距離T-DNA邊界區(qū)25 bp處），另外1條引物錨定在載體骨架區(qū)（圖1）。對(duì)于所用的PCR分子檢測(cè)引物優(yōu)先在空白對(duì)照、陰性對(duì)照以及陽(yáng)性對(duì)照中進(jìn)行擴(kuò)增，選擇在空白對(duì)照和陰性對(duì)照中無(wú)擴(kuò)增產(chǎn)物，僅在陽(yáng)性對(duì)照中能擴(kuò)增出預(yù)期大小的片段的引物進(jìn)行轉(zhuǎn)基因植株的檢測(cè)。載體骨架序列左右邊界檢測(cè)的引物分別命名為TL-L/R和TR-L/R，TL-L/R和TR-L/R擴(kuò)增的片段大小分別為872 bp和525 bp。引物序列由華大基因合成。

PCR總反應(yīng)體系為15 μL，組成為10×Taq Buffer 1.5 μL，25 mmol/L Mg2+ 1.2 μL；10 mmol/L dNTPs 0.3 μL；5 U/μL Taq酶0.15 μL；10 μmol/L正反向引物各0.3 μL，上述DNA溶液2 μL，滅菌水補(bǔ)足至15 μL。PCR反應(yīng)參數(shù)為：94 ℃預(yù)變性3 min；94 ℃變性30 s，引物TL-L/R 67 ℃復(fù)性40 s（引物TR-L/R 66 ℃復(fù)性40 s），引物TL-L/R 72 ℃延伸1 min（引物TR-L/R 72 ℃延伸40 s），35個(gè)循環(huán)；72 ℃終延伸10 min；4 ℃保存。PCR產(chǎn)物用1%的瓊脂糖凝膠電泳檢測(cè)。PCR檢測(cè)時(shí)，分別以水做空白對(duì)照，未經(jīng)遺傳轉(zhuǎn)化的橡膠樹(shù)葉片基因組DNA為陰性對(duì)照，含有目的基因的表達(dá)載體質(zhì)粒為陽(yáng)性對(duì)照。

1.2.3 PCR產(chǎn)物的回收 TL-L/R和TR-L/R兩對(duì)引物擴(kuò)增出轉(zhuǎn)基因單株的PCR產(chǎn)物用1%的瓊脂糖凝膠電泳檢測(cè)，選擇轉(zhuǎn)基因株系中的一個(gè)單株擴(kuò)增的目標(biāo)條帶進(jìn)行挖膠，用E.Z.N.A.TM Gel Extraction Kit（OMEGA）對(duì)目的產(chǎn)物回收，并將回收產(chǎn)物取2 μL進(jìn)行瓊脂糖凝膠電泳檢測(cè)。

1.2.4 PCR產(chǎn)物的克隆與序列分析 將PCR產(chǎn)物回收純化后與pMDR-19T vector（TaKaRa）連接，方法參照載體試劑盒說(shuō)明書(shū)。轉(zhuǎn)化大腸桿菌DH5α，鑒定并挑選至少兩個(gè)陽(yáng)性克隆送至華大基因測(cè)序，其中TL-L/R鑒定的轉(zhuǎn)基因植株進(jìn)行雙向測(cè)序，而TR-L/R鑒定的轉(zhuǎn)基因植株進(jìn)行單向測(cè)序即可。測(cè)序獲得的序列采用LASERGENE 6.0中的Seqman軟件去除載體序列后進(jìn)行片段拼接。利用MUSCLE 軟件（http：//www.ebi.ac.uk/Tools/msa/muscle/）進(jìn)行序列數(shù)據(jù)的比對(duì)[13]，并用GeneDoc軟件可視化輸出。

2 結(jié)果與分析

2.1 PCR分子檢測(cè)

上述62份材料的葉片DNA擴(kuò)增結(jié)果表明，TL-L/R引物在15份轉(zhuǎn)基因陽(yáng)性植株（代表5個(gè)株系）中擴(kuò)增出872 bp的片段大小，占轉(zhuǎn)基因陽(yáng)性株系的33.3%，而非陽(yáng)性單株和其余22份陽(yáng)性單株中無(wú)擴(kuò)增產(chǎn)物；TR-L/R引物在6份轉(zhuǎn)基因陽(yáng)性植株（代表2個(gè)株系）中擴(kuò)增出525 bp的片段大小，占轉(zhuǎn)基因陽(yáng)性株系的13.3%，而在非陽(yáng)性單株及其余31份陽(yáng)性單株中無(wú)擴(kuò)增產(chǎn)物。PCR分析的結(jié)果顯示，含有左側(cè)邊界序列的陽(yáng)性植株相對(duì)較多，而含有右側(cè)邊界序列的轉(zhuǎn)基因植株同時(shí)也含有左側(cè)的邊界序列，在15個(gè)轉(zhuǎn)基因陽(yáng)性株系中有2個(gè)株系均含有左右邊界的序列。部分樣品檢測(cè)的結(jié)果如圖2所示，全部植株檢測(cè)結(jié)果如表2所示。

2.2 生物信息學(xué)分析結(jié)果

將TL-L/R檢測(cè)到的轉(zhuǎn)基因植株P(guān)T55-1、Tj200863-1、Tp0828-1、Tj2008171-1和T2009401-1進(jìn)行測(cè)序，并與pCAMBIA2301載體序列進(jìn)行比對(duì)發(fā)現(xiàn)序列完全一致（圖3-A），該比對(duì)結(jié)果說(shuō)明這些轉(zhuǎn)基因植株中確實(shí)存在左邊界的序列整合到植物基因組中。同理，將TR-L/R檢測(cè)到的兩株轉(zhuǎn)基因植株P(guān)T55-1和T2009401-1進(jìn)行測(cè)序，與pCAMBIA2301載體序列的右邊界區(qū)域完全一致（圖3-B），說(shuō)明這兩株植株中，右邊界的序列也隨著T-DNA區(qū)的轉(zhuǎn)移而轉(zhuǎn)移至植物基因組中。但在這些植株中，邊界序列插入到基因組中的片段大小目前還在進(jìn)一步的研究中。

3 討論

載體骨架序列整合進(jìn)植物染色體中是轉(zhuǎn)基因植物中的一種普遍現(xiàn)象，并無(wú)作物特異性。在擬南芥[3，5，14]、水稻[15]、玉米[16]、煙草[4]、馬鈴薯[17]、棉花[8]和矮牽牛[18]中均有報(bào)道。研究結(jié)果表明，在農(nóng)桿菌介導(dǎo)的遺傳轉(zhuǎn)化過(guò)程中，載體骨架序列與植物基因組之間的整合頻率一般在20%～50%之間，有時(shí)甚至高達(dá)75%以上[4，19-20]。本研究左右邊界整合的頻率分別為33.3%和13.3%，相比他人的研究報(bào)道來(lái)講，骨架序列存在的比例偏低，可能是由于本研究中利用PCR方法結(jié)合測(cè)序的方法分析的轉(zhuǎn)基因陽(yáng)性株系數(shù)目有限，僅為15個(gè)株系，因此還需要更多的轉(zhuǎn)基因陽(yáng)性株系來(lái)評(píng)估橡膠樹(shù)中載體骨架序列的轉(zhuǎn)移比率。

目前，對(duì)于載體骨架序列整合進(jìn)植物染色體的原因還沒(méi)有得到一致的結(jié)論，T-DNA鏈從LB處起始覆蓋載體骨架亦或是在RB處開(kāi)始形成的T-鏈在載體的LB處終止失敗而最終造成通讀（read-through），這兩方面的推測(cè)均有相關(guān)的報(bào)道[3，6-7，21-22]。

本研究中發(fā)現(xiàn)轉(zhuǎn)基因陽(yáng)性植株中檢測(cè)到左邊界序列的株系明顯高于右邊界序列（分別為33.3%和13.3%），這與De Buck等[6]的研究結(jié)果較為一致，載體骨架DNA多與左側(cè)邊界序列而非右側(cè)邊界串聯(lián)，且有很高比率的轉(zhuǎn)化子含有與左右邊界均相連的載體DNA序列。而本研究還發(fā)現(xiàn)左右邊界序列均有部分整合在植物染色體上的現(xiàn)象，且含有右邊界序列的株系均同時(shí)含有左邊界序列，這與Wenck 等[5]的研究結(jié)果一致，右邊界載體骨架序列的整合必然會(huì)伴隨著左邊界序列的整合，進(jìn)而形成整個(gè)載體的全部整合。但由于本研究檢測(cè)過(guò)程中所使用的引物只是距離兩側(cè)邊界序列較近的引物，因此對(duì)于這些轉(zhuǎn)基因植株中是否含有整個(gè)的載體序列，還將進(jìn)一步設(shè)計(jì)相對(duì)較遠(yuǎn)的引物或者利用Tail-PCR的方法確定這些檢測(cè)到的轉(zhuǎn)基因陽(yáng)性株系插入骨架載體片段的大小。目前已經(jīng)獲得少數(shù)植株側(cè)翼序列的結(jié)果（利用Tail-PCR法獲得的PT55-1、T2009401-1左右序列，與本研究中的PCR結(jié)果和測(cè)序分析的結(jié)果一致，均為載體骨架區(qū)序列；獲得的HbCBF1-3左右序列經(jīng)分析發(fā)現(xiàn)均為非載體骨架序列，將獲得的側(cè)翼序列與橡膠樹(shù)基因組比較發(fā)現(xiàn)為橡膠樹(shù)基因組序列，表明該植株僅為T-DNA區(qū)插入到橡膠樹(shù)染色體上，并根據(jù)生物信息學(xué)分析獲得該插入基因全長(zhǎng)，為功能未知基因），其余植株的Tail-PCR仍在進(jìn)行中。待這些側(cè)翼序列均獲得后，可以對(duì)所有的插入位點(diǎn)側(cè)翼序列進(jìn)行生物信息學(xué)分析，初步尋找橡膠樹(shù)基因組序列中外源基因整合位點(diǎn)的特征，初步分析整合位點(diǎn)的純合性（目前，根據(jù)側(cè)翼序列以及基因特異性引物組合進(jìn)行PCR擴(kuò)增，推測(cè)獲得的轉(zhuǎn)基因植株HbCBF1-3為整合位點(diǎn)雜合子）、T-DNA插入是否存在偏愛(ài)性（預(yù)示基因組該段DNA序列易于外源基因插入）或是隨機(jī)整合性，為外源基因在橡膠樹(shù)轉(zhuǎn)基因植株中整合機(jī)制的研究和轉(zhuǎn)基因效率的提高具有重大的意義。

正因?yàn)檩d體骨架序列存在可能會(huì)導(dǎo)致一系列的問(wèn)題存在[6，9-11]，因此對(duì)其成功的剔除顯得尤為重要。目前己有或者正在嘗試多種方法來(lái)剔除載體骨架序列，如對(duì)表達(dá)載體進(jìn)行改造，即設(shè)計(jì)和使用較小的雙元表達(dá)載體用于遺傳轉(zhuǎn)化[23]或篩選只含有T-DNA區(qū)的轉(zhuǎn)基因植株，即構(gòu)建含非T-DNA區(qū)致死基因的載體，篩選或富集只含有T-DNA序列的轉(zhuǎn)基因后代等[24-25]。不過(guò)大多數(shù)減少載體非必需DNA序列的方法主要用于煙草和擬南芥等模式植物中，真正普遍應(yīng)用于轉(zhuǎn)基因植物還需要很長(zhǎng)一段時(shí)間[26]。對(duì)于橡膠樹(shù)而言，載體骨架序列的發(fā)現(xiàn)，要求我們必須對(duì)其進(jìn)行改造剔除，目前本研究室已經(jīng)開(kāi)展設(shè)計(jì)和使用較小的載體進(jìn)行遺傳轉(zhuǎn)化，還可以嘗試參照模式植物T-DNA區(qū)致死基因的載體研究的思路進(jìn)行優(yōu)化，以期獲得僅含有T-DNA區(qū)的轉(zhuǎn)基因植株。

參考文獻(xiàn)

[1] Zambryski P C. Chronicles from the Agrobacterium-plant cell DNA transfer story[J]. Annu Rev Plant Physiol Plant MOI Biol， 1992， 43： 465-490.

[2] Martineau B， Voelker T A， Sanders R A. On defining T-DNA[J]. The Plant Cell， 1994， 6（8）： 1 032-1 033.

[3] van der Graaff E， den Dulk-Ras A， Hooykaas P J J. Deviating T-DNA transfer from Agrobacterium tumefaciens to plants[J]. Plant Mol Biol， 1996， 31： 677-681.

[4] Kononov M E， Bassuner B， Gelvin S B. Integration of T-DNA binary vector‘backbone sequences into the tobacco genome： evidence for multiple complex patterns of integration[J]. The Plant Journal， 1997， 11（5）： 945-957.

[5] Wenck A， Czakó M， Kanevski I， et al. Frequent collinear long transfer of DNA inclusive of the whole binary vector during Agrobacterium-mediated transformation[J]. Plant molecular biology， 1997， 34（6）： 913-922.

[6] De Buck S， Wilde C D， Montagu M V， et al. T-DNA vector backbone sequences are frequently integrated into the genome of transgenic plants obtained by Agrobacterium-mediated transformation[J]. Molecular Breeding， 2000， 6： 459-468.

[7] Abdal-Aziz S A， Pliego-Alfaro F， Quesada M A， et al. Evidence of frequent integration of non-T-DNA vector backbone sequences in transgenic strawberry plant[J]. Journal of bioscience and bioengineering， 2006， 101（6）： 508-510.

[8] 馬小波， 王芙蓉， 張傳云， 等. 多位點(diǎn)整合的轉(zhuǎn)基因棉花后代分子鑒定[J]. 棉花學(xué)報(bào)， 2008， 20（1）： 9-13.

[9] Iglesias V A， Moscone E A， Papp I， et al. Molecular and cytogenetic analyses of stably and unstably expressed transgene loci in tobacco[J]. Plant Cell， 1997， 9： 1 251-1 264.

[10] Stoger E， Williams S， Keen D， et al. Molecular characteristics of transgenic wheat and the effect on transgene expression[J]. Transgenic Res， 1998， 7： 463-471.

[11] 王利華， 蘇 喬， 包永明. 轉(zhuǎn)基因植物中載體骨架序列的安全性隱患及解決方案[J]. 中國(guó)生物工程雜志， 2004， 24（5）： 38-42.

[12] 李 季， 黃天帶， 蔡海濱， 等. 橡膠樹(shù)轉(zhuǎn)基因植株遺傳穩(wěn)定性分析[J]. 熱帶作物學(xué)報(bào)， 2013， （4）： 591-595.

[13] Edgar R C. MUSCLE： multiple sequence alignment with high accuracy and high throughput[J]. Nucleic Acids Res， 2004， 32： 1 792-1 797.

[14] De Greve H， Nguyen V K， Deboeck F， et al. T-DNA tagging of the translation initiation factor eIF-4Al of Arabidopsis thaliana[J]. Plant Sci， 2001， 161（4）： 685-693.

[15] Zhu Q H， Ramm K， Eamens A L， et al. Transgene structures suggest that multiple mecha-nisms are involved in T-DNA integration in plants[J]. Plant Sci， 2006， 171（3）： 308-322.

[16] Shou H， Frame B R， Whitham S A， et al. Assessment of transgenic maize events produced by particle bombardment or Agrobacterium-mediated transformation[J]. Mol Breed， 2004， 13（2）： 201-208.

[17] Wolters Am-A， Trindade L M， Jacobsen E， et al. Fluorescence in situ hybridization on extended DNA fibres as a tool to analyse complex T-DNA loci in potato[J]. Plant J， 1998， 13： 837-847.

[18] Cluster P D， O'Dell M， Metzlaff M， et al. Details of T-DNA structural organization from a transgenic Petunia population exhibiting co-suppression[J]. Plant Mol Biol， 1996， 32： 1 197-1 203.

[19] Kuraya Y， Ohta S， Fukuda M， et al. Suppression of transfer of non-T-DNA ‘vector backbone sequences by multiple left border repeats in vectors for transformation of higher plants mediated by Agrobacterium tumefaciens[J]. Molecular Breeding， 2004， 14： 309-320.

[20] 成善漢， 謝從華， 柳 俊， 等. 馬鈴薯轉(zhuǎn)基因株系載體骨架的整合分析[J]. 海南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2006， 24（3）： 305-308.

[21] Sallaud C， Meynard D， van Boxtel J， et al. Highly efficient production and characterization of T-DNA plants for rice （Oryza sativa L.）functional genomics[J]. Theor Appl Genet， 2003， 106： 1 396-1 408.

[22] van Haaren M J J， Pronk J T， Schilperoort R A， et al. Functional analysis of the Agrobacterium tumefaciens octopine Ti-plasmid left and right T-region border fragments[J]. Plant Mol Biol， 1987， 8： 95-104.

[23] Xiang C B， Han P， Lutziger， et al. A mini binary vector series for plant transformation[J]. Plant Molecular Biology， 1999， 40（4）： 711-717.

[24] Hanson B， Engler D， Moy Y， et al. A simple method to enrich an Agrobacterium-transformed population for plants containing only T-DNA sequences[J]. The Plant Journal， 1999， 19（6）： 727-734.

[25] 董志峰， 馬榮才， 彭于發(fā)， 等. 轉(zhuǎn)基因植物中外源非目的基因片段的生物安全研究進(jìn)展[J]. 植物學(xué)報(bào)， 2001， 43（7）： 661-672.

[26] 馬小波. 農(nóng)桿菌介導(dǎo)獲得的轉(zhuǎn)基因棉花基因組中外源基因結(jié)構(gòu)與整合方式[M]. 濟(jì)南： 山東師范大學(xué)圖書(shū)館， 2007.