Cloning and function analysis of LsMYB7 gene in Lycoris sprengeri

2019年8—9月于浙江农林大学石蒜属植物种质资源圃采集换锦花花瓣，取换锦花小花苞(2.0~3.5 cm)、大花苞(4.5~6.0 cm)、盛花期和败花期4个不同花发育时期(图1A)和不同花色无性系H1、H2、H3、H4和H5盛花期花瓣(图1B)，液氮速冻后于−80 ℃冰箱保存备用。

Figure 1.  Different stages of flower development (A) and peak flowering stage in different clonal plants (B)

采用RNAiso plus法提取换锦花花瓣总RNA^[15]，参照PrimeScript 1st Strand cDNA Synthesis Kit试剂盒(Takara)方法合成cDNA 第1链。

根据转录组测序信息设计LsMYB7基因 cDNA序列特异引物LsMYB7-F：CAAGCAGTGGTCTCAACA和LsMYB7-R：AGAACAGCACTACTAAAGGT，参照Premix PrimeSTAR HS的PCR扩增体系扩增目的基因cDNA序列，PCR扩增反应程序为94 ℃变性30 s；58 ℃ 退火30 s；72 ℃延伸90 s，32个循环，72 ℃延伸10 min，质量浓度为1.0%琼脂糖凝胶电泳检测PCR扩增产物，于凝胶成像系统(Bio-RAD)观察拍照；采用Hingene琼脂糖凝胶回收试剂盒(杭州麦克德勒科技有限公司)回收目的DNA，于−20 ℃保存备用。

参照pEASY-Blunt Zero Cloning Kit说明书将目的基因LsMYB7序列连接到pEASY-Blunt载体，转化大肠埃希菌Escherichia coli DH5α感受态细胞；在含有氨苄青霉素(Amp)、异丙基-β-D-硫代半乳糖苷(IPTG)和5-溴-4-氯-3-吲哚 β-D-半乳糖苷(X-Gal)的LB培养基(Luria-Bertani medium)上培养过夜，挑取阳性克隆于LB液体(含50 mg·L⁻¹ Amp)培养基中，37 ℃振荡培养，PCR检测呈阳性的菌液送浙江有康生物科技有限公司测序；采用质粒DNA提取试剂盒(杭州创试生物科技有限公司)提取目的序列重组质粒DNA，于−20 ℃保存。

采用Primer 5.0设计LsMYB7基因和花色形成相关基因(LsCHS、LsF3H、LsANS、LsUFGT1和LsUFGT2)的RT-qPCR引物(表1)，使用PrimeScriptTM RT reagent Kit with gDNA Eraser试剂盒(Takara)方法合成cDNA 第1链。参照SYBR® Premix Ex TaqTMⅡ(Takara)方法，于CFX96TM荧光定量 PCR 仪(BIO-RAD)进行RT-qPCR验证。RT-qPCR体系(20.0 uL)：cDNA(＜100 ng) 1.6 uL，正反引物各0.8 uL，TB Green Premix Ex Taq Ⅱ 10.0 uL，RNase Free ddH₂O 6.8 uL。RT-qPCR程序为95 ℃ 30 s，95 ℃ 3 s，60 ℃ 30 s，循环40次。以LsGADPH为内参基因^[15]，2^−ΔΔCt方法计算实时荧光各基因的相对表达量，重复3次。

参照刘跃平等^[17]方法提取并用高效液相色谱(HPLC)测定换锦花花瓣花色苷质量浓度。精密称量换锦花花瓣干粉0.040 0 g，加入2 mL体积分数为1%甲醇/HCl提取溶液，振荡混匀；20 ℃超声提取30 min；12 000 r·min⁻¹，离心15 min，取上清液；0.45 μm滤膜过滤；梯度洗脱：流动相为甲醇(A)-体积分数为1%的甲酸水(B)，0~20 min(A体积分数从5%升至60%)，20~25 min(A体积分数从60%升至100%)，25~30 min(A体积分数保持100%)，流速1 mL·min⁻¹，检测波长280 nm，柱温30 ℃，进样量10 uL；以矢车菊素-3-O-葡萄糖苷、天竺葵素-3-O-葡萄糖苷、飞燕草素-3-O-葡萄糖苷等3个花色苷为标准品。每样3个生物学重复。

采用XbaⅠ和BamHⅠ分步酶切法酶切亚细胞定位载体pAN580，再采用ClonExpress Ⅱ One Step Cloning Kit (杭州霆喜生物科技有限公司)将LsMYB7基因序列连接到pAN580上，并转化大肠埃希菌DH5α感受态细胞，在含有50 g·L⁻¹Amp的LB培养基上进行筛选，挑选阳性克隆进行菌液PCR检测，成功构建pAN580-LsMYB7亚细胞定位载体，采用无内毒素质粒提取试剂盒(AxyPrep)提取pAN580-LsMYB7重组载体质粒DNA；聚乙二醇法(PEG 4000)转化烟草Nicotiana tabacum原生质体，用激光共聚焦荧光显微镜观察拍照。

