Bifurcation analysis of a two–dimensional p53 gene regulatory network without and with time delay

1.   Introduction

Neural networks (NNs) are computational models that represent the information processing aspect of the human brain's neural network. NNs have been the focus of extensive studies across various domains in recent decades, including but not limited to automatic control ^[1], transmission of biological information ^[2], medical image analysis ^[3], system identification ^[4], and associate memory ^[5]. Synchronization is an essential tool for secure communication and biocomputing. It is the process of maintaining consistency in a system's dynamic signals as time passes. Synchronization, a common collective behavior, has emerged as a significant concern in the field of NNs. Researchers have developed various synchronization approaches, including pinning synchronization ^[6], finite time synchronization ^[7], and exponential synchronization ^[8]. In contrast to solitary NNs, coupled NNs (CNNs) exhibit more intricate and imperceptible behaviors. In recent years, researchers have made considerable progress in studying synchronization of CNNs ^[9,10,11,12].

Due to the limitations caused by the speed at which signals may travel, specifically the speed at which neuron amplifiers in NNs can switch and transmit, the system trajectory is contingent not only on the previous state but also on the current one. The delay category, also known as neutral-type delay, has the potential to induce instability or other undesirable dynamic behaviors ^[13,14,15]. Furthermore, the time delays observed in CNNs may be subject to random occurrences due to the temporal signals and synaptic voltage fluctuations transmitted by the transmitters. This can result in extremely large time delay values, despite the extremely low probability that they will occur. A degree of conservatism results if information regarding the range of time delay variations is the only factor considered, without taking their probability into consideration. Thus, a succession of research findings concerning CNNs with probabilistic time-varying latencies are available ^[16,17,18]. In addition to time delay effects, real NNs are susceptible to impulsive effects when system states undergo instantaneous disturbances or abrupt state changes at specific instants. These effects can have a similar impact on systems' dynamical behaviors as time-delay effects do. Thus far, considerable interest has been devoted to the synchronization or stability of the coupled neutral NNs (CNNNs) ^[19,20,21] and impulsive coupled NNs (ICNNs) ^[22,23].

Nerve signals are conveyed through chaotic electrical pulses in biological nervous systems, which are susceptible to stochastic disturbances and random noises. Frequently used to model these stochastic factors are Gaussian white noise or Brownian motion, considered reasonable approximations. As of now, a multitude of dynamic behaviors exhibited by coupled stochastic NNs, such as synchronization and stability, have been examined through the modeling of randomness using Brownian motion ^{[14,16,20,24,25]}. Nevertheless, rare instances are not uncommon among genuine biological neurons. Neuronal impulses discharge and chemical processes in neuronal synaptic receivers generate jump-type noise ^[26]. Therefore, the Lévy process, which expands Brownian motion to include jump-diffusion, is a more suitable model for these circumstances compared to Brownian motion ^{[27,28,29,30]}. Notwithstanding the augmented mathematical intricacy, the findings of this study suggest that examining synchronization in impulsive coupled neutral stochastic NNs (ICNSNNs) propelled by Lévy noise is more significant.

The majority of the previously published findings on the examination of NNs synchronization necessitate the activation function to adhere to Lipschitz continuity. Nevertheless, in actual situations, the Lipschitz conditions frequently impose excessively stringent requirements, with certain conditions proving exceedingly challenging to fulfill ^[31]. There has been considerable scholarly activity aimed at easing restrictions on activation functions, including but not limited to one-sided Lipschitz conditions and local Lipschitz conditions ^[32,33,34]. Accordingly, it is critical, in light of these findings, to derive lower-limit synchronization conditions for ICNSNNs with Lévy noise.

To attain synchronization among driver-response systems, it is customary to devise control protocols that specifically target the desired control outcome. CNNs are capable of implementing a wide range of control schemes, including adaptive control ^[13,20], impulsive control ^[30,35], and event-triggered control ^[27,28]. It is widely acknowledged that adaptive controllers possess numerous advantages and find utility across diverse domains. Adaptive control is appealing and intriguing due to its robustness and ability to adjust autonomously in response to various updating laws, suitable for systems characterized by strong nonlinearity and discontinuity on the righthand side.

Based on the preceding discussion, we consider the problem of exponential synchronization in ICNSNNs with Lévy noise, specifically in the presence of non-Lipschitz conditions. The subsequent text provides a concise overview of the main advancements in contrast to the existing body of research.

