外源基因插入位點對抗蟲棉Bt蛋白含量的影響

中圖分類號：S562；Q786 文獻標志碼：A 文章編號：1008-0864（2025）07-0044-10

Effect of Exogenous Gene Insertion Site on Bt Protein Content in Insect-resistantCotton

ZHAO Guiyuan ^1，2 ，WANG Yongqiang2，LIU Jianguang2，GENG Zhao2， ZHANG Hanshuang2 WULiqiang1，WANG Xingfen1*，ZHANG Guiyin1*

（1.StateKeyLaboratoryofNorthChinaCropImprovemenandRegulation，KeyLaboratoryforCropGermplasmResourcesof Hebei，ColegeofAgronomy，Hebei AgriculturalUniversity，HebeiBaodingO7Oo，China;2.KeyLaboratoryofCottonBiology andGeneticBreedinginHuanghuaihaiSemiaridAreaofMinistryofAgricultureandRuralAfairs，InstituteofCoton，Hebei AcademyofAgricultureand forestrySciences，ShijiazhuangO5oO51，China）

Abstract：The insertion sites of exogenous genes can affct the gene expression and character expression of recipient materials.Inorder todeterminetheefectof exogenousgeneinsertionsitesonBtproteincontentof insect-resistant coton，29 insect-resistant coton lines were selected for Bt protein content determination and wholegenome resequencing analysis.The results showed thatthere were significant diferences in Bt protein content among diffrent lines，and the K3 line was the highest.Whole-genome resequencing analysis and PCR identification revealedthat the Btgene wasa single copyin 29 lines，but the insertion sites were different. The Bt gene of K3 line was located in the 2 791 303＼～2 791 335 bp of A02 chromosome，while the Bt gene of other lines were located in D12 chromosome. To further verifytheefectof insertionsiteon Btproteincontent，theK3 line was hybridization withcontrol Xinmian 33B，and F₂ population was constructed. The results showed that the Bt protein content of F₂ individual plants with theinsertion siteof Bt gene in AO2 chromosome was significantly higher than that D12 type individual plants.In summary，the Bt protein content was affected by the insertion site of Bt gene，and its location on chromosome A02 was more conducive to the efficient expression of Bt gene.Above results could provide germplasm resources and theoretical support for breeding new varieties of insect-resistant cotton.

Keywords：insect-resistant cotton；exogenous gene；insertion site；Bt protein

棉花（Gossypium spp.）是集纖維、油料、飼料于一身的重要經濟作物。棉鈴蟲（Helicoverpaarmigera）是棉花生長過程中面臨的主要害蟲之一，每年給我國棉花產業造成的損失超過總產量的 10%^[2-3] 。近年來，隨著 Bt 基因成功植入棉花，棉鈴蟲危害得到了顯著控制[4-5]。棉花已成為全球第三大轉基因作物，2022年轉基因抗蟲棉的種植面積達到2540萬 hm² ，占棉花總種植面積的 80% 以上。

盡管轉基因抗蟲棉在棉花生產中發揮了重要作用，但不同品種間Bt蛋白含量的差異導致了抗蟲效果的不一致性，進而影響了產量的進一步提高。目前，國內外學者對這一問題進行了大量研究，但多集中在外界環境因素的影響方面，如溫度、王壤鹽度和肥料等。周冬生等和陳源等研究顯示，溫度是影響轉 Bt 棉抗蟲性的重要環境因素，尤其是低溫對抗蟲棉葉片Bt蛋白含量的影響最為顯著。同時，土壤鹽度較高也會顯著影響抗蟲棉在不同生長階段的Bt蛋白含量]。氮素代謝與抗蟲棉Bt蛋白含量之間呈顯著正相關，直接影響抗蟲性的高低]。此外，空氣中的二氧化碳和臭氧含量等環境因素也會影響轉 Bt 基因棉的抗蟲性[12-13]

