基于PTN系统分析不同种植地转基因水稻种子可培养内生真...

奥鹏网院作业

基于PTN系统分析不同种植地转基因水稻种子可培养内生真菌菌群的多样性
赵艳，王天圻，朱军莉

（浙江工商大学食品与生物工程学院，杭州 310018）

摘要：【目的】比较不同种植地、不同品种来源的转基因（genetically modified，GM）水稻近等基因系种子内生真菌的丰富度和多样性，探讨变异发生的影响因素，为研究GM水稻种子内生真菌种群结构的非预期变异提供科学基础。【方法】收集转2mG2-epsps抗草甘膦水稻株系（transgenic line，T），及相应的对照样本亲本品种（parent variety，P）和非转基因组培再生株系（non-transgenic regeneration line from tissue culture，NR），建立亲本对照-转基因株系-非转基因对照（parent control-transgenic plant line-non-transgenic control，PTN）近等基因系。粳稻品种日本晴亲本（P1）及其相应的T16和NR25组成P1近等基因系（P1 near-isogenic line，P1L），粳稻品种PJ574亲本（P2）及其相应的T23和NR18组成P2近等基因系（P2 near-isogenic line，P2L）。在海南省（Hainan province，H）和浙江省富阳市（Fuyang，Zhejiang province，F）两地种植并收获水稻种子，采用组织培养法分离内生真菌，通过形态学和分子生物学相结合的方法对菌株分类鉴定。以分离率（isolation rate，IR）、分离频率（isolation frequency，IF）、丰富度指数Margalef（D）、多样性指数Shannon-Wiener（H'）和均匀度指数Evenness（E）反映真菌的结构和分布，以Sorenson相似性系数（Cs）和Fisher精确检验（Fisher's exact test）结果描述水稻样本之间的内生真菌菌群组成差异。【结果】从海南省和浙江省富阳市两地的P1L和P2L近等基因系水稻种子样本中共分离121株内生真菌，鉴定为15个属，其中，Curvularia、Dendryphiella、Epicoccum、Fusarium为优势菌属，Fusarium为两地共有的优势菌属。F地水稻种子内生真菌总IR（4.61%）是H地样本（0.83%）的5.05倍。F地样本内生真菌类群的丰富度指数（2.29），Shannon-Wiener（H'）多样性指数（1.63）和均匀度指数（0.66）均显著高于H地样本。H地种植的P1L和P2L种子内生真菌菌属组成的相似性系数（Cs）为0.615，F地种植时二者的Cs=0.737，均为中等相似，Fisher精确检验分析表明二者无显著性差异（P＞0.05）。除HT16外的GM株系与相应亲本种子内生真菌类群的相似性系数为0.500—0.667，均为中等相似。但与相应亲本相比，GM水稻株系种子内生真菌的IF和菌属数存在显著性非预期变异，其中IF变异方向和幅度因种植地不同而波动。而且FT16增加2个菌属Phaeosphaeria和Nigrospora，HT23增加了1个菌属Letendraea helminthicola，FT23增加了Curvularia和Cladosporium 2个菌属，这些GM株系种子内生真菌菌属增多的变异来源于转基因插入突变，需要关注其安全性。而GM株系种子内生真菌种类减少的变异来源于组织培养无性系变异，较为安全。水稻种子真菌IF变异幅度及排序：P1和P2品种两地之间差异（30.58%）＞P1和P2品种之间差异（27.28%）＞NR株系变异（23.14%）＞GM株系变异（22.32%）。种子内生真菌总属数变异幅度及排序：P1和P2品种两种植地之间差异（9）＞P1和P2品种之间差异（8）=NR株系变异（8）=GM株系变异（8）。【结论】水稻种子内生真菌丰富多样，可培养内生真菌绝大多数为子囊菌。不同种植地域水稻种子内生真菌菌群组成存在差异，镰刀菌（Fusarium）为海南省和浙江省富阳市两地水稻种子的共有优势内生菌属。转基因技术和组织培养技术均显著影响水稻种子内生真菌类群的结构，但其非预期变异效应低于种植地域影响和品种差异，其中GM水稻种子内生真菌菌属种类增多的变异主要来源于转基因插入突变效应，需安全性评估。

关键词：转基因水稻种子；可培养内生真菌；非预期变异；PTN系统；多样性

0 引言
【研究意义】水稻是最重要的主粮作物之一，世界上一半以上的人口以稻米为食。基因工程已被用于改善包括水稻、玉米和小麦在内的多种作物的农艺性状和营养价值。评估转基因植物的安全性对于其商业化至关重要，实质等同性是基本原则[1]，公众对转基因（genetically modified，GM）作物非预期效应影响人类健康和生态环境的担忧是当前GM作物商品化应用的障碍[2]。研究GM水稻近等基因系的种子内生真菌多样性为进一步评估有关非预期效应奠定科学基础，对农业生物技术的健康发展意义重大。【前人研究进展】GM作物的非预期效应一直是国内外研究热点。Jiao等[3]研究表明GM水稻种子某些营养成分的含量存在非预期变异。赵艳等[2]利用PTN样本系统研究发现GM水稻种子蛋白质非预期变异主要来源于组织培养变异，基因插入的影响较小。Park等[1]研究证实水稻种子的营养品质不受基因插入的影响，并且表明环境因素对稻谷营养成分的影响比遗传背景的影响更大。GM植物一些代谢物的非预期变异可能影响植物体内环境进而影响植物内生微生物的定植。寄居在植物组织和器官内部，并且不会对宿主造成损害的真菌和细菌等微生物被称为内生菌或内生微生物[4]。前人研究表明转Bt植物可直接或间接地影响农作物周围环境微生物[5]。Pan等[6]报道内生真菌的多样性与玉米植物的基因型有关。刘月廉等[7]发现与亲本相比，GM水稻根部内生真菌分离率、种类多样性以及真菌侵染率和侵染指数无显著差异，转基因对水稻根部内生真菌分布影响不显著。【本研究切入点】将种植地域环境和转基因因素结合探讨GM水稻种子内生菌类群的研究仍鲜见报道。植物内生菌可在植物系统间进行内部传播，水稻种子既是内生菌的携带者和传播者，又是主食稻米的直接原料，内生真菌菌群变异可能涉及病害菌和产毒菌，从而影响水稻的种植生产和食用安全。【拟解决的关键问题】本研究以粳稻品种日本晴（P1）和PJ574（P2）为材料，收集亲本品种（parent variety，P）、转基因株系（transgenic line，T）、及非转基因组培再生株系（non-transgenic regeneration line from tissue culture，NR）等样本，分别建立水稻PTN近等基因系。分析不同种植地、不同品种来源的PTN近等基因系种子内生真菌群落结构特征和生物多样性，比较样本间内生真菌种群结构的差异，探究种植地域环境、品种遗传背景对水稻种子内生真菌多样性的影响贡献大小，为从种子内生菌分布角度客观评估GM水稻的非预期效应奠定基础。