1) Compared to earlier neutral stochastic NN models, the model in this study incorporates impulse effects and probabilistic delays into the coupled neutral stochastic NNs, making it more versatile and beneficial in engineering practice.

2) The requirements necessary for synchronization analysis, which do not need Lipschitz conditions, have been created, therefore relaxing the limitations imposed by Lipschitz conditions.

3) Due to the lack of satisfaction of the chain rule by the Itô-type stochastic integral, the Dini differential method ^{[16,22,23,35]} presents challenges in solving the proof difficulties arising from the combination of impulsive component and neutral delay. Our technique is resistant to the aforementioned limitations and can efficiently tackle these challenges.

4) Under the stochastic perturbations, impulses, and various delays, by resorting to an adaptive controller, some sufficient conditions have been established to make ICNSNNs with Lévy noise achieve exponential synchronization.

Notations: Please refer to the following Table 1 for specific symbols.

2.   Models and preliminaries

We consider the neutral NNs with probabilistic delays of the form

where the variable $\tilde{x}(t)\in \mathbb{R}^{n}$ reflects the state of neutral NNs that are connected to $n$ neurons, and $\tilde{H}(\tilde{x}(t))\in \mathbb{R}^{n}$ is the activation function of the neurons. $\mathfrak{G}(t)$ is time-varying delay and satisfies conditions $\tilde{\mathfrak{G}} \geq\mathfrak{G}(t)\geq 0$ and $1 > \bar{\mathfrak{G}}\geq\dot{\mathfrak{G}}(t)$. $A$ is the neutral parameter matrix. $B = diag(b_{1}, b_{2}, \ldots, b_{n}) > 0$. $\tilde{E}\in \mathbb{R}^{n}$. The connection weight matrix is denoted as $C$, whereas the delay connection weight matrix is denoted as $D$. The delay $\mathfrak{P}(t)$ in system (2.1) is bounded and meets condition $0\leq\mathfrak{P}(t)\leq \tilde{\mathfrak{P}}_{2}$. Practically, there is a constant $\tilde{\mathfrak{P}}_{1}$ that satisfies $0\leq\tilde{\mathfrak{P}}_{1}\leq \tilde{\mathfrak{P}}_{2}$. Moreover, $\mathfrak{P}(t)$ assumes values within the intervals $[0, \tilde{\mathfrak{P}}_{1}]$ and $(\tilde{\mathfrak{P}}_{1}, \tilde{\mathfrak{P}}_{2}]$ with a specific probability, as determined by its probability distribution, i.e., $P\{\mathfrak{P}(t)\in[0, \tilde{\mathfrak{P}}_{1}]\} = \mathfrak{K}_{0}$ and $P\{\mathfrak{P}(t)\in(\tilde{\mathfrak{P}}_{1}, \tilde{\mathfrak{P}}_{2}]\} = 1-\mathfrak{K}_{0}$, $0\leq \mathfrak{K}_{0}\leq 1$.

Then

where $\mathfrak{P}_{k}(t)$ satisfies $1 > \bar{\mathfrak{P}_{k}}\geq\dot{\mathfrak{P}_{k}}(t), k = 1, 2.$ The stochastic variable $\mathfrak{K}(t)$ follows a Bernoulli distribution, denoted by $P\{\mathfrak{K}(t) = 1\} = P\{\mathfrak{P}(t)\in [0, \tilde{\mathfrak{P}}_{1}]\}$ and $P\{\mathfrak{K}(t) = 0\} = P\{\mathfrak{P}(t)\in (\tilde{\mathfrak{P}}_{1}, \tilde{\mathfrak{P}}_{2}]\}$. Then

By utilizing the new functions $\mathfrak{P}_{1}(t)$ and $\mathfrak{P}_{2}(t)$, together with the stochastic variables $\mathfrak{K}(t)$, system (2.1) may be reformulated as

The system (2.2) is considered the drive system, with the state variable represented as $\tilde{\chi}(t)$. The response system, on the other hand, is described by the following ICNSNN, with the state variable denoted as $\tilde{y}_{k}(t)$. $e_{k}(t) = \tilde{y}_{k}(t)-\tilde{\chi}(t)$ specifies the following error vector:

