磁流作用下网络神经元中簇放电的模式.pdf

呼吸节呼吸节律是哺乳动物重要的生理性活动，其产生机制一直受到广泛的关注。关于正常呼吸节律的形成机制，律是哺乳动物重要的生理性活动，其产生机制一直受到广泛的关注。关于正常呼吸节律的形成机制，目前主要为起步细胞学说和神经网络学说。其中，神经元网络学说认为依赖于呼吸神经元之间的复杂相互目前主要为起步细胞学说和神经网络学说。其中，神经元网络学说认为依赖于呼吸神经元之间的复杂相互作用。本文，我们主要研究磁流对吸气周期和呼气周期的影响，当呼吸网络中考虑磁流的影响时，网络中作用。本文，我们主要研究磁流对吸气周期和呼气周期的影响，当呼吸网络中考虑磁流的影响时，网络中的呼吸节律将会表现出许多特性。由于哺乳动物的脑干呼吸中枢模式发生器的呼吸节律将会表现出许多特性。由于哺乳动物的脑干呼吸中枢模式发生器(CPG)具有许多对生存至关重具有许多对生存至关重要的功能，而要的功能，而pre-Btzinger复合体作为一种复合体作为一种CPG，是呼吸节律产生的关键部位。本文研究在磁流作用下，是呼吸节律产生的关键部位。本文研究在磁流作用下，pre-Btzinger复合体网络系统的节律转迁，并从动复合体网络系统的节律转迁，并从动力学的角度研究节律转迁的机制。结果表明磁流反馈力学的角度研究节律转迁的机制。结果表明磁流反馈系数系数k0对呼吸网络系统的节律模式和周期都具有重要影响。结果对呼吸节律的研究提供有益的参考。对呼吸网络系统的节律模式和周期都具有重要影响。结果对呼吸节律的研究提供有益的参考。关键词关键词 呼吸中枢模式发生器呼吸中枢模式发生器，分岔分岔，外激励外激励，吸气周期吸气周期，呼气周期呼气周期 The Pattern of Bursting in Network Neurons under Magnetic Hao Xu*,Xueli Chen School of Science,North China University of Technology,Beijing Received:Apr.10th,2023;accepted:Aug.21st,2023;published:Aug.28th,2023 Abstract Respiratory rhythm is an important physiological activity in mammals,and its production me-chanism has been widely concerned.Regarding the normal respiratory rhythm formation me-chanism,the main theory is the pacemaker theory and the network theory.Among them,the theory of neuronal networks believes that it relies on complex interactions between respiratory neurons.In this paper,the effects of magnetic flow are considered in the respiratory network,the *通讯作者。徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 53 生物物理学 respiratory rhythm in the network will exhibit many properties.As the brainstem respiratory cen-tral pattern generator(CPG)in mammals has many functions that are essential for survival,and the pre-Btzinger complex,as one type of CPG,is a key site of respiratory rhythm generation.In this paper,we focus on the effect of magnetic currents on the inspiratory and expiratory cycles,and when the effect of magnetic currents is taken into account in a respiratory network,the rhythmic transition of the pre-Btzinger complex network system under the action of magnetic currents is studied,and the mechanism of the rhythmic transition is investigated from a kinetic point of view.The results suggest that magnetic flow feedback coefficient k0 has an important influence on both the rhythmic pattern and the period of the respiratory network system.The results provide a bene-ficial reference for the study of respiratory rhythm in neural system.Keywords CPG,Bifurcation,External Excitation,Inspiratory Period,Expiratory Period Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 呼吸节律是由位于延髓中起步型神经元的节律性兴奋引起的，已有的研究成果表明，前包钦格复合体是呼吸节律活动产生最主要的部位1 2。Lieske 等人发现脑干切片中的 pre-Btzinger 复合体可以产生两种与正常呼吸极为相似的节律，该节律是多种呼吸模式的基础3。最近研究表明，在适当的外部电流刺激下，神经元的簇放电节律会发生改变，例如，磁流可以直接影响神经元的动态特性4 5 6，这种节律的控制模式对治疗节律相关的精神类疾病有一定的参考价值7 8 9。马军等人发现，通过电荷和磁通可以分别描述电场与电磁感应的变化，并以此为基础提出了几个神经元控制模型10 11。对于许多物种的脑干和脊髓中枢模式发生器(CPG)中的神经元网络，也都会表现出有规律有节奏的活动。其中，CPG 是指产生节律运动活动的神经环路，在许多中枢模式发生器中，这些模式由不同神经元群的顺序激活组成，而这些神经元群通过突触连接相互作用12。虽然有许多研究已经探索了呼吸系统中pre-Btzinger 复合体神经元放电模式的产生机制13 14，并引入磁流和电流条件下对 pre-Btzinger 复合体放电模式的产生进行了讨论，但仍有一些关键问题值得进一步研究，如：可放电的神经元嵌入到通过突触相互作用的完整网络中时，对引入磁流之后的放电模式会产生怎样的影响？若将 pre-Btzinger 神经元直接嵌入到网络中，它的放电模式是否会影响这些神经元嵌入后的呼吸网络中的放电机制？这些问题还需要进一步探究。呼吸节律的鲁棒性是指在各种条件下保持输出模式，或者在表示这些条件的模型参数或特征的变化之间保持输出模式。而且可以出现一种更强的鲁棒性的形式，其中不仅保留输出模式，而且还会保留产生该模式的动力机制。根据已有的实验表明，在甲壳类动物中，CPG 的节律性具有显著的鲁棒性，它会对温度等外在因素以及离子电流等内在因素的变化产生一定的消化节律15 16 17。本文中，我们运用呼吸 CPG 的高度简化模型，由代表重要神经元群体18的少量单位组成。在该模型网络中会产生三阶段呼吸节律，包括吸气神经元，吸气后神经元和呼气神经元19 20 21。我们引入磁流，对呼吸网络中具有放电模式的神经元进行建模，该微分方程支持已建立的簇放电产生机制，而且已在这些神经元中得到很好的实验记录22 23 24。我们发现网络的三阶段呼吸节律在网络中会受到磁流的严重影响。虽然外部刺激对神经元中混合簇放电的产生已有一定的研究，但是电流和磁流等外激励Open AccessOpen Access徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 54 生物物理学 对混合簇放电的影响，其现象和动力学机制都极其复杂，值得进一步探究。本文结构如下：第二部分给出了 pre-Btzinger 复合体兴奋性神经元和呼吸 CPG 网络的简化常微分方程模型；第三部分，研究磁流调控对神经元放电模式的影响，通过相平面分析，给出了网络系统的节律的产生和转迁的动力学机制；第四部分给出结论和展望。本文中的相平面分析是通过 XPPAUT 24软件计算得到。2.模型介绍模型介绍 pre-Btzinger 复合体兴奋性神经元和呼吸 CPG 网络的简化常微分方程模型如方程(1)(6)。该呼吸CPG 网络包括吸气前神经元(pre-I)，早期吸气神经元(early-I)，后期吸气神经元(post-I)，和增强型呼气神经元(aug-E)20 22。在这个网络中，我们假设每个神经元都能在静息与放电状态之间正常转换。对于编号为第一个神经元的 pre-I，引入磁流后的神经元模型的微分方程表达式如下：()()()()()()11111111dd,NaPKLsyn Isyn EFBCVtIV hIV nIVIVIVIC=+(1)()()()1111dd,hhthVhV=(2)其中，兴奋性神经元泄漏电流的方程表达式是()()1,1LL excNaIVgVE=，持续钠电流为()()()11,111,NaPNaP excNaIV hgmVhVE=，钾电流为()()()411,KK exciiKIV ngnVVE=；而模型中的其他神经元分别编号为：第二个神经元为 early-I，第三个神经元为 post-I，第四个神经元为 aug-E，且具体微分方程表达式如下：()()()()()(),dd,iNaP iiiK iiiL iisyn I iisyn E iiVtIV hIV nIVIVIVC=(3)()()()dd,iiihihthVhV=(4)其中2,3,4i。ih和in分别对应持续钠电导和钾电导的门控变量，抑制性神经元泄漏电流的方程表达式为()(),L iiL inhiNaIVgVE=，各单元的持续钠电流为()()(),NaP iiiNaP inhiiiNaIV hgmVhVE=，钾电流为()(),K iiiK inhiiKIV ngpVE=，对于抑制性单元，具有专用的门控变量 pi由电压相关的激活函数()(),ddiiout inhiip iptd fVp=表示。此外，动作电位跨此网络所有神经元时需要通过突触传递信号，抑制性突触电流(),syn E iiiIV h和兴奋性突触电流(),syn I iiiIV h的方程如下：()()(),inhisyn I iisyn Iisyn Ijiout inhjj NIVgVEb fV=(5)()()()44,1211,excisyn E iisyn Eisyn Ejiout excjiij NiiIVgVEa fVcc=+(6)exciN和inhiN分别表示突触到神经元单元 i 的兴奋性和抑制性单元集，对于 pre-I 兴奋性单元电流表现为；其中，常数 bji，aji分别表示抑制性和兴奋性连接的强度，而常数 cji，是对每个单元的驱动强度的参数。磁流激励的表达式是：()()203FBCeIkvv=+，k0为方程中的反馈系数，模型中涉及的其他函数和参数见附录。3.主要结果主要结果 3.1.IFBC对对呼吸节律周期呼吸节律周期的影响的影响 我们探索了磁流 IFBC对整个网络节律周期(T)、吸气阶段的持续时间(TI)和呼气阶段的持续时间(TE)的影响，即磁流强度的反馈系数 k0对 pre-I 单元输出周期的影响。由突触连接的呼吸网络放电模式如图 1 所示，黑色曲线表示 pre-I，红色曲线表示 early-I，蓝色曲线表示 post-I，粉色曲线表示 aug-E。这些网络放电模式在图 1(a)(c)显示了整个网络节律周期(T)、吸气阶段的持续时间(TI)和呼气阶段的持续时间(TE)，对于组成徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 55 生物物理学 网络的 pre-I，early-I，post-I 和 aug-E 各单元的输出，在 pre-I 阶段和 early-I 阶段的活动水平突然剧增之前，pre-I 阶段和 early-I 阶段的输出都会缓慢上升，这与吸气阶段的出现相对应。Pre-I 和 early-I 然后逐渐下降，使得 post-I 单元放电活动激活，并对这两个单元进行相应的抑制，标志着此时已经切换到呼气期。在这种转换之后，post-I 输出下降，对于 aug-E 的抑制减少，使其输出上升，产生所谓的呼气期，aug-E 输出并抑制但不完全抑制 post-I。随着磁流的反馈系数 k0的增加，整个网络节律周期(T)以及吸气(TI)和呼气(TE)阶段输出的周期会发生相应的改变。接下来，我们探索反馈系数 k0对包含 pre-I 单元的的网络中节律的频率以及周期的影响。当k0分别是 0，0.5，1 时，pre-I 神经元输出的整个网络节律周期(T)与呼气周期(TE)的持续时间变小，而吸气周期(TI)的持续时间不变，由突触连接的呼吸网络的簇放电个数随 k0的增加而增加。徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 56 生物物理学 Figure 1.Time series of membrane potential V for the entire network rhythm at a variation of 0k.(a)00k=;(b)00.5k=;(c)01k=图图 1.0k变化时整个网络节律的膜电位 V 时间序列。(a)00k=；(b)00.5k=；(c)01=k 当磁流的反馈系数 k0继续增加时，发现整个网络节律又出现了新的现象，周期(T)以及吸气(TI)和呼气(TE)阶段输出的周期与上述时间序列有所不同，如图 2 所示。当反馈系数 k0分别是 1.2，1.22，1.5 时，pre-I 神经元输出的整个网络节律周期(T)与吸气周期(TI)的持续时间变长，而的呼气周期(TE)持续时间不变，由突触连接的呼吸网络的簇放电个数随 k0的增加呈现先增加后减少的变换趋势。徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 57 生物物理学 Figure 2.Time series of membrane potential V for the entire network rhythm at a variation of 0k.(a)01.2k=;(b)01.22k=;(c)01.25k=图图 2.0k变化时整个网络节律的膜电位 V 时间序列。(a)01.2k=；(b)01.22k=；(c)01.25=k 整个网络节律周期 T(黑色)、吸气持续时间 TI(蓝色)和呼气持续时间 TE(红色)与在网络节律内反馈系数 k0之间的关系，如图 3 所示。随着磁流的反馈系数 k0的增加，整个网络节律先是通过调节呼气 TE来控制整个网络的呼吸周期。结果表明，由于磁流在呼吸节律性中起到了抑制作用，即周期 T 随着 k0的增加而减少。当 pre-Btzinger 复合体受到磁流刺激作用后，TE持续时间缩短，从最大值约 3000 ms 到最小值约 1000 ms，如图 3 所示；而 pre-I 单元随着磁流的反馈系数 k0从 1.0 增加到 1.25 时，TE持续时间处于平缓状态，而 TI持续时间开始急剧增加，此时磁流刺激是通过影响吸气 TI的持续时间来控制网络节律周期 T 的呼吸周期，增加 k0延长吸气时间，导致呼吸节律的周期增加，如图 3 所示。徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 58 生物物理学 Figure 3.The feedback function k0 varies through the network rhythm period(T)as well as the durations of the inspiratory phase(TI)and the expiratory phase(TE)图图 3.反馈系数 k0变化时整个网络节律周期(T)以及吸气阶段的持续时间(TI)和呼气阶段的持续时间(TE)3.2.IFBC对网络放电节律影响的动力学分析对网络放电节律影响的动力学分析 由前一个神经元释放信号，通过两个神经元接触的结构，后一个神经元的树突前端的轴突接收信号。因此，每个单元的膜电位方程取决于另一个单元的输出。在相平面中可以用可视化的动态模型单元来模拟网络中的每个神经元，用几何形式来表示网络的振荡机制。考虑投影到各二维相平面的网络输出，在这些相平面内，两个变量 hi，Vi的对应的零等值线，如图 4 所示。其中，此模型输出的 pre-I 单元的(h1,V1)相平面如图 4(a)所示，全系统的轨线(黑色曲线)也叠在此相平面上，蓝色实线是指 V1零等值线，右下方的零等值线对应于对 pre-I 的没有抑制作用的情况，而最左侧的零等值线对应于 pre-I 接收的最大抑制时的情况，而中间的零等值线是 pre-I 单元接收到较低的抑制水平，此时 pre-I 单元开始转至活动阶段。其中，黑色实心圆和方形分别对应 pre-I 单元抑制最低值和最高值的点。此外，红色曲线为 h1零等值线，与抑制水平无关。徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 59 生物物理学 Figure 4.Phase plane analysis of the neurons in the network at feedback coefficient 00k=.(a)(h1,V1);(b)(h3,V3);(c)(h4,V4)图图 4.反馈系数00=k时网络中各单元的相平面分析。(a)(h1,V1)平面；(b)(h3,V3)平面；(c)(h4,V4)平面 另外，对于 post-I 单元的(h3,V3)相平面上，最右侧的 Vpre-I零等值线对应于 post-I 的最小抑制，如图 4(a)中所示。图 4(b)中的曲线与图 4(a)类似的，而沿轨线获得 post-I 单元的零等值线 V3最大值时(黑色方块)，此时恰好与图 4(a)中的时间点相同。接下来，当对 post-I 单元达到最左侧零等值线数值，即达到对 post-I单元的最大抑制时，pre-I 也开始转变为低抑制水平(与图 4(a)中相同时间点用黑实心圈标记)，如图 4(b)所示。网络中全系统的轨线(黑色曲线)也叠在 aug-E 单元的相平面(h4,V4)上，如图 4(c)所示，从右到左，V4 零等值线依次对应于对 aug-E 的最小抑制以及到 pre-I 单元开始被激活时时对 aug-E 的抑制，此时，也与其他图相同时间点时用黑实心圈表示。其中，最右侧的零等值线 V1表示为达到峰值时对 aug-E 的抑制的曲线，而最左侧的零等值线表示为 V3达到峰值时(用黑色方块标记)对 aug-E 的抑制的曲线值。通过对网络呼吸节律中的周期内的几个关键时刻绘的相平面分析，对网络放电模式作出了解释。这些轨迹投影通常位于恰当的电压零等值线上。而我们在讨论中省略了 early-I 单元，这是由于 early-I 单元缺乏被激活向吸气阶段过渡时所需的自我调节能力。随着磁流的反馈系数 k0的增加，呼吸网络中的 pre-I 单元出现不同的放电状态，如图 5 所示，因此，也出现了非常丰富的神经元网络的动力学机制。不同的磁流度使得网络各单元都会出现放电振荡机制。与图 4(a)具有相似的曲线，右下方的零等值线对应于对 pre-I 的没有抑制作用的情况，而最左侧的零等值线对应于 pre-I 接收的最大抑制时的情况，而在 pre-I 单元的相平面中间的零等值线通过 h1I与 V1交点后，是其 pre-I 单元接收到较小的抑制水平，让其可以开始变成到活动阶段。引入磁流激励作用后，pre-I 单元被激活时的零等值线 Vpre-I的立方形状呈现单调的状态。Pre-I 单元的开始振荡唯一可能起作用的时刻是在被激活开始的时候，这取决于零等值线 h1与 V1交点处。然而，不同反馈系数 k0却使得 pre-I 单元具有相似的相平面，如图 5 所示。在磁流作用下，零等值线 V1在激活开始的时刻时是呈现单调状态的，而转变为呼气活动的机制与网络呼吸周期的持续时间对应(图 5(a)(e)。也就是说，即使在缺乏抑制的情况下，零等值线 Vpre-I在这一变化中呈现单调状态，使足够强的抑制可以产生较长的静息状态。当对 pre-I 单元的抑制在一段时间后开始衰减时，INap的失活可以允许 pre-I 单元不受抑制，这种机制是构成整个网络的呼吸节律的。当 pre-Btzinger 复合体嵌入网络后振荡时，整个网络周期 T 随着 k0的增加出现两种变化，一种是网络周期 T 缩短到一定值然后保持稳定，从初始状态 4.2 秒降到 2.1 秒左右，然后再增加到 6 秒；TE缩短到一定值然后保持稳定，从最大值 3.2 秒降到 1.3 秒左右保持稳定；另一种变化是，随着磁流的反馈系数 k0的增加，TI保持稳定值，当 k0增加到一定值时，TI突然增大，即 TI从保持稳定 1.8 秒上升到 4.5 秒左右。徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 60 生物物理学 Figure 5.Phase plane analysis of pre-I neurons in networks.(a)00.5k=;(b)01k=;(c)01.2k=;(d)01.22k=;(e)01.25k=图图 5.网络中 pre-I 单元的相平面分析。(a)00.5=k；(b)01=k；(c)01.2=k；(d)01.22=k；(e)01.25=k 4.结论结论 中枢模式发生器(CPG)对生存相关的节律性运动具有许多至关重要的功能，作为一种 CPG，徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 61 生物物理学 pre-Btzinger 复合体对哺乳动物的呼吸节律的输出，特别是在各种代谢和环境条件下持续存在是至关重要的。因此，即使 CPG 网络中的特定部分表现出输出驱动和神经自我调节的内在节律性，它们的放电模式也可能随着这些水平的变化而改变。所以，内在动力学变化下网络输出和调节的不变性可能是在各种条件下保持 CPG 稳定性能的关键特征。在呼吸网络的情况下，这种不变性确保了稳定地输出动态节律，这也是如此重要的稳态 CPG 的关键属性。本文主要研究了 pre-Btzinger 复合体作为 CPG 的一类，在磁流的作用下网络系统的节律改变及其动力学研究。本文中显示了磁流的反馈系数 k0在保持网络节律性的范围内变化时，周期、吸气持续时间和呼气持续时间与磁流的反馈系数 k0的值之间的关系，即随着磁流的反馈系数 k0的增加，pre-I 神经元输出的整个网络节律周期(T)先缩短到一定值后快速增加，呼气周期(TI)持续时间呈现出由稳定转为快速增加，而呼气周期(TE)的持续时间减少到一定值后趋于平缓，而且由突触连接的呼吸网络中簇的个数随 k0的增大而增加。参考文献参考文献 1 Smith,J.C.,Ellenberger,H.H.,Ballanyi,K.,Richter,D.W.and Feldman,J.L.(1991)Pre-Btzinger Complex:A Brains-tem Region That May Generate Respiratory Rhythm in Mammals.Science,254,726-729.https:/doi.org/10.1126/science.1683005 2 Rekling,J.C.and Feldman,J.L.(1998)PreBtzinger Complex and Pacemaker Neurons:Hypothesized Site and Kernel for Respiratory Rhythm Generation.Annual Review of Physiology,60,385-405.https:/doi.org/10.1146/annurev.physiol.60.1.385 3 Lieske,S.,Thoby-Brisson,M.,Telgkamp,P.and Ramirez,J.M.(2000)Reconfiguration of the Neural Network Con-trolling Multiple Breathing Patterns:Eupnea,Sighs and Gasps.Nature Neuroscience,3,600-607.https:/doi.org/10.1038/75776 4 Wang,H.and Chen,Y.(2016)Spatiotemporal Activities of Neural Network Exposed to External Electric Fields.Non-linear Dynamics,85,881-891.https:/doi.org/10.1007/s11071-016-2730-4 5 Zhan,F.and Liu,S.(2017)Response of Electrical Activity in an Improved Neuron Model under Electromagnetic Radiation and Noise.Frontiers in Computational Neuroscience,11,Article 107.https:/doi.org/10.3389/fncom.2017.00107 6 Ge,M.,Jia,Y.,Xu,Y.and Yang,L.(2017)Mode Transition in Electrical Activities of Neuron Driven by High and Low Frequency Stimulus in the Presence of Electromagnetic Induction and Radiation.Nonlinear Dynamics,91,515-523.https:/doi.org/10.1007/s11071-017-3886-2 7 Gluckman,B.J.,Nguyen,H.,Weinstein,S.L.and Schiff,S.J.(2001)Adaptive Electric Field Control of Epileptic Sei-zures.Journal of Neuroscience,21,590-600.https:/doi.org/10.1523/JNEUROSCI.21-02-00590.2001 8 Berzhanskaya,J.,Gorchetchnikov,A.and Schiff,S.J.(2007)Switching between Gamma and Theta:Dynamic Network Control Using Subthreshold Electric Fields.Neurocomputing,70,2091-2095.https:/doi.org/10.1016/j.neucom.2006.10.124 9 Durand,D.M.and Bikson,M.(2001)Suppression and Control of Epileptiform Activity by Electrical Stimulation:A Review.Proceedings of the IEEE,89,1065-1082.https:/doi.org/10.1109/5.939821 10 Liu,Y.,Ma,J.,Xu,Y.and Jia,Y.(2019)Electrical Mode Transition of Hybrid Neuronal Model Induced by External Stimulus and Electromagnetic Induction.International Journal of Bifurcation and Chaos,29,Article ID:1950156.https:/doi.org/10.1142/S0218127419501566 11 Zhang,Y.,Xu,Y.,Yao,Z.and Ma,J.(2020)A Feasible Neuron for Estimating the Magnetic Field Effect.Nonlinear Dynamics,102,1849-1867.https:/doi.org/10.1007/s11071-020-05991-y 12 Marder,E.and Calabrese,R.L.(1996)Principles of Rhythmic Motor Pattern Generation.Physiological Reviews,76,687-717.https:/doi.org/10.1152/physrev.1996.76.3.687 13 Del Negro,C.A.,Funk,G.D.and Feldman,J.L.(2018)Breathing Matters.Nature Reviews Neuroscience,19,351-367.https:/doi.org/10.1038/s41583-018-0003-6 14 Richter,D.W.and Smith,J.C.(2014)Respiratory Rhythm Generation in Vivo.Physiology,29,58-71.https:/doi.org/10.1152/physiol.00035.2013 15 Prinz,A.A.,Bucher,D.and Marder,E.(2004)Similar Network Activity from Disparate Circuit Parameters.Nature Neuroscience,7,1345-1352.https:/doi.org/10.1038/nn1352 徐浩，陈雪丽 DOI:10.12677/biphy.2023.113005 62 生物物理学 16 Marder,E.,Goeritz,M.L.and Otopalik,A.G.(2015)Robust Circuit Rhythms in Small Circuits Arise from Variable Cir-cuit Components and Mechanisms.Current Opinion in Neurobiology,31,156-163.https:/doi.org/10.1016/j.conb.2014.10.012 17 Tang,L.S.,Goeritz,M.L.,Caplan,J.S.,et al.(2010)Precise Temperature Compensation of Phase in a Rhythmic Motor Pattern.PLOS Biology,8,e1000469.https:/doi.org/10.1371/journal.pbio.1000469 18 Rubin,J.E.,Shevtsova,N.A.,Ermentrout,G.B.,et al.(2009)Multiple Rhythmic States in a Model of the Respiratory Central Pattern Generator.Journal of Neurophysiology,101,2146-2165.https:/doi.org/10.1152/jn.90958.2008 19 Richter,D.W.,Ballanyi,K.and Ramirez,J.M.(1996)Respiratory Rhythm Generation.In:Miller,A.D.,Bianchi,A.L.and Bishop,B.P.,Eds.,Neural Control of the Respiratory Muscles,CRC Press,Boca Raton.20 Richter,D.W.(1982)Generation and Maintenance of the Respiratory Rhythm.The Journal of Experimental Biology,100,93-107.https:/doi.org/10.1242/jeb.100.1.93 21 Del Negro,C.A.,Morgado-Valle,C.,Hayes,J.A.,et al.(2005)Sodium and Calcium Current-Mediated Pacemaker Neurons and Respiratory Rhythm Generation.Journal of Neuroscience,25,446-453.https:/doi.org/10.1523/JNEUROSCI.2237-04.2005 22 Yamanishi,T.,Koizumi,H.,Navarro,M.A.,et al.(2018)Kinetic Properties of Persistent Na+Current Orchestrate Oscil-latory Bursting in Respiratory Neurons.Journal of General Physiology,150,1523-1540.https:/doi.org/10.1085/jgp.201812100 23 Del Negro,C.A.,Koshiya,N.,Butera Jr.,R.J.and Smith,J.C.(2002)Persistent Sodium Current,Membrane Prop