1 材料与方法
1.1 供试材料
水稻转基因株系（T）以及非转基因对照样本（P、NR）的来源：水稻亲本（P）→愈伤组织的诱导培养和感受态细胞制备→外源基因2mG2-epsps转化愈伤组织→草甘膦抗性愈伤组织的筛选→再生植株→转基因株系（T）和非转基因组培再生株系（NR）。具体包括粳稻日本晴品种（P1）及其相应的GM株系T16和组培对照株系NR25组成的P1近等基因系（P1 near-isogenic line，P1L）；粳稻品种PJ574（P2）及其相应的GM株系T23和组培对照株系NR18组成的P2近等基因系（P2 near-isogenic line，P2L）。

供试水稻样本于2017年12月种植于海南省（Hainan province，H），2018年5月收获。2018年6月种植于浙江省富阳市（Fuyang，Zhejiang province，F），11月收获。收获种子自然晾干后，4℃保存，种子收获后2个月内完成微生物分离试验，尽量减少储存时间的影响。

1.2 水稻种子可培养内生真菌的分离纯化
水稻种子表面消毒：水稻种子脱壳后，参照CUI等[8]方法对水稻种子表面消毒，用无菌水冲洗种子5次，洗去表面灰尘杂质，再依次使用70%乙醇浸泡3 min、5%次氯酸钠浸泡5 min、70%乙醇清洗30 s、最后无菌水冲洗6次。取最后一次清洗液100 μL涂布在PDA培养基，28℃培养7 d，作为检验种子表面彻底消毒的对照。

内生真菌的分离：采用组织块分离法[8]，将表面消毒的种子放到垫有滤纸的无菌平板上，超净台上稍晾干，切成两段，切面朝下半埋接种于PDA培养基。在28℃恒温培养箱中培养6—10 d，待组织块周围长出菌丝体后，挑取菌丝接到PDA培养基上纯化培养。

1.3 水稻种子可培养内生真菌的鉴定
无菌条件下，将纯化后的真菌续接到PDA固体培养基，28℃培养7 d后观察菌落形态并拍照记录。刮取足量的菌丝体到灭菌洁净研钵中，用液氮充分研磨。采用CTAB法提取菌体DNA[9]，以DNA为模板，以ITS1（5′-TCCGTAGGTGAACCTGCGG-3′）和ITS4（5′-TCCTCCGCTTATTGATATGC-3′）为引物进行PCR扩增[8]，PCR反应体系（50 μL）：5 μL 10×PCR buffer（含Mg2+）、1 μL上游引物、1 μL下游引物、3 μL DNA模板、4 μL dNTP、1 μL Taq酶（最后加），dH2O补足50 μL。PCR扩增程序为94℃ 5 min；94℃ 45 s，54℃ 30 s，72℃ 90 s，32个循环；72℃ 10 min。PCR产物经1.0%琼脂糖凝胶电泳检测。

PCR产物送生工生物工程（上海）股份有限公司测序。用Blast软件将测序结果与GenBank中已知菌株的ITS rDNA序列进行同源性比对鉴定。

1.4 数据处理与统计分析
对稻米内生真菌的分离率（isolation rate，IR）、分离频率（isolation frequency，IF）进行统计，用Margalef指数（D）分析真菌类群的丰富度，Shannon-Weiner多样性指数（H'）、Evenness均匀度指数（E）及Sorenson相似性系数（Cs）分析内生真菌的多样性及分布特征[10]。

分离率IR（%）=（稻米样本长菌组织块数/稻米全部样本组织块数）×100。分离率是衡量植物组织中内生真菌的定植丰度和每个组织块受侵染的发生频率。

分离频率IF（%）=（Ni/N）×100。其中：Ni为某种稻米内生真菌的菌株数；N为全部内生真菌菌株数。

IF＞10%的内生真菌为优势菌，1%＜IF≤10%的内生真菌为常见菌，IF≤1%为稀有菌[11]。当N取值为所有样本内生真菌合集时，某种真菌相应占比用IF表示；当N取值为部分特定样本组合时某种真菌相应占比称为相对分离频率（relative frequency，RF），便于样本间的直观比较。

Margalef指数D=（S-1）/ln N。其中：S为类群总数；N为内生真菌总数目。

Shannon-Wiener多样性指数H'=-width=22,height=19（Ni/N）×ln（Ni/N）。其中，k为所有真菌种类的数目；Ni为某种真菌的数目；N为内生真菌总数目。

均匀度指数E=H'/lnS。其中：H'为Shannon-Wiener指数，S为群落中内生真菌的类群总数。

Sorenson相似性指数Cs=2j/（a+b）。式中：j为两种植地（或品系或不同基因型）共有的种数或属数；a为一个种植地（或品系或不同基因型）内生真菌的种数或属数；b为另一种植地（或品系或不同基因型）内生真菌的种数或属数。当Cs为0.00—0.25时，极不相似；当Cs为0.25—0.50时，中等不相似；当Cs为0.50—0.75时，中等相似；当Cs为0.75—1.00时，极相似。