In system (2.3), the control input vector, denoted as $\Psi_{k} = (\varphi_{k1}, \varphi_{k2}, \ldots, \varphi_{kn})^{T} \in \mathbb{R}^{n}$, $k = 1, ..., N$, represents the control inputs for a system. $\tilde{\varrho} > 0$ is coupling strength. The configuration matrix $O = \{o_{kj}\}_{N\times N}$ denotes the topological structure of systems. A connection from node $k$ to $j$ is defined as $o_{kj} > 0$, otherwise $o_{kj} = 0$. Additionally, $o_{kk} = -\sum\limits_{j = 1, j\neq k}^{N}o_{kj}$. The positive definite diagonal matrix $\Upsilon$ represents the internal connection strength between two interconnected NNs. $e_{k}(\sigma_{l}) = e_{k}(\sigma_{l}^{+}) = \lim\limits_{t\rightarrow \sigma_{l}^{+}}e_{k}(t)$ and $e_{k}(\sigma_{l}^{-}) = \lim\limits_{t\rightarrow \sigma_{l}^{-}}e_{k}(t)$. $\tilde{J}_{l}(\cdot, \cdot)\in \mathbb{R}^{n}$. The impulsive time instants, denoted as $\sigma_{l}$, follow the conditions $0 = \sigma_{0} < \sigma_{1} < \sigma_{2} < \cdots < \sigma_{l} < \cdots$, and $\lim\limits_{k\rightarrow \infty}\sigma_{l} = \infty$. $\mu: \mathbb{R}^{n}\times \mathbb{R}^{n}\times \mathbb{R}^{n}\times \mathbb{R}^{n}\rightarrow \mathbb{R}^{n \times n}$ and $q: \mathbb{R}^{n}\times \mathbb{R}^{n}\times \mathbb{R}^{n}\times \mathbb{R}^{n}\times \mathbb{V} \rightarrow \mathbb{R}^{n}$ indicate Lévy noise intensity functions. A vector Brownian motion is denoted by $\omega(t)$, which is defined on the probability space $(\Omega, \mathcal{F}, \{\mathcal{F}_{t}\}_{t\geq0}, P)$, where $\{\mathcal{F}_{t}\}_{t\geq0}$ meets the standard conditions. $N(dt, d\nu)$ is a Poisson counting measure utilizing a characteristic measure $\pi$ on a measurable subset $\mathbb{V}$ of $\mathbb{R}$, whereas $N(t, \nu)$ represents a Poisson process. $\tilde{N}(dt, d\nu) = N(dt, d\nu)-\pi(d\nu)dt$. The research operates under the assumption that stochastic processes $\omega$ and $N$ are independent.

According to systems (2.2) and (2.3) and the quality of matrix $O$,

the error dynamics system can be obtained such that

in which $H(e_{k}(t)) = \tilde{H}(\tilde{y}_{k}(t))-\tilde{H}(\tilde{\chi}(t))$, $H(e_{k}(t-\mathfrak{P}_{1}(t))) = \tilde{H}(\tilde{y}_{k}(t-\mathfrak{P}_{1}(t)))-\tilde{H}(\tilde{\chi}(t-\mathfrak{P}_{1}(t)))$, and $H(e_{k}(t-\mathfrak{P}_{2}(t))) = \tilde{H}(\tilde{y}_{k}(t-\mathfrak{P}_{2}(t)))-\tilde{H}(\tilde{\chi}(t-\mathfrak{P}_{2}(t)))$, where $\tilde{\mathfrak{G}}_{0} = \max(\tilde{\mathfrak{G}}, \tilde{\mathfrak{P}}_{2})$. To establish synchronization, we assume control is defined as

where $\aleph$ represents a nonnegative real number. The adaptive feedback gain is denoted as $\varphi_{k}(t) = diag(\varphi_{k1}(t), \varphi_{k2}(t), ..., \varphi_{kn}(t))(k = 1, ..., N)$.

To determine the synchronization requirements for systems (2.2) and (2.3), we provide the following assumptions.

$\textbf{A}_{1}$ There is a constant $\varrho > 0$ that fulfills

for $\forall \theta_{1}, \theta_{2} \in \mathbb{R}^{n}$, $\tilde{H}(0) = 0$, where $\varpi(\cdot): \mathbb{R}^{+}\rightarrow \mathbb{R}^{+}$ is a concave increasing continuous function with $\varpi(0) = 0$ and $\int_{0}^{1}\frac{du}{\varpi(u)} = \infty$. Additionally, it is supposed that there exist a positive constant $\varpi_{0}$ and a nonnegative function $\kappa(t)$ with finite upper bound, such that

for $\forall x\geq 0$. Here, $\kappa(t)$ meets $\int_{0}^{t}\exp(\vartheta s)\kappa(s)ds < \infty$ for $\forall\vartheta > 0$.

