下载:
下载:
-
中药丹参是唇形科鼠尾草属药用植物丹参（Salvia miltiorrhiza Bunge）的干燥根及根茎，具有活血祛瘀、通经止痛、清心除烦、凉血消痈的功效[1]。水溶性的丹酚酸类物质和脂溶性的丹参酮类物质是丹参发挥药效的重要物质基础[2-4]。丹参主要来源于人工栽培，但随着种植规模的扩大，品质退化现象严重[5]。因此，阐明丹参活性成分形成规律对保证其临床用药安全尤为重要。近年来，科研工作者从转录水平对丹参品质形成机制开展了研究，鉴定出多个参与品质调控的转录因子[6-7]。但随着研究的深入，逐渐发现酶的稳定性在药材品质形成中发挥着不可或缺的作用。
泛素-26S蛋白酶体系统通过介导蛋白泛素化降解影响靶蛋白的稳定性[8-11]。其中，F-box蛋白作为Skip-Cullin-F-box（SCF）复合物的核心组分，通过其N端的F-box基序与SKP1相互作用形成骨架，并通过C端的多种蛋白质相互作用结构域选择性识别底物蛋白，这些结构域包括WD 40重复序列、TUB 结构域、Kelch 结构域和富含亮氨酸的重复序列（LRR）等[12]。基于这些不同的结构域，F-box蛋白家族被进一步细分为多个亚家族[13-16]。近年来，F-box-LRR （FBXL）家族蛋白在植物应对生物及非生物逆境中的关键作用日益受到关注。例如，拟南芥茉莉酸受体COI1的C端为18个串联的LRRs结构域[17]，能够特异性识别转录抑制子JAZ，并通过26S蛋白酶体途径降解，解除了JAZ对茉莉素途径转录激活因子MYC2的抑制作用，从而激活茉莉素信号途径下游信号防御通路[18-20]。水稻OsCOI1是拟南芥COI1的同源基因，通过与阻遏蛋白JAZ及E3泛素连接酶复合体SCF-COI1的相互作用，诱导一系列防御及生长发育相关基因的表达[21]。拟南芥中的MAX2蛋白，其C端富含LRR结构域，通过调控气孔开合，有效防御丁香假单胞菌和胡萝卜果胶菌的入侵[22]。TIR1/AFBs生长素受体蛋白具有高度保守的F-box-LRR结构域[23]，外源病原菌侵染后，miR393过表达导致TIR1水平下降，特异性提高了拟南芥对丁香假单胞杆菌的抗性[24-25]。拟南芥F-box蛋白AFBA1通过脱落酸（ABA）信号途径正调控植物抗旱反应[26]。F-box-Nictaba的表达在拟南芥中受热胁迫、丁香假单胞菌和水杨酸（SA）的诱导, 并通过SA途径介导植物抗病性反应[27]。 拟南芥AT5G15710 基因编码的F-box蛋白参与了对重金属Cu2+或Cd2+的胁迫过程[28]。目前，F-box家族在植物中的功能研究主要集中在拟南芥、水稻等模式植物中，在丹参这一重要药用植物中的研究尚显不足。本研究从丹参基因组中筛选出104条FBXL 基因，并进行生物信息学和表达模式分析，旨在为SmFBXL的后续研究提供依据。
-
从CNCB网站（https://www.cncb.ac.cn/）中下载丹参基因组（GWHAOSJ00000000）数据及gff注释信息。从拟南芥信息网站（TAIR，https://www.arabidopsis.org/）中下载拟南芥FBXL 基因的氨基酸序列。使用TBtools v2.101中Blast Compare Two Seqs模块，将E值设置为10−5进行本地blast。此外，通过Pfam数据库（http://pfam.xfam.org/）查找丹参FBXL 基因家族特征结构域（PF00646、PF12937、PF13013），通过Simple HMM Search进行丹参FBXL 氨基酸序列的筛选。取两次筛选的序列ID的交集，最后使用在线生物信息学工具国家生物技术信息中心CDD搜索数据库（https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi）、InterPro数据库（https://www.ebi.ac.uk/interpro/）和SMART网站（https://smart.embl.de）, 去除不含F-box和LRRs的序列。
使用ProtParam工具（https://web.expasy.org/protparam/）计算了预测丹参FBXL 基因蛋白的理化特征，包括氨基酸序列长度、理论等电点（pI）、分子量（Mw）、不稳定指数、脂肪族指数和总平均亲水性指数（GRAVY）等指标。使用CELLO在线网站（http://cello.life.nctu.edu.tw/）进行丹参FBXL 基因亚细胞定位预测。
-
从PudMed（https://www.ncbi.nlm.nih.gov/）网站中下载大豆FBXL基因的氨基酸序列，连同丹参和拟南芥FBXL 基因的氨基酸序列进行系统发育关系分析。使用Clustal进行多序列比对，比对完成后使用MEGA11软件Neighbor-Joining（NJ），参数Bootstrap值设为
1000 ，进行系统发育树构建。使用iTOL（https://itol.embl.de/）对进化树进行可视化及美化。 -
使用在线MEME软件（http://www.OMI csclass. com/article/67）分析丹参FBXL 家族成员的基序结构，其中基序的最大数量设置为20，基序的最小宽度设置为6，基序的最大宽度设置为50。从丹参的基因组注释文件中提取丹参FBXL 基因的结构，包括CDS和UTR。最后，使用TBtools v2.101 进行可视化分析。
-
从丹参的基因组注释文件中提取丹参FBXL 基因染色体位置信息，并使用Tbtools软件中Gene Location Visualize from GTF/GFF模块进行可视化。BLASTP程序用于鉴定丹参中的同源FBXL 基因，e值阈值设置为< e−5。使用MCScanX默认参数分析丹参FBXL 基因之间的共线性关系。使用TBtools进行可视化。为了进行共线性分析，从Phytozome v13 网站（https://phytozome-next.jgi.doe.gov/）下载拟南芥基因组序列和注释文件，从CNCB网站南丹参基因组序列和注释文件（GWHASIU00000000）。利用MCScanX程序分析了丹参与拟南芥、南丹参之间的共线性。
-
每个SmFBXL 基因起始密码子（ATG）上游的2 000 bp区域被认为是启动子序列。使用TBtools提取启动子序列并用PlantCARE（https://bioinformatics.PSB.ugent.be/web tools/plant care/html/）预测顺式调控元件。使用TBtools中HeatMap模块对与胁迫反应、植物生长和发育、植物激素响应和光响应的顺式元件进行可视化和总结。
-
丹参植株生长于海军医大学药学系温室，为唇形科植物丹参（Salvia miltiorrhiza Bunge），于丹参生长一月后取其根、叶组织，立即液氮冷冻, 于−80℃冰箱保存，进行转录组测序。使用TBtools的HeatMap模块进行热图可视化，为避免基因表达量差异过大，影响可视化，对数据取Log2（以2为底差异表达倍数的对数）进行标准化，并进行列归一化和聚类分析。
-
经CDD和Pfam在线分析后最终获得104个SmFBXL序列，均含有F-box和LRRs保守结构域，将所得到的SmFBXL基因序列进行编号，即 SmFBXL1~SmFBXL104。通过ProtParam工具预测理化性质，发现SmFBXL基因长度范围为645 ~ 3 012 bp，编码的蛋白氨基酸长度差异亦较大，为 214~2 362 aa，平均为 494 aa；相对分子质量为 24 333~268 357，与氨基酸长度成正比；等电点为4.41~9.36，其中43个为碱性蛋白质，61个为酸性蛋白质。亚细胞定位预测结果表明，42个 SmFBXLs蛋白定位在细胞核，表明其在细胞核发挥相关调控的功能；9个 SmFBXL蛋白定位在细胞质，表明其在细胞质发挥调控的功能；3 个 SmFBXL 蛋白定位在叶绿体，推测其可能参与叶绿体的基因表达调控；此外，还有1个成员定位于线粒体，37个定位于细胞质膜，12个定位于细胞外间质。
-
为明确 SmFBXL基因家族成员间的亲缘关系及其与生物学功能的关系，采用NJ法构建丹参、拟南芥和大豆的FBXL 基因家族成员的系统发育进化树（图1）。SmFBXL 基因家族共分为8个亚族（亚族Ⅰ~Ⅷ）。其中，亚族Ⅰ的所有成员大部分来自大豆，仅有一个AtFBXL 家族成员，其他7个亚类中均有SmFBXL家族成员分布。其中，亚族Ⅱ中SmFBXL成员分布最多，有52个SmFBXL成员，其次是亚族Ⅲ、亚族Ⅳ、亚族Ⅷ，分别有17、10、8个，亚族Ⅴ、亚族Ⅶ均有6个SmFBXL成员，亚族Ⅵ仅有5个SmFBXL成员。另外，亚族Ⅱ中仅有SmFBXL和AtFBXL 家族成员，推测SmFBXL在进化过程中可能出现了功能特异蛋白。系统进化上与拟南芥直系同源基因比较相近的FBXL 基因，在其结构与功能上可能比较相似，可以作为推测丹参中FBXL 基因功能的参考依据。
图 1 丹参FBXL 基因的系统发育树
-
为了进一步分析SmFBXL基因家族成员的结构和潜在功能，利用GSDS、MEME分析SmFBXL基因的外显子与内含子的分布、保守基序。SmFBXL基因家族的蛋白序列中，基序3出现在大多数亚族中，出现频率最高，其与F-Box-LRR保守结构域相关。同一亚族中的大部分SmFBXL基序组成相似，此外，有些基序仅出现在特定亚族中，如亚族Ⅱ成员中含有独有的基序2、基序6，亚族Ⅲ特有的基序5，亚族Ⅷ特有的基序8，推测这些亚族中的FBXL基因可能具有其他重要的生物学功能（图2）。对SmFBXL基因家族成员的结构分析发现，其内含子和外显子数量存在高度的差异性，基本都含有内含子1~10个，多数SmFBXL基因含有2个内含子，同一亚族的SmFBXL转录因子的基因结构相对保守（图2）。
图 2 丹参FBXL 编码蛋白基序分析及基因结构
-
对丹参基因组中获得的SmFBXLs基因进行染色体定位分析（图3A），大部分基因分布在8条染色体的上下臂，每条染色体上的SmFBXL 基因数量为3~15个，基因在Chr01、Chr02、Chr03、Chr06、Chr07、Chr08染色体上分布的较少，集中在Chr04、Chr05染色体的较多，这可能与染色体结构和大小的差异有关。基因家族的形成与扩张很大程度上是由基因复制产生的，为了探究丹参中FBXL 家族的基因扩张情况，本研究对SmFBXL基因家族成员进行了基因复制分析。种内共线性分析（图3B）表明，丹参中FBXL 成员主要通过片段重复进行基因复制，其中涉及到片段重复的基因共有38个，占比36.5%。丹参同拟南芥、南丹参的基因组共线性分析（图3C）表明，丹参中36个FBXL 家族成员同拟南芥FBXL 具有共线性，丹参5号染色体上的基因大多对应拟南芥 3号染色体上；丹参中87个FBXL 家族成员同南丹参FBXL 具有共线性，说明丹参中FBXL 家族成员同南丹参FBXL 具有较高同源性。
图 3 丹参FBXL 基因染色体定位及共线性分析
-
位于基因上游的启动子区域中存在的顺式作用元件可以通过结合转录因子调控基因表达，并对不同环境条件做出反应，这可能为基因功能推测和分化研究提供重要依据。对 SmFBXL 启动子进行顺式作用元件预测，观察到大量不均匀分布在所有基因上的顺式作用元件，将其分为4类：应激响应、光响应、激素响应和生长发育（图4）。统计显示，应激响应类别占比最高，占所有元件的33.7%，包括10种类型的元件，前4个主要元件分别是MYC（低温胁迫）、STRE（渗透压胁迫）、ARE（抗氧化响应）和AS-1（砷胁迫）。第二大类别是光响应元件，占所有元件的26.9%，所有SmFBXL基因家族成员启动子区均含有光响应元件，其中，Box 4的频率最高（光响应元件的34.5%），其次是G-box元件（光响应元件的21.4%）。激素响应占所有元件的25.1%，包括脱落酸（ABA）响应（ABRE）、水杨酸（SA）响应（TCA-element）、茉莉酸甲酯（MeJA）响应（CGTCA-motif和TGACG-motif）、生长素（IAA）响应（TGA-element和AuxRR-core）、赤霉素（GA）响应（P-box和GARE-motif）以及植物雌激素响应（ERE）。ABRE（激素响应元件的26.1%）是最常见的激素响应元件。生长发育元件占比最低（占所有元件的14.3%），包括：负责分生组织表达的CAT-box、CCGTCC-box和CCGTCC motif、负责玉米谷蛋白代谢调控的O2-site、负责次生木质部发育的AAGAA-motif、负责胚乳表达的GCN4-motif、负责细胞周期调节的MSA-like以及负责生长发育的HD-Zip 1等。综上所述，SmFBXL 基因家族的启动子含有低温、渗透压、抗氧化、砷、光响应、SA、ABA、GA、IAA、MeJA和分生组织等响应元件，暗示SmFBXL可能在植物应激防御、激素响应和生长发育中发挥潜在作用。
图 4 丹参FBXL 基因启动子顺式作用元件及其数量分布
-
通过转录组测序，分析104个 SmFBXL 基因在不同组织中的表达模式（图5）。结果表明，104个SmFBXL在丹参不同组织中的表达水平差异很大，其中SmFBXL20、SmFBXL10、SmFBXL89、SmFBXL39、SmFBXL7、SmFBXL17、SmFBXL78、SmFBXL25、SmFBXL6、SmFBXL14、SmFBXL24、SmFBXL72、SmFBXL19 这13个基因在叶和根中表达量均相对较高；SmFBXL40、SmFBXL41 2个基因仅在叶片中表达量较高；SmFBXL44、SmFBXL1、SmFBXL91 3个基因仅在根中的表达水平较高；其中44个基因的在叶和根中的表达量都很低，FPKM 值<1。以上结果表明SmFBXL在根和叶中均能发挥功能，其中有13个SmFBXL基因表达量明显高于其他家族基因，可作为后续研究SmFBXL基因家族的候选基因，其原因可能是SmFBXL 基因功能分化所导致的。
图 5 丹参FBXL 基因表达模式分析
-
丹参酮和丹酚酸类物质是丹参中主要生物活性次生代谢产物。因此，明确丹参次生代谢产物的生物合成及其调控机制是丹参质量控制的基础。大量的研究集中在酚酸生物合成的转录调控上，但对控制酚酸生物合成的蛋白质翻译后修饰知之甚少。次生代谢物的产生受环境刺激的影响很大。当植物遭受胁迫时，相应的次生代谢物会迅速合成和积累，同时消耗大量的能量和营养[29]，这种生产过剩可能会随着逆境的消除而停止。因此，开启和关闭生物合成途径在植物的生命周期中发挥着同等重要的作用。泛素-蛋白酶体系统作为蛋白质翻译后修饰方式之一，通过介导转录因子或生物合成酶等靶蛋白的降解，在关闭生物合成途径中发挥着重要作用。例如，寒冷胁迫中，花青素可有效清除活性氧造成的氧化损伤，增强耐寒性。低温时，苹果中MdMYB23可激活MdANR的转录，促进花青素的生物合成；常温下，MdBT2介导的UPS可以降解MdMYB23，抑制MdANR，从而终止花青素的生物合成[30]，可以看到泛素化修饰在植物抗逆和次生代谢调控中发挥关键作用。
FBXL蛋白通过蛋白-蛋白相互作用域选择性募集靶蛋白，从而控制 SCF泛素连接酶复合物的特异性。FBXL 蛋白分布广泛，迄今为止，在拟南芥、水稻、大豆、玉米、苹果中分别鉴定出29、61、19、16、34 个FBXL 基因[31-33]。然而，关于丹参的FBXL 基因家族的研究还很少，有待进一步挖掘并进行功能鉴定。本研究在丹参基因组中鉴定了104个SmFBXL基因，比多数物种中的FBXL 基因数目多，其大小的差异可能是由于每个物种在进化和分化过程中全基因组多倍体化事件的发生。理化分析表明SmFBXL的氨基酸长度、相对分子质量和等电点差异较大，但基因结构较为保守，基本都含有内含子1~10个，多数SmFBXL基因含有2个内含子，类似的结果亦在水稻[31]、大豆[32]等植物中发现。保守基序对于蛋白质行使生物功能非常重要， 基序3在SmFBXL基因家族出现频率最高，其与F-Box结构域相关，有些基序仅出现在特定亚族中，推测这些亚族中的 FBXL 基因可能具有其他重要的生物学功能。这一现象与油菜[34]的结果相似。SmFBXL 蛋白具有较高的基序保守性，同一亚族的SmFBXL蛋白基序组成基本相似，但也有少数聚在一支的SmFBXL 蛋白基序构成并不相同，推测其可能与SmFBXL 家族成员的功能分化有关。
本研究采用NJ法构建了丹参、拟南芥和大豆的FBXL 基因家族成员的系统发育进化树，将丹参FBXL 基因家族分为了7个亚族。蛋白序列的相似性往往决定其功能的相似性，研究发现拟南芥FBXL 蛋白AT2G42620.1通过调控气孔开合，有效防御丁香假单胞菌和胡萝卜果胶菌的入侵[35]，同源进化分析发现，SmFBXL36与AT2G42620.1的同源关系较近，组织表达分析发现在叶中表达量较高，启动子分析发现含有低温响应元件LTR，推测可能与其同源蛋白具有相似的生物学功能。拟南芥AT3G60350.1和AT2G44900.1编码蛋白可以促进侧根发育，与野生型幼苗相比，突变株形成的侧根较少，而高表达株产生的侧根较多[36-37]。系统进化分析发现SmFBXL86、SmFBXL79 与AT3G60350.1、AT2G44900.1同源关系较近，SmFBXL79含有分生组织表达元件CCGTCC 基序和CCGTCC-box，暗示其可能在调控丹参侧根生长中起重要作用。AT5G01720.1（又称作RAE1）与SmFBXL82同属一个分支，是 SCF 型 E3 连接酶复合物的 F-box 蛋白成分。它是明矾诱导调节环路的一部分：STOP1 会诱导 RAE1 的活性，反过来RAE1又会泛素化 STOP1，然后将其降解[38]。AT4G33210编码 SLOMO（SLOW MOTION）蛋白，可以平衡生长素和启动嫩枝分生组织，用于调节温度介导的下胚轴生长[39]， SmFBXL40、SmFBXL11和AT4G33210同属一个分支，并且均含有生长素响应元件TGA-element，推测可能具有生长素调控功能。AT3G54650.1和SmFBXL15同源关系较近，可能参与细胞周期基因调控[40]。AT4G07400.1（VFB3）和SmFBXL5、SmFBXL32同源关系较近，与 VFB1、VFB2 和 VFB4 似乎存在功能冗余，当所有 4 个基因的水平都降低时，植株会表现出生长缺陷[41]。综上所述，结合启动子顺式作用元件和组织表达分析，推测SmFBXL家族在丹参逆境胁迫和生长发育中可能发挥重要作用。
基因表达是其发挥生物学功能的前提。SmFBXL基因呈现不同的组织表达特征，104个在丹参不同组织中的表达水平差异很大，部分家族基因仅在叶或根中表达，44个基因在根和叶中表达量都很低。SmFBXL20、SmFBXL10、SmFBXL89、SmFBXL39、SmFBXL7、SmFBXL17、SmFBXL78、SmFBXL25、SmFBXL6、SmFBXL14、SmFBXL24、SmFBXL72、SmFBXL19 这13个基因在叶和根中表达量均相对较高，表明这些SmFBXL可能在不同组织的转录调节中起协助作用，可能是SmFBXL基因功能分化所导致的。部分基因在丹参根中高水平表达，提示这些SmFBXL可能在丹参根发育中发挥重要作用。
当前研究关于SmFBXL基因在丹参中的作用尚待深入探索。该基因家族是否能增强丹参的抗逆性，诱导丹参酮及酚酸生物合成关键酶基因表达，以及其具体作用机制尚不明确。本研究对丹参FBXL家族进行了系统性鉴定，为中药模式植物丹参中FBXL介导的泛素化修饰的生物合成与调控机制解析提供了数据支持，有望为提升丹参的抗逆性及其品质改良提供分子育种上的新策略与参考依据。未来的研究将聚焦于这些核心问题，以期推动丹参种质资源的优化与发展。
Identification and expression pattern analysis of FBXL gene family in Salvia miltiorrhiza
-
摘要:
目的 基于基因组数据鉴定丹参F-box-LRR （FBXL ）基因家族，并对其进行生物信息学与表达模式分析，为进一步深入阐明其基因功能提供依据。 方法 从丹参基因组数据库中鉴定出SmFBXL基因，运用生物信息学方法及在线工具分析其基因结构特征，启动子顺式作用元件，编码的蛋白理化性质、系统进化、组织表达等。 结果 从丹参的基因组中共鉴定出104个SmFBXL基因（SmFBXL1~SmFBXL104），不均等分布于8条染色体上，上游启动子含有与植物抗逆、生长发育和激素应答等相关的顺式作用元件。构建丹参、拟南芥和大豆的FBXL 家族成员的系统发育树，将104个SmFBXL基因分为7个亚族。通过同源进化分析，猜测SmFBXL36可能参与防御病原菌入侵，SmFBXL86、SmFBXL79可能在调控丹参侧根生长中起重要作用，SmFBXL11、SmFBXL40可能调节下胚轴生长。转录组数据显示SmFBXL基因在丹参不同组织中差异表达，其中13个SmFBXL基因在根和叶中的表达水平较高，可作为后续研究SmFBXL基因家族的候选基因。 结论 研究结果为进一步解析SmFBXL基因在丹参逆境响应及次生代谢产物生物合成中的调控机制提供了参考。 -
关键词:
- 丹参 /
- F-box-LRR基因家族 /
- 生物信息学 /
- 功能分析
Abstract:Objective To identify and analyze the bioinformatics and expression patterns of the F-box-LRR（FBXL） gene family of Salvia miltiorrhiza based on genomic data, and provide a foundation for further elucidating its gene functions. Methods The SmFBXL gene was identified from the Salvia miltiorrhiza genomic database. Its gene structure features, promoter cis-acting elements, physicochemical properties of encoded proteins, evolutionary relationships, and tissue expression were analyzed by bioinformatics methods and online tools. Results A total of 104 SmFBXL genes were identified from the Salvia miltiorrhiza genome, unevenly distributed on 8 chromosomes, with upstream promoters containing cis-acting elements related to plant stress resistance, growth and development, and hormone response. A phylogenetic tree of the FBXL family members of Salvia miltiorrhiza, Arabidopsis thaliana, and Glycine max was constructed, dividing the 104 SmFBXL genes into 7 subfamilies. Through homologous evolution analysis, it was speculated that SmFBXL36 might be involved in defense against pathogen invasion, SmFBXL86 and SmFBXL79 might play important roles in regulating lateral root growth in Salvia miltiorrhiza, and SmFBXL11 and SmFBXL40 might regulate hypocotyl growth. Transcriptome data showed differential expression of SmFBXL genes in different tissues of Salvia miltiorrhiza, with 13 SmFBXL genes showing higher expression levels in roots and leaves, serving as candidate genes for further research on the SmFBXL gene family. Conclusion The research results provided a reference for further elucidating the regulatory mechanisms of SmFBXL genes in stress response and secondary metabolite biosynthesis in Salvia miltiorrhiza. -
Key words:
- Salvia miltiorrhiza /
- F-box-LRR gene family /
- Bioinformatics /
- Functional analysis
-
METRNL（Meteorin-like）是一个新发现的分泌蛋白，为神经营养调节因子Meteorin的同源蛋白。2014年，本实验室首次报道METRNL是一个新的脂肪因子，由于其在皮下白色脂肪组织中表达很丰富，故也称为Subfatin[1]。10年来，我们对METRNL的功能进行了多方面的探索，并不断扩展METRNL的研究工具与平台。我们的研究已发现，该蛋白参与调节机体多种病理生理过程，比如：脂肪细胞METRNL可促进白色脂肪分化、脂质代谢并抑制脂肪炎症，从而抵抗高脂饮食诱导的胰岛素抵抗[2]；肠上皮细胞METRNL参与调节肠道抗菌肽的平衡[3]，且肠道METRNL缺乏会加重溃疡性结肠炎[4]；此外，METRNL促进小鼠皮肤创伤愈合[5]，并能对抗D-半乳糖诱导的衰老小鼠的认知功能障碍[6]。最近，我们新报道了血液METRNL的主要分泌来源是血管内皮细胞，并发现内皮细胞METRNL对维持血管内皮正常功能和对抗动脉粥样硬化具有重要作用[7]。在这些研究中，我们构建METRNL基因的全身性和各种组织特异性的敲除小鼠，以及多种双基因敲除小鼠，并在体外细胞实验中充分利用METRNL重组蛋白探索相关治疗学意义。除了本实验室，全球其他多个实验室也展开了对METRNL的功能探索，并发现METRNL在能量代谢[8-9]、炎症[10-11]、心脏疾病[12-13]等多种病理生理过程中发挥积极作用。
尽管目前有很多关于METRNL的研究结果提示，该蛋白具有非常好的临床治疗潜力，但有关整体METRNL治疗学探索不多，尤其是长期治疗研究几乎没有。主要原因之一是市场METRNL重组蛋白价格昂贵，而且我们前期研究结果发现：对C57BL/6J 小鼠单次静脉注射1.75 µg METRNL重组蛋白后，血清METRNL在15 min后急剧升高（226 ng/ml），接着在4 h内迅速下降约90%。虽然在注射后24 h仍明显高于基础水平，但此时血中METRNL浓度已下降约97%[2]。因此，以重组蛋白给药方式在动物整体水平研究METRNL的治疗学作用，尤其是长期治疗学作用将产生巨大经济成本。另一方面，对于一些已经体现METRNL治疗潜力的疾病（如动脉粥样硬化），疾病发展缓慢，短期给予METRNL重组蛋白很难起到治疗作用。
因此，本研究旨在构建一株长期稳定高表达METRNL的小鼠作为METRNL的治疗学研究工具，并对该小鼠高表达METRNL的情况进行验证。
1. 材料与方法
1.1 实验试剂和仪器
鼠尾DNA提取试剂盒（CW2094S）购自北京康伟试剂生物科技有限公司；5 × PrimeScript RT Master Mix（Takara 公司）；小鼠Tubulin抗体（ AT819，碧云天公司）；通用型RNA提取试剂盒Ⅱ（AG21022，艾瑞克生物科技）；Human Meteorin-like/METRNL DuoSet ELISA试剂盒（DY7867-05，R&D system公司）；山羊抗兔 IgG（ab175471）、山羊抗小鼠 IgG（ab216772）、抗METRNL抗体（ab235775）购自Abcam公司。
LightCycler96实时荧光定量PCR仪（Roche公司）；TP600PCR仪（Takara公司）；5200S化学发光分析系统（Tanon公司）。
1.2 实验动物
人METRNL基因条件性过表达（R26-LSL-METRNL+/-）小鼠为实验室前期构建所得。SPF级8周龄 C57BL/6J 小鼠和Dppa-Cre小鼠购自上海南方模式生物技术有限公司。
所有实验小鼠均饲养在独立通气笼盒（IVC）系统中，温度（24±2）℃，相对湿度为40%~60%，饲养期间笼盒内保持清洁，小鼠在笼内自由活动、进食及饮水，动物房内照明系统为自动控制（12 h照明、12 h黑暗）。动物实验标准均依照国家《实验动物护理使用卫生指南》，并经过海军军医大学医学研究伦理委员会批准指导。
1.3 实验方法
1.3.1 动物基因型鉴定
将剪刀消毒后剪取小鼠尾尖约3 mm，剪碎，按照DNA提取试剂盒的说明书方法提取DNA之后，对目的基因进行PCR扩增，各引物序列见表1。
表 1 PCR扩增实验中的引物序列基因名称 上游引物（5’→3’） 下游引物（5’→3’） R26-WT TCAGATTCTTTTATAGGGGACACA TAAAGGCCACTCAATGCTCACTAA R26-L-METRNL AAAGTCCCGGAAAGGAGCTG GAGGCTCCATCCAGCAAGTT R26-Stop GGGCAACGTGCTGGTTATTG ACTTGCCCCTTGCTCCATAC 内参基因 TGGGTTGGGTGTCTGTTTCATTGT GATCCACCTGTCTCTGCCTTCC Dppa-Cre TGGGTTGGGTGTCTGTTTCATTGT GACCTTGCATTCCTTTGGCGAGAG 将PCR产物进行1.2% 琼脂糖凝胶电泳，上样量为每孔6 μl，电泳条件为100 V，30 min，结束后进行拍照、分析。
1.3.2 小鼠取材
将小鼠称重后，腹腔注射1%戊巴比妥钠溶液（100 mg/kg），待小鼠处于深度麻醉后，打开其胸腔，自上下腔静脉汇合处缓慢抽取血液，并转移至1.5 ml EP管静置于室温。迅速剪取小鼠心、肝、脾、肺、肾、脑、白色脂肪和肌肉组织，放入组织冻存管扔进液氮速冻，待取材结束后及时转入−80 ℃超低温冰箱储存。血液于室温静置2 h后离心：4 ℃，3000×g，15 min，分离血清，储存至−80 ℃超低温冰箱。
1.3.3 实时荧光定量PCR实验
使用RNA提取试剂盒提取组织RNA，将得到的RNA进行浓度测定与吸光度测定后进行逆转录，得到cDNA用于实时荧光定量PCR实验，各引物序列见表2。
表 2 实时荧光定量PCR实验中的引物序列基因名称 上游引物（5′→3′） 下游引物（5′→3′） 人 METRNL ACCAGCGACTTCGTAATTCAC CAGCTCCACGTCATGGGTG 小鼠 Gapdh GTATGACTCCACTCACGGCAAA GGTCTCGCTCCTGGAAGATG 1.3.4 蛋白免疫印迹实验
取适量组织至2 ml高速离心管，加入蛋白裂解液后，使用高通量匀浆仪匀浆240 s。取出高速离心管，离心：12000 × g，20 min。将上清液转移至另一干净1.5 ml EP管中，进行蛋白浓度测定，剩余样品加入5 × 蛋白上样缓冲液，97 ℃变性10 min得到蛋白样品。
使用10% SDS-PAGE凝胶进行电泳，电泳条件为：150 V，60 min。使用PVDF膜进行转膜，转膜条件为100 V，60 min。
转膜结束后，使用快速封闭液封闭15 min，之后使用1×TBST缓冲液洗膜，5 min × 3次。加入一抗（1∶1000稀释）4 ℃孵育过夜。次日去除一抗孵育液，使用1 × TBST缓冲液洗膜，5 min × 3次。加入二抗孵育液（1∶2 000稀释）常温孵育1 h，用1 × TBST缓冲液洗去二抗，5 min × 4次，结束后即可进行扫膜。
1.3.5 血清酶联免疫吸附实验
使用酶联免疫吸附实验试剂盒（DY7867-05）测定小鼠血清中的METRNL水平，具体操作按照试剂盒说明书进行。
1.3.6 统计学分析
实验数据使用GraphPad Prism 8.0进行统计分析。两组间的比较使用双尾t检验，P<0.05视为差异有统计学意义。
2. 结果
2.1 全身过表达人METRNL基因小鼠的构建和基因型鉴定
本研究基于前期构建好的人METRNL基因条件性过表达小鼠（简称为R26-LSL-METRNL+/-小鼠）和Cre-LoxP技术，最终获得全身过表达人METRNL基因小鼠（简称为R26-L-METRNL+/-小鼠），具体构建策略如图1所示。
R26-LSL-METRNL+/-小鼠是前期通过CRISPR/Cas9基因编辑技术构建而成，即在其中一个Rosa26基因位点定点插入了CAG-LoxP-Stop-LoxP-METRNL-3XFlag-Wpre-pA表达框，且该表达框的终止密码子Stop两侧插有同向LoxP位点，可基于Cre-loxP系统在Cre酶的作用下，将LoxP位点之间的序列切除，只留下一个LoxP位点，最终达到人METRNL基因过表达的目的。在该小鼠的名称“R26-LSL-METRNL+/-”中，“+”表示有外源基因表达框的插入，“-”表示无外源基因表达框插入。
Dppa3-Cre小鼠是由Dppa3基因启动子介导Cre重组酶在全身表达的工具鼠，将R26-LSL-METRNL+/-小鼠和Dppa3-Cre小鼠杂交，可获得R26-L-METRNL+/-小鼠，具体繁殖方法如图2所示。其中，野生对照小鼠简写为R26-WT，表示在Rosa26位点没有外源人METRNL基因表达框的插入。
在繁殖过程中进行基因型鉴定时，需确认外源人METRNL基因表达框、终止密码子Stop和Dppa-Cre基因的存在情况。采用相应基因上下游引物分别进行鼠尾基因型鉴定，外源性表达框阳性条带为699 bp，对应野生型序列条带为996 bp，终止密码子Stop阳性条带为408 bp，Cre基因阳性条带为100 bp。如图3所示，泳道1为R26-LSL-METRNL+/-小鼠，泳道2为R26-L-METRNL+/-小鼠，泳道3为R26-WT小鼠，泳道4为R26-L-METRNL+/-Cre小鼠，泳道5为R26-WT;Cre小鼠。
2.2 R26-L-METRNL+/-小鼠各组织人METRNL mRNA表达
为验证R26-L-METRNL+/-小鼠是否存在人METRNL基因过表达，本研究首先利用实时荧光定量PCR技术检测了该小鼠各组织中人METRNL mRNA的表达情况。如图4所示，以 R26-WT 小鼠白色脂肪的人 METRNL mRNA表达量为1，肝、脾、肺、肾、白色脂肪和脑组织的相对表达量分别为94 008.6、618.1、88 537.3、68 897.9、32 386.3和24 816.5。该结果说明，R26-L-METRNL+/-小鼠在组织mRNA水平上实现了人METRNL基因过表达。
2.3 R26-L-METRNL+/-小鼠各组织人METRNL蛋白表达
接着，本研究提取各组织的蛋白，使用蛋白免疫印迹实验方法验证R26-L-METRNL+/-小鼠各组织中人METRNL蛋白表达情况。在该实验中，所用METRNL抗体可同时抗人和小鼠的METRNL蛋白。如图5所示，METRNL蛋白在R26-L-METRNL+/-小鼠的心、肝、脾、肺和肌肉组织中的含量明显高于R26-WT小鼠，而在肾组织中METRNL抗体的结合效果不佳，且并未发现两组小鼠肾METRNL蛋白存在明显差异。该结果提示，R26-L-METRNL+/-小鼠在组织蛋白水平上实现了人METRNL基因过表达。
图 5 R26-L-METRNL+/-小鼠各组织的METRNL 蛋白表达情况A.心脏组织；B.肝组织；C.脾组织；D.肺组织；E.肾组织；F.肌肉组织 （n=3，$ \bar x \pm s $ 2.4 R26-L-METRNL+/-小鼠血清中人METRNL蛋白水平
由于METRNL为分泌性蛋白，本研究最后使用ELISA方法测定R26-L-METRNL+/-小鼠及其对照小鼠血清中人METRNL蛋白水平。如图6所示，可成功检测到R26-L-METRNL+/-小鼠血液中的人METRNL，浓度在1 100.3～1 579.3 pg/ml，而在R26-WT小鼠血液中检测不到人METRNL。该结果说明，R26-L-METRNL+/-小鼠的血液中存在大量人METRNL蛋白，提示R26-L-METRNL+/-小鼠实现了人METRNL基因过表达。
图 6 R26-L-METRNL+/-小鼠及其对照小鼠的血清METRNL蛋白水平（n=5，$ \bar x \pm s $ ）A.“METRNL浓度-吸光度”标准曲线（R2=0.992 9）；B.R26-L-METRNL +/-小鼠血清METRNL浓度 *** P <0.001，与 R26-WT 组比较。3. 讨论和总结
本研究利用Dppa-Cre小鼠和实验室前期构建的R26-LSL-METRNL+/-小鼠进行杂交繁殖，最终获得R26-L-METRNL+/-小鼠。该小鼠与R26-LSL-METRNL+/-小鼠都在Rosa26基因位点含有外源性基因表达框，不同的是在R26-L-METRNL+/-小鼠外源表达框中的终止密码子Stop被Cre酶成功切除，因此R26-L-METRNL+/-小鼠可实现全身细胞过表达人METRNL。
本研究从mRNA水平、组织蛋白水平和血清蛋白水平考察了R26-L-METRNL+/-小鼠过表达METRNL的情况。由于该小鼠插入的外源METRNL基因是人METRNL基因，因此在进行实时荧光定量PCR实验时，使用的是人METRNL引物和小鼠Gapdh引物。类似地，在ELISA实验中，使用的是检测人METRNL的试剂盒。由于没有特异性抗人的METRNL抗体，因此在蛋白免疫印迹实验中使用的是可同时抗人和小鼠的METRNL抗体，而内参使用的是抗鼠的Tubulin抗体。本研究结果显示，R26-L-METRNL+/-小鼠的各组织存在人METRNL mRNA高表达，虽然各组织的表达量有所波动，但都比R26-L-WT小鼠的表达量高几千倍甚至是数十万倍。蛋白免疫印迹实验结果显示，在R26-L-METRNL+/-小鼠的心、肝、脾、肺、肌肉组织存在明显升高的METRNL蛋白水平，但该实验中检测的METRNL蛋白量升高倍数不多，可能与使用的METRNL抗体可以同时抗人和小鼠的METRNL有关。另外，我们也注意到两组小鼠肾组织的METRNL蛋白水平并没有明显差异，并且两组条带都非常微弱，为确证实验结果，我们进行了重复实验，得出相同结果。经过分析可能是因为该METRNL抗体对肾组织蛋白的亲和力不是很高，导致检测出的蛋白绝对量都太低而使两组之间很难出现差异。在血清水平，本研究结果提示R26-L-METRNL+/-小鼠的血液中存在大量人METRNL蛋白，而R26-L-WT小鼠血中检测不到该蛋白，说明R26-L-METRNL+/-小鼠成功过表达人METRNL，并可以成功分泌至血液中。同时，该实验提示，本研究所使用的人METRNL ELISA试剂盒特异性比较好，可以清晰区别人和鼠来源的METRNL蛋白。
在小鼠培育过程当中，尚未发现R26-L-METRNL+/-小鼠和同窝对照WT小鼠在体重、形态等方面有何差异。根据图2中R26-L-METRNL+/-小鼠的培育方式，采用R26-L-METRNL+/-小鼠与C57BL/6J 小鼠杂交进行扩大繁殖和保种，基于孟德尔遗传定律，后代鼠中R26-L-METRNL+/-小鼠理论得率为50%，但实际中，我们发现该小鼠的得率仅为15%左右，远低于理论值。这一现象非常有趣，因为此前本实验室构建过多种基因工程动物模型，均没有发现类似偏离孟德尔遗传定律的现象，而这一现象是否与METRNL蛋白的全身性过表达有关，值得我们进一步研究。
-
图 3 丹参FBXL 基因染色体定位及共线性分析
A. 丹参FBXL 基因染色体定位；B. 丹参FBXL 基因的共线性分析：FBXL 家族片段重复间用红色连接线连接；C. 丹参、拟南芥与南丹参物种间FBXL 家族基因共线性分析：橘黄色代表丹参染色体，绿色代表拟南芥染色体，蓝色代表南丹参染色体，红色连接线代表有不同物种间FBXL 基因存在的共线性
-
[1] 马晓晶, 杨健, 马桂荣, 等. 中药丹参的现代化研究进展[J]. 中国中药杂志, 2022, 47(19):5131-5139. [2] 陈雨萌, 李倩, 刘维海, 等. 丹参活性成分治疗心血管疾病的药理作用、临床应用及不良反应研究进展[J]. 药学研究, 2023, 42(12):1028-1034. [3] 焦育强, 刘文斌, 袁夏, 等. 丹参素及其衍生物心血管作用机制的研究进展[J]. 药学实践杂志, 2015, 33(5):389-391,405. [4] 温萍, 张俊平. 隐丹参酮及其衍生物抗肿瘤活性研究进展[J]. 药学实践与服务, 2023, 41(4):207-211. [5] 杨彬, 赵文博, 张海燕, 等. 丹参资源的遗传多样性及其保护利用[J]. 寒旱农业科学, 2023, 2(11):1002-1008. doi: 10.3969/j.issn.2097-2172.2023.11.004 [6] YANG N, ZHOU W, SU J, et al. Overexpression of SmMYC2 Increases the Production of Phenolic Acids in Salvia miltiorrhiza[J]. Front Plant Sci, 2017, 8:1804. doi: 10.3389/fpls.2017.01804 [7] DENG C, HAO X, SHI M, et al. Tanshinone production could be increased by the expression of SmWRKY2 in Salvia miltiorrhiza hairy roots[J]. Plant Sci, 2019, 284:1-8. doi: 10.1016/j.plantsci.2019.03.007 [8] SADANANDOM A, BAILEY M, EWAN R, et al. The ubiquitin-proteasome system: central modifier of plant signalling[J]. New Phytol, 2012, 196(1):13-28. doi: 10.1111/j.1469-8137.2012.04266.x [9] BUCKLEY D L, CREWS C M. Small-molecule control of intracellular protein levels through modulation of the ubiquitin proteasome system[J]. Angew Chem Int Ed Engl, 2014, 53(9):2312-2330. doi: 10.1002/anie.201307761 [10] SANG Y, YAN F, REN X. The role and mechanism of CRL4 E3 ubiquitin ligase in cancer and its potential therapy implications[J]. Oncotarget, 2015, 6(40):42590-42602. doi: 10.18632/oncotarget.6052 [11] XU G, MA H, NEI M, et al. Evolution of F-box genes in plants: different modes of sequence divergence and their relationships with functional diversification[J]. Proc Natl Acad Sci U S A, 2009, 106(3):835-840. doi: 10.1073/pnas.0812043106 [12] LECHNER E, ACHARD P, VANSIRI A, et al. F-box proteins everywhere[J]. Curr Opin Plant Biol, 2006, 9(6):631-638. doi: 10.1016/j.pbi.2006.09.003 [13] CUI H R, ZHANG Z R, LV W, et al. Genome-wide characterization and analysis of F-box protein-encoding genes in the Malus domestica genome[J]. Mol Genet Genomics, 2015, 290(4):1435-1446. doi: 10.1007/s00438-015-1004-z [14] JIA F, WU B, LI H, et al. Genome-wide identification and characterisation of F-box family in maize[J]. Mol Genet Genomics, 2013, 288(11):559-577. doi: 10.1007/s00438-013-0769-1 [15] GUPTA S, GARG V, KANT C, et al. Genome-wide survey and expression analysis of F-box genes in chickpea[J]. BMC Genomics, 2015, 16(1):67. doi: 10.1186/s12864-015-1293-y [16] WANG G M, YIN H, QIAO X, et al. F-box genes: Genome-wide expansion, evolution and their contribution to pollen growth in pear(Pyrus bretschneideri)[J]. Plant Sci, 2016, 253:164-175. doi: 10.1016/j.plantsci.2016.09.009 [17] YAN J, ZHANG C, GU M, et al. The Arabidopsis CORONATINE INSENSITIVE1 protein is a jasmonate receptor[J]. Plant Cell, 2009, 21(8):2220-2236. doi: 10.1105/tpc.109.065730 [18] SHEARD L B, TAN X, MAO H, et al. Jasmonate perception by inositol-phosphate-potentiated COI1-JAZ co-receptor[J]. Nature, 2010, 468(7322):400-405. doi: 10.1038/nature09430 [19] CHINI A, BOTER M, SOLANO R. Plant oxylipins: COI1/JAZs/MYC2 as the core jasmonic acid-signalling module[J]. FEBS J, 2009, 276(17):4682-4692. doi: 10.1111/j.1742-4658.2009.07194.x [20] ZHANG C, LEI Y, LU C, et al. MYC2, MYC3, and MYC4 function additively in wounding-induced jasmonic acid biosynthesis and catabolism[J]. J Integr Plant Biol, 2020, 62(8):1159-1175. doi: 10.1111/jipb.12902 [21] LEE S H, SAKURABA Y, LEE T, et al. Mutation of Oryza sativa CORONATINE INSENSITIVE 1b(OsCOI1b)delays leaf senescence[J]. J Integr Plant Biol, 2015, 57(6):562-576. doi: 10.1111/jipb.12276 [22] PIISILA M, KECELI M A, BRADER G, et al. The F-box protein MAX2 contributes to resistance to bacterial phytopathogens in Arabidopsis thaliana[J]. BMC Plant Biol, 2015, 15:53. doi: 10.1186/s12870-015-0434-4 [23] PARRY G, CALDERON-VILLALOBOS L I, PRIGGE M, et al. Complex regulation of the TIR1/AFB family of auxin receptors[J]. Proc Natl Acad Sci U S A, 2009, 106(52):22540-22545. doi: 10.1073/pnas.0911967106 [24] NAVARRO L, DUNOYER P, JAY F, et al. A plant miRNA contributes to antibacterial resistance by repressing auxin signaling[J]. Science, 2006, 312(5772):436-439. doi: 10.1126/science.1126088 [25] SI-AMMOUR A, WINDELS D, ARN-BOULDOIRES E, et al. miR393 and secondary siRNAs regulate expression of the TIR1/AFB2 auxin receptor clade and auxin-related development of Arabidopsis leaves[J]. Plant Physiol, 2011, 157(2):683-691. doi: 10.1104/pp.111.180083 [26] KIM Y Y, CUI M H, NOH M S, et al. The FBA motif-containing protein AFBA1 acts as a novel positive regulator of ABA response in Arabidopsis[J]. Plant Cell Physiol, 2017, 58(3):574-586. doi: 10.1093/pcp/pcx003 [27] STEFANOWICZ K, LANNOO N, ZHAO Y, et al. Glycan-binding F-box protein from Arabidopsis thaliana protects plants from Pseudomonas syringae infection[J]. BMC Plant Biol, 2016, 16(1):213. doi: 10.1186/s12870-016-0905-2 [28] REMANS T, SMEETS K, OPDENAKKER K, et al. Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations[J]. Planta, 2008, 227(6):1343-1349. doi: 10.1007/s00425-008-0706-4 [29] WATERMAN P G. Roles for secondary metabolites in plants[J]. Ciba Found Symp, 1992, 171:255-269. [30] AN J P, LI R, QU F J, et al. R2R3-MYB transcription factor MdMYB23 is involved in the cold tolerance and proanthocyanidin accumulation in apple[J]. Plant J, 2018, 96(3):562-577. doi: 10.1111/tpj.14050 [31] JAIN M, NIJHAWAN A, ARORA R, et al. F-box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress[J]. Plant Physiol, 2007, 143(4):1467-1483. doi: 10.1104/pp.106.091900 [32] JIA Q, XIAO Z X, WONG F L, et al. Genome-wide analyses of the soybean f-box gene family in response to salt stress[J]. Int J Mol Sci, 2017, 18(4):818. [33] KURODA H, TAKAHASHI N, SHIMADA H, et al. Classification and expression analysis of Arabidopsis F-box-containing protein genes[J]. Plant Cell Physiol, 2002, 43(10):1073-1085. doi: 10.1093/pcp/pcf151 [34] 左蓉, 吴姗, 刘杰, 等. 油菜F-box-LRR基因全基因组鉴定与核盘菌诱导应答分析[J]. 中国油料作物学报, 2022, 44(3):503-514. [35] CHANG W, QIAO Q, LI Q, et al. Non-transcriptional regulatory activity of SMAX1 and SMXL2 mediates karrikin-regulated seedling response to red light in Arabidopsis[J]. Mol Plant, 2024, 17(7):1054-1072. doi: 10.1016/j.molp.2024.05.007 [36] NIBAU C, GIBBS D J, BUNTING K A, et al. ARABIDILLO proteins have a novel and conserved domain structure important for the regulation of their stability[J]. Plant Mol Biol, 2011, 75(1-2):77-92. doi: 10.1007/s11103-010-9709-1 [37] GIBBS D J, VOSS U, HARDING S A, et al. AtMYB93 is a novel negative regulator of lateral root development in Arabidopsis[J]. New Phytol, 2014, 203(4):1194-1207. doi: 10.1111/nph.12879 [38] DING Z J, XU C, YAN J Y, et al. The LRR receptor-like kinase ALR1 is a plant aluminum ion sensor[J]. Cell Res, 2024, 34(4):281-294. doi: 10.1038/s41422-023-00915-y [39] ZHU S, PAN L, VU L D, et al. Phosphoproteome analyses pinpoint the F-box protein SLOW MOTION as a regulator of warm temperature-mediated hypocotyl growth in Arabidopsis[J]. New Phytol, 2024, 241(2):687-702. doi: 10.1111/nph.19383 [40] PAN T, GAO S, CUI X, et al. APC/CCDC20 targets SCFFBL17 to activate replication stress responses in Arabidopsis[J]. Plant Cell, 2023, 35(2):910-923. doi: 10.1093/plcell/koac360 [41] van den BURG H A, TSITSIGIANNIS D I, ROWLAND O, et al. The F-box protein ACRE189/ACIF1 regulates cell death and defense responses activated during pathogen recognition in tobacco and tomato[J]. Plant Cell, 2008, 20(3):697-719. doi: 10.1105/tpc.107.056978 -
点击查看大图
计量
- 文章访问数: 2126
- HTML全文浏览量: 757
- PDF下载量: 16
- 被引次数: 0