利用软件spss17.0进行Fisher精确检验（Fisher's exact test）分析比较样本间真菌菌群组成的差异显著性。

2 结果
2.1 H和F种植地水稻种子可培养内生真菌的组成
接种后在28℃培养7 d，洗涤液对照无菌落长出（图1-a），说明试验稻米材料表面已彻底消毒。待组织块在PDA中长出较明显菌落时（图1-b—图1-e）续接纯化培养并进行菌落形态观测。从不同种植地（H、F）PTN系统样本组成的2个近等基因系2 400个稻米组织块中共分离内生真菌121株，总分离率（IR）为5.04%（表1），根据形态结构合并为43种真菌（图2）。以提取的真菌DNA为模版，ITS1和ITS4为引物的PCR扩增条带约在550 bp，测序后进行BLAST分析，将相似度＞97%的鉴定为同属，＞99%的鉴定为同种[12]，因此，将43个菌株归为3纲、4目、7科、15属（表1）。

width=418.5,height=97.9
a：对照；b、c、d、e：菌落分离图 a: control; b, c, d, e: colony separation diagram

图1 部分水稻种子内生真菌分离图

Fig. 1 Partial separation mapping of endophytic fungi from rice seeds

width=479.1,height=579.35
上：培养皿正面；下：培养皿反面

Upper: Petri dish front; Bottom: Petri dish reverse side

图2 水稻种子内生真菌在PDA平板上的菌落形态

Fig. 2 Colony morphology of endophytic fungi from rice seeds on PDA plates

2.2 不同种植地水稻种子内生真菌的多样性比较
2.2.1 海南省和浙江省富阳市两地水稻样本间的综合比较 海南省种植（H）的水稻样本共分离内生真菌20株，归为9个属，其中8个属为子囊菌纲，1个属为半知菌纲。浙江省富阳市（F）种植地分离内生真菌101株，归为12个属，全部属于子囊菌纲。可见水稻种子可培养内生真菌大多为子囊菌纲。

表1 水稻种子可培养内生真菌的鉴定结果

Table 1 Identification results of culturable endophytic fungi in rice seeds




H地水稻种子内生真菌的总IR（0.83%）显著低于F地样本（4.61%）。以两种植地所有水稻样本内生真菌菌株总数作为N值，计算不同样本中各内生真菌属的IF（表2），分析可知H和F两地内生真菌菌群的组成既有共性，又存在一定的差异。两地的共有真菌包括Curvularia、Phaeosphaeria、Setophoma、Nigrospora、Cladosporium和Fusarium6个属，共46株，分别占H和F地总菌株数的75.00%和30.70%，占两地总菌株38.02%。H地特有Alternaria、Letendraea helminthicola和Pleosporales 3个菌属共5株，占两地总菌数4.13%。F地特有真菌包括6个属共70株，分别为Ophiosphaerella、Dendryphiella、Epicoccum、Bipolaris、Pithomyces和Leptosphaeria，占比高达F地总菌株的69.3%和两地总菌数的57.85%。表明水稻种子内生真菌种群具有明显的地理分布特性。

H地6个样本内生真菌的IF和所含真菌类群的Margalef（D）、Shannon-Wiener（H'）以及均匀度指数（E）均比对应F地的低（表2和表3），表明F地内生真菌要比H地更加丰富和多样。两地水稻种子样本内生真菌类群属水平上Cs为0.57，为中等相似（表4），并且Fisher精确检验结果显示二者分析两地真菌菌群组成差异显著（P＜0.05）（表5），表明水稻种子内生真菌群落的组成在一定程度上受地理环境影响。

2.2.2 亲本品种及其PTN近等基因系种子内生真菌分布的地域差异比较 按地域分组计算RF值比较两地样本的内生优势菌属（RF＞10%，图3），Fusarium是两地水稻种子共有优势菌属。H地样本有3个优势菌属，其中Curvularia菌属在F地样本中丰度降为常见菌属，Alternaria在F地未检出。F地独有优势菌属Dendryphiella和Epicoccum广泛存在于6个样本中（表2），可见不同种植地域水稻种子内生真菌的优势菌群组成存在差异。

表2 不同种植地水稻种子内生真菌分离频率比较（IF）（%）

Table 2 Comparison of isolation frequency (IF) (%) of endophytic fungi in rice seeds from different planting sites (%)


HP1L和FP1L：海南省和浙江省富阳市种植的P1近等基因系，包括亲本品种P1及其相应的转基因株系T16和非转基因组培对照株系NR25；HP2L和FP2L：海南省和浙江省富阳市种植的P2近等基因系，包括亲本品种P2及其相应的转基因株系T23和组培对照株系NR18；—：未检出。下同

HP1L and FP1L: P1 near-isogenic line planted in Hainan and Fuyang, Zhejiang, containing parent variety P1, and its corresponding genetically modified line T16 and non-transgenic control line NR25 regenerated from tissue culture; HP2L and FP2L: P2 near-isogenic line planted in Hainan and Fuyang, Zhejiang, containing parent variety P2, and its corresponding genetically modified line T23 and non-transgenic control line NR18 regenerated from tissue culture; —: no detection. The same as below

width=472.1,height=182.85
图3 海南省（H）和浙江省富阳市（F）种植水稻种子内生真菌的相对分离频率（RF）

Fig. 3 Relative isolation frequency ( RF ) of endophytic fungi from rice seeds grown in Hainan (H) and Fuyang, Zhejiang (F)

表3 海南省（H）和浙江省富阳市（F）种植的PTN近等基因系水稻种子内生真菌菌群多样性指数

Table 3 The endophytic fungal diversity indexes of rice seeds from PTN near-isogenic lines in Hainan (H) and Fuyang, Zhejiang (F) cultivation sites


P1L：P1近等基因系，包括亲本品种P1及其相应的转基因株系T16和非转基因组培对照株系NR25；P2L：P2近等基因系，包括亲本品种P2及其相应的转基因株系T23和组培对照株系NR18；—：未检出。下同