$\textbf{A}_{2}$ There exist constants $\Lambda_{1} > 0$, $\Lambda_{2} > 0$, $\Lambda_{3} > 0$, $\Gamma_{1} > 0$, $\Gamma_{2} > 0$, and $\Gamma_{3} > 0$ such that

and

hold for $\forall \theta_{1}, \theta_{2}, \theta_{3}, \theta_{4} \in \mathbb{R}^{n}$. In addition, $\mu(0, 0, 0) = 0$ and $q(0, 0, 0, 0) = 0$.

$\textbf{A}_{3}$ There exists a constant $\varsigma \in (0, 1)$, which ensures that $\rho (A)$ fulfills $\rho (A)\leq \varsigma$.

$\textbf{A}_{4}$ The function $\tilde{J}_l(\cdot, \cdot)$ is assumed to meet the specified conditions.

where $\Xi_{l} > 0$.

Remark 2.1. If $\varpi(x) = x$, then the Lipschitz criteria are satisfied, indicating that the assumption $\textbf{A}_{1}$ is less stringent than the prerequisite of the Lipschitz continuity. Certain activation functions do not meet the criteria for Lipschitz continuity, yet they do exhibit non-Lipschitz continuity in practical applications. For example, if we are examining a function $\tilde{H}(x) = x\sin x$ that fails to meet the criteria of Lipschitz continuity, we may construct the concave nondecreasing function

for small enough $\epsilon \in (0, 1)$^[36], so that $\tilde{H}(x)$ encounters the assumption $\textbf{A}_{1}$.

According to the literature ^[36,37,38], for any initial value $\{e_{k}(\epsilon): -\tilde{\mathfrak{G}}_{0} \leq \epsilon \leq 0\} = \phi^{0}_{k}\in \mathbb{L}^{2}_{\mathcal{F}_{0}}([-\tilde{\mathfrak{G}}_{0}, 0]; \mathbb{R}^{n})$, the assumptions $\textbf{A}_{1}$-$\textbf{A}_{4}$ are believed to demonstrate that the solution symbolized by $e_{k}(t, \phi^{0}_{k})$ on $t \geq 0$ to the system (2.4) is both existent and unique. Evidently, the system (2.4) allows for $e_{k}(t, 0) = 0$. The expressions $e_{k}(t, \phi^{0}_{k})$, $e_{k}(t-\mathfrak{G}(t))$, $e_{k}(t-\mathfrak{P}_{1}(t))$, and $e_{k}(t-\mathfrak{P}_{2}(t))$ will be referred to as $e_{k}(t)$, $e_{k\mathfrak{G}(t)}$, $e_{k\mathfrak{P}_{1}(t)}$, and $e_{k\mathfrak{P}_{2}(t)}$, respectively, for simplicity.

Lemma 2.1. ^[30,39] Consider the function $\delta(t)$, which is piecewise continuous and nonnegative. It fulfills the condition

where $t_{i}$ are the points of discontinuity of the first type for $\delta(t)$. Additionally, we have $\tilde{\lambda}_{1}\geq0$, $\tilde{\lambda}_{2}\geq0$, $\tilde{\lambda}_{3}\geq0$. Then, for $t_{i} < t\leq t_{i+1}$,

Here, the notation $i(t_{0}, t)$ signifies the number of points $t_{i}$ that lie inside the interval $[t_{0}, t)$.

Definition 2.1. ^[20] The drive system (2.2) and the response system (2.3) are considered to be exponentially synchronized in mean square (ESMS) if the error system (2.4) is exponentially stable in mean square, meaning that for all $\phi^{0}_{k} \in \mathbb{L}^{2}_{\mathcal{F}_{0}}([-\tilde{\mathfrak{G}}_{0}, 0]; \mathbb{R}^{n})$,

3.   Exponential synchronization analysis

This section will provide the general requirements for the ESMS in systems (2.2) and (2.3) based on the assumptions $\textbf{A}_{1}$–$\textbf{A}_{4}$ mentioned above.

Theorem 3.1. Assume that conditions $\textbf{A}_{1}$–$\textbf{A}_{4}$ are satisfied. If

where

then when the adaptive controller (2.5) with update rule

is applied, where $\upsilon_{kj} > 0, k = 1, ..., N, j = 1, 2, \ldots, n,$ the systems (2.3) and (2.2) are ESMS.

Proof. Define

