• 中国科学引文数据库(CSCD)来源期刊
  • 中文核心期刊
  • 中国农林核心期刊(A类)
  • 中国高校百佳科技期刊

版纳微型猪近交系CPN10克隆、亚细胞定位、组织表达和生物信息学分析

宋雪, 秦彩艳, 霍金龙, 王配, 李罗刚, 张霞, 王淑燕, 王雪飞

宋雪, 秦彩艳, 霍金龙, 等. 版纳微型猪近交系CPN10克隆、亚细胞定位、组织表达和生物信息学分析[J]. 云南农业大学学报(自然科学), 2019, 34(2): 255-262. DOI: 10.12101/j.issn.1004-390X(n).201803013
引用本文: 宋雪, 秦彩艳, 霍金龙, 等. 版纳微型猪近交系CPN10克隆、亚细胞定位、组织表达和生物信息学分析[J]. 云南农业大学学报(自然科学), 2019, 34(2): 255-262. DOI: 10.12101/j.issn.1004-390X(n).201803013
Xue SONG, Caiyan QIN, Jinlong HUO, et al. Cloning, Subcellular Localization, Tissue Expression and Bioinformatics Analysis of CPN10 Gene in Banna Mini-pig Inbred Line[J]. JOURNAL OF YUNNAN AGRICULTURAL UNIVERSITY(Natural Science), 2019, 34(2): 255-262. DOI: 10.12101/j.issn.1004-390X(n).201803013
Citation: Xue SONG, Caiyan QIN, Jinlong HUO, et al. Cloning, Subcellular Localization, Tissue Expression and Bioinformatics Analysis of CPN10 Gene in Banna Mini-pig Inbred Line[J]. JOURNAL OF YUNNAN AGRICULTURAL UNIVERSITY(Natural Science), 2019, 34(2): 255-262. DOI: 10.12101/j.issn.1004-390X(n).201803013

版纳微型猪近交系CPN10克隆、亚细胞定位、组织表达和生物信息学分析

基金项目: 国家自然科学基金项目(31660650,31160439,31460580,31660637)
详细信息
    作者简介:

    #对本文贡献等同,为并列第一作者。宋雪(1995—),女,山西长治人,在读硕士研究生,主要从事动物分子遗传学研究。E-mail:604855210@qq.com

    秦彩艳(1988—),女,山西临县人,硕士,主要从事动物分子遗传学研究。E-mail:907788347@qq.com

    通信作者:

    霍金龙(1975—),男,山西五寨人,博士,副教授,硕士生导师,主要从事动物分子遗传学研究。E-mail:jinlonghuo973@163.com

  • 中图分类号: S 828.2

摘要:
目的以版纳微型猪近交系(BMI)为材料克隆CPN10基因,从亚细胞、mRNA水平及生物信息学方面初步探索其结构及功能特性,为该基因在猪—人异种器官移植免疫排斥中的作用研究奠定基础。
方法采用RT-PCR方法获得CPN10基因cDNA序列,应用分子克隆技术构建目的基因和报告基因EGFP的重组真核表达载体pEGFP-C1-CPN10,将其转染猪肾细胞PK15进行瞬时表达,通过EGFP示踪CPN10在PK15细胞中的表达和定位。进一步应用荧光定量技术明确其mRNA的多组织表达特征。最后利用生物信息学分析CPN10蛋白质的结构特征。
结果得到BMI CPN10 CDS 309 bp序列(GenBank登录号:KM098149),成功构建了BMI CPN10基因的绿色荧光蛋白融合表达载体pEGFP-C1-CPN10,转染PK15细胞且主要在细胞质中检测到绿色荧光蛋白表达。多组织荧光定量表达分析表明:BMI CPN10 mRNA在皮肤中具有最高表达量,在肝和肾上腺中的表达量较高,在其余各组织中呈中低度表达或几乎不表达。生物信息学分析表明:CPN10编码102个氨基酸,其蛋白分子量为10.93 ku,等电点为8.89,含1个CPN10保守结构域和5个磷酸化位点,二级结构以延伸链结构为主,无跨膜螺旋结构和信号肽。
结论CPN10基因定位在细胞质中,在皮肤中高表达,该研究结果为进一步研究该基因在免疫排斥方面的功能奠定了理论基础。

 