轉基因作物的分子特征包括外源基因的表達載體、T-DNA片段特征、拷貝數、插入位點和穩定性等遺傳信息[4]，其中拷貝數和插入位點尤為關鍵[15]。郭旺珍等研究發現，抗蟲棉品系中 Bt 基因拷貝數的增多并沒有使其抗性增強，反而使外源基因的表達受到一定抑制。但也有研究發現，擁有多拷貝外源基因的轉基因作物能夠持續的表達和遺傳[17]。夏蘭芹等[18]利用轉基因方法以多拷貝方式將 Bt 基因整合并轉入相關植物體中，并且檢測到穩定的蛋白含量，從而說明外源基因能夠在受體植株中持續穩定的表達。外源基因隨機插入宿主基因組，其整合位置可能是染色體的特定位點或染色體上的任意位置[9]。外源基因的表達水平受其插入位點的影響，從而在植物體上產生特定的外在表現。研究表明，轉基因植株間外源基因表達的顯著差異往往由其在宿主基因組中的不同插入位點引起[20-21]。Pedersen等[22]探討了不同麥類轉化體的插入整合規律，并強調了外源基因插入位點對其穩定表達的重要性。李亞麗等[23]對轉基因水稻的插入位點分析發現，T-DNA表現出一定的位置效應。

盡管關于外源基因插入位點的研究較多，但關于不同插入位點對受體材料影響的研究卻相對較少。本研究團隊前期培育出一個Bt蛋白含量極高的轉基因抗蟲棉純合品系K3，其平均蛋白含量是傳統轉基因抗蟲棉新棉33B的1.8倍以上。因此，本研究對不同轉基因抗蟲棉品系的Bt蛋白含量進行比較，并通過全基因組重測序技術分析了不同抗蟲棉品系的 Bt 基因插入位點；利用Bt蛋白含量超高但 Bt 基因插入位點不同的K3品系與對照新棉33B構建雜交群體，鑒定分析F群體中不同類型單株的Bt蛋白含量差異規律，從而進一步驗證插入位點對抗蟲棉Bt蛋白含量的影響，旨在從外源基因分子水平揭示抗蟲棉Bt蛋白差異機制，研究結果將為轉基因棉高抗蟲育種研究提供重要參考。

1材料與方法

1.1供試材料

試驗材料分為品系材料和群體材料。品系材料包括29個轉 Bt 基因抗蟲棉品系，均為高代純合品系，由來源不同的優勢非轉 Bt 基因ND系列材料與新棉33B（33B）雜交獲得，以新棉33B為陽性對照，以非轉基因抗蟲棉中棉所12（CCRI12）為陰性對照（表1）。群體材料為品系K3與新棉33B雜交后的 F₁ 和 F₂ 群體。品系材料用于分析不同品系間Bt蛋白含量和 Bt 基因插入位點差異；群體材料用于鑒定不同插入位點群體Bt蛋白含量差異。所涉及的 Bt 基因均為 CryIAc 基因。所有供試材料均由河北農業大學棉花遺傳育種研究室選育與創制。

試驗于2019—2023年在河北省農林科學院棉花研究所小安舍試驗站（石家莊市新華區，38^°05^′N，114^°25^′E） 進行。試驗地為中壤土，土壤肥力中等，有嚴格的轉基因隔離防護帶。于2019—2022年種植29份品系材料，采用隨機區組排列，每個品系種植4行區，3次重復，兩邊設保護行。于2021年配置品系K3與新棉33B雜交組合，2022年種植 F_r 群體，2023年種植 F₂ 群體及父母本。所有材料均為4月下旬田間播種，等行距種植，行距 76cm ，行長 7m ，株距 20cm 。采用常規棉花田間栽培水肥管理，整個生育期不施用防治棉鈴蟲的藥劑。

1.2試驗方法

1.2.1Bt蛋白含量檢測試驗材料種植期間于蕾期取頂部第1片完全展開葉片測定Bt蛋白含量。樣品用液氮速凍后放入干冰中，迅速帶回實驗室檢測或放入 -80^°C 超低溫冰箱中保存并盡快檢測。

采用酶聯免疫吸附測定試劑盒（enzymes linked immuno sorbent assays，ELISA）測定Bt蛋白 含量[]

1.2.2全基因組重測序采用Paterson等24方法提取各品系葉片DNA，采用瓊脂糖凝膠電泳和Nanodrop檢測DNA質量；將質量檢測合格的DNA樣品用于建庫，送至安諾優達基因科技（北京）有限公司進行測序。

建庫流程如下：檢驗合格的DNA樣品通過Covaris破碎機隨機打斷成長度為 350bp 的片段[25]。采用 TruSeq Library Construction Kit進行建庫，嚴格使用說明書推薦的試劑和耗材。DNA片段經末端修復、加ployA尾、加測序接頭、純化、PCR擴增等步驟完成整個文庫制備。文庫構建完成后進行illumina NovaSeq 600O 測序[26]。

1.2.3插入位點分析與鑒定方法將所測序列與TM-1基因組序列進行合并得到參考文件ref.fa。將過濾后cleanreads通過BWA比對到ref.fa基因組得到bam文件[27]。用SAMtools對比對結果進行排序；再用Picard標記重復序列（markduplicatereads）[28]。從bam文件中篩選出比對到的外源序列reads。將篩選出的reads與TM-1基因組進行比對，預測插入位點。

根據預測出的插入位點側翼序列和 Bt 基因序列，利用PrimerPremier5.0軟件設計特異擴增側翼序列的PCR引物（表2），并由生工生物工程（上海）股份有限公司合成。引物設計的原則是跨拼接區域：左側鑒定引物的上游位于參考基因組上，下游位于載體序列上（LeftBorder，即LB引物）；右側鑒定引物的上游位于載體序列，下游位于參考基因組（RightBorder，RB引物）29。通過瓊脂糖凝膠電泳，能擴增到目的長度片段的即可驗證插入位點。

PCR體系： 2× SuperPfxMasterMix 12.5μL ，正引物 1μL ，反引物 1μL ，DNA 1μL μL，ddH₂O9.5μL_c PCR程序： 98^°C 預變性 5min;98^°C 變性 30s，55～ 65^°C 退火 30s，72° 延伸30s，34個循環； 72^°C 終延伸 10min 。PCR產物用 1% 瓊脂糖凝膠電泳檢測。

1.2.4群體純合型單株鑒定于2023年5月苗期取 F₂ 群體單株上部 1cm² 鮮嫩葉片置于 2mL 離心管，冰上保存帶回實驗室。采用Paterson等24方法提取DNA，根據表2引物序列對其進行瓊脂糖凝膠電泳鑒定。對 F₂ 單株鑒定過程中，擴增到A02型條帶且未擴增到D12型條帶的稱為A02型單株；擴增到D12型條帶且未擴增到A02型條帶的稱為D12型單株。盡快鑒定完并標記純合型單株，用于蕾期不同類型單株Bt蛋白含量的檢測。

1.3 數據分析

使用Excel軟件進行試驗數據的整理與分析，采用Origin軟件制圖。

2 結果與分析

2.1不同抗蟲棉品系Bt蛋白含量差異

對29份抗蟲棉品系的Bt蛋白含量比較（圖1）表明，不同品系的平均Bt蛋白含量存在差異。盡管同一品系在不同年份會因氣候條件、栽培措施等因素出現波動，但整體趨勢保持一致，表明Bt蛋白含量主要受品系自身的遺傳背景和內在分子機制影響。

進一步分析表明，Bt蛋白含量在 以上的只有K3品系，極顯著高于其他品系；含量在 1000～2000ng^?g^-1 有2個，為K12品系和陽性對照新棉33B（K26）；其余品系的Bt蛋白含量低于 1 000ng?g^-1 。29個品系中，K3品系的Bt蛋白含量極顯著高于對照新棉33B，K12品系與新棉33B差異不顯著，有7個品系顯著低于新棉33B，有19個品系極顯著低于新棉 33B 。綜上表明，各品系間Bt蛋白含量存在顯著差異。

2.2 插入位點分析和PCR鑒定

對29個品系進行全基因組重測序，共獲得總堿基數超 700Gb 的數據，測序深度超過 10×，Q30 值均超過 92%，Q20 值均為 100% 。其中品系K3的原始reads數為 145 574 824 條，原始Bases為23.7G，Q30為 95.69% ，其中測序深度 10× 的比例為76.75% 。新棉33B的原始reads數為158525832條，原始Bases為23.6G，Q30為 93.75% ，其中測序深度 10× 的比例為 76.69% 。