By adopting the general It$\hat{o}$'s formula ^[40] for the interval $t\in (\sigma_{l}, \sigma_{l+1})$ and doing integration on both sides from $\sigma_{l}$ to $t$, we derive

where

Through employing Young's inequality, one may infer that

Based on the given assumptions $\textbf{A}_{1}$, we may conclude that

and

Comparable to (3.3), it is possible to get

Besides, we have

Following the assumption $A_{2}$, there exist

and

When $t = \sigma_{l}$, it produces

where $\tilde{\Xi} = \max\limits_{l\in \mathbb{N}_{+}}\{1+\Xi_{l}\}$. By replacing (3.3)–(3.9) into (3.2) and subsequently calculating the mathematical expectation for both sides of (3.2), it yields

Therefore, for any values of $t$ inside the interval $[0, \sigma_{l+1})$, it follows that

Since

and

it can be deduced that

In addition, employing the inequality in reference ^[36], it can be yielded that

and we have

where $\hbar > 0$ is sufficiently large for $\tilde{\hbar} = \varsigma^{2}(1+\hbar)\hbar^{-1}\exp(\vartheta \tilde{\mathfrak{G}}) < 1$. Then,

Thus, Lemma 2.1 subsequently results in

Based on the inequality $m\leq\inf\limits_{l\in \mathbb{N}}(\sigma_{l+1}-\sigma_{l})$, we ultimately get

With respect to criterion (3.1), it may be concluded that systems (2.2) and (2.3) are ESMS.

Remark 3.1. Based on our comprehension, certain current CNNN and neutral NN models ^{[13,14,15,19,20]} do not take into account the impact of impulses. This research explores a wider model by incorporating impulsive effects. The Dini differential technique^{[16,22,23,35]} is not applicable to our primary problem since it involves a combination of neutral delay and impulsive factors, and Itô-type stochastic integrals do not adhere to the chain rule. The proof strategy outlined in Theorem 3.1 offers a method for resolving the difficulty.

Remark 3.2. Activation functions are required to provide Lipschitz continuity in several studies on the synchronization or stability of CNNs, CNNNs, and neutral NNs, as indicated by the references ^{[13,20,22,27,35]}. Theorem 3.1 provides sufficient criteria for achieving synchronization in the drive-response system without the need for the Lipschitz situation, therefore relaxing the limitations imposed by Lipschitz conditions. As a result, this study extends the current findings (see ^{[13,20,22,27,35]}).

When Lévy jump $q(e_{k}(t), e_{k}(t-\mathfrak{G}(t)), e_{k}(t-\mathfrak{P}_{1}(t)), e_{k}(t-\mathfrak{P}_{2}(t)), \nu) = 0$, the coupled error system (2.4) reduced to the following form:

Under these circumstances, we ascertain the subsequent outcomes.

Corollary 3.1. If conditions $\textbf{A}_{1}$–$\textbf{A}_{4}$ are fulfilled, for

where

then when the adaptive controller (2.5) with update rule

is applied, where $\upsilon_{kj} > 0, \; k = 1, ..., N, j = 1, 2, \ldots, n$, the systems (2.3) and (2.2) are ESMS.

Proof. Define

Employing Theorem 3.1 readily proves the result of Corollary 3.2, hence the proof is omitted in this context.

4.   Numerical simulations

This section presents numerical simulations to demonstrate the efficacy of our suggested strategy.

Example 4.1. Consider the following two-dimensional neutral NNs with probabilistic delays, where

$\tilde{\chi}(t) = (\tilde{\chi}_{1}(t), \tilde{\chi}_{2}(t))^{T}$, $P\{\mathfrak{K}(t) = 0\} = 0.2$, $E = [0.1, 0.1]^{T}$. The initial values are $\tilde{\chi}(t) = [2.3, 2.6]^{T}$. The parameters of the corresponding response system could be expressed as follows.

$\tilde{y}_{k}(t) = (\tilde{y}_{k1}(t), \tilde{y}_{k2}(t))^{T}$. $\tilde{y}_{1}(t) = [-1.2, -1.4]^{T}$, $\tilde{y}_{2}(t) = [-1.5, -1.8]^{T}$, $\tilde{y}_{3}(t) = [-2.2, -2.4]^{T}$. $\varphi_{kj}(t) = 0, k = 1, 2, 3, j = 1, 2.$ $\varrho = 0.5$. The approaches provided in references ^{[13,14,15,19,20]} are not applicable to this particular instance. Verifying that the parameters in Example 4.1 fulfill the constraints of Theorem 3.1 is straightforward. Consequently, Theorem 3.1 may be utilized to establish that the response system (2.3) should be considered ESMS in theoretical analysis, together with the driving system (2.2). The simulated findings depicted in Figures 1–4 unequivocally validate our conclusion. We can find that the error systems in Figure 3 (a)–(b) converge to zero as time progresses. Figure 4 depicts the evolutionary trajectory of the adaptive feedback gain $\varphi_{kj}(t) (k = 1, 2, 3, j = 1, 2)$.

