鸡骨骼肌中天然反义lncRNA VGLL2-AS的鉴定及其与VGLL2的关系研究

李文雅，牛欣然，任团辉，蔡含芳，韩瑞丽，田亚东，刘小军，康相涛，李转见

(河南农业大学动物科技学院，郑州 450046)

生物体内在自然状态下形成的，能够与其正义链RNA互补的转录产物叫做自然反义链转录本(NATs)。以往对NATs的研究结果表明，NATs大多分布在其正义链基因的启动子区域与外显子区域，既包含了具有编码功能的mRNA，也包含了非编码RNA(如lncRNA)，根据转录的方向不同，又可以分为cis-NATs和trans-NATs[1]。

早期有研究者认为NATs属于一种转录异常副产物，因而其本身没有调控任何基因表达的功能[2]，近年来的研究则推翻了这一观点。研究发现，NATs能通过与其正义链转录本序列以重叠或互补配对，以此而实现可以在包括转录或者转录后等多个水平参与到对正义链转录本基因的表达调控[2-3]。已逐步证实NATs广泛的存在于人类[4-5]、小鼠[6]、猪[7]、鸡[8]等多个常见物种的转录组中，暗示着它们可能在参与这些物种的基因表达以及在调控功能等方面具有重要的作用。而与人类和小鼠NATs的研究相比，在家禽中相关的研究相对滞后。

鸡的VGLL2基因(vestigial like family member 2,VGLL2)位于3号染色体，在GenBank数据库中，鸡VGLL2 mRNA有4个外显子，而X1、X2和X3均因内含子3的部分滞留而多形成一个外显子。在分析的所有脊椎动物物种中，VGLL2表达主要定位在早期胚胎及以后的骨骼肌的体细胞分裂中，这与在C2C12成肌细胞的肌肉分化期间观察到VGLL2表达的上调相一致[9-10]。这些数据表明，VGLL2在骨骼肌生长发育过程中可能发挥着非常重要的作用，但VGLL2在家禽中的表达和调控还尚不清楚。

本课题组基于前期对6 周龄地方鸡(固始鸡)和商业肉鸡(AA鸡)肌肉组织的转录组测序(RNA-Seq)鉴定出一个在鸡肌肉组织中特异性表达，由VGLL2基因序列逆向转录而来的lncRNA，命名为VGLL2-AS(GenBank登录号: KY126094)，猜测VGLL2-AS可能在鸡的肌肉发育中发挥重要作用，但具体机制尚不清楚[11]。为此，本研究对VGLL2-AS的生物学功能展开进一步的探索，以期为NATs调控家禽骨骼肌生长提供理论基础。

1 材料与方法

1.1 试验材料

本试验所用动物均来自于河南农业大学家禽种质资源场，包括11胚龄固始鸡胚蛋若干枚、1日龄(1 d)、6周龄(6 W)、16周龄(16 W)、22周龄(22 W)、30周龄(30 W)固始鸡各6只。

本试验所用的核酸染料(DNA Green)以及2×Taq PCR Master Mix分别购自TIANDZ公司和CWBIO公司；试验用到的DNA Marker，包括Ladder 100 和D2000，分别购自上海莱枫以及北京中科瑞泰生物科技有限公司；琼脂糖凝胶回收试剂盒、Trizol、RNA反转录试剂盒、SYBR Premix和DH5α感受态细胞均购自日本TaKaRa公司；试验中用到的有机溶剂包括无水乙醇、三氯甲烷等均购自天津市富宇精细化工有限公司；PBS、PRMI-1640细胞培养基和质粒小提试剂盒购自Promega公司；细胞培养试验中用到的细胞消化液(0.25% EDTA)、青链霉素混合液和胎牛血清(FBS)购自于Gibco公司；放线菌素D和PARISTM Kit分别购自于Sigma和Invitrogen公司。

1.2 引物设计

根据NCBI在线网站所提供的基因序列，利用NCBI在线设计引物网站(https://www.ncbi.nlm.nih.gov/)和Sangon Biotech公司分别设计和合成表1所示的引物序列。

表1 本研究所使用的引物Table 1 The primers used in this study

1.3 RNA抽提、反转录和实时荧光定量PCR

利用组织(细胞)裂解液Trizol来提取不同时期固始鸡不同组织样品(包括肺、心、肾、肝、皮脂、腹脂、腿肌、胸肌和下丘脑)或者细胞样品的RNA，按照反转录试剂盒使用说明书对质量合格的RNA进行反转录，以反转录而来的cDNA为模板进行实时荧光定量PCR(real-time quantitative PCR, RT-qPCR)，具体操作及分析参考任团辉[11]的方法。

1.4 RT-PCR验证VGLL2-AS与VGLL2基因的双向转录

分别以VGLL2的正义链第3个外显子设计引物GF1和GR1，以VGLL2-AS的反义链第1个外显子设计引物DL-F和DL-R，GAPDH作为内参基因。引物序列见表1。

用于特异基因逆转录cDNA的合成步骤如下：以固始鸡胸肌总RNA(1 μg)为模板，去除模板中的基因组DNA污染后，于反转录反应中分别以DL-R和GAPDH-R或GR1和GAPDH-R特异引物进行反转录，得到的链特异cDNA为接下来Touchdown PCR扩增的模板。反应体系为：cDNA 1 μL，上、下游引物(10 μmol·L-1)各0.5 μL、ddH2O 3 μL，2×SYBR Green Mix 5 μL。反应程序为： 95 ℃ 4 min；95 ℃ 30 s，退火65 ℃～50 ℃(-1 ℃/cycle) 30 s，72 ℃ 40 s；95 ℃ 30 s，50 ℃ 30 s，72 ℃ 40 s，分别为15和26个循环。扩增所得产物于核酸染料预染的1.5%琼脂糖凝胶中进行电泳(120 V，110 mA，30 min)，最后在凝胶成像系统下成像。

1.5 细胞培养

鸡原代成肌细胞的分离和培养详见戴巍等[12]的研究。

1.6 核酸酶保护试验

向固始鸡胸肌总RNA(1 μg)中加入2 μL DNase Ⅰ和1 μL RNase A以消除基因组的DNA及单链RNA，用反转录试剂盒将消化产物反转录，反转录产物作为模板进行下一步的PCR检测。反应体系为：cDNA 1 μL，上、下游引物(10 μmol·L-1)各0.5 μL、ddH2O 3 μL，2×Taq PCR Master Mix 5 μL，反应条件同RT-qPCR。

1.7 放线菌素D添加试验分析VGLL2-AS与VGLL2半衰期

将鸡原代成肌细胞接种至12孔板中，待其细胞汇合度达到70%左右时，收取3个孔的细胞培养品，设置为0 h样品；将最终浓度为2 μg·mL-1的放线菌素D加入细胞中，之后每1 h收集3个孔细胞样品，提取细胞RNA后检测靶基因的表达。

1.8 鸡原代成肌细胞VGLL2和 VGLL2-AS 亚细胞定位

鸡成肌细胞细胞核与细胞质的分离使用PARISTM Kit试剂盒进行，细胞质和细胞核的内参基因分别为GAPDH和U6基因，通过RT-qPCR方法检测成肌细胞中VGLL2-AS和VGLL2基因的表达，进而确定两者的细胞定位。

1.9 统计分析与作图

本试验采用相对定量2-ΔΔCt计算方法对所得qRT-PCR的原始数据进行分析，以2-ΔΔCt值来代表不同基因在不同样品中的相对表达水平。2-ΔΔCt方法具体计算过程如下：

ΔCt= Ct(目的基因)-Ct(内参基因)；

ΔΔCt=ΔCt-ΔCt(对照组均值)；

目的基因在某一样品中的相对表达量=2-ΔΔCt。

本研究所有的定量数据均是通过罗氏实时荧光定量软件LightCycler 96进行收集，再经Excel等软件进行整理；用SPSS 20.0对整理后的数据进行显著性分析，采用单因素方差分析进行显著性检验，多重比较采用新复极差法(SSR)，设置显著性水平为0.05。最后用作图软件GraphPad Prism 6进行制图，数据表示的形式为“平均值±标准差”[10]。

2 结 果

2.1 鸡肌肉组织中VGLL2-AS的验证

为了验证VGLL2-AS是否真实存在，在VGLL2-AS的特异区(与VGLL2相比，VGLL2-AS缺少第3外显子，故选取VGLL2-AS跨第2外显子与第4外显子区域作为VGLL2-AS特异区)设计引物F1/R1并进行PCR扩增和测序检测。结果如图1所示，特异区PCR产物测序结果与目的基因VGLL2-AS序列完全一致，证明VGLL2-AS在鸡的转录组中真实存在。

A. VGLL2-AS的琼脂糖凝胶电泳图；B. VGLL2-AS的测序图A. Agarose gel electrophoresis of VGLL2-AS; B. DNA sequencing map of VGLL2-AS图1 VGLL2-AS的验证Fig.1 Verification of VGLL2-AS

2.2 VGLL2与VGLL2-AS均在肌肉中特异表达

为了验证VGLL2与VGLL2-AS是否在肌肉中高表达，qRT-PCR检测VGLL2与VGLL2-AS在固始鸡不同时期的肺、心、肾、肝、皮脂、腹脂、腿肌、胸肌和下丘脑组织之间的表达差异。结果如图2所示，可以看出VGLL2与VGLL2-AS表达趋势一致，即均在肌肉组织中特异高表达，而在其余组织中几乎不表达。且在6 w和16 w肌肉组织中的表达量与1 d相比，表达量均显著降低(P<0.05)，说明VGLL2和VGLL2-AS参与了鸡的肌肉发育进程。

A. VGLL2时空表达谱；B. VGLL2-AS的时空表达谱。经单因素方差分析检验，*表示差异显著(P<0.05)A. The spatio-temporal expression pattern of VGLL2; B. The spatio-temporal expression pattern of VGLL2-AS. Tested by one-way analysis of variance, * means significant difference(P<0.05)图2 VGLL2与VGLL2-AS的表达规律Fig.2 Expression patterns of VGLL2 and VGLL2-AS

2.3 VGLL2-AS与VGLL2可以进行双向转录

为了验证VGLL2-AS与VGLL2是否可以进行双向转录，通过设计特异引物，以胸肌cDNA进行Touchdown PCR检测。结果如图3所示，VGLL2和VGLL2-AS都可以扩增出目的条带，表明VGLL2-AS与VGLL2均可以进行双向转录。

图3 链特异引物鉴定鸡VGLL2-AS与VGLL2双向转录Fig.3 Identification of bidirectional transcription of chicken VGLL2-AS and VGLL2 by strand specific primers

2.4 VGLL2-AS与VGLL2可以形成双链RNA

为了探究VGLL2-AS的作用机制，利用核酸酶保护试验在VGLL2-AS与VGLL2的重叠区(第二外显子)设计了引物F2/R2(图4A)，之后扩增重叠区的目的片段以检测VGLL2-AS与VGLL2之间是否形成了RNA二聚体。如图4B所示，可以看到在F2与R2之间可以扩增到大小132 bp的目的片段，证明VGLL2-AS与VGLL2间形成了RNA二聚体。

A. VGLL2和VGLL2-AS之间重叠和非重叠区域的引物示意图；B. 消化产物扩增的琼脂糖凝胶电泳A. Schematic illustration of primers in overlapping and non-overlapping region between VGLL2 and VGLL2-AS; B. Agarose gel electrophoresis of amplification of digestion products图4 核酸酶保护试验Fig.4 Ribonuclease protection assay

2.5 VGLL2-AS与VGLL2的半衰期检测

既然VGLL2-AS与VGLL2可以形成 RNA双链，而两者之间形成的双链是否可以增加其在鸡原代成肌细胞增殖过程中的稳定性还不可知。为此，本研究进行了放线菌素D处理成肌细胞试验。结果如图5所示，与VGLL2的2 h左右的半衰期相比较，VGLL2-AS的半衰期可以达到5 h左右，表明VGLL2-AS半衰期较VGLL2长，证明其稳定性高于VGLL2基因。

A. VGLL2半衰期检测; B. VGLL2-AS半衰期检测A. The stability assay of VGLL2; B. The stability assay of VGLL2-AS图5 VGLL2和VGLL2-AS的半衰期检测Fig.5 The stability assay of VGLL2 and VGLL2-AS

2.6 VGLL2及其反义链转录本VGLL2-AS亚细胞定位

为了更进一步探究VGLL2-AS的作用机制，对VGLL2-AS和VGLL2进行亚细胞定位。分离的RNA产物用GAPDH(细胞质内参基因)和U6(细胞核内参基因)进行RT-qPCR检测，确定核质分离是否成功。如图6所示，VGLL2-AS和VGLL2基因都主要表达于细胞质中，说明在成肌细胞中，VGLL2-AS和VGLL2主要定位于鸡成肌细胞的细胞质中。

GAPDH和U6分别为细胞质和细胞核内参基因。经单因素方差分析检验，***表示差异极显著(P<0.001)GAPDH and U6 are cytoplasmic and nuclear reference genes, respectively. Tested by one-way analysis of variance, *** means extremely significant difference (P<0.001)图6 鸡原代成肌细胞中VGLL2-AS与VGLL2亚细胞定位Fig.6 Subcellular localization of hepatic VGLL2-AS and VGLL2

2.7 VGLL2-AS与VGLL2各转录本结构分析

利用NCBI在线网站查询VGLL2基因发现，VGLL2具有4个转录本，分别为：VGLL2-mRNA、VGLL2-X1、VGLL2-X2和VGLL2-X3。为了阐明VGLL2-AS与VGLL2各转录本的异同，首先利用PCR扩增及单克隆测序对VGLL2的转录本进行验证。如图7所示，利用VGLL2各转录本的共有引物VBF和VBR进行扩增后发现胸肌组织中VGLL2实际不存在VGLL2-X1转录本，只存在VGLL2-mRNA、VGLL2-X2和VGLL2-X3转录本。

图7 VGLL2不同转录本的凝胶电泳Fig.7 The agarose gel electrophoresis of VGLL2 transcripts

通过NCBI和USCS(http://genome.ucsc.edu/)数据库以及验证得到的各转录本序列信息，对VGLL2各转录本序列进行分析。结果如图8所示，VGLL2-X2和VGLL2-X3均含有5个外显子，而VGLL2-mRNA和VGLL2-AS分别含有4个和3个外显子；观察各个转录本之间的序列信息发现，与VGLL2的3种形式转录本相比，VGLL2-AS缺少外显子E3；此外，相比于VGLL2-X2和VGLL2-X3转录本，VGLL2-AS缺少VGLL2-X2的外显子I3a和VGLL2-X3的外显子I3b。

双斜线代表起始和终止密码子The double diagonal represents the initial and terminal codons图8 VGLL2-AS与VGLL2各剪切体结构分析Fig.8 Schematic comparison of transcripts of the VGLL2-AS and VGLL2 genes

2.8 VGLL2-AS与VGLL2在胸肌中的时序表达特征

前面的试验证明VGLL2-AS与VGLL2都在肌肉组织中高表达，但两者之间的相关性尚不清楚。为此，本研究利用RT-qPCR检测了它们在固始鸡胸肌组织中的时空表达规律(由于VGLL2-X2的表达量在各个组织中均较低且在其序列中没有找到合适的引物进行检测，因此本试验中只对VGLL2-mRNA和VGLL2-X3进行了检测)。结果如图9所示，VGLL2-mRNA和VGLL2-X3与VGLL2-AS表达趋势一致，且VGLL2-mRNA、VGLL2-X3和VGLL2-AS的表达呈现极强的正相关(r分别为0.943和0.935)，且差异极显著(P<0.01)。

图9 VGLL2-AS和VGLL2-mRNA (A)、VGLL2-X3 (B)在胸肌中的表达Fig.9 Expression patterns of VGLL2-AS, VGLL2-mRNA (A) and VGLL2-X3 (B) in breast muscle

3 讨 论

NATs可以直接与其正义链转录本相互作用，通过转录和表观遗传学水平来调节正义链基因表达，从而参与到各种生物学过程，对生物体的生长发育具有重要意义[13-16]。目前，在哺乳动物中关于NATs与其正义链转录本之间的研究较多，如DHRS4-AS[17]、p15AS[18]，而在家禽中的相关研究较少，如bFGF[19]、IGF-II[20]、Collagen[21]和Pdcd2[22]。因此，对NATs在家禽肌肉发育过程中的调控作用进行解析，既能补充lncRNAs在动物肌肉中的调控作用理论，又能为鸡肌肉发育和开发家禽的分子标记提供理论基础。

研究报道，半衰期可以衡量NATs的稳定性[23-24]。此外，与正义链转录本相比，大部分NATs的半衰期更长，这也意味着大多数NATs的稳定性更高，表明NATs有可能比其正义链转录本承担着更多的生物学功能[25-26]，而VGLL2-AS半衰期较VGLL2长，证明其稳定性高于VGLL2基因。同时，lncRNAs的功能与所处的细胞定位相关，而不同细胞定位的lncRNAs发挥的调节机制也不一样。根据细胞定位，lncRNAs又可分类为细胞质lncRNAs以及细胞核lncRNAs[27-30]。研究发现，大多数存在于细胞质中的lncRNAs可以通过与双链RNA之间的相互作用，进而参与到对mRNA稳定性和表达的调节。例如，主要存在于细胞质中的lncRNA MD1可以通过ceRNA的调控机制从而参与骨骼肌的分化进程[31]。而位于细胞核内的lncRNAs可能通过与靶基因的结合从而激活或者抑制靶基因表达，除此之外，还可能通过参与表观调控作用介导基因表达的调节[11, 32-34]。此外，已有研究表明，位于细胞核内的lncRNAs比细胞质中的lncRNAs的稳定性弱[35]。VGLL2-AS主要表达于鸡原代成肌细胞的细胞质中，暗示VGLL2-AS可能通过转录水平以及转录后水平来调控VGLL2。

基因的可变剪切能够改变转录本的结构和其编码的蛋白，同时也可以使基因产生多种结构不同的mRNA以及多肽序列，是基因表达多样性的一种重要的表现形式[36-39]。反义lncRNAs可以通过与其正义链转录本形成双链RNA以覆盖其顺式元件，从而直接调节前体mRNA的可变剪切。例如，NATs Zeb2补充了其正义链转录本的剪接位点，并导致其内含子的滞留，最终促进Zeb2的表达以及翻译[40]。内含子滞留是基因形成不同剪切体的方式之一，而分析VGLL2不同的剪切体信息，推测可能是VGLL2-AS与VGLL2通过碱基互补配对形成双链RNA，造成了VGLL2的内含子滞留，使之产生了3种剪切体，但VGLL2-AS是否影响VGLL2基因的可变剪切，仍待进一步的研究。

转录起始点之间的距离在1 kb碱基范围内的相邻的两个转录方向相反的基因被称为双向转录基因对，原核生物和真核生物的基因组中都存在双向转录基因，它们的表达大多是正相关的，并且在进化上保守[41-42]。VGLL2-AS和VGLL2-mRNA、VGLL2-X3在不同时期的胸肌组织中的表达具有很强的正相关，这与前期证明VGLL2-AS与VGLL2可以双向转录的结果一致。而VGLL2-AS的表达水平随着胸肌的发育而下调，这个基因是否负向调控肌肉发育，或者只在肌肉发育前期起作用有待进一步解析。

4 结 论

综上所述，VGLL2-AS作为VGLL2反义链编码的lncRNA定位于细胞质中，可能通过与VGLL2形成的双链RNA之间的相互作用，然后参与调节VGLL2的表达并维持其稳定性，最终在鸡的早期肌肉发育中发挥重要作用。本研究的结果扩展了鸡中关于NATs的研究，并为鸡VGLL2基因与其天然反义链转录本VGLL2-AS在鸡骨骼肌发育中的生物学功能奠定了基础，对于提高禽类的生长性能具有一定的意义。