P1L: P1 near-isogenic line, containing parent variety P1, and its corresponding genetically modified line T16 and non-transgenic control line NR25 regenerated from tissue culture; P2L: P2 near-isogenic line, containing parent variety P2, and its corresponding genetically modified line T23 and non-transgenic control line NR18 regenerated from tissue culture ; —: no detection. The same as below

将两地种植的2个PTN近等基因系共6个供试水稻样本种子可培养内生真菌菌群组成和多样性指数汇总于表3。观察相同品种水稻种植于不同地点时的种子内生真菌类群丰富度和多样性。对P1品种而言，HP1和FP1样本表现为“无”和“有”的极端差异，可见天然亲本品种P1种子内生真菌类群分布受种植地域影响很大。将观测样本群体扩大至P1近等基因系P1L，除均匀度指数（E）外，浙江省富阳市种植的FP1L近等基因系分离的真菌株数、菌属数、丰富度指数（D）和多样性指数（H'）都高于相应海南省种植样本（表3），这与上述H和F两种植地所有水稻样本的内生真菌多样性变化趋势一致。对P2品种而言，虽然FP2分离的内生真菌株数和属数比HP2多，但多样性指数均低于相应的HP2样本；将观测样本群体扩大至P2近等基因系P2L，FP2L分离的内生真菌株数和菌属数仍比HP2L多，Margalef（D）指数1.52也大幅高于海南省种植HP2L群体的0.30。

2.3 不同品种水稻种子内生真菌的多样性比较
综合表2和表3结果可分析不同样本间种子内生真菌的组成和分布差异，探测水稻种子非预期变异的特征和变异来源。以亲本品种P1和P2为比对组，分析天然水稻品种之间种子内生真菌菌群结构和多样性变异程度。

海南省种植地样本中，HP1未检出真菌，HP2分离到内生真菌6株，归为Curvularia、Phaeosphaeria、Nigrospora、Alternaria和Fusarium5个属（表2）。可见在海南省种植时P2品种内生真菌类群比P1丰富得多。浙江省富阳市种植地样本中，FP1分离37株内生真菌归为9个属，FP2获得10株内生真菌归为6个属，FP2的IF（8.26%）显著低于FP1（30.58%）。Curvularia、Setophoma、Ophiosphaerella和Pithomyces 4个为FP1特有菌属，Nigrospora为FP2特有菌属；FP1和FP2共有菌属5个，包括Dendryphiella、Bipolaris、Epicoccum、Leptosphaeria和Fusarium。FP2内生真菌的多样性指数Margalef（D）、Shannon-Wiener（H'）以及均匀度指数（E）均低于FP1，下降幅度分别达27.08%、39.39%和25.00%，说明浙江省富阳市种植的P1品种内生真菌类群比P2更丰富多样。

综合两地，在属水平，P1和P2之间种子内生真菌菌群的相似性系数（Cs）分别为0（H）和0.67（F）（表4），表现极不相似和中等相似；近等基因系P1L和P2L之间内生真菌菌群的Cs分别为0.62（H）和0.74（H），为中等相似，Fisher' exact test检验显示二者差异均不显著（P＞0.05）（表5）。

在内生真菌类群组成上，HP2和FP2样本都分离出了Nigrospora和Fusarium，说明这两个菌属是P2品种的固有菌属，只是Fusarium菌属丰度因地域而表现不同，FP2中Fusarium菌属的IF（1.65%）是HP2的2倍（表2）。综合表2和表3，可知HP2中的3个菌属Curvularia、Phaeosphaeria、Alternaria和FP2中的Dendryphiella等4个菌属在水稻种子内定植主要受地域影响。对P1品种而言，只在F地分离到9个属真菌，种植地域影响更为突出。

2.4 GM水稻种子内生真菌的非预期变异特征
海南省种植的GM水稻株系T16（HT16）种子共分离7株内生真菌，总IF为5.79%，归为Curvularia、Fusarium、Pleosporales、Cladosporium和Alternaria 5个属（表2），相应亲本HP1未检出真菌，二者相似性系数为0，完全不等同。浙江省富阳市种植时，GM水稻株系FT16种子内生真菌总IF（14.05%）极显著低于相应亲本FP1（30.58%），丰富度指数（D）、多样性指数（H'）和均匀度指数（E）也都低于FP1样本。菌群结构上，FT16与FP1样本共有Ophiosphaerella、Dendryphiella、Epicoccum和Fusarium 4个内生真菌属，未检出FP1携带的Bipolaris、Leptosphaeria、Curvularia、Setophoma和Pithomyces菌属，却额外增加了2个特有菌属Phaeosphaeria和Nigrospora。FT16与FP1内生真菌类群的相似性系数0.53，为中等相似（表4），且二者内生真菌的组成不具有显著差异（P=0.23＞0.05，表5）。

海南省种植的GM水稻株系T23（HT23）种子内生真菌总IF（2.48%）显著低于亲本对照HP2（4.96%）。菌群组成上，与亲本HP2相比，HT23携带的3个菌属中Letendraea helminthicola为新增特有菌属，Curvularia和Fusarium与HP2所共有，此外未检出HP2特有菌属Phaeosphaeria、Nigrospora和Alternaria，二者的相似性系数为0.50，中等相似，且不具有显著差异（表5）。浙江省富阳市种植的FT23样本内生真菌总分离率13.22%显著高于FP2。FT23分离获得6个属共16株真菌，其中Dendryphiella、Epicoccum、Leptosphaeria和Fusarium4个菌属与亲本FP2共有，Curvularia和Cladosporium为新增特有菌属，比FP2减少了2个菌属Bipolaris和Nigrospora。FT23种子内生真菌类群与FP2的相似性系数0.67，中等相似，二者内生真菌的组成差异不显著（P=0.05）。

综合两地样本，GM水稻品系T16种子固有的内生真菌只有Fusarium属，在H和F两地的IF分别为1.65%和4.13%；T23种子固有的内生真菌包括Curvularia和Fusarium。Nigrospora和Fusarium是亲本品种P2种子固有内生真菌（表2）。