Example 4.2. Consider the following two-dimensional neutral NNs with probabilistic delays, where

$\tilde{\chi}(t) = (\tilde{\chi}_{1}(t), \tilde{\chi}_{2}(t))^{T}$, $P\{\mathfrak{K}(t) = 0\} = 0.2$, $E = [0.1, 0.1]^{T}$. The initial values are $\tilde{\chi}(t) = [3.2, 3.4]^{T}$. The parameters of the corresponding response system can be given as

$k = 1, 2, 3.$ $\tilde{y}_{k}(t) = (\tilde{y}_{k1}(t), \tilde{y}_{k2}(t))^{T}$. $\tilde{y}_{1}(t) = [2.2, 2.4]^{T}$, $\tilde{y}_{2}(t) = [-1.3, -1.1]^{T}$, $\tilde{y}_{3}(t) = [-2.6, -2.8]^{T}$. $\varphi_{kj}(t) = 0, k = 1, 2, 3, j = 1, 2.$ $\varrho = 0.5$. By Remark 2.1, the activation function $\tilde{H}(\cdot)$ in this instance does not meet the Lipschitz continuity. When $\epsilon = 0.001$, $\varpi_{0} = 13$, and $\kappa(t) = \exp(-2\vartheta t)$ are used, then the parameters in this example fulfill the circumstances of Theorem 3.1, implying that the systems (2.2) and (2.3) should be ESMS in the theoretical analysis. Figures 5–8 display the results of the simulation, which amply supports our conclusion. Figure 6 shows the trajectories of the system (2.2) and the system (2.3) among them, and it is evident that as time goes on, their trajectories become consistent. The error's trajectory converges to zero, as illustrated in Figure 7. Figure 8 illustrates the evolutionary trajectory of the feedback gain $\varphi_{kj}(t)(k = 1, 2, 3, j = 1, 2)$.

5.   Conclusions

Through adaptive controller, we address the problem of exponential synchronization in ICNSNNs with Lévy noise and probabilistic time delays in cases where the conditions are non-Lipschitz. We additionally discover sufficient criteria for achieving ESMS. In this research, we enhance the practicality of the examined model by incorporating impulse and probabilistic time-varying delays in coupled neutral stochastic NNs, as compared to previous findings^{[13,14,15,19,20]}. The activation function's Lipschitz continuity is not necessary for our investigation, suggesting that our criteria acquired are less limiting than existing ones and can be used in a broader range of situations. Furthermore, the findings of this study propose an approach to the challenge of analyzing synchronization, which arises from the presence of both a neutral delay term and an impulse term. Finally, we present two numerical illustrations that confirm the theoretical discoveries.

In the non-Lipschitz case, developing the discriminant criteria for the exponential synchronization of ICNSNNs with Lévy noise is still a challenging issue. In addition, more discussion is required on the relaxation of the time delay condition. Lately, there has been a significant focus on the dynamic characteristics of systems that utilize event-triggered control and sliding-mode control. This is evident in the relevant literature ^{[41,42,43,44,45]}. Such issues deserve extra investigation and analysis.

Author contributions

S. M.: Conceptualization; S. M. and JM. L.: Writing-original draft; S. M. and Q. L.: Writing-review & editing; JM. L. and RN. L.: Visualization. All authors have read and agreed to the published version of the manuscript.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The work was supported in part by the Fundamental Research Funds for the Central Universities, North Minzu University under Grant No.2020KYQD17, and in part by the National Natural Science Foundation of China under Grant No. 12202005, and in part by the High-Level Talent Research Foundation of Anhui Agricultural University under Grant No. RC382106, and in part by the Philosophy and Social Science Foundation of Universities of Anhui Province under Grant No. 2023AH050970, and in part by the Nature Science Foundation of Jiangsu Province under Grant No. BK20220233, and in part by the Ningxia Outstanding Talent Support Program.

Conflict of interest

The authors declare no conflict of interest.