
Citation: Kerri Glassner, Eamonn MM Quigley, Lissa Franco, David W Victor III. Autoimmune liver disease and the enteric microbiome[J]. AIMS Microbiology, 2018, 4(2): 334-346. doi: 10.3934/microbiol.2018.2.334
[1] | Erxi Zhu . Time-delayed feedback control for chaotic systems with coexisting attractors. AIMS Mathematics, 2024, 9(1): 1088-1102. doi: 10.3934/math.2024053 |
[2] | Ning Cui, Junhong Li . A new 4D hyperchaotic system and its control. AIMS Mathematics, 2023, 8(1): 905-923. doi: 10.3934/math.2023044 |
[3] | Binhao Hong, Chunrui Zhang . Bifurcations and chaotic behavior of a predator-prey model with discrete time. AIMS Mathematics, 2023, 8(6): 13390-13410. doi: 10.3934/math.2023678 |
[4] | Kexin Zhang, Caihui Yu, Hongbin Wang, Xianghong Li . Multi-scale dynamics of predator-prey systems with Holling-IV functional response. AIMS Mathematics, 2024, 9(2): 3559-3575. doi: 10.3934/math.2024174 |
[5] | A. E. Matouk . Chaos and hidden chaos in a 4D dynamical system using the fractal-fractional operators. AIMS Mathematics, 2025, 10(3): 6233-6257. doi: 10.3934/math.2025284 |
[6] | Sukono, Siti Hadiaty Yuningsih, Endang Rusyaman, Sundarapandian Vaidyanathan, Aceng Sambas . Investigation of chaos behavior and integral sliding mode control on financial risk model. AIMS Mathematics, 2022, 7(10): 18377-18392. doi: 10.3934/math.20221012 |
[7] | Sulasri Suddin, Fajar Adi-Kusumo, Mardiah Suci Hardianti, Gunardi . Bifurcation analysis of a diffuse large b-cell lymphoma growth model in germinal center. AIMS Mathematics, 2025, 10(5): 12631-12660. doi: 10.3934/math.2025570 |
[8] | Xiaoyan Hu, Bo Sang, Ning Wang . The chaotic mechanisms in some jerk systems. AIMS Mathematics, 2022, 7(9): 15714-15740. doi: 10.3934/math.2022861 |
[9] | Hany A. Hosham, Thoraya N. Alharthi . Bifurcation and chaos in simple discontinuous systems separated by a hypersurface. AIMS Mathematics, 2024, 9(7): 17025-17038. doi: 10.3934/math.2024826 |
[10] | Na Min, Hongyang Zhang, Xiaobin Gao, Pengyu Zeng . Impacts of hunting cooperation and prey harvesting in a Leslie-Gower prey-predator system with strong Allee effect. AIMS Mathematics, 2024, 9(12): 34618-34646. doi: 10.3934/math.20241649 |
Chaos theory is a field of study that involves the mathematical analysis of complex deterministic systems exhibiting unstable aperiodic behavior. In chaotic systems, even a tiny difference in the initial conditions can lead to significantly different outcomes over time. As a way of understanding and describing the long-term behavior of a chaotic system, a chaotic attractor has at least one positive Lyapunov exponent [1]. Bifurcation diagrams can provide insights into the mechanisms behind the chaotic behavior. After transients decay, bifurcation diagrams can be constructed by recording either the Poincaré section intersections or the local peaks of the system's trajectory. By understanding the underlying mechanisms, it is possible to control and predict the behavior of the system. Amplitude control refers to the ability to regulate the magnitude of a signal, typically in the context of chaotic systems. It was shown in [2] that symmetric chaotic systems can be both multistable and also have an independent amplitude control parameter. Chaos has found applications across various fields, including economics [3], cryptography [4], and robotics [5], among others.
In their work on Chua circuits, Leonov, Kuznetsov et al. [6,7,8] introduced the classification of attractors as hidden or self-excited, and they proposed an effective procedure for localizing hidden attractors based on homotopy and numerical continuation. The classification is a crucial concept, which represents a contemporary extension of Andronov's theory of oscillations and provides a fundamental understanding of the existence and characteristics of hidden oscillatory behavior in chaotic models [9,10].
A self-excited attractor refers to an attractor whose basin of attraction intersects an open neighborhood of an equilibrium. It can be numerically found by plotting the trajectory of the system with random initial conditions near the equilibria. On the other hand, a hidden attractor is an attractor whose basin of attraction does not intersect with small neighborhoods of any equilibrium points. This implies that it cannot be numerically found using random initial conditions and demands special analytical-numerical procedures to locate its basin of attraction. In practice, there is a need to control the hidden attractors, because some structures (such as bridges or airplane wings) can display catastrophic responses under disturbances [11,12]. According to [13], a saddle focus of index 2 is defined by an equilibrium with three eigenvalues: one real $ \gamma $ and a pair of complex conjugates $ \sigma \pm j\omega $, where $ \gamma < 0 $ and $ \sigma > 0 $. Shilnikov homoclinic and heteroclinic theorems provide sufficient conditions for the existence of self-excited chaotic attractors in 3D systems with unstable hyperbolic equilibria. Saddle-foci with an index of 2 are frequently encountered in such systems. As noted by Zhou and Chen [14], the chaotic attractors of Lorenz, Rössler, Chua, Chen, and many other systems are in the sense of Shilnikov, but not all chaotic attractors are of Shilnikov sense. Beyond the scope of Shilnikov theorems, many chaotic flows with various types of non-hyperbolic equilibria were introduced [15,16,17]. Li and Hai [18] introduced a three-dimensional robust chaotic system with two non-hyperbolic equilibrium points, exhibiting features of amplitude and position modulation. It is shown that the presence of linear terms is not essential in producing complex behaviors [19].
In a review paper, Jafari et al. [20] noted that dissipative hidden chaotic flows can be categorized into three types: those with a stable equilibrium [21,22], those without equilibrium [23], and those with an infinite number of equilibria [24]. Molaie et al. [25] proposed twenty-three chaotic flows with a single stable equilibrium point. Jafari and Sprott [26] proposed nine chaotic flows with hidden attractors that have a line equilibrium. Kumarasamy et al. [27] investigated how saddle-node bifurcations contribute to the emergence of hidden attractors, both in systems with stable equilibrium points and those without equilibrium points. Antimonotonicity refers to the intertwined creation and destruction of periodic orbits in nonlinear dynamical systems, leading to inevitable reversals of bifurcation cascades [28]. For a family of jerk systems with a hyperbolic tangent function, Li et al. [29] illustrated that antimonotonicity is important in elucidating the generation mechanism of hidden attractors with a stable equilibrium. Recently, Dong et al. [30] investigated a chaotic system that has no equilibrium through numerical simulations and circuit implementations, and discussed engineering applications such as synchronization, offset boosting control, and image encryption.
A chameleon system is a chaotic system that, as parameter values vary, can exhibit various types of chaotic attractors, including self-excited attractors and the three types of hidden attractors. This feature provides flexibility in applications needing different types of chaotic signals. In a review paper, Sprott [31] categorized chaotic flows from a quadratic family based on equilibrium points, in which both self-excited and hidden attractors are involved. Jafari et al. [32] introduced a chameleon system with quadratic nonlinearities, contributing to an enhanced understanding of chaotic attractors, especially those that are hidden. Recall that a hyperchaotic system is defined as a dynamical system that exhibits chaotic behavior, characterized by the presence of multiple positive Lyapunov exponents. Motivated by [32], Rajagopal et al. [33] presented a hyperchaotic chameleon system, contributing to a deeper understanding of hidden chaotic flows of higher dimensions. Natiq et al. [34] presented a chameleon system with a trigonometric nonlinear controller. Cang et al. [35] investigated the hidden and self-excited coexisting attractors in the Yang system [36]. Signing et al. [37] introduced a chameleon cryptosystem that combines chameleon chaotic systems with dynamic DNA coding to enhance image encryption security through the use of chaotic sequences. Fan et al. [38] investigated the two-parameter space bifurcation of the Wei chameleon system. Recently, Tiwari et al. [39] introduced a class of quadratic chameleon systems that allow for the coexistence of hidden chaotic attractors and conservative tori.
A Hopf bifurcation is a local bifurcation where an equilibrium of a dynamical system changes stability as a parameter is varied. This phenomenon occurs when the equilibrium possesses a pair of purely imaginary eigenvalues, accompanied by specific transversality conditions and an absence of zero eigenvalues. The bifurcation can be classified as either supercritical or subcritical, resulting in either a stable or an unstable limit cycle within an invariant two-dimensional manifold, respectively. The latter situation presents a potential risk, as stable, large-amplitude limit cycles can coexist with the stable equilibrium point [40]. Stankevich et al. [41] and Zhao et al.[42] studied the dynamics of the Chua systems, establishing a connection between subcritical Hopf bifurcations of equilibria and the emergence of hidden attractors. Using a family of quadratic jerk systems, Liu et al. [43] also showed that subcritical Hopf bifurcations are closely related to the emergence of hidden attractors. Yang et al. [44] proposed a 7D hyperchaotic system with five positive exponents and investigated Hopf bifurcation to understand how complex dynamics emerge in this system. Li and Chen [45] demonstrated the complex neuromorphic behaviors of a novel current-controlled Chua corsage memristor-based neuron circuit through Hopf bifurcation and edge-of-chaos analysis. Liu et al. [46] explored a 4D hyperchaotic system through Hopf bifurcation analysis and the linear control method.
The design of a chaotic system with independent non-bifurcation parameters faces a significant challenge due to the broadband random-like nature of chaos [47,48]. Signal polarity refers to whether a signal oscillates around zero (bipolar) or around a non-zero reference level (unipolar). Offset boosting involves adjusting the average value of a state variable by introducing a constant term into the system model, resulting in signal offset and enabling control of signal polarity. Li et al. [49] proposed an offset boosting-based strategy for the identification of multistability in dynamical systems, employing nonbifurcation operations for diagnosing multistability. In chaotic and hyperchaotic systems [50,51], offset boosting controls have been implemented, enabling direct offset regulation through the use of a single constant.
The paper is organized as follows. In Section 2, we introduce a novel chameleon system and conduct a detailed study on its symmetry, dissipation, stability analysis, Hopf bifurcation, and chaotic dynamics with $ \mu $ varying. In Section 3, we perform a bifurcation analysis for scenarios 1 to 4 with $ \mu = 0 $. We investigate the system's behavior under various parameter values, focusing on the emergence of different dynamical phenomena. In Section 4, we investigate hidden chaotic bistability with offset boosting through the variation of $ \mu $. In the final section, we synthesize our key findings and discuss their implications.
Consider the quadratic system
$ {dxdt=y,dydt=−x+y(z+μ),dzdt=f(x,y,z), $
|
(2.1) |
where
$ f(x,y,z)=a−bz+cy2−0.7x2−3xy. $
|
It is symmetric with respect to 180-degree rotation about the $ z $-axis. This means that the system is invariant with respect to the transformation $ (x, y, z)\rightarrow (-x, -y, z) $, corresponding to a 180-degree rotation about the $ z $-axis.
Let $ h = (h_{1}, h_{2}, h_{3}) $ be the vector field of system (2.1). Thus the divergence function is
$ ∇⋅V=h1x+h2y+h3z=−b+(z+μ). $
|
It is clear that system (2.1) is not always dissipative as the divergence depends on $ z $. However, in specific regions of the phase space where the condition $ \nabla \cdot {V} < 0 $ holds on average, we have a dissipative flow. Due to the symmetry, system (2.1) can have a symmetrical attractor or a symmetric pair of coexisting attractors.
If $ a = b = 0 $, system (2.1) has a line equilibrium along the $ z $-axis. This means that every point $ (0, 0, z_{0}) $ on the axis is an equilibrium point.
If $ a\neq 0 $ and $ b = 0 $, we have $ f(0, 0, z) = a\neq 0 $. Hence, there is no equilibrium point. When $ x = y = 0 $, we have $ \dfrac{dx}{dt} = \dfrac{dy}{dt} = 0 $ and $ \dfrac{dz}{dt} = a\neq 0 $. This indicates that the $ z $-axis is invariant under the system's flow. The direction of the flow along the $ z $-axis is determined by the sign of $ a $: if $ a > 0 $, the flow is in the forward direction; if $ a < 0 $, the flow is in the backward direction.
If $ b\neq 0 $, there exists a unique equilibrium $ P: (0, 0, a/b) $. The Jacobian matrix of system (2.1) at $ P $ is
$ A=[010−1μ+ab000−b]. $
|
The corresponding characteristic polynomial is
$ φA(λ)=λ3+p1λ2+p2λ+p3=(λ+b)(λ2−(μ+ab)λ+1). $
|
(2.2) |
Proposition 2.1. Consider system (2.1) with $ b\neq 0 $. It is locally asymptotically stable around $ P $ if
$ b>0anda<−bμ. $
|
(2.3) |
Proof. From (2.2) and Routh-Hurwitz criterion, the equilibrium $ P $ of system (2.1) is asymptotically stable if $ b > 0 $ and $ \mu+\dfrac{a}{b} < 0 $. Such constraints are equivalent to the inequalities stated in (2.3). Thus the conclusion holds.
Corollary 2.1. Consider system (2.1) with $ b > 0 \; \mathit{\text{and}} \; a < -b\mu $, and let
$ Δ=(μ+ab)2−4. $
|
We have the following conclusions.
(i) If $ \Delta < 0 $, then the equilibrium $ P $ is a stable node-focus.
(ii) If $ \Delta\geq 0 $, then the equilibrium $ P $ is a stable node.
Proof. (ⅰ) Assume that $ \Delta < 0 $. From (2.2), we know that the characteristic polynomial has a negative root and a pair of complex conjugate roots with a negative real part. Thus the equilibrium $ P $ is a stable node-focus.
(ⅱ) Assume that $ \Delta\geq 0 $. From (2.2), we know that the characteristic polynomial has three negative roots, counted with multiplicity. Therefore, the equilibrium $ P $ is a stable node.
Assume $ a, b\neq 0 $ and $ \mu = 0 $. Then system (2.1) can assume all types of hyperbolic equilibrium points at $ P:\; (0, 0, a/b) $.
Corollary 2.2. Consider system (2.1) with $ b\neq 0 $ and $ \mu = 0 $. The following statements hold.
(i) It is locally asymptotically stable around $ P $ if $ a < 0 $ and $ b > 0 $. In this case, if $ a > -2b $, then the equilibrium $ P $ is a stable node-focus; if $ a\leq -2b $, then the equilibrium is a stable node.
(ii) If $ a > 0 $ and $ b > 0 $, the equilibrium $ P $ is unstable. In this case, if $ a < 2b $, the equilibrium $ P $ is a saddle-focus of index 2; if $ a\geq 2b $, the equilibrium $ P $ is a saddle of index 2.
(iii) If $ b < 0 $, the equilibrium $ P $ is unstable. In this case, if $ 0 < a < -2b $, the equilibrium $ P $ is a saddle-focus of index 1; if $ a\geq -2b $, the equilibrium $ P $ is a saddle of index 1; if $ 2b < a < 0 $, the equilibrium $ P $ is an unstable node-focus; and if $ a\leq 2b $, the equilibrium $ P $ is an unstable node.
Proof. (ⅰ) By setting $ \mu = 0 $, the statements follow from Proposition 2.1 and Corollary 2.1.
(ⅱ) Setting $ \mu = 0 $, then the characteristic polynomial (2.2) becomes
$ φA(λ)=(λ+b)(λ2−abλ+1). $
|
(2.4) |
Since $ a > 0 $ and $ b > 0 $, the quadratic factor of (2.4) can have a root with a positive real part. Thus the equilibrium $ P $ is unstable.
Assume that $ a, b > 0 $. If $ a < 2b $, the cubic polynomial (2.4) has a negative root and a pair of complex conjugate roots with a positive real part. Therefore, the equilibrium $ P $ is a saddle-focus of index 2. If $ a\geq 2b $, the cubic polynomial (2.4) has a negative root and two positive roots. Therefore, the equilibrium $ P $ is a saddle of index 2.
(ⅲ) If $ b < 0 $, the characteristic polynomial (2.4) has a positive root $ \lambda = -b $. Thus the equilibrium $ P $ is unstable. The rest of the proof for classifying the equilibrium is similar to that of case (ⅱ).
When $ a = 0 $, system (2.1) becomes
$ {dxdt=y,dydt=−x+y(z+μ),dzdt=−bz+cy2−0.7x2−3xy. $
|
(2.5) |
It has a unique equilibrium point at the origin. Here we assume that $ b > 0 $.
Proposition 2.2. When $ \mu = 0 $, a Hopf bifurcation occurs at the origin for system (2.5).
Proof. The characteristic polynomial of this system is
$ φA(λ)=λ3+p1(μ)λ2+p2(μ)λ+p3(μ), $
|
where
$ p1=−μ+b,p2=−bμ+1,p3=b. $
|
Note that
$ (p1p2−p3)(μ)=bμ2−(1+b2)μ. $
|
So
$ (p1p2−p3)(0)=0, $
|
(2.6) |
$ (p1p2−p3)′(0)=−(1+b2)<0. $
|
(2.7) |
Thus, a Hopf bifurcation occurs from the origin at $ {\mu} = 0 $.
When $ \mu = 0 $, system (2.5) becomes
$ {dxdt=y,dydt=−x+F2,dzdt=−bz+F3, $
|
(2.8) |
where
$ F2=yz,F3=cy2−0.7x2−3xy. $
|
According to [52], the first Lyapunov number $ a_{1} $ at the Hopf bifurcation can be obtained as
$ a1=−216λ3F2yz(F3xx+F3yy)+F2yz16(4ω2+λ23)(λ3(F3xx−F3yy)+4ωF3xy), $
|
(2.9) |
where $ \lambda_{3} = -b $ and $ \omega = -1 $.
By a simple computation, we have
$ a1=(30c−7)b2+60b+(80c−56)80b(b2+4). $
|
(2.10) |
Let us consider the following system:
$ {dxdt=y,dydt=−x+y(z+μ),dzdt=−bz+y2−0.7x2−3xy, $
|
(2.11) |
where $ b > 0 $.
Corollary 2.3. When $ \mu = 0 $, a subcritical Hopf bifurcation occurs at the origin for system (2.11), which generates one unstable limit cycle when $ {\mu} $ is negative but not too small. For $ \mu = 0 $, the origin is unstable.
Proof. Note that $ b > 0 $ and system (2.11) is a special case of (2.5) with $ c = 1 $.
By setting $ c = 1 $ in (2.10), we have
$ a1=23b2+60b+2480b(b2+4)>0, $
|
(2.12) |
where $ a_{1} $ is the first Lyapunov number. According to Proposition 2.2, conditions (2.6) and (2.7), and the sign of (2.12), we have the conclusions.
The following results give us insights into the stability of the system (2.5).
Proposition 2.3. Consider system (2.5) with $ \mu\neq 0 $ and $ b > 0 $. The stability properties of the origin are characterized as follows:
(i) If $ -2 < \mu < 0 $, the origin is a stable node-focus.
(ii) If $ 0 < \mu < 2 $, the origin is a saddle-focus of index 2.
(iii) If $ \mu\leq -2 $, the origin is a stable node.
(iv) If $ \mu\geq 2 $, the origin is a saddle of index 2.
Proof. The Jacobian matrix of the system at the origin is
$ A=[010−1μ000−b]. $
|
The corresponding characteristic polynomial is
$ φA(λ)=(λ+b)(λ2−μλ+1). $
|
(2.13) |
By analyzing of the roots of (2.13), we can draw the four conclusions.
Consider the following three systems:
$ {dxdt=y,dydt=−x+y(z+μ),dzdt=y2−0.7x2−3xy, $
|
(2.14) |
$ {dxdt=y,dydt=−x+y(z+μ),dzdt=0.001+y2−0.7x2−3xy, $
|
(2.15) |
$ {dxdt=y,dydt=−x+y(z+μ),dzdt=−0.01z+y2−0.7x2−3xy, $
|
(2.16) |
where $ \mu\in [-0.08, 0.02] $. These systems can be obtained from system (2.1) by setting $ (a, b, c) = (0, 0, 1) $, $ (a, b, c) = (0.001, 0, 1) $, and $ (a, b, c) = (0, 0.01, 1) $, respectively. Figure 1 shows the bifurcation diagrams for these systems with respect to the parameter $ \mu $.
Let $ (L_1^{(k)}, L_2^{(k)}, L_3^{(k)}), \; k = 1, 2, 3 $, be the Lyapunov exponent spectra of systems (2.14), (2.15), and (2.16), respectively. For each $ k $, the exponents satisfy $ L_1^{(k)} > 0 $, $ L_2^{(k)} = 0 $, and $ L_3^{(k)} < 0 $. The variations of $ L_{1}^{(k)} $ and $ L_{3}^{(k)} $ with respect to the parameter $ \mu $ for $ 1\leq k\leq 3 $ are shown in Figure 2.
Based on the two figures above and the nature of equilibrium points, we classify the dynamical behaviors of these three systems into four scenarios:
Scenario 1: Hidden chaotic dynamics with a line equilibrium Note that system (2.14) has a line equilibrium $ x = y = 0 $. Thus, from panel (a) of Figure 1, we see that system (2.14) exhibits a symmetric pair of hidden chaotic attractors for $ \mu\in [-0.08, 0.02] $, coexisting with a line equilibrium.
Scenario 2: Hidden chaotic dynamics without equilibrium It is clear that there is no equilibrium in system (2.15). Thus, from panel (b) of Figure 1, we see that system (2.15) exhibits a symmetric pair of hidden chaotic attractors without an equilibrium point for $ \mu\in [-0.08, 0.02] $.
Scenario 3: Hidden chaotic dynamics with a stable equilibrium By Proposition 2.3, we know that system (2.16) has a stable equilibrium at the origin when $ \mu\in [-0.08, 0) $. Thus, from panel (c) of Figure 1, we see that system (2.16) exhibits a symmetrical hidden chaotic attractor for $ \mu\in [-0.08, 0) $.
Scenario 4: Self-excited chaotic dynamics with an unstable equilibrium By Corollary 2.3 and Proposition 2.3, we know that system (2.16) has an unstable equilibrium at the origin when $ \mu\in [0, 0.02] $. Thus, from panel (c) of Figure 1, we see that system (2.16) exhibits a symmetrical self-excited chaotic attractor for $ \mu\in [0, 0.02] $.
Using two fixed initial conditions and varying parameters $ c $ and $ a $, we demonstrate the system's period-doubling routes to chaos when $ \mu = 0 $. The analysis reveals that system (2.1) can exhibit either self-excited or hidden chaotic attractors under different parameter settings.
Let
$ f1(x,y)=y2−0.7x2−3xy, $
|
(3.1) |
$ f2(x,y)=1.008y2−0.7x2−3xy. $
|
(3.2) |
When $ \mu = 0 $ and $ f(x, y, z) = f_{1}(x, y) $, system (2.1) becomes
$ {dxdt=y,dydt=−x+yz,dzdt=y2−0.7x2−3xy. $
|
(3.3) |
It has a symmetric pair of hidden chaotic attractors with a line equilibrium, which is shown in Figure 3 by initiating with $ (x, y, z) = (\pm{0.8}, \pm{0.8}, 0) $. Figure 4 shows various projections of the coexisting hidden chaotic attractors. The corresponding Poincaré maps for the two attractors are shown in Figure 5. The set of dense points covering a certain region on the plane $ y = 0 $ shows that each attractor is chaotic. When the initial conditions are set to $ (0.8, 0.8, 0) $, system (3.3) has Lyapunov exponents of $ (L_{1}, L_{2}, L_{3}) = (0.0548, 0, -0.6914) $ and a Kaplan-Yorke dimension of $ D_{\mathrm{KY}} = 2-L_{1} / L_{3} \approx 2.079 $. The sum of the three exponents $ L_{1}+L_{2}+L_{3} = -0.6366 $ signifies that the system exhibits dissipative behavior, characterized by a chaotic attractor.
When $ \mu = 0 $ and $ f(x, y, z) = f_{2}(x, y) $, system (2.1) becomes
$ {dxdt=y,dydt=−x+yz,dzdt=1.008y2−0.7x2−3xy. $
|
(3.4) |
It has a symmetric hidden chaotic attractor with respect to the $ z $-axis, which is visualized in Figure 6 by starting from the initial conditions $ (x, y, z) = (0.8, 0.8, 0) $. Moreover, the line equilibrium $ z $-axis is denoted by a green line. The Poincaré maps of the attractor are shown in Figure 7, in which forward and backward directions are considered, respectively. When the initial conditions are set to $ (0.8, 0.8, 0) $, system (3.4) has Lyapunov exponents of $ (L_{1}, L_{2}, L_{3}) = (0.1179, 0, -0.7641) $ and a Kaplan-Yorke dimension of $ D_{\mathrm{KY}} = 2-L_{1} / L_{3} \approx 2.154 $.
Consider the following system:
$ {dxdt=y,dydt=−x+yz,dzdt=cy2−0.7x2−3xy, $
|
(3.5) |
where $ c\in [0.8, 1.008] $ is the coefficient of the nonlinear term $ y^2 $. It contains systems (3.3) and (3.4) as special cases.
For the two initial conditions $ (\pm{0.8}, \pm{0.8}, 0) $, the bifurcation diagram is presented in Figure 8(a) for $ c\in [0.94, 1.008] $. It shows the development of the hidden chaotic motions (with a line equilibrium) via period-doubling cascades. The periodic orbit undergoes a symmetry-breaking bifurcation around $ c \approx 0.967 $, transitioning into a saddle orbit, and giving rise to a pair of asymmetric stable periodic orbits. For each initial data, the first period-doubling bifurcation occurs at $ c\approx 0.991 $. The period and amplitude of the two asymmetric attractors grow with increasing $ c $ until they touch the saddle orbit and merge into a broader symmetrical hidden chaotic attractor at $ c\approx 1.0027 $ in an attractor merging crisis. Recall that an attractor merging crisis refers to two or more chaotic attractors merging to form a single attractor as the critical parameter value is passed [53]. For the initial data $ (x(0), y(0), z(0)) = (0.8, 0.8, 0) $, the variations of Lyapunov exponents and their sum with respect to the parameter $ c $ are shown in Figure 8(b). The Lyapunov exponent spectrum depicted in Figure 8(b) aligns with the bifurcation diagram illustrated in Figure 8(a). In Figure 8(b), we also plotted the variation of $ L_\mathrm{sum} = L_{1}+L_{2}+L_{3} $ with respect to $ c $. For each value of $ c $, the sum $ L_\mathrm{sum} $ is negative, indicating that the system is dissipative. Hence, the system is always dissipative for $ c $ within the interval $ [0.94, 1.008] $.
Letting $ b = \mu = 0 $ in system (2.1), we have the following system:
$ {dxdt=y,dydt=−x+yz,dzdt=a+(y2−0.7x2−3xy). $
|
(3.6) |
When $ a\neq 0 $, it has no equilibrium but does exhibit an invariant straight line along the $ z $-axis. Recall that when $ a = 0 $, the system can have a symmetric pair of hidden chaotic attractors with a line equilibrium ($ z $-axis).
For the two initial conditions $ (\pm{0.8}, \pm{0.8}, 0) $, the bifurcation diagram is presented in Figure 9(a) for $ a\in [-0.02, 0.004] $. It shows the development of the hidden chaotic motions (with no equilibrium) via period-doubling cascades. The left panel of Figure 9 with the change of $ a $ is qualitative similar to Figure 8(a) with the change of $ c $. For the initial data $ (x(0), y(0), z(0)) = (0.8, 0.8, 0) $, the variation of Lyapunov exponents with respect to the parameter $ a $ is shown in Figure 9(b). Note that at $ a\approx 0.0013 $, a symmetric pair of hidden chaotic attractors combine into a symmetrical hidden chaotic attractor. At $ a = 0.004 $, there exists a symmetrical hidden chaotic attractor with a Lyapunov exponent spectrum of $ (L_1, L_2, L_3) = (0.1162, 0, -0.7621) $ and a Kaplan-Yorke dimension of 2.152.
Letting $ b = 0.01, \mu = 0 $ in system (2.1), we have the following system:
$ {dxdt=y,dydt=−x+yz,dzdt=(a−0.01z)+(y2−0.7x2−3xy), $
|
(3.7) |
which has a single equilibrium at $ P:\; (0, 0,100a) $. According to Corollary 2, the equilibrium $ P $ is locally asymptotically stable for $ a < 0 $.
For the two initial conditions $ (\pm{0.8}, \pm{0.8}, 0) $, the bifurcation diagram is presented in Figure 10(a) for $ a\in [-0.03, 0] $. It shows the development of chaotic motions via period-doubling cascades. For the initial data $ (x(0), y(0), z(0)) = (0.8, 0.8, 0) $, the variation of Lyapunov exponents with respect to the parameter $ a $ is show in Figure 10(b). From Figure 10, there is a narrow range near $ a = 0 $ with $ a < 0 $ that corresponds to hidden chaotic dynamics with a stable node-focus.
For two initial conditions $ (\pm{0.8}, \pm{0.8}, 0) $, the bifurcation diagram is presented in Figure 11 for $ a\in [0, 2\times 10^{(-4)}] $. When $ a = 0 $, from (2.12), we know that the first Lyapunov number $ a_{1} $ of system (3.7) is positive, thus the origin is unstable. Recalling from Corollary 2.2 that the equilibrium $ P $ is a saddle-focus of index 2 for $ a\in (0, 2\times 10^{(-4)}] $, hence, we have that the chaotic range $ [0, 2\times 10^{(-4)}] $ corresponds to self-excited dynamics.
With two initial conditions $ (\mp{0.8}, \mp{0.8}, -\mu) $, the variations of the average values of the state variables in systems (2.14) and (2.15) with respect to $ \mu $ are shown in Figures 12 and 13, respectively. Since the difference between the two system is caused by a small constant of 0.001, the figures are almost the same. For both systems, the average value of $ z $ decreases smoothly as the parameter $ \mu $ increases.
Let $ (L_1^{(k)}, L_2^{(k)}, L_3^{(k)}), \; k = 1, 2 $, be the Lyapunov exponent spectra of systems (2.14) and (2.15), respectively. For each $ k $, we have $ L_1^{(k)} > 0 $, $ L_2^{(k)} = 0 $, and $ L_3^{(k)} < 0 $. For $ 1\leq k\leq 2 $, the variations of $ L_{1}^{(k)} $ and $ L_{3}^{(k)} $ with respect to the parameter $ \mu $ are shown in Figure 14.
Figure 14 indicates that as the initial conditions vary with the parameter $ \mu $, both of these systems exhibit chaotic dynamics. For system (2.14), there exists a symmetric pair of hidden chaotic attractors accompanied by an infinite number of equilibrium points. As for system (2.15), there also exists a symmetric pair of hidden chaotic attractors, but there are no equilibrium points. As shown in panel (a) of Figure 14, the largest Lyapunov exponent $ L_1^{(2)} $ (red curve) of system (2.15) consistently exceeds the largest Lyapunov exponent $ L_1^{(1)} $ (blue curve) of system (2.14), indicating that system (2.15) demonstrates stronger chaotic dynamics than system (2.14).
For systems (2.14) and (2.15) with $ \mu = -5, 5 $, the coexistence of hidden chaotic attractors is illustrated in Figures 15 and 16, respectively.
In this paper, we have conducted a comprehensive investigation of a novel chameleon system, examining its symmetry, dissipation, stability analysis, Hopf bifurcation, and various chaotic dynamics. The parameter $ \mu $ serves as both the Hopf bifurcation parameter and the offset boosting parameter. The other parameters are also critical control parameters, with variation of which, the system undergoes period-doubling bifurcations leading to various types of chaos, including hidden chaos and self-excited chaos, with hidden chaos being predominant. By offset boosting the variable $ z $, we also investigated the coexistence of hidden chaotic attractors. The results of this study enrich our understanding of the dynamical behavior of chameleon systems, providing important insights for research in related fields such as information security, cryptography, and nonlinear dynamics. In the future, it is hoped that there will be studies about chameleon systems subjected to small random perturbations.
Jie Liu: Writing-first draft. Bo Sang: Conceptualization, Writing-second and third drafts. Lihua Fan: Software. Chun Wang: Software. Xueqing Liu: Software. Ning Wang: Conceptualization, Software. Irfan Ahmad: Conceptualization. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research is supported by the Shandong Provincial Natural Science Foundation (No. ZR2018MA025).
The authors declare that they have no competing interests in this paper.
[1] |
Huttenhower C, Gevers D, Knight R, et al. (2012) Structure, function and diversity of the healthy human microbiome. Nature 486: 207–214. doi: 10.1038/nature11234
![]() |
[2] |
Backhed F (2011) Programming of host metabolism by the gut microbiota. Ann Nutr Metab 58: 44–52. doi: 10.1159/000328042
![]() |
[3] | Cani P, Everard A (2015) Talking microbes: when gut bacteria interact with diet and host organs. Mol Nutr Food Res 60: 58–66. |
[4] |
Cotillard A, Kennedy S, Kong L, et al. (2013) Dietary intervention impact on gut microbial gene richness. Nature 500: 585–588. doi: 10.1038/nature12480
![]() |
[5] | Wu G, Compher C, Chen E, et al. (2014) Comparative metabolomics in vegans and omnivores reveal constraints on diet-dependent gut microbiota metabolite production. Gut 65: 63–72. |
[6] | Lim M, Rho M, Song Y, et al. (2014) Stability of gut enterotypes in korean monozygotic twins and their association with biomarkers and diet. Sci Rep 4: 7348. |
[7] |
Org E, Parks B, Joo J, et al. (2015) Genetic and environmental control of host-gut microbiota interactions. Genome Res 25: 1558–1569. doi: 10.1101/gr.194118.115
![]() |
[8] |
Tilg H, Cani P, Mayer E, et al. (2016) Gut microbiome and liver diseases. Gut 65: 2035–2044. doi: 10.1136/gutjnl-2016-312729
![]() |
[9] | Tabibian J, Varghese C, LaRusso N, et al. (2015) The enteric microbiome in hepatobiliary health and disease. Liver Int 36: 480–487. |
[10] |
Lazaridis K, LaRusso N (2016) Primary sclerosing cholangitis. New Engl J Med 375: 2500–2502. doi: 10.1056/NEJMc1613273
![]() |
[11] |
Trivedi P, Adams D (2016) Gut-liver immunity. J Hepatol 64: 1187–1189. doi: 10.1016/j.jhep.2015.12.002
![]() |
[12] | Kummen M, Holm K, Anmarkrud J, et al. (2016) The gut microbial profile in patients with primary sclerosing cholangitis is distinct from patients with ulcerative colitis without biliary disease and healthy controls. Gut 66: 611–619. |
[13] |
Bajer L, Kverka M, Kostovcik M, et al. (2017) Distinct gut microbiota profiles in patients with primary sclerosing cholangitis and ulcerative colitis. World J Gastroentero 23: 4548. doi: 10.3748/wjg.v23.i25.4548
![]() |
[14] |
Sabino J, Vieira-Silva S, Machiels K, et al. (2016) Primary sclerosing cholangitis is characterised by intestinal dysbiosis independent from IBD. Gut 65: 1681–1689. doi: 10.1136/gutjnl-2015-311004
![]() |
[15] | Tabibian J, O-Hara S, Trussoni C, et al. (2015) Absence of the intestinal microbiota exacerbates hepatobiliary disease in a murine model of primary sclerosing cholangitis. Hepatology 63: 185–196. |
[16] |
Tabibian J, O'Hara S, Splinter P, et al. (2014) Cholangiocyte senescence by way of n-ras activation is a characteristic of primary sclerosing choalngitis. Hepatology 59: 2263–2275. doi: 10.1002/hep.26993
![]() |
[17] | Tabibian J, O'Hara S, Lindor K (2014) Primary sclerosing cholangitis and the microbiota: current knowledge and perspectives on etiopathogenesis and emerging therapies. Scand J Gastroentero 49: 901–908. |
[18] |
Eaton J, Juran B, Atkinson E, et al. (2015) A comprehensive assessment of environmental exposures among 1000 north american patients with primary sclerosing cholangitis, with and without inflammatory bowel disease. Aliment Pharm Therap 41: 980–990. doi: 10.1111/apt.13154
![]() |
[19] |
Lichtman S, Keku J, Schwab J, et al. (1991) Hepatic injury associated with small bowel bacterial overgrowth in rats is prevented by metronidazole and tetracycline. Gastroenterology 100: 513–519. doi: 10.1016/0016-5085(91)90224-9
![]() |
[20] | Tabibian J, Talwalkar J, Lindor K (2013) Role of the microbiota and antibiotics in primary sclerosing cholangitis. Biomed Res Int 2013: 1–7. |
[21] |
Tabibian J, Weeding E, Jorgensen R, et al. (2013) Randomised clinical trial: vancomycin or metronidazole in patients with primary sclerosing cholangitis-a pilot study. Aliment Pharm Therap 37: 604–612. doi: 10.1111/apt.12232
![]() |
[22] |
Mistilis S, Skyring A, Goulston S (1965) Effect of long-term tetracycline therapy, steroid therapy and colectomy in pericholangitis associated with ulcerative colitis. Australas Ann Med 14: 286–294. doi: 10.1111/imj.1965.14.4.286
![]() |
[23] |
Färkkilä M, Karvonen A, Nurmi H, et al. (2004) Metronidazole and ursodeoxycholic acid for primary sclerosing cholangitis: a randomized placebo-controlled trial. Hepatology 40: 1379–1386. doi: 10.1002/hep.20457
![]() |
[24] |
Silveira M, Torok N, Gossard A, et al. (2009) Minocycline in the treatment of patients with primary sclerosing cholangitis: results of a pilot study. Am J Gastroenterol 104: 83–88. doi: 10.1038/ajg.2008.14
![]() |
[25] |
Rankin J, Sydney M, Boden R, et al. (1959) The liver in ulcerative colitis; treatment of pericholangitis with tetracycline. Lancet 274: 1110–1112. doi: 10.1016/S0140-6736(59)90098-4
![]() |
[26] | Mathew K (1983) Metronidazole in primary cholangitis. J Indian Med Assoc 80: 31–33. |
[27] |
Kozaiwa K, Tajiri H, Sawada A, et al. (1998) Case report: three paediatric cases of primary sclerosing cholangitis treated with ursodeoxycholic acid and sulphasalazine. J Gastroen Hepatol 13: 825–829. doi: 10.1111/j.1440-1746.1998.tb00740.x
![]() |
[28] |
Cox KL, Cox KM (1998) Oral vancomycin: treatment of primary sclerosing cholangitis in children with inflammatory bowel disease. J Pediatr Gastr Nutr 27: 580–583. doi: 10.1097/00005176-199811000-00015
![]() |
[29] | Broccoletti T, Ciccimarra E, Spaziano M, et al. (2002) Refractory primary sclerosing cholangitis becoming responsive after sulphasalazine treatment of an underlying silent colitis. Ital J Pediatr 28: 515–517. |
[30] |
Tada S, Ebinuma H, Saito H, et al. (2006) Therapeutic benefit of sulfasalazine for patients with primary sclerosing cholangitis. J Gastroenterol 41: 388–389. doi: 10.1007/s00535-005-1758-x
![]() |
[31] |
Boner A, Peroni D, Bodini A, et al. (2007) Azithromycin may reduce cholestasis in primary sclerosing cholangitis: a case report and serendipitous observation. Int J Immunopath Ph 20: 847–849. doi: 10.1177/039463200702000423
![]() |
[32] |
Davies Y, Cox K, Abdullah B, et al. (2008) Long-term treatment of primary sclerosing cholangitis in children with oral vancomycin: an immunomodulating antibiotic. J Pediatr Gastr Nutr 47: 61–67. doi: 10.1097/MPG.0b013e31816fee95
![]() |
[33] |
Pohl J, Ring A, Stremmel W, et al. (2006) The role of dominant stenoses in bacterial infections of bile ducts in primary sclerosing cholangitis. Eur J Gastroen Hepat 18: 69–74. doi: 10.1097/00042737-200601000-00012
![]() |
[34] |
Hiramatsu K, Harada K, Tsuneyama K, et al. (2000) Amplification and sequence analysis of partial bacterial 16S ribosomal RNA gene in gallbladder bile from patients with primary biliary cirrhosis. J Hepatol 33: 9–18. doi: 10.1016/S0168-8278(00)80153-1
![]() |
[35] |
Olsson R, Bjornsson E, Backman L, et al. (1998) Bile duct bacterial isolates in primary sclerosing cholangitis: a study of explanted livers. J Hepatol 28: 426–432. doi: 10.1016/S0168-8278(98)80316-4
![]() |
[36] |
Sasatomi K, Noguchi K, Sakisaka S (1998) Abnormal accumulation of endotoxin in biliary epithelial cells in primary biliary cirrhosis and primary sclerosing cholangitis. J Hepatol 29: 409–416. doi: 10.1016/S0168-8278(98)80058-5
![]() |
[37] |
Vleggaar F, Monkelbaan J, Van Erpecum K (2008) Probiotics in primary sclerosing cholangitis: a randomized placebo-controlled crossover pilot study. Eur J Gastroen Hepat 20: 688–692. doi: 10.1097/MEG.0b013e3282f5197e
![]() |
[38] | Allegretti J, Kassam Z, Carrellas M, et al. (2017) Fecal microbiota transplantation improves microbiome diversity and liver enzyme profile in primary sclerosing cholangitis. World Congress of Gastroenterology at ACG2017 Meeting Abstracts. |
[39] | David L, Maurice C, Carmody R, et al (2013) Diet rapidly and reproducibly alters the human gut microbiome. Nature 505: 559–563. |
[40] |
Andersen I, Tengesdal G, Lie B, et al. (2014) Effects of coffee consumption, smoking, and hormones on risk for primary sclerosing cholangitis. Clin Gastroenterol H 12: 1019–1028. doi: 10.1016/j.cgh.2013.09.024
![]() |
[41] |
Lammert C, Juran B, Schlicht E, et al. (2014) Reduced coffee consumption among individuals with primary sclerosing cholangitis but not primary biliary cirrhosis. Clin Gastroenterol H 12: 1562–1568. doi: 10.1016/j.cgh.2013.12.036
![]() |
[42] | Torres J, Palmela C, Bao X, et al. (2016) The gut microbiota in primary sclerosing cholangitis and inflammatory bowel disease: correlations with disease and diet. Microbiology 150: S583–S584. |
[43] |
Boonstra K, Beuers U, Ponsioen C (2012) Epidemiology of primary sclerosing cholangitis and primary biliary cirrhosis: a systematic review. J Hepatol 56: 1181–1188. doi: 10.1016/j.jhep.2011.10.025
![]() |
[44] |
Dyson J, Hirschfield G, Adams D, et al. (2015) Novel therapeutic targets in primary biliary cirrhosis. Nat Rev Gastro Hepat 12: 147–158. doi: 10.1038/nrgastro.2015.12
![]() |
[45] |
Mattner J (2016) Impact of microbes on the pathogenesis of primary biliary cirrhosis (PBC) and primary sclerosing cholangitis (PSC). Int J Mol Sci 17: 1864. doi: 10.3390/ijms17111864
![]() |
[46] |
Kaplan M, Gershwin E (2005) Primary biliary cirrhosis. New Engl J Med 353: 1261–1273. doi: 10.1056/NEJMra043898
![]() |
[47] |
Selmi C, Bowlus C, Gershwin E, et al. (2011) Primary biliary cirrhosis. Lancet 377: 1600–1609. doi: 10.1016/S0140-6736(10)61965-4
![]() |
[48] |
Chassaing B, Etienne-Mesmin L, Bonnet R, et al. (2013) Bile salts induce long polar fimbriae expression favouring crohn's disease-associated adherent-invasive escherichia coli interaction with peyer's patches. Environ Microbiol 15: 355–371. doi: 10.1111/j.1462-2920.2012.02824.x
![]() |
[49] |
Kakiyama G, Pandak W, Gillevet P, et al. (2013) Modulation of the fecal bile acid profile by gut microbiota in cirrhosis. J Hepatol 58: 949–955. doi: 10.1016/j.jhep.2013.01.003
![]() |
[50] |
Lv L, Fang D, Shi D, et al. (2016) Alterations and correlations of the gut microbiome, metabolism and immunity in patients with primary biliary cirrhosis. Environ Microbiol 18: 2272–2286. doi: 10.1111/1462-2920.13401
![]() |
[51] |
Sasatomi K, Noguchi K, Sakisaka S, et al. (1998) Abnormal accumulation of endotoxin in biliary epithelial cells in primary biliary cirrhosis and primary sclerosing cholangitis. J Hepatol 29: 409–416. doi: 10.1016/S0168-8278(98)80058-5
![]() |
[52] |
Maldonado R, Sa-Correia I, Valvano M (2016) Lipopolysaccharide modification in gram-negative bacteria during chronic infection. FEMS Microbiol Rev 40: 480–493. doi: 10.1093/femsre/fuw007
![]() |
[53] |
Wang A, Migita K, Ito M, et al. (2005) Hepatic expression of toll-like receptor 4 in primary biliary cirrhosis. J Autoimmun 25: 85–91. doi: 10.1016/j.jaut.2005.05.003
![]() |
[54] |
Benias P, Gopal K, Bodenheimer H, et al. (2012) Hepatic expression of toll-like receptors 3, 4, and 9 in primary biliary cirrhosis and chronic hepatitis c. Clin Res Hepatol Gas 36: 448–454. doi: 10.1016/j.clinre.2012.07.001
![]() |
[55] |
Seki E, Brenner D (2008) Toll-like receptors and adaptor molecules in liver disease: update. Hepatology 48: 322–335. doi: 10.1002/hep.22306
![]() |
[56] |
Wahlström A, Sayin S, Marschall H, et al. (2016) Intestinal crosstalk between bile acids and microbiota and its impact on host metabolism. Cell Metab 24: 41–50. doi: 10.1016/j.cmet.2016.05.005
![]() |
[57] |
Visschers R, Luyer M, Schaap F, et al. (2013) The gut-liver axis. Curr Opin Clin Nutr 16: 576–581. doi: 10.1097/MCO.0b013e32836410a4
![]() |
[58] |
Lindor K, Gershwin E, Poupon R, et al. (2009) Primary biliary cirrhosis. Hepatology 50: 291–308. doi: 10.1002/hep.22906
![]() |
[59] |
Beuers U, Boberg K, Chapman R (2009) EASL clinical practice guidelines: management of cholestatic liver diseases. J Hepatol 51: 237–267. doi: 10.1016/j.jhep.2009.04.009
![]() |
[60] |
Lee J, Arai H, Nakamura Y, et al. (2013) Contribution of the 7β-hydroxysteroid dehydrogenase from ruminococcus gnavus N53 to ursodeoxycholic acid formation in the human colon. J Lipid Res 54: 3062–3069. doi: 10.1194/jlr.M039834
![]() |
[61] | Tang R, Wei Y, Che W, et al. (2017) Gut microbial profile is altered in primary biliary cholangitis and partially restored after UDCA therapy. Gut 67: 534–541. |
[62] |
Hopf U, Moller B, Stemerowicz R, et al. (1989) Relation between Escherichia coli R(rough)-forms in gut, lipid A in liver, and primary biliary cirrhosis. Lancet 334: 1419–1422. doi: 10.1016/S0140-6736(89)92034-5
![]() |
[63] | Bogdanos D, Baum H, Grasso A, et al. (2004) Microbial mimics are major targets of crossreactivity with human pyruvate dehydrogenase in primary biliary cirrhosis. J Hepatol 40: 31–39. |
[64] | Olafsson S, Gudjonsson H, Selmi C, et al. (2004) Antimitochondrial antibodies and reactivity to N. aromaticivorans proteins in icelandic patients with primary biliary cirrhosis and their relatives. Am J Gastroenterol 99: 2143–2146. |
[65] |
Mattner J, Savage P, Leung P, et al. (2008) Liver autoimmunity triggered by microbial activation of natural killer T cells. Cell Host Microbe 3: 304–315. doi: 10.1016/j.chom.2008.03.009
![]() |
[66] |
Smyk D, Bogdanos D, Kriese S, et al. (2012) Urinary tract infection as a risk factor for autoimmune liver disease: from bench to bedside. Clin Res Hepatol Gas 36: 110–121. doi: 10.1016/j.clinre.2011.07.013
![]() |
[67] | Rashid T, Ebringer A (2012) Autoimmunity in rheumatic diseases is induced by microbial infections via crossreactivity or molecular mimicry. Autoimmune Dis: 1–9. |
[68] |
Stemerowicz R, Hopf U, Moller B, et al. (1988) Are antimitochondrial antibodies in primary biliary cirrhosis induced by R(rough)-mutants of enterobacteriaceae? Lancet 332: 1166–1170. doi: 10.1016/S0140-6736(88)90235-8
![]() |
[69] |
Gossard A, Lindor K (2012) Autoimmune hepatitis: a review. J Gastroenterol 47: 498–503. doi: 10.1007/s00535-012-0586-z
![]() |
[70] |
Boberg K, Aadland E, Jahnse J, et al. (1998) Incidence and prevalence of brimary biliary cirrhosis, primary sclerosing cholangitis, and autoimmune hepatitis in a norwegian population. Scand J Gastroentero 33: 99–103. doi: 10.1080/00365529850166284
![]() |
[71] |
Lee Y, Teo E, Ng T, et al. (2001) Autoimmune hepatitis in singapore: a rare syndrome affecting middle-aged women. J Gastroen Hepatol 16: 1384–1389. doi: 10.1046/j.1440-1746.2001.02646.x
![]() |
[72] |
Hurlburt K, McMahon B, Deubner H, et al. (2002) Prevalence of autoimmune liver disease in Alaska natives. Am J Gastroenterol 97: 2402–2407. doi: 10.1111/j.1572-0241.2002.06019.x
![]() |
[73] |
Primo J, Merino C, Fernandez J, et al. (2004) Incidence and prevalence of autoimmune hepatitis in the area of the Hospital de Sagunto (spain). Gastroent Hepat-Barc 27: 239–243. doi: 10.1016/S0210-5705(03)70452-X
![]() |
[74] |
Ngu J, Bechly K, Chapman B, et al. (2010) Population-based epidemiology study of autoimmune hepatitis: a disease of older women? J Gastroen Hepatol 25: 1681–1686. doi: 10.1111/j.1440-1746.2010.06384.x
![]() |
[75] | Delgado J, Vodonos A, Malnick S, et al. (2013) Autoimmune hepatitis in southern Israel: a multicenter study of 15 years. J Digest Dis 14: 611–618. |
[76] |
Czaja A (2016) Factoring the intestinal microbiome into the pathogenesis of autoimmune hepatitis. World J Gastroentero 22: 9257–9278. doi: 10.3748/wjg.v22.i42.9257
![]() |
[77] |
Czaja A (1998) Frequency and nature of the variant syndromes of autoimmune liver disease. Hepatology 28: 360–365. doi: 10.1002/hep.510280210
![]() |
[78] | Czaja A (2013) The overlap syndromes of autoimmune hepatitis. Digest Dis Sci 58: 326–343. |
[79] | Lin R, Zhou L, Zhang J, et al. (2015) Abnormal intestinal permeability and microbiota in patients with autoimmune hepatitis. Int J Clin Exp Patho 8: 5153–5160. |
[80] |
Cai W, Ran Y, Li Y, et al. (2017) Intestinal microbiome and permeability in patients with autoimmune hepatitis. Best Pract Res Cl Ga 31: 669–673. doi: 10.1016/j.bpg.2017.09.013
![]() |
[81] |
Yuksel M, Wang Y, Tai N, et al. (2015) A novel humanized mouse model for autoimmune hepatitis and the association of gut microbiota with liver inflammation. Hepatology 62: 1536–1550. doi: 10.1002/hep.27998
![]() |
[82] |
Reumaux D, Duthilleul P, Roos D (2004) Pathogenesis of diseases associated with antineutrophil cytoplasm autoantibodies. Hum Immunol 65: 1–12. doi: 10.1016/j.humimm.2003.09.013
![]() |
[83] |
Terjung B, Sohne J, Lechtenberg B, et al. (2010) p-ANCAs in autoimmune liver disorders recognise human beta-tubulin isotype 5 and cross-react with microbial protein FtsZ. Gut 59: 808–816. doi: 10.1136/gut.2008.157818
![]() |