Loading [MathJax]/jax/output/SVG/jax.js
Research article

Potential plant benefits of endophytic microorganisms associated with halophyte Glycyrrhiza glabra L.

  • Received: 11 April 2024 Revised: 14 August 2024 Accepted: 11 September 2024 Published: 30 September 2024
  • In this study, bacteria associated with licorice (Glycyrrhiza glabra L.) were characterized through 16S rRNA gene analysis. Profiling of endophytic bacteria isolated from Glycyrrhiza glabra tissues revealed 18 isolates across the following genera: Enterobacter (4), Pantoea (3), Bacillus (2), Paenibacillus (2), Achromobacter (2), Pseudomonas (1), Escherichia (1), Klebsiella (1), Citrobacter (1), and Kosakonia (1). Furthermore, the beneficial features of bacterial isolates for plants were determined. The bacterial isolates showed the capacity to produce siderophores, hydrogen cyanide (HCN), indole-3-acetic acid (IAA), chitinase, protease, glucanase, lipase, and other enzymes. Seven bacterial isolates showed antagonistic activity against F. culmorum, F. solani, and R. solani. According to these results, licorice with antimicrobial properties may serve as a source for the selection of microorganisms that have antagonistic activity against plant fungal pathogens and may be considered potential candidates for the control of plant pathogens. The selected bacterial isolates, P. polymyxa GU1, A. xylosoxidans GU6, P. azotoformans GU7, and P. agglomerans GU18, increased root and shoot growth of licorice and were able to colonize the plant root. They can also serve as an active part of bioinoculants, improving plant growth.

    Citation: Gulsanam Mardonova, Vyacheslav Shurigin, Farkhod Eshboev, Dilfuza Egamberdieva. Potential plant benefits of endophytic microorganisms associated with halophyte Glycyrrhiza glabra L.[J]. AIMS Microbiology, 2024, 10(4): 859-879. doi: 10.3934/microbiol.2024037

    Related Papers:

    [1] Yan’e Wang , Zhiguo Wang, Chengxia Lei . Asymptotic profile of endemic equilibrium to a diffusive epidemic model with saturated incidence rate. Mathematical Biosciences and Engineering, 2019, 16(5): 3885-3913. doi: 10.3934/mbe.2019192
    [2] Yoichi Enatsu, Yukihiko Nakata . Stability and bifurcation analysis of epidemic models with saturated incidence rates: An application to a nonmonotone incidence rate. Mathematical Biosciences and Engineering, 2014, 11(4): 785-805. doi: 10.3934/mbe.2014.11.785
    [3] Pan Yang, Jianwen Feng, Xinchu Fu . Cluster collective behaviors via feedback pinning control induced by epidemic spread in a patchy population with dispersal. Mathematical Biosciences and Engineering, 2020, 17(5): 4718-4746. doi: 10.3934/mbe.2020259
    [4] Jummy F. David, Sarafa A. Iyaniwura, Michael J. Ward, Fred Brauer . A novel approach to modelling the spatial spread of airborne diseases: an epidemic model with indirect transmission. Mathematical Biosciences and Engineering, 2020, 17(4): 3294-3328. doi: 10.3934/mbe.2020188
    [5] Chang-Yuan Cheng, Shyan-Shiou Chen, Xingfu Zou . On an age structured population model with density-dependent dispersals between two patches. Mathematical Biosciences and Engineering, 2019, 16(5): 4976-4998. doi: 10.3934/mbe.2019251
    [6] Yicang Zhou, Zhien Ma . Global stability of a class of discrete age-structured SIS models with immigration. Mathematical Biosciences and Engineering, 2009, 6(2): 409-425. doi: 10.3934/mbe.2009.6.409
    [7] Zongwei Ma, Hongying Shu . Viral infection dynamics in a spatial heterogeneous environment with cell-free and cell-to-cell transmissions. Mathematical Biosciences and Engineering, 2020, 17(3): 2569-2591. doi: 10.3934/mbe.2020141
    [8] Yujie Sheng, Jing-An Cui, Songbai Guo . The modeling and analysis of the COVID-19 pandemic with vaccination and isolation: a case study of Italy. Mathematical Biosciences and Engineering, 2023, 20(3): 5966-5992. doi: 10.3934/mbe.2023258
    [9] Zhen Jin, Guiquan Sun, Huaiping Zhu . Epidemic models for complex networks with demographics. Mathematical Biosciences and Engineering, 2014, 11(6): 1295-1317. doi: 10.3934/mbe.2014.11.1295
    [10] Qian Yan, Xianning Liu . Dynamics of an epidemic model with general incidence rate dependent on a class of disease-related contact functions. Mathematical Biosciences and Engineering, 2023, 20(12): 20795-20808. doi: 10.3934/mbe.2023920
  • In this study, bacteria associated with licorice (Glycyrrhiza glabra L.) were characterized through 16S rRNA gene analysis. Profiling of endophytic bacteria isolated from Glycyrrhiza glabra tissues revealed 18 isolates across the following genera: Enterobacter (4), Pantoea (3), Bacillus (2), Paenibacillus (2), Achromobacter (2), Pseudomonas (1), Escherichia (1), Klebsiella (1), Citrobacter (1), and Kosakonia (1). Furthermore, the beneficial features of bacterial isolates for plants were determined. The bacterial isolates showed the capacity to produce siderophores, hydrogen cyanide (HCN), indole-3-acetic acid (IAA), chitinase, protease, glucanase, lipase, and other enzymes. Seven bacterial isolates showed antagonistic activity against F. culmorum, F. solani, and R. solani. According to these results, licorice with antimicrobial properties may serve as a source for the selection of microorganisms that have antagonistic activity against plant fungal pathogens and may be considered potential candidates for the control of plant pathogens. The selected bacterial isolates, P. polymyxa GU1, A. xylosoxidans GU6, P. azotoformans GU7, and P. agglomerans GU18, increased root and shoot growth of licorice and were able to colonize the plant root. They can also serve as an active part of bioinoculants, improving plant growth.



    Population migration is a common phenomenon. With the migration of population, infectious diseases can easily spread from one area to another, so it is meaningful to consider population migration when studying the spread of infectious diseases [1,2,3,4,5,6,7].

    Wang and Mulone [2] established an SIS infectious disease model with standard incidence based on two patches. It is proved that the basic reproduction number is the threshold of the uniform persistence and disappearance of the disease. The dispersal rate of the population will make the infectious disease persist or disappear in all patches. There will be no the phenomena that infectious diseases persists in one patch but disappears in the other.

    Sun et al. [3] put forward an SIS epidemic model with media effect in a two patches setting. Under the assumption that the migration matrix is irreducible, it is proved that if the basic reproduction number is greater than 1 then the system persists and solutions converge to an endemic equilibrium and that if the basic reproduction number is less than 1 then solutions tend to an equilibrium without disease.

    Gao et al. [5,6] studied an SIS multi-patch model with variable transmission coefficients. Their results show that the basic reproduction number R0 is a threshold parameter of the disease dynamics.

    All the patch models referenced above assume that the migration matrix is irreducible. The studies in which the migration matrix is reducible are few. Therefore, based on the case of two patches, we consider that the individuals can only migrate from one patch to the other. In this case, the migration matrix is reducible. It can characterize the phenomenon that individuals migrate in one direction between two regions, such as, from the rural patch to the urban one [8] and from a small community hospital to a large teaching hospital [4].

    It is well known that the incidence rate plays an important role in the modeling of infectious disease. Considering the saturation phenomenon for numerous infected individuals, Capasso and Serio [9] first introduce a nonlinear bounded function g(I) to form the interaction term g(I)S in 1978. It can characterize the behavioral changes of individuals, such as wearing masks or reducing their social activities and direct contact with others with the increase of infectious individuals. After that, the saturation incidence rate has attracted much attention and various nonlinear types of incidence rate are employed. The most commonly used types are Holling type Ⅱ λSI1+αI [10,11,12] and βSI1+αS [13], Monod-Haldane type λSI1+αI2 [14], Beddington-DeAngelis type λSI1+αS+βI [15,16,17] and Crowley-Martin type λSI(1+αS)(1+βI) [18,19].

    In this paper, we consider infectious disease transmission models with saturation incidence rate. The rest of this paper is organized as follows: In Section 2, we establish a two-patch SIS model with saturating contact rate and one-directing population dispersal. We discuss the existence of disease-free equilibrium, boundary equilibrium and endemic equilibrium and prove the global asymptotic stability of the equilibriums in Section 3. In Section 4, we perform simulations to illustrate the results and analyze the effect of the contact rate and population migration on epidemic transmission. Finally, we discuss in Section 5.

    In the two patches, the population is divided into two states: susceptible and infective. Thus we can establish a two-patch SIS model with saturating contact rate and one-directing population dispersal

    {dS1(t)dt=A1d1S1β1S1I11+α1I1mS1+γ1I1,dS2(t)dt=A2d2S2β2S2I21+α2I2+mS1+γ2I2,dI1(t)dt=β1S1I11+α1I1d1I1mI1γ1I1,dI2(t)dt=β2S2I21+α2I2d2I2+mI1γ2I2, (2.1)

    where Si is the number of susceptible population in patch i (i=1,2), Ii is the number of infective population in patch i (i=1,2), Ai is the recruitment into patch i (i=1,2), di is the natural mortality rate, γi is the recovery rate of an infective individual in patch i (i=1,2), m is the migration rate form patch 1 to patch 2. Since the individuals can migrate from the first patch to the second, patch 1 is the source patch and patch 2 is the sink patch. The initial conditions is

    Si(0)>0,  Ii(0)0,  i=1,2,  I1(0)+I2(0)>0. (2.2)

    Denote the population in patch i by Ni. Then Ni=Si+Ii. From system (2.1), the differential equations governing the evolution of N1 and N2 are

    {dN1(t)dt=A1(d1+m)N1,dN2(t)dt=A2d2N2+mN1. (2.3)

    Obviously, system (2.3) has a unique equilibrium (N1,N2)=(A1d1+m,A2d2+mA1d2(d1+m)) which is globally asymptotically stable for (2.3). So (2.1) is equivalent the following system

    {dN1(t)dt=A1(d1+m)N1,dN2(t)dt=A2d2N2+mN1,dI1(t)dt=β1(N1I1)I11+α1I1d1I1mI1γ1I1,dI2(t)dt=β2(N2I2)I21+α2I2d2I2+mI1γ2I2. (2.4)

    Because limtNi(t)Ni (i=1,2), system (2.4) leads to the following limit system

    {dI1(t)dt=β1(N1I1)I11+α1I1d1I1mI1γ1I1,dI2(t)dt=β2(N2I2)I21+α2I2d2I2+mI1γ2I2. (2.5)

    Let Ω={(I1,I2)|0I1N1,0I2N2}. Then Ω is invariant region for system (2.5).

    Define the basic reproduction number in the two patches respectively by R10=β1A1(d1+γ1+m)2=β1d1+γ1+mN1, R20=β2(d2+γ2)N2. The basic reproduction number R10 gives the expected secondary infections in the source patch produced by a primary infected individual in the source patch when the population is supposed to be in the disease-free equilibrium. The basic reproduction number R20 gives the expected secondary infections in the sink patch produced by a primary infected individual in the sink patch when the population is supposed to be in the disease-free equilibrium. Then we have the following theorem.

    Theorem 3.1. For the system (2.5), we have

    (i) The disease-free equilibrium E0:=(0,0) always exists;

    (ii) The boundary equilibrium E1:=(0,β2N2d2γ2(d2+γ2)α2+β2) exists if R20>1;

    (iii) There is a unique epidemic equilibrium E if R10>1.

    Proof. (ⅰ) can be easily proved.

    Let

    β1(N1I1)I11+α1I1d1I1mI1γ1I1=0, (3.1)
    β2(N2I2)I21+α2I2d2I2+mI1γ2I2=0. (3.2)

    From Eq (3.1), we can have I1=0 always satisfies Eq (3.1). When I1=0, from Eq (3.2), we have

    I2=β2N2d2γ2(d2+γ2)α2+β2.

    If R20>1, then I2=β2N2d2γ2(d2+γ2)α2+β2>0. So The boundary equilibrium E1:=(0,β2N2d2γ2(d2+γ2)α2+β2) exists if R20>1. The conclusion (ⅱ) is proved.

    If R10>1, Eq (3.1) has a positive solution I1=β1N1(d1+m+γ1)(d1+m+γ1)α1+β1. Solve Eq (3.2), we have

    I2=(β2N2d2γ2+mα2I1)±(β2N2d2γ2+mα2I1)2+4[(d2+γ2)α2+β2]mI12[(d2+γ2)α2+β2]. (3.3)

    Substituting I1 into Eq (3.3), we have

    I2=(β2N2d2γ2+mα2I1)±(β2N2d2γ2+mα2I1)2+4[(d2+γ2)α2+β2]mI12[(d2+γ2)α2+β2].

    Since I20 is meaning only, we take

    I2=(β2N2d2γ2+mα2I1)+(β2N2d2γ2+mα2I1)2+4[(d2+γ2)α2+β2]mI12[(d2+γ2)α2+β2].

    So if R10>1, there is a unique epidemic equilibrium E=(I1,I2), where I1=β1N1(d1+m+γ1)(d1+m+γ1)α1+β1 and I2=(β2N2d2γ2+mα2I1)+(β2N2d2γ2+mα2I1)2+4[(d2+γ2)α2+β2]mI12[(d2+γ2)α2+β2]. The conclusion (ⅲ) is proved.

    This completes the proof of the theorem.

    From the above analysis, we have the following theorem.

    Theorem 3.2. For the system (2.5), we have

    (i) If R10<1 and R20<1, there is the disease-free equilibrium E0 only;

    (ii) If R10<1 and R20>1, there are the disease-free equilibrium E0 and the boundary equilibrium E1;

    (iii) If R10>1 and R20<1, there are the disease-free equilibrium E0 and the epidemic equilibrium E;

    (iv) If R10>1 and R20>1, there are the disease-free equilibrium E0, the boundary equilibrium E1 and the epidemic equilibrium E.

    Remark 3.1. Define the basic reproduction number R0 of the system (2.5) by the spectral radius of the next generation matrix [20], we have

    R0=ρ(β1d1+γ1+mN10mβ2N2(d2+γ2)(d1+γ1+m)β2(d2+γ2)N2),

    where ρ(A) denotes the spectral radius of a matrix A. So from the above analysis, we know that R0=max{R10,R20}.

    The next, we shall discuss the local stability of the disease-free equilibrium firstly. Then we discuss the global asymptotical stability.

    Theorem 3.3. For the system (2.5), we have

    (i) If R10<1 and R20<1, the disease-free equilibrium E0 is locally asymptotically stable;

    (ii) If R10>1 or R20>1, the disease-free equilibrium E0 is unstable.

    Proof. The linearized system of (2.5) at the equilibrium E0 is

    {dI1(t)dt=(β1N1d1mγ1)I1,dI2(t)dt=(β2N2d2γ2)I2+mI1. (3.4)

    The associated characteristic equation of the linearized system of (3.4) at the equilibrium E0 is

    F(λ)=|λ(β1N1(d1+m+γ1))0mλ(β2N2d2γ2)|=0 (3.5)

    It is easy to see that the two eigenvalues of characteristic Eq (3.5) are

    λ1=β1N1(d1+m+γ1)=(R101)(d1+m+γ1)

    and

    λ2=β2N2d2γ2=(R201)(d2+γ2).

    So, when R10<1 and R20<1, the disease-free equilibrium E0 is locally asymptotically stable; However, if R10>1 or R20>1, the disease-free equilibrium E0 is unstable.

    Remark 3.2. From Theorem 3.3, we know that for the system (2.5), if R0<1 the disease-free equilibrium E0 is locally asymptotically stable; if R0>1, the disease-free equilibrium E0 is unstable.

    Theorem 3.4. For the system (2.5), if R10<1 and R20<1, the disease-free equilibrium E0 is globally asymptotically stable.

    Proof. Since Ii1+αiIiIi for i=1,2, from system (2.5), we can obtain that

    {dI1(t)dt(β1N1d1mγ1)I1,dI2(t)dt(β2N2d2γ2)I2+mI1. (3.6)

    Define an auxiliary linear system using the right hand side of (3.6) as follows

    {dI1(t)dt=(β1N1d1mγ1)I1,dI2(t)dt=(β2N2d2γ2)I2+mI1.

    It can be rewritten as

    (I1I2)=(β1N1d1mγ10mβ2N2d2γ2)(I1I2). (3.7)

    if R10<1 and R20<1, we can solve (3.7) and know that limtI1(t)=0 and limtI2(t)=0. By the comparison principle [21], we can conclude that when R10<1 and R20<1, all non-negative solutions of (2.5) satisfy limtIi(t)=0 for i=1,2. So the disease-free equilibrium E0 is globally asymptotically stable.

    In this subsection, we will discuss the local stability of the boundary equilibrium firstly. Then discuss the global asymptotical stability.

    Theorem 3.5. For the system (2.5), if R10<1 and R20>1, the boundary equilibrium E1 is globally asymptotically stable.

    Proof. The Jacobian matrix at the boundary equilibrium E1 of system (2.5) is

    J=(β1N1d1mγ10m(d2+γ2)(1R20)(β2α2+α22)((d2+γ2)(R201)(d2+γ2)α2+β2)2(1+α2(d2+γ2)(R201)(d2+γ2)α2+β2)2).

    The two eigenvalues of the Jacobian matrix are

    λ1=β1N1(d1+m+γ1)=(R011)(d1+m+γ1)

    and

    λ2=(d2+γ2)(1R02)(β2α2+α22)((d2+γ2)(R21)(d2+γ2)α2+β2)2(1+α2(d2+γ2)(R21)(d2+γ2)α2+β2)2.

    So, when R10<1 and R20>1, λ1<0 and λ2<0. That is the boundary equilibrium E1 is locally asymptotically stable.

    For every (I1(0),I2(0))Ω, assume the solution of the system (2.5) with initial value (I1(0),I2(0)) is (I1(t),I2(t)). Since

    dI1(t)dt=(β1+(d1+γ1)α1+mα1)(β1N1d1mγ1β1+(d1+γ1)α1+mα1I1)I11+α1I1,

    if R10<1, dI1(t)dt<0, then I1(t) is positive and decreasing and limtI1(t)=0. So for sufficiently small positive number ϵ1, there exists a T, such that I1(T)=ϵ1 and I1(t)<ϵ1 when t>T.

    The following, we prove that for any ϵ>0, there exists a T>T such that |I2(T)β2N2d2γ2β2+(d2+γ2)α2|<ϵ. And because E1=(0,β2N2d2γ2β2+(d2+γ2)α2) is locally asymptotically stable, we have E1 is globally asymptotically stable.

    Since

    dI2(t)dt=(β2+(d2+γ2)α2)(β2N2d2γ2β2+(d2+γ2)α2I2)I21+α2I2+mI1,

    if I2(T)<β2N2d2γ2β2+(d2+γ2)α2, then dI2(t)dt>0 for t>T. So I2(t) is increasing and there exists T1 such that |I2(T1)β2N2d2γ2β2+(d2+γ2)α2|<ϵ;

    if I2(T)>β2N2d2γ2β2+(d2+γ2)α2, there are two cases:

    ⅰ) I2(t) is decreasing for t>T. In this case, there exists T2>T, such that |I2(T2)β2N2d2γ2β2+(d2+γ2)α2|<ϵ;

    ⅱ) There exists T1>T, such that dI2(T1)dt>0. That is

    (β2+(d2+γ2)α2)(β2N2d2γ2β2+(d2+γ2)α2I2(T1))I2(T1)1+α2I2(T1)+mI1(T1)>0.

    Since I1(t)<ϵ1 for t>T, we have

    (β2+(d2+γ2)α2)(β2N2d2γ2β2+(d2+γ2)α2I2(T1))I2(T1)1+α2I2(T1)+mϵ1>0.

    Since I2(T)>β2N2d2γ2β2+(d2+γ2)α2, we have

    (β2+(d2+γ2)α2)(β2N2d2γ2β2+(d2+γ2)α2I2(T1))α2+β2+(d2+γ2)α2β2N2d2γ2+mϵ1>0.

    So

    I2(T1)β2N2d2γ2β2+(d2+γ2)α2<α2+β2+(d2+γ2)α2β2N2d2γ2(d2+γ2)α2+β2mϵ1.

    If only ϵ1<(d2+γ2)α2+β2(α2+β2+(d2+γ2)α2β2N2d2γ2)mϵ, then |I2(T2)β2N2d2γ2β2+(d2+γ2)α2|<ϵ. It is completed.

    In this subsection, we will discuss the local stability of the epidemic equilibrium firstly and then discuss the global asymptotical stability.

    Theorem 3.6. For the system (2.5), if R10>1, the epidemic equilibrium E is locally asymptotically stable.

    Proof. The Jacobian matrix at the epidemic equilibrium E of system (2.5) is

    J=(β1I11+α1I1+β1(N1I1)(1+α1I1)2d1mγ10mβ2N22β2I22β2α2I22(d2+γ2)(1+α2I2)2(1+α2I2)2)=((d1+m+γ1)1R10(β1α1d1+m+γ1+α21)I21(1+α1I1)20mβ2N22β2I22β2α2I22(d2+γ2)(1+α2I2)2(1+α2I2)2).

    The two eigenvalues of the Jacobian matrix are

    λ1=(d1+m+γ1)1R10(β1α1d1+m+γ1+α21)I21(1+α1I1)2

    and

    λ2=β2N22β2I22β2α2I22(d2+γ2)(1+α2I2)2(1+α2I2)2.

    It is easy to see that if R10>1, λ1<0. The next, we need only prove the second eigenvalue λ2<0 if R10>1. Since (1+α2I2)2>0, we need only prove β2N22β2I22β2α2I22(d2+γ2)(1+α2I2)2<0. Let

    G(I2)=β2N22β2I22β2α2I22(d2+γ2)(1+α2I2)2.

    Since I2(β2N2d2γ2+mα1I1)(d2+γ2)α2+β2, so

    G(I2)=β2N22β2I22β2α2I22(d2+γ2)(1+α2I2)2β2N2(β2+(d2+γ2)α2)(β2N2d2γ2+mα1I1)(d2+γ2)α2+β2d2γ2=mα1I1<0.

    This completes the proof.

    Theorem 3.7. For the systeml (2.5), if R10>1, the epidemic equilibrium E is globally asymptotically stable.

    Proof. Since E is stable when R10>1, we need only prove E is globally attractive.

    Consider the equation

    dI1(t)dt=I1(β1(N1I1)1+α1I1d1mγ1).

    Let

    f1(I1)=β1(N1I1)1+α1I1d1mγ1.

    Then f1(I1)=β11N1α1(1+α1I1)2<0. So f1(I1) is a monotonic decreasing function for all I1>0. Furthermore, f1(0)>0, f1(N1)<0 and f1(I1)=0 when R10>1. That means if I1(0,I1), f1(I1)>0 and dI1(t)dt>0; if I1(I1,N1), f1(I1)<0 and dI1(t)dt<0. Hence limtI1(t)=I1. By Eq (3.3), limtI2(t)=I2. Thus E is globally asymptotically stable.

    The results about the existence and stability of equilibria are summarized in Table 1.

    Table 1.  Existence and stability of equilibria.
    Conditions E0 E1 E
    R10<1 and R20<1 Yes (GAS) No No
    R10<1 and R20>1 Yes (Unstable) Yes (GAS) No
    R10>1 and R20<1 Yes (Unstable) No Yes (GAS)
    R10>1 and R20>1 Yes (Unstable) Yes (Unstable) Yes (GAS)

     | Show Table
    DownLoad: CSV

    Remark 3.3. From Theorems 3.5 and 3.7, we know that for the system (2.5), if R0>1, the infectious disease is uniformly persistent. However, the infectious disease is not always uniformly persistent in every patch. If R0>1, but R10<1, the disease is uniformly persistent in the sink patch, but is extinct in the source patch. If R10>1, the disease is always uniformly persistent in every patch. This is a different conclusion resulted by the reducible migration matrix.

    Remark 3.4. From Theorem 3.5, in the case that R10<1 and R20>1, the infection does not persist in the source patch but is able to persist in the sink patch. So, in the early stage of the spread of infectious disease, the sink patch should assess the reproduction number R20 reasonably and take control measures timely to prevent the epidemic.

    In this section, we carry on numerical simulations to verify the theoretical conclusions, reveal the influence of the migration rate form patch 1 to patch 2 on the basic reproduction number, the transmission scale and transmission speed, and discuss the influence of the parameters α1 and α2 that measure the inhibitory effect on the basic reproduction number, the transmission scale and transmission speed.

    To numerically illustrate the theoretical results, we need to choose some parameter values (see Table 2).

    Table 2.  Description and values of parameters.
    Parameter Description Value
    A1 the recruitment rate of the population in patch 1 0.018 (Figures 1 and 2)
    0.03 (Figure 3)
    A2 the recruitment rate of the population in patch 2 0.0005 (Figure 1)
    0.004 (Figures 2 and 3)
    β1 the transmission rate in patch 1 0.00001
    β2 the transmission rate in patch 2 0.00005
    α1 the parameter that measure the inhibitory effect in patch 1 0.02
    α2 the parameter that measure the inhibitory effect in patch 1 0.02
    d1 the death rate in patch 1 0.0003
    d2 the death rate in patch 2 0.0003
    γ1 the death rate in patch 1 0.0001
    γ2 the death rate in patch 2 0.0001
    m the migration rate form the patch 1 to the patch 2 0.00005

     | Show Table
    DownLoad: CSV

    We verify the theoretical conclusions firstly. Denote the density of the infective individuals in patch 1 by i1(t)=I1(t)N1. Denote the density of the infective individuals in patch 2 by i2(t)=I2(t)N2. Figure 1 shows the evolution of the density of infective individuals in the two patches when R10=0.8889 and R20=0.7500. As predicted by the analytic calculation, the infectious disease in the two patches will disappear eventually. Figure 2 shows the evolution of the density of infective individuals in the two patches when R10=0.8889 and R20=1.8750. We can see the infectious disease will be endemic in patch 2 and the infectious disease in patch 1 will disappear eventually. Figure 3 shows the evolution of the density of infective individuals in the two patches when R10=1.4815 and R20=2.2917. We can see the infectious disease will be endemic in the two patches. And we can see the infectious disease will be endemic in the two patches if R10>1 from the subfigures (c) and (d) of Figure 4.

    Figure 1.  When R10=0.8889 and R20=0.7500, the evolution of the density of the infective individuals in the two patches.
    Figure 2.  When R10=0.8889 and R20=1.8750, the evolution of the density of the infective individuals in the two patches.
    Figure 3.  When R10=1.4815 and R20=2.2917, the evolution of the density of the infective individuals in the two patches.
    Figure 4.  Phase portraits for (a) R10<1 and R20<1; (b) R10<1 and R20>1; (c) R10>1 and R20<1; (d) R10>1 and R20>1.

    Second, we reveal the influence of the migration rate m on the transmission in Figure 5. With the increasing of m, the density of infective individuals in patch 1 i1 is decreasing, however the density of infective individuals in patch 2 i2 is increasing.

    Figure 5.  When m is increasing, (a) the density of infective individuals in patch 1 i1 is decreasing; (b) the density of infective individuals in patch 2 i2 is increasing.

    Third, we reveal the parameters α1 and α2 on the transmission scale and transmission speed. We can see that when α1 is increasing, the density of infective individuals in patch 1 i1 is decreasing from Figure 6 and When α2 is increasing, the density of infective individuals in patch 2 i1 is decreasing from Figure 7.

    Figure 6.  When α1 is increasing, the density of infective individuals in patch 1 i1 is decreasing.
    Figure 7.  When α2 is increasing, the density of infective individuals in patch 2 i1 is decreasing.

    Many scholars have studied infectious disease transmission with population migration [1,2,3,4,5,6,7], assuming that the migration matrix is irreducible, and found that the propagation dynamics of infectious diseases is determined by the basic reproduction number of the system. When the basic reproduction number is less than 1, the infectious disease eventually becomes extinct; when the basic reproduction number is larger than 1, the infectious disease is epidemic eventually. Since the migration matrix is irreducible, all patches are a connected whole. In all patches, infectious diseases are either extinct or epidemic. That is there is not the phenomenon that infectious diseases are extinct in some patches but epidemic in the others.

    Because the studies about the spread of infectious diseases with reducible migration matrix are rare, in this paper, we proposed a two-patch SIS model with saturating contact rate and one-directing population dispersal, discussed the global asymptotic stability of the disease-free equilibrium, the boundary equilibrium and the endemic equilibrium respectively, and revealed the influence of saturating contact rate and migration rate on basic reproduction number and the transmission scale. We have the following main conclusions:

    1) If R10>1 then the system tends to a global endemic equilibrium in which infected individuals are present in both patches provided initially there were infected individuals in the source patch; If R10<1 and R20>1 then the system converges to an equilibrium with infected individuals only in the sink patch; If R10<1 and R20<1 then the system converges to the disease-free equilibrium.

    2) When migration rate is increasing, the density of infective individuals in the source patch is decreasing; but the density of infective individuals in the sink patch is increasing;

    3) With the increasing of the parameter αi (i=1,2) in saturating contact rate, the density of infective individuals in patch i (i=1,2) is decreasing.

    The similar conclusions can be obtained for the two patch SI model

    {dS1(t)dt=A1d1S1β1S1I11+α1I1mS1,dS2(t)dt=A2d2S2β2S2I21+α2I2+mS1,dI1(t)dt=β1S1I11+α1I1d1I1mI1,dI2(t)dt=β2S2I21+α2I2d2I2+mI1.

    We can generalize the current model in many aspects to increase realism. For instance, the infection rate can be given by β(I)SI. We can give the properties on function β(I) such that β(I) is decreasing and tends to 0 when I tends to infinity. The mortality rates of the susceptible and infected individuals are the same in the current model. In fact, the disease-induced death rate can not be neglected sometimes. So the disease-induced death rate can be considered. It is also significant to consider heterogeneous number of contacts for each individual on complex network. There are many paper on this topic [22,23]. One can investigate the multi-patch epidemic model with reducible migration matrix.

    This work is supported by the National Natural Sciences Foundation of China (Nos.12001501, 12071445, 11571324, 61603351), Shanxi Province Science Foundation for Youths (201901D211216), the Fund for Shanxi '1331KIRT'.

    All authors declare no conflicts of interest in this paper.


    Acknowledgments



    This research was supported by the grant funded by Innovation Development Agency of Uzbekistan.

    Conflict of interest



    The authors declare no conflict of interest.

    Author contributions



    DE designed the experiment. GM conducted the laboratory experiments. VS and FE analyzed the results of experiments. DE, VS, and FE wrote the manuscript. All authors read and approved the manuscript.

    [1] Santos MLD, Berlitz DL, Wiest SLF, et al. (2018) Benefits associated with the interaction of endophytic bacteria and plants. Braz Arch Biol Technol 61: 1-11. https://doi.org/10.1590/1678-4324-2018160431
    [2] Egamberdieva D, Alimov J, Shurigin V, et al. (2022) Diversity and plant growth-promoting ability of endophytic, halotolerant bacteria associated with Tetragonia tetragonioides (Pall.) Kuntze. Plants 11: 49. https://doi.org/10.3390/plants11010049
    [3] Egamberdieva D, Ma H, Reckling M, et al. (2022) Interactive effects of biochar and N and P nutrients on the symbiotic performance, growth, and nutrient uptake of soybean (Glycine max L.). Agronomy 12: 27. https://doi.org/10.3390/agronomy12010027
    [4] Weyens N, van der Lelie D, Taghavi S, et al. (2009) Phytoremediation: Plant–endophyte partnerships take the challenge. Curr Opin Biotechnol 20: 248-254. https://doi.org/10.1016/j.copbio.2009.02.012
    [5] Shurigin V, Alaylar B, Davranov K, et al. (2021) Diversity and biological activity of culturable endophytic bacteria associated with marigold (Calendula officinalis L.). AIMS Microbiol 7: 336-353. https://doi.org/10.3934/microbiol.2021021
    [6] Shurigin V, Egamberdieva D, Li L, et al. (2020) Endophytic bacteria associated with halophyte Seidlitzia rosmarinus Ehrenb. ex Boiss. from saline soil of Uzbekistan and their plant beneficial traits. J Arid Land 12: 730-740. https://doi.org/10.1007/s40333-020-0019-4
    [7] Marui A, Nagafuchi T, Shinogi Y, et al. (2012) Soil physical properties to grow the wild licorice at semi-arid area in Mongolia. J Arid Land Studies 22: 33-36.
    [8] Lewis G, Schrire B, Mackinder B, et al. (2005) Legumes of the world. London: Royal Botanic Gardens, Kew Publishing.
    [9] Hayashi H, Hattori S, Inoue K, et al. (2003) Field survey of Glycyrrhiza plants in Central Asia (1). Characterization of G. uralensis, G. glabra and the putative intermediate collected in Kazakhstan. Biol Pharm Bull 26: 867-871. https://doi.org/10.1248/bpb.26.867
    [10] Patil SM, Patil MB, Sapkale GN (2009) Antimicrobial activity of Glycyrrhiza glabra Linn. roots. Int J Chem Sci 7: 585-591.
    [11] Sharma V, Agrawal RC, Pandey S (2014) Phytochemical screening and determination of anti-bacterial and anti-oxidant potential of Glycyrrhiza glabra root extracts. J Envir Res Devel 7: 1552-1558.
    [12] Chao H, Wang W, Hou J (2019) Plant growth and soil microbial impacts of enhancing licorice with inoculating dark septate endophytes under drought stress. Front Microbiol 10: 2277. https://doi.org/10.3389/fmicb.2019.02277
    [13] Li L, Sinkko H, Montonen L, et al. (2012) Biogeography of symbiotic and other endophytic bacteria isolated from medicinal Glycyrrhiza species in China. FEMS Microbiol Ecol 79: 46-68. https://doi.org/10.1111/j.1574-6941.2011.01198.x
    [14] Cao XM, Cai J, Li SB, et al. (2013) Fusarium solani and Fusarium oxysporum associated with root rot of Glycyrrhiza uralensis in China. Plant Dis 97: 1514. https://doi.org/10.1094/PDIS-12-12-1111-PDN
    [15] Parray J, Egamberdieva D, Abd_Allah E, et al. (2023) Editorial: Soil microbiome metabolomics: A way forward to sustainable intensification. Front Sustain Food Syst 7: 1251054. https://doi.org/10.3389/fsufs.2023.1251054
    [16] Abdullaeva Y, Mardonova G, Eshboev F, et al. (2024) Harnessing chickpea bacterial endophytes for improved plant health and fitness. AIMS Microbiol 10: 489-506. https://doi.org/10.3934/microbiol.2024024
    [17] Egamberdieva D, Wirth S, Alqarawi AA, et al. (2017) Phytohormones and beneficial microbes: Essential components for plants to balance stress and fitness. Front Microbiol 8: 2104. https://doi.org/10.3389/fmicb.2017.02104
    [18] Rezaei-Chiyaneh E, Mahdavikia H, Subramanian S, et al. (2021) Co-inoculation of phosphate-solubilizing bacteria and mycorrhizal fungi: Effect on seed yield, physiological variables, and fixed oil and essential oil productivity of ajowan (Carum copticum L.) under water deficit. J Soil Sci Plant Nutr 21: 3159-3179. https://doi.org/10.1007/s42729-021-00596-9
    [19] Pawlik M, Cania B, Thijs S, et al. (2017) Hydrocarbon degradation potential and plant growth-promoting activity of culturable endophytic bacteria of Lotus corniculatus and Oenothera biennis from a long-term polluted site. Environ Sci Pollut Res 24: 19640-19652. https://doi.org/10.1007/s11356-017-9496-1
    [20] Chamkhi I, Sbabou L, Aurag J (2023) Improved growth and quality of saffron (Crocus sativus L.) in the field conditions through inoculation with selected native plant growth-promoting rhizobacteria (PGPR). Ind Crops Prod 197: 116606. https://doi.org/10.1016/j.indcrop.2023.116606
    [21] Egamberdieva D, Shurigin V, Alaylar B, et al. (2020) Bacterial endophytes from horseradish (Armoracia rusticana G. Gaertn., B. Mey. & Scherb.) with antimicrobial efficacy against pathogens. Plant Soil Environ 66: 309-316. https://doi.org/10.17221/137/2020-PSE
    [22] Egamberdieva D, Shurigin V, Alaylar B, et al. (2020) The effect of biochars and endophytic bacteria on growth and root rot disease incidence of Fusarium infested narrow-leafed lupin (Lupinus angustifolius L.). Microorganisms 8: 496. https://doi.org/10.3390/microorganisms8040496
    [23] Egamberdieva D, Wirth S, Behrendt U, et al. (2017) Antimicrobial activity of medicinal plants correlates with the proportion of antagonistic endophytes. Front Microbiol 8: 199. https://doi.org/10.3389/fmicb.2017.00199
    [24] Nejatzadeh-Barandozi F (2013) Antibacterial activities and antioxidant capacity of Aloe vera. Bioorganic Med Chem Lett 3: 1-8. https://doi.org/10.1186/2191-2858-3-5
    [25] Bafana A, Lohiya R (2013) Diversity and metabolic potential of culturable root-associated bacteria from Origanum vulgare in sub-Himalayan region. World J Microbiol Biotechnol 29: 63-74. https://doi.org/10.1007/s11274-012-1158-3
    [26] Phetcharat P, Duangpaeng A (2012) Screening of endophytic bacteria from organic rice tissue for indole acetic acid production. Procedia Eng 32: 177-183. https://doi.org/10.1016/j.proeng.2012.01.1254
    [27] Katoch M, Pull S (2017) Endophytic fungi associated with Monarda citriodora, an aromatic and medicinal plant and their biocontrol potential. Pharm Biol 55: 1528-1535. https://doi.org/10.1080/13880209.2017.1309054
    [28] Farhaoui A, El Alami N, Khadiri M, et al. (2023) Biological control of diseases caused by Rhizoctonia solani AG-2-2 in sugar beet (Beta vulgaris L.) using plant growth-promoting rhizobacteria (PGPR). Physiol Mol Plant Pathology 124: 101966. https://doi.org/10.1016/j.pmpp.2023.101966
    [29] Rustamova N, Wubulikasimu A, Mukhamedov N, et al. (2020) Endophytic bacteria associated with medicinal plant Baccharoides anthelmintica diversity and characterization. Curr Microbiol 77: 1457-1465. https://doi.org/10.1007/s00284-020-01924-5
    [30] Mora-Ruiz MDR, Font-Verdera F, Díaz-Gil C, et al. (2015) Moderate halophilic bacteria colonizing the phylloplane of halophytes of the subfamily Salicornioideae (Amaranthaceae). Syst Appl Microbiol 38: 406-416. https://doi.org/10.1016/j.syapm.2015.05.004
    [31] Dashti AA, Jadaon MM, Abdulsamad AM, et al. (2009) Heat treatment of bacteria: A simple method of DNA extraction for molecular techniques. Kuwait Med J 41: 117-122.
    [32] Lane DJ (1991) 16S/23S rRNA Sequencing. Nucleic Acid Techniques in Bacterial Systematic. New York: John Wiley and Sons 115-175.
    [33] Jinneman KC, Wetherington JH, Adams AM, et al. (1996) Differentiation of Cyclospora sp. and Eimeria spp. by using the polymerase chain reaction amplification products and restriction fragment length polymorphisms. Food and Drug Administration Laboratory Information Bulletin LIB no 4044 .
    [34] Tamura K, Nei M, Kumar S (2004) Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA 101: 11030-11035. https://doi.org/10.1073/pnas.0404206101
    [35] Kumar S, Stecher G, Li M, et al. (2018) MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol 35: 1547-1549. https://doi.org/10.1093/molbev/msy096
    [36] Castric PA (1975) Hydrogen cyanide, a secondary metabolite of Pseudomonas aeruginosa. Can J Microbiol 21: 613-618. https://doi.org/10.1139/m75-088
    [37] Schwyn B, Neilands JB (1987) Universal chemical assay for the detection and determination of siderophores. Anal Biochem 160: 47-56. https://doi.org/10.1016/0003-2697(87)90612-9
    [38] Brown MRW, Foster JHS (1970) A simple diagnostic milk medium for Pseudomonas aeruginosa. J Clin Pathol 23: 172-177. https://doi.org/10.1136/jcp.23.2.172
    [39] Walsh GA, Murphy RA, Killeen GF, et al. (1995) Technical note: Detection and quantification of supplemental fungal b-glucanase activity in animal feed. J Anim Sci 73: 1074-1076. https://doi.org/10.2527/1995.7341074x
    [40] Malleswari D, Bagyanarayan G (2017) In vitro screening of rhizobacteria isolated from the rhizosphere of medicinal and aromatic plants for multiple plant growth promoting activities. J Microbiol Biotechnol Res 3: 84-91.
    [41] Howe TG, Ward JM (1976) The utilization of tween 80 as carbon source by Pseudomonas. J Gen Microbiol 92: 234-235. https://doi.org/10.1099/00221287-92-1-234
    [42] Bano N, Musarrat J (2003) Characterization of a new Pseudomonas aeruginosa strain NJ-15 as a potential biocontrol agent. Curr Microbiol 46: 324-328. https://doi.org/10.1007/s00284-002-3857-8
    [43] Egamberdieva D, Kucharova Z (2009) Selection for root colonising bacteria stimulating wheat growth in saline soils. Biol Fertil Soils 45: 561-573. https://doi.org/10.1007/s00374-009-0366-y
    [44] Egamberdieva D, Wirth SJ, Shurigin VV, et al. (2017) Endophytic bacteria improve plant growth, symbiotic performance of chickpea (Cicer arietinum L.) and induce suppression of root rot caused by Fusarium solani under salt stress. Front Microbiol 8: 1887. https://doi.org/10.3389/fmicb.2017.01887
    [45] Simons M, van der Bij AJ, Brand I, et al. (1996) Gnotobiotic system for studying rhizosphere colonization by plant growth-promoting Pseudomonas bacteria. Mol Plant Microbe Interact 9: 600-607. https://doi.org/10.1094/mpmi-9-0600
    [46] Makuwa SC, Motadi LR, Choene M, et al. (2023) Bacillus dicomae sp. nov., a new member of the Bacillus cereus group isolated from medicinal plant Dicoma anomala. Int J Syst Evol Microbiol 73: 006112. https://doi.org/10.1099/ijsem.0.006112
    [47] Khan MS, Gao J, Chen X, et al. (2020) Isolation and characterization of plant growth-promoting endophytic bacteria Paenibacillus polymyxa SK1 from Lilium lancifolium. Biomed Res Int 27: 8650957. https://doi.org/10.1155/2020/8650957
    [48] Anandan K, Vittal RR (2019) Endophytic Paenibacillus amylolyticus KMCLE06 extracted dipicolinic acid as antibacterial agent derived via dipicolinic acid synthetase gene. Curr Microbiol 76: 178-186. https://doi.org/10.1007/s00284-018-1605-y
    [49] Shurigin V, Alimov J, Davranov K, et al. (2022) The diversity of bacterial endophytes from Iris pseudacorus L. and their plant beneficial traits. Curr Res Microb Sci 3: 100133. https://doi.org/10.1016/j.crmicr.2022.100133
    [50] Egamberdieva D, Kamilova F, Validov S, et al. (2008) High incidence of plant growth-stimulating bacteria associated with the rhizosphere of wheat grown in salinated soil in Uzbekistan. Environ Microbiol 19: 1-9. https://doi.org/10.1111/j.1462-2920.2007.01424.x
    [51] Lim JA, Lee DH, Heu S (2014) The interaction of human enteric pathogens with plants. Plant Pathol J 30: 109-116. https://doi.org/10.5423/PPJ.RW.04.2014.0036
    [52] Cho ST, Chang HH, Egamberdieva D, et al. (2015) Genome analysis of Pseudomonas fluorescens PCL1751: A rhizobacterium that controls root diseases and alleviates salt stress for its plant host. PLoS ONE 10: e0140231. https://doi.org/10.1371/journal.pone.0140231
    [53] Koberl M, Ramadan EM, Adam M, et al. (2013) Bacillus and Streptomyces were selected as broad-spectrum antagonists against soilborne pathogens from arid areas in Egypt. FEMS Microbiol Lett 342: 168-178. https://doi.org/10.1111/1574-6968.12089
    [54] Egamberdiyeva D, Hoflich G (2004) Effect of plant growth-promoting bacteria on growth and nutrient uptake of cotton and pea in a semi-arid region of Uzbekistan. J Arid Environ 56: 293-301. https://doi.org/10.1016/S0140-1963(03)00050-8
    [55] Medison RG, Tan L, Medison MB, et al. (2022) Use of beneficial bacterial endophytes: A practical strategy to achieve sustainable agriculture. AIMS Microbiol 8: 624-643. https://doi.org/10.3934/microbiol.2022040
    [56] Ali S, Duan J, Charles TC, et al. (2014) A bioinformatics approach to the determination of genes involved in endophytic behavior in Burkholderia spp. J Theor Biol 343: 193-198. https://doi.org/10.1016/j.jtbi.2013.10.007
    [57] Goryluk A, Rekosz-Burlaga H, Blaszczyk M (2009) Isolation and characterization of bacterial endophytes of Chelidonium majus L. Pol J Microbiol 58: 355-361.
    [58] Akinsanya MA, Goh JK, Lim SP, et al. (2015) Metagenomics study of endophytic bacteria in Aloe vera using next-generation technology. Genom Data 6: 159-163. https://doi.org/10.1016/j.gdata.2015.09.004
    [59] Liu Y, Mohamad OAA, Salam N, et al. (2019) Diversity, community distribution and growth promotion activities of endophytes associated with halophyte Lycium ruthenicum Murr. 3 Biotech 9: 144. https://doi.org/10.1007/s13205-019-1678-8
    [60] Rana KL, Kour D, Yadav AH (2019) Endophytic microbiomes: Biodiversity, ecological significance and biotechnological applications. Res J Biotechnol 14: 142-162.
    [61] Michelsen CF, Stougaard P (2012) Hydrogen cyanide synthesis and antifungal activity of the biocontrol strain Pseudomonas fluorescens In5 from Greenland is highly dependent on growth medium. Can J Microbiol 58: 381-390. https://doi.org/10.1139/w2012-004
    [62] Wozniak M, Gałazka A, Tyskiewicz R, et al. (2019) Endophytic bacteria potentially promote plant growth by synthesizing different metabolites and their phenotypic/physiological profiles in the Biolog GEN III MicroPlateTM Test. Int J Mol Sci 20: 1-24. https://doi.org/10.3390/ijms20215283
    [63] Musa Z, Ma J, Egamberdieva D, et al. (2020) Diversity and antimicrobial potential of cultivable endophytic actinobacteria associated with medicinal plant Thymus roseus. Front Microbiol 11: 191. https://doi.org/10.3389/fmicb.2020.00191
    [64] Glick BR (2014) Bacteria with ACC deaminase can promote plant growth and help to feed the world. Microbiol Res 169: 30-39. https://doi.org/10.1016/j.micres.2013.09.009
    [65] Yadav A, Yadav K (2019) Plant growth-promoting endophytic bacteria and their potential to improve agricultural crop yields. Microbial Interventions in Agriculture and Environment. Singapore: Springer 143-169. https://doi.org/10.1007/978-981-32-9084-6_7
    [66] Fouda A, Eid AM, Elsaied A, et al. (2021) Plant growth promoting endophytic bacterial community inhabiting the leaves of Pulicaria incisa (Lam.) DC inherent to arid regions. Plants 10: 76. https://doi.org/10.3390/plants10010076
    [67] Shurigin V, Li L, Alaylar B, et al. (2024) Plant beneficial traits of endophytic bacteria associated with fennel (Foeniculum vulgare Mill.). AIMS Microbiol 10: 449-467. https://doi.org/10.3934/microbiol.2024022
    [68] Sudarshna, Sharma N (2024) Endophytic bacteria associated with critically endangered medicinal plant Trillium govanianum (Wall ex. Royle) and their potential in soil nutrition alleviation. Plant Stress 11: 100349. https://doi.org/10.1016/j.stress.2024.100349
    [69] Deepa N, Chauhan Sh, Singh A (2024) Unraveling the functional characteristics of endophytic bacterial diversity for plant growth promotion and enhanced secondary metabolite production in Pelargonium graveolens. Microbiol Res 283: 127673. https://doi.org/10.1016/j.micres.2024.127673
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(988) PDF downloads(85) Cited by(0)

Figures and Tables

Figures(5)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog