PEG模拟干旱胁迫对冬樱花幼苗叶片生理特性的影响
Effect of Simulated Drought Stress on the Physiological Characteristics of Cerasus cerasoides Leaves
-
Keywords:
- Cerasus cerassoides (D. Don) Sok. /
- seeding /
- PEG-8000 /
- drought stress /
- leaf /
- physiological characteristics
-
雄性增强的抗原-1基因(the male-enhanced antigen-1 gene,MEA1)最初是从用抗雄性小鼠血清细胞构建的原核表达cDNA文库中分离获得[1]。分析其结构,发现MEA1基因序列侧面有2个基因:Ppp2r5d和Peas,使得Peas-MEA1-Ppp2r5d 成为1个基因复合体[2]。有研究报道MEA1基因位于人6号染色体,位于小鼠17号染色体,位于牛23号染色体,在哺乳动物系统发育过程中是比较保守的[3]。通过原位杂交,揭示了MEA1在生精管上皮细胞的粗线期高表达[4]。免疫组化分析显示该基因在精子生成后期发生转录[4]。在小鼠的研究中表明:MEA1在睾丸中明显高表达,在心脏、肝脏、大脑中微量表达[5]。Northern-blot结果显示:该基因编码大约1.5 kb的单个转录本,当葡萄糖、棉子糖或果糖存在的情况下,有表达量[6]。MEA2被报道在小鼠的精母细胞中表达[7]。通过分析该基因的SNP位点及其与生精过程中各个激素指标的相关性,其中两个位点可作为选育性欲较强的雄性动物的分子标记[8]。
版纳微型猪近交系(Banna mini-pig inbred line,BMI)作为一个大型哺乳动物近交系,其遗传背景较为清楚,生理、病理特征与人类相似,在临床医学、人类疾病模型、异种器官移植等方面有一定的价值[9-13]。但是在正常繁育BMI的37年来,有多个亚系断代,前期研究发现是部分公猪不育造成的,这些不育公猪已成为阻碍其群体进一步扩大的障碍,因此对版纳微型猪近交系雄性不育的研究很有必要。由于MEA1基因在哺乳动物精子生成后期发挥着重要作用,本研究通过克隆版纳微型猪近交系MEA1基因,确定其序列特征;进而对蛋白质进行功能生物信息学分析预测其功能,通过mRNA的多组织表达谱确定其发挥功能的重要组织,对该基因在BMI公猪精子生成方面的功能研究奠定一定的基础。
1. 材料与方法
1.1 材料
屠宰云南农业大学版纳微型猪近交系实验猪场的成年公猪,取睾丸和附睾(性腺),前列腺、精囊腺和尿道球腺(副性腺),5种雄性生殖重要相关组织;取心脏、肝脏、脾脏、肺脏、肾脏、胃脏、脑、肌肉、十二指肠(小肠)、结肠(大肠)等常规组织;共15种组织样品。
1.2 方法
1.2.1 总RNA提取和第一链cDNA合成
参照说明书利用柱式RNA提取试剂盒(TaKaRa:9767)分别提取各组织总RNA,检测所提RNA的浓度、纯度及完整性,合格样品按反转录试剂盒(TaKaRa:6210A)操作说明书反转录为第一链cDNA保存备用。
1.2.2 引物设计及合成
根据GenBank下载的猪(NM_001110429)、牛(NM_174648)、小鼠(NM_001277309)、人(NM_014623)等物种的MEA1 mRNA序列,结合猪表达序列标签EST序列(CF176017、CJ008523、CK452977、FS670642、HX242707),利用Oligo 7软件设计F1/R1特异引物(F1: AGGCTAAAGCTGGGCCTAAT,R1:AACCAGGATGCGTTCAGTCC)来扩增MEA1基因665 bp,其中包括全长CDS及部分5′和3′侧翼序列;设计F2/R2特异引物(F2: CCAGGATCGAATCCAGGCCTTG,R2: CGTGTTCTGGGTCCATGGGAAT)扩增125 bp,并以GAPDH基因为内参(F: CCTTCATTGACCTCCACTACATGGT,R: CCACAACATACGTAGCACCAGCATC)扩增183 bp,利用荧光定量法分析各组织的表达水平。
1.2.3 RT-PCR分离MEA1基因和生物信息学分析
PCR反应总体系25 µL,采用TaKaRa Ex TaqE体系,琼脂糖凝胶电泳确定扩增的MEA1基因片段与引物设计时预期的片段大小吻合,胶回收后测序。利用Lasergene软件中的SeqMan程序进行组装、拼接、并人工校对,确定BMI MEA1基因编码序列并提交GenBank数据库,同时推导出其编码的氨基酸序列并进行生物信息学分析。
1.2.4 MEA1基因多组织表达分析
以50 ng/μL的睾丸cDNA为模板,5倍倍比连续稀释获得7个标准品,以F2/R2为引物,制作标准曲线。以各组织cDNA为模板,以GAPDH为内参进行多组织qPCR表达水平检测。
2. 结果与分析
2.1 MEA1基因RT-PCR级序列分析结果
以BMI睾丸组织cDNA为模板,扩增的MEA1基因长度为665 bp (图1)。将F1/R1扩增产物测序结果使用DNAStar软件中的SeqMan程序进行组装,获得BMI MEA1全长编码序列525 bp,GenBank登录号为KU705641,推测其包含174个氨基酸,对应的氨基酸登录号为AOC89059,含有1个完整的保守结构域MEA1,属于MEA1超家族(图2)。
2.2 MEA1基因组结构分析
为了获得BMI MEA1基因座的基因组DNA,利用获得的BMI MEA1 CDS序列作为种子序列,在NCBI数据库中进行BLAST搜索,找到1个含有MEA1 CDS全长的猪种混合染色体序列(Sus scrofa chromosome 7,GenBank登录号:NC_010449.5)。利用在线的ensembl分析BMI MEA1基因的外显子/内含子结构,结果表明:BMI MEA1基因长2 368 bp,含有4个外显子和3个内含子,外显子—内含子剪接均符合GT-AG规则(图3)。
2.3 MEA1蛋白质生物信息学分析结果
2.3.1 MEA1蛋白质基本信息
BMI MEA1蛋白质基本信息见表1。
表 1 MEA1蛋白质的基本信息Table 1. Primary information of MEA1 protein一级结构特性
primary structure characteristics预测结果
prediction results氨基酸数the number of amino acids 174 分子量/ku molecular weight 18.77 等电点isoelectric point 3.98 分子式formula C812H1249N219O285S4 正电荷残基 positively charge residues 8 负电荷残基 negative charge residues 38 疏水性hydropathicity −0.705 不稳定系数instability index 60.84 2.3.2 MEA1蛋白质二级结构信息
BMI MEA1蛋白质二级结构包含54个氨基酸构成的α-螺旋结构(31.03%),95个氨基酸构成的无规则卷曲结构(54.60%),14个氨基酸构成的延伸链结构(8.05%),11个氨基酸构成的β-转角结构(6.32%)。
2.3.3 MEA1蛋白质的疏水性
BMI MEA1蛋白质N端疏水,C端亲水,且在第7位置有最大疏水值1.798,在第44位置有最小疏水值−3.289。
2.3.4 MEA1蛋白质功能活性位点分析
BMI MEA1蛋白质包含N-糖基化位点、酪蛋白激酶II磷酸化位点和N-豆蔻酰化位点3种功能位点(表2)。
表 2 MEA1蛋白质3种功能活性位点Table 2. The three functional active sites of MEA1 protein活性位点名称active site 序列号sequence No. 氨基酸位置amino acids position N-糖基化位点N-glycosylation sites PS00001 19~22:NQTE;116~119:NHSS 酪蛋白激酶II磷酸化位点casein kinase II phosphorylation sites PS00006 38~41:SseE;103~106:SedE N-豆蔻酰化位点N-myristoylation site PS00008 50~55:GTgpAG 2.3.5 MEA1蛋白质N端信号肽及跨膜螺旋
MEA1蛋白质不包含N端信号肽序列,不存在跨膜螺旋,是非跨膜蛋白。
2.3.6 MEA1蛋白质序列多物种系统进化分析
利用MEGA 7.0软件将BMI MEA1的蛋白质序列与绵羊(XP_004018879)、人(NP_001305871)、长臂猿(XP_012351797)、黄牛(NP_777073)和牦牛(XP_005900181)等蛋白质序列进行了相似性比对,结果显示相似度分别为96.6%、96.5%、96.0%、94.8%和94.8%。并构建6个物种的分子系统进化树(图4)。
2.4 MEA1基因表达分析
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
图 5 MEA1基因多组织相对表达水平注:1. 心;2. 肝;3. 脾;4. 肺;5. 肾;6. 胃;7. 脑;8. 肌肉;9. 十二指肠;10. 结肠;11. 精囊腺;12. 前列腺;13. 尿道球腺;14. 睾丸;15. 附睾。Figure 5. Relative expression level of MEA1 gene in multi-tissuesNote: 1. heart; 2. liver; 3. spleen; 4. lung; 5. kidney; 6. stomach; 7. brain; 8. muscle; 9. duodenum; 10. colon; 11. seminal vesicle gland; 12. prostate; 13. urethral ball glands; 14. testis; 15. epididymis.3. 讨论
鉴于部分公猪不育已成为阻碍其群体数目继续扩大的主要瓶颈,为了尽可能繁育更多亚系,对其精子发生相关基因的研究尤为重要。近年来对于MEA1基因的结构与功能的研究越来越深入,通过SCOTT等[14]的研究发现:MEA1可能是血清学检测的雄性候选抗原之一,也可能是睾丸特异表达基因,还可能是位于白细胞抗原和组织相容性抗原之间的重叠抗原基因,但在版纳微型猪近交系的研究中尚属空白。因此,通过生物信息学分析MEA1的氨基酸构成、C端和N端氨基酸的特性,以及这些氨基酸形成的功能位点、功能域,还有卷曲折叠构成的二级结构等是很有必要的。
蛋白质的结构与功能之间存在必然联系[15],功能生物信息学分析显示BMI MEA1蛋白质二级结构中无规则卷曲和α螺旋含量较高,分别为54.60%和31.03%,无规卷曲是构成酶重要活性部位和蛋白质特异功能部位的二级结构元件,而α螺旋是蛋白质二级结构中最常见、含量也较为丰富的结构元件[15]。MEA1蛋白含有1个完整的保守结构域,属于MEA1超家族。
BMI MEA1蛋白无N端信号肽序列,表明其属于非分泌蛋白,可直接进入到细胞质基质中发挥其功能。MEA1蛋白无跨膜螺旋表明其属于非跨膜蛋白。化学修饰是蛋白质翻译后加工的重要内容,N-糖基化在蛋白质折叠、运输等过程中扮演着重要角色,蛋白质磷酸化可以大大提高酶蛋白的活性,它们是蛋白质化学修饰的两种重要方式,发现MEA1蛋白存在N-糖基化功能位点、酪蛋白激酶II磷酸化功能位点、N-豆蔻酰化位点,为今后深入进行MEA1蛋白的功能奠定了基础。
通过比对BMI猪、绵羊、人、长臂猿、黄牛和牦牛等物种的MEA1氨基酸序列,发现其相似性均在90%以上,说明该基因保守性较高。这几个物种MEA1蛋白质序列构建的系统进化树表明,牛和羊聚为一类,牛和羊均属于偶蹄目牛科动物,牛和羊是人类驯化较早的家养动物;人和长臂猿聚为一类,然后与猪聚为一类,MEA1基因编码的氨基酸与人类同源性更高,保守性更好,该基因在BMI中发挥的作用可能与人类更相似,有望对人类疾病的探索提供一定的理论依据。
实时荧光定量多组织表达谱分析表明:MEA1基因在BMI性腺睾丸中高表达,在其他组织中微量表达,由此推断该基因主要在猪的睾丸中发挥重要功能。睾丸是很重要的雄性生殖器官之一,可产生精子,分泌雄激素,产生睾丸液,维持精子的生成及存活[16]。据报道,MEA1基因已在小鼠[4]中进行扩增并对多组织的表达情况进行了实验鉴定,该基因在睾丸中明显高表达。我们研究发现MEA1基因在猪睾丸中也是高表达,其结果一致,说明该基因在睾丸组织中起着很重要的作用。
-
表 1 干旱胁迫下冬樱花叶片电解质外渗率的比较
Table 1 Comparison of electrolyte leakage rate of C. cerassoides leaves under drought stress
w (PEG-8000)/% 胁迫时间stress time 0 d 2 d 4 d 6 d 8 d 0 0.251±0.012 Aa 0.262±0.076 Aa 0.273±0.010 Aa 0.295±0.038 Aa 0.300±0.011 Aa 5 0.243±0.032 Aa 0.297±0.003 Ab 0.319±0.002 Bbc 0.343±0.008 Bcd 0.361±0.010 Bd 10 0.245±0.032 Aa 0.293±0.026 Ab 0.336±0.006 Cc 0.401±0.021 Cd 0.453±0.003 Ce 20 0.253±0.011 Aa 0.353±0.013 Bb 0.412±0.003 Dc 0.531±0.006 Dd 0.732±0.005 De 注:同列中不同大写字母表示同一测定时间不同处理间差异显著(P<0.05);同行中不同小写字母表示同一处理下不同测定时间差异显著(P<0.05);下同。
Note:Different capitals in the same columns indicate significant difference of different treatments in the same measurement time, and there is a significant difference among different treatments when P<0.05; different lowercase in the same lines indicate significant difference in the time of measurement under the same treatment(P<0.05); the same as below.表 2 干旱胁迫下冬樱花叶片丙二醛含量的比较
Table 2 Comparison of malondialdehyde content in leaves of C. cerassoides under drought stress
w (PEG-8000)/% 胁迫时间stress time 0 d 2 d 4 d 6 d 8 d 0 4.751±0.643 Aa 5.126±0.038 Aab 5.532±0.048 Abc 5.877±0.057 Acd 6.157±0.130 Ad 5 5.041±0.035 Aa 5.453±0.112 Bb 6.174±0.384 Bc 6.975±0.050 Bd 7.853±0.131 Be 10 4.805±0.097 Aa 6.143±0.057 Cb 7.374±0.032 Cc 8.132±0.039 Cd 9.561±0.095 Ce 20 4.895±0.101 Aa 7.535±0.034 Db 9.327±0.097 Dc 10.763±0.933 Dd 12.005±0.520 De 表 3 干旱胁迫下冬樱花叶片可溶性糖含量的比较
Table 3 Comparison of soluble sugar content in leaves of C. cerassoides under drought stress
w (PEG-8000)/% 胁迫时间stress time 0 d 2 d 4 d 6 d 8 d 0 2.398±0.018 Aa 2.416±0.008 Aa 2.271±0.036 Aa 2.553±0.024 Ab 2.563±0.030 Ab 5 2.403±0.010 Aa 2.563±0.018 Bb 2.669±0.030 Bc 2.874±0.029 Bd 3.154±0.033 Be 10 2.418±0.007 Aa 2.699±0.087 Cb 3.112±0.006 Cc 3.546±0.036 Cd 4.012±0.017 Ce 20 2.405±0.032 Aa 2.796±0.061 Cb 3.401±0.061 Dc 3.992±0.004 Dd 4.516±0.055 De 表 4 干旱胁迫下冬樱花叶片相对含水量的比较
Table 4 Comparison of relative moisture content of C. cerassoides leaves under drought stress
w (PEG-8000)/% 胁迫时间stress time 0 d 2 d 4 d 6 d 8 d 0 90.35±0.08 Aa 88.85±0.03 Ab 86.25±0.02 Ac 87.91±0.02 Ad 85.46±0.03 Ae 5 90.33±0.17 Aa 87.05±0.14 Bb 84.33±0.01 Bc 81.59±0.02 Bd 80.77±0.03 Be 10 90.39±0.15 Aa 85.44±0.04 Cb 81.25±0.02 Cc 74.33±0.03 Cd 62.59±0.02 Ce 20 90.36±0.01 Aa 82.79±0.01 Db 78.19±0.02 Dc 70.32±0.02 Dd 52.75±0.02 De 表 5 干旱胁迫下冬樱花叶片电解质外渗率、MDA、可溶性糖含量和相对含水量的相关性分析
Table 5 Correlation analysis of electrolyte leakage rate, malondialdehyde content, soluble sugar content and relative water content of C. cerassoides leaves under drought stress
相关系数
correlation coefficient电解质外渗率
electrolyte leakage rate丙二醛含量
malondialdehyde content可溶性糖含量
soluble sugar content叶片相对含水量
relative water content电解质外渗率
electrolyte leakage rate1.000 丙二醛含量
malondialdehyde content0.933** 1.000 可溶性糖含量
soluble sugar content0.933** 0.964** 1.000 叶片相对含水量
relative water content−0.944** −0.933** −0.970** 1.000 注:“**”表示在0.01水平显著相关。
Note:“**” means the correlation on 0.01 significance level. -
[1] 荣秀连, 王波, 刘刊, 等. PEG-6000模拟干旱胁迫对冷季型草坪种子萌发特性影响[J]. 北方园艺, 2010(8): 80 [2] 尹晗, 李耀辉. 我国西南干旱研究最新进展综述[J]. 干旱气象, 2013(1): 182 [3] 鲁松. 干旱胁迫对植物生长及其生理的影响[J]. 江苏林业科技, 2012(4): 51. DOI: 10.3969/j.issn.1001-7380.2012.04.015 [4] 王宏哲. 我国城市缺水的原因与对策[J]. 吉林建筑工程学院学报, 2001(2): 33. DOI: 10.3969/j.issn.1009-0185.2001.02.009 [5] RYAN M G. Tree responses to drought[J]. Tree Physiology, 2011, 31(3): 237. DOI: 10.1093/treephys/tpr022.
[6] 何磊. 干旱区园林绿化设计中的节水措施[J]. 绿化与生活, 2013(12): 23 [7] 洪震, 练发良, 刘术新, 等. 3种乡土园林地被植物对干旱胁迫的生理响应[J]. 浙江农林大学学报, 2016(4): 636 [8] 段晓梅, 王自辉, 樊国盛. 云南冬樱花的群落特征及地理分布[J]. 南京林业大学学报(自然科学版), 2004, 28(6): 83. DOI: 10.3969/j.issn.1000-2006.2004.06.022 [9] 关文灵. 冬日里的浪漫花木—冬樱花[J]. 云南林业, 2006(6): 27 [10] 杨明艳, 李兴明, 杨发军, 等. 冬樱花嫁接繁殖试验[J]. 农业研究与应用, 2012(2): 17. DOI: 10.3969/j.issn.2095-0764.2012.02.006 [11] 杨文良, 关文灵, 毛昆明, 等. 冬樱花不同时期播种试验[J]. 中国园艺文摘, 2010(12): 7. DOI: 10.3969/j.issn.1672-0873.2010.12.003 [12] 和凤美, 朱永平, 杨晓红, 等. 冬樱花愈伤组织诱导和抑制褐化初探[J]. 中国农学通报, 2010(12): 130 [13] 杨明艳, 李兴明, 杨发军, 等. 冬樱花扦插生根处理研究[J]. 热带农业科学, 2011(12): 20. DOI: 10.3969/j.issn.1009-2196.2011.12.006 [14] 段晓梅. 冬樱花扦插繁殖研究[J]. 西南林学院学报, 2003, 23(1): 43. DOI: 10.3969/j.issn.2095-1914.2003.01.012 [15] 谢凤云, 戴明仲, 张飞, 等. 冬樱花的引种栽培及其在园林中的应用[J]. 现代农业科技, 2017(3): 140. DOI: 10.3969/j.issn.1007-5739.2017.03.095 [16] 赵大克, 郑丽. 云南冬樱花及其在园林中的应用[J]. 云南农业大学学报(自然科学), 2009, 24(5): 778. DOI: 10.3969/j.issn.1004-390X.2009.05.027 [17] 陈剑英, 董琼, 杨超本. 冬樱花生态地理学研究[J]. 西南林学院学报, 1999, 19(1): 1. DOI: 10.3969/j.issn.2095-1914.1999.01.001 [18] 李丽, 范眸天, 杨德, 等. 昆明地区冬樱花开花生物学习性的初步研究[J]. 云南农业大学学报, 2005, 20(3): 448. DOI: 10.3969/j.issn.1004-390X.2005.03.034 [19] 孙存华, 李扬, 贺鸿雁, 等. PEG6000渗透胁迫对藜幼苗叶片渗透调节物质的影响(英文)[J]. Agricultural Science & Technology, 2007(Z1): 20. DOI: 10.3969/j.issn.1001-3601.2007.z1.006 [20] 李合生. 植物生理生化实验原理和技术[M]. 北京: 高等教育出版社, 2000. [21] 郝再彬, 苍晶, 徐仲. 植物生理实验[M]. 哈尔滨: 哈尔滨工业大学出版社, 2004. [22] 孙程旭, 刘立云, 李荣生, 等. PEG处理对蛇皮果幼苗的生理响应及酶变化研究[J]. 云南农业大学学报(自然科学), 2011, 26(5): 673. DOI: 10.3969/j.issn.1004-390X(n).2011.05.017 [23] 吴庆贵, 杨敬天, 邹利娟, 等. 珙桐幼苗生理生态特性对土壤干旱胁迫的响应[J]. 江苏农业科学, 2014, 42(2): 119. DOI: 10.3969/j.issn.1002-1302.2014.02.040 [24] 崔江慧, 李霄, 常金华. PEG模拟干旱胁迫对高粱幼苗生理特性的影响[J]. 中国农学通报, 2011, 27(9): 160 [25] 苏寒之, 金建邦, 祝遵凌. 干旱胁迫对北美红栎幼苗生理特性的影响[J]. 东北林业大学学报, 2014(8): 34. DOI: 10.3969/j.issn.1000-5382.2014.08.008 [26] 丁菲, 杨帆, 杜天真. 干旱胁迫对构树幼苗抗氧化酶活性变化的影响[J]. 江西农业大学学报, 2008, 30(4): 680. DOI: 10.3969/j.issn.1000-2286.2008.04.022 [27] 李书平. 福建山樱花和日本樱花对干旱胁迫的生理响应研究[D]. 福州: 福建农林大学, 2014. [28] 姜英淑, 陈书明, 王秋玉, 等. 干旱胁迫对2个欧李种源生理特征的影响[J]. 林业科学, 2009, 45(6): 6. DOI: 10.3321/j.issn:1001-7488.2009.06.002 [29] 张尧. 梭梭幼苗干旱死亡过程中的水力结构变化与碳分配动态[D]. 北京: 中国科学院大学, 2016. [30] 郭春芳, 孙云. 干旱胁迫下植物的渗透调节及脯氨酸代谢研究进展[J]. 福建教育学院学报, 2015, 16(1): 114. DOI: 10.3969/j.issn.1673-9884.2015.01.032 [31] 刘忠霞. 干旱胁迫下苹果树苗光合、耗水特性及抗旱性研究[D]. 咸阳: 西北农林科技大学, 2011. [32] 蒲光兰, 袁大刚, 胡学华, 等. 土壤干旱胁迫对3个杏树品种生理生化特性的影响[J]. 浙江林学院学报, 2005, 22(4): 375. DOI: 10.3969/j.issn.2095-0756.2005.04.004 [33] 胡学华. 李、桃、杏和梨四种经济树木对干旱胁迫的生理响应及其抗旱性的研究[D]. 雅安: 四川农业大学, 2005. [34] 赵福庚, 何龙飞, 罗庆云. 植物逆境生理生态学[M]. 北京: 化学工业出版社, 2004. [35] 唐承财, 钟全林, 王健. 林木抗旱生理研究进展[J]. 世界林业研究, 2008(1): 20 [36] 段胜芝, 马士祝, 牟凤娟. 裸子植物水分特征与抗旱机理研究进展[J]. 陕西林业科技, 2016(3): 83. DOI: 10.3969/j.issn.1001-2117.2016.03.022 [37] 李云飞. 土壤干旱胁迫对李属彩叶植物抗旱生理及叶色的影响[D]. 保定: 河北农业大学, 2008. [38] 孙宗国. PEG-6000模拟干旱胁迫下裸果木幼苗的生理响应[J]. 甘肃农业科技, 2014(12): 41. DOI: 10.3969/j.issn.1001-1463.2014.12.015 [39] 仇云峰, 李亚光, 李青山. 干旱胁迫对小叶杨幼苗生长的影响[J]. 水土保持通报, 2015, 35(2): 42
计量
- 文章访问数: 2062
- PDF下载量: 21