表4 不同种植地水稻PTN近等基因系样本种子内生真菌菌群的相似性系数（Cs）

Table 4 Similarity coefficients (Cs) of endophytic fungi of rice seeds from PTN near-isogenic line samples in different cultivation places


表5 不同种植地水稻PTN近等基因系样本种子内生真菌菌群组成的差异性分析（P）

Table 5 Analysis the differences (P) composition of endophytic fungi of rice seeds from PTN near-isogenic line samples in different cultivation places


2.5 GM水稻种子内生真菌类群的非预期变异技术来源解析
两地种植的P1L和P2L近等基因系中，与相应亲本相比，2个GM株系种子内生真菌总分离频率（IF）在海南省和浙江省富阳市的变异方向相反，但在相同种植地时又分别与相应组培对照NR株系变异方向相同（表2）。以T16为例，HT16种子内生真菌总IF（5.79%）显著高于HNR25（1.65%），变异方向一致只是大小不同，其中GM水稻品HT16种子内生真菌分离频率增加的变异效应中转基因插入突变效应＞组织培养无性系变异效应。但浙江省富阳市种植时，FT16种子内生真菌总IF（14.05%）比亲本FP1降低的幅度显著低于FNR25（9.09%），因此，推测HT16种子内生真菌分离频率降低的变异效应中转基因插入突变效应＜组织培养无性系变异效应。对GM水稻品系T23而言，海南省种植的HT23和组培对照HNR18的种子内生真菌总IF均显著低于亲本对照HP2（4.96%），二者变异方向一致而且HNR18样本变异幅度更大，似乎这种内生真菌总分离率降低的变异效应来源于组织培养无性系变异。但浙江省富阳市种植的FT23样本内生真菌总IF13.22%极显著高于FP2（8.26%），组培对照FNR18与FP2亲本对照相同（表2），又可将该变异效应归因于转基因插入突变。综合2个种植地变异方向相反且技术来源不甚明确的分析结果，认为GM水稻品系种子内生真菌总IF的波动变异受种植地域环境影响较大，非预期变异效应较小。

菌群组成和结构上，HT16种子携带内生真菌Curvularia等5个属，HNR25却只有Setophoma 1个属（表2）；FT16额外增加的2个菌属Phaeosphaeria和Nigrospora不存在于其相应组培对照样本FNR25中，说明T16种子定植真菌类群增加的变异主要为转基因插入突变效应。与FP1相比，FT16种子内生真菌减少了Bipolaris等5个菌属，组培对照样本FNR25却减少了7个菌属，而且其中5个菌属与FT16的一致，由此可知GM水稻株系T16种子内生真菌类群的消减变异效应为组织培养无性系变异。同理，HT23与亲本HP2相比增加了1个菌属Letendraea helminthicola并且该菌属只出现在该样本中，浙江省富阳市种植的FT23比亲本FP2增加了Curvularia和Cladosporium 2个菌属，这种内生真菌菌属增多的变异并未在相应组培参照样本中出现，可判断变异来源于转基因插入突变，需要重点关注。海南省种植的HT23比HP2减少3个菌属Phaeosphaeria、Nigrospora和Alternaria，HNR18比HP2减少4个菌属Curvularia、Phaeosphaeria、Nigrospora和Fusarium，浙江省富阳市种植的FT23与组培参照样本FNR18一样也减少了Bipolaris和Nigrospora 2个菌属，可知GM水稻品系T23种子内生真菌减少的变异主要受组织培养无性系变异效应的影响。

2.6 不同样本水稻种子内生真菌分布的聚类分析及影响因素主次评价
在属水平上，根据种子可培养内生真菌的种类和丰度进行聚类分析（图4），可见两地种植的12个水稻样本中，聚类距离最近的样本组包括FNR25和FNR18、FT16和HP2、HP1和HNR25。分析样本具体来源可知，HP1和HNR25为海南省种植的P1品种及其组培品系，遗传关系相近外；FNR25和FNR18对应P1和P2品种来源的组培品系之间；FT16和HP2分别属于不同近等基因系、不同种植地域、GM水稻品系和天然亲本P2，实际遗传关系较远，可见种子内生真菌类群分布受品种、种植地、转基因技术和组织培养技术改良的综合影响。在根据样本聚类的3个大组中，除FP1内生真菌类群和丰度最大单独归为一个大组外，FP2/FT23/FNR25/FNR18/FT16/HP2归为一组，HT16/HT23/HNR18/HP1/HNR25归为另一组，后2个大组近似于种植地域的样本集合，说明种植地域的影响较大。

width=375.35,height=373.95
图4 海南省（H）和浙江省富阳市（F）种植的6个水稻样本种子内生真菌组成的聚类分析

Fig. 4 The cluster analysis of endophytic fungi composition of rice seeds from the 6 samples planted in Hainan (H) and Fuyang, Zhejiang (F)

综合分析两地种植的2个PTN近等基因系共12个样本的变异情况，将变异幅度峰值采取同向取高，反向绝对值相加的计算方法，得出最大变异幅度并根据大小进行排序，真菌总分离频率（IF）排序结果：P1和P2品种两种植地之间变异幅度3.30%—30.58%（30.58%）＞P1和P2品种之间-22.32%—4.96%（27.28%）＞NR株系变异幅度-21.49%—1.65%（23.14%）＞GM株系变异幅度-16.53%—5.79%（22.32%），可见种植地域差异对水稻种子可培养内生真菌分离频率的影响最大，品种影响其次，组培技术的变异效应较小，转基因技术的变异效应最小。从携带内生真菌的总属数变异幅度排序结果：P1和P2品种两种植地之间变异幅度1—9（9）＞P1和P2品种之间-3—5（8）=NR株系变异幅度-7—1（8）=GM株系变异幅度-3—5（8），可见在属水平上，种植地域差异对水稻种子可培养内生真菌类群组成的影响最大，其他3个因素的影响较小且相当。

3 讨论
“实质等同性”原则一直是转基因作物安全性评价的重要指标，尚未见基于种子内生菌评估GM作物非预期变异的研究报道。植物内生菌与植物生长发育密切相关，进而可能影响农作物的产量以及食用安全。有研究表明Curvularia lunata和Nigrospora oryzae能抑制胚芽生长，并且减少洋葱根分生组织的有丝分裂指数[13]。

本研究发现水稻种子可培养内生真菌绝大多数为子囊菌纲，并且Fusarium为海南省和浙江省富阳市两种植地的优势菌属（图3）。Tian等[14]也发现Fusarium为稻米中常见菌株。Fusarium是水稻恶苗病的病原体[15]，也可能产生威胁食品安全的真菌毒素[16]。本研究发现地域对内生真菌的影响较为明显，H和F两地内生真菌菌群的组成既有共性，又存在一定的差异。Alternaria作为H地特有优势菌株未在F地中分离得到，Dendryphiella和Epicoccum作为在F地特有的优势菌株存在于6个稻米样本中（表2），表明内生真菌的分布具有地域的专化性，优势菌群也随着地理环境的变化而改变。这与李盼盼的研究结果一致[11]。并且F地6个样本中的IF、丰富度指数（D）、多样性指数（H'）以及均匀度指数（E）均高于H地，表明F地稻米内生真菌更加丰富，薛庆婉等[17]研究也发现地理环境是影响植物内生真菌组成、丰富度和优势菌群的重要因素。本研究发现H和F地水稻种子内生真菌的组成存在显著性差异，差异的原因可能涉及土壤化学成分组成[14]，环境条件和气候的影响[13]。

理论上，GM水稻种子的非预期变异效应包括公认安全的组织培养无性系变异和风险性较大的转基因插入突变，技术来源分别为GM株系培育过程中的组织培养和转基因插入整合。同等条件下，本研究PTN样本系统中NR的性状表现可作为解析GM水稻品系非预期变异技术来源的有效参照[2]。在H和F种植地中，某些菌属表现明显的品种特异性，如Nigrospora为P2特有菌株（表2），原因可能与水稻品种的遗传背景相关，也可能受感染种子从母代到子代的垂直遗传[18]。GM水稻T16和T23种子内生真菌分离频率的变异效应即包括组织培养无性系变异，又包括转基因插入突变，但变异幅度显著小于P1品种的种植地域差异，也小于浙江省富阳市种植地的FP1和FP2的品种间变异（表2）。可推定GM水稻品系种植内生真菌IF的变异幅度在安全范围内。与相应亲本相比，FT16额外增加2个菌属Phaeosphaeria和Nigrospora，HT23与亲本HP2相比增加了1个独有菌属Letendraea helminthicola，FT23比亲本FP2增加了Curvularia和Cladosporium 2个菌属，基于PTN系统分析可知这些内生真菌菌属种类增多的变异均来源于转基因插入突变效应，需要重点关注其安全性。

与前人关于转基因技术对包括水稻在内的植物根际土壤真菌群落、根际细菌群落和根内生真菌群落结构的影响不明显[19-22]的报道不同，本研究发现转基因技术和组织培养技术均能显著影响水稻种子内生真菌类群的结构。两地种植的PTN系统样本的最大变异幅度排序结果可知种植地域对水稻种子可培养内生真菌丰度和类群组成的影响最大，品种影响次之，组织培养无性系变异效应影响较小，转基因插入突变效应的影响最小。这与前人对GM马铃薯和GM水稻根际微生物群落分布的研究结果类似，与种植地点、土壤类型、植物生长阶段、植物品种相比，转基因插入突变对马铃薯根际细菌群落的影响较小，其中种植地域影响要大的多[22-25]。结合聚类分析（图4）和P1品种在海南省和浙江省富阳市种植时种子内生真菌多样性和丰富度的极端差异表现（表2），可明确种植地域对水稻种子内生真菌类群的影响最为显著。

4 结论
水稻种子内生真菌丰富多样，可培养内生真菌绝大多数为子囊菌纲。不同种植地域水稻种子内生真菌菌群组成存在差异，Fusarium为海南省和浙江省富阳市两种植地的水稻种子的共有优势内生菌属。转基因技术和组织培养技术均显著影响水稻种子内生真菌类群的结构，其中GM水稻种子内生真菌菌属种类增多的变异主要来源于转基因插入突变效应，需评估其安全性。对水稻种子可培养内生真菌丰度和类群组成的影响大小排序为种植地域＞品种差异＞组织培养无性系变异效应＞转基因插入突变效应。

References

[1] Park S Y, Kim J K, Jang J S, Lee S Y, Oh S, Lee S M, Yang C I, Yeo Y. Comparative analysis of nutritional composition between the disease-resistant rice variety OsCK1 and conventional comparators. Food Science and Biotechnology, 2015, 24(1): 225-231.

[2] 赵艳, 李燕燕. 利用PTN系统快速解析转基因水稻种子可溶性蛋白非预期变异的来源. 中国农业科学, 2015, 48(22): 4397-4407.

Zhao Y, Li Y Y. Rapid identify the sources of unintended variations of seed soluble proteins in genetically modified rice by using PTN system. Scientia Agricultura Sinica, 2015, 48(22): 4397-4407. (in Chinese)

[3] Jiao Z, Si X X, Li G K, Zhang Z M, Xu X P. Unintended compositional changes in transgenic rice seeds (Oryza sativa L.) studied by spectral and chromatographic analysis coupled with chemometrics methods. Journal of Agricultural and Food Chemistry, 2010, 58(3): 1746-1754.

[4] Almeida T T, Orlandelli R C, Azevedo J L, Pamphile J A. Molecular characterization of the endophytic fungal community associated with Eichhornia azurea (Kunth) and Eichhornia crassipes (Mart.) (Pontederiaceae) native to the Upper Paraná River floodplain, Brazil. Genetics & Molecular Research Gmr, 2015, 14(2): 4920-4931.

