Rhizosphere colonization of Bacillus subtilis biocontrol strain SEM-9 and the effect on microbial diversity in rhizosphere soil
-
摘要:目的
研究枯草芽孢杆菌生防菌株SEM-9在作物根际的定植规律及对根际土壤微生物多样性的影响。
方法通过自然转化法对SEM-9菌株进行绿色荧光蛋白标记,利用倒置荧光显微镜观察该菌株在作物根际土壤、根表面及根系内组织的定植情况;并以土传病害土壤为试验对象,通过高通量测序初步分析SEM-9菌株处理后根际土壤中微生物多样性的变化。
结果构建了稳定表达绿色荧光蛋白的重组菌株SEM-9-pGFP22。荧光显微镜观察结果表明,SEM-9-pGFP22重组菌株可以在根际土壤和根表面定植,但是无法在根内组织或细胞内定植。SEM-9菌液处理显著降低了黄瓜土传病害的发病率,提高了根际土壤真菌微生物的多样性。
结论成功建立了SEM-9菌株的GFP标记方法,明确了该菌株的根际定植规律及对黄瓜土传病害的防治效果,为后期开发替代型微生物肥料奠定了基础。
Abstract:ObjectiveIn order to study the colonization rule of Bacillus subtilis biocontrol strain SEM-9 in the rhizosphere of crops and its influence on the microbial diversity of rhizosphere soil.
MethodThe strain SEM-9 was labeled with green fluorescent protein by natural transformation method, and the colonization in rhizosphere soil, on root surface and in root tissue were observed by inverted fluorescence microscope. The changes of microbial diversity in rhizosphere soil treated with the strain SEM-9 were analyzed by high-throughput sequencing using soil with soilborne diseases as test material.
ResultThe recombinant strain SEM-9-pGFP22 stably expressing green fluorescent protein was constructed. The observations of fluorescence microscopy showed that SEM-9-pGFP22 could colonize on the rhizosphere soil and root surface, but not in root tissue or cell. After treated with the SEM-9 suspension, the incidence rate of cucumber soilborne disease significantly reduced, and the fungal diversity in rhizosphere soil increased.
ConclusionThe GFP labeling method of SEM-9 strain was successfully established, and the rhizosphere colonization rule of the strain and the control effect on cucumber soilborne diseases were clarified, which lays a foundation for the later development of alternative microbial fertilizer.
-
藏猪是少有的高原型地方猪种,是我国宝贵的地方品种资源[1]。据调查,藏猪产仔数并不低,母猪的营养水平低下和乳腺发育不佳以及生存条件较恶劣可能是导致仔猪死亡率高的主要原因[2]。乳腺的良好发育是正常泌乳的前提,仔猪的存活率与母猪的乳腺发育密不可分。因此,研究藏猪妊娠期乳腺发育状况对于判断乳腺是否正常发育以及提高藏猪繁殖能力具有重要意义。
妊娠期是母猪乳腺发育的关键时期,特别是妊娠后1/3阶段,即妊娠75 d后,乳腺快速发育,其质量快速增加,乳腺结构由怀孕初期的以脂肪细胞为主转化为怀孕后期以导管和腺泡结构为主[3];妊娠期乳腺发育受到雌二醇(Estradiol,E2)、孕酮(Progesterone,P)、催乳素(Prolactin,PRL)[4]等激素的调控。E2对于乳腺导管的伸长和分支具有重要作用[5];P与E2相似,由卵巢分泌,调控乳腺组织的导管分枝、腺泡形成[6];PRL促进乳腺腺泡的发育及乳汁的分泌[7-8]。此外,在信号通路方面,PI3K/Akt是细胞内重要的信号转导通路,在乳腺细胞的增殖、分化、凋亡等活动中发挥重要的生物学功能[9-11]。Jak2/STAT5信号通路对乳腺腺泡的生成和多种乳汁蛋白基因的转录有着重要的调控作用[12-13]。但是,目前对于藏猪妊娠期乳腺的发育情况及激素和信号通路调控尚不清楚。
本试验以藏猪为对象,选取妊娠期不同时间点,在研究乳腺发育形态的基础上,进一步探索不同时间点血清中E2、P、PRL的水平,乳腺中激素受体的表达及乳腺发育关键信号通路PI3K/Akt和Jak2/STAT5的变化。研究旨在初步探究妊娠期藏猪乳腺发育过程及其潜在调控机制,为日后藏猪乳腺发育规律的揭示和地方品种的保护提供科学依据。
1. 材料与方法
1.1 试验动物
妊娠藏猪选取4个时间点(妊娠33、50、75和90 d)进行屠宰采样,采集血后离心取血清,采集第3、4对乳腺组织提取蛋白质,取第4对靠近乳头部乳腺进行石蜡切片染色。
1.2 试验材料
雌二醇、孕酮和催乳素ELISA试剂盒购于南京建成生物工程有限公司;催乳素受体(Prolactin receptor,PRLR)(货号:382057)、雌激素受体(Estrogen receptor,ER)(货号:220467)、孕酮受体(Progesterone receptor, PR)(货号:220124)抗体购于正能生物有限公司,蛋白酪氨酸激酶2(Janus kinase 2, Jak2)(货号:3230)、磷酸化蛋白酪氨酸激酶2 (p-Jak2)(货号:3771)、信号转导及转录激活因子5 (Signal transducers and activators of transduction 5, STAT5)(货号:9359)、磷酸化信号转导及转录激活因子5 (p-STAT5)(货号:4322)、磷脂酰肌醇三激酶(Phosphatidylinositol 3-kinase, PI3K)(货号:4249)、磷酸化磷脂酰肌醇三激酶(p-PI3K)(货号:4228)、蛋白激酶B (Protein kinase B, AKT)(货号:9272)和磷酸化蛋白激酶B (p-AKT)(货号:4060)抗体购于Cell Signaling Technology公司;苏木素染液、伊红染液、苏木素分化液和苏木素返蓝液购于塞维尔生物公司;BCA蛋白定量试剂盒购自白泰克生物技术有限公司(北京);ECL化学发光液购自上海雅酶生物医药科技有限公司。
1.3 测定指标及方法
乳腺采集和HE染色:切下右侧乳腺腹腺体(第4对乳腺),体积分数为4%的多聚甲醛溶液固定24 h,石蜡包埋,切片,进行HE染色,显微镜下观察并拍照。
蛋白质免疫印迹(Western blot):按照每10 μg乳腺组织加入100 μL裂解液进行匀浆,试剂盒抽提法提取蛋白质,按照BCA蛋白定量试剂盒进行蛋白浓度测定,调整蛋白浓度并用5×loading buffer制样,按照每孔20 μg总蛋白上样电泳,经电转至聚偏二氟乙烯(PVDF) 膜上后封闭,孵育一抗过夜,并用TBST缓冲液洗净孵育二抗,TBST缓冲液洗净后按照ECL化学发光底物说明书1∶1配置工作液,使PVDF膜与其充分反应30 s,置于曝光仪中曝光显色并拍照。
血清激素检测:按照南京建成ELISA试剂盒说明书进行检测;加入准备好的样品、标准品和生物素抗原,37 ℃条件下反应30 min;洗板5次,加入亲和素−HRP,37 ℃条件下反应30 min;洗板5次,加入显色液A、B,37 ℃条件下显色10 min;加入终止液;10 min之内读取D450 nm,计算浓度。
1.4 统计分析
数据结果用平均值±标准误表示,统计分析采用SigmaPlot 12.5软件分析,采用单因素方差分析,并用Duncan’s法对各组进行多重比较分析。
2. 结果与分析
2.1 藏猪妊娠期乳腺形态变化
乳腺HE染色结果如图1所示,妊娠33 d时,藏猪乳腺中主要是导管结构;50 d时,乳腺中出现少量腺泡结构;75 d时,乳腺中腺泡结构快速增多;至90 d时,乳腺中主要是腺泡结构。
2.2 藏猪妊娠期乳腺发育标志蛋白的表达
利用Western blot方法检测了藏猪妊娠期不同时间点乳腺发育标志蛋白Elf-5和PLIN2的表达,结果发现,在妊娠50、75和90 d时,Elf-5蛋白表达水平显著高于33 d;PLIN2的蛋白水平在75和90 d时显著高于33 d (图2、图3)。
图 3 藏猪妊娠期不同时间点乳腺发育标志蛋白相对表达量相同标志蛋白柱子上方的不同小写字母表示差异显著(P<0.05,Duncan’s法)Figure 3. Relative expressions of marker proteins for mammary gland development at different time points during gestation in Tibetan pigsDifferent lowercase letters on bars of the same marker protein indicate significant differences (P<0.05,Duncan’s method)2.3 藏猪妊娠期血清激素水平
如表1所示,妊娠期不同时间点藏猪血清中E2、P和PRL的水平随着妊娠的进行呈升高趋势。其中E2水平逐渐升高,到妊娠90 d达到最高水平,为42.82 ng/L;P水平在75和90 d显著高于33和50 d,90 d P水平达到36.76 μg/L;PRL水平在50 d升高,75 d显著高于33 d但与50 d无显著差异,90 d时达到最高,为66.53 μg/L,显著高于其他时间点。
表 1 藏猪妊娠期不同时间点的血清激素水平1)Table 1. Serum hormone levels at different time points during gestation in Tibetan pigst妊娠/d
Days of gestationρ(E2)/
(ng·L−1)ρ(P)/
(μg·L−1)ρ(PRL)/
(μg·L−1)33 14.56±0.82a 31.55±1.15a 44.86±1.36a 50 22.30±0.71b 32.36±0.62a 54.12±2.73b 75 27.56±0.91c 35.19±0.90b 52.91±1.31b 90 42.82±2.25d 36.76±0.94b 66.53±2.87c 1)同列数据后的不同小写字母表示差异显著(P<0.05,Duncan’s 法)
1) Different lowercase letters in the same column indicate significant differences (P<0.05,Duncan’s method)2.4 藏猪妊娠期激素受体的蛋白表达模式
利用Western blot方法检测了藏猪妊娠期不同时间点乳腺中激素受体的蛋白表达。结果显示,妊娠50 d时PRLR表达水平显著高于33 d,在90 d时达到最高;ER表达水平在50 d时显著增加,在50、75和90 d时水平相当;在75 d时PR表达水平显著高于33和50 d,在90 d时达到更高(图4、图5)。
图 5 藏猪妊娠期不同时间点乳腺中激素受体蛋白相对表达量相同激素受体柱子上方的不同小写字母表示差异显著(P<0.05,Duncan’s法)Figure 5. Relative expressions of the hormone receptors at different time points during gestation in mammary glands of Tibetan pigsDifferent lowercase letters on bars of the same hormone receptor indicate significant differences (P<0.05,Duncan’s method)2.5 藏猪妊娠期乳腺发育相关信号通路变化
利用Western blot方法检测藏猪妊娠期不同时间点乳腺发育相关信号通路Jak2/STAT5和PI3K/AKT的激活情况。由图6、图7可知,妊娠75 d时,Jak2、STAT5、PI3k和AKT的磷酸化水平显著升高,90 d时,Jak2、STAT5和PI3k的磷酸化水平显著升高,提示Jak2/STAT5和PI3K/AKT信号通路被显著激活。
图 7 藏猪妊娠期不同时间点乳腺中Jak2/STAT5和PI3K/AKT信号通路的蛋白相对表达量相同信号通路柱子上方的不同小写字母表示差异显著(P<0.05,Duncan’s法)Figure 7. Relative expressions of proteins from Jak2/STAT5 and PI3K/AKT signaling pathways at different time points during gestation in mammary glands of Tibetan pigsDifferent lowercase letters on bars of the same signaling pathway indicate significant differences (P<0.05,Duncan’s method)3. 讨论与结论
3.1 藏猪妊娠期乳腺形态变化及发育标志蛋白的表达模式
妊娠期是母猪乳腺发育的重要时期。我们的研究发现,妊娠33 d时,藏猪乳腺中主要是导管结构,50 d时出现少量腺泡结构,75 d时腺泡快速增多,90 d时乳腺中主要是腺泡结构。与我们的研究结果类似的是,Ji等[14]研究发现,母猪乳腺在妊娠45 d仅316 g,75 d乳腺质量达1606 g,90 d达到2357 g;此外,Kensinger等[15]研究表明,母猪乳腺在90 d时,腺泡数量达到最大,在90~105 d,乳腺腺泡开始分泌并蓄积大量乳汁,泌乳活动即将开始等。高慧杰等[16]在奶山羊的妊娠前期也发现乳腺并没有进入快速增殖分化阶段,而是代谢和呼吸作用增强,妊娠中期有大量细胞增殖分化。
Elf-5在妊娠期和哺乳期对乳腺腺泡的增殖和分化有重要作用,是调节乳腺发育中必不可少的调控因子[17-18]。PLIN2是调控乳脂生成的关键分子[19]。我们的研究结果显示,Elf-5和PLIN2在妊娠50 d后表达量显著升高,结合乳腺的形态和Elf-5、PLIN2蛋白水平,说明50 d乳腺开始发育出腺泡,75 d乳腺进入快速发育的阶段,90 d达到更高的发育程度,其中主要是腺泡结构。
3.2 藏猪妊娠期血清激素水平、乳腺激素受体及关键信号通路的蛋白表达模式
妊娠期乳腺发育受到多种激素的调节,其中,E2、P是调节妊娠期乳腺发育的主要激素,PRL是调节泌乳期乳腺发育及泌乳的的主要激素[5-6]。激素通过与其受体结合发挥作用,若敲除其受体,则乳腺无法正常发育[20]。方莉莉等[21]研究发现,在牦牛的妊娠早中期,PR在乳腺组织中表达量较少,与本文藏猪妊娠早中期PR蛋白表达较少一致。本文研究结果显示,E2、P和PRL等激素及其受体在妊娠期呈升高趋势,与乳腺的发育程度相吻合。类似地,Horigan等[22]通过体外给卵巢切除并抑制PRL分泌的猪注射E2、P、E2+PRL、E2+PRL+P,结果表明,注射E2+PRL+P这3种激素的组合方式才能最大程度地促进乳腺导管和腺泡产生,说明E2、P和PRL之间的相互作用对母猪乳腺发育起着关键作用。
PI3K/Akt和Jak2/STAT5信号通路在乳腺发育中起着重要调控作用。其中,PI3K/Akt是细胞内参与细胞信号转导的重要通路,参与细胞生长、增殖及分化等细胞过程[9, 23]。Meng等[10-11]研究表明,PI3K/Akt信号通路对乳腺发育和乳腺细胞的增殖具有重要作用。此外,JAK2/STAT5信号通路对乳腺腺泡的生成和多种乳汁蛋白基因的转录有着重要的调控作用[12-13]。我们的研究结果表明,PI3K/Akt和Jak2/STAT5信号通路在妊娠75 d后被显著激活,这与乳腺的高度发育及泌乳活动的开始有关。Palin等[24]的研究表明,梅山猪与大白猪乳腺组织实质中STAT5A和STAT5B的表达水平存在差异,STAT5发生磷酸化后易位至细胞核,与产乳靶基因启动子结合,激活并维持泌乳,梅山猪妊娠期乳腺组织中STAT5A和STAT5B的表达水平更高,能够生成更多的磷酸二聚体易位至核,同时梅山猪乳腺发育情况更好,具有更高的泌乳力。
综上所述,本文研究了藏猪妊娠期乳腺形态和乳腺发育标志蛋白、相关激素及信号通路的变化。结果发现,在藏猪妊娠过程中,其乳腺在妊娠50 d开始腺泡发育,75 d乳腺进入腺泡快速发育期,90 d发育程度更高,同时伴随着血清中乳腺发育相关激素(E2、P和PRL)和乳腺中激素受体表达的显著升高,以及乳腺发育相关通路PI3K/AKT和Jak2/STAT5的激活。研究结果为认识藏猪的乳腺发育和繁殖功能奠定了科学依据,为保护藏猪资源奠定了理论基础。
-
图 2 SEM-9菌株质粒提取结果
M为DNA marker DL2000,1为菌株中提取的pGFP22质粒,2为阳性转化子提取的pGFP22质粒,3、4为阳性转化子繁殖3(约200代)和4 d(约280代)提取的质粒
Figure 2. Plasmid extraction results of SEM-9 strain
M is DNA marker DL2000, Lane 1 is the pGFP22 plasmid extracted from the preserved strain, lane 2 is the pGFP22 plasmid extracted from the positive transformant with green fluorescence, and lanes 3 and 4 are the plasmids extracted from the positive transformant breeding for three (about 200 generations) and four days (about 280 generations)
图 4 SEM-9-pGFP22菌株的根际定植分析
图中绿色荧光为SEM-9-pGFP22菌落,其中A、B分别为SEM-9-pGFP22菌株在根际土壤和根表面的定植;C、D分别为空白对照和SEM-9-pGFP22菌株在根内的定植
Figure 4. Analysis of rhizosphere colonization of the recombinant strain SEM-9-pGFP22
The green fluorescence in the figure shows the colony of the recombinant strain SEM-9-pGFP22. A and B are the recombinant strains SEM-9-pGFP22 colonizing in rhizosphere soil and on the root surface, respectively; C and D are blank control and the recombinant strain SEM-9-pGFP22 colonizing in roots, respectively
表 1 SEM-9处理后根际土壤微生物多样性分析1)
Table 1 Analysis of microbial diversity in rhizosphere soil after treated with the strain SEM-9
病原菌 Pathogen 处理 Treatment Chao Ace Shannon Simpson 细菌
BacteriaCK 1973±207.4a 1940±218.3a 5.84±0.389a 0.0091±0.006a SEM-9 2112±37.2a 2114±68.2a 6.08±0.196a 0.0068±0.003a SJ 2064±121.4a 2035±121.8a 6.06±0.266a 0.0059±0.003a 真菌
FungiCK 176±15.9a 182±2.8a 2.57±0.321a 0.24±0.08a SEM-9 206±7.2b 206±6.2b 2.93±0.158b 0.11±0.02b SJ 196±23.1ab 197±25.3ab 2.65±0.363a 0.19±0.05a 1) 同种病原菌同列数据后的不同小写字母表示不同处理差异显著(P<0.05, LSD 法)
1) Different lowercase letters after the data of the same pathogen of the same column indicate significant differences among different groups (P<0.05, LSD method)表 2 SEM-9处理后根际土壤中镰刀菌在种属水平上的相对丰度变化1)
Table 2 Analysis of relative abundance of Fusarium in rhizosphere soil at species and genus level after treated by the suspension of the strain SEM-9
属 Species CK SEM-9 SJ 镰刀菌属
Fusarium0.73±0.131a 0.52±0.049b 0.51±0.104b 尖孢镰刀菌
F.oxysporum0.14±0.025a −0.25±0.282b 0.26±0.182a 1) 表中数据为相对丰度的对数值,同行数据后的不同小写字母表示不同处理差异显著(P<0.05, LSD 法)
1) The data in the table are the logarithm values of relative abundances, and different lowercase letters after the data of the same row indicate significant differences among different groups (P<0.05, LSD method) -
[1] LUGTENBERG B J J, DEKKERS L, BLOEMBERG G V. Molecular determinants of rhizosphere colonization by Pseudomonas[J]. Annual Review of Phytopathology, 2001, 39(0): 461-490. doi: 10.1146/annurev.phyto.39.1.461
[2] COMPANT S, CLÉMENT C, SESSITSCH A. Plant growth-promoting bacteria in the rhizo and endosphere of plants: Their role, colonization, mechanisms involved, and prospects for utilization[J]. Soil Biology and Biochemistry, 2010, 42(5): 669-678. doi: 10.1016/j.soilbio.2009.11.024
[3] GHIRARDI S, DESSAINT F, MAZURIER S, et al. Identification of traits shared by rhizosphere competent strains of fluorescent pseudomonads[J]. Microbial Ecology, 2012, 64(3): 725-737. doi: 10.1007/s00248-012-0065-3
[4] HALL-STOODLEY L, COSTERTON J W, STOODLEY P. Bacterial biofilms: From the natural environment to infectious diseases[J]. Nature Reviews Microbiology, 2004, 2(2): 95-108. doi: 10.1038/nrmicro821
[5] BAIS H P, FALL R, VIVANCO J M. Biocontrol of Bacillus subtilis against infection of Arabidopsis roots by Pseudomonas syringae is facilitated by biofilm formation and surfactin production[J]. Plant Physiology, 2004, 134(1): 307-319. doi: 10.1104/pp.103.028712
[6] CHEN Y, CAO S G, CHAI Y R, et al. A Bacillus subtilis sensor kinase involved in triggering biofilm formation on the roots of tomato plants[J]. Molecular Microbiology, 2012, 85(3): 418-430. doi: 10.1111/j.1365-2958.2012.08109.x
[7] CHU F, KEARNS D B, BRANDA S S, et al. Targets of the master regulator of biofilm formation in Bacillus subtilis[J]. Molecular Microbiology, 2006, 59(4): 1216-1228. doi: 10.1111/j.1365-2958.2005.05019.x
[8] KEARNS D B, CHU F, BRANDA S S, et al. A master regulator for biofilm formation by Bacillus subtilis[J]. Molecular Microbiology, 2005, 55(3): 739-749.
[9] TASAKI S, NAKAYAMA M, SHOJI W. Morphologies of Bacillus subtilis communities responding to environmental variation[J]. Development, Growth & Differentiation, 2017, 59(5): 369-378. doi: 10.1111/dgd.12383
[10] HAMON M A, LAZAZZERA B A. The sporulation transcription factor Spo0A is required for biofilm development in Bacillus subtilis[J]. Molecular Microbiology, 2001, 42(5): 1199-1209.
[11] VLAMAKIS H, CHAI Y, BEAUREGARD P, et al. Sticking together: Building a biofilm the Bacillus subtilis way[J]. Nature Reviews Microbiology, 2013, 11(3): 157-168. doi: 10.1038/nrmicro2960
[12] LI Q R, LIAO S T, ZHI H Y, et al. Characterization and sequence analysis of potential biofertilizer and biocontrol agent Bacillus subtilis strain SEM-9 from silkworm excrement[J]. Canadian Journal of Microbiology, 2019, 65(1): 45-58. doi: 10.1139/cjm-2018-0350
[13] 梁肇均, 林毓娥, 何晓明, 等. 黄瓜土传性病害的发生与防治技术[J]. 长江蔬菜, 2019(7): 52-54. [14] COOMBS J T, FRANCO C M M. Visualization of an endophytic Streptomyces species in wheat seed[J]. Applied and Environmental Microbiology, 2003, 69(7): 4260-4262. doi: 10.1128/AEM.69.7.4260-4262.2003
[15] WEYENS N, BOULET J, ADRIAENSEN D, et al. Contrasting colonization and plant growth promoting capacity between wild type and a GFP-derative of the endophyte Pseudomona sputida W619 in hybrid poplar[J]. Plant and Soil, 2012, 356(1/2): 217-230. doi: 10.1007/s11104-011-0831-x
[16] WANG X J, LI M J, YAN Q, et al. Across genus plasmid transformation between Bacillus subtilis and Escherichia coli and the effect of Escherichia coli on the transforming ability of free plasmid DNA[J]. Current Microbiology, 2007, 54(6): 450-456. doi: 10.1007/s00284-006-0617-1
[17] 李瑞芳, 薛雯雯, 黄亮, 等. 枯草芽孢杆菌感受态细胞的制备及质粒转化方法研究[J]. 生物技术通报, 2011(5): 227-230. [18] ZBORALSKI A, FILION M. Genetic factors involved in rhizosphere colonization by phytobeneficial Pseudomonas spp.[J]. Computational and Structural Biotechnology Journal, 2020, 18: 3539-3554. doi: 10.1016/j.csbj.2020.11.025
[19] 崔晓辰. 根际微生物与土壤植物关系的研究进展[J]. 现代农业研究, 2021, 27(5): 34-35. doi: 10.3969/j.issn.1674-0653.2021.05.015 [20] 刘京伟, 李香真, 姚敏杰. 植物根际微生物群落构建的研究进展[J]. 微生物学报, 2021, 61(2): 231-248. [21] 刘王锁, 李海泉, 何毅, 等. 根际微生物对植物与土壤交互调控的研究进展[J]. 中国土壤与肥料, 2021(5): 318-327. doi: 10.11838/sfsc.1673-6257.20292 [22] 郑婷婷. 生防菌GHt_q6对黄瓜根系土壤微生态及根结线虫的影响[D]. 晋中: 山西农业大学, 2019. -
期刊类型引用(2)
1. 韩江涛,张帅博,秦雅蕊,韩硕洋,张雅康,王吉庆,杜清洁,肖怀娟,李猛. 甜瓜β-淀粉酶基因家族的鉴定及对非生物胁迫的响应. 生物技术通报. 2025(03): 171-180 . 百度学术
2. 梅玉琴,刘意,王崇,雷剑,朱国鹏,杨新笋. 甘薯PHB基因家族的全基因组鉴定和表达分析. 作物学报. 2023(06): 1715-1725 . 百度学术
其他类型引用(3)