Analysis of bacterial community diversity in rhizosphere soil of continuous cropping Andrographis paniculata based on high-throughput sequencing
-
摘要:目的
研究穿心莲Andrographis paniculata连作对土壤细菌多样性及种群结构的影响。
方法通过16S rRNA 基因高通量测序对连作组(连作5年穿心莲的根际土壤)和对照组(连作5年穿心莲周边相同土质的未耕作自然土壤)样品进行测序分析，比较2组土壤细菌的丰富度、多样性和群落结构差异。
结果α 多样性指数分析结果显示，连作组土壤中细菌丰富度指数Chao1指数、Observed species指数显著低于对照组(P<0.05)，连作组土壤中的多样性指数Shannon指数和Simpson指数均低于对照组，但差异不显著。主坐标分析及分子方差分析进行组间群落结构差异显著性分析结果显示，穿心莲连作显著改变了土壤细菌群落结构(P<0.05)。从2个处理组土壤样品中共检测到2769个可操作分类单元，分属于47门885属，穿心莲连作明显改变了土壤细菌在门和属水平上的结构分布。在门分类水平上，酸杆菌门Acidobacteria、变形菌门Proteobacteria、拟杆菌门Bacteroidetes和Patescibacteria等11个细菌门是2个处理组的主要细菌类群，连作组土壤中变形菌门和酸杆菌门的相对丰度较对照组分别增加了21.97%和46.33% (P<0.05)，而Patescibacteria较对照组减少了68.99% (P<0.05)；在属分类水平上，有10个优势菌属的相对丰度变化显著(P<0.05)，其中连作组土壤中有益生防细菌黄杆菌属Flavobacterium和Haliangium的相对丰度较对照组分别减少了78.40%和54.55%，而植物病原细菌伯克氏菌属Burkholderia-Caballeronia-Paraburkholderia的相对丰度较对照组增加了804.17%。
结论穿心莲5年连作降低了土壤细菌丰富度和多样性水平，使土壤微生态系统稳定性降低，同时有益菌相对丰度显著降低，而病原菌相对丰度显著增加造成细菌群落结构改变。上述变化打破了原有的土壤微生态平衡，可能是穿心莲连作障碍的重要原因之一。
Abstract:ObjectiveTo investigate the effects of continuous cropping of Andrographis paniculata on soil bacterial diversity and community structure.
MethodSequencing analysis of samples from continuous cropping group (rhizosphere soil of A. paniculata continuous cropping for five years) and control group (uncultivated natural soil with the same texture as the soil around A. paniculata continuous cropping for five years) were performed using 16S rRNA gene high-throughput sequencing technology. The differences in richness, diversity and community structure of soil bacteria from two groups were analyzed.
ResultThe analysis results of α diversity index showed that the bacterial richness indices (Chao1 index and Observed species index) in the soil of continuous cropping group were significantly lower than those of the control group (P<0.05), and the diversity indices (Shannon index and Simpson index) in the soil of continuous cropping group were lower than those of the control group, but the differences were not significant. Principal coordinates analysis and molecular variance analysis of community structure between groups showed that continuous cropping ofA. paniculata significantly changed soil bacterial community structure (P<0.05). A total of 2769 operational taxonomic units were detected from the soil samples of two groups, belonging to 47 phyla and 885 genera. Continuous cropping ofA. paniculata obviously changed the structure distribution of soil bacteria at the phylum and genus level. At the phylum classification level, 11 bacterial phyla including Acidobacteria, Proteobacteria, Bacteroidetes, Patescibacteria and etc. were the main bacterial communities in two groups. Compared with the control group, the relative abundance of Proteobacteria and Acidobacteria in the continuous cropping group increased by 21.97% and 46.33% respectively (P<0.05), while the relative abundance of Patescibacteria reduced by 68.99% (P<0.05). At the genus classification level, the relative abundance of 10 dominant bacterial genera changed significantly (P<0.05). Among them, compared with the control group, the relative abundance of the beneficial bacteriaFlavobacterium and Haliangiumin in the soil of continuous cropping group reduced by 78.40% and 54.55% respectively, while the relative abundance of the plant pathogenic bacterium Burkholderia-Caballeronia-Paraburkholderia in the soil of continuous cropping group increased by 804.17%.
ConclusionAfter five years of A. paniculata continuous cultivation, the level of soil bacterial richness and diversity decreased. At the same time, the relative abundance of beneficial bacteria significantly reduced, while the relative abundance of pathogenic bacteria significantly increased, causing a change in the bacterial community structure. The above-mentioned changes broke the original soil microecological balance, which may be one of the important reasons for A. paniculata continuous cropping obstacles.
-
赤桉Eucalyptus camaldulensis是桃金娘科Myrtaceae桉属Eucalyptus比较广泛栽培的树种，原产于澳大利亚，与其他桉树相比，赤桉具有高抗寒性，是研究桉树抗寒分子机理的理想材料[1]。
miRNA是由microRNA基因所编码、II型RNA聚合酶(Pol II)和DCL1酶(Dicer-like 1)等加工形成的长度约18~24 nt的单链非编码RNA[2]。研究证实，miRNA通过对靶基因的剪切或抑制来实现转录后水平上的负调控，在植物的生长发育、形态建成和信号转导以及逆境胁迫响应中都发挥着重要的调控作用[3-5]。
低温是一种常见的环境因子，它限制植物只能在特定的地理位置和季节生长。植物为适应低温胁迫，在生理生化和分子水平上形成了一系列的抵御机制，包括瞬时Ca2+和脱落酸(ABA)含量的提高、植物体液组成的改变、抗氧化物含量的增加，以及一些渗透调节物质的集聚[6]。最早在拟南芥Arabidopsis thaliana中发现低温胁迫下miRNA的表达丰度会发生变化[7]，之后许多研究报道了木本植物低温胁迫下差异表达的miRNA[4,8]。华南农业大学桉树育种课题组前期对赤桉耐寒关键性转录因子ICE1(Inducer of CBF expression 1)[9]、SIZ1(SAP and Miz)苏素化[10]和HOS1(High expression of osmotically responsive gene 1)泛素化[11]通路进行了一些研究，但在赤桉ICE1通路及其上下游中是否有miRNA参与调节，仍然未知。本研究拟通过高通量测序和生物信息学软件相结合的方式，挖掘、分析涉及低温胁迫应答的赤桉miRNA，以进一步弄清赤桉低温胁迫应答的完整分子调控网络，并为后续的桉树抗寒分子育种提供一些理论依据和技术支撑。
1. 材料与方法
1.1 植物材料与处理
赤桉无性系CV103无菌组培苗取自广东省森林植物种质创新与利用重点实验室。在MS培养基(pH=5.8)附加0.5 mg·L−1 6-BA、0.1 mg·L−1 NAA、30 g·L−1蔗糖和6 g·L−1琼脂，并接种赤桉无菌苗，25 ℃条件下16 h光照、8 h黑暗光周期培养3~4个月。低温处理组于4 ℃低温条件下处理24 h，对照组(CK)不进行处理。分别选取长势健壮无菌苗的茎尖组织，液氮速冻后磨碎，用Trizol(上海生工生物工程技术服务有限公司)提取总RNA。
1.2 Small RNA文库构建与测序
用Qubit 2.0 (Thermo, 美国)对总RNA进行定量，使用T4 RNA连接酶(NEB, 美国)将RNA 3′端和5′端与接头连接，随后用M-MuLV Reverse Transcriptase(NEB, 美国)合成cDNA单链，经PCR扩增后通过12 g·L−1聚丙烯酰胺凝胶电泳回收140~150 bp左右的PCR产物条带，通过质检后得到满足Illumina平台测序的文库，并委托上海生工生物工程技术服务有限公司进行高通量测序。
1.3 测序数据处理
测序下机的原始序列(Raw reads)经过去接头，去低质量、冗余序列等处理后得到纯净序列(Clean reads)，使用Bowtie(参数为-v 1)将纯净序列匹配到巨桉E. grandis基因组上，得到的序列再与Rfam14.1数据库对比，过滤其他非编码RNA。
1.4 miRNA的预测
以miRBase21.0库中植物和巨桉[12]的miRNA序列为参考源，使用blastn软件预测已知miRNA。具体参数为：-task blastn-short-word_size 16-evalue 0.01。将错配数≤1、e＜0.01的序列作为潜在miRNA。同时，将这些潜在的miRNA 用Bowtie(参数为 -v 2)匹配回巨桉基因组后，通过miREAP和RNAfold软件预测所有已知miRNA及其二级结构。未匹配上的序列，继续使用miRDeep2(v2.0.0.8)[13]软件预测新miRNA，采用默认参数。
对不同长度miRNA的首位碱基进行统计并绘制直方图，分析赤桉miRNA碱基的偏好性。本文除了miRNA二级结构采用RNAfold软件生成，其他所有图形均使用R语言绘制[14]。
1.5 差异表达miRNA的挖掘与分析
低温处理下赤桉差异表达的miRNA采用Bioconductor中的DEGSeq[15]包分析。首先将2组数据使用TPM进行归一化处理，归一化表达量=(Readcount
$ \times $ 1 000 000)/libsize，其中：Readcount为miRNA序列数目，libsize为样品所有miRNA 序列数目之和。设置DEGSeq的算法为MARS，进行差异分析。定义表达量变化2倍以上且P <0.001的为显著差异表达miRNA。1.6 miRNA靶基因预测与富集分析
使用psRNATarget[16]软件对miRNA进行靶基因预测，采用默认参数。同时对差异表达miRNA的靶基因分别进行Gene ontology[17] ( http://www.geneontology.org/)和KEGG[18] ( http://www.genome.jp/kegg/)富集分析。GO注释使用在线网站agriGO v2.0[19] ( http://systemsbiology.cau.edu.cn/agriGOv2/index.php)，并定义P<0.05为显著差异的类别。KEGG使用KOBAS3.0[20]软件，该软件以KEGG通路为单位，应用超几何检验进行分析：
$$ p=1-\displaystyle\sum\limits_{i=0}^{m-1}\dfrac{\left(\dfrac{M}{i}\right)\left(\dfrac{N-M}{n-i}\right)}{\left(\dfrac{N}{n}\right)},$$ (1) 式中，p为KEGG通路富集概率，i为求和中的特定项，N为所有基因中具有通路注释的基因数目，n为N中候选靶基因的数目，M为所有基因中注释为某特定通路的基因数目，m为注释为某特定通路的候选靶基因数目。同时，以拟南芥的低温胁迫响应网络和华南农业大学桉树育种课题组前期研究为参考，进一步挖掘涉及ICE1-CBFs-COR通路的miRNA。
2. 结果与分析
2.1 测序数据分析
对照组和低温处理组的小RNA高通量测序分别得到 29725524和31104120条原始序列，过滤掉低质量、两端接头和重复序列后，最后分别得到15813414和14456077条纯净序列。质控数据显示，低温处理组的Q30为96.47%，对照组的Q30为96.46%，结果可靠性高。此外，纯净序列的长度统计结果(图1)显示：2组小RNA长度的分布特征基本一致，以21 nt分布最多，其次是20和24 nt。
![]()
图 1 赤桉小RNA纯净序列的长度分布图Figure 1. Length distribution of clean reads of small RNA in Eucalyptus camaldulensis2.2 赤桉已知miRNA和新miRNA的预测
以RFam数据库对纯净序列进行注释，结果如表1所示，除了其他类别外，核糖体RNA(rRNA)占比最高，在对照和低温处理组中分别占比30.27%和26.08%，而核仁小RNA (snoRNA)占比最少，分别占比0.21%和0.17%。过滤掉rRNA、snRNA、snoRNA和tRNA序列后，以巨桉和miRBase21中所有植物的成熟miRNA为参考组，进行miRNA预测，结果如图2所示：2组数据一共预测到392个已知miRNAs(隶属于54个家族)，其中对照组282个，低温处理组329个，它们共同的有219个；预测到97个新miRNA，其中对照组65个，低温处理组51个，它们共同的有19个。在预测到的所有miRNA中，长度均以21 nt为主，其次，已知miRNA为20和22 nt，新miRNA为24 nt(图3)。图4展示了部分预测miRNA前体的二级结构。
表 1 赤桉小RNA分类统计Table 1. Classification statistics of small RNAs in Eucalyptus camaldulensis种类
Type对照 CK 低温处理 Low temperature treatment 数量 Count 占比/% Percentage 数量 Count 占比/% Percentage 核糖体RNA rRNA 4 787 452 30.27 3 769 781 26.08 核内小RNA snRNA 55 640 0.35 32 977 0.23 核仁小RNA snoRNA 33 402 0.21 24 413 0.17 转运RNA tRNA 828 159 5.24 368 918 2.55 其他 Other 10 108 761 63.93 10 259 988 70.97 总计 Total 15 813 414 100.00 14 456 077 100.00 ![]()
图 3 赤桉预测miRNA长度分布Figure 3. The length distribution of the predicted miRNAs in Eucalyptus camaldulensis![]()
图 4 赤桉部分预测miRNA前体二级结构A: eca-miR164b-3p, B: eca-miR390b-5p, C: eca-miR395f-3p, D: eca-miR-n40, E: eca-miR-n45, F: eca-miR-n51; A~C: 已知miRNA, D~F: 新miRNAFigure 4. The secondary structure of some predicted miRNAs in Eucalyptus camaldulensisA: eca-miR164b-3p, B: eca-miR390b-5p, C: eca-miR395f-3p, D: eca-miR-n40, E: eca-miR-n45, F: eca-miR-n51; A−C: Known miRNAs; D−F: Novel miRNAsmiRNA是由DCL1酶剪切得到的，DCL1酶在剪切pre-miRNA时，对其5′端第1个碱基具有强烈的偏好性，尤其体现在对尿嘧啶(U)的偏好上。miRNA首位碱基偏好性结果表明，对照和低温处理组均呈现出相似规律，具体为：长度为18 nt的miRNA其主要为腺嘌呤(A)，19~22 nt的miRNA其主要为尿嘧啶(U)，而23~26 nt的miRNA其以腺嘌呤(A)或鸟嘌呤(G)为主(图5)。
![]()
图 5 赤桉不同长度miRNA首位碱基偏好性Figure 5. Bias of the first base of different length miRNAs in Eucalyptus camaldulensis2.3 差异表达miRNA分析
以TPM标准化后的序列数为源数据，使用DEGSeq包分析对照组和低温处理组2个样品间差异表达的miRNA。共筛选出80个显著差异表达的miRNA(图6)，包括55个显著上调(含10个新miRNA)和25个显著下调的miRNA(含16个新miRNA)。上调的miRNA归属于miR395、miR167、miR399、miR8282、miR477、miR156、miR171和miR482(eca-miR482a-3p, eca-miR482d-5p)等8个家族，而下调的miRNA归属于miR862、miR482(eca-miR482b-5p)、miR6300、miR397和miR10506等5个家族。
![]()
图 6 赤桉差异表达miRNA火山图Figure 6. Volcanic map of differentially expressed miRNAs in Eucalyptus camaldulensis2.4 靶基因预测和基因富集分析
默认参数下，psRNATarget软件共预测到17 614个靶基因，其中低温胁迫下差异表达miRNA的靶基因有 5983个。这些差异表达miRNA的靶基因GO注释结果(图7)显示，定位到87个功能群中，包括分子功能(Molecular function, MF)49个、生物学过程(Biological process, BP)31个和细胞组成(Cell component)7个。这些靶基因中最多的一类是分子功能，主要集中在结合作用；在生物学过程中，主要集中在单一生物过程、生物调节和信号转导；在细胞成分中，靶基因主要与细胞膜有关。KEGG富集分析结果表明，这些靶基因参与的通路有23条，其中代谢通路(Metabolic pathway)和次生代谢物的生物合成(Biosynthesis of secondary metabolite)最多，分别占到了36.7%和17.8%(图8)。
![]()
图 7 赤桉差异表达miRNA靶基因的GO注释Figure 7. GO annotation of the target genes of differentially expressed miRNAs in Eucalyptus camaldulensis![]()
图 8 赤桉差异表达miRNA靶基因KEGG富集分析Figure 8. KEGG pathway analysis of the target genes of differentially expressed miRNAs in Eucalyptus camaldulensis2.5 赤桉ICE1-CBFs-COR通路相关miRNA的筛选
ICE1-CBFs-COR通路是植物低温胁迫应答的重要信号转导途径[21]，本试验进一步筛选到25个与该通路相关的miRNA(表2)。这些靶基因主要编码蛋白激酶、E3泛素化连接酶和转录因子，具体包含：ICE1，直接调控ICE1的MPK6(Mitogen-activated protein kinase 6)、OST1(Open stomata 1)、HOS1、SIZ1，与调控ICE1下游的CBF1~CBF3有关的BTF3(Basic transcription factor 3)、CAMTA3(Calmodulin-binding transcription activator 3)、BZR1(Brassinazole-resistant 1)、EIN3(Ethylene in sensitive 3)、PIF3(Phytochrome-interacting factor 3)和SOC1(Suppressor of constans overexpression 1)等。其中，仅eca-miR-n60有2个靶基因，其余均有1个靶基因；多个miRNA靶向同一个基因，比如eca-miR5780b-5p、eca-miR390b-5p、eca-miR-n51和eca-miR-n60均靶向OST1。这些miRNA长度主要为20~21 nt，包括8个保守型miRNA(miR156、miR159、miR164、miR171和miR482等家族)、10个非保守型miRNA以及7个新miRNA。此外，有3个miRNA(eca-miR-n51，eca-miR-n41和eca-miR862a-5p)经过4 ℃低温处理24 h后显著差异表达，其中，eca-miR-n51靶向OST1，而eca-miR-n41和eca-miR862a-5p均靶向EBF1，其余的miRNA则差异不显著。
表 2 赤桉ICE1-CBFs-COR通路相关miRNATable 2. The miRNAs associated with ICE1-CBFs-COR pathway in Eucalyptus camaldulensismiRNA1) 长度/nt
Length序列(5′→3′)
Sequence靶基因
Target gene基因ID
Gene ID蛋白质特征
Protein characteristiceca-miR-n33 21 ACGGAAUUGUUCGAGCCGACU ICE1 Eucgr.G01938 转录因子 Transcription factor eca-miR171g-3p 19 UGAGCCGGACCAAUAUCAC MPK6 Eucgr.L00026 蛋白激酶 Protein kinase eca-miR171j-3p 22 GAUGAGCCGGACCAAUAUCACG MPK6 Eucgr.L00026 蛋白激酶 Protein kinase eca-miR5780b-5p 23 UCCAGUCUCUGAUCAAUUUUGAC OST1 Eucgr.E00345 蛋白激酶 Protein kinase eca-miR390b-5p 21 GGCGCUAUCCCUCCUGAGCUU OST1 Eucgr.I00977 蛋白激酶 Protein kinase eca-miR-n51↓ 21 GAAUGUCUCCAAUCUGCCCGA OST1 Eucgr.H04745 蛋白激酶 Protein kinase eca-miR-n60 20 AGCUCAUCCAUCUGUAAGAG OST1 Eucgr.D02135 蛋白激酶 Protein kinase BZR1 Eucgr.H01239 转录因子 Transcription factor eca-miR156m-3p 20 UGCUCUCUCUCUUCUGUCAA BZR1 Eucgr.F01541 转录因子 Transcription factor eca-miR156o-3p 20 UGCUCUCUAUCUUCUGUCAA SOC1 Eucgr.A02846 转录因子 Transcription factor eca-miR156j-5p 21 UUGACAGAAGAGAGAGAGCAC SOC1 Eucgr.D02427 转录因子 Transcription factor 续表 2 Continued table 2 miRNA1) 长度/nt
Length序列(5′→3′)
Sequence靶基因
Target gene基因ID
Gene ID蛋白质特征
Protein characteristiceca-miR159k-3p 19 UUUGGAUUGAAUGGAGUCU SOC1 Eucgr.K00208 转录因子 Transcription factor eca-miR94a-3p 21 UCCCGGGAACAGAAUCAUUAC EIN3 Eucgr.J00631 转录因子 Transcription factor eca-miR845c-3p 20 CCUACAAUUGGUAUCAGAGC PIF3 Eucgr.B01825 转录因子 Transcription factor eca-miR-n38 21 AGGUGAAUUCUUAUAGAUCCA PIF3 Eucgr.B01825 转录因子 Transcription factor eca-miR482f-3p 21 UCUUUCCUAUUCCUCCAUUCC SIZ1 Eucgr.B02470 E3苏素化连接酶 E3 SUMO ligase eca-miR23a-5p 25 UGAGAGUGAGUGUAGAGUAGGGAAU HOS1 Eucgr.E00402 E3泛素化连接酶 E3 ubiquitin ligase eca-miR-n2 21 GCUCCCCAAACUGACUACCAA HOS1 Eucgr.E00402 E3泛素化连接酶 E3 ubiquitin ligase eca-miR-n41↓ 22 UCGGAAGUCUUUGAGGGAGAGA EBF1 Eucgr.C01723 E3泛素化连接酶 E3 ubiquitin ligase eca-miR862a-5p↓ 21 AGUUUCCUUGAAGACAUCCAA EBF1 Eucgr.C01524 E3泛素化连接酶 E3 ubiquitin ligase eca-miR845a-5p 20 AGCUCUGAUACCAAUUGUUG EBF1 Eucgr.C02778 E3泛素化连接酶 E3 ubiquitin ligase eca-miR396a-3p 21 AAGCUCAAGAAAGCUGUGGGA EBF1 Eucgr.C02778 E3泛素化连接酶 E3 ubiquitin ligase eca-miR7782a-5p 19 AGUGGUAUCAGAGCAGGUU BTF3 Eucgr.K02308 NAC蛋白β亚基 β-subunit of NAC protein) eca-miR7782b-5p 23 AGUGGUAUCAGAGCAGGUCGUCG BTF3 Eucgr.K02308 NAC蛋白β亚基 β-subunit of NAC protein) eca-miR827b-5p 22 UUUUGUUGAUGGCCAUCUAAUC CAMTA3 Eucgr.H04783 转录激活子 Transcription activator eca-miR164b-3p 20 UGGAGAAGCAGGGCACGUAA PhyB Eucgr.A00380 光感受器 Photoreceptor 1)“↓”表示在4 ℃低温处理24 h后显著下调
1)“↓” shows significant down-regulation after 4 ℃ low temperature treatment for 24 h3. 讨论与结论
本试验通过高通量测序共预测到隶属于54个家族的392个已知miRNA和97个新miRNA，并通过软件预测到17614个靶基因。与对照组相比，低温处理组在4 ℃处理24 h后，预测到80个显著差异表达的miRNA，包括55个上调miRNA和25个下调miRNA。GO基因功能注释表明，低温胁迫下miRNA的靶基因主要与结合作用、单一生物过程、生物调节和信号转导以及细胞膜有关；KEGG富集分析表明，差异表达miRNA靶基因可能在代谢通路和次生代谢物的生物合成中起关键作用。同时，本试验还分析了可能参与到赤桉ICE1-CBFs-COR通路有关的miRNA。
在miRNA预测过程中，发现Lin等[12]的巨桉成熟miRNA序列存在一些序列无法匹配到巨桉基因组的现象，因此我们将Lin等[12]的miRNA序列与其他植物的miRNA整合后，作为miRNA参考组，再与测序下机后的纯净序列匹配，得到潜在的已知miRNA。随后，将这些潜在已知miRNA匹配回巨桉基因组，再使用miREAP和RNAFold软件寻找前体序列。但是，由于赤桉和巨桉属于近缘种，基于巨桉基因组预测的miRNA有可能存在假阳性，后续仍需试验进一步验证。
在植物中，行使重要功能的miRNA，往往是保守的。本试验中，miR156、miR167、miR395和miR399家族的多个成员在低温处理后均显著差异表达。Zhou等[22]指出，miR156家族广泛响应植物低温胁迫，比如低温胁迫会诱导杨树miR156表达量下调[23-24]，而甘蔗抗寒品种‘FN39’在低温胁迫下miR156表达量上调[25]。研究发现，miR167和ARF(Auxin response factors)有关，并且miR167通过调控生长素水平抵御冷胁迫[26]。薄维平等[27]报道，低温胁迫会诱导木薯miR395的上调表达，其中miR395abcd的上调有利于减轻木薯的低温伤害。Gao等[28]发现在番茄中过表达拟南芥ath-miR399d可以提高番茄的耐寒性。除了这些保守型miRNA，本试验中其他显著差异表达miRNA在其他植物低温胁迫中也有过类似报道[29]。
ICE1作为关键性抗寒转录因子，其主导的ICE1-CBFs-COR途径是植物响应低温胁迫的重要通路。在该通路中，ICE1通过与MYB15[30]或HOS1[31]的互作来负调控CBFs及其下游基因，但其通过与SIZ1[32]或OST1[33]的互作来实现正调控。最近，Li等[34]报道MPK3/MPK6可与ICE1相互作用并将其磷酸化，降低ICE1的稳定性及其转录活性，从而负调控拟南芥CBF表达及其抗寒力。本试验进一步筛选到25个与该通路相关的miRNA，其中，eca-miR-n33靶向ICE1，可正调控ICE1-CBFs-COR通路的miRNA有靶向OST1的eca-miR5780b-5p、eca-miR390b-5p、eca-miR-n51和eca-miR-n60以及靶向SIZ1的eca-miR482f-3p，而负调控该通路的miRNA有靶向HOS1的eca-miR23a-5p和eca-miR-n2以及靶向MPK6的eca-miR171g-3p和eca-miR171j-3p。此外，Li等[35]报道BZR1可以正调控CBF1/2，增强植物抵御低温胁迫的能力。Jiang等[36]报道与光形态建成有关的转录因子PIF3可与CBF启动子直接结合，实现对低温胁迫的负调控；而EBF1(EIN3-BINDING F-BOX 1)可直接靶向PIF3进行26S蛋白酶体介导的降解，提高植物的抗寒力。综上可知，上述的这些miRNA可为后续深入研究其参与调控赤桉低温胁迫应答的分子网络奠定基础。
-
![]()
图 1 穿心莲连作土壤与对照土壤的PCoA (Unweighted Unifrac)分析
Figure 1. PCoA (Unweighted Unifrac) analysis of soil of Andrographis paniculata continuous cropping and control
![]()
图 2 穿心莲连作土壤与对照土壤细菌在门水平上的UPGMA聚类与柱状图组合分析
S1~S5为对照组，C1~C5为连作组；图中左侧是UPGMA 聚类树结构，右侧是各样品在门水平上的物种相对丰度分布图
Figure 2. UPGMA clustering and histogram combination analysis of bacteria in soil of Andrographis paniculata continuous cropping and control at the phylum level
S1−S5 belong to the control group, C1−C5 belong to the continuous cropping group; The left side of the figure is the UPGMA cluster tree structure, and the right side is the relative abundance distribution map of each sample at the phylum level
![]()
图 3 穿心莲连作土壤与对照土壤中相对丰度前30的细菌群落在属水平的热图
S1~S5为对照组，C1~C5为连作组；热图中着色方格对应的值为样品在属分类上的Z值，即为样品在该分类上的相对丰度和所有样品在该分类的平均相对丰度的差除以所有样品在该分类上的标准差所得到的值；方格颜色越红，说明该样品相对丰度越高
Figure 3. Heat map of bacterial communities with top 30 relative abundances in soil of Andrographis paniculata continuous cropping and control at the genus level
S1−S5 belong to the control group, C1−C5 belong to the continuous cropping group; The value corresponding to the colored square in the heatmap is the Z score of a sample in the genus classification, which is the difference between the relative abundance of the sample in the category and the average relative abundance of all samples in the category divided by the standard deviation of all samples in the category; The redder the color of the box, the higher the relative abundance of the sample
表 1 穿心莲连作土壤与对照土壤细菌的α多样性1)
Table 1 α diversity of bacteria in soil with continuous cropping of Andrographis paniculata and control soil
组别
GroupChao1指数
Chao1 indexObserved species指数
Observed species indexShannon指数
Shannon indexSimpson指数
Simpson index对照组 Control group 2798.60±72.31 1994.60±83.74 9.05±0.40 0.992±0.0064 连作组 Continuous cropping group 2546.67±67.61 1782.40±45.94 8.92±0.12 0.991±0.0022 P 0.000*** 0.001** 0.517 0.877 1)表中数据为平均值 ± 标准差；“**”和“***”分别表示差异达到0.01和0.001的显著水平(t检验)
1) Data in the table are means ± standard deviations; “**”and“***” indicate differences at 0.01 and 0.001 significance levels respectively (t test)
下载: 导出CSV
表 2 穿心莲连作土壤与对照土壤主要细菌门的相对丰度比较1)
Table 2 Comparison of relative abundance of major bacterial phylums in soil of Andrographispaniculata continuous cropping and control
门
Phylum相对丰度/%
Relative abundanceP 对照组
Control group连作组
Continuous cropping group变形菌门 Proteobacteria 22.53±3.40 27.48±2.42 0.029* 酸杆菌门 Acidobacteria 19.62±5.60 28.71±3.56 0.015* 拟杆菌门 Bacteroidetes 17.75±4.57 12.14±1.05 0.050 绿弯菌门 Chloroflexi 10.96±1.81 11.62±3.33 0.705 浮霉菌门 Planctomycetes 6.26±0.91 5.43±0.60 0.129 疣微菌门 Verrucomicrobia 4.50±2.52 4.24±0.74 0.834 放线菌门 Actinobacteria 3.18±0.65 2.86±0.23 0.325 芽单胞菌门 Gemmatimonadetes 3.53±0.93 2.48±0.17 0.064 蓝藻细菌门 Cyanobacteria 4.81±8.69 0.10±0.04 0.293 Patescibacteria 1.58±0.27 0.49±0.15 0.000*** 硝化螺旋菌门 Nitrospirae 1.09±0.26 0.89±0.17 0.184 1) 表中数据为平均值 ± 标准差；“*”和“***”分别表示差异达到0.05和0.001的显著水平(t检验)
1) Data in the table are means ± standard deviations; “*”and“***” indicate differences at 0.05 and 0.001 significance levels respectively (t test)
下载: 导出CSV
表 3 穿心莲连作土壤与对照土壤细菌主要细菌属相对丰度比较1)
Table 3 Comparison of relative abundance of major bacterial genera in soil of Andrographis paniculata continuous cropping and control
属
Genus相对丰度/%
Relative abundanceP 对照组
Control group连作组
Continuous cropping groupFlavobacterium 8.52±4.19 1.84±0.41 0.023* HSB_OF53-F07 1.94±0.50 1.81±0.49 0.672 Bryobacter 1.57±0.35 2.15±0.23 0.014* ADurb.Bin063-1 1.16±0.66 1.65±0.38 0.182 Candidatus_Solibacter 1.03±0.33 1.76±0.17 0.005** Burkholderia-Caballeronia-Paraburkholderia 0.24±0.05 2.17±0.13 0.000*** Nitrospira 1.08±0.26 0.88±0.17 0.172 1921-2 0.87±0.20 1.00±0.18 0.319 MND1 1.46±0.28 0.37±0.16 0.000*** Candidatus_Udaeobacter 1.07±0.63 0.68±0.11 0.216 Gemmatimonas 1.18±0.31 0.57±0.09 0.010* Mucilaginibacter 0.38±0.15 1.37±0.18 0.000*** Acidibacter 0.28±0.09 1.43±0.23 0.000*** Ellin6067 1.19±0.52 0.41±0.14 0.026* Haliangium 1.10±0.18 0.50±0.08 0.000*** 1) 表中数据为平均值±标准差；“*”、“**”和“***”分别表示差异达到0.05、0.01和0.001的显著水平(t检验)
1) Data in the table are means±standard deviations; “*”，“**”and“***”indicate differences at 0.05, 0.01 and 0.001 significance levels respectively (t test)
下载: 导出CSV
-
[1] 国家药典委员会. 中华人民共和国药典: 1部: 2015年版[M]. 北京: 中国医药科技出版社, 2015: 268. [2] LI J R, CHEN X Z, ZHAN R T, et al. Transcriptome profiling reveals metabolic alteration in Andrographis paniculatain response to continuous cropping[J]. Industrial Crops and Products, 2019, 137: 585-596. doi: 10.1016/j.indcrop.2019.05.067
[3] 王瑞, 董林林, 徐江, 等. 基于病虫害综合防治的人参连作障碍消减策略[J]. 中国中药杂志, 2016, 41(21): 3890-3896. [4] 蒋景龙, 余妙, 李丽, 等. 西洋参根腐病发生与根际土壤细菌群落结构变化关系研究[J]. 中草药, 2018, 49(18): 4399-4407. doi: 10.7501/j.issn.0253-2670.2018.18.027 [5] QIAO Y J, ZHANG J J, SHANG J H, et al. GC-MS-based identification and statistical analysis of liposoluble components in the rhizosphere soils of Panax notoginseng[J]. RSC Advances, 2019, 9(36): 20557-20564. doi: 10.1039/C9RA02110H
[6] 张亚琴, 陈雨, 雷飞益, 等. 药用植物化感自毒作用研究进展[J]. 中草药, 2018, 49(8): 1946-1956. doi: 10.7501/j.issn.0253-2670.2018.08.032 [7] 檀国印, 杨志玲, 袁志林, 等. 药用植物连作障碍及其防治途径研究进展[J]. 西北农林科技大学学报(自然科学版), 2012, 40(4): 197-204. [8] 曾令杰, 刘意, 禇晨亮, 等. 穿心莲化感作用与GAP栽培规范的研究[J]. 现代中药研究与实践, 2011, 25(3): 5-7. [9] 李俊仁, 陈秀珍, 刘凯频, 等. 穿心莲连作土壤提取物对其种子的化感作用[J]. 广州中医药大学学报, 2016, 33(3): 389-395. [10] 李玲梅, 李明. 不同处理方法对穿心莲种子萌发影响的初步研究[J]. 广东药学院学报, 2011, 27(4): 371-374. doi: 10.3969/j.issn.1006-8783.2011.04.010 [11] 黎韵琪, 李明, 唐堃, 等. 穿心莲营养体的化感自毒作用研究[J]. 北方园艺, 2014(21): 157-160. [12] ZHENG Y B, LI M. Autotoxicity of phenolic acids in root exudates of Andrographis paniculata (Burm. f.) Nees[J]. Allelopathy Journal, 2018, 45(2): 153-162. doi: 10.26651/allelo.j./2018-45-2-1183
[13] 李珍, 陈义三, 陈荣珠, 等. 生物炭对连作穿心莲根际土壤微生物的影响[J]. 上饶师范学院学报, 2019, 39(6): 50-55. [14] BRUMFIELD K D, HUQ A, COLWELL R R, et al. Microbial resolution of whole genome shotgun and 16S amplicon metagenomic sequencing using publicly available NEON data[J]. PLoS One, 2020, 15(2): e0228899. doi: 10.1371/journal.pone.0228899.eColletion2020.
[15] MENDES R, GARBEVA P, RAAIJMAKERS J M. The rhizosphere microbiome: Significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms[J]. FEMS Microbiology Reviews, 2013, 37(5): 634-663. doi: 10.1111/1574-6976.12028
[16] 沈菊培, 张丽梅, 贺纪正. 几种农田土壤中古菌、泉古菌和细菌的数量分布特征[J]. 应用生态学报, 2011, 22(11): 2996-3002. [17] 胡江春, 薛德林, 王书锦, 等. 大豆连作障碍研究Ⅲ. 海洋放线菌MB-97促进连作大豆增产机理[J]. 应用生态学报, 2002, 13(9): 1095-1098. doi: 10.3321/j.issn:1001-9332.2002.09.010 [18] 寇智瑞, 周鑫斌. 不同连作年限黄壤烟田土壤细菌群落的差异[J]. 植物营养与肥料学报, 2020, 26(3): 511-521. [19] 谭勇, 崔尹赡, 季秀玲, 等. 三七种植前后土壤细菌群落结构与多样性分析[J]. 昆明理工大学学报(自然科学版), 2016, 41(6): 92-99. [20] 张庚. 穿心莲(Androyraphis paniculata(Brum. f)Nees)连作的自毒作用研究[D]. 广州: 华南农业大学, 2010. [21] 母容, 潘开文, 王进闯, 等. 阿魏酸、对羟基苯甲酸及其混合液对土壤氮及相关微生物的影响[J]. 生态学报, 2011, 31(3): 793-800. [22] WANG Y Z, XU X M, LIU T M, et al. Analysis of bacterial and fungal communities in continuous-cropping ramie (Boehmerianivea L. Gaud) fields in different areas in China[J]. Scientific Reports, 2020, 10(1): 3264. doi: 10.1038/s41598-020-58608-0
[23] 张重义, 陈慧, 杨艳会, 等. 连作对地黄根际土壤细菌群落多样性的影响[J]. 应用生态学报, 2010, 21(11): 2843-2848. [24] JOHNSTON-MONJE D, LUNDBERG D S, LAZAROVITS G, et al. Bacterial populations in juvenile maize rhizospheres originate from both seed and soil[J]. Plant and Soil, 2016, 405(1/2): 337-355.
[25] 何玮, 郭琳微, 樊鹏辉, 等. 黄花棘豆在腐解过程中的化感作用及其土壤细菌群落结构分析[J]. 草业学报, 2015, 24(7): 21-29. doi: 10.11686/cyxb2014504 [26] FREY B, RIME T, PHILLIPS M, et al. Microbial diversity in European alpine permafrost and active layers[J]. FEMS Microbiology Ecology, 2016, 92(3): fiw018. doi: 10.1093/femsec/fiw018.
[27] BHATTACHARYYA P N, JHA D K. Plant growth-promoting rhizobacteria (PGPR): Emergencein agriculture[J]. World Journal of Microbiology and Biotechnology, 2012, 28(4): 1327-1350. doi: 10.1007/s11274-011-0979-9
[28] 陈三凤, 李季伦, 裘维蕃. 关于抑制植物病原真菌几丁质酶来源及效应的研究强作用黄杆菌的分离和鉴定[J]. 植物病理学报, 1992(4): 37-41. [29] KUNDIM B A, ITOU Y, SAKAGAMI Y, et al. New haliangicin isomers, potent antifungal metabolites produced by a marine myxobacterium[J]. The Journal of Antibiotics, 2003, 56(7): 630-638. doi: 10.7164/antibiotics.56.630
[30] 张琦梦. 伯克氏菌属植物病原菌DNA条形码检测技术研究[D]. 南京: 南京农业大学, 2016. [31] ZHAO J, ZHANG D, YANG Y Q, et al. Dissecting the effect of continuous cropping of potato on soil bacterial communities as revealed by high-throughput sequencing[J]. PLoS One, 2020, 15(5): e0233356. doi: 10.1371/journal.pone.0233356.
[32] DU L, HUANG B, DU N, et al. Effects of garlic/cucumber relay intercropping on soil enzyme activities and the microbial environment in continuous cropping[J]. HortScience, 2017, 52(1): 78-84. doi: 10.21273/HORTSCI11442-16
[33] 王悦, 杨贝贝, 王浩, 等. 不同种植模式下丹参根际土壤微生物群落结构变化[J]. 生态学报, 2019, 39(13): 4832-4843. [34] 王彩云, 武春成, 曹霞, 等. 生物炭对温室黄瓜不同连作年限土壤养分和微生物群落多样性的影响[J]. 应用生态学报, 2019, 30(4): 1359-1366. -
期刊类型引用(2)
1. 周文灵，陈迪文，吴启华，方界群，敖俊华. 耕作措施对宿根甘蔗产量、碳排放及经济效益的影响. 中国生态农业学报(中英文). 2024(09): 1481-1491 . 
百度学术
2. 凤晓岗，汪永诚，刘书雅，马横宇，王寅，董召荣，车钊. 有机肥施用量对燕麦–青贮玉米产量和N_2O排放的影响. 草业科学. 2024(11): 2675-2684 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 1133
- HTML全文浏览量: 28
- PDF下载量: 2538
- 被引次数: 6