對全基因組重測序數據的進一步分析表明，29個品系中 Bt 基因的拷貝數均為1，但插入位點有2種，分別位于A02染色體和D12染色體（表3）。K3品系中 Bt 基因插入到A02染色體2791303＼～2791335處，插入位點處有 32bp 的堿基缺失，將這種類型稱為A02型；其余品系均插入到D12染色體，除K13和K23外的26份品系均插入49127447＼～49127527bp處，且插入位點處有80bp 的堿基缺失，將這種類型稱為D12型。

為進一步明確各品系 Bt 基因插入位點，分別利用D12和A02染色體特異性引物對其進行PCR鑒定。D12染色體特異性引物PCR鑒定表明，除K3外的其他品系都擴增出相應的目的條帶（圖2）；A02染色體特異性引物PCR鑒定表明，只有K3品系擴增出目的條帶，其他品系均未擴增出條帶（圖3）。綜合重測序結果表明， Bt 基因插入到K3品系A02染色體相應位置，而其余品系則插入到D12染色體相應位置。

2.3群體純合型單株鑒定與Bt蛋白含量比較

由圖4可知，K3品系的Bt蛋白含量極顯著高于新棉 33B 。2019—2023年新棉33B的Bt蛋白含量為 1151.00～1543.70ng?g^-1， ，而K3的蛋白含量為 2140.93～2765.84ng^?g^-1 ，較新棉33B顯著提高65.17%～96.76% ，二者差異達極顯著水平。

因K3品系的Bt蛋白含量極顯著高于對照新棉33B，且二者Bt基因插入位點亦不相同，故利用K3與新棉33B雜交后的 F₁ 和 F₂ 群體單株Bt蛋白含量分布規律進一步驗證插入位點對Bt蛋白含量的影響。2022年獲得K3與新棉33B的 F₁ 單株369株，2023年獲得 F₂ 單株921株。通過特異性引物對其進行基因型鑒定表明，A02型和D12型單株的數量分別為42和53株。

2022年部分 F₁ 代單株Bt蛋白含量測定（表4）表明，其含量為 1567.73～2762.57ng^?g^-1. ，不同單株間差異較大，但介于2個親本之間。通過對A02型和D12型2種純合型 F₂ 單株Bt蛋白含量研究（圖5）表明，A02純合型單株的蛋白含量為2311.40～2771.43ng^?g^-1， ，均極顯著高于D12型。由此表明， Bt 基因插入位點位于A02染色體的 F₂ 單株的蛋白含量顯著高于D12型單株。綜合分析表明，Bt蛋白含量受 Bt 基因插入位點影響，當其插入A02染色體時更利于 Bt 基因的高效表達。

3討論

隨著對轉基因作物的不斷研究，外源基因的不穩定性和表達多樣性等因素增加了基因高效表達的難度，從而制約了轉基因育種的研究與利用。因此，尋找和創制高表達且穩定的轉基因材料就顯得尤為迫切。本研究發現，K3品系的 Bt 基因插入到A02染色體后極大地提高了抗蟲棉的Bt蛋白含量，是其他品系的1.8＼～5.3倍，且通過K3品系構建的后代群體實現了超高Bt蛋白的遺傳和表達。研究結果為培育高效抗蟲棉新品種提供了優異種質資源。

轉基因技術作為農業領域中發展較快、應用較廣泛的新技術之一，已成為打破物種間隔離障礙的有效手段[30]。外源基因在受體植物基因組中的插入具有隨機性，因此，明確外源基因的插入位點對于轉基因作物的環境釋放和安全評估至關重要3]。隨著下一代測序技術（next-generationsequencing，NGS）的成熟發展和廣泛商業化運用，利用全基因組測序技術分析轉基因作物的插入位點已成為可能，這種方法不僅精準而且便捷[32]。研究證實，傳統的Southern雜交與NGS檢測結果之間具有較高的一致性[33]。利用重測序分析鑒定插入位點的方法已在多種轉基因動植物基因組中成功應用，實現了對外源基因的精確定位[26.3436]

本研究表明，不同抗蟲棉品系的Bt蛋白含量顯示出不同程度差異。其中，K3品系由于Bt基因插人位點的不同，Bt蛋白含量顯著提高。而在插入位點相同的情況下，關于外源基因表達差異的研究多集中在甲基化方面。甲基化包括DNA甲基化和RNA甲基化，是植物體內重要的表觀遺傳修飾方式，能夠顯著調控基因的表達和系統發育[37-40]。DNA甲基化對基因的調控多發生在分子水平上，在一定程度上影響植物的生長和發育[4u]。DNA甲基化常常與啟動子區的甲基化聯系在一起[42]。研究發現，抗蟲棉在生殖生長階段發生的甲基化能夠直接影響Bt蛋白的表達，從而降低抗蟲棉的抗蟲性；并且Bt蛋白含量的高低跟35S啟動子甲基化直接相關43]。涉及啟動子的DNA甲基化可發生在植物的多種代謝途徑中，對基因的不表達或高效表達起關鍵作用44。作為轉錄后水平的一種調控方式，RNA甲基化對真核生物基因的表達同樣起著重要的調節作用。RNA轉錄后修飾已成為表觀遺傳學領域的研究熱點[45]。 N⁶. 甲基腺苷 （ N⁶ -methyladenosine， m⁶A 是RNA腺嘌呤核苷酸在第6位氮原子上發生的甲基化修飾現象。這種修飾在mRNA中極為常見，并具有重要的轉錄后調控功能。從哺乳動物到病毒再到植物， m⁶A 修飾都發揮著不可忽視的作用[46-47]。對于具有相同插入位點的品系間Bt蛋白含量差異的分子機制，未來可進一步利用DNA和RNA甲基化相關分子技術進行探究。

將外源基因成功插入受體材料并且實現高效表達，意味著該區域是理想的轉錄激活區域，可以作為候選的轉基因友好位點。因此，未來通過多種方式的外源基因轉化事件，結合基因編輯技術，能夠將感興趣的特定基因精確地插入這些友好位點，從而獲得穩定遺傳和表達的轉基因材料，這將避免隨機整合所引發的諸多問題[35]。

參考文獻

