Geographical distribution and floristic significance of Pyrola in China
-
摘要:
本研究通过查阅植物志和标本，获取了中国鹿蹄草属（Pyrola）的地理分布信息。结果显示：（1）我国所有省份皆有鹿蹄草属分布，其中吉林省和四川省分布最多，达到了10种；（2）鹿蹄草属为北温带分布型，但种一级以中国特有分布型最多；（3）从海平面以上至海拔4 100 m皆有鹿蹄草属分布，主要集中分布于海拔1 500 ~ 2 500 m；（4）我国鹿蹄草属分布区与针叶林高度重叠，最多见于寒温性针叶林，具有植被恢复的指示意义。
Abstract:This study presents comprehensive data on the geographical distribution of the genus Pyrola across China, gathered through extensive literature review and specimen analysis. Results showed that: (1) Pyrola is ubiquitously distributed throughout all Chinese provinces, with Jilin and Sichuan recording the highest species counts, each hosting 10 species. (2) The Pyrola genus exhibits a predominant north temperate distribution type, but at the species level, the Chinese endemic distribution type is most common. (3) The altitudinal range of Pyrola spans from sea level to 4 100 m, primarily between 1 500–2 500 m. (4) There is a significant overlap between Pyrola distribution areas and coniferous forests in China, predominantly within boreal coniferous forests, indicating their significance for vegetation restoration.
-
Keywords:
- Pyrola /
- Geographical distribution /
- Flora /
- Boreal coniferous forest
-
VQ（Valine-glutamine）是植物特有的一类蛋白，因其包含一个特有的保守基序FxxxVQxLTG而得名（F：苯丙氨酸，V：缬氨酸，Q：谷氨酰胺，L：亮氨酸，T：苏氨酸，G：甘氨酸，x 代表任意氨基酸）[1-3]。首个VQ蛋白SIB1（SIGMA FACTOR-BINDING PROTEIN1）在拟南芥（Arabidopsis thaliana (L.) Heynh.）中被发现[4]，随后，从拟南芥中又陆续鉴定出33个VQ蛋白。目前，科研人员已在茶（Camellia sinensis (L.) O. Ktze. ）[2]、大豆（Glycine max (L.) Merr.）[5]、小麦（Triticum aestivum L.）[6]及烟草（Nicotiana tabacum L.）[3]中分别鉴定出25、75、113和59个VQ基因。
VQ蛋白在植物应答生物与非生物胁迫中发挥着重要作用[7, 8]，其作用机制主要有两种。一是通过与WRKY转录因子相互作用发挥功能，如AtVQ23和AtVQ16通过与WRKY33的WRKY结构域结合，增强后者与靶基因的结合能力，从而影响植物的抗病能力[9]。而AtVQ9与WRKY8互作抑制了WRKY8与W-box的结合活性，从而负调控植物的抗盐性[10]。在番茄（Solanum lycopersicum L.）叶片中，SlVQ7与SlWRKY37互作可提高后者的稳定性，并转录激活其下游靶基因，从而正调控茉莉酸及黑暗诱导的叶片衰老[11]。二是作为丝裂原活化蛋白激酶（MAPK）家族的磷酸化底物发挥作用[12-14]。如拟南芥MKS1（AtVQ21）作为AtMPK4的底物，与WRKY25和WRKY33相互作用，有助于激活MPK4调节的病原体防御反应[3]。
多年生黑麦草（Lolium perenne L.）为禾本科早熟禾亚科黑麦草属植物，是世界上广泛种植的冷季型草坪草及牧草[15]，具有分蘖多、成坪速度快、粗蛋白含量高、耐牧性好等优点。多年生黑麦草在生长发育过程中常面临高温、盐碱、荫蔽等非生物胁迫，严重影响草坪质量和牧草品质。因此，鉴定多年生黑麦草VQ基因，探究其在应答非生物胁迫中的作用十分必要。然而，目前尚未见相关报道。本研究系统鉴定了多年生黑麦草VQ蛋白家族成员，并对其保守结构域、进化关系、共线性、亚细胞定位及其在高温、盐碱、黑暗等非生物胁迫下的表达模式进行了分析，研究结果旨在为进一步探讨VQ蛋白在多年生黑麦草生长发育及应答非生物逆境胁迫中的作用奠定基础。
1. 材料与方法
1.1 实验材料
本研究以多年生黑麦草品种‘百灵鸟’为实验材料。黑麦草种子于营养土（基质∶蛭石=1∶1）中萌发，萌发后转移至光照培养箱，光周期为16 h / 8 h（光 / 黑暗），温度为25 ℃ / 22 ℃（光 / 黑暗）。待植株生长1个月后进行黑暗、高温和NaCl处理。黑暗处理的条件为：25 ℃ / 22 ℃，24 h黑暗，分别于处理0、2、8 d后取样。高温处理的条件为：42 ℃，光周期16 h / 8 h，分别于处理0、1、24 h后取样。NaCl处理的条件为：25 ℃ / 22 ℃，光周期16 h / 8 h，使用200 mL 200 mmol/L的NaCl溶液进行浇灌，分别于盐处理0、3、38 h后取样。
1.2 多年生黑麦草VQ基因家族成员的鉴定
以Pfam数据库中VQ结构域（PF05678）的隐马尔可夫模型（HMM）图谱作为查询，通过HMMER搜索多年生黑麦草基因组，鉴定了可能的VQ蛋白序列。使用Pfam数据库（http://pfam.xfam.org/）进一步确认了这些序列。在拟南芥基因组数据库（TAIR）中下载拟南芥全基因组序列及其注释信息。
1.3 多年生黑麦草VQ家族成员的理化性质及亚细胞定位
使用在线网站（https://web.expasy.org/compute_pi/）分析多年生黑麦草VQ蛋白的理化性质。使用BIOINF（http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-）在线工具对多年生黑麦草VQ家族成员进行亚细胞定位预测。
1.4 多年生黑麦草VQ家族基因系统发育树构建
使用MEGAX（http://www.megasoftware.net）软件的ML法构建系统发育树，进行系统发育分析，并用iTOL（https://itol.embl.de/）网站进行美化。
1.5 多年生黑麦草VQ保守基序及基因结构分析
利用MEME（https://meme-suite.org/meme/tools/meme）网站对VQ蛋白的氨基酸序列进行保守基序分析。结合植物基因组序列和注释文件，用TBtools提取VQ基因的基因组序列及CDS序列，用TBtools软件工具中的Gene Structure View进行可视化。
1.6 VQ基因染色体分布及种内共线性分析
结合植物基因组序列和注释文件，使用TBtools软件工具中的Amazing Gene Location From GFF3/GTF File进行可视化，使用TBtools软件工具中的one step MCScanX进行种内基因组比对，并用Multiple Synreny Plot进行可视化。
1.7 胁迫条件下多年生黑麦草VQ基因的表达模式分析
利用Trizol试剂（Takara，大连）提取多年生黑麦草在不同非生物胁迫处理后的总RNA。使用反转录试剂盒Hifair® Ⅲ 1st Strand cDNA Synthesis SuperMix for qPCR（gDNA digester plus）（Yeasen，上海）合成cDNA。利用Hieff® qPCR SYBR Green Master Mix（High Rox Plus）试剂盒（Yeasen，上海）在Step One Plus Real-Time PCR Systems系统上进行qRT-PCR反应。反应体系总体积为20 μL，包括2 μL cDNA，上、下游引物各0.8 μL，10 μL SYBR，用ddH2O补足体积。反应程序为：95 ℃预变性5 min；95 ℃变性10 s，60 ℃退火20 s，72 ℃延伸20 s，40个循环，系统自带熔解曲线。引物序列见附表1
1 ，由擎科公司（擎科，武汉）合成，多年生黑麦草eEF1A(s)为内参基因。采用2−△△Ct方法计算候选基因的相对表达量，实验设3次重复。利用TBtools绘制热图，进行表达模式分析。2. 结果与分析
2.1 LpVQ家族成员鉴定
为了系统分析多年生黑麦草的VQ基因，以VQ结构域（PF05678）为参考，从全基因组水平挖掘其VQ家族基因。结果显示，多年生黑麦草中共存在52个VQ基因。对这些基因的编码产物LpVQ家族成员进行理化性质分析，发现LpVQs的氨基酸长度介于118~427 aa，蛋白分子量范围为11.67~43.84 kD，等电点在5.01~11.6。通过BIOINF在线网站对LpVQ家族成员进行亚细胞定位预测，发现48个蛋白定位在细胞核中，LpVQ22为叶绿体与细胞核共定位，LpVQ33为细胞壁与细胞核共定位，LpVQ37为细胞膜与细胞核共定位，而LpVQ45预测定位在细胞膜中（附表2
1 ）。2.2 LpVQ家族系统进化树分析
为了分析LpVQs的进化关系，利用多年生黑麦草的52个LpVQ蛋白和拟南芥的34个AtVQ蛋白共同构建系统进化树。结果显示，拟南芥与多年生黑麦草的VQ家族成员可划分为7个亚家族（Ⅰ~Ⅶ），分别包含5、8、6、11、12、6和4个LpVQ蛋白。亚家族Ⅰ中只包含LpVQs，亚家族Ⅳ、Ⅴ和Ⅵ中含有较多的LpVQ家族成员，而AtVQ家族成员则更多地归类于亚家族Ⅱ、Ⅲ和Ⅶ（图1）。进一步分析发现，多年生黑麦草与拟南芥的VQ家族间存在4对直系同源基因，分别是LpVQ6 / AtVQ7、LpVQ11 / AtVQ9、LpVQ50 / AtVQ32和LpVQ5 / AtVQ31。
![]()
图 1 拟南芥和多年生黑麦草VQ家族成员系统进化树分析Figure 1. Phylogenetic analysis of VQ family members in Arabidopsis thaliana and Lolium perenne2.3 LpVQ家族成员保守基序和基因结构分析
为了研究蛋白的结构多样性，利用MEME在线软件对52个LpVQ蛋白的保守基序进行分析，结果发现17%的LpVQs含有1个Motif，48%的LpVQs含有两个Motif，19%的LpVQs含有3个Motif，6%的LpVQs含有4个Motif，10%的LpVQs含有5个Motif。所有家族成员均具有Motif1，表明它是LpVQ家族的核心保守结构域（图2）。
此外，本研究用TBtools工具对LpVQs进行了基因结构分析。结果发现，46个LpVQ基因没有内含子，5个LpVQ基因含有1个内含子（LpVQ24、LpVQ29、LpVQ31、LpVQ32、LpVQ52），仅有1个基因（LpVQ161）含有两个内含子（附图1
1 ）。2.4 LpVQ家族成员染色体定位以及共线性分析
染色体定位分析结果表明，LpVQ基因分布在多年生黑麦草的7条染色体上，其中分布在chr1上的有7个，占13%，分布在chr2上的最多，为11个，占21%，分布在chr3上的为5个，占10%，分布在chr4、chr6和chr7上的基因均各有9个，各占17%，而分布在chr5上的最少，仅有2个，占4%（图3）。
为了了解LpVQ基因之间的进化关系，本研究对VQ家族成员进行了种内共线性分析。结果表明，3对VQ基因之间存在共同进化关系，分别是LpVQ16与LpVQ36、LpVQ34与LpVQ38以及LpVQ14与LpVQ19（图4）。
2.5 LpVQ基因响应非生物胁迫的表达模式分析
从7个亚家族中分别随机选取3个LpVQ基因，检测其在黑暗、高温以及盐胁迫下多年生黑麦草中的表达量变化。结果显示，黑暗处理后，只有4个LpVQ基因的表达持续上调。大部分LpVQ基因的表达量在黑暗处理2 d后下降，其中，LpVQ2、LpVQ8、LpVQ14、LpVQ22及LpVQ38的表达量在黑暗处理8 d后又显著上升（图5：A）。
![]()
图 5 LpVQ基因响应黑暗、高温及NaCl胁迫的表达模式分析A：黑暗处理；B：高温处理；C：NaCl处理。红色、蓝色和黄色分别代表上调、下调和没有表达差异的基因。Figure 5. Expression analyses of LpVQs under dark, high temperature and NaCl treatmentsA: Dark treatment; B: Heat treatment; C: NaCl treatment. Red, blue, and yellow represent up-regulated, down-regulated, and no differential genes, respectively.42 ℃高温处理后，LpVQ的表达模式被分为3类。LpVQ18、LpVQ8、LpVQ28、LpVQ2、LpVQ7、LpVQ4及LpVQ12受1 h高温处理的显著诱导，而高温处理24 h后，其表达又显著下调。LpVQ9、LpVQ22、LpVQ15、LpVQ3、LpVQ42、LpVQ6、LpVQ10、LpVQ1、LpVQ38的表达量在1 h和24 h高温处理后均显著下降。LpVQ5、LpVQ27、LpVQ34、LpVQ14和LpVQ44的表达量只在24 h高温处理后下降（图5：B）。
NaCl处理3 h后，13个LpVQ基因显著上调表达，其中LpVQ4的表达上升了13倍。而NaCl处理38 h后，大部分基因的表达明显下调。仅有3个基因的表达受盐诱导持续上调，分别是 LpVQ8、LpVQ10和LpVQ38（图5：C）。
3. 讨论
近年来，VQ基因已经在多种植物中被发掘鉴定，包括拟南芥、水稻（Oryza sativa L.）、大豆和葡萄（Vitis vinifera L.）等[16-19]，而关于多年生黑麦草VQ家族基因的鉴定和分析还鲜见报道。本研究通过全基因组鉴定，得到52个多年生黑麦草LpVQ基因，并系统分析了其理化特性、蛋白结构、基因结构、共线性、进化关系以及响应逆境胁迫的表达模式。
基因结构分析发现，88.46%的LpVQ基因无内含子，与之类似，在拟南芥[16]、水稻[17]、玉米（Zea mays L.）[20]和马铃薯（Solanum tuberosum L.）[21]中，无内含子的VQ基因分别占88.2%、92.5%、88.5%和91.89%。而苔藓（Bryophyta）[22]中无内含子的VQ基因比例则远低于以上物种，仅为28.0%，说明VQ基因可能在长期进化过程中丢失了内含子。
系统进化分析结果显示，多年生黑麦草和拟南芥的VQ家族成员可被分为7个亚家族，且LpVQs与AtVQs亲缘关系较远，只存在4对直系同源基因，分别是LpVQ6 / AtVQ7、LpVQ11 / AtVQ9、LpVQ50 / AtVQ32和LpVQ5 / AtVQ31。根据拟南芥VQ基因的功能，可以推测多年生黑麦草中其直系同源VQ基因的功能。有研究指出，AtVQ9作为AtWRKY8的抑制子，能提高拟南芥的耐盐性[9]，由此推测，其直系同源基因LpVQ11可能也在多年生黑麦草抵御盐胁迫中发挥重要作用。
VQ基因在植物响应非生物胁迫中发挥着重要作用。如，过表达PtVQ1基因显著增强了转基因拟南芥的耐盐性[23]；而过表达MdVQ37则提高了转基因苹果（Malus pumila Mill.）对高温胁迫的敏感性[24]。本研究对黑暗、高温及盐胁迫处理下多年生黑麦草的VQ表达量进行了分析。结果显示，LpVQ基因广泛响应黑暗、高温及盐胁迫，且大部分LpVQ基因均在胁迫初期即做出快速响应，这与西葫芦（Cucurbita pepo L.）[8]、黄瓜（Cucumis sativus L.）[25]等物种中的研究结果相类似。然而，LpVQ基因的表达模式在不同胁迫处理下也存在差异。如，LpVQ1在黑暗胁迫处理后持续下调表达，而在高温和盐胁迫下持续上调表达，推测该基因可能在不同胁迫中发挥着不同的作用。尽管LpVQ可能参与多年生黑麦草对黑暗、高温及盐等胁迫响应的调控，但其具体功能还需要进一步验证。
综上所述，本研究对多年生黑麦草VQ基因家族成员进行了系统鉴定与分析，为后续验证VQ家族成员在多年生黑麦草应答生物或非生物胁迫中的功能提供了基础。
1 1) 如需查阅附表内容请登录《植物科学学报》网站（http://www.plantscience.cn）查看本期文章。 -
![]()
图 1 中国鹿蹄草属种数与海拔梯度关联图
Figure 1. Correlation map between number of Pyrola species and altitudinal gradient in China
![]()
图 2 中国鹿蹄草属种数与生境关联图
Figure 2. Correlation map between number of Pyrola species and habitat in China
表 1 中国鹿蹄草属分布区类型统计
Table 1 Distribution area types of Pyrola in China
编号
No.分布型
Distribution pattern种数
Species比例
Proportion8-2 北极、高山间断分布 1 0.04 10 旧世界温带分布 2 0.08 11 温带亚洲分布 2 0.08 14-SH 中国-喜马拉雅分布 2 0.08 14-SJ 中国-日本分布 4 0.15 15 中国特有分布 15 0.58 
下载: 导出CSV
-
[1] Qin HN, Stevens PF. Pyrola[M]//Wu WZ, Raven P, Hong DY, eds. Flora of China: Vol. 14. Beijing: Science Press; Missouri: Missouri Botanical Garden Press, 2005: 8-15.
[2] Freudenstein JV. Relationships and character transformation in Pyroloideae (Ericaceae) based on its sequences,morphology,and development[J]. Syst Bot,1999,24(3):398−408. doi: 10.2307/2419697
[3] Liu ZW,Zhou J,Liu ED,Peng H. A molecular phylogeny and a new classification of Pyrola (Pyroleae,Ericaceae)[J]. Taxon,2010,59(6):1690−1700. doi: 10.1002/tax.596004
[4] Liu ZW,Wang ZH,Zhou J,Peng H. Phylogeny of Pyroleae (Ericaceae):implications for character evolution[J]. J Plant Res,2011,124(3):325−337. doi: 10.1007/s10265-010-0376-8
[5] Liu ZW,Jolles DD,Zhou J,Peng H,Milne RI. Multiple origins of circumboreal taxa in Pyrola (Ericaceae),a group with a tertiary relict distribution[J]. Ann Bot,2014,114(8):1701−1709. doi: 10.1093/aob/mcu198
[6] 李东,杨松,宋宝安,姜林锟,胡德禹,薛伟. 鹿蹄草属植物的研究进展[J]. 贵州大学学报(自然科学版),2008,25(2):188−194. Li D,Yang S,Song BA,Jiang LK,Hu DY,Xue W. Research advances on Pyrola L.[J]. Journal of Guizhou University (Natural Sciences)
,2008,25(2):188−194. [7] 马玉心,蔡体久,崔大练. 兴安鹿蹄草分布规律的研究[J]. 北方园艺,2007(5):180−181. Ma YX,Cai TJ,Cui DL. Study on the distribution pattern of Pyrola dahurica[J]. Northern Horticulture,2007(5):180−181.
[8] 余贝贝,康永祥,张世姣,苏悦,曾珵. 太白山雅美鹿蹄草群落植物种间关联性[J]. 东北林业大学学报,2012,40(11):74−77. doi: 10.3969/j.issn.1000-5382.2012.11.019 Yu BB,Kang YX,Zhang SJ,Su R,Zeng C. Interspecific association among species of Pyrola decorata community in Taibai Mountain[J]. Journal of Northeast Forestry University,2012,40(11):74−77. doi: 10.3969/j.issn.1000-5382.2012.11.019
[9] 苏悦,康永祥,张世姣,余贝贝,薛利艳. 太白山雅美鹿蹄草和鹿蹄草种群分布格局[J]. 东北林业大学学报,2013,41(1):89−93. Su Y,Kang YX,Zhang SJ,Yu BB,Xue LY. Distribution pattern of Pyrola decorata and Pyrola calliantha population in Taibai Mountain[J]. Journal of Northeast Forestry University,2013,41(1):89−93.
[10] Selosse MA,Roy M. Green plants that feed on fungi:facts and questions about mixotrophy[J]. Trends Plant Sci,2009,14(2):64−70. doi: 10.1016/j.tplants.2008.11.004
[11] Suetsugu K,Matsuoka S,Shutoh K,Okada H,Taketomi S,et al. Mycorrhizal communities of two closely related species,Pyrola subaphylla and P. japonica,with contrasting degrees of mycoheterotrophy in a sympatric habitat[J]. Mycorrhiza,2021,31(2):219−229. doi: 10.1007/s00572-020-01002-5
[12] Logacheva MD,Schelkunov MI,Shtratnikova VY,Matveeva MV,Penin AA. Comparative analysis of plastid genomes of non-photosynthetic Ericaceae and their photosynthetic relatives[J]. Sci Rep,2016,6:30042. doi: 10.1038/srep30042
[13] Yu JJ,Wang CB,Gong X. Degeneration of photosynthetic capacity in mixotrophic plants,Chimaphila japonica and Pyrola decorata (Ericaceae)[J]. Plant Divers,2017,39(2):80−88. doi: 10.1016/j.pld.2016.11.005
[14] 朱华,马友鑫,阎丽春,胡华斌. 中国种子植物属的地理成分分布格局及其与气候和地理的关系(英文)[J]. 植物分类学报,2007,45(2):134−166. Zhu H,Ma YX,Yan LC,Hu HB. The relationship between geography and climate in the generic-level patterns of Chinese seed plants[J]. Acta Phytotaxonomica Sinica,2007,45(2):134−166.
[15] 沈泽昊,杨明正,冯建孟,李新辉,彭培好,郑智. 中国高山植物区系地理格局与环境和空间因素的关系[J]. 生物多样性,2017,25(2):182−194. Shen ZH,Yang MZ,Feng JM,Li XH,Peng PH,Zheng Z. Geographic patterns of alpine flora in China in relation to environmental and spatial factors[J]. Biodiversity Science,2017,25(2):182−194.
[16] 吴征镒. 中国种子植物属的分布区类型[J]. 云南植物研究, 1991（S4）: 1-139. Wu ZY. The distribution types of seed plant genera in China[M]. Acta Botanica Yunnanica, 1991（S4）: 1-139.
[17] 吴征镒, 周浙昆, 孙航, 李德铢, 彭华. 种子植物分布区类型及其起源和分化[M]. 昆明: 云南科技出版社, 2006: 60-64. [18] 吴征镒. 中国植被[M]. 北京: 科学出版社, 1980: 159-210. [19] 刘增力,方精云,朴世龙. 中国冷杉、云杉和落叶松属植物的地理分布[J]. 地理学报,2002,57(5):577−586. Liu ZL,Fang JY,Piao SL. Geographical distribution of species in genera Abies,Picea and Larix in China[J]. Acta Geographica Sinica,2002,57(5):577−586.
[20] Vincenot L,Tedersoo L,Richard F,Horcine H,Kõljalg U,Selosse MA. Fungal associates of Pyrola rotundifolia,a mixotrophic Ericaceae,from two Estonian boreal forests[J]. Mycorrhiza,2008,19(1):15−25. doi: 10.1007/s00572-008-0199-9
[21] Malysheva EF,Malysheva VF,Voronina EY,Kovalenko AE. Diversity of fungal communities associated with mixotrophic pyroloids(Pyrola rotundifolia,P. media and Orthilia secunda) in their natural habitats[J]. Bot Pac,2018,7(2):31−39.
-
其他相关附件
-
PDF格式
董洪进附表1 点击下载(140KB)
-
计量
- 文章访问数: 0
- HTML全文浏览量: 0
- PDF下载量: 0