[5] Wang X, Cai M Y, Zhou Y. Biological influence of cry1Ab gene insertion on the endophytic bacteria community in transgenic rice. Turkish Journal of Biology, 2018, 42(3): 231-239.

[6] Pan J J, Baumgarten A M, May G. Effects of host plant environment and Ustilago maydis infection on the fungal endophyte community of maize (Zea mays). New Phytologist, 2008, 178(1): 147-156.

[7] 刘月廉, 郭建夫, 袁红旭, 李玥仁, 李伟, 黄永相. 两个转抗真菌基因水稻根部内生真菌群落变化的分析. 生物过程, 2012, 2: 105-110.

Liu Y L, Guo J F, Yuan H x, Li Y R, Li W, Huang Y X. Analysis of endophytic fungal community changes in roots of two transgenic antifungal genes. Biology Process, 2012, 2: 105-110. (in Chinese)

[8] Cui J L, Guo T T, Ren Z X, Zhang N S, Wang M L. Diversity and antioxidant activity of culturable endophytic fungi from alpine plants of Rhodiola crenulata, R. angusta, and R. sachalinensis. PLoS ONE, 2015, 10(3): e0118204.

[9] Weising K, Nybom H, Pfenninger M, Wolff K, Meyer W.DNA Fingerprinting in Plants and Fungi.Florida: CRC press, 1994.

[10] 柴新义, 柴钢青, 向玉勇, 张微微, 殷培峰. 青檀叶片内生和附生真菌组成及生态分布. 生态学报, 2016, 16: 5163-5172.

Chai X Y, Chai G Q, Xiang Y Y, Zhang W W, Yin P F. Composition and ecological distribution of endophytic and epiphytic fungi in leaves of P. chinensis. Chinese Journal of Ecology, 2016, 16: 5163-5172. (in Chinese)

[11] 李盼盼, 袁晓龙, 李金海, 申国明, 高林, 李青城, 高加明, 张鹏. 湖北烟草内生真菌生物多样性和种群结构分析. 微生物学报, 2018, 58(10): 169-179.

Li P P, Yuan X L, Li J h, Shen G M, Gao L, Li Q C, Gao J M, Zhang P. Biodiversity and population structure analysis of endophytic fungi from tobacco in Hubei Province. Acta Microbiologica Sinica, 2018, 58(10): 169-179. (in Chinese)

[12] 王梦亮, 贾岩, 崔晋龙, 王俊宏.锁阳及其寄主白刺内生真菌的生态分布及遗传关系. 应用生态学报, 2017, 28(3): 815-820.

Wang M L, Jia Y, Cui J L, Wang J H. Ecological distribution and genetic relationship of endophytic fungi from Cynomorium and its host white thorn. Chinese Journal of Applied Ecology, 2017, 28(3): 815-820. (in Chinese)

[13] Rajagopal K, Meenashree B, Binika D, Joshila D, Tulsi P S, Arulmathi R, Kathiravan G, Tuwar A. Mycodiversity and biotechnological potential of endophytic fungi isolated from hydrophytes. Current Research in Environmental and Applie Myclogy, 2018, 8(2): 172-182.

[14] Tian X L, Cao L X, Tan H M, Zeng Q G, Jia Y Y, Han W Q, Zhou S N. Study on the communities of endophytic fungi and endophytic actinomycetes from rice and their antipathogenic activities in vitro. World Journal of Microbiology & Biotechnology, 2004, 20(3): 303-309.

[15] 王拱辰, 陈鸿逵. 水稻恶苗病病原菌的研究. 植物病理学报, 1990, 2: 93-97.

Wang G C, Chen H K. Studies on pathogens of rice blight disease. Chinese Journal of Plant Pathology, 1990, 2: 93-97. (in Chinese)

[16] 史建荣, 刘馨, 仇剑波, 祭芳, 徐剑红, 董飞, 殷宪超, 冉军舰. 小麦中镰刀菌毒素脱氧雪腐镰刀菌烯醇污染现状与防控研究进展. 中国农业科学, 2014, 47(18): 3641-3654.

Shi J R, Liu X, Qiu J B, JI F, XU J H, DONG F, YIN X C, RAN J J. Research status and prevention and control of Fusarium toxin deoxynivalenol contamination in wheat. China Agricultural Science, 2014, 47(18): 3641-3654. (in Chinese)

[17] 薛庆婉, 牛永春, 邓晖. 常见瓜类植物的内生真菌多样性. 菌物学报, 2015, 34(2): 196-203.

Xue Q W, Niu Y C, Deng H. Endophytic fungal diversity of common melon plants. Journal of the Physiology, 2015, 34(2): 196-203. (in Chinese)

[18] 贾彤, 任安芝, 魏茂英, 尹立佳, 高玉葆. 不同传播方式的内生真菌感染对羽茅的生理生态影响. 植物生态学报, 2015, 39(1): 72-80.

Jia T, Ren A Z, Wei M Y, Yin L J, Gao Y B. Physiological and ecological effects of endophytic fungal infections with different propagation patterns on P. chinensis. Chinese Journal of Plant Ecology, 2015, 39(1): 72-80. (in Chinese)

[19] Danielsen L, Thurmer A, Meinicke P, Buee M, Morin E, Martin F, Pilate G, Daniel R, Polle A, Reich M. Fungal soil communities in a young transgenic poplar plantation form a rich reservoir for fungal root communities. Ecology and Evolution, 2012, 2(8):1935-1948.

[20] Heuer H, Kroppenstedt R M, Lottmann J, Berg G, Smalla K. Effects of T4 lysozyme release from transgenic potato roots on bacterial rhizosphere communities are negligible relative to natural factors. Applied and Environmental Microbiology, 2002, 68(3): 1325-1335.

[21] Silva R D A, Dourado M N, Fávaro L C L, Mendes R, Ferreira A, Araujo W L. Diversity of cultivated fungi associated with conventional and transgenic sugarcane and the interaction between endophytic trichoderma virens and the host Plant. PLoS ONE, 2016, 11(7): e0158974.