[1］夏薈菁.轉 CryIC^* 基因棉花材料的抗蟲性鑒定與評價[D]. 武漢：華中農業大學，2023. XIA _HJ ，Insect-resistant identification and evaluation of transgenic CryIC^* cotton [D].Wuhan：Huazhong Agricultural University，2023.

[2]文學，張寶紅.轉基因抗蟲棉研究現狀與展望[J].農業生物 技術學報，2000，8（2）：194-199. WENX，ZHANG BH.Progressand prospects for insect-resistant transgeniccotton[J].J.Agric.Biotechnol.，20OO，8（2）：194-199.

[3]張京飛.轉Bt基因抗蟲棉抗蟲性研究與新品系篩選[D].保 定：河北農業大學，2014. ZHANG JF.Identification of insect resisitance and new variety screening in Bt transgenic cotton [D].Baoding：Hebei Agricultural University，2014.

[4]鄒奎.我國轉基因抗蟲棉品種現狀分析[J].中國棉花，2003， 30（8）：2-4.

[5] 徐澤俊，聶以春，張獻龍，等.轉雙價抗蟲基因棉花的主要農 藝性狀的遺傳變異[J].植物遺傳資源學報，2011，12（1）： 125-130. XUZJ，NIEYC，ZHANGXL，etal..Genetic variation of mainagronomic traits of transgenic insect resistant linesin cotton[J].J.PlantGenet.Resour.，2011，12（1）：125-130.

[6]李浩輝，劉彩月，張海文，等.2022年度全球轉基因作物產業 化發展現狀及趨勢分析[JJ.中國農業科技導報，2023，25（12）： 6-16. LI HH，LIU CY，ZHANGHW，et al.. Global genetically modified crop industrialization trendsin2022[J].J.Agric.Sci. Technol.，2023，25（12）：6-16.

[7]谷淇深，李志坤，王省芬，等.Bt抗蟲棉新品系毒蛋白表達差 異分析[J].棉花學報，2016，29（2）：128-137. GUQS，LI ZK，WANGXF，et al.Differential expression analysis of Bt protein in new Bt cotton lines [J].Cott.Sci.， 2016，29（2）：128-137.

[8]周冬生，吳振延，王學林，等.轉 Bt 基因棉的抗棉鈴蟲性及其 生理作用研究進展[J].安徽農業科學，2000，28（1）：65-68. ZHOU D S，WU Z Y， WANG X L，et al. Progress and prospectsof the resistance to cotton boll wormand physiological action of Bt transgenic cotton [J]. J. Anhui Agric. Sci.，2000，28（1）：65-68.

[9]陳源，顧超，王桂霞，等.蕾期低溫及濕度脅迫對Bt棉殺蟲蛋 白表達量的影響[J].作物學報，2013，39（1）：184-189. CHENY，GU C，WANG G X， et al..Effect on stresses of 18°℃ and different relative humidities on Bt protein ex-pression at squaring stage inBt cotton [J].Acta Agron.Sin.，2013，39（1）： 184-189.

[10]JIANG L J， DUAN L S，TIAN X L，et al. NaCl salinity stress decreased Bacilus thuringiensis（Bt） protein content of transgenic Bt cotton （Gossypium hirsutum L.） seedlings [J]. Environ.Exp.Bot.，2006，55（3）：315-320.

[11］趙紅霞，王士杰，朱繼杰，等.不同遺傳背景轉基因抗蟲棉Bt 蛋白表達與氮代謝關系研究[J].棉花學報，2018，30（6）： 498-504. ZHAO HX， WANG S J，ZHUJJ， et al.Relationship between Bt protein expression and nitrogen metabolism in insectresistant transgenic cotton lineswith different genetic backgrounds [J].Cott.Sci.，2018，30（6）：498-504.

[12]HIMANEN SJ，NERGA M，NISSINEN A，et al..Elevated atmospheric ozone increases concentration of insecticidal Bacillus thuringiensis （Bt）CryAcprotein inBtBrasicanapus and reduces feeding ofa Bt target herbivore on the nontransgenic parent [J].Environ.Pollut.，20o9，157（1）：181-185.

[13］王家寶，王留明，沈法富，等.環境因素對轉Bt基因棉Bt殺蟲 蛋白表達量的影響[J].山東農業科學，2000，32（6）：4-6. WANGJB，WANGL M，SHENFF，et al.Effect of environment elements on Bt-protein content in transgenic Bt cotton [J]. Shandong Agric.Sci.，20oo，32（6）：4-6.

[14]LIR，QUAN S，YAN XF，et al..Molecular characterizationof genetically-modified crops： challengesand strategies[J]. Biotechnol. Adv.，2017，35（2）：302-309.

[15］王翠云，劉艷，劉允軍.外源基因在轉基因玉米中的整合位點 分析[J].生物技術通報，2019，35（3）：1-5. WANG C Y，LIUY，LIUYJ.Analysis of the integration site of exogenous gene in transgenic maize [J]. Biotechnol.Bull.， 2019，35（3）：1-5.

[16]郭旺珍，孫敬，郭玉芳，等.轉基因抗蟲棉Bt基因不同劑量的 聚合與抗蟲性表現[J].遺傳學報，2001，28（7）：668-676. GUO W Z，SUNJ， GUO YF，etal.. Investigation of different dosages of inserted Bt genes and their insect-resistance in 668-676.

[17］趙梅香.轉基因大麥的 gfp 基因遺傳表達及整合位點結 構[D].揚州：揚州大學，2007. ZHAO M X. Integration loci structure and genetic expression of green fluorescent proteingfp gene in transgenic barley （Horsdeum vulgareL.）[D].Yangzhou：Yangzhou University，2007.

[18］夏蘭芹，王遠，郭三堆.外源基因在轉基因植物中的表達與穩 定性[J].生物技術通報，2000（3）：8-12. XIALQ，WANG Y，GUO S D. The stability of the expression of foreign genes in transgenic plants [J].Biotechnol. Bull.，2000 （3）：8-12.

[19］董美.NGS在生物技術作物分子特征解析中的應用[D].北 京：中國農業科學院，2018. DONG M.Application of NGS for molecular characterization of biotech crops [D].Beijing：Chinese Academy of Agricultural Sciences，2018.

[20]陳天子，凌溪鐵，楊郁文，等.轉GbVe1基因在棉花基因組中的 整合與定位分析[J].棉花學報，2019，31（1）：1-11. CHENTZ，LINGXT，YANGYW，etal.. The integration and insertion siteof GbVel gene in thegenome of transgenic coton （Gossypium hirsutum） [J].Cott.Sci.，2019，31（1）：1-11.

[21］金永梅，馬瑞，于志晶，等.轉 cryIC 基因抗蟲水稻吉生粳3號 外源基因整合分析與品系特異性檢測[J].生物技術通報， 2019，35（3）：6-12. JIN Y M， MA R， YU Z J， et al. Integration analysis of exogenous gene and line-specific detection in the insectresistant cryIC -transgenic riceJishengjing3[J].Biotechnol. Bull.，2019，35（3）：6-12.

[22]PEDERSEN C， ZIMNY J，BECKER D，et al..Localization of introduced genes on the chromosomes of transgenic barley， wheat and Triticale by fluorescence insitu hybridization [J]. Theor.Appl. Genet.，1997，94（6）：749-757.

[23］李亞麗，劉中來，周潔，等.轉mCherry基因水稻的遺傳分析及 T-DNA整合位點的研究[J].分子植物育種，2012，10（2）： 121-130. LI YL，LIU ZL， ZHOU J， et al.. Genetic analysis of mCherry transgenic rice and molecular characterization of T-DNA integration sites in the rice genome[J]. Mol. Plant Breed.， 2012， 10（2）：121-130.

[24]PATERSON A H， BRUBAKER C L， WENDEL JF.A rapid method for extraction of cotton （Gossypium spp.） genomic DNA suitable for RFLPor PCR analysis[J].Plant Mol. Biol. Rep.， 1993，11（2）：122-127.

[25]MA Z Y， HE S P， WANG X F， et al. Resequencing a core collection of upland cotton identifies genomic variation and loci influencing fiber quality and yield [J].Nat.Genet.，2018， 50（6）： 803-813.

[26]KERSTEN B， LEITE MONTALVAO A P， HOENICKA H， t al. Sequencing of two transgenic early-flowering poplar lines confirmed vector-free single-locus T-DNA integration [J]. Transgenic Res.，2020，29（3）：321-337.

[27］魏嘉，王秀東，卜顯峰，等.基于重測序鑒定耐鹽轉基因大豆 事件外源T-DNA整合位點及特異性檢測[J].南京農業大學 學報.2023.46（4）：710-717. WEIJ，WANGXD，BUXF，etal..IdentificationofT-DNA insertion and specific detection of transgenic soybean event based on re-sequencing technology [J].J. Nanjing Agric. Univ.， 2023，46（4）：710-717.

[28］張偉.玉米超表達突變體生成系統的完善與突變體篩選鑒 定[D].武漢：華中農業大學，2023. ZHANG W. Improvement of maize overexpression mutant generation system and identification of mutants [D]. Wuhan： Huazhong Agricultural University，2023.

[29］周星魯.轉雙Bt基因107楊T-DNA整合及對基因表達的影 響[D].保定：河北農業大學，2021. ZHOU X L.T-DNA integration information of double Bt gene populusXeuramericana“Neva” and its response to Bt gene [D]. Baoding：Hebei Agricultural University，2021.

[30]程俊凌，趙亮，徐劍文，等.基于重測序鑒定GbTMEM214基因 在陸地棉基因組的插入位點[J].棉花學報，2024，36（4）： 285-295. CHENGJL，ZHAO L， XUJW. Identification of GbTMEM214 insertion site in Gossypium hirsutum genome based on resequencing[J]. Cott. Sci.，2024，36（4）：285-295.

[31］張美冬，孫玲，熊秋芳.轉基因作物的安全性及其評價[J].湖 北農業科學，2015，54（5）：1025-1030. ZHANGMD，SUNL，XIONGQF.Safetyand evaluation of transgenic crops [J].Hubei Agric.Sci.，2015，54（5）：1025-1030.

[32]CADE R，BURGIN K，SCHILLING K，et al. Evaluation of wholegenomesequencingandaninsertionsite characterization method for molecular characterization of GM maize[J]. J.Regul.Sci.，2018，6（1）：1-14.

[33]PARKD，PARK SH，BANYW，et al..A bioinformatics approach for identifying transgene insertion sites using whole genome sequencing data [J/OL]. BMC Biotechnol.，2017，17（1）： 67[2024-09-20]. https：//doi.org/10.1186/s12896-017-0386-x.

[34]SCHOUTEN HJ，VANDE GEEST H，PAPADIMITRIOU S， et al..Re-sequencing transgenic plants revealed rearrangements at T-DNA inserts，and integration of a short T-DNA fragment， but no increase of small mutations elsewhere [J].Plant Cell Rep.，2017，36（3）：493-504.

[35］魏強，奧巖，楊漫漫，等.利用全基因組重測序技術鑒定五指 山豬GHR突變體轉基因插入位點[J].遺傳，2021，43（12）： 1149-1158. WEI Q，AOY，YANG M M，et al. Identification of genomic insertion of dominant-negative GHR mutation transgenes in Wuzhishan pig using whole genome sequencing method [J]. Hereditas，2021，43（12）：1149-1158.

[36]徐紀明，胡晗，毛文軒，等.利用重測序技術獲取轉基因植物 T-DNA插入位點[J].遺傳，2018，40（8）：676-682. XUJM，HU H，MAO W X，et al.. Identifying T-DNA insertion site（s）of transgenic plants bywhole-genome resequencing[J]. Hereditas，2018，40（8）：676-682.

[37]TIRNAZ S，BATLEY J.DNA methylation：toward crop disease resistance improvement [J]. Trends Plant Sci.，2019，24（12）： 1137-1150.

[38]GALLEGO-BARTOLOME J.DNA methylation in plants： mechanisms and tools fortargeted manipulation[J].New Phytol.，2020，227（1）：38-44.

[39]DELERIS A，HALTER T，NAVARRO L.DNA methylation anddemethylation in plantimmunity [J].Annu.Rev. Phytopathol.，2016，54：579-603.

[40]PERTEA M，KIM D，PERTEAG M，etal..Transcript-level expressionanalysisof RNA-seq experiments with HISAT， StringTie and Ballgown [J]. Nat.Protoc.，2016，11（9）：1650-1667.

[41］孟彩鳳，周穎.DNA甲基化對植物生長發育的調控研究[J] 安徽農學通報，2015，21（10）：29-32. MENGCF，ZHOUY.Regulation effectofDNA methylation on the growth and development of plant [J]. Anhui Agric. Sci. Bull.，2015，21（10）：29-32.

[42]ZHANGXY，YAZAKI J，SUNDARESANA，etal..Genome wide high-resolution mapping and functional analysis of DNA methylation inArabidopsis [J].Cell，2006，126（6）：1189-1201.

[43］夏蘭芹，徐瓊芳，郭三堆.抗蟲棉生長發育過程中Bt殺蟲基 因及其表達的變化[J].作物學報，2005，31（2）：197-202. XIALQ，XUQF，GUOSD.Btinsecticidal gene and its temporal expression in transgenic cotton plants [J].Acta Agron.Sin.，2005，31（2）：197-202.

[44]CHEN ZY，RIU E，HE CY，et al.. Silencing of episomal transgene expressionin liverbyplasmid bacterialbackbone DNA is independent of CpG methylation [J]. Mol. Ther.， 2008， 16（3）：548-556.

[45］朱婷.m6A甲基化修飾調控YS型小麥溫敏雄性不育系花藥 育性的分子機制[D].楊凌：西北農林科技大學，2023. ZHUT. The molecularregulation mechanismof m⁶A modification on anther fertility of YS-type thermo-sensitive male sterile line in wheat （Triticum aestivumL.）[D]. Yangling： NorthwestAamp;FUniversity，2023.

[46]AKHTAR J，LUGOBONI M，JUNION G. m＼～6A RNA modification in transcription regulation [J]. Transcription， 2021，12（5）：266-276.

[47]SONG P Z， TAYIER S， CAI ZH， et al..RNA methylation in mammalian development and cancer [J].Cell Biol. Toxicol.， 2021，37（6）：811-831.