|
黄瓜ERF基因家族鉴定及其在雌花芽分化中的表达分析黄瓜ERF基因家族鉴定及其在雌花芽分化中的表达分析 潘健，温海帆，何欢乐，连红莉，王刚，潘俊松，蔡润 （上海交通大学农业与生物学院，上海 200240） 摘要：【目的】通过以黄瓜9930_V2版本基因组为参照进行生物信息学分析，对ERF基因家族在基因组中的数量、结构以及表达特征进行分析，为研究ERF转录因子在黄瓜雌花分化与发育中的作用提供数据支持。【方法】根据已报道的拟南芥ERF，利用黄瓜基因组数据库中9930_V2版本基因组进行BLAST比对，通过MEGA、MEME、TBtools、ExPASy等工具进行生物信息学分析。采用qRT-PCR方法检测不同性型黄瓜材料、雌花发育初期不同阶段中ERF基因家族成员的表达水平。采用酵母单杂交方法验证家族成员与乙烯响应元件GCC-box的互作。【结果】从黄瓜材料9930基因组中鉴定得到138个ERF基因家族成员，共分为10个亚族，编码氨基酸长度介于126—745。按照基因家族成员在染色体上的位置分布，将其命名为CsERF1-CsERF138。多序列比对和motif分析结果表明，黄瓜ERF基因家族均具有AP2/ERF结构域，其中4个成员具有B3结构域。表达分析结果显示，在不同性型材料中共有19个ERF家族成员差异表达，其中9个在FFMMAA基因型中高表达，10个在ffMMAA基因型中高表达。通过雌花芽发育初期ERF家族成员的表达趋势分析，发现31个ERF随子房发育表达上调，30个表达下调。初步证明CsERF9和CsERF31具有结合GCC-box元件的功能。【结论】从黄瓜基因组中鉴定出138个ERF基因家族成员，均拥有1个或多个AP2/ERF结构域；其中部分成员在不同性型材料中差异表达，并可能参与雌花分化初期的基因表达调控；部分成员具有结合保守元件GCC-box调控下游基因表达的功能。 关键词：黄瓜;乙烯信号;雌花发育;ERF基因；性别决定 0 引言【研究意义】黄瓜（Cucumis sativus L.）是重要的园艺作物之一，其栽培面积仅次于番茄和洋葱[1]，为世界第三大蔬菜。雌花子房是黄瓜经济价值的体现，其分化与发育决定了果实的产量和质量。黄瓜果实的形成经历了花芽性别分化和子房发育两个阶段，植物激素乙烯在这一过程中起至关重要的作用[2]。乙烯响应因子（ethylene response factor，ERF）家族是一类植物特异性转录因子家族，多个ERF家族成员在不同物种中参与到乙烯信号的响应与调控[3-5]。乙烯是影响黄瓜雌花分化与子房发育的重要因素，系统地对黄瓜基因组中ERF基因家族进行鉴定与分析，可以为深入探究黄瓜雌花分化与发育的分子机制提供数据基础。【前人研究进展】植物ERF基因家族成员数量众多，且在各种生物和非生物胁迫下发挥重要功能[6]。在一些茄科植物中，LeERF2/TERF2通过与GCC-box 结合，诱导乙烯的生物合成[3]。在香蕉（Musa acuminate）中，MaERF9通过结合MaACS1的启动子激活其转录，而MaERF11通过结合MaACS1和MaACO1的启动子，抑制其转录[4]。类似的，苹果（Malus domestica）中MdERF3和MdERF2通过结合MdACS1启动子，分别激活和抑制其转录[5]。这些研究结果说明ERF在乙烯合成过程中起重要作用。ERF1是较早被克隆并进行功能研究的ERF之一[7]。后续研究中，ERF1被报道通过与茉莉酸、生长素等互作，参与植物防卫反应[8-9]，也有报道显示ERF1通过调控下胚轴细胞壁增厚和伸长参与植物的光形态建成[10]。在黄瓜性别分化研究中，已报道CsERF110和CsERF31具有分别结合A基因（CsACS11）[11]和M基因（CsACS2）的启动子，激活其转录的功能[12]。【本研究切入点】在已有研究中，黄瓜ERF基因家族已基于9930_V1版本基因组进行过鉴定分析[13]，但由于早期版本基因组拼接和注释的质量有限，且无法将基因定位于染色体，目前研究人员普遍使用的是9930_V2版本基因组[14]。所以，基于更通用版本基因组对ERF基因家族进行鉴定，将会提供更全面和准确的鉴定结果。【拟解决的关键问题】通过基于9930_V2版本基因组对ERF基因家族进行鉴定与表达分析，明确ERF在雌花分化与发育阶段的特异表达情况，为解释这一生物学过程的分子机制奠定理论基础。 1 材料与方法试验于2017年春季在上海交通大学农业与生物学院进行。 1.1 黄瓜材料本研究试验材料为华北类型黄瓜近等基因系：全雌株系419-19（FFMMAA基因型）和主茎强雄株系419-2（ffMMAA）。材料种植培养环境为28℃ 18 h昼，18℃ 6 h夜。顶芽选取四叶一心时期，不同时期雌花芽选取四叶一心时期顶芽内不同大小的雌花芽。所取材料采用3个独立生物学重复，液氮速冻后，置于-80℃超低温冰箱备用。 应用型本科教育对于满足中国经济社会发展，对高层次应用型人才需要以及推进中国高等教育大众化进程起到了积极的促进作用。最近几年，随着外资企业的不断引进和中国社会经济的发展，对商务英语人才的需求越来越多。商务英语人才质量的高低对于外资企业的各种业务开展起到举足轻重的作用。当前，作为培养商务英语人才的应用型本科院校商务英语教师面临着专业化发展上的各种障碍。很明显，社会对商务英语人才的迫切需求与商务英语教师专业化发展困难之间出现了不均衡的矛盾。因此，加强应用型本科院校商务英语教师专业化发展，是应用型本科院校目前师资队伍建设的一个重要任务。 1.2 黄瓜ERF基因家族成员鉴定根据文献和TAIR数据库获得拟南芥（Arabidopsis thaliana）ERF基因家族的ID和氨基酸序列信息。黄瓜基因组数据信息均获取自CuGenDB（http:// cucurbitgenomics.org）。采用双向BLAST方法获取黄瓜基因组中ERF基因ID。双向BLAST方法简述如下：采用Geneious软件将拟南芥ERF基因家族成员比对至黄瓜9930_V2版本基因组，获取基因ID，再将所获取基因ID比对至Swiss-Port数据库，确认获取ID为ERF家族成员。通过Geneious软件进行序列多重比对，鉴定保守结构域。同时，通过NCBI网站的CDD工具预测ERF基因的AP2/ERF结构域，确定每个基因均包含该结构域。采用TBtools软件[15]绘制黄瓜ERF基因家族成员在染色体上的物理位置。 1.3 黄瓜ERF基因家族进化、基因结构和GO富集分析利用MEGA7软件对黄瓜ERF基因家族成员进行系统进化树分析，分析序列为氨基酸，采用Neighbor Joining（NJ）法进行建树，Bootstrap值设置为1 000。采用EMBOSS软件进行ERF基因家族成员的信号肽和跨膜域预测。通过MEME（ http://meme-suite.org）网站工具进行motif预测分析。将Motif预测、CDD结构域鉴定以及进化树分析结果，通过TBtools软件整合并制图。VENN交集分析采用VENNY2.1[16]，GO富集分析在葫芦科基因组网站（http:// cucurbitgenomics. org）进行。 将植物材料研磨，采用康为世纪公司试剂盒RNA抽提、反转录获取cDNA。qRT-PCR试验采用罗氏公司FastStart Essential DNA Green Maste试剂，程序参见试剂盒说明书。以 CsActin（ Csa6M484600.1）作为内参，所用引物均查询自qPrimerDB数据库（ https://biodb.swu.edu.cn/qprimerdb/）。采用3个生物学重复与3个技术重复。试验数据用Excel软件进行处理，采用2-ΔΔCT方法进行相对定量分析，显著性分析采用 t-test公式计算。在表达趋势分析中，将不同时期花芽中表达的CsERF基因家族成员的qRT-PCR相对定量结果导入Omicshare趋势分析工具（ www.omicshare.com/tools），采用 t-test 检验两个时期之间表达差异大于2倍（ P value＜0.05）的基因，定义为表达存在显著差异的基因。 在生工生物公司合成包含3个串联重复GCC-box的启动子DNA序列“TCCAAGGAGCCGCCGGG”（来源于CsERF39的启动子区，下划线为酶切位点），将合成序列通过诺唯赞ClonExpress II重组试剂盒与pLacZ载体融合，构建为pLacZ-3xGCC-box载体。采用TaKaRa的PrimerStar GXL DNA聚合酶克隆CsERF9和CsERF31的开放和阅读框，将其与pB42AD载体融合，构建为pB42AD-CsERF9和pB42AD- CsERF31载体。酵母单杂交分析中，采用EGY48a菌株，以空载体为负对照，在SD/-Trp/-Ura缺陷型培养基进行筛选培养，并在加有X-gal的SD/-Trp/-Ura缺陷型定性培养基上进行互作鉴定[17]。 人工智能取代劳动者导致大规模失业，这是人工智能在劳动领域里引发的最大忧虑，有关人工智能影响的论述无不提及这一点。《新一代人工智能发展规划》中也明确指出，人工智能的发展可能带来就业结构的改变。机器取代人，这是人工智能进化发展的宿命，抑或只是理论上的臆想？对可能造成的劳动者失业劳动法上应如何应对？ 2 结果2.1 黄瓜ERF基因家族成员的鉴定与相关序列信息通过双向BLAST分析，以拟南芥AP2/EREBP基因家族为参考，对黄瓜基因组中ERF基因家族进行鉴定，共鉴定出138个具有AP2/ERF结构域的基因（图1），根据其在9930_V2基因组中的ID，分别命名为CsERF1—CsERF138。长度介于126—745个氨基酸。通过BLAST比对，将9930_V2基因组ERF家族成员比对到拟南芥TAIR10基因组中，获取对应基因名称（表1）。相比根据9930_V1版本所鉴定的103个ERF基因[13]，本研究新鉴定出25个ERF家族成员，并将全部138个ERF定位于染色体。 基于扩展序列运算的含不确定性需求响应电力系统优化调度//徐箭,曹慧秋,唐程辉,魏聪颖,江海燕,廖思阳//(13):152 信道建模方法主要包括统计信道建模、确定性信道建模和半确定性信道建模.统计性建模主要依赖于信道测量,文献[2]对高速铁路进行了无线信道测量,但基于实测的建模灵活性及移植性较差.确定性建模如射线跟踪法可对地铁隧道内的多径信道建模[3],但算法复杂、消耗资源较多.二者结合的半确定性信道建模[4]能更好地符合实际环境,其中使用较多的是几何建模方法.  图1 黄瓜ERF家族成员AP2/ERF保守结构域序列 Fig. 1 The AP2/ERF conserved domain in ERF members of cucumber 通过EMSS软件对138个CsERF家族成员进行信号肽和跨膜域预测分析，结果显示有40个CsERF蛋白预测含有信号肽（表1中标注*），其中CsERF55和CsERF116预测含有跨膜结构域（表1中标注#）。 将138个CsERF成员根据物理位置定位到7条染色体上，结果如图2所示。其中，CsaUNG003730和CsaUNG031640未定位于7条染色体中，通过两基因比对至9930_V3版本基因组，发现其定位于5号染色体。同时，通过将138个CsERF家族成员的CDS序列与转录组拼接结果比较，发现Csa2G001000和Csa2G307850两个基因的转录本存在注释错误，依照转录组数据修正后的两基因用于后期生物信息学分析。 2.2 黄瓜ERF基因家族进化分析为了解CsERF家族成员的进化特性，对138个基因构建系统进化树（图3）。CsERF家族成员共分为10个亚家族（Class I—X），其中成员最多的为Class IX，包含37个基因。成员最少的是Class V，包含4个具有B3结构域的基因。 2.3 黄瓜ERF基因家族基因结构分析CsERF基因家族的结构分析结果如图4所示，共鉴定出7个motif，其中motif 1在所有家族中均出现，而在除了Class VI外的其他家族中，大部分成员均包含按顺序排列的motif 2-3-1-4。在Class VI中，大部分成员具有按相同顺序排列的motif 2-3/6-1-4-5-2-7- 1-4，说明这一结构域在该亚家族中较为保守。在Class IX中，具有一个特殊的motif 8与motif 4相邻。综合motif分析结果，虽然不同亚族中包含motif的种类以及排列顺序不尽相同，但在同一亚族中，均包含相同或相似的motif排序，这可能与该亚家族的分子功能相关。 表1 CsERF基因家族相关信息 Table 1 Information of CsERF gene family members 
 名称Name黄瓜基因组基因编号ID in cucumber (9930_V2)注释到拟南芥成员Annotated to Arabidopsis成员名称Name in Arabidopsis所属亚族Class氨基酸个数Number of amino acids CsERF1*Csa1G000550AT5G61890ERF1142216 CsERF2Csa1G042290AT1G68550ERF1182336 CsERF3Csa1G075030AT1G28360ERF124164 CsERF4Csa1G075060AT5G44210ERF94211 CsERF5Csa1G269870AT1G50680-5339 CsERF6Csa1G275920AT1G16060ADAP6363 CsERF7Csa1G340430AT1G15360WIN17223 CsERF8Csa1G364460AT5G11590DREB39250 CsERF9Csa2G177210AT3G23240ERF1B1228 CsERF10Csa1G597730AT1G19210ERF0179203 CsERF11*Csa1G605670AT2G05990MOD15745 CsERF12*Csa2G001000AT5G11590DREB39252 CsERF13Csa2G006270AT1G15360WIN17246 CsERF14*Csa2G092800AT5G17430BBM6650 CsERF15Csa2G138780AT4G13620ERF06210400 CsERF16Csa1G423190AT5G57390AIL56466 CsERF17Csa2G177220AT3G23230ERF0981126 CsERF18Csa2G279250AT4G36920AP26483 CsERF19*Csa2G279260AT2G23340ERF0089159 CsERF20*Csa2G297760AT2G33710ERF1122285 CsERF21*Csa2G307850AT1G78080RAP2-410362 CsERF22Csa2G349090AT1G12980ESR14358 CsERF23*Csa2G354000AT4G37750ANT6696 CsERF24Csa2G356600AT1G46768RAP2-19142 CsERF25*Csa2G363010AT1G75490DREB2D8212 CsERF26Csa2G374590AT5G11590DREB39231 CsERF27Csa2G382540AT5G25190ERF0037175 CsERF28*Csa2G382550AT5G25190ERF0037164 CsERF29*Csa2G401240AT1G22810ERF0199148 CsERF30*Csa3G017320AT4G17500ERF1A1279 CsERF31Csa3G135120AT3G23240ERF1B1228 CsERF32Csa3G018320AT4G17490ERF61344 CsERF33Csa3G019400AT2G44940ERF0349214 CsERF34Csa3G042390AT5G52020ERF0259170 CsERF35Csa3G073900AT3G20310ERF74214 CsERF36*Csa3G114470AT5G65510AIL76283 CsERF37Csa3G116720AT2G44840ERF131144 CsERF38Csa3G116730AT2G44840ERF131157 CsERF39Csa3G389850AT3G23240ERF1B1224 CsERF40*Csa3G120390AT4G27950CRF43283 续表1 Continued table 1 名称Name黄瓜基因组基因编号ID in cucumber (9930_V2)注释到拟南芥成员Annotated to Arabidopsis成员名称Name in Arabidopsis所属亚族Class氨基酸个数Number of amino acids CsERF41Csa3G124760AT5G19790RAP2-113243 CsERF42Csa3G133130AT1G68550ERF1182213 CsERF43Csa3G135620AT3G23230ERF0981158 CsERF44Csa3G135630AT3G23230ERF0981133 CsERF45Csa3G152130AT5G18560ERF0864326 CsERF46*Csa3G164580AT3G16770RAP2-32225 CsERF47Csa3G180260AT4G25470DREB1C9213 CsERF48*Csa3G207390AT4G34410ERF1092274 CsERF49Csa3G357110AT4G27950CRF43326 CsERF50*Csa3G118010AT2G44940ERF0349230 CsERF51Csa3G630280AT4G13040-10234 CsERF52Csa3G646560AT3G14230RAP2-22370 CsERF53*Csa3G652380AT1G72570AIL16490 CsERF54Csa3G736760AT2G28550RAP2-76527 CsERF55*#Csa3G751430AT2G36450ERF0249183 CsERF56*Csa3G751440AT4G25490DREB1B9204 CsERF57Csa3G812170AT3G57600DREB2F8333 CsERF58*Csa3G822440AT2G41710-6432 CsERF59Csa3G827310AT3G54320WRI16337 CsERF60Csa3G865380AT1G75490DREB2D8184 CsERF61Csa3G878210AT5G25390SHN37194 CsERF62Csa3G895680AT4G32800ERF0439227 CsERF63Csa4G001970AT3G16770RAP2-32272 CsERF64Csa4G004960AT1G15360WIN17190 CsERF65Csa4G006210AT5G25190ERF0037180 CsERF66Csa4G007070AT1G28160ERF0874227 CsERF67Csa4G007650AT1G71450ERF0219179 CsERF68Csa4G023020AT5G61600ERF1041194 CsERF69Csa4G051360AT4G23750CRF23319 CsERF70Csa4G192030AT5G19790RAP2-113225 CsERF71Csa4G268100AT5G19790RAP2-113211 CsERF72*Csa4G290800AT1G51190PLT26615 CsERF73Csa4G292470AT2G28550RAP2-76456 CsERF74Csa4G314390AT2G40340DREB2C8372 CsERF75Csa4G370550AT1G21910ERF0129232 CsERF76Csa4G630010AT3G15210ERF44193 CsERF77Csa4G641590AT3G50260ERF0119149 CsERF78*Csa4G644740AT4G37750ANT6573 CsERF79Csa4G649620AT1G74930ERF0189186 CsERF80Csa4G652640AT5G25190ERF0037203 CsERF81Csa5G139630AT4G27950CRF43268 续表1 Continued table 1 名称Name黄瓜基因组基因编号ID in cucumber (9930_V2)注释到拟南芥成员Annotated to Arabidopsis成员名称Name in Arabidopsis所属亚族Class氨基酸个数Number of amino acids CsERF82Csa5G146300AT5G19790RAP2-113210 CsERF83Csa5G150420AT1G64380ERF06110304 CsERF84Csa5G151530AT4G23750CRF23318 CsERF85Csa5G155560AT5G52020ERF0259252 CsERF86*Csa5G155570AT4G25470DREB1C9225 CsERF87Csa5G165850AT4G17500ERF1A1250 CsERF88Csa5G167110AT5G07580ERF1061160 CsERF89*Csa5G167120AT5G51190ERF1051235 CsERF90Csa5G174570AT5G51990DREB1D9200 CsERF91Csa5G175970AT2G28550RAP2-76441 CsERF92Csa5G598600AT1G24590ESR24341 CsERF93Csa5G608380AT1G68840RAV25344 CsERF94Csa5G609620AT1G68550ERF1182309 CsERF95*Csa5G612310AT3G54320WRI16439 CsERF96*Csa5G637750AT5G19790RAP2-113395 CsERF97Csa5G647260AT5G25810TINY9163 CsERF98Csa5G649870AT1G33760ERF0229163 CsERF99Csa5G649890AT1G71450ERF0219192 CsERF100*Csa6G011730AT2G40220ABI48305 CsERF101Csa6G012810AT2G40340DREB2C8410 CsERF102Csa6G167230AT1G50640ERF34219 CsERF103*Csa6G040610AT1G21910ERF0129201 CsERF104Csa6G042450AT2G46310CRF53243 CsERF105*Csa6G055940AT5G61890ERF1142404 CsERF106Csa6G091830AT5G13330ERF1132191 CsERF107Csa6G104640AT4G32800ERF0439213 CsERF108*Csa6G124180AT2G40340DREB2C8179 CsERF109Csa6G133770AT4G27950CRF43281 CsERF110*Csa6G017030AT5G61890ERF1142310 CsERF111*Csa6G296960AT2G28550RAP2-76497 CsERF112Csa6G318160AT3G16770RAP2-32231 CsERF113Csa6G361330AT1G19210ERF0179213 CsERF114Csa6G404260AT3G54320WRI16366 CsERF115Csa6G421660AT1G50640ERF34214 CsERF116#*Csa6G450420AT2G36450ERF0249224 CsERF117*Csa6G486790AT5G61890ERF1142338 CsERF118Csa6G490860AT4G34410ERF1092193 CsERF119Csa6G491020AT4G36920AP26537 CsERF120*Csa6G491030AT5G67190ERF0109173 CsERF121*Csa6G496390AT4G37750ANT6629 CsERF122Csa6G500550AT1G78080RAP2-410271 续表1 Continued table 1 名称Name黄瓜基因组基因编号ID in cucumber (9930_V2)注释到拟南芥成员Annotated to Arabidopsis成员名称Name in Arabidopsis所属亚族Class氨基酸个数Number of amino acids CsERF123Csa6G518040AT3G15210ERF44246 CsERF124Csa6G518290AT5G47220ERF21177 CsERF125Csa6G518300AT4G18450ERF0911277 CsERF126Csa7G047400AT2G20880ERF05310397 CsERF127Csa7G049230AT5G44210ERF94209 CsERF128*Csa7G073700AT4G34410ERF1092256 CsERF129Csa7G352440AT4G13040-10231 CsERF130Csa7G375820AT3G23240ERF1B1231 CsERF131Csa7G375830AT3G23230ERF0981148 CsERF132*Csa7G431330AT1G16060ADAP6367 CsERF133Csa7G432080AT5G25190ERF0037196 CsERF134Csa7G432130AT1G80580ERF0843183 CsERF135Csa7G447150AT5G11590DREB39201 CsERF136*Csa7G448110AT1G53910RAP2-122389 CsERF137CsaUNG003730AT5G13910LEP4228 CsERF138*CsaUNG031640AT1G13260RAV15317