采用双酶切法构建VIGS沉默载体pTRV2-LsMYB7：利用Primer 5.0设计长为400 bp的LsMYB7基因cDNA序列插入片段引物pTRV-LsMYB7-F: TACCGAATTCTCTAGAGAGGACAACATGCTACAATCCC和pTRV-LsMYB7-R：GCTCGGTACCGGATCCGCTCGAATCGCTAACATCCG；用限制性内切酶XbaI和BamHI分别双酶切VIGS载体(pTRV2)与LsMYB7基因的插入序列，T4 DNA连接酶连接载体与目的序列，转化大肠埃希菌DH5α，并在含有50 g·L⁻¹ 卡那霉素(Kan)的LB固体培养基筛选，提取pTRV2-LsMYB7重组质粒DNA。

采用微量注射法转化换锦花花苞^[12]，将重组质粒pTRV2-LsMYB7、pTRV2和pTRV1质粒DNA转化根癌农杆菌Agrobacterium tumefaciens GV3101，于YEP液体培养基(50 mg·L⁻¹ Rif和50 mg·L⁻¹ Kan)中，28 ℃，220 r·min⁻¹振荡培养12~14 h；4 000 r·min⁻¹, 离心10 min，去上清液；加入侵染液(10 mmol·L⁻¹ MES + 200 μmol·L⁻¹ AS + 10 mol·L⁻¹ MgCl₂，pH 5.6±0.03)重悬后的菌液D(600)=0.8~1.0；将pTRV1与pTRV2空载体、pTRV2-LsMYB7农杆菌液按照体积比1∶1混合均匀，室温黑暗静置4~6 h；使用1 mL注射器，吸取1 mL混合农杆菌菌液。选取3~4 d的换锦花花苞，在花苞基部缓慢注射混合根癌农杆菌菌液，充分侵染换锦花花苞，侵染5 d后采集花瓣用于LsMYB7和花色形成相关结构基因的表达分析(见1.2.4和表1)。

采用Excel 2020 统计和整理数据，采用SPSS 20.0 对数据进行差异显著性分析，采用Graphpad 8.0作图。

以换锦花花瓣的cDNA为模板，采用RT-qPCR技术扩增获得长度为951 bp 的LsMYB7 cDNA序列(图2)。该cDNA序列含有1个825 bp的开放阅读框(ORF)，编码274个氨基酸，分子量为6.84 kD，蛋白分子式为C₂₄₀₂H₃₉₈₁N₈₂₅O₉₈₂S₂₅₈，理论等电点(pI)为5.04，不稳定系数为43.43，总体亲水性(GRAVY)为0.936。

Figure 2.  Full length of RT-qPCR

通过与美国国家生物技术信息中心(NCBI)搜索的其他植物MYB氨基酸同源序列进行比对，发现换锦花LsMYB7蛋白含有2个R2R3-MYB典型的SANT结构域R2和R3^[18](图3)，属于R2R3-MYB类转录因子。与茶Camellia sinensis KAF5956278.1、狭叶油茶Camellia lanceoleosa KAI8015846.1、粗柄象腿蕉Ensete ventricosum RRT35038.1/RWV93016.1、椰子Cocos nucifera XP_KAG1363931.1、番红花Crocus sativus QBF29477.1、香根鸢尾Iris pallida KAJ6804060.1/KAJ6829988.1、河岸葡萄Vitis riparia XP_034676575.1和拟南芥Arabidopsis thaliana等其他植物的同源性为44.52%~60.57%，其中，与茶KAF5956278.1的同源性最高达60.57%，其次是香根鸢尾 KAJ6804060.1/KAJ6829988.1(59.32%)，与拟南芥AtMYB44(Q9FDW1.1)、AtMYB70(AEC07437.1)、AtMYB73(O23160.1)和AtMYB77(Q98N12.1)的同源性为44.52%~47.44%，且同源部分均集中在N端的R2R3 DNA结合结构域，而C端的同源性较低。

Figure 3.  Comparison of R2R3-MYB amino acid homologous sequences between LsMYB7 and other plants, R2R3-MYB

从NCBI数据库下载拟南芥(AtMYB 16条)、换锦花(LsMYB 2条)和石蒜L. radiata (LrMYB 2条) R2R3-MYB氨基酸序列进行比对并构建系统进化树(图4)。参考拟南芥MYB基因家族的分类方法^[19]，LsMYB7与拟南芥R2R3-MYB S22亚家族的AtMYB44、AtMYB70、AtMYB73和AtMYB77 聚为一类，由于拟南芥S22亚家族参与调控拟南芥的生长和发育过程，响应高盐、干旱低温等非生物胁迫反应，同一亚族基因的功能相似，结合花色基因表达分析，推测LsMYB7可能通过响应干旱和高温等调控换锦花花瓣花色苷的形成。

Figure 4.  Phylogenetic tree between LsMYB7 and R2R3-MYB from A. thaliana

利用HPLC分别测定换锦花4个花发育时期和5个不同花色无性系(H1、H2、H3、H4和H5)盛花期花瓣的花色苷质量浓度(图5A和5B)，结果表明：换锦花花瓣花色苷的主要成分为矢车菊素，含有少量的天竺葵素和飞燕草素，说明花色苷的种类和质量浓度决定换锦花花色的多样性。随着换锦花花器官的生长发育，花瓣花色苷总质量浓度和矢车菊素质量浓度呈逐渐下降的趋势，小花苞时期矢车菊素质量浓度大约是败花期的3倍，说明换锦花花瓣花色苷的积累可能主要在花瓣发育的早期完成(图1A和5A)；在浅色换锦花无性系H4和H5中，矢车菊素质量浓度明显低于深色换锦花无性系H1、H2和H3(图1B和5B)。因此，从换锦花不同花发育时期和不同花色无性系花瓣颜色和花色苷质量浓度来看，矢车菊素对换锦花花瓣的颜色影响最大，矢车菊素质量浓度越高，换锦花花瓣的颜色就越深。

Figure 5.  Anthocyanin contents in petals of different stages of flower development (A) and different clones (B) of L. sprengeri