Cloning, Subcellular Localization, Tissue Expression and Bioinformatics Analysis of CPN10 Gene in Banna Mini-pig Inbred Line

Abstract:
PurposeIn this study, the CPN10 gene was cloned from Banana mini-pig inbred line (BMI) and its structural and functional properties were preliminary explored from the subcellular, mRNA levels and bioinformatics, which laid the foundation for the research on the role of this gene in immune rejection of pig-human xenotransplantation.
MethodsThe cDNA sequence of BMI CPN10 was obtained by RT-PCR, and the recombinant eukaryotic expression vector pEGFP-C1-CPN10 of target gene and reporter gene EGFP was constructed by molecular cloning technology. Then, the recombinant vector was transfected transiently into porcine kidney cells (PK15) to confirm CPN10 localization in cells by tracing EGFP. Further, mRNA expression characteristics of multi-tissue were examined by qPCR techniques. Finally, structural characteristics of CPN10 protein were analyzed by bioinformatics.
ResultThe CDS sequence of BMI CPN10 was 309 bp with GenBank accession number KM098149. The fusion expression vector pEGFP-C1-CPN10 of green fluorescence protein and CPN10 gene was successfully constructed and transfected into PK15 cells. The target proteins were mainly detected in the cytoplasm by tracing the green fluorescent protein. BMI CPN10 mRNA had the highest expression in the skin, higher expression in the liver and adrenal glands, and was moderate, low or hardly expressed in the remaining tissues. Bioinformatics analysis indicated the CPN10 encoded 102 amino acids with molecular weight of 10.93 ku and isoelectric point of 8.89, containing a conserved CPN10 domain and 5 phosphorylation sites. The secondary structure of the CPN10 protein is mainly extended chain structure without transmembrane helical structure and signal peptide.
ConclusionThe CPN10 gene was located in the cytoplasm and highly expressed in the skin, which provided a theoretical basis for further study on the function of CPN10 gene in rejection of xenotransplantation.

 

  • 百脉根(Lotus corniculatus L.)是豆科百脉根属多年生牧草,又名牛角花、五叶草,是一种与紫花苜蓿相似的豆科牧草,有人称之为“瘠地苜蓿”[1]。百脉根原产欧亚大陆温带地区,北美、大洋洲、非洲北部也有分布。其粗蛋白含量高于其他豆科牧草,而粗纤维含量相对较低,是一种具有较高营养价值和栽培价值的优良牧草[2-3],且叶量丰富,产量高,再生力强,用途广,世界各国都进行了引种栽培,多年来形成了很多地区品种[4-6]。云南省草地动物科学研究院于1983年从国外引种栽培了8个品种,但由于气候、环境的差异,都未在生产中大面积推广应用[7-8]。百脉根在云南省是一个广布种,云南省野生百脉根资源丰富。但目前,针对云南野生百脉根种质资源的评价、筛选研究较少[9],因此,加大对野生百脉根种质资源的评价和筛选可以有效推进百脉根新品种的选育工作。

    本研究对云南野生资源中甸百脉根和引进品种——帝国百脉根在云南迪庆州高寒地区的区域适应性和品种生产潜力进行评价,以期为百脉根的品种选育和生产种植提供参考。

    供试野生百脉根种子于2009年10月采于云南迪庆州香格里拉县小中甸,帝国百脉根种子购于云南绿盛美地草业有限公司。

    试验地设在迪庆州香格里拉县小中甸,位于东经99°50′,北纬27°90′,海拔为3 200~3 220 m,年平均气温5.8 ℃,年均降雨量849.8 mm,无霜期120 d。土壤为褐土,有机质含量4.5%,pH 7.3。

    试验材料于2010年4月栽种,试验小区采用随机区组设计,小区面积均为5 m×3 m,组间距为1 m,区间距为50 cm,条播,播种量一般为11.25 kg/hm2,播种深度0.5 cm,行距30 cm,并在试验地四周设置保护行。2013年11月在试验地进行株高、叶长、叶宽的测定,用0.5 m×0.5 m样方选取3个重复,齐地面刈割,测定每个品种的鲜草产量。随后全部置于烘箱于65 ℃下烘24 h,回潮恒重测定鲜草水分含量、茎叶比。并采集材料带回实验室进行生理指标的测定。形态指标重复10次,生理指标重复3次。

    测定供试材料的自然植株高度、拉直植株高度、叶长和叶宽,重复10次。

    取1 m2样区材料,测定地上部分、地下部分生物量及茎叶比,重复5次。

    测定供试材料叶片相对电导率、丙二醛、可溶性蛋白含量以及超氧化化物歧化酶(SOD)和过氧化物酶(POD)活性;测定叶片和根中可溶性糖、脯氨酸含量[10-11]

    采用SPSS 16.0和Excel 2003对数据进行统计分析和作图。

    中甸百脉根和帝国百脉根植株自然高度分别为14.4、20.3 cm,拉直高度为27.8、36.2 cm,叶长分别为1.65、2.22 cm,叶宽分别为0.84、1.36 cm。从图1可见:帝国株高、叶长、叶宽显著高于中甸(P<0.05)。总体来看,中甸百脉根植株个体较矮小,叶片狭小。

    图  1  帝国与中甸百脉根的形态特征对比分析
    注:小写字母表示两个地区之间存在显著差异(P<0.05); 下同。
    Figure  1.  Comparison of morphological characteristics between Empire and Zhongdian
    Note:Different letters mean significant difference between two areas(P<0.05); the same as below.

    图2可见:中甸百脉根和帝国百脉根地下部分的鲜重分别为13.2、69.7 g,干重为3.0、20.1 g;地上部分的鲜重分别为46.7、86.6 g,干重为10.4、16.4 g,帝国百脉根单位面积生物量显著高于中甸百脉根(P<0.05)。同时,帝国百脉根茎叶比(1.22)显著低于中甸百脉根(2.43) (P<0.05),可利用营养价值更高。总体上看帝国百脉根生产性能和产量价值更高。

    图  2  帝国与中甸百脉根生物量对比分析
    Figure  2.  Comparison of biomass between Empire and Zhongdian

    图3可见:帝国百脉根叶片相对电导率、可溶性蛋白含量略高于中甸百脉根,但差异不显著(P>0.05)。而供试帝国和中甸叶片丙二醛含量分别为5.37、12.33 μmol/g,帝国显著低于中甸(P<0.05)。

    图  3  帝国与中甸百脉根相对电导率、丙二醛、可溶性蛋白含量对比分析
    Figure  3.  Comparison of relative conductivity, MDA and soluble protein content between Empire and Zhongdian

    图4可见:供试材料叶片可溶性糖含量差异不显著(P>0.05),而根中可溶性糖含量帝国显著高于中甸(P<0.05)。帝国百脉根叶片脯氨酸含量略高于中甸百脉根,但差异不显著(P>0.05)。而根中脯氨酸含量帝国(36.89 μg/g)显著低于中甸(58.40 μg/g)(P<0.05)。综合分析认为:帝国百脉根可通过降低叶片膜脂过氧化产物丙二醛含量、提高根中可溶性糖含量增强其在高寒地区的抗寒性。

    图  4  帝国与中甸百脉根叶片和根中可溶性糖与脯氨酸含量对比分析
    Figure  4.  Comparison of soluble sugar and proline content in leaves and root content between Empire and Zhongdian

    图5可知:帝国和中甸叶片超氧化物歧化酶(SOD)活性分别为993.6、1 403.4 U/(g·h),过氧化物酶(POD)活性分别为85.0、119.0 μg/(g·min),帝国均显著高于中甸(P<0.05)。

    图  5  帝国与中甸百脉根SOD、POD活性对比分析
    Figure  5.  Comparison of SOD and POD activities between Empire and Zhongdian

    株高不仅反映了品种的种群高度,而且在一定程度上反映了牧草生长能力的强弱,是衡量牧草高产的一个重要指标。崔亚飞等[12]研究表明:百脉根的自然高度影响百脉根的单株产量并且在提高单株草产量的育种工作中首先要考虑到自然高度这一影响因素。本研究结果表明:帝国百脉根在植株的自然高度、拉直高度都显著高于中甸百脉根,且其叶长和叶宽也显著大于中甸百脉根,表明帝国百脉根的生长能力和经济性状要优于中甸野生百脉根。同时,产量的测定也表明:帝国百脉根的地上部分的干重湿重、地下部分的干重湿重都显著高于中甸百脉根。所以在迪庆州香格里拉县格咱乡试验基地或者在与其相似的地域和气候环境中选择种植帝国百脉根其产量和经济性能都优于中甸百脉根。

    当植物在逆境条件下时往往会发生膜脂过氧化作用,丙二醛作为膜脂过氧化产量之一,其含量往往能反应植物细胞膜的受伤害程度。从试验结果来看中甸百脉根中丙二醛含量显著高于帝国,暗示中甸百脉根在受到相同低温胁迫的时候细胞膜受伤害程度更大。植物体内,可溶性糖可以提高细胞液浓度,降低水势,从而使冰点下降,提高植物抗寒能力[13]。本研究中帝国百脉根根中的可溶性糖含量显著高于中甸百脉根中的含量,表明帝国百脉根可通过提高其根部内含物的寒冷增强其的抗寒性。植物受到逆境胁迫时体内的脯氨酸含量会显著增加。中甸百脉根的根中的脯氨酸含量显著高于帝国,表明中甸百脉根对寒冷反应更加敏感。植物内源保护酶系超氧化物歧化酶(SOD)和过氧化物酶(POD)是植物细胞抵御活性氧伤害的重要保护酶[14]。从本研究的结果可以看出:帝国百脉根超氧化物歧化酶(SOD)和过氧化物酶(POD)的活性都显著高于中甸百脉根。表明帝国百脉根可通过提高保护酶活性来加强对活性氧伤害的抵御能力。同时,两种百脉根中相对电导率差异、可溶性蛋白含量差异不显著,说明百脉根作为冷季型牧草,对寒冷都有一定的抵御能力。

    本研究表明:帝国百脉根植株更高大,牧用经济性能更高。在迪庆州香格里拉县小中甸试验基地或者在与其相似的地域和气候环境中种植百脉根则帝国百脉根的产量和经济价值优于中甸百脉根。

    对供试材料的相对电导率、丙二醛、可溶性糖、可溶性蛋白、脯氨酸以及超氧化物歧化酶(SOD)和过氧化物歧化酶(POD)活性分析表明,帝国百脉根的抗寒性优于中甸百脉根。帝国百脉根能适应中甸地区高寒生境,且在相对寒冷的地区选择种植帝国百脉根比选择中甸百脉根更好。

    总之,在迪庆州香格里拉县小中甸试验基地或在与其相似的地域和气候环境中或者更寒冷的地区帝国百脉根的区域适应性优于中甸百脉根。

  • 图  1   CPN10基因PCR结果

    注:M. DNA相对分子质量标准DL2000;1. 引物F1/R1扩增的PCR产物。

    Figure  1.   The PCR result of CPN10

    Note: M. Marker-DL2000; 1. PCR product amplified by primers F1/R1.

    图  2   BMI CPN10真核表达PCR产物电泳结果

    注:M. DNA相对分子质量标准DL2000;1. 引物F3/R3扩增的PCR产物。

    Figure  2.   The agarose gel electrophoresis result of eukaryocyte expression PCR products of BMI CPN10

    Note: M. Marker-DL2000; 1. PCR product amplified by primers F3/R3.

    图  3   pEGFP-C1-CPN10重组质粒在PK15细胞中亚细胞定位

    注:a)绿色荧光蛋白;b) Hoechst33342对细胞核的染色;c) Mito Tracker对线粒体的染色;d)绿色荧光蛋白与细胞核的叠加;e)绿色荧光蛋白与线粒体的叠加;f)绿色荧光蛋白、细胞核和线粒体的叠加。

    Figure  3.   Subcellular localization of pEGFP-C1-CPN10 recombinant plasmids in PK15 cell

    Note: a) green fluorescent protein; b) nucleus staining by Hoechst33342; c) mitochondria staining by Mito Tracker; d) superposition of green fluorescent protein with nucleus; e) superposition of green fluorescent protein with mitochondria; f) superposition of green fluorescent protein, nucleus and mitochondria.

    图  4   BMI CPN10基因27个组织的mRNA表达谱

    Figure  4.   The mRNA expression profile 27 tissues of BMI CPN10 gene

    图  5   BMI CPN10基因编码序列及其对应的氨基酸序列

    注:ATG为起始密码子;*为终止密码子;方框为保守结构域:8~100 AA (CPN10)。

    Figure  5.   The coding sequence of BMI CPN10 gene and its corresponding amino acid sequence

    Note: ATG. start codon; *. stop codon; the conserved domain is boxed: 8-100 AA (CPN10).

    图  6   预测的BMI CPN10蛋白质保守结构域

    Figure  6.   Predicted conserved domain of BMI CPN10 protein

    图  7   BMI CPN10蛋白三级结构分析

    Figure  7.   The tertiary structure of BMI CPN10 protein

    表  1   PCR引物

    Table  1   List of PCR primers

    引物名称primer name 引物序列(5′→3′) primer sequences 产物大小/bp product size 用途usage
    CPN10-F1 CATGGCAGGACAAGCATTTA 343 序列克隆sequence cloning
    CPN10-R1 GCTTCATGTGATGCCATTAGAC
    CPN10-F2 CGTAGTAGCTGTTGGATCAGG 167 荧光定量qPCR
    CPN10-R2 CGTATTTTCCAAGAATGTCACC
    CPN10-F3 CCGGAATTCATGGCAGGACAAGCATTTAG 326 真核表达eukaryotic expression
    CPN10-R3 CGGGATCCTCAGTCTACGTATTTTCCAAG
    GAPDH-F CCTTCATTGACCTCCACTACATGGT 183 内参internal reference
    GAPDH-R CCACAACATACGTAGCACCAGCATC
    下载: 导出CSV

    表  2   预测的BMI CPN10氨基酸磷酸化位点

    Table  2   The prediction results of BMI CPN10 phosphorylated sites

    基因
    gene
    磷酸化氨基酸
    phosphorylated amino acids
    位置
    position
    分值
    score
    CPN10 丝氨酸serine 21 0.723
    64 0.988
    53 0.507
    苏氨酸threonine 79 0.563
    酪氨酸tyrosine 100 0.565
    下载: 导出CSV

    表  3   BMI CPN10功能预测

    Table  3   Prediction of CPN10 function in BMI

    功能分类functional category 概率probability 基因本体分类gene ontology category 概率probability
    氨基酸合成amino acid biosynthesis 0.011 信号转导signal transducer 0.093
    辅因子生物合成biosynthesis of cofactors 0.061 结构蛋白structural protein 0.018
    细胞被膜cell envelope 0.034 转运transporter 0.025
    细胞过程cellular processes 0.110 离子通道ion channel 0.009
    中间代谢中枢central intermediary 0.060 阳离子通过cation channel 0.010
    能量代谢energy metabolism ≥0.334 电压门控离子通道voltage gated ion channel 0.004
    脂肪酸代谢fatty acid metabolism 0.032 应激反应stress response 0.024
    嘌呤和嘧啶purines and pyrimidines 0.090 免疫反应immune response 0.052
    调节功能regulatory functions 0.019 生长因子growth factor 0.007
    复制和转录replication and transcription 0.113 转录调节transcription regulation 0.029
    翻译translation 0.116 转录transcription 0.029
    转运和结合transport and binding 0.021
    下载: 导出CSV
  • [1] 胡为民, 李幼平, 卢晓风, 等. 猪−人移植细胞性排斥反应相关候选基因pOSR1的克隆和鉴定[J]. 现代免疫学, 2004, 24(4): 268. DOI: 10.3969/j.issn.1001-2478.2004.04.002
    [2] 马学军, 侯云德. 分子伴侣基因的克隆及表达[J]. 中华实验和临床病毒学杂志, 1995, 9(2): 107
    [3]

    HU W M, CHENG J Q, LI Y P, et al. Identification of the candidate genes associated with cellular rejection in pig-to-human xenotransplantation[J]. Progress in Natural Science, 2002, 12: 710.

    [4] 胡为民. 猪到人异种移植排斥反应相关候选基因的克隆和鉴定[D]. 成都: 四川大学, 2003.
    [5] 胡为民, 程惊秋, 曾令宇, 等. 猪内皮细胞与血清作用后表达上调基因的初步鉴定[J]. 川北医学院学报, 2004, 19(1): 7. DOI: 10.3969/j.issn.1005-3697.2004.01.004
    [6]

    YU P, ZHANG L, LI S, et al. Screening and analysis of porcine endogenous retrovirus in Chinese Banna mini-pig inbred line[J]. Transplant Proceedings, 2004, 36(8): 2485. DOI: 10.1016/j.transproceed.2004.07.055.

    [7]

    ZENG R, ZENG Y Z. Molecular cloning and characterization of SLA-DR genes in the 133-family of the Banna mini-pig inbred line[J]. Animal Genetics, 2005, 36(3): 267. DOI: 10.1111/j.1365-2052.2005.01277.x.

    [8]

    HUO J L, WANG P, ZHAO Y, et al. Molecular cloning, mRNA expression and characterization of a novel FAIM1 gene from Chinese Banna mini-pig inbred line (BMI)[J]. Journal of Animal and Veterinary Advances, 2012, 11(8): 1080.

    [9] 胡为民, 程惊秋, 曾令宇, 等. 猪伴侣蛋白10基因克隆和鉴定[J]. 川北医学院学报, 2004, 19(2): 7. DOI: 10.3969/j.issn.1005-3697.2004.02.004
    [10]

    CHEN D H, LUKE K, ZHANG J J, et al. Location and flexibility of the unique C-terminal tail of Aquifex aeolicus co-chaperonin protein 10 as derived by cryo-electron micros-copy and biophysical techniques[J]. Journal of Molecular Biology, 2008, 381(3): 707. DOI: 10.1016/j.jmb.2008.06.021.

    [11]

    KULKARNI O P, RYU M, KANTNER C, et al. Recombinant chaperonin 10 suppresses cutaneous lupus and lupus nephritis in MRL-(Fas)lpr mice[J]. Nephrol Dial Transplant, 2012, 27: 1358. DOI: 10.1093/ndt/gfr544.

    [12]

    CAVANAGH A C, MORTON H. The purification of early-pregnancy factor to homogeneity from human platelets and identification as chaperonin 10[J]. European Journal of Biochemistry, 1994, 222: 551. DOI: 10.1111/j.1432-1033.1994.tb18897.x.

    [13]

    CORRAO S, CAMPANELLA C, ANZALONE R, et al. Human Hsp10 and early pregnancy factor (EPF) and their relationship and involvement in cancer and immunity: current knowledge and perspectives[J]. Life Sciences, 2010, 86: 145. DOI: 10.1016/j.lfs.2009.11.004.

    [14] 陈一曲. pIRES2-AcGFP1-CD真核表达载体的构建及其在MSCs中的表达[D]. 大连: 大连理工大学, 2007.
    [15]

    NING Z Y, PENG Y Z, HAO W B, et al. Generation of recombinant Orf virus using an enhanced green fluorescent protein reporter gene as a selectable marker[J]. BMC Veterinary Research, 2011, 7: 80. DOI: 10.1186/1746-6148-7-80.

    [16]

    ZHANG G H, GURTU V, KAIN S R. An enhanced green fluorescent protein allows sensitive detection of gene transfer in mammalian cells[J]. Biochemical and Biophysical Research Communications, 1996, 227(3): 707. DOI: 10.1006/bbrc.1996.1573.

    [17] 李凤娥, 郝峰, 藏雨轩, 等. 重组EGFP-aquaporin-4融合蛋白真核表达载体的构建及其在FRT细胞中的表达和定位[J]. 吉林大学学报(医学版), 2013, 39(2): 273. DOI: 10.7694.jidxyxb20130218.
    [18] 刘稳, 李杨, 段新源. 豆壳过氧化物酶的序列联配、二级结构预测及疏水性分析[J]. 纤维素科学与技术, 2001, 9(1): 8. DOI: 10.16561/j.cnki.xws.2001.01.002
    [19] 鞠红. 蛋白质二级结构特征分析与相互作用预测[D]. 哈尔滨: 哈尔滨工业大学, 2009.
  • 期刊类型引用(2)

    1. 刘伟,赵越,杨龙,孟翔,颜安,谢开云,崔荷婷,褚皓清. 补播优良牧草对新疆昭苏退化草地生产性能和牧草品质的影响. 草地学报. 2025(01): 231-240 . 百度学术
    2. 欧阳青,袁福锦,钟声,吴文荣,王馨,李乔仙,徐驰,李世平,高月娥,余梅,吴梦霞. 2种豆科牧草生产性能与营养价值比较研究. 草学. 2024(04): 33-37 . 百度学术

    其他类型引用(1)

图(7)  /  表(3)
计量
  • 文章访问数:  2068
  • PDF下载量:  30
  • 被引次数: 3
出版历程
  • 通信作者:  霍金龙 jinlonghuo973@163.com
  • 收稿日期:  2018-03-05
  • 修回日期:  2018-07-23
  • 网络首发日期:  2019-02-28

目录

/

返回文章
返回