*：有信号肽；#：有跨膜结构域 *: Signal peptide; #: Transmembrane domain  图2 黄瓜ERF基因家族在染色体上的位置 Fig. 2 The chromosome location of the ERF gene family in cucumber  图3 黄瓜ERF家族进化树（罗马数字I-X代表10个亚家族） Fig. 3 The phylogenetic tree of ERF family in cucumber (Roman numerals I-X indicate 10 sub-families) 结构域分析表明，所有CsERF家族成员均具有保守的AP2/ERF结构域，其中4个基因具有B3结构域。在CsERF家族成员中，多数亚家族成员无内含子结构，但Class II、Class VI和Class VII亚家族中的大部分成员，基因内具有多个内含子结构。其中Class VI亚家族较为特殊，其成员均包含2个AP2/ERF结构域，且多数内含子分布在AP2/ERF结构域中，推测该亚家族在进化过程中AP2/ERF结构域发生了复制，造成这一结构域被内含子间隔的特殊现象。 970 智能预警系统在脑卒中患者院前院内衔接中的应用 张洪剑，刘团结，王文安，徐建华，解 炯，蒋 超，徐中杰，张永巍，杨鹏飞，邓本强，朱勤忠，刘建民 2.4 黄瓜ERF在不同性型材料中的特异表达分析ERF基因家族在多种植物中被报道参与乙烯信号传导过程，在黄瓜花芽性别分化过程中，乙烯信号是决定雌花分化的关键。为鉴定黄瓜ERF家族成员与性别决定与雌花分化过程的关系，通过RT-qPCR方法对全雌株系419-19（FFMMAA）和419-2主茎强雄株系M12（ffMMAA）植株顶芽转录本进行检测。表达存在显著差异（p value＜0.05）的基因如图5所示，CsERF9、CsERF30、CsERF31、CsERF37、CsERF87、CsERF98、CsERF105、CsERF106和CsERF110这9个基因在全雌株系中高表达（图5-A），CsERF10、CsERF48、CsERF68、CsERF69、CsERF71、CsERF76、CsERF81、CsERF89、CsERF127和CsERF128这10个基因在强雄株系中高表达（图5-B）。在ERF基因家族中，上述差异表达基因主要集中于Class I和Class II亚家族，且具有成簇分布的趋势，如CsERF9和CsERF31相邻并在全雌株系中高表达，CsERF30、CsERF87和CsERF37也在全雌株系中高表达。 红色字体代表在FFMMAA基因型材料顶芽中高表达，蓝色字体代表在ffMMAA基因型材料顶芽中高表达，红色三角号代表雌花芽发育初期表达呈上调趋势，蓝色三角号代表雌花芽发育初期表达呈下调趋势 Genes in red indicate highly express in FFMMAA genotype line, genes in blue indicate highly express in ffMMAA genotype line, red and blue triangles indicate genes show an up- and down-regulated trend in female bud development, respectively 图4 黄瓜ERF基因家族的基因结构分析 Fig. 4 Gene structure analysis of ERF family members in cucumber 2.5 黄瓜ERF基因家族成员在雌花发育初期的表达模式分析有研究表明，乙烯信号在促进雌花分化的同时，对果实发育也存在调控作用，如CsSF1通过结合M基因（CsACS2）介导其泛素化，进而精细调控M基因的表达，控制果实的发育[18]。M基因的特异表达时期是雌花分化的1—10 mm时期[19]，所以本研究分别选取这一时期2、5和10 mm的雌花芽，对ERF基因家族成员和部分ACO和ACS的表达变化趋势进行检测，其中5个ACO命名与前人一致[20]。将不同时期间存在显著表达变化（P＜0.05）且差异大于2倍以上的基因，定义为表达存在显著差异的基因。结果如图6所示，共有61个ERF在雌花芽发育初期表达趋势发生变化，其中29个ERF在2—5 mm阶段表达上调，而5—10 mm阶段表达无显著变化，2个ERF在2—10 mm阶段均呈上升趋势。同时，25个ERF在2—5 mm阶段表达下调，而5—10 mm阶段表达无显著变化，5个ERF在2—10 mm阶段均呈下降趋势。此外，A和CsACO2在5—10 mm时期特异高表达，而M和CsACO1在2—5 mm阶段特异高表达。 值得注意的是，ERF基因家族在雌花芽发育初期的表达趋势具有明显的结构特异性。Class I和Class II亚家族中，表达变化的ERF均呈现下降趋势。同时，Class VI和Class VII中，表达变化的ERF多数呈上升趋势。此外，一些在系统进化树中相邻的基因，也表现出相似的变化趋势，如Class VI中的CsERF16、CsERF36、CsERF53、CsERF23、CsERF78和CsERF121，以及Class IX中的CsERF33、CsERF107、CsERF135和CsERF50。这说明一些进化上趋于保守的基因家族成员具有相似的表达模式，在雌花芽发育初期所发挥的功能可能相似。  图5 在FFMMAA（A）和ffMMAA（B）基因型材料顶芽中显著高表达的CsERF家族成员 Fig. 5 The significant differential expression CsERF in FFMMAA (A) and ffMMAA (B) cucumber lines shoot apical (t test, P value＜0.05)  图6 CsERF家族成员在雌花芽发育初期表达存在显著变化的成员 Fig. 6 The significant differences in the expression of CsERF in early stages of female bud development 很多ERF转录因子通过与其下游应答元件GCC-box的相互作用来调控基因表达。通过对黄瓜基因组注释基因启动子区进行序列分析，共检测到882个基因启动子具有GCC-box（AGCCGCC）。为验证黄瓜ERF是否具有保守的结合GCC-box功能，选取拟南芥ERF1同源基因CsERF9和CsERF31，通过酵母单杂交系统验证其与GCC-box的结合作用，发现两基因均可以直接结合GCC-box（图7）。这说明部分黄瓜ERF家族成员具有保守结合GCC-box元件的功能，同时也为从转录组数据中高效筛选下游直接调控靶基因提供了分子依据。 为进一步研究雌花芽发育初期受ERF转录因子调控的基因表达变化，从国家基因库生物大数据平台（CNGBdb）的CNP0000657项目中提取2、5和10 mm雌花芽转录组数据（CNS0110220、CNS0110221、CNS0110222、CNS0110223、CNS0110224、CNS0110225）进行表达趋势分析。共筛选出2 969个基因与CsERF9和CsERF31表达趋势相似，随着雌花芽分化表达呈下降趋势，同时筛选出7 051个基因的表达随雌花芽分化呈上升趋势。通过将这些基因与882个启动子具有GCC-box的基因进行VENN交集分析，发现有148个基因的表达随着雌花芽分化下调，而232个基因上调，且这些基因启动子均具有GCC-box（图7-B）。将上调和下调基因进行分子功能（molecular function）GO富集分析，发现上调基因中富集最显著的是GO: 0005667途径，共有10个基因富集其中，注释与转录复合体（transcription factor complex）相关。在下调基因中富集最显著的是GO:0003723途径，共有14个基因富集其中，注释与RNA结合（RNA binding）相关。 山东地区气候温和，属于暖温带季风气候，四季分明，年平均降水量为600～750 mm，雨量集中。每年7—8月份葡萄果实成熟，却大量降雨，造成葡萄病害发生严重，裂果、烂果率高[5]。笔者对不同模式避雨棚搭建成本进行分析，筛选出低成本、寿命长的避雨棚搭建模式，以期为葡萄避雨栽培生产提供理论依据。  图7 酵母单杂交系统验证CsERF31和CsERF9对GCC-box的互作（A），雌花芽发育初期表达成下降趋势的基因与启动子具有GCC-box的基因VENN分析（B），交集中表达下调（左）和上调（右）基因的GO富集分析（C） Fig. 7 Yeast one-hybrid assay of interaction between CsERF31, CsERF39 and GCC-box (A), the two sets of down-regulated and GCC-box in promoter genes analyzed using the VENN method (B), and GO enrichment analysis of down-regulated (left) and up-regulated (right) genes with GCC-box in promoters (C) 自1994年在拟南芥中分离出第一个与花发育相关基因AP2[21]以来，已报道大量对AP2/ERF转录因子的基因家族分析，包括玉米[22]、粟[23]、葡萄[24]和橙[25]等。本研究主要关注黄瓜ERF基因家族在黄瓜雌花芽的分化与发育过程中的表达趋势变化，进而为深入研究这一重要生物学过程的分子机制提供数据基础。根据前人对黄瓜花芽性别分化时期的划分，6—7时期（Stage 6—7）是花芽性别分化时期[26]。在6—7时期，花芽长度小于1 mm，而在2 mm阶段，是花芽由性别分化向雌器官起始发育的转变期。2 mm阶段的雌花芽，子房已经起始发育，而在5—10 mm阶段，已发育出具有明显形态特征的子房，各雌器官均已分化完成。在ERF基因家族中，大部分基因的表达变化集中于雌花芽2—5 mm阶段，说明这一阶段可能发生了复杂的基因表达调控，使雌性器官得以正确的分化。 3.2 Class I和Class II类基因可能在雌花的分化与发育过程中发挥功能在黄瓜ERF基因家族的10个亚家族中，Class I和Class II成员的表达模式较为特殊。值得注意的是在雌花芽分化初期，Class I和Class II中所有表达发生变化的基因，均呈现下降趋势（图4中蓝色箭头标注基因）。表达量呈下降趋势暗示这些基因可能在花芽发育最初阶段发挥功能，随着子房发育，其表达被抑制而趋于稳定。 此外，选取不同性别类型植株和不同雌花芽发育阶段进行表达检测，是将雌花的性别分化与雌性器官发育假定为两个生物学过程。因此，图4中“红色字体，蓝色三角号”基因（即在全雌株顶芽高表达，并随雌花芽发育表达下降）推测与性别分化更相关，“蓝色字体，红色三角号”基因（即在强雄株顶芽高表达，并随雌花芽发育表达上升）则可能与雌花芽发育更相关，“蓝色字体，蓝色三角号”基因（即在强雄株顶芽高表达，并随雌花芽发育表达下降）更有可能与雄花分化相关。值得注意的是，在所有具有表达变化的ERF中，未发现在全雌株系中高表达，且随着雌花芽分化表达上调的基因（即图4中“红色字体，红色三角号”）。据此推测ERF参与性别分化与雌花芽发育的过程受到严格调控，在不同性别花芽、不同花芽分化阶段发挥不同的功能。 3.3 部分黄瓜ERF的调控机制与模式植物相似在黄瓜ERF基因家族中，CsERF9、CsERF31、CsERF39以及CsERF130在拟南芥中的直系同源基因均为ERF1B。拟南芥ERF1B在乙烯信号传导通路中，受上游EIN3/EIL1等基因调控，通过激活下游基因参与多种生物学过程[27]。在上述4个基因中，CsERF31已被报道具有正调控M基因（CsACS2）的功能，这说明其很可能通过乙烯信号传导途径参与花芽的性别分化过程。此外，CsERF105和CsERF110在拟南芥中的同源基因均为ERF114，且两基因具有与CsERF31和CsERF9相似的表达趋势，其中CsERF110已被报道具有正调控A基因（CsACS11）的功能[10]，而CsERF31被报道可能参与M基因自反馈调控[28-30]。有报道显示，F（CsACS1G）、M（CsACS2）和A（CsACS11）的不同基因型材料，其乙烯释放速率存在显著差异[31]，上述结果暗示，这类基因很可能通过乙烯信号传导途径介导黄瓜雌花的分化。 在Class VI中，CsERF23、CsERF78和CsERF121在拟南芥中的同源基因均为ANT。ANT在拟南芥中调控花芽等侧生器官发育，其突变体在育性、子房发育和侧生器官分化方面均存在缺陷[32]。本研究中，CsERF23、CsERF78和CsERF121随雌花芽发育表达上调，这暗示其可能具有调控雌性器官发育的功能。CsERF69和CsERF84在拟南芥中的同源基因为CRF2。CRF2在拟南芥中受细胞分裂素诱导，转移至细胞核调控下游基因[33]。CsERF69和CsERF84的表达随雌花芽分化上调，推测其同样受到细胞分裂素等激素的调控，在子房发育过程中行使功能。 此外，拟南芥ERF1同源蛋白CsERF9和CsERF31具有结合乙烯响应元件GCC-box的功能，说明黄瓜ERF家族基因可能存在与模式植物相似的分子调控方式影响下游基因的表达。通过进一步的生物信息学分析，鉴定出大量在雌花芽发育初期表达变化的基因，对启动子具有GCC-box的基因进行GO聚类分析，发现富集最显著的是转录因子复合体和RNA结合两个途径。转录因子复合体和RNA结合蛋白均为调控RNA转录或转录后翻译的影响因子，其在雌花芽发育初期表达的集中变化，说明这一过程可能涉及转录激活调控、转录后调控以及翻译后修饰等一系列调控过程，近期报道的CsSF1介导M（CsACS2）蛋白泛素化调控子房发育也为这一过程提供了新的分子证据[18]。综上所述，推测ERF在雌花芽发育初期的表达变化，直接或间接影响了下游基因的表达，使雌花芽中各器官得以正确分化，进而完成果实的发育。 （141）尖瓣扁萼苔 Radula apiculata Sande Lac.ex.Steph.余夏君等（2018） 4 结论本研究从黄瓜9930_V2基因组中鉴定出138个具有AP2/ERF结构域的ERF基因家族成员，可分为10个亚族。在花芽性别分化和雌花芽发育不同阶段中，ERF具有明显的表达变化，并且其表达变化与基因结构存在相关性。一些ERF可以结合保守元件GCC-box，推测其通过结合GCC-box来调控下游靶基因的表达。 References [1] UN Food and Agriculture Organization, Corporate Statistical Database (FAOSTAT). Production of cucumbers and gherkins in 2016[OL], 2017. [2] Malepszy S, Niemirowicz-Szczytt K. Sex determination in cucumber (Cucumis sativus) as a model system for molecular biology. Plant Science, 1991, 80(1/2): 39-47. [3] Zhang Z J, Zhang H W, Quan R D, Wang X C, Huang R F. Transcriptional regulation of the ethylene response factor LeERF2 in the expression of ethylene biosynthesis genes controls ethylene production in tomato and tobacco. Plant Physiology, 2009, 150(1): 365-377. [4] Xiao Y y, Chen J y, Kuang J f, Shan W, Xie H, Jiang Y m, Lu W J. Banana ethylene response factors are involved in fruit ripening through their interactions with ethylene biosynthesis genes. Journal of Experimental Botany, 2013, 64(8): 2499-2510. [5] Li T, Jiang Z Y, Zhang L C, Tan D M, Wei Y, Yuan H, Li T L, Wang A D. Apple (Malus domestica) MdERF2 negatively affects ethylene biosynthesis during fruit ripening by suppressing MdACS1 transcription. The Plant Journal, 2016, 88(5): 735-748. [6] Zhu X L, Qi L, Liu X, Cai S B, Xu H J, Huang R F, Li J R, Wei X N, Z, Zhang Z Y. The wheat ethylene response factor transcription factor PATHOGEN-INDUCED ERF1 mediates host responses to both the necrotrophic pathogen rhizoctoniacerealis and freezing stresses. Plant physiology, 2014, 164(3): 1499. [7] Ohme-Takagi M, Shinshi H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. ThePlant Cell, 1995, 7: 173-182. [8] Lorenzo O, Piqueras R, Sánchez-Serrano J J, Solano R. ETHYLENE RESPONSE FACTOR1 integrates signals from ethylene and jasmonate pathways in plant defense. The Plant Cell, 2003, 15: 165-178. [9] Mao J L, Miao Z Q, Wang Z, Yu L H, Cai X T, Xiang C B. Arabidopsis ERF1 mediates cross-Talk between ethylene and auxin biosynthesis during primary root elongation by regulating ASA1 expression. PLoS Genet, 2016, 12: e1006076. [10] Zhong S W, Shi H, Xue C, Wei N, Guo H W, Deng X W. Ethylene-orchestrated circuitry coordinates a seedling’s response to soil cover and etiolated growth. Proceeding of the National Academy of Sciences of the United States of America, 2014, 111: 3913-3920. [11] Tao Q Y, Niu H H, Wang Z Y, Zhang W H, Wang H, Wang S H, Zhang X, Li Z. Ethylene responsive factor ERF110 mediates ethylene-regulated transcription of a sex determination-related orthologous gene in two Cucumis species. Journal of Experimental Botany, 2018, 69: 2953-2965. [12] Pan J, Wang G, Wen H F, Du H, Lian H L, He H, Pan J S, Cai R. Differential gene expression caused by the F and M loci provides insight into ethylene-mediated female flower differentiation in cucumber. Fronts in Plant Science, 2018, 9: 1091. [13] Hu L f, Liu S Q. Genome-wide identification and phylogenetic analysis of the ERF gene family in cucumbers. Genetic Molecular Biology, 2011, 34: 624-633. [14] Huang S W, Li R q, Zhang Z H, Li L, Gu X f, Fan W, Lucas W J, Wang X w, Xie B y, Ni P X, et al. The genome of the cucumber, Cucumis sativus L. Nature Genetics, 2009, 41: 1275-1281. [15] Chen C j, Xia R, Chen H, He Y h. TBtools, a Toolkit for Biologists integrating various HTS-data handling tools with a user-friendly interface. bioRxiv, 2018, 3: 289660. doi: https://doi.org/ 10.1101/289660. [17] Yan T x, Chen M h, Shen Q, Li L, Fu X q, Pan Q F, TANG Y L, SHI P, LV Z Y, JIANG W M, MA Y N, HAO X L, SUN X F, Tang K X. HOMEODOMAIN PROTEIN 1 is required for jasmonate- mediated glandular trichome initiation in Artemisia annua. New Phytologist, 2017, 213: 1145-1155. [18] Xin T x, Zhang Z, Li S, Zhang S, Li Q, Zhang Z h, Huang S W, Yang X Y. Genetic regulation of ethylene dosage for cucumber fruit elongation. The Plant Cell, 2019.doi: https://doi.org/10.1105/ tpc.18.00957. [19] Saito S, Fujii N, Miyazawa Y, Yamasaki S, Matsuura S, Mizusawa H, Fujita Y, Takahashi H. Correlation between development of female flower buds and expression of the CS-ACS2 gene in cucumber plants. Journal of Experimental Botany, 2007, 58(11): 2897-2907. [20] Chen H m, Sun J j, Li S, Cui Q z, Zhang H m, Xin F j, Wang H s, Lin T, Gao D l, Wang S H, LI X, WANG D H, ZHANG Z H, XU Z H, HUANG S W. An ACC oxidase gene essential for cucumber carpel development. Molecular Plant, 2016, 9: 1315-1327. [21] Jofuku K D, DEN Boer B G w, van montagu m, Okamuro j k. Control of Arabidopsis flower and seed development by the homeotic geneapetala2. The Plant Cell, 1994, 6(9): 1211-1225. [22] Zhou M L, Tang Y X, Wu Y M. Genome-wide analysis of AP2/ERF transcription factor family in zea mays. Current Bioinformatics, 2012, 7(3): 324-332. [23] Lata C, Mishra A K, Muthamilarasan M, Bonthala V S, Khan Y, Prasad M. Genome-wide investigation and expression profiling of AP2/ERF transcription factor superfamily in foxtail millet (Setaria italica L.). PLoS ONE, 2014, 9: e113092. [24] Licausi F, Giorgi F M, Zenoni S, Osti F, Pezzotti M, Perata P. Genomic and transcriptomic analysis of the AP2/ERF superfamily in Vitis vinifera. BMC Genomics, 2010, 11: 719. [25] Ito T M, Polido P B, Rampim M C, Kaschuk G, Souza S G H. Genome-wide identification and phylogenetic analysis of the AP2/ERF gene superfamily in sweet orange (Citrus sinensis). Genetic Molecular Research, 2014, 13: 7839-7851. [26] Bai S L, Peng Y B, Cui J X, Gu H T, Xu L Y, Li Y Q, Xu Z H, Bai S N. Developmental analyses reveal early arrests of the spore-bearing parts of reproductive organs in unisexual flowers of cucumber (Cucumis sativus L.). Planta, 2004, 220: 230-240. [27] 张存立, 郭红卫. 乙烯信号转导通路研究. 自然杂志, 2012, 34(4): 219-228. Zhang C L, Guo H W. Study on ethylene signal transduction pathway. China Journal of Nature, 2012, 34(4): 219-228. (in Chinese) [28] Li Z, Huang S W, Liu S q, Pan J S, Zhang Z H, Tao Q Y, Shi Q x, Jia Z q, Zhang W w, Chen H m, si l t, zhu l h, Cai R. Molecular isolation of the M gene suggests that a conserved-residue conversion induces the formation of bisexual flowers in cucumber plants. Genetics, 2009, 182: 1381-1385. [29] Li Z, Pan J s, Guan Y, Tao Q Y, He H L, Si L T, Cai R. Development and fine mapping of three co-dominant SCAR markers linked to the M/m gene in the cucumber plant (Cucumis sativus L.). Theoretical and Applied Genetics, 2008 117: 1253-1260. [30] Li Z, Wang S, Tao Q T, Pan J S, Si L T, Gong Z h, Cai R. A putative positive feedback regulation mechanism in CsACS2 expression suggests a modified model for sex determination in cucumber (Cucumis sativus L.). Journal of Experimental Botany, 2012, 63: 4475-4484. [31] 崔清志, 陈宸, 田云, 刘晓虹, 陈惠明. 不同基因型黄瓜性别与乙烯释放速率的关系. 中国蔬菜, 2016(3): 37-42. Cui Q Z, Chen C, Tian Y. Liu X H, Chen H M. Relationship between cucumber sex of different genotypes and ethylene release rate. China Vegetables, 2016(3): 37-42. [32] Krizek B A, Lewis M W, Fletcher J C. RABBIT EARS is a second-whorl repressor of AGAMOUS that maintains spatial boundaries in Arabidopsis flowers. Plant Journal, 2006, 45(3): 369-383. [33] Jeon J. Cho C, Lee M R, Van Binh N, Kim J. CYTOKININ RESPONSE FACTOR2 (CRF2) and CRF3 regulate lateral root development in response to cold stress in Arabidopsis. ThePlant Cell, 2016, 28: 1828-1843. Genome-Wide Identification of Cucumber ERF Gene Family and Expression Analysis in Female Bud DifferentiationPAN Jian, WEN HaiFan, HE HuanLe, LIAN HongLi, WANG Gang, PAN JunSong, CAI Run (School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai 200240) Abstract:【Objective】The objectives of this research were to identify the Ethylene Response Factor (ERF) family genes from cucumber (Cucumis sativus L.) genome and to know the profile of ERF family such as gene number, gene structure and expression characters in cucumber, so as to provide a theoretical basis for exploring what roles the ERF transcription factors played in female flowers development. 【Method】The ERF genes in cucumber genome were identified by BLAST software in the 9930_V2 genome database based on Arabidopsis ERF genes. EMBOSS, MEME, TBtools, ExPASy and MEGA 7.0 software were used for carrying out various bioinformatics analysis of ERF. The qRT-PCR was used to detect the expression of cucumber ERF gene family in female buds’ differentiation. 【Result】Total of 138 ERF genes were identified from cucumber genome, which could be divided into 10 classes. These 138 ERF genes were named from CsERF1 to CsERF138, and the number of amino acids of these ERF genes was between 126 and 745. Multiple sequence alignments and Motif analysis showed that the CsERF gene family had two conserved domains, namely AP2/ERF domain and B3 domain. There were 19 ERF genes showed different expression between FFMMAA and ffMMAA genotype, 9 of them up-regulated in FFMMAA genotype and the other 10 up-regulated in ffMMAA genotype. The expression trend analysis showed 31 ERF genes up-regulated and 30 genes down-regulated in the initial-phase of female buds’ differentiation. Furthermore, it was demonstrated that CsERF9 and CsERF31 could bind the GCC-box. 【Conclusion】138 ERF gene family members were identified from the cucumber 9930_V2 genome. These entire ERF gene shared AP2/ERF domain. Some of the ERF family members were related to sex determination and female flower development and could bind GCC-box to regulate downstream genes expression. Key words: cucumber; ethylene signal; female flower development; ERF, sex determination doi: 10.3864/j.issn.0578-1752.2020.01.013 开放科学（资源服务）标识码（OSID）： 收稿日期：2019-05-28； 接受日期：2019-10-29 基金项目：国家自然科学基金（31772308，31272185，31471879） （责任编辑 赵伶俐）
| 
|Archiver|手机版|小黑屋|www.openhelp100.com
( 冀ICP备19026749号-1 )
发表于 2021-2-23 21:49:51
