时间:2024-05-23
黄 宁 惠乾龙 方振名 李姗姗 凌 辉 阙友雄 袁照年,*
1 福建农林大学国家甘蔗工程技术研究中心, 福建福州 350002; 2 玉林师范学院农学院, 广西玉林 537000
分蘖或分枝是水稻、小麦、大麦和甘蔗等禾本科植物生长发育过程中普遍存在的一种发育生物学现象[1-2], 也是一种影响作物产量的重要农艺性状,通常受植物激素[3-4]、生物胁迫[5-6]及非生物胁迫调控[7-8]。甘蔗是热带和亚热带地区重要的糖料作物,蔗茎是糖分储藏和收获的重要部位, 获得更多有效分蘖茎, 能够提高甘蔗产量[9]。甘蔗有效分蘖茎是由主茎地下基部未伸长节间的侧芽发育而来的成熟个体[2]。当受黑穗病菌侵染后, 甘蔗分蘖增多、蔗茎变细, 蔗梢长出黑鞭, 引起其有效茎数及蔗茎中蔗糖含量严重下降, 最终导致甘蔗产量降低和品质变劣[10]。
1966年, Cook等[11]首次从棉花根系中分离得到一类萜类小分子化合物—独脚金内酯(strigolactones,SLs), 能够抑制植物分蘖或分枝[12-13], 调控植物对各种生物和非生物胁迫的响应[14]。β-胡萝卜素异构酶(D27)可以产生独脚金内酯前体或中间产物, 调控独脚金内酯生物合成[15]。据报道, 水稻[16]和拟南芥[15]中D27基因突变体均表现为分枝增多。小麦TaD27基因沉默突变体植株分蘖数增多, 而过表达小麦TaD27突变体植株分蘖数减少[17]。
迄今, 关于甘蔗D27基因, 只有1篇基因克隆及其组织特异性和非生物胁迫下表达模式的报道[18],尚未见该基因家族鉴定在甘蔗中的研究报道。本研究首先在甘蔗原始亲本割手密种基因组中鉴定D27基因家族, 并进行生物信息学分析; 其后, 基于植物生长调节剂处理相关和受黑穗病菌侵染相关甘蔗栽培种转录组数据, 对甘蔗栽培种中与割手密种D27基因家族同源基因的表达模式进行分析; 最后,筛选并克隆出1个同时参与分蘖和响应黑穗病菌胁迫的甘蔗栽培种D27基因, 并对该基因进行生物信息学分析, 多种胁迫下的表达定量和亚细胞定位分析。研究旨在为解析D27在甘蔗分蘖过程及与黑穗病菌互作过程中的作用提供理论依据。
Waters等[15]发现, 在91个真核生物D27蛋白序列中, 均包含1个功能未知的保守结构域DUF4033,使用NCBI中的conserved domains工具检索发现包含该结构域的蛋白, 通过将反式-β-胡萝卜素异构成9-顺式-β-胡萝卜素, 参与 SL 的生物合成, 可使用DUF4033检索 D27蛋白家族。以甘蔗原始亲本割手密种基因组的基因序列、蛋白序列及染色体位置信息为基础, 结合甘蔗栽培种转录组序列, 构建甘蔗基因家族鉴定模型。通过在线数据库(http://www.life.illinois.edu/ming/downloads/Spontan eum_genome/)[19]下载甘蔗割手密种基因组相关信息。使用 TBtools软件[20]中的 HMM 检索(Simple HMM Search)工具, 以DUF4033为探针在割手密种基因组的蛋白序列中进行检索, 获得割手密种D27基因家族成员所编码的蛋白序列。同时利用DNAMAN软件中的多序列比对(multiple alignment)工具, 计算割手密种D27基因家族成员间氨基酸序列的相似性。最后分别通过 TBtools软件[20]中的基因位置可视化(Gene Location Visualize)工具和热图(Heatmap)工具绘制割手密种D27基因家族成员的染色体分布图和基因家族成员间氨基酸序列相似性热图。
1.2.1 系统进化树构建及结构域分析 使用interpro在线软件(http://www.ebi.ac.uk/interpro/)检索DUF4033, 分别下载拟南芥(Arabidopsis thaliana)、水稻(Oryza sativasubsp.)、高粱(Sorghum bicolor)、小米(Setaria italica)和玉米(Zea mays)中的 D27s氨基酸序列, 同时利用MEGA7.0软件对甘蔗割手密种及其他5个物种中的D27s氨基酸序列进行多序列比对, 并使用最大似然法(Maximum Likelihood Method)构建系统进化树 (Bootstrap = 1000)。使用 Smart search在线工具 (http://smart.embl-heidelberg.de/smart/set_mode.cgi?NORMAL=1) 检索 6个物种中D27s氨基酸序列的结构域, 并通过TBtools软件[20]中的基因结构查看(gene structure viewer)工具进行展示。
1.2.2 顺式元件预测 基于甘蔗割手密种基因组染色体位置信息, 使用TBtools软件[20]中的Gtf/Gff3序列提取(Gtf/Gff3 sequences Extract)工具, 提取甘蔗割手密种D27基因家族成员CDS上游2000 bp的启动子序列, 将提取到的序列提交至 PlantCare(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)进行顺式元件预测。同时, 使用TBtools软件[20]中的生物序列查看(BioSequence viewer)工具对甘蔗割手密种D27基因家族成员的顺式元件进行可视化。
1.3.1 植物材料及处理 参考 Ling等[22]的方法培养甘蔗主要栽培品种 ROC22组培苗, 分别使用 ABA (100 μmol L-1)、MeJA (100 μmol L-1)、SA(5 mmol L-1)和H2O2(500 mmol L-1) 水溶液处理长势一致的组培苗, 并于0 h和6 h取整株组培苗作为样品, 保存于-80°C。
1.3.2 总 RNA提取和 cDNA合成 采用TRIzol法提取甘蔗总RNA后, 使用DNaseI对RNA样品进行处理, 同时使用 PrimeScript RT-PCR Kit (Perfect for Real Time)合成cDNA。
1.3.3 差异表达ScD27.1基因筛选 项目组前期研究发现, 甘蔗受6-BA (6-benzylaminopurine, 2 mg L-1)和DA-6 (Diethylaminoethylhexanoate, 10 mg L-1)共同诱导后分蘖增多, 同时前人研究表明, 黑穗病菌侵染导致甘蔗分蘖增多[10], 说明植物生长调节剂(6-BA和 DA-6)和黑穗病菌均引起了甘蔗分蘖相关代谢途径调控基因的变化。因此, 本研究以甘蔗割手密种D27基因家族所编码的氨基酸序列为探针,分别在6-BA (2 mg L-1)和DA-6 (10 mg L-1)共同诱导的甘蔗栽培种 ROC22转录组 (数据未公布)和黑穗病菌胁迫的甘蔗栽培种ROC22和YC05-179转录组[21]数据库中进行tblastn比对, 提取与甘蔗割手密种D27基因家族同源的甘蔗栽培种D27基因家族在2个转录组中的表达数据, 筛选2个转录组中均差异表达的甘蔗割手密种D27基因。
1.3.4ScD27.1基因的克隆和分析 以在 6-BA(2 mg L-1)和DA-6 (10 mg L-1)共同诱导(数据未公布)和黑穗病菌胁迫的甘蔗栽培种转录组[21]中均差异表达的甘蔗割手密种D27基因 Sspon.06G0012830-1A的序列作为参考设计PCR扩增引物ScD27.1-F/R (表1), 通过ExTaq酶(TaKaRa, 中国大连)在甘蔗栽培种ROC22整株组培苗cDNA中进行扩增, 扩增程序如下: 94℃ 4 min; 94℃ 30 s, 58℃ 30 s, 72℃ 1 min 30 s, 35个循环; 72℃ 10 min。目的片段PCR产物经过Gel Extraction Kit (Omega, 中国北京) 纯化后连接至pMD19-T载体, 并转化大肠杆菌DH5α进行TA克隆及测序。同时, 使用 ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/)及 blastp (https://blast.ncbi.nlm.nih.gov/Blast.cgi)分析测序结果。
1.4.1 生物信息学分析 分别使用在线工具ExPaSy (http://web.expasy.org/protparam/)和 Prabi(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?pa ge=/NPSA/npsa_hnn.html)预测 ScD27.1一级结构和二级结构; 同时使用 Plant-PLoc (http://www.csbio.sjtu.edu.cn/bioinf/plant/)、ProtComp 9.0 (http://www.softberry.com/berry.phtml?topic=protcomppl&group=programs&subgroup=proloc)和WoLFPSORT (https://wolfpsort.hgc.jp/)工具预测 ScD27.1亚细胞定位; 分别使用 TargetP-2.0 (http://www.cbs.dtu.dk/services/TargetP/index.php)、UbPred (http://www.ubpred.org/)和NetPhos 3.1 Server (http://www.cbs.dtu.dk/services/NetPhos/)对 ScD27.1进行靶向肽类型、泛素化位点及磷酸化位点分析。
1.4.2 qRT-PCR表达分析 参考甘蔗割手密种D27基因家族顺式元件预测结果, 分析参与植物生长发育和逆境响应的植物激素信号分子 ABA、MeJA、SA[32-33]和H2O2[34]如何调控甘蔗D27基因的转录表达。使用 qRT-PCR技术, 以 q-ScD27.1-F/R(表1)为引物,CAC和CUL为内参基因[22], 分析ScD27.1基因在ABA、MeJA、SA和H2O2胁迫下的表达模式。为了获得准确基因相对表达分析结果,在3倍梯度稀释的ROC22 cDNA样品中对ScD27.1基因的定量引物进行特异性扩增检测和扩增效率检测。基因的相对表达量计算及差异显著性分析分别使用 2-ΔΔCT算法[23]和新复极差法 DPS (7.05)。qRT-PCR反应体系及扩增程序均参照 SYBR Green PCR Master Mix试剂盒说明书。
1.4.3 亚细胞定位分析 参考 Gateway BP Clonase II Enzyme mix (Invitrogen, 美国加利福尼亚)说明书, 以 pMD19-T-ScD27.1质粒为模板,G-ScD27.1-F/R (表1)为引物, 构建入门载体pDONR-Zeo-ScD27.1。根据Gateway LR Clonase II Enzyme mix (Invitrogen, 美国加利福尼亚)说明书,通过 LR反应构建亚细胞定位载体pFAST-R06-ScD27.1。将 pFAST-R06和pFAST-R06-ScD27.1质粒分别转入农杆菌 GV3101并涂板(含50 μg mL-1壮观霉素和 35 μg mL-1利福平), 28℃倒置培养 2 d后, 进行菌液 PCR鉴定。分别将含有pFAST-R06 和 pFAST-R06-ScD27.1质粒的农杆菌注射侵染烟草叶片表皮细胞(操作参考 Su等[24]),48 h 后, 用 16.5 μmol L-1的渥曼青霉素(Wortmannin,Wort)处理烟草叶片, 并在0 h和0.5 h后[25]通过激光共聚焦显微镜(Leica, 德国韦茨拉尔)观察烟草叶片表皮细胞[26]。
通过HMM共检索到5个具有完整DUF4033保守结构域的甘蔗割手密种D27序列, 分别为 Sspon.06G0012830-1A、Sspon.06G0012830-3C、Sspon.06G0016140-2C、Sspon.07G0014180-2B和 Sspon.07G0014180-3C, 分别定位于6A、6C、6C、7B、7C染色体(图1-A)。多序列比对结果显示, 甘蔗割手密种D27基因家族成员间的氨基酸序列一致性为 52.36%, 其中Sspon.07G0014180-2B与Sspon.07G0014180-3C一致性最高, 为 90.09%, Sspon.06G0016140-2C与 Sspon.07G0014180-2B一致性最低, 为20.41% (图1-B), 说明甘蔗割手密种D27基因家族氨基酸序列间差异较大。
表1 引物列表Table 1 List of primers
2.2.1 系统进化树构建及结构域分析 使用甘蔗割手密种的5个D27s氨基酸序列与其他5个物种(拟南芥、水稻、高粱、小米及玉米)的62个D27s氨基酸序列共同构建系统进化树, 结果如图2所示。6个物种的D27s蛋白聚类成7个系统发育分支, 定位于甘蔗割手密种7B染色体的Sspon.07G0014180-2B和7C染色体的Sspon.07G0014180-3C分布于分支I, 定位于6A染色体的Sspon.06G0012830-1A和6C染色体的 Sspon.06G0012830-3C分布于分支 V, 定位于6C染色体的Sspon.06G0016140-2C分布于分支VII,并且甘蔗割手密种D27s与高粱D27s 在系统发育分支中距离最近。Smart search检索到的结构域包括Pfam: DUF4033、Pfam: COesterase和Pfam: Abhydrolase_3等, 说明D27基因家族除了包含典型的 Pfam:DUF4033结构域外, 还可能包含其他结构域, 但这些结构域在进化过程中不保守。
2.2.2 顺式元件预测 为了解甘蔗割手密种D27基因家族成员的调控功能, 使用PlantCare对甘蔗割手密种D27基因家族成员CDS上游2000 bp的启动子序列进行分析, 同时提取与胁迫响应/诱导相关的14类顺式元件(图3)。其中与植物激素响应相关的顺式元件包括 ABA响应元件、MeJA响应元件、GA响应元件和 SA响应元件。与植物生长发育相关的顺式元件包括玉米醇溶蛋白代谢调控元件、分生组织表达调控元件、种子特异性调控元件和胚乳表达调控元件。与胁迫响应相关的顺式元件包括低温响应元件、厌氧诱导元件、光响应元件、干旱诱导元件、防御和应激反应元件和伤害响应元件。甘蔗割手密种D27基因家族所有成员的启动子区均含有厌氧诱导元件、光响应元件、MeJA响应元件和 GA响应元件。Sspon.06G0012830-3C和Sspon.06G0016140-2C的启动子区分别含有特有的种子特异性调控元件和胚乳表达调控元件。推测厌氧诱导元件、光响应元件、MeJA响应元件和GA响应元件相关的反式作用因子可能可以结合到甘蔗割手密种D27基因家族启动子区, 从而受植物激素响应和胁迫诱导表且部分成员在特定组织器官中特异性表达。
2.3.1 差异表达ScD27.1基因筛选 经tblastn比对后, 甘蔗割手密种D27基因家族在植物生长调节剂诱导甘蔗栽培种ROC22分蘖转录组数据中存在3条相似性较高的序列, 分别为 c116894.graph_c0、c112688.graph_c0和 c161113.graph_c0; 在甘蔗栽培种ROC22和YC05-179应答黑穗病菌侵染转录组数据中存在 3条相似性较高的序列, 分别为Sugarcane_Unigene.66941、Sugarcane_Unigene.58014和Sugarcane_Unigene.87626。D27s序列间相似性详见表2。
6-BA和 DA-6诱导甘蔗栽培种 ROC22分蘖转录组中,c112688.graph_c0(P= 0.006)和c116894.graph_c0(P= 0.02)经6-BA和DA-6处理后的表达量,相对于对照组均显著性下调表达, 且在处理20 d时下调表达幅度最大, 分别为对照组的-0.54和-0.42倍(图4-A)。受黑穗病菌侵染后,Sugarcane_Unigene.58014在甘蔗抗黑穗病品种YC05-179和甘蔗感黑穗病品种 ROC22中均无显著差异表达;尽管Sugarcane_Unigene.66941在甘蔗抗黑穗病品种YC05-179受侵染后48 h和120 h的下调表达并不显著, 但其在感黑穗病品种 ROC22受侵染后 48 h和120 h的上调表达均达到显著水平(P= 0.01) (图4-B)。此外, 由于c161113.graph_c0和Sugarcane_Unigene.87626在转录组中的 FPKM 值均小于 0.2,属于极低表达, 不符合数据分析要求。
表2 甘蔗割手密种与栽培种D27s间序列相似性Table 2 Sequence similarity between D27s from Saccharum spontaneum and sugarcane cultivar
综上所述, 与甘蔗割手密种D27基因家族中Sspon.06G0012830-1A相似性分别为 97.93%和98.34%的甘蔗栽培种c116894.graph_c0和Sugarcane_Unigene.66941(表2)分别在 6-BA 和DA-6共同诱导甘蔗栽培种ROC22分蘖转录组和甘蔗栽培种 ROC22应答黑穗病菌侵染转录组中显著差异表达, 推测与甘蔗割手密种Sspon.06G0012830-1A同源的甘蔗栽培种D27基因同时参与甘蔗分蘖调控及黑穗病菌响应, 可用于下一步的研究。
2.3.2ScD27.1基因克隆 根据琼脂糖凝胶电泳结果可知, PCR扩增获得1条约1127 bp的特异条带(图5-A)。PCR产物经连接, 转化后进行测序分析。测序结果经 ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/)分析表明, 克隆所得的基因长1127 bp, 编码 266个氨基酸(图5-B)。经 blastp(https://blast.ncbi.nlm.nih.gov/Blast.cgi)比对后发现,其氨基酸序列与高粱 β-胡萝卜素异构酶 (GenBank登录号为 XP_002444917.2)相似性为 91.82%, 因此将该基因命名为甘蔗 β-胡萝卜素异构酶基因(ScD27.1, GenBank登录号为MT499895)。
2.4.1 生物信息学分析 生物信息学分析结果显示, 甘蔗栽培种ScD27.1编码266个氨基酸, 等电点为8.91, 分子量为30.00 kD, 属于不稳定蛋白。二级结构主要包括α-螺旋和无规则卷曲(表3)。甘蔗栽培种 ScD27.1具有叶绿体转运肽可能定位于叶绿体,包含4个泛素化位点和18个磷酸化位点(表3)。
2.4.2 qRT-PCR表达分析ScD27.1基因引物扩增效率为 1.026, 在 0.95~1.05之间, 扩增曲线具有很好的相关性(R2= 0.9991), 说明引物具有良好的扩增特异性(图6-A)。受MeJA胁迫后,ScD27.1基因的表达量相对于对照组下调表达; 受ABA、SA及H2O2胁迫后,ScD27.1基因的表达量相对于对照组上调表达, 并且在受ABA和H2O2胁迫后显著高于对照组,分别为对照组的 5.37倍和 4.60倍(图6-B)。推测ScD27.1基因受 ABA及 H2O2显著诱导表达, 但对MeJA、SA胁迫响应不明显。
表3 生物信息学分析结果Table 3 Result of bioinformatics analysis
2.4.3 亚细胞定位分析 烟草叶片表皮细胞亚细胞定位结果如图7所示, 空载35S::GFP定位于细胞膜及细胞核, 35S::ScD27.1::GFP定位在细胞膜上,胞内也出现了点状荧光, 如叶绿体。Wort能够促进细胞内液泡前体或者运输小泡发生融合形成更大的膜泡结构[25,27]。经过Wort处理0.5 h后观察发现, 空载35S::GFP的定位未发生改变, 烟草叶片表皮细胞中未观察到具有绿色荧光的膜泡结构。然而, 尽管35S::ScD27.1::GFP的定位也未发生改变, 但在烟草叶片表皮细胞中可以观察到含有绿色荧光的膜泡结构, 因此, ScD27.1可能定位于细胞膜和叶绿体, 并参与细胞内膜泡运输或由液泡前体、胞内运输小泡完成其在在胞质内的分拣运输过程。
SLs是一种广泛存在、能够抑制植物分蘖或分枝的植物激素[12-13]。目前, 在拟南芥[15]、水稻[16]和小麦[17]中均有关于SLs合成相关基因d27突变体中分蘖增多的报道, 说明D27基因能够通过抑制 SLs合成来减少植物分蘖或分枝。本研究中,ScD27.1在分蘖相关转录组中的转录本 c116894.graph_c0, 经6-BA和DA-6 (具有促进分蘖功能)处理后相对于对照组, 显著下调表达, 与前人研究结果相同[15-17]。
SLs除了在生长发育过程中起作用外, 还被认为能够响应病原菌侵染[28]。番茄SLs合成相关基因Slccd8突变体对灰葡萄球菌和链格孢菌均更敏感[29]。拟南芥 SLs信号转导相关基因Atmax2突变体对果胶杆菌和假丁香单胞菌的敏感性增加, 但对灰葡萄球菌的敏感性不增加[30]。拟南芥SLs生物合成(max1、max3和max4)和信号转导(max2)突变体对带化红球菌更加敏感[31]。甘蔗 SLs合成相关基因ScD27.1的同源基因Sugarcane_Unigene.66941在抗病品种中受黑穗病菌抑制表达, 在感病品种中受黑穗病菌诱导表达, 说明ScD27.1响应黑穗病菌侵染,但其在甘蔗分蘖及甘蔗与黑穗病菌互作过程中的作用还需进一步试验验证。
ABA、MeJA和SA作为重要的植物激素信号分子, 不仅能够调控植物生长发育, 同时参与植物逆境反应[32-33]。活性氧主要包括 O2-和 H2O2, 也是植物中一种重要的信号分子, 在植物受病原体侵染或其生长发育过程中发挥作用[34]。ScD27.1基因受ABA、SA及 H2O2胁迫后上调表达, 受 MeJA胁迫后下调表达, 但SA及MeJA胁迫下差异表达不显著,推测ScD27.1基因受 ABA和H2O2诱导表达, 不受MeJA和SA调控。Abuauf等[35]在ABA处理拟南芥侧根后对AtD27进行表达分析, 得到了相似的结果。苜蓿根系在营养缺乏条件下, 激发活性氧生成的NADPH氧化酶活性增强, SLs合成相关基因MtCCD7、MtCCD8、MtD27基因上调表达[36], 说明活性氧与 SLs合成相关基因可能存在协同作用, 这与ScD27.1基因受H2O2诱导表达的研究结果相似。真核生物中由转录因子等反式作用因子结合基因启动子区启动基因表达转录, 然而一般情况下反式作用因子常常形成同源或异源复合体行使其功能, 基因的表达往往不受单一顺式元件调控。因此,ScD27.1基因的割手密同源基因启动子区域虽然具有SA和MeJA响应元件, 但ScD27.1也有可能在某些条件下不响应SA和MeJA的诱导表达。
D27同源基因存在于藻类(蓝藻)和高等植物,但不存在于动物或真菌中, 是一种植物特有蛋白[16]。本研究首次从甘蔗割手密种基因组中鉴定出D27基因家族的5个成员, 并从甘蔗栽培种ROC22中克隆出1个同时参与分蘖和黑穗病菌胁迫响应的割手密种同源D27基因, 该甘蔗栽培种ScD27.1基因编码266个氨基酸。该蛋白等电点为 8.91, 分子量为30.00 kD, 具有DUF4033结构域, 与已报道的编码288个氨基酸, 等电点为 5.04, 分子量为 71.58 kD,具有2个锌指蛋白结构域的甘蔗D27基因[18]同源性为31.40%, 推测2个甘蔗D27基因同属于甘蔗D27基因家族, 可能具有不同的生物学功能。目前关于D27亚细胞定位的报道较少, Waters等[16]发现, 拟南芥D27定位于质体, 水稻D27定位于叶绿体, 本研究中甘蔗 ScD27.1定位于细胞膜和叶绿体。随着越来越多的D27s被发现, D27的亚细胞定位研究将会更加深入。甘蔗割手密种D27基因家族顺式元件预测结果显示,D27s参与植物激素响应、生长发育调控及胁迫响应。其中关于D27参与ABA、MeJA和SA激素响应已在本研究中验证, 干旱胁迫[37]、光胁迫[38]和磷饥饿胁迫[39]调控已有相关报道。进一步的研究中, SL在低温胁迫[40]、种子萌发[41]和果实发育[42]过程中的作用可作为D27s相关调控功能研究的参考。
本研究首次鉴定出5个甘蔗割手密种D27基因家族成员, 并对其进行染色体定位、系统进化、结构域及顺式元件分析。其后, 筛选并克隆到 1个在甘蔗分蘖转录组和甘蔗应答黑穗病侵染转录组中均差异表达的甘蔗栽培种ScD27.1基因。该基因编码266个氨基酸, 为等电点8.91、分子量30.00 kD的不稳定蛋白, 二级结构主要包括 α-螺旋和无规则卷曲。表达分析显示, 该基因受 ABA和 H2O2诱导表达, 不受 MeJA和 SA调控。亚细胞定位明确了ScD27.1蛋白定位于细胞膜和叶绿体。以上结果不仅为解析甘蔗ScD27.1基因在信号传导机制及胁迫响应机制中的功能作用提供参考, 还为深入了解甘蔗受黑穗病菌侵染后的抗病机制与分蘖表型之间的关系研究提供新思路。
我们致力于保护作者版权,注重分享,被刊用文章因无法核实真实出处,未能及时与作者取得联系,或有版权异议的,请联系管理员,我们会立即处理! 部分文章是来自各大过期杂志,内容仅供学习参考,不准确地方联系删除处理!