[22] Lee S H, Kim C G, Kang H. Temporal dynamics of bacterial and fungal communities in a genetically modified (GM) rice ecosystem. Microbial Ecology, 2011, 61(3): 646-659.

[23] Rasche F, Hodl V, Poll C, Kandeler E, Gerzabek M H,Elsas J D V, Sessitsch A. Rhizosphere bacteria affected by transgenic potatoes with antibacterial activities compared with the effects of soil, wild-type potatoes, vegetation stage and pathogen exposure. Fems Microbiology Ecology, 2010, 56(2): 219-235.

[24] Emilia H S, Wietse D B, Johannes V V, Kevin M C. A 3-year study reveals that plant growth stage, season and field site affect soil fungal communities while cultivar and GM-trait have minor effects. PLoS ONE, 2012, 7(4): e33819.

[25] Liu W, Lu H H, Wu W, Wei Q K, Chen Y X, Thies J E. Transgenic Bt rice does not affect enzyme activities and microbial composition in the rhizosphere during crop development. Soil Biology & Biochemistry, 2008, 40(2): 475-486.

Diversity of Endophytic Fungi in Transgenic Rice Seeds from Different Planting Sites Based on PTN System
ZHAO Yan, WANG Tianqi, ZHU JunLi

(College of Food Science and Biotechnology, Zhejiang Gongshang University, Hangzhou 310018)

Abstract: 【Objective】To provide scientific bases for study of the unintended variation in seed endophytic fungi community structure of genetically modified (GM) rice and explore the causing factors of variation, comparative analysis of the richness and biodiversity of seed endophytic fungi were conducted, using GM rice relative near-isogenic lines of different varieties in different cultivation sites as materials. 【Method】Collect transgenic rice line which harboring the glyphosate resistant 2mG2-epsps gene (T)and its corresponding parent variety (P) and non-transgenic tissue culture regeneration control line (NR) to make parent control-transgenic plant line-non-transgenic control (PTN) near-isogenic line. The parent japonica rice Nipponbare (P1) and its corresponding transgenic line T23 and the NR control line NR18 formed the P1 near-isogenic line (P1L), the P2 near-isogenic line (P2L) composed of parent japonica rice PJ574 (P2) and its corresponding transgenic line T23 and the NR control line NR18. All rice samples were planted in two plantations including Hainan province (H) and Fuyang city of Zhejiang province (F) and the resulting seeds were harvested. The rice endophytic fungi were isolated by tissue separation method, strains were classified and identified with morphology and molecular biology methods. The isolation rate (IR), isolation frequency (IF), richness Margalef index (D), diversity Shannon-Wiener index (H') and Evenness index (E) were used to reflect the structure and distribution of rice seed endophytic fungi, and Sorenson similarity coefficient (Cs) and Fisher’s exact test were employed to describe the composition difference of endophytic fungi between rice samples. 【Result】A total of 121 endophytic fungi strains were isolated from rice seed samples of P1L and P2L near-isogenic lines that harvested in Hainan (H) and Zhejiang (F) plantations, they were identified as 15 genera, of which Curvularia, Dendryphiella, Epicoccum, Fusarium were confirmed as the dominant flora, with Fusarium as the common dominant genus of both H and F plantations. The total RF (4.61%) of endophytic fungi from Zhejiang-grown rice seeds is 5.05 times than that of Hainan-grown rice samples (0.83%). The richness Margalef index (D=2.29), Shannon-Wiener diversity index (H'=1.63) and Evenness index (E=0.66) of endophytic fungi flora from rice seeds in F plantation were higher than that (D=1.67, H'=0.63, E=0.29) of samples in H plantation. The seed endophytic fungi communities of P1L and P2L in H plantation had a similarity coefficient of 0.615, and that was 0.737 in F plantation, showing moderate similarity, and the Fisher’s exact test analysis suggested that there were no significant difference (P＞0.05) between them. The GM lines except HT16 showed moderate similarity to their corresponding parent controls referring seed endophytic fungi communities, with similarity coefficient ranged between 0.500 and 0.667. But compared with their corresponding parents, the GM rice lines showed notably unintended variations in IF and genus numbers of endophytic fungi, the variation direction and amplitude of IF varied between different plantations. Also, the GM line FT16 increased 2 additional fungi genera Phaeosphaeria and Nigrospora, HT23 and FT23 respectively added 1 genus Letendraea helminthicola and 2 genera including Curvularia and Cladosporium. These genera increasing variations of endophytic fungi in GM rice seeds were derived from transgenic insertion mutation, whose safety needs to be focused on. While the genera decreasing variations of endophytic fungi in GM rice seeds were derived from somaclonal variation of tissue culture, which were safer. The variationamplitudeorder oftotal IF of rice endophytic fungi were as follows. Difference between H and F plantations referring to P1 and P2 (30.58%)＞Difference between P1 and P2 varieties (27.28%)＞Variation of NR lines (23.14%)＞Variation of GM lines (22.32%). The genus numbervariations of rice endophytic fungi ranked as follows. Difference between H and F plantations referring to P1 and P2 (9)＞Difference between P1 and P2 varieties (8) = Variation of NR lines (8) = Variation of GM lines (8). 【Conclusion】Rice seeds have abundant and diverse endophytic fungi, with majority of the culturable stains belong to Ascomycetes. Composition of rice seed endophytic fungi community shows geographical differences, and Fusarium is very common dominant genus in rice seeds grown both in Hainan and Zhejiang. The structure of endophytic fungi flora in rice seeds are affected by transgenic manipulation as well as tissue culture technology, while their unintended variation effects are less than that of rice growing locations and variety differences. The genera increasing variations of endophytic fungi in GM rice seeds are derived from transgenic insertion mutation, and the safety needs to be assessed.

Key words: transgenic rice seed; culturable endophytic fungi; unintended variation; PTN system; diversity