1. 引言
近日节律(周期接近于24小时的节律)广泛存在于各种生物中,包括单细胞生物和哺乳动物等 [1]。在自然光照周期下,生物的生理活动和行为活动能够精确地与外界的24小时周期一致。令人有趣的是,在没有周期光照的情况下(如持续黑暗),生物仍然有显示出周期近似为24小时的节律的能力 [2] [3]。这种近似于24小时的周期被称为自振荡(或内源性)周期。但不同物种的自振荡周期存在一定的差异,主要分布在22至28小时之间 [4] [5]。从进化上看,如果生物的自振荡周期为24小时,将与地球的自转周期保持一致。依常理应是一个非常有利的生存机制,因为自振荡周期约接近24小时,越容易被外界周期所同步化以便于生物适应外界环境,从而更好地生存。但是为什么不是所有物种的自振荡周期都精确等于24小时?这种差异是如何产生的?科学家们对如何解释这种差异的存在引起了兴趣。
一方面从生物自身出发,Gu等人发现生物钟神经元之间的耦合强度的离散性能够导致这种差异性 [3] [6]。另一方面从生物外部因素出发,Daido通过模型为这种差异性给出了一个理论解释 [2]。该模型有如下假设条件:种群增长率取决于环境以及各物种之间的竞争。假设(1)物种和环境存在相位差,在最佳相位处(相位差为0时)种群增长率达到最大;假设(2)物种之间的竞争会导致增长率变小,两物种之间的竞争取决于他们的相位差。结合这两个因素,如果物种的自振荡周期都倾向于与外界周期一致,固然有利于得到最佳相位,但物种之间的相位差变小了从而使得竞争性变强,因此自振荡周期与外界周期一致有可能并不是最优的选择。为了得到最优的选择,在选择与外界环境周期一致和物种竞争之间应该存在某种平衡。
最近的实验发现,除了物种之间的自振荡周期不同之外,同一物种内部的各子群体还表现出自振荡周期不一样。例如人种、性别的不同会导致自振荡周期不一样。女性的自振荡周期为24.09 ± 0.2小时,男性的为24.19 ± 0.2小时 [7]。黑种人比其他人种的自振荡周期要短12.6分钟 [8]。即使同一个个体,自振荡周期也不是不变的。例如随着年龄的增加,自振荡周期会变小 [8] [9]。另外外界光照条件的转换等因素也会使自振荡周期发生改变 [8] [10]。然而迄今为止,尚未有论文解释同一物种内部的自振荡周期的差异性的、同一个个体自振荡周期不稳定的原因。
本文尝试着在Daido模型的基础上引入噪声来为这种差异性和不稳定提供一种解释。我们假设在物种的层面上来看的话,这种物种内部的差异性和个体不稳定可以看作是一种物种自振荡周期的噪声。基于此,全文的安排如下。第二部分是模型介绍部分,该模型在Daido模型的基础上增加了一项噪声项。第三部分使用数值模拟的方法来寻找噪声对种群规模的的影响。第四部分是结论和讨论部分。
2. 加噪声项的Daido模型介绍
Daido模型中每个物种有2个变量,包含种群密度
和时钟相位
[2]。含有N个物种的模型为:
(1)
(2)
其中i代表第i个物种,
,
和
,这里
。参数
都是正的;
是生物的自振荡频率,对应的自振荡周期为
。这里的环境频率
取为1。可以看出,该模型是Lotka-Volterra类型的方程(1) [11] 和Kuramoto类型的方程(2) [12] 的结合。本文假设,这里的任二个物种之间都相互耦合,即全联通网。
根据文献 [2] 的描述,对生物种群增长有影响的有两个因素,一是外界环境的影响,二是物种之间的竞争。外界的影响体现在:首先,位于方程(1)的A项用以描述假设(1),其中当
时与外界环境的相位一致从而被认为为最有利状态。参数a测度了这种环境作用的影响强度。其次,方程(2)右边的第二项是由于外部周期环境驱动,来促使种群相位朝最佳相位移动。这种驱动强度我们用r表示。物种之间的竞争影响也可以从两个方面体现:一个是来自方程(1)中的B项。该项用以描述假设(2),这种竞争效应由参数b来测度。当两物种的相位越接近时,竞争越激烈。另一个是方程(2)的最后一项,表示物种相位的排斥耦合。然而,不存在已知的实验数据来证明这种相互作用,因此
被设为零。为了简便,如果我们使用
来代替
,那么模型于就是一个关于
和
的方程,如下:
(3)
(4)
根据物种是否被外界所导引,所有的物种可分为两大类:能被外界导引的物种(ES,Entrained Species)和没被导引的物种(NES, Nonentrained Species) [2]。前者包括所有满足
的物种,即物种随着时间
都与环境建立了固定的相位关系
。其余的物种构成后者,他们的
都是一直变化的。为了方便起见,Daido把ES在
时的物种叫做中央导引物种(CES, Central ES),在导引的边缘也就是
最接近
或者
时叫做边缘导引物种(BES, Boundary ES)。如果a很大,说明环境很重要,CES的生存指数就会超过其他的物种。但是随着a值的降低,竞争会逐步变得相对重要,CES的生存指数会下降。
在Daido模型的基础上,本文在方程(2)中引入了噪声项
(方程(1)不变):
(5)
其中
是高斯白噪声,用于描述物种的自振荡周期不是稳定的、会有随机变化。
满足
和
,其中g是关键参数为噪声强度。
与文献 [2] 的参数取值一致,本文取
,
,
和
。假设物种的编号和其自振荡频率成正比,即对物种i来讲满足
。为了简便,下文我们直接用物种的序号i来代表其自振荡频率。根据能否被外界周期所导引,得到物种
为ES (其中两个BES和一个CES的序号为
),剩下的物种全都是NES。用
来估计每个物种的种群数量,其中
代表对时间求平均。由于s测量了物种的生存情况,因此我们可以把s的取值叫做生存指数:
代表着完美生存(没有物种比其生存指数更高);相反,
代表着此物种无法生存,即灭绝。生存指数取决于初始状态,本文通过1000次模拟来求每一物种的s的平均值,其中初值随机选取并满足
和
。
3. 数值模拟结果
在Daido模型中,研究了物种的自振荡频率(序号i)和该物种的生存指数(
)之间的关系受外界影响强度的值a的影响。与Daido论文一致,本文考察了a值分别为0.2,0.6,0.8和1.0的情况。为了描述的方便,将物种分为二类,即ES(
)或者NES(
和
)。其中ES进一步可以分为BES(
和
)和CES(
)。
3.1. 噪声强度对生存指数的影响
图1的每个子图(a)-(d)显示的是在给定的a值情况下,噪声强度g对自振荡频率(i)和生存指数(s(i))之间关系的影响。这里选取了四种不同的噪声强度g = 0,0.03,0.05和0.07,其中当g = 0时本文的模型退化为Daido模型。
当a很小为0.2时(图1(a)),在
的条件下(即无噪声的情况下),s(i)的值在NES处大于0但很小,在BES处s(i)约等于1,在ES的其他位置为0。因此,在BES处的物种占绝对的生存优势,NES和ES(不含BES)的物种不占优势。可以明显看出,所有的线都几乎重叠在一起,因此噪声强度g对s(i)的值影响不大,也就是说随着噪声强度的增加并不能改变BES的优势或者其他物种的劣势。
当a增加到0.6时(图1(b)),在
的条件下,s(i)在NES处等于0,在BES处略小于1,在ES的其他位置大于0。随着噪声强度g的增加,s(i)在NES处并没有明显的变化,在BES处略增加,在CES处变化不明显,ES的其他区域略增加。因此在a = 0.6时,噪声强度的增加会使ES处的s(i)的值略有增加(不含CES),但不影响
的定性关系。
当a增加到0.8时(图1(c)),在
的条件下,s(i)在NES处等于0,在BES处略大于0,在ES的其他位置大于0。随着噪声强度g的增加,s(i)在NES处并没有明显的变化,在ES处有所增加,其中在CES处增加最明显。因此在a = 0.8时,噪声强度的增加会使ES处的s(i)的值有所增加,但不影响
的定性关系。
当a增加到1时(图1(c)),在
的条件下,s(i)在NES处等于0,在BES处等于0,在CES处等于1,在ES的其他位置为0。随着噪声强度g的增加,s(i)在NES和BES处并没有明显的变化仍然为0,在CES处仍然等于1,但在靠近CES的位置上大于0。因此在a = 1时,噪声强度的增加会使靠近CES处的s(i)值有所增加,但不影响
的定性关系。
综上所述,噪声对s(i)值的影响取决于a值。当a等于0.2时,噪声对s(i)值没有影响;当a等于0.6时,噪声与ES处(除CES)的s(i)值成正相关;当a等于0.8时,噪声与ES处的s(i)值成正相关;当a等于1时,只有邻近CES的位置,噪声的增加才能使其s(i)值增加。
3.2. 噪声强度对标准方差的影响
从直觉上来说,如果噪声强度g的增加能够增加s(i)的值,那么s(i)对应的标准方差也应该按照比例增加。图2检验了这一直觉是否正确。令人惊讶地,除了在a = 0.2情况下(a),其他子图(b)-(d)显示当s(i) > 0时,噪声强度g的增加能够明显地减少生存指数的标准方差d(i)。这就说明噪声强度的增加,能够增加数值模拟的稳定性或者说能够减少生存指数对于初值的依赖性。
Figure 1. Effect of noise intensity g on species Survival Index s(i)
图1. 噪声强度g对物种生存指数s(i)的影响
Figure 2. Effect of noise intensity g on Standard Variance d(i) of species Survival Index
图2. 噪声强度g对物种生存指数的标准方差d(i)的影响
4. 总结与讨论
本文在Daido模型的基础上,考虑自振荡频率的扰动(白噪声)对各个物种生存的影响。我们假设用白噪声来代表个体自振荡周期的不稳定性和同一物种内不同子群体之间的自振荡周期差异性,并发现白噪声能够有利于(外界影响强度a的值适中,a = 0.6或0.8)或无害于(a = 0.2或1)物种的生存。实验表明,所有物种的自振荡周期是分布在24小时左右的正态分布,因此a = 0.8的情况可能接近于真实的情况。所以,我们建议现实中可能存在的白噪声有利于物种的生存。此外,本文发现随着噪声强度的增加会使生存指数对初值的依赖性减低。这里给本文的主要结论给出一个定性的解释。自振荡周期的白噪声的存在会使物种之间的相位差变得不固定,因此导致物种之间的竞争性变低,从而导致物种的生存指数提高。本文的研究有助于理解自然界中为什么会存在个体自振荡周期的不稳定性和同一物种内不同子群体之间的自振荡周期差异性。
基金项目
感谢上海理工大学科研项目课程(课程代码:13002100)和上海理工大学大学生创新创业项目“基于小世界网络的哺乳动物近日节律模型”(项目号:XJ2019141)。