-
长链非编码 RNA (lncRNA)是一类长度大于200个核苷酸的调控RNA，不具有编码蛋白质能力[1−2]。根据 lncRNA 在基因组上的位置，通常被分为长内含子非编码 RNA (intronic lncRNA)、长基因间非编码RNA (lincRNAs) 和长非编码自然反义转录本 (lncNATs)[3−4]。lncRNA 具有顺式调控和反式调控2种调控模式。顺式调控是指 lncRNA 对距离较近的 mRNA 的一种转录激活与表达调控方式。lncRNA 主要以序列互补为原理结合在未解链的染色质大沟中，或者解链的 DNA 链上。反式调控是指 lncRNA 通过与启动子和增强子或与这些位点结合并受影响的其他蛋白质相互作用来调节非连锁基因的表达、染色质状态和 RNA 聚合酶活性[5]。
lncRNA 可以起到支架的作用，组织蛋白复合物，充当亚细胞结构，从而行使功能。lncRNA 以剂量依赖形式在细胞质或者细胞核内调控结合的蛋白或者 RNA 的活性[6]。lncRNA 还可以通过与 RNA 或 DNA 的序列互补性或同源性，在结构上形成大分子复合物组装的分子框架和支架[7]。lncRNA 可作为与其互作分子的招募者、引导者和信号分子等，通过表观遗传调控、转录调控、转录后调控和翻译调控而发挥其功能[8−12]。研究表明：lncRNA 在毛竹Phyllostachys edulis次生壁的合成过程中发挥重要作用[13]。
转座子是生物体的重要组成成分。研究表明：许多 lncRNA 来源于转座子序列，称为转座子来源的lncRNA (TE-lncRNA1)[14−15]。23%的拟南芥Arabidopsis thaliana、50%的水稻Oryza sativa、51%的玉米Zea mays的 TE-lncRNA 主要来源于逆转座子序列[16−17]。它主要参与作物的生长发育、非生物胁迫应答等不同的生物过程[18]。转座子来源的 lncRNA 受胁迫诱导表达，且表达量较高，具有较高的研究价值[19]。
毛竹分布范围广，种植面积大，是中国极为重要的非木材森林资源[20]。吴佳军[21]构建了低温、高温、紫外、高盐等4种胁迫毛竹全转录组数据库。在本研究中，我们利用4种胁迫全转录组数据库筛选并鉴定了在紫外胁迫下差异表达的1个lncRNA，它来源于 LARD 转座子序列，命名为Ph-TElncRNA1。本研究以毛竹Ph-TElncRNA1为研究对象，鉴定其调控的靶基因，分析了其与靶基因的表达协同性，以及Ph-TElncRNA1和靶基因在叶片着色过程中的表达趋势，初步探究了Ph-TElncRNA1的功能。
-
毛竹种子采自广西省桂林市灵川县同一母竹，发芽后放置在室温人工气候室中进行培养。紫外胁迫处理：取五叶一心期实生苗，30 W紫外灯分别处理2、4、6、8 h[22]。对照为正常的日常光。每个处理做3个重复，处理后立即取样，每个样取5张叶片和心。处理好的样品进行液氮速冻，存于−80 ℃冰箱，用于RNA提取。
不同着色程度毛竹叶片采集：选取长势良好且一致的毛竹实生幼苗，待其在暗培养箱里长到三叶一心时，取黄色叶片的幼苗，并将其中3盆转移到光照培养箱中；培养第2天取淡黄色叶片；培养第4天，取淡绿色叶片；培养第6天，取绿色的叶片(图1)。每个类型的叶片取3张，进行液氮速冻，存于−80 ℃冰箱，用于RNA提取。
Figure 1. Leaves with yellow, light yellow, light green and green color
-
以前期转录组测序后的低温、高温、紫外、高盐等4种胁迫全转录组为基础，在吴佳军[21]构建的毛竹 lncRNA 数据库中，筛选到在紫外胁迫处理中差异表达，且靶基因差异表达，缺乏开放阅读框的 lncRNA。使用以下4个的软件：CNCI (CodingNon-CodingIndex)、CPC2 (Coding Potential Calculator)、PLEK (predictor of long non-coding RNAs and messenger RNAs based on an improved k-mer scheme)和Pfam (PfamScan) 判断转录本的编码情况，确定典型的 lncRNA。
根据其在基因组的位置和转座子位置关系，确定其为来源转座子的 lncRNA，命名为Ph-TElncRNA1。
在毛竹基因组中提取Ph-TElncRNA1的序列，进行引物设计(F：AGGGGCGCTTCCGGTTCTG； R：AAGAGATATATTTTACATCCCAGA)，PCR 扩增，测序确认。
-
使用2种方法预测lncRNA的靶基因。第1种是顺式 (cis) 靶基因，根据lncRNA与基因的位置关系，将lncRNA附近100 kb范围以内的基因作为其靶基因，通过使用perl脚本查找这类靶基因[23]。第2种是根据lncRNA和mRNA之间的碱基互补配对关系，利用在线软件LncTar计算lncRNA和mRNA配对位点的自由能和标准化自由能，将低于标准化自由能阈值的基因认定为lncRNA的靶基因，称作反式 (trans) 靶基因。通过转录组数据筛选与对应TE-lncRNA变化趋势一致的靶基因，且差异表达的基因作为研究对象。
-
提取经胁迫处理后毛竹实生苗RNA，以反转录后的cDNA作为定量模板，进行实时荧光定量PCR (RT-qPCR)验证，分析Ph-TElncRNA1和靶基因表达趋势的关联性，剔除相关性不强的靶基因。使用IBM SPSS Statistics软件对Ph-TElncRNA1和靶基因的相对表达量进行皮尔逊相关性分析，确定Ph-TElncRNA1和靶基因的相关情况。
-
提取4种颜色(黄色、淡黄、淡绿、绿色)毛竹实生苗叶片(图1)RNA，以反转录后的cDNA作为定量模板，进行RT-qPCR验证，分析Ph-TElncRNA1和靶基因psbA在叶片着色过程中的表达模式。
-
以吴佳军[21]构建的 lncRNA 数据库为基础，筛选到紫外胁迫处理中差异表达，且靶基因差异表达的 lncRNA。通过 CNCI、Pfam、CPC2、PLEK 等4个软件分析lncRNA 的真实性。 CNCI 和Pfam软件均表明筛选的 lncRNA 不具有编码蛋白的能力，CPC2[编码可能性(coding probability)为0.144 099，菲克特分数(Fickettscore)为0.465880)]和 PLEK (软件评分为−2.282730，该软件评分小于 0 为非编码 RNA)也表明筛选的 lncRNA 编码蛋白的可能性较低。由此可以判定筛选的 lncRNA 是一个典型的非编码 RNA。
分析表明：该lncRNA 位于一个反转录转座子大片段(large retrotransposons derivatives, LARD)内部(图2)，为转座子来源的lncRNA，命名为Ph-TElncRNA1。Ph-TElncRNA1全长为 342 bp，其中 185 个碱基来源于外显子，157 个碱基来源于内含子(图3)。
Figure 2. 2Positional relationship between the LARD transposon and Ph-TElncRNA1
Figure 3. Composition and sequence of Ph-TElncRNA1
利用在线软件LncTar计算lncRNA和mRNA配对位点的自由能和标准化自由能，低于标准化自由能的靶基因认定为是该 lncRNA 的靶基因。进一步通过转录组数据筛选与对应Ph-TElncRNA1变化趋势一致的靶基因，且差异表达的基因作为研究对象。得到该 lncRNA 的3个靶基因，分别为PH010-02523G0060、PH01002523G0140 和 PH010025-23G0050。
以紫外处理2 h的毛竹叶片 cDNA 为模板，通过 PCR 扩增出相应的片段(图4)，测序结果表明成功克隆到Ph-TElncRNA1 (图5)。
Figure 4. Electropherogram of Ph-TElncRNA1
Figure 5. Alignment of sequencing results of Ph-TElncRNA1
-
通过实时荧光定量 PCR 验证，发现 PH01002523G0060 和PH01002523G0050 在紫外胁迫处理 2、4、6、8 h 后，相对表达量的趋势杂乱，与Ph-TElncRNA1不一致，说明 PH01002523G0060 和 PH01002523G0050 这2个靶基因与Ph-TElncRNA1的关联性不强，从而剔除掉。最终得到靶基因 PH01002523G0140 与Ph-TElncRNA1关联性较强，且相对表达量呈现完全一致的趋势。靶基因 PH01002523G0140 是psbA同源基因，编码 photosystem Ⅱ protein D1。与对照相比，Ph-TElncRNA1和靶基因psbA的相对表达量在紫外胁迫处理初期表达量上升显著，随着紫外灯照射时间的延长逐渐降低。Ph-TElncRNA1和psbA在紫外胁迫处理 2 h 时，相对表达量达到最高值(图6)，Ph-TElncRNA1与psbA呈显著正相关(P=0.000，皮尔逊相关系数为0.898)。
Figure 6. Relative expression of Ph-TElncRNA1 and psbA under UV treatment
-
Ph-TElncRNA1的靶基因psbA编码光系统Ⅱ (PSⅡ) 反应中心蛋白 D1，是叶绿体的重要组成部分。通过对叶片着色的4个时期样本(黄色、淡黄、淡绿、绿色)进行RT-qPCR发现：Ph-TElncRNA1和靶基因psbA在4个时期的相对表达量呈正相关(图7)。根据图7中的表达趋势可知：在淡黄时期，也就是将毛竹放入培养箱2 d时，表达量最高，在这一时期叶绿体正在形成，而光系统Ⅱ (PSⅡ) 反应中心蛋白 D1作为重要的蛋白亚基，也处于表达高峰时期。Ph-TElncRNA1和靶基因psbA在叶绿体形成过程中表达量最高，在叶绿体形成后呈现逐渐下降的趋势。两者在叶片着色过程中，相对表达量的变化趋势呈显著正相关(P=0.000，皮尔逊相关系数为0.926)。
Figure 7. Relative expression of Ph-TElncRNA1 and psbA in four periods of leaves colorization
-
由于lncRNA具有低表达和低序列保守性的特点，曾被认为是RNA聚合酶Ⅱ转录的副产物[24]。随着人们对lncRNA研究的不断深入，逐渐发现其在动植物等真核生物生长发育的过程中扮演着至关重要的作用[25]。目前，对lncRNA的研究大多集中于治疗人体各类疾病的研究上，如胃癌[26]、口腔鳞状细胞癌[27]、变应性鼻炎[28]等。相比于动物中lncRNA的研究，对植物lncRNA的探索仍处于生物信息学分析的初级阶段，研究主要集中在拟南芥[29]、烟草Nicotiana tabacum[30]及水稻[31]等植物中。研究表明：大量非编码RNA来源于转座子序列，如lncRNA、miRNA、siRNA等。转座子来源的lncRNA参与了真核生物生长发育等活动中[16]，但大多数转座子来源的lncRNA的功能尚未研究。
从实验室前期转录组数据分析可知：毛竹63.59%的lncRNA是转座子起源的，毛竹中受非生物胁迫诱导的 TE-ncRNA 非常丰富[21]。本研究通过对Ph-TElncRNA1和psbA在紫外胁迫下相对表达量的分析发现：Ph-TElncRNA1与psbA的相对表达量呈完全相同的趋势，在紫外胁迫下这两者是协同表达的。Ph-TElncRNA1与psbA相对表达量显著相关，进一步说明psbA可能是Ph-TElncRNA1靶基因。
psbA基因编码光系统Ⅱ(PSⅡ)反应中心蛋白D1[32]。D1蛋白是PSⅡ反应中心复合物中的一个非常重要的蛋白亚基，D1蛋白在叶绿体中的合成效率决定PSⅡ复合体的修复效率，进而决定植物的光合作用效率。研究表明：光可以诱导D1蛋白的表达，实现精准控制[33]。本研究对毛竹实生苗叶片着色过程Ph-TElncRNA1和psbA相对表达量的分析表明：Ph-TElncRNA1与psbA在叶绿体正在形成时期达到表达峰值，呈正相关，说明lncRNA可以通过靶基因参与调控毛竹叶片的着色过程。本研究为进一步深入分析lncRNA与靶基因的互作机制奠定基础，对于研究lncRNA的生物学功能也有一定的意义。
Identification of Ph-TElncRNA1 in Phyllostachys edulis and its regulation of target genes
doi: 10.11833/j.issn.2095-0756.20220396
- Received Date: 2022-06-14
- Accepted Date: 2022-10-17
- Rev Recd Date: 2022-10-10
- Available Online: 2023-04-03
- Publish Date: 2023-04-20
-
Key words:
- Ph-TElncRNA1 /
- leaf coloration /
- D1 protein /
- Phyllostachys edulis
Abstract:
| Citation: | ZHAO Jiamin, YU Lu, DING Yiqian, ZHOU Mingbing. Identification of Ph-TElncRNA1 in Phyllostachys edulis and its regulation of target genes[J]. Journal of Zhejiang A&F University, 2023, 40(2): 314-320. doi: 10.11833/j.issn.2095-0756.20220396
|
DownLoad: