Regulation effects of microRNA-1285 and its target DDX3X on Senecavirus A infected PK-15 cells
-
摘要:目的
探究MicroRNA-1285(miR-1285)及其靶标DDX3X在猪塞内卡病毒(Senecavirus A，SVA)感染PK-15细胞中的调控作用。
方法利用qRT-PCR、双荧光素酶活性及Western blot等方法研究miR-1285和DDX3X对I型干扰素(IFN-β)分泌及RIG-I信号通路的作用，分析miR-1285及DDX3X对SVA 3C蛋白基因表达的影响。
结果SVA感染PK-15细胞后，miR-1285表达量显著升高，并且miR-1285与DDX3X存在负靶向关系，二者可促进IFN-β转录及蛋白水平的表达。miR-1285通过靶向DDX3X对RIG-I信号通中的MAVS、TRAF3信号分子起调控作用。对于SVA 3C蛋白基因，DDX3X可以显著抑制其转录，并且可以逆转miR-1285所诱导的上调趋势。
结论SVA感染PK-15细胞后，宿主miR-1285及其靶标DDX3X对IFN-β及病毒3C蛋白的表达具有调控作用，研究结果将为明确miRNAs调控SVA感染的分子机制奠定基础，并为SVA的防控和诊断提供新的科学依据。
-
关键词:
- MicroRNA-1285 /
- DDX3X /
- 猪塞内卡病毒 /
- IFN-β /
- RIG-I信号通路
Abstract:ObjectiveTo explore the regulation roles of microRNA-1285 (miR-1285) and its target DDX3X in Senecavirus A (SVA) infected PK-15 cells.
MethodBy qRT-PCR, double luciferase activity and Western blot, the effects of miR-1285 and its target DDX3X on IFN-β secretion and the RIG-I signaling pathway were studied, and their effects on the expression of SVA 3C protein gene were analyzed.
ResultIn SVA infected PK-15 cells, the expression of miR-1285 increased significantly, and there was a negative targeting relationship between miR-1285 and DDX3X. Both miR-1285 and DDX3X promoted the transcription and protein expression of IFN-β. MiR-1285 regulated MAVS and TRAF3 signaling molecules in the RIG-I signaling pathway by targeting DDX3X. For SVA 3C protein, DDX3X significantly inhibited the transcription of 3C and reversed the up-regulation trend induced by miR-1285.
ConclusionAfter infecting PK-15 cells with SVA, host miR-1285 and its target DDX3X can regulate the expression of IFN-β and the viral 3C protein, which will lay a foundation for clarifying the molecular mechanism of miRNAs regulating SVA infection, and provide a new scientific basis for the prevention, control and diagnosis of SVA.
-
Keywords:
- MiRNA-1285 /
- DDX3X /
- Senecavirus A /
- IFN-β /
- RIG-I signal pathway
-
柑橘Citrusspp.是一种重要的经济作物，主要分布在我国的广东、广西和四川等地[1-3]。我国的柑橘种植面积和年产量已经位居世界首位，并且连续多年不断增长[4-5]。目前我国的柑橘灌溉由传统的漫灌逐渐升级到滴灌，能够较好地实现节水灌溉[6]。但是滴灌作业主要还是由人为主观意识根据气象状况来进行判断，缺乏客观依据[7-8]。柑橘生长需要适宜的土壤含水量和电导率，过低或者过高的土壤含水量都会对根系生理活动产生负面影响[9]，而土壤电导率则直接反映土壤的含盐量[10]，过高的盐胁迫会导致柑橘产生盐毒害，过低的盐含量导致柑橘难以吸收所需的无机盐离子。因此，在柑橘种植作业中，及时获取柑橘果园的气象、土壤数据，以及建立果园的土壤含水量和电导率预测模型，对于制定和指导果园灌溉施肥管理具有重要意义。
近年来，物联网技术在农业中得到了广泛的应用，其结合相关传感器能够实现果园生长环境的远程实时监测。杨伟志等[11]建立了基于物联网和人工智能的柑橘灌溉系统，该研究利用传感器实时监测柑橘果园土壤水分等信息，整合物联网系统与专家系统知识库的灌溉规则来制定灌溉控制措施，该系统主要依靠已有的控制规则，部分参数的设置需要依靠经验来确定。李中良等[12]利用物联网和空气温湿度传感器、土壤温湿度传感器设计了实时监测柑橘环境的系统，综合了无线传感器网络、通用无线分组服务(General packet radio service，GPRS)网络与远程服务器网络来实现田间数据的远程传输，但该系统没有针对果园环境数据进行进一步的分析。罗党等[13]利用灰色神经网络，研究降雨量、温度、湿度、风速等因子对土壤含水量的影响，该模型说明了神经网络在土壤含水量预测方面是可行的，具有不错的精度。蔡亮红等[14]设计了基于极限学习机(Extremelearning machine，ELM)算法的土壤含水量预测模型，采用定点采集土壤样本并获取光谱数据的方法，经过小波变换和特征工程处理后，利用ELM算法建立非线性模型，实现土壤含水量的预测，但该系统没有利用物联网等实时传输方法，而且加入的训练因子较少，在复杂环境条件下，还需要进一步研究。张武等[15]利用反向传播(Back propagation，BP)神经网络和长短期记忆(Long short-term memory，LSTM)深度学习模型研究了不同采样间隔对茶园土壤墒情预测模型性能的影响，结果表明在30 min的采样间隔下，模型的性能表现最好，但是该研究的数据采集时间跨度较短，模型也未给出更多可解释性指标来描述其模型性能。
针对以上研究现状，本文设计一种基于物联网的柑橘园环境信息实时采集与传输系统，并建立基于LSTM的土壤含水量和电导率预测模型，为果园灌溉施肥管理提供指导。
1. 柑橘园环境信息物联网传输系统设计
1.1 系统总体架构
物联网传输系统包括远程服务器、数据库、无线通信模块、信息采集节点、气象站等[16]。物联网系统架构如图1所示。信息采集节点负责采集柑橘园中的土壤含水量和电导率信息，气象站负责采集果园中的气象信息，包括空气温度、湿度、风速、降雨量等，通过无线通信模块传输至物联网网关与远程服务器，传输的数据最终存储在数据库中。
1.2 物联网传输系统硬件
硬件系统包括STM32F103C8T6处理器模块、电源模块、GPRS模块、ZigBee模块、空气温湿度传感器、风速传感器、降雨量传感器和土壤温度、含水量、电导率三合一传感器。传感器模块用于直接采集相关果园环境信息，ZigBee模块负责各个节点与网关数据传输，GPRS模块将各个节点数据远程传输至远程服务器。硬件框图如图2所示。
1.3 物联网传输系统软件
1.3.1 ZigBee组网协议
常见的ZigBee组网协议主要包括树型组网、星型组网和网状型组网等[17-19]。本文信息采集节点和气象站布置于柑橘果园中，需要考虑节点数据传输性能和续航能力。由于星型组网方式下，各个节点需要与网关直连以传输数据，在果园中存在着数据无法传输的可能性。网状型组网方式需要经过多轮数据转发，容易引起额外的能量损耗和丢包。因此本文采用树型组网方式，能够利用少量的节点中继实现数据传输，也能够尽可能减少路由转发引起的丢包问题。
1.3.2 节点软件
系统节点主要负责数据采集与转换，包含2种工作模式，即正常模式和休眠模式[20-21]。在正常工作模式下，节点负责开启传感器并采集数据；而在休眠模式下，节点仅保留时钟唤醒机制，定时开启并采集数据，以实现节能的目的。节点软件工作流程如图3所示。
2. 材料与方法
2.1 区域概况与数据采集预处理
试验地点设在华南农业大学试验基地，该基地面积约13 hm2，海拔约为200 m，属于山地地形，土壤主要为含沙土质。该地区夏季高温多雨，年降雨量为1700 mm左右，但是冬季降雨较少，易干旱。昼夜温差较大，适合柑橘等作物的糖分积累。在果园设置了5个节点，气象站和网关位于果园中央位置；使用激光测距仪测量得到单节点之间最长距离为344 m，节点使用高度为2.5 m的支架进行固定以减少果树对无线信号传播的干扰。其布置示意图如图4所示。节点数据与气象站数据分别传输至网关，由网关负责将数据打包传输至远程服务器，完成数据的传输。
本文的系统搭建完成后正式数据采集时间为2019年10月1日至2020年6月16日，节点数据传输间隔设置为每10 min传输1个数据包，单个节点理论上每天的数据包量为144个。数据采集期间平均每个节点接收数据包总数为36326个，总平均丢包率为2.6%。利用汇总后的数据包统计当天的最高空气温度、最低空气温度、平均空气温度、空气湿度、降雨量、土壤含水量、土壤电导率和土壤温度。土壤含水量和电导率一般与蒸散量、降雨量、温度、湿度等相关，根据气象参数和彭曼公式计算蒸散量[22]：
$$ {\rm{E}}{{\rm{T}}}_{0}=\dfrac{0.408{\Delta }\left({R}_{n}-G\right)+{\gamma }\dfrac{900}{T+273}{u}_{2}\left({e}_{{\rm{s}}}-{e}_{{\rm{a}}}\right)}{{\Delta }+{\gamma }\left(1+0.34{u}_{2}\right)} ,$$ (1) 式中，ET0为日参考蒸散量，mm·d−1；Rn为净辐射量，MJ·m−2·d−1；G为土壤热通密度，MJ·m−2·d−1；T为离地2 m高的平均空气温度，℃；u2为离地2 m高的风速，m·s−1；es为饱和蒸汽压，kPa；ea为实际蒸汽压，kPa；Δ为蒸汽压力曲线斜率，kPa·℃−1；γ为温湿度常数，kPa·℃−1。采用估算公式(2)计算辐射数据[23]：
$$ {R}_{{\rm{s}}}=\left({a}_{{\rm{s}}}+{b}_{{\rm{s}}}\frac{t}{{t}_{{\rm{max}}}}\right){R}_{{\rm{a}}} ,$$ (2) 式中，Rs和Ra分别为太阳辐射和地外辐射，MJ·m−2·d−1；t为日照时长，h；tmax为最大日照时长，h；
${a}_{{\rm{s}}}$ 和${b}_{{\rm{s}}}$ 通常取常数0.25和0.50。由于各个环境参数具有不同的数量级，因此需要在数据预处理环节进行归一化处理。各个节点预处理后的数据取平均值后，分别建立各个环境参数与土壤含水量和电导率的相关性分析[24]，结果如表1所示。
表 1 土壤含水量和电导率与环境参数的相关系数Table 1. Correlation coefficients of soil moisture, electrical conductivity and environmental parameters土壤因子
Soil factor空气温度 Air temperature 降雨量
Precipitation空气湿度
Air
humidity土壤温度
Soil
temperature蒸散量
Evapotranspiration最高值
Max.
最低值
Min.
平均值
Mean
含水量 Moisture content 0.490 0.598 0.644 0.123 0.430 0.652 0.158 电导率 Electrical conductivity 0.467 0.311 0.231 0.300 0.168 0.251 0.148 由表1可见，气象环境因子与土壤含水量和电导率存在一定的相关性，具备输入模型参与拟合的条件。
2.2 GRNN模型
GRNN模型最早是由Specht[25]提出的。GRNN模型能够利用连续的非线性数据进行训练，并收敛于底层回归曲面，该算法能够应用在学习植物模型的动力学中。GRNN建立在非参数回归基础上，以样本数据为后验条件，执行Parzen非参数估计，依据最大概率原则计算网络的输出并传输给下一级。GRNN的基础是径向基网络，具有较好的非线性拟合功能，但在训练过程中需要选定合适的学习参数，否则会容易产生过拟合现象。GRNN的网络结构分为4层，分别是输入层、模式层、求和层和输出层。其基础结构图如图5所示。
GRNN的网络输入为矩阵X=
$[ $ x1，x2，…，xn$]^\top$ ，输出为矩阵Y=$[ $ y1，y2，…，yk$]^\top$ ，k为输出向量的维数。GRNN的模式层的神经元数目等于输入样本数目n，模式层的传输函数(pi)为公式(3)：$$ {p}_{i}=\mathrm{exp}\left[-\frac{{\left({\boldsymbol{X}}-{{\boldsymbol{X}}}_{i}\right)}^\top \left({\boldsymbol{X}}-{{\boldsymbol{X}}}_{i}\right)}{2{\mathrm{\delta }}^{2}}\right],{\text{其中}},i=\mathrm{1,2},\cdots ,n ,$$ (3) 式中，i表示模式层中第i个神经元，Xi为第i个神经元对应的学习样本，
$ \delta $ 为平滑参数，改变$ \delta $ 即可改变网络学习率。GRNN网络求和层包含2种求和方式，即直接求和与加权求和，计算公式分别如公式(4)和(5)所示：
$${S_{\rm{D}}} = \sum\limits_{i = 1}^n {{p_i}} ,$$ (4) $${S_{{\rm{N}}j}} = \sum\limits_{i = 1}^n {{y_{ij}}} {P_i},{\text{其中}},j = 1,2, \cdots ,k,$$ (5) 式中，
${S}_{{\rm{D}}}$ 为求和层神经元输出，${S}_{{\rm{N}}j}$ 为求和层加权神经元输出，j表示输出层中第j个神经元，$ {P}_{i} $ 为模式层第i个神经元的权值。GRNN网络的输出层中神经元数目等于输出向量的维数，如公式(6)所示：
$$ {y}_{j}=\frac{{S}_{{\rm{N}}j}}{{S}_{{\rm{D}}}},{\text{其中}},j={1,2},\cdots ,k{\text{。}} $$ (6) 2.3 LSTM模型
实践证明，神经网络能够有效地解决复杂的非线性问题，但传统的神经网络无法利用时间序列数据特征预测未来的情况，特别是对市场销售、股票价格、pm2.5等时间序列数据的时序建模。时间序列数据不仅与当前条件相关，而且与最后的数据具有很高的相关性，而传统的神经网络难以对时间序列进行建模。近年来，递归神经网络[26-28](Recursive neural network，RNN)在自然语言处理、机器翻译和文本预测等方面得到了成功的应用，其效果明显优于神经网络。不同于一般的前馈网络，RNN包含隐藏状态提取特征，然后转化为输出数据。另一方面，RNN在所有计算步骤中共享权值、偏差等，使其与所连接的节点相依赖。然而，标准RNN也有不足之处，如要求输入和输出数据大小相同。另外，RNN通过对前一层隐状态的变换,经过激活函数后得到新的隐状态，但当RNN进行反向传播时容易产生梯度消失或者梯度爆炸问题。因此，RNN有时不能很好地处理长期依赖问题。
作为一种特殊的RNN网络，LSTM[29]包含了1个控制应该丢弃哪些信息的门，即遗忘门(Forget gate)。其次，为了避免梯度消失问题，在计算隐层状态时，用求和代替梯度导数。因此，理论上LSTM具有解决长期问题的能力。在本文中，LSTM网络是由3个部分组成的，即1个输入层用于处理数据，50个隐藏层用于训练和建模，1个输出层用于输出建模结果。模型的核心部分即LSTM节点[30]，主要由3个部分组成，分别是输入门、输出门和遗忘门。其结构如图6所示。
其中，遗忘门控制丢弃的信息(ft)，如公式(7)所示：
$${{{f}}_{{t}}} = \sigma \left( {{{{W}}_{\rm{f}}} \left[ {{{{h}}_{{{t}} - 1}},{{{x}}_{{t}}}} \right] + {{{b}}_{\rm{f}}}} \right),$$ (7) 式中，
$ {f}_{t} $ 的范围为[0,1]，$ \mathrm{\sigma } $ 函数使得方程非线性化，更接近实际值从而提高预测精度。${{W}}_{\mathrm{f}}$ 和${{b}}_{\mathrm{f}}$ 分别为遗忘门的权重和偏差参数。${{h}}_{{t}-1}$ 和${{x}}_{{t}}$ 为上一状态值和当前输入值。输入门的功能是控制输入到当前单元中的信息，其状态更新如公式(8)和(9)所示：
$${{{i}}_{{t}}} = \sigma \left( {{{{W}}_{\rm{i}}} \left[ {{{{h}}_{{{t}} - 1}},{{{x}}_{{t}}}} \right] + {{{b}}_{\rm{i}}}} \right),$$ (8) $$\widetilde {{{{C}}_{{t}}}} = {\rm{tanh}}\left( {{{{W}}_{\rm{C}}} \left[ {{{{h}}_{{{t}} - 1}},{{{x}}_{{t}}}} \right] + {{{b}}_{\rm{C}}}} \right),$$ (9) 式中，
${{i}}_{{t}}{\text{由}}{\sigma }$ 激活函数控制输入数据，如果其值为0，表示输入数据被完全阻塞，如果值为1，则表示完全允许数据进入门中。细胞状态由$\widetilde {{{C}}_{{t}}}$ 负责，它的激活函数是tanh。Wi和bi表示激活函数的权重和偏差参数，WC和bC表示细胞状态方程的权重和偏差参数。输入门经过运算的数据，需要更新到细胞状态中，更新公式为：
$${{{C}}_{{t}}} = {{{f}}_{{t}}}{{{C}}_{{{t}} - 1}} + {{{i}}_{{t}}}\widetilde {{{{C}}_{{t}}}},$$ (10) 式中，
${f}_{t} {C}_{t-1}$ 代表遗忘的信息，${i}_{t} \widetilde {{C}_{t}}$ 代表需要更新到当前细胞状态的信息，两者之和共同决定向下一个状态传递的信息量。LSTM网络结构的最后一部分是输出门，其更新公式为：
$${{{o}}_{{t}}} = \sigma \left( {{{{W}}_{\rm{o}}}\left[ {{{{h}}_{{{t}} - 1}},{{{x}}_{{t}}}} \right] + {{{b}}_{\rm{o}}}} \right),$$ (11) $${{{h}}_{{t}}} = {{{o}}_{{t}}}{\rm{tanh}}\left( {{{{C}}_{{t}}}} \right),$$ (12) 式中，
$o_t $ 表示激活函数的输出，由σ函数决定，同时更新后的细胞状态值$C_t $ 在tanh函数作用下，与传递函数$o_t $ 的乘积作为向下一个单元输出的内容，Wo和bo表示激活函数的权重和偏差参数。经过上述LSTM单元的运算后，输出的信息能够被存储，也能够进入下一个节点进行运算，从而在训练过程中提取出信息。
2.4 模型性能测试标准
本文模型的性能表现使用均方根误差(RMSE)来表示，计算模型预测值与真实测量值的决定系数(R2)[31]来表征模型可解释性。相关公式为：
$${\rm{RMSE}} = \sqrt {\frac{{\displaystyle\sum _{m = 1}^{{N}}{{\left( {{{{y}}_{{\rm{pred}}}} - {{{y}}_{{\rm{cal}}}}} \right)}^2}}}{{{N}}}}, $$ (13) $${R^2} = \frac{{\displaystyle\sum _{m = 1}^N{{\left( {{{{y}}_{{\rm{pred}}}} - \overline {{{{y}}_{{\rm{cal}}}}} } \right)}^2}}}{{\displaystyle\sum _{m = 1}^N{{\left( {{{{y}}_{{\rm{cal}}}} - \overline {{{{y}}_{{\rm{cal}}}}} } \right)}^2}}},$$ (14) 式中，
$ {{y}}_{\rm{pred}} $ 表示预测值，$ {{y}}_{\rm{cal} }$ 表示实测值，$ \overline{{{y}}_{\rm{pred}}} $ 表示预测值均值，$ \overline{{{y}}_{\rm{cal}} }$ 表示实测值均值，N表示样本数。3. 结果与分析
3.1 LSTM和GRNN模型训练设置
试验中，使用了果园气象站气象数据和5个信息采集节点采集到的土壤墒情数据，现场布置的节点如图7所示。由于无线传感器网络传输过程中存在复杂因素的干扰，产生了个别异常(如空数据)与丢包数据，使得数据包不完整或整体丢失，本文将异常数据进行剔除，最终得到符合模型需求的数据。将70%数据划分为训练集，30%数据划分为测试集，5个节点数据划分结果如表2所示。训练参数设置方面，LSTM模型包含2层，每层128个LSTM单元，使用Relu激活函数，损失函数使用均方误差(MSE)，优化算法使用Adam，共设置72个epoch和128个batch size。GRNN模型包含4层，模式层与输入样本相同，求和层使用128个神经元，激活函数、损失函数和优化算法、epoch和batch size等参数与LSTM模型保持一致。
表 2 节点数据信息Table 2. Datasets of nodes节点
Node数据总数
Total number训练集
Training dataset测试集
Testing dataset1 216 151 65 2 241 168 73 3 241 168 73 4 196 137 59 5 203 142 61 气象站
Weather station241 168 73 3.2 模型性能表现
基于本文所提出的模型和预处理的输入样本经过训练后的性能测试标准如表3所示。
表 3 模型性能表现Table 3. Performance of models节点
Node模型
Model土壤含水量
Soil moisture
content土壤电导率
Soil electrical
conductivityRMSE R2 RMSE R2 1 LSTM 8.65 0.811 6.88 0.850 GRNN 13.30 0.126 12.30 0.106 2 LSTM 7.30 0.906 6.68 0.727 GRNN 14.70 0.323 13.70 0.185 3 LSTM 6.74 0.814 8.22 0.648 GRNN 7.01 0.369 8.26 0.268 4 LSTM 7.93 0.760 8.50 0.767 GRNN 10.27 0.256 7.85 0.167 5 LSTM 8.40 0.768 7.78 0.769 GRNN 8.80 0.206 7.60 0.132 从表3中可以看出，在土壤含水量的预测模型表现中，LSTM模型的RMSE值范围为6.74~8.65，低于GRNN模型的7.01~14.70。LSTM模型中，节点2和3的RMSE在5个节点中较低，这可能是因为这2组训练样本数较多。节点4的RMSE比节点5更低，但其训练样本数比后者少，因此RMSE除了受训练样本数量影响外，还与数据集自身样本情况有关。GRNN模型的RMSE与训练样本数量没有呈现明显的关系或趋势。RMSE用来衡量离散程度，与原样本具有相同的量纲，而R2作为一种无量纲的性能指标，能够反映模型的可解释性。从表3可以看出，对于土壤含水量，虽然相较于GRNN模型，LSTM模型的RMSE没有表现出明显优势，但是其R2范围为0.760~0.906，远大于GRNN模型的0.126~0.369，这意味着LSTM模型预测的土壤含水量呈现出了较高的可解释性，意味着该模型给出的预测值能够较好地反映土壤含水量。表3中土壤电导率模型预测表现与土壤含水量呈现类似的趋势，同样表明LSTM模型在给定的试验样本下，适合用来对土壤电导率进行预测。
3.3 模型拟合分析
LSTM模型的预测精度较高、性能稳定，故而使用该模型观察预测效果。利用模型和气象数据得到预测值，与节点采集到的实测值进行对比，结果如图8和图9的a1~e1所示，之后利用最小二乘法进行回归拟合，其结果如图8和图9的a2~e2所示。其中，5个节点分别布置在果园5个分区中，由于果园面积较大，同时不同区域有高低落差，因此5个节点数据分别代表果园各个分区的土壤墒情。
![图 8 LSTM模型在不同节点的土壤含水量实测值与预测值的比较]()
图 8 LSTM模型在不同节点的土壤含水量实测值与预测值的比较a1~e1中，蓝线代表实测值，红线代表预测值Figure 8. Comparison of measured and predicted soil moisture of LSTM model in different nodesIn a1−e1, blue line represents measured value, and red line represents predicted value![图 9 LSTM模型在不同节点的土壤电导率实测值与预测值的比较]()
图 9 LSTM模型在不同节点的土壤电导率实测值与预测值的比较a1~e1中，蓝线代表实测值，红线代表预测值Figure 9. Comparison of measured and predicted soil electrical conductivity of LSTM model in different nodesIn a1−e1, blue line represents measured value, and red line represents predicted value节点1的土壤含水量(图8a1、a2)预测值相对于实测值偏小，而拟合得到的回归方程中，实测的土壤含水量数据点主要集中在56~76 mm区间，这个趋势在图8a1中的折线图中无法反映出来，这可能是由于节点1处在相对低洼的位置，加上外界环境变化如大量降雨积水或者人工灌溉导致的短时间内土壤含水量升高，在图9a1中可以观察到此时的土壤电导率在土壤含水量升高后迅速下降的现象。节点2土壤含水量(图8b1、8b2)的回归拟合决定系数高达0.906，折线图中，预测值很好地贴近了实测的土壤含水量。节点3(图8c1、8c2)的拟合效果与节点2接近，说明模型在节点2和3的数据集中能够获得较好的训练结果。节点4和5(图8d2、8e2)模型拟合的R2分别为0.760和0.768，低于0.8，根据回归方程图可以看到，主要是由于个别数据相对于其他样本偏差较大导致的。这2个节点对应模型的RMSE并不是最大的，这也再次说明了仅根据RMSE指标不能证明模型的性能高低。
从图9可以看出，节点3的土壤电导率(图9c1、9c2)数据点多集中在30~55 μS.cm−1区间，从折线图中可以看到预测值与实测值曲线比较接近，而其回归拟合R2仅为0.648，说明模型预测性能会受到数据分布的影响，较多的数据量可能会减少此种影响。虽然使用了同样的气象数据，但土壤电导率的预测模型性能总体上稍逊于土壤含水量预测模型，这可能是由于土壤电导率除了受气象因子、柑橘果树蒸散量影响外，还受到人工施肥、生物无机盐代谢等多种因素的影响，这可能是需要进一步研究的方向。
4. 结论
本文提出了利用物联网来实现远程实时获取柑橘果园土壤含水量和土壤电导率的方法，并搭建气象站获取柑橘果园气象信息，利用LSTM和GRNN模型对土壤墒情进行建模，实现对柑橘园土壤含水量与土壤电导率的预测。部分已有的对土壤墒情进行建模的研究，以BP神经网络或光谱反演方法为基础，存在着实时性不足、长时间预测精度较低的问题。本文利用物联网及相关传感器搭建远程获取果园环境数据的系统，在果园中设置5个土壤墒情采集节点和1个气象站，实现了实时采集和传输果园环境数据与气象数据的目的。
利用采集到的气象数据与节点土壤墒情数据在LSTM模型和GRNN模型上的训练结果表明：LSTM模型训练的土壤含水量和土壤电导率的RMSE范围分别为6.74~8.65和6.68~8.50，R2范围分别为0.76~0.906和0.648~0.850，相比于GRNN模型的RMSE和R2(7.01~14.70和7.60~13.70，0.126~0.369和0.132~0.268)，表现出了较好的模型性能与预测效果。在LSTM模型的土壤含水量和土壤电导率拟合分析中可以发现，模型的RMSE与R2之间并没有明确的相关性，即较小的RMSE并不意味着较高的R2，因此只用RMSE来衡量模型性能可能是片面的。
LSTM模型的预测值与实测值的拟合回归方程和对应的折线图说明模型的预测效果较好，能够满足预测的需求。但是部分节点如节点3的土壤电导率拟合决定系数较低，说明数据样本分布可能会对模型效果产生影响。本文建立的LSTM预测模型在土壤电导率的预测性能方面稍逊于土壤含水量模型，这可能是由于土壤电导率受到环境中其他因素影响更多，需要更进一步地深入研究。另外，果园的土壤含水量和土壤电导率在不同的种植和生理条件下，存在一定的差异，根据生产过程中的灌溉与施肥情况，利用本文提出的方法能够为柑橘果园的灌溉和施肥提供一定的指导，也可以根据LSTM预测模型，结合未来一段时间的气象数据与当前田间土壤墒情数据，预测未来一段时间的土壤墒情数据，及时掌握土壤墒情变化趋势，方便果园管理者提前做好施肥灌溉准备工作，能够帮助实现科学精细的果园灌溉与施肥管理，带来较好的管理效益。
-
![]()
图 1 不同SVA感染时间(A)和感染剂量(B)条件下PK-15细胞中miR-1285的表达量
“*”“**”分别表示处理与对照在P < 0.05和P < 0.01水平差异显著(Duncan’s法)
Figure 1. Expression of miR-1285 in PK-15 cells infected by SVA at different infection time (A) and dosages (B)
“*” and “**” represented statistical difference in comparison with control group at P < 0.05 and P < 0.01 levels respectively (Duncan’s method)
![]()
图 3 转染miR-1285 mimics、inhibitor至PK-15细胞后DDX3X蛋白的表达
Figure 3. DDX3X protein expression after the transfection of miR-1285 mimics, inhibitor into PK-15 cells
![]()
图 4 转染miR-1285 mimics、inhibitor至PK-15细胞后DDX3X的mRNA相对表达量
“**”表示处理与对照在P < 0.01水平差异显著 (Duncan’s 法)
Figure 4. The relative expression of DDX3X mRNA after transfection of miR-1285 mimics and inhibitor into PK-15 cells
“**” represented statistical difference in comparison with control group at P < 0.01 level (Duncan’s method)
![]()
图 5 转染不同DDX3X重组载体质粒至PK-15细胞后miR-1285双荧光素酶活性
“*”表示处理与对照在P < 0.05水平差异显著 (Duncan’s 法)
Figure 5. The relative dual-luciferase activity of miR-1285 after transfection of different DDX3X recombinant vector plasmids into PK-15 cells
“*” represented statistical difference in comparison with control group at P < 0.05 level (Duncan’s method)
![]()
图 6 miR-1285及其靶标DDX3X对IFN-β mRNA表达的调控作用
“*”和“**”分别表示处理与对照在P < 0.05和P < 0.01水平差异显著(Duncan’s法)
Figure 6. Regulation effects of miR-1285 and its target DDX3X on the mRNA expression of IFN-β
“*” and “**” represented statistical difference in comparison with control group at P < 0.05 and P < 0.01 levels respectively (Duncan’s method)
![]()
图 7 miR-1285及其靶标DDX3X对IFN-β 蛋白表达的调控作用
Figure 7. Regulation effects of miR-1285 and its target DDX3X on the protein expression of IFN-β
![]()
图 8 miR-1285对RIG-I通路信号转导分子的影响
“*”表示处理与对照在P < 0.05水平差异显著(Duncan’s法)
Figure 8. Effects of miR-1285 on signal transduction molecules of the RIG-I pathway
“*” represented statistical difference in comparison with control group at P < 0.05 level respectively (Duncan’s method)
![]()
图 9 DDX3X沉默对RIG-I通路信号转导分子的影响
“*”和“**”分别表示处理与对照在P < 0.05和P < 0.01水平差异显著(Duncan’s法)
Figure 9. Effects of DDX3X silencing on signal transduction molecules of the RIG-I pathway
“*” and “**” represented statistical difference in comparison with control group at P < 0.05 and P < 0.01 levels respectively (Duncan’s method)
![]()
图 10 miR-1285及其靶标DDX3X对SVA 3C表达的影响
“**”表示处理与对照在P < 0.01水平差异显著(Duncan’s法)
Figure 10. Effects of miR-1285 and its target DDX3X on the expression of SVA 3C
“**” represented statistical difference in comparison with control group at P < 0.01 level respectively (Duncan’s method)
表 1 基因引物序列
Table 1 Primer sequences of genes
基因名称 Gene name 引物/探针序列(5′→3′) Primer/Probe sequence θ退火/ ℃ Annealing temperature 产物大小/bp Product size 文献 Reference RIG-I F: ATCCCAGCAACGAGAA 60 188 [36] R: GCCACGTCCAGTCAAT MDA5 F: GAGGAATCAGCACGAGGAA 58 73 [37] R: GTCAGTAATCCACTGGGA MAVS F: ATAGCCAGCCTTTCTCGG 60 237 [36] R: TAGCCTCAGTCTTGACCTCTTC TRAF3 F: GTGTCAAGAAGGCATCG 60 164 [36] R: CCTCAAACTGGCAATCA TANK F: GGACGCCTTGAACTACCTGT 60 119 R: GCCTGCCGAAAGGCTTCATA TBK1 F: GCCTTTCTCGGGGTCTTCAA 60 74 R: ACACTTTTCCTGATCCGCCT IRF3 F: CCAGTGGTGCCTACACTCCT 61 191 [38] R: AGAGGTGTCTGGCTCAGGAA IRF7 F: CGCCTCCTGGAAAACCAA 60 76 [37] R: CCCTGAGTTGTCCTGCAACA IFN-β F: GCTAACAAGTGCATCCTCCAAA 60 77 [39] R:AGCACATCATAGCTCATGGAAAGA GAPDH F: ACATGGCCTCCAAGGAGTAAGA 60 106 [40] R: GATCGAGTTGGGGCTGTGACT SVA-3C F: GAGCTTCAATCTCCTAGA 59 115 R: GTGTCATCATTCTCGTTAG 探针 Probe CAGACATTCGAGCCAAGCAACAA 69 
下载: 导出CSV
-
[1] YATES L A, NORBURY C J, GIBERT R J C. The long and short of microRNA[J]. Cell, 2013, 153(3): 516-519. doi: 10.1016/j.cell.2013.04.003
[2] CHENG J, WU R, LONG L, et al. MiRNA-451a targets IFN regulatory factor 8 for the progression of systemic lupus erythematosus[J]. Inflammation, 2017, 40(2): 676-687. doi: 10.1007/s10753-017-0514-8
[3] WANG L, ZHOU L, HU D, et al. Porcine reproductive and respiratory syndrome virus suppresses post-transcriptionally the protein expression of IFN-β by upregulating cellular microRNAs in porcine alveolar macrophages in vitro[J]. Experimental and Therapeutic Medicine, 2018, 15(1): 115-126.
[4] MORIN R D, O’CONNOR M D, GRIFFITH M, et al. Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells[J]. Genome Research, 2008, 18(4): 610-621. doi: 10.1101/gr.7179508
[5] ANTHON C, TAFER H, HAVGAARD J H, et al. Structured RNAs and synteny regions in the pig genome[J]. BMC Genomics, 2014, 15: 459. doi: 10.1186/1471-2164-15-459.
[6] GAO X, WANG Y, ZHAO H, et al. Plasma miR-324-3p and miR-1285 as diagnostic and prognostic biomarkers for early stage lung squamous cell carcinoma[J]. Oncotarget, 2016, 7(37): 59664-59675. doi: 10.18632/oncotarget.11198
[7] BORRELLI N, DENARO M, UGOLINI C, et al. MiRNA expression profiling of ‘noninvasive follicular thyroid neoplasms with papillary-like nuclear features’ compared with adenomas and infiltrative follicular variants of papillary thyroid carcinomas[J]. Modern Pathology, 2017, 30(1): 39-51. doi: 10.1038/modpathol.2016.157
[8] LIU J, YAN J, ZHOU C, et al. MiR-1285-3p acts as a potential tumor suppressor miRNA via down regulating JUN expression in hepatocellular carcinoma[J]. Tumor Biology, 2015, 36: 219-225. doi: 10.1007/s13277-014-2622-5
[9] HUANG H, XIONG G, SHEN P, et al. MicroRNA-1285 inhibits malignant biological behaviors of human pancreatic cancer cells by negative regulation of YAP1[J]. Neoplasma, 2017, 64(3): 358-366. doi: 10.4149/neo_2017_306
[10] LIAO J, LI Q, HU Z, et al. Mitochondrial miR-1285 regulates copper-induced mitochondrial dysfunction and mitophagy by impairing IDH2 in pig jejunal epithelial cells[J]. Journal of Hazardous Materials, 2022, 422: 126899. doi: 10.1016/j.jhazmat.2021.126899.
[11] VENKATARAMANNAN S, GADEK M, CALVIELLO L, et al. DDX3X and DDX3Y are redundant in protein synthesis[J]. RNA, 2021, 27(12): 1577-1588. doi: 10.1261/rna.078926.121
[12] VAN VOSS M R H, KAMMERS K, VESUNA F, et al. Global effects of DDX3 inhibition on cell cycle regulation identified by a combined phosphoproteomics and single cell tracking approach[J]. Translational Oncology, 2018, 11(3): 755-763. doi: 10.1016/j.tranon.2018.04.001
[13] HEATON S M, BORG N A, DIXIT V M. Ubiquitin in the activation and attenuation of innate antiviral immunity[J]. Journal of Experimental Medicine, 2016, 213(1): 1-13. doi: 10.1084/jem.20151531
[14] HATHAICHOTI S, VISITNONTHACHAI D, NGAMSIRI P, et al. Paraquat induces extrinsic pathway of apoptosis in A549 cells by induction of DR5 and repression of anti-apoptotic proteins, DDX3 and GSK3 expression[J]. Toxicology in Vitro, 2017, 42: 123-129. doi: 10.1016/j.tiv.2017.04.016
[15] LIN T C. DDX3X multifunctionally modulates tumor progression and serves as a prognostic indicator to predict cancer outcomes[J]. International Journal of Molecular Sciences, 2019, 21(1): 281. doi: 10.3390/ijms21010281.
[16] SOULAT D, BŰRCKSTŰMMER T, WESTERMAYER S, et al. The DEAD-box helicase DDX3X is a critical component of the TANK-binding kinase 1-dependent innate immune response[J]. EMBO Journal, 2008, 27(15): 2135-2146. doi: 10.1038/emboj.2008.126
[17] SCHRŐDER M, BARAN M, BOWIE A G. Viral targeting of DEAD box protein 3 reveals its role in TBK1/IKKepsilon-mediated IRF activation[J]. EMBO Journal, 2008, 27(15): 2147-2157. doi: 10.1038/emboj.2008.143
[18] OSHIUMI H, SAKAI K, MATSUMOT M, et al. DEAD/H BOX 3 (DDX3) helicase binds the RIG-I adaptor IPS-1 to up-regulate IFN-β-inducing potential[J]. European Journal of Immunology, 2010, 40(4): 940-948. doi: 10.1002/eji.200940203
[19] SZAPPANOS D, TSCHISMAROV R, PERLOT T, et al. The RNA helicase DDX3X is an essential mediator of innate antimicrobial immunity[J]. PLoS Pathogens, 2018, 14(11): e1007397. doi: 10.1371/journal.ppat.1007397
[20] ADAMS M J, LEFKOWITZ E J, KING A M Q, et al. Ratification vote on taxonomic proposals to the international committee on taxonomy of viruses[J]. Archives of Virology, 2016, 161(10): 2921-2949. doi: 10.1007/s00705-016-2977-6
[21] VANNUCCI F A, LINHARES D C L, BARCELLOS D E S N, et al. Identification and complete genome of Seneca Valley virus in vesicular fluid and sera of pigs affected with idiopathic vesicular disease, Brazil[J]. Transboundary Emerging Disease, 2015, 62: 589-593. doi: 10.1111/tbed.12410
[22] CANNING P, CANON A, BATES J L, et al. Neonatal mortality, vesicular lesions and lameness associated with Senecavirus A in a U. S. sow farm[J]. Transboundary and Emerging Disease, 2016, 63(4): 373-378. doi: 10.1111/tbed.12516
[23] ZHANG X L, ZHU Z X, YANG F, et al. Review of Seneca Valley virus: A call for increased surveillance and research[J]. Frontiers in Microbiology, 2018, 9: 940. doi: 10.3389/fmicb.2018.00940.
[24] FERNANDES M H V, MAGGIOLI M F, OTTA J, et al. Senecavirus A 3C protease mediates host cell apoptosis late in infection[J]. Frontiers in Immunology, 2019, 10: 363. doi: 10.3389/fimmu.2019.00363.
[25] LEME R A, OLIVEIRA T E S, ALFIERI A F, et al. Pathological, immunohistochemical and molecular findings associated with Senecavirus A-induced lesions in neonatal piglets[J]. Journal of Comparative Pathology, 2016, 155(2/3): 145-155.
[26] HAUSE B M, MYERS O, DUFF J, et al. Senecavirus A in pigs, United States, 2015[J]. Emerging Infectious Diseases, 2016, 22(7): 1323-1325. doi: 10.3201/eid2207.151591
[27] XU W, HOLE K, GOOLIA M, et al. Genome wide analysis of the evolution of Senecavirus A from swine clinical material and assembly yard environmental samples[J]. PLoS One, 2017, 12(5): e0176964. doi: 10.1371/journal.pone.0176964
[28] WU Q, ZHAO X, BAI Y, et al. The first identification and complete genome of Senecavirus A affecting pig with idiopathic vesicular disease in China[J]. Transboundary and Emerging Disease, 2017, 64(5): 1633-1640. doi: 10.1111/tbed.12557
[29] SAENG-CLUTO K, RODTIAN P, TEMEEYASEN G, et al. The first detection of Senecavirus A in pigs in Thailand, 2016[J]. Transboundary and Emerging Disease, 2018, 65(1): 285-288. doi: 10.1111/tbed.12654
[30] SUN D, VANNUCCI F, KNUTSON T P, et al. Emergence and whole-genome sequence of Senecavirus A in Colombia[J]. Transboundary and Emerging Diseases, 2017, 64(5): 1346-1349. doi: 10.1111/tbed.12669
[31] ARZT J, BERTRAM M R, VU L T, et al. First detection and genome sequence of Senecavirus A in Vietnam[J]. Microbiology Resource Announcements, 2019, 8(3). doi: 10.1128/MRA.01247-18.
[32] SUN Y, CHENG J, WU R T, et al. Phylogenetic and genome analysis of 17 novel Senecavirus A isolates in Guangdong Province, 2017[J]. Frontiers in Veterinary Science, 2018, 5: 314. doi: 10.3389/fvets.2018.00314.
[33] QIAN S, FAN W, QIAN P, et al. Isolation and full-genome sequencing of Seneca Valley virus in piglets from China, 2016[J]. Virology Journal, 2016, 13: 173. doi: 10.1186/s12985-016-0631-2.
[34] ZHU Z, YANG F, CHEN P, et al. Emergence of novel Seneca Valley virus strains in China, 2017[J]. Transboundary and Emerging Diseases, 2017, 64(4): 1024-1029. doi: 10.1111/tbed.12662
[35] WANG H, LI C, ZHAO B, et al. Complete genome sequence and phylogenetic analysis of Senecavirus A isolated in Northeast China in 2016[J]. Archives of Virology, 2017, 162(10): 3173-3176. doi: 10.1007/s00705-017-3480-4
[36] ZHANG J, MIAO J, HOU J, et al. The effects of H3N2 swine influenza virus infection on TLRs and RLRs signaling pathways in porcine alveolar macrophages[J]. Virology Journal, 2015, 12: 61. doi: 10.1186/s12985-015-0284-6.
[37] 谢立兰. 伪狂犬病毒和猪传染性胃肠炎病毒诱导β干扰素产生的分子机制研究[D]. 武汉: 华中农业大学, 2011. [38] ISLAM M A, GROßE-BRINKHAUS C, PRŐLL M J, et al. Deciphering transcriptome profiles of peripheral blood mononuclear cells in response to PRRSV vaccination in pigs[J]. BMC Genomics, 2016, 17: 641. doi: 10.1186/s12864-016-2849-1.
[39] 王荡. 口蹄疫病毒Lpro和3Cpro调控宿主抗病毒天然免疫反应的分子机制研究[D]. 武汉: 华中农业大学, 2011. [40] 王国庆. 口蹄疫病毒2B蛋白拮抗RIG-I抗病毒作用研究[D]. 兰州: 甘肃农业大学, 2015. [41] DERRICK T, ROBERTS C H, RAJASEKHAR M, et al. Conjunctival microRNA expression in inflammatory trachomatous scarring[J]. PLoS Neglected Tropical Diseases, 2013, 7(3): e2117. doi: 10.1371/journal.pntd.0002117
[42] KAUR S, KRISHN S R, RACHAGANI S, et al. Significance of microRNA-based biomarkers for pancreatic cancer[J]. Annals of Translational Medicine, 2015, 3(18): 277.
[43] LIU Y, XU X, XU X, et al. MicroRNA-193a-3p inhibits cell proliferation in prostate cancer by targeting cyclin D1[J]. Oncology Letters, 2017, 14(5): 5121-5128.
[44] MORGUL M H, KLUNK S, ANASTASIADOU Z, et al. Diagnosis of HCC for patients with cirrhosis using miRNA profiles of the tumor-surrounding tissue: A statistical model based on stepwise penalized logistic regression[J]. Experimental and Molecular Pathology, 2016, 101(2): 165-171. doi: 10.1016/j.yexmp.2016.07.014
[45] POIRIER J T, DOBROMILSKAVA I, MORIARTY W F, et al. Selective tropism of Seneca Valley virus for variant subtype small cell lung cancer[J]. Journal of the National Cancer Institute, 2013, 105(14): 1059-1065. doi: 10.1093/jnci/djt130
[46] ZHOU Z H, SUN Y, YAN X L, et al. Swine acute diarrhea syndrome coronavirus (SADS-CoV) antagonizes interferon-β production via blocking IPS-1 and RIG-I[J]. Virus Research, 2020, 278: 197843. doi: 10.1016/j.virusres.2019.197843.
[47] ANDREJEVA J, CHILDS K S, YOUNG D F, et al. The V proteins of paramyxoviruses bind the IFN-inducible RNA helicase, mda-5, and inhibit its activation of the IFN-beta promoter[J]. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(49): 17264-17269. doi: 10.1073/pnas.0407639101
[48] CAO L, GE X, GAO Y, et al. Porcine epidemic diarrhea virus inhibits dsRNA-induced interferon-beta production in porcine intestinal epithelial cells by blockade of the RIG-I-mediated pathway[J]. Virology Journal, 2015, 12: 127. doi: 10.1186/s12985-015-0345-x.
[49] CHEN Q, LIU Q, LIU D, et al. Molecular cloning, functional characterization and antiviral activity of porcine DDX3X[J]. Biochemical and Biophysical Research Communications, 2014, 443(4): 1169-1175. doi: 10.1016/j.bbrc.2013.12.098
[50] ZHAO X, WU Q, BAI Y, et al. Phylogenetic and genome analysis of seven Senecavirus A isolates in China[J]. Transboundary and Emerging Diseases, 2017, 64(6): 2075-2082. doi: 10.1111/tbed.12619
[51] THOMPSON S R, SARNOW P. Regulation of host cell translation by viruses and effects on cell function[J]. Current Opinion in Microbiology, 2000, 3(4): 366-370. doi: 10.1016/S1369-5274(00)00106-5
[52] QIAN S, FAN W, LIU T, et al. Seneca Valley virus suppresses host type I interferon production by targeting adaptor proteins MAVS, TRIF, and TANK for cleavage[J]. Journal of Virology, 2017, 91(16): e00823-17.
[53] LIU T, LI X, WU M, et al. Seneca Valley virus 2C and 3Cpro induce apoptosis via mitochondrion-mediated intrinsic pathway[J]. Frontiers in Microbiology, 2019, 10: 1202. doi: 10.3389/fmicb.2019.01202.
[54] XUE Q, LIU H, ZHU Z, et al. Seneca Valley Virus 3C protease negatively regulates the type I interferon pathway by acting as a viral deubiquitinase[J]. Antiviral Research, 2018, 160: 183-189. doi: 10.1016/j.antiviral.2018.10.028
-
期刊类型引用(9)
1. 吴超，周紫静，黄锦铧，许啸寅，邱洪，彭业萍. 基于长短时记忆的农作物生长环境数据预测. 深圳大学学报(理工版). 2024(05): 563-573 . 
百度学术
2. 谢妍，贾志军，王二英，贾永国，张旭东. 基于IPSO-BiLSTM和SVM的灌区农田土壤含水率预测研究. 节水灌溉. 2024(12): 63-71 . 百度学术
3. 刘嘉豪，刘浩. 基于LSTM神经网络模型的石油单井产量预测. 石油化工应用. 2023(03): 38-41+52 . 百度学术
4. 王登亮，马创举，曹锦萍，吴雪珍，吴群，刘丽丽，陈骏，孙建城. 柑桔生产数字技术研究及应用进展. 中国南方果树. 2023(03): 228-235 . 百度学术
5. 刘丽贞，黄琪，迟殿委，方朝阳，楚明航. 基于随机森林回归模型和高频数据的鄱阳湖子湖电导率预测. 水电能源科学. 2023(10): 50-53 . 百度学术
6. 代秋芳，熊诗路，李震，宋淑然，陈梓蔚，王元. 基于LSTM的柑橘幼苗蒸发量预测. 湖南农业大学学报(自然科学版). 2023(06): 743-747 . 百度学术
7. 彭旭，饶元，乔焰. 基于宽度卷积神经网络的异常农情数据检测方法. 华南农业大学学报. 2022(02): 113-121 . 本站查看
8. 高鹏，杨明欣，周平，孙道宗，谢家兴，陆健强，王卫星. 基于物联网和Deep-LSTM的茶树净光合速率动态预测模型. 农业工程学报. 2022(04): 159-168 . 百度学术
9. 王武，杨海建，洪林，杨蕾，王敏，李霜. 柑橘需水规律研究进展及趋势展望. 南方农业. 2022(19): 141-147+152 . 百度学术
其他类型引用(4)
