DOI: 10.3724/SP.J.1006.2022.11043
长链非编码RNA (LncRNA)在印度梨形孢促进大麦根部生长发育中的调控作用 郭楠楠刘天策史硕胡心亭牛亚丹李亮*
河北工业大学化工学院，天津300130
摘要：印度梨形孢定殖植物会促进植物生物产量提高，其分子机制有待于深入挖掘。LncRNA是一类长链非编码RNA，在植物生长发育过程中具有重要调控作用。然而，目前我们还不清楚大麦中LncRNA是否对印度梨形孢的定殖有响应。本研究发现，印度梨形孢定殖大麦会诱导大麦根系迅速发育，促进较多根分枝。采用全基因组高通量测序RNA-seq和生物信息学方法鉴定LncRNA，发现在P. indica定殖后3 d和7 d分别有752和932个差异表达的LncRNA，7 d相对于3 d有70个差异表达的LncRNA。其中在P. indica定殖后3 d有375个LncRNA表达上调，377个LncRNA表达下调；在P. indica定殖后7 d有459个LncRNA表达上调，473个LncRNA表达下调；7 d相对于3 d组中，有39个LncRNA表达上调，31个LncRNA表达下调。qPCR验证LncRNA的表达与RNA-seq结果一致。GO和KEGG分析表明，在P. indica定殖大麦促生过程中，部分LncRNA参与了激素信号途径的转录调控。该工作对于进一步理解LncRNA与靶基因的相互作用以及其对靶基因的调控功能提供了新的理论基础和实验依据，并为以LncRNA为靶点，进行作物性状改良提供新的思路和方向。 t型铝型材
关键词：LncRNA；大麦；印度梨形孢；RNA-seq；转录因子；细胞周期
Regulation of long non-coding RNA (LncRNA) in barley roots in response to Piriformospora indica colonization
GUO Nan-Nan, LIU Tian-Ce, SHI Shuo, HU Xin-Ting, NIU Ya-Dan, and LI Liang *
College of Chemical Engineering and Technology, Hebei University of Technology, Tianjin 300130, China
Abstract: The molecular mechanism of biomass enhancement by Piriformospora indica in colonization plants needs to be further explored. LncRNA is a kind of long-chain non-coding RNA, which plays an important role in the regulation of plant growth and development. However, it remains unclear whether barley LncRNAs are responsive to Piriformospora indica colonization. It was found that barley roots exhibited fast development and large roots branched after P. indica colonization. Genome-wide high throughput RNA-seq and bioinformatical analysis showed that 752 and 932 differentially expressed LncRNAs were detected in responsive to P. indica at 3-day and 7-day after colonization, respectively. And 70 differentially expressed LncRNAs were found at 7-day compared to 3-day. Among these, 375 were up-regulated and 377 were down-regul ated after 3 days’ colonizati
on, and 459 were up-regulated and 473 were down-regulated after 7 days’ colonization, 39 were up-regulated and 31 were down-regulated in 7-day to 3-day comparison group. The qPCR results verified the validity of LncRNAs data in RNA-seq. GO and KEGG analysis indicated that a few LncRNAs might be involved in the molecular functions, cellular components, and biological processes upon P. indica colonization. This study provides a new theoretical basis and experimental basis for further understanding of the interaction
本研究由国家自然科学基金项目(31801948)，河北省重点研发计划项目(19226505D)和河北省自然科学基金项目(C2021202005)资助。
This study was supported by National Natural Science Foundation of China (31801948), the Key Research & Development Projects in Hebei Province (19226505D) and the Natural Science Foundation of Hebei Province (C2021202005).
*通信作者(Corresponding author): 李亮，E-mail: liangli@
第一作者: E-mail: 150********@
Received (收稿日期): 2021-04-13; Accepted (接受日期): 2021-10-19; Published online (网络出版日期): 2021-11-02.
URL: /kcms/detail/11.1809.S.20211101.1150.006.html
between LncRNAs and coding sequences and regulatory functional networks, and provides new ideas and directions for crop shape improvement based on LncRNAs.
Keywords: LncRNA; barley; Piriformospora indica; RNA-seq; transcription factor; cell cycle
研究人员曾经认为长链非编码RNA (LncRNA)是转录过程中产生的“无用物”，最近发现其在植物的诸多调控过程中发挥重要作用[1]。尽管大多数LncRNA功能未知，但已鉴定的LncRNA表现出组织特异性表达并参与蛋白定位以及调控相关疾病的特点[2-3]。顺式(或反式) LncRNA具有调节蛋白质活性或定位[4-5]的作用，并可作为亚细胞结构的组织框架[2-3,6]。此外，某些LncRNA可作为竞争性内源RNA (ceRNA)与小分子RNA (microRNA)配对[7-9]。
目前已对部分植物在遭受冷、热、干旱、盐和氮等一系列胁迫反应中的LncRNA进行了鉴定，如毛果杨[10]、蒺藜苜蓿[11]、棉花[12]、玉米[13]和小麦[14-15]等植物。此外，在对拟南芥的研究中发现，冷诱导的LncRNA COOLAIR在春化早期起作用，基因3′末端的反义转录机制可能是以条件/阶段依赖的方式调控基因转录，LncRNA COLDAIR可调节春化介导的表观遗传沉默从而应对冷胁迫[16-17]；LncRNA LDMAR在杂交稻中起着调节光周期敏感雄性不育的作用[18]；LncRNA DRIR能正向调控拟南芥对干旱和盐胁迫的响应[19]。另外，Zou等[20]和Karlik等[21]在大麦全基因组范围内发现了大约8000
个LncRNA，并揭示了它们在过量硼处理下的表达模式。此外，与LncRNA共表达编码转录本的功能注释揭示了离子转运、定位建立和叶对刺激响应的分子功能[22]。LncRNA HvCesA6调控大麦细胞壁的合成[23]。基于LncRNA-mRNA相互作用表明，LncRNA介导的调控丝氨酸/苏氨酸蛋白激酶SMG1信号转导途径可能与野生大麦XZ5的耐旱性有关[24]。Karlik[25]等报道了盐胁迫下大麦通过抑制LncRNAAK370814的表达反向促进AK372815PLR基因表达，并进一步明确了LncRNA-AK363461和LncRNAAAK370506在细胞中的定位。这些发现增进了我们对LncRNA作为抗旱重要调节因子的认识。然而，LncRNA在更广泛的环境条件下的生物学功能还需深入探索。
大功率控制器内生担子菌印度梨形孢(P. indica)属蜡壳耳目，梨形孢属[26]。P. indica与大麦、拟南芥、小麦、玉米、烟草、番茄[2,27-31]等多种植物形成了互惠菌根共生关系。P. indica具有促进植物生长、增强抗逆性和抗病性、促进营养物质和生物活性物质积累等多种功能[32-37]。分子生物学研究表明，接种P. indica可诱导与激素信号转导、细胞壁代谢、碳水化合物运输代谢和根系形成等相关基因表达[38-39]。然而，在P. indica促进大麦生长中，LncRNA的响应及调控机制仍未见有报道，需进行深入探索。
1 材料与方法
核桃剥壳机
1.1 植株生长
将春大麦品种福大麦1号(162粒)用70%乙醇表面杀菌3 min，然后在次氯酸钠(5%活性氯)中杀菌20 min。在无菌条件下，用无菌水(pH为3)清洗1次，无菌蒸馏水冲洗6次。将灭菌种子放在1/2 MS培养基的方形培养皿中，每个培养皿中种6粒种子，每个处理组(对照组即Mock、印度梨形孢定殖后3 d、印度梨形孢定殖后7 d)进行3个重复。在20℃的培养箱中培养，条件如下：光照8 h (荧光冷白，Toshiba FL40SSW/37，180 µmoL m‒2 s‒1光子通量密度)/黑暗16 h，22℃/18℃，相对湿度60%。为了保证所有根系都能在培养基表面生长，将含有种子的培养皿垂直放置。培养结束后，挑选每个处理组长势良好的大麦各36颗，对其根部进行取样，样品送天津诺禾致源生物信息科技有限公司进行测序。
1.2 植物根系接种印度梨形孢
P. indica在CM培养基上生长3~4周，可制备孢子悬浮液。在CM培养基中加入0.05%吐温-20的灭菌水，用涂布棒轻轻刮擦平板表面，收集孢子，悬浮液通过miracloth (中国鼎国)过滤以除去菌丝体。将悬浮液在1359 ×g下离心7 min收集孢子。用灭菌的吐温水清洗孢子3次。显微镜下，血球计数板测定孢子密度。将孢子用无菌吐温水稀释浓度至105个孢子mL‒1。接种时，用移液将3 mL孢子悬浮液滴到方形培养皿中的植物根表面。
1.3 RNA提取及测序定量与鉴定
发电机集电环用TRIzol (Invitgen，#15596-018)法将对照组和处理组的根(0.5 g)在液氮中研磨成粉末。在1%琼脂糖
凝胶上监测RNA降解和污染。分别使用NanoPhotometer分光光度计和Qubit 2.0 Flurometer中的Qubit RNA 试剂盒检测RNA纯度与浓度。纯化的RNA用Revert-Aid第一链cDNA合成试剂盒进行反转录。根据制造商的说明，使用Ribo-Zeror RNA Removal Kit (Epicentre，BioTechnoLogies)和NEB Next RULtraTMRNA Library Prep Kit for Illumina (New EngLand BioLabs)构建RNA-Seq文库。根据末端配对法，使用Illumina HiSeqTM2000对得到的文库进行测序。每个样本设置3个生物学重复。
1.4 差异表达基因的GO和KEGG富集分析
GO (Gene Ontology)是描述基因功能的综合性数据库，可分为分子功能(Molecular Function)，生物过程(biological process)和细胞组成(cellular component) 3个部分。利用clusterProfiler R软件包与GO网站(/)对差异表达的LncRNA进行GO富集分析。KEGG (Kyoto Encyclopedia of Genes and Genomes)是一种数据库资源，用于从分子水平信息，尤其是基因组测序等高通量实验技术产生的大规模分子数据集(jp/kegg/)了解生物系统(如细胞、生物体和生态系统)的高级功能和实用性。使用clusterProfiler R软件包测试KEGG通路中差异表达的LncRNA的统计富集情况。GO富集和KEGG通路富集以p_value或p_adj小于0.05为显著富集。
1.5 实时定量聚合酶链反应(qRT-PCR)分析
分别用对照和处理组的总RNA (1 µg)与M-MLV Reverse转录酶来制备cDNA。经DNase Ⅰ处理后，将c
DNA作为qRT-PCR的模板，使用LncRNA特异性引物和靶mRNA特异性引物对所选择的LncRNA和mRNA进行定量分析。定量RT-PCR检测各基因的表达水平(SYBR绿荧光法测定)。qPCR实验在Light Cyker96快速实时PCR系统(Roche)上进行，该反应溶液含有2×ULtra SYBR 混合物10 µL，cDNA模板100 ng ，正反向引物10 µmol L‒1。以HvUBIQUTIN为对照，进行3次重复。扩增程序按以下步骤进行：初始活化步骤在95℃下进行5 min，然后进行30个循环(95℃ 20 s，56℃ 35 s，72℃ 35 s，65℃ 20 s)。在每个循环结束时，分别测定熔融曲线，以保证单个PCR产物的扩增。
2 结果与分析
2.1 P. indica的定殖提高了大麦根的产量全自动文具盒
大麦幼苗在P. indica处理的第3 d和第7 d的根长均比对照组显著增长(图1-A、B)。对3 d和7 d处理组与Mock组相比，根长度分别增加了75%和58%。此外，根在P. indica定殖7 d后分枝显著，根数在3 d 和7 d比Mock组分别增加了25%和42.85% (图1-C, D)。该结果表明，P. indica定殖能促进大麦根系生长和根数增加。
2.2 大麦根系对P.indica定殖反应的差异表达的LncRNA谱
为揭示大麦LncRNA对P. indica定殖的分子响应机制，从12个大麦根部样本中提取RNA，并进行全转
录组序列测定。LncRNA鉴定的高通量测序流程和结果见图2。每个样本的RNA序列产生了1.042~1.285亿次raw reads和1.038~1.252亿次clean reads。Raw reads提交至NCBI (PRJNA71930)。利用Hisat将组装好的clean reads映射到大麦基因组上，所有样本的平均映射率为74.13% (表S1)。利用RSeQC测定RNA-seq 中LncRNA的饱和度见表S2。
图1 大麦在P. indica定殖后3 d和7 d形态、根长与根数统计结果图
Fig. 1 Barley root physio morphology, root length and root number at 3 day and 7 day after P. indica colonization
A：处理组与对照组在3 d的根长和根数，Scale bar = 1.0 cm；B：处理组与对照组在3 d根的分枝情况，标尺为1.0 cm；C：处理组与对照组在3 d和7 d的根长统计；D：处理组与对照组在3 d和7 d的根数统计。进行3个生物学重复实验和统计分析。*差异显著(P < 0.05)。A: root length and root number were compared between treatment and Mock at 3 days. Bar = 1.0 cm; B: branched roots were compared between treatment and Mock at 3 days. Bar: 1.0 cm; C: chart showing the statistical results of the root length of barley at 3 days and 7 days; D: chart showing the statistical results of the root number of barley at 3 days and 7 days. The experiments and statistical analysis were performed in three biological
replicates. *: P < 0.05.
图2 LncRNA鉴定的高通量测序流程图
Fig. 2 High-throughput sequencing flowchart of LncRNA identification
层次聚类热图(图3)表明3个处理组(即3 d vs Mock、7 d vs Mock和7 d vs 3 d)的差异表达的LncRNA 表达丰度不同，说明这些基因位于不同的代谢途径及信号通路中。由图4可知，在3 d vs Mock组中，共有752个差异表达的LncRNA，其中，有375个LncRNA表达上调，377个LncRNA表达下调；在7 d vs Mock
组中，共有932个差异表达的LncRNA，其中上调表达的LncRNA有459个，下调表达的LncRNA有473
个；在7 d vs 3 d组中，共有70个差异表达的LncRNA，其中39个LncRNA表达上调，31个LncRNA表达下调。
按摩毯
图3 差异表达的LncRNA的层次聚类热图
Fig. 3 Hierarchical clustering heat map of differentially expressed LncRNA
横坐标为样本名称，纵坐标为差异表达的LncRNA，左侧根据表达相似程度对基因进行聚类，上方根据表达谱的相似程度对每个样本进行聚类，由蓝至红表达量逐渐上调，数字为均一化后的相对表达量。
The abscissa is the sample, the ordinate is the differentially expressed LncRNA. On the left side, the genes are clustered according to the degree of similarity of expression, and the upper part of the sample is clustered according to the degree of similarity of the expression profile. The relative expression level is gradually up-regulated from blue to red, and the number is the relative expression after standardization.
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议