利用RT-qPCR对LsMYB7和换锦花花色形成相关基因(LsC4H、LsCHS、LsF3H、LsF3'H、LsANS和LsUFGT1、LsUFGT2)在换锦花不同发育时期花瓣中的表达进行分析，结果(图6)表明：LsMYB7与LsCHS和LsF3'H基因的表达随着换锦花花苞发育呈逐渐上升趋势，LsC4H、LsCHS、LsUFGT1、LsUFGT2、LsF3'H和LsMYB7在败花期大量表达，LsF3H则在花苞发育前期几乎无表达，盛花期开始大量表达，而在败花期表达量开始下降；LsANS、LsUFGTs等花色苷合成的后期基因在盛花期的表达最低。其中LsCHS和LsF3'H基因的表达与花色苷总质量浓度和矢车菊素质量浓度正好相反(图5A)。

Figure 6.  Expression of LsMYB7 and structure genes related to anthocyanin biosythesis in different flower development stages

LsMYB7和换锦花花色形成关键基因LsC4H、LsCHS、LsF3H、LsF3'H、LsANS、LsUFGT1和LsUFGT2在换锦花不同花色无性系中的表达差异显著(图7)。2个花色苷合成的前期基因LsC4H和LsCHS及后期基因LsANS在淡色的无性系换锦花花瓣中表达量高。LsF3H在蓝色为主的H3无性系中表达量达到最高。LsMYB7在H1和H4中表达量比较高。LsMYB7基因的表达与LsCHS和LsF3'H基因的表达正好相反，而与不同花色无性系花色苷总质量浓度和矢车菊素质量浓度基本一致，这一结果与不同花发育时期的基因表达结果一致。因此，LsMYB7可能对LsCHS和LsF3'H基因的表达有一定的调控作用。

Figure 7.  Expression of LsMYB7 and structure genes related to anthocyanin biosythesis in different flower clones

通过双酶切法构建亚细胞定位载体pAN580-LsMYB7，转化烟草原生质体，于激光共聚焦荧光显微镜观察，LsMYB7荧光信号定位在细胞核中(图8)，说明LsMYB7基因为转录因子基因，在细胞核中起转录调控的作用。

Figure 8.  Subcellular localization of LsMYB7

构建LsMYB7的VIGS基因沉默载体pTRV2-LsMYB7，通过微量注射法转化换锦花花苞，LsMYB7 基因沉默后，换锦花同朵花中一半花瓣明显变短，且颜色变深(图9A)；对LsMYB7和花色苷形成相关基因在换锦花花瓣中的表达进行分析(图9B和C)，与ck (空载)相比，LsMYB7基因换锦花花瓣中的表达量明显下降，同时，花色苷形成相关基因LsCHS、LsF3'H、LsANS、LsUFGT1和LsUFGT2的表达量极显著下降(P＜0.01)，而LsF3H基因的表达量反而上升，LsC4H基因表达变化不大，说明LsMYB7转录因子可能参与调控花瓣的生长发育，且对LsCHS、LsF3'H、LsANS、LsUFGT1和LsUFGT2的表达起正调控作用，而对LsF3H基因的表达起负调控作用。

Figure 9.  Phenotype and relative expression levels of LsMYB7 and genes related to anthocyanin biosynthesis after LsMYB7 silenced

植物R2R3-MYB转录因子广泛参与调控次生代谢、细胞形态发生、激素刺激、环境胁迫应答、分生组织形成和细胞周期等过程^[20]。根据C-末端的不同，拟南芥R2R3-MYB基因家族可分成22个不同的亚族，其中，拟南芥S4、S5、S6和S7亚族基因参与花青素和类黄酮类化合物生物合成途径的结构基因的转录调控^{[18, 21−22]}，换锦花LsMYB4、LsMYB5和拟南芥S4亚族基因聚为一类。通过花青素形成相关基因表达和VIGS基因沉默技术研究表明，LsMYB4和LsMYB5对花青素生物合成基因的表达有负调控的作用^[10−12]。

S22亚族基因AtMYB44、AtMYB77、AtMYB73和AtMYB70主要参与拟南芥响应高盐、干旱、低温等非生物胁迫反应，其中，AtMYB73基因启动子中含有ABA响应元件ABRE及干旱胁迫和热胁迫顺式作用元件。AtMYB73突变体atmyb73在干旱胁迫处理下，ABA 下游基因ABI2、ABI5 的表达水平均较野生型明显增强。外源 ABA 处理野生型和突变体种子、幼苗，获得了与干旱胁迫处理类似的结果^[23−24]。本研究结果表明：LsMYB7蛋白含有2个R2R3-MYB典型的SANT结构域R2和R3，属R2R3-MYB类转录因子，且R2和R3的保守结构域与其他植物高度同源，系统进化树分析结果显示LsMYB7和S22亚族基因AtMYB44、AtMYB77、AtMYB73和AtMYB70聚为一类，推测LsMYB7可能与S22亚族基因有相似的功能。

本研究中，LsMYB7基因的表达与换锦花花色苷形成相关基因LsCHS、Ls4CL2和LsUFGT2的表达趋势一致，推测LsMYB7可能参与换锦花花瓣花色苷的生物合成。而LsMYB7属于拟南芥S22亚族，未见该亚族基因AtMYB44、AtMYB77、AtMYB73和AtMYB70参与花色苷的生物合成的转录调控。已有报道认为AtMYB73/44能够参与调控拟南芥对干旱胁迫的响应^[25−26]，而干旱胁迫可诱导植物细胞合成和积累花色苷。花色苷的光化学性质、亚细胞积累位点及在植物器官、组织中的空间分布决定了花色苷能强化植物的耐旱性，其中，花色苷提高植物细胞在干旱胁迫下的抗氧化能力可能是花色苷强化植物耐旱性的主要原因^[27]。有研究表明：干旱胁迫可激活紫麦Triticum aestioum ‘Guizi 1’、甘薯Ipomoes batats等植物花色苷合成相关基因表达，通过提高花青素含量抵御干旱胁迫^[28−29]。本研究中LsMYB7 基因沉默后，花色苷形成相关基因LsCHS、LsF3'H、LsANS、LsUFGT1和LsUFGT2的表达明显下降，说明LsMYB7可能直接或间接调控花色苷的积累，而LsMYB7与S22亚族基因聚为一类，由于换锦花开花季为夏末秋初的8—9月，此时多为高温和干旱季节，因此推测LsMYB7可能通过对干旱胁迫响应来调控花色苷形成相关基因的表达，大量积累花色苷。这一研究结果与云南文山辣椒Capsicum annuum、番茄Lycopersicon esculentum的研究^[30−31]相似，因此，植物可以通过干旱胁迫激活与花色苷合成相关基因的表达水平促进花色苷积累，进而提高抗旱能力。

本研究通过RT-qPCR获得长为951 bp 的LsMYB7 cDNA序列，开放阅读框(ORF) 825 bp，编码274个氨基酸，LsMYB7为R2R3-MYB转录因子家族，定位于细胞核，与其他植物R2R3-MYB有较高的同源性。系统进化分析表明LsMYB7和参与调控干旱等非生物胁迫响应的S22亚族基因聚为一类；LsMYB7基因主要在败花期和花色苷含量较高的H1无性系中表达，与花色苷合成相关基因的表达趋势一致；LsMYB7基因沉默后，换锦花部分花瓣明显变短，颜色变深，LsCHS、LsF3'H、LsANS、LsUFGT1和LsUFGT2等花色苷形成相关基因的表达显著下调，因此，LsMYB7参与调控换锦花花瓣的生长发育，且通过对LsCHS、LsF3'H、LsANS、LsUFGT1和LsUFGT2花色苷形成相关基因的表达调控花色苷的积累。此外，LsMYB7可能通过响应干旱胁迫激活花青素合成相关基因表达促进花色苷积累，进而提高换锦花的抗旱能力，LsMYB7调控花色苷积累抵御干旱的机制需要进一步研究。