Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js
Review Special Issues

A review of molecular biology detection methods for human adenovirus

  • Received: 03 December 2022 Revised: 05 February 2023 Accepted: 15 February 2023 Published: 10 March 2023
  • In humans particularly in children, adenovirus is one of the most common viruses that cause respiratory illnesses. Knowing how to detect adenovirus proficiently and rapidly will help reinforce surveillance of adenovirus infections, detect epidemic situations in real-time, and understand the trend of virus epidemics, which will allow effective actions to be taken quickly. The rapid detection of antiviral antibodies or viral antigens in clinical samples can be achieved by molecular diagnostic techniques like PCR, Real-Time PCR, LAMP, mPCR-RLB, PCR-ELISA, Tem-PCR, Gene Chip, and so on. Some of the molecular diagnostic methods are relatively economical, exceedingly sensitive and explicit. There are several commercially accessible molecular diagnostic techniques that enable their use in clinical laboratories all over the world. In this review, the principles, characteristics, and applications of molecular biology surveillance methods commonly used in labs and clinics for the detection of human adenoviruses are examined and highlighted.

    Citation: Neelabh Datta. A review of molecular biology detection methods for human adenovirus[J]. AIMS Biophysics, 2023, 10(1): 95-120. doi: 10.3934/biophy.2023008

    Related Papers:

    [1] Lin Zhang, Yongbin Ge, Zhi Wang . Positivity-preserving high-order compact difference method for the Keller-Segel chemotaxis model. Mathematical Biosciences and Engineering, 2022, 19(7): 6764-6794. doi: 10.3934/mbe.2022319
    [2] Sunwoo Hwang, Seongwon Lee, Hyung Ju Hwang . Neural network approach to data-driven estimation of chemotactic sensitivity in the Keller-Segel model. Mathematical Biosciences and Engineering, 2021, 18(6): 8524-8534. doi: 10.3934/mbe.2021421
    [3] Thierry Colin, Marie-Christine Durrieu, Julie Joie, Yifeng Lei, Youcef Mammeri, Clair Poignard, Olivier Saut . Modeling of the migration of endothelial cells on bioactive micropatterned polymers. Mathematical Biosciences and Engineering, 2013, 10(4): 997-1015. doi: 10.3934/mbe.2013.10.997
    [4] Shangbing Ai, Zhian Wang . Traveling bands for the Keller-Segel model with population growth. Mathematical Biosciences and Engineering, 2015, 12(4): 717-737. doi: 10.3934/mbe.2015.12.717
    [5] Tong Li, Zhi-An Wang . Traveling wave solutions of a singular Keller-Segel system with logistic source. Mathematical Biosciences and Engineering, 2022, 19(8): 8107-8131. doi: 10.3934/mbe.2022379
    [6] Maghnia Hamou Maamar, Matthias Ehrhardt, Louiza Tabharit . A nonstandard finite difference scheme for a time-fractional model of Zika virus transmission. Mathematical Biosciences and Engineering, 2024, 21(1): 924-962. doi: 10.3934/mbe.2024039
    [7] Qigang Deng, Fugeng Zeng, Dongxiu Wang . Global existence and blow up of solutions for a class of coupled parabolic systems with logarithmic nonlinearity. Mathematical Biosciences and Engineering, 2022, 19(8): 8580-8600. doi: 10.3934/mbe.2022398
    [8] Z. Jackiewicz, B. Zubik-Kowal, B. Basse . Finite-difference and pseudo-spectral methods for the numerical simulations of in vitro human tumor cell population kinetics. Mathematical Biosciences and Engineering, 2009, 6(3): 561-572. doi: 10.3934/mbe.2009.6.561
    [9] Benjamin Wacker, Jan Christian Schlüter . A non-standard finite-difference-method for a non-autonomous epidemiological model: analysis, parameter identification and applications. Mathematical Biosciences and Engineering, 2023, 20(7): 12923-12954. doi: 10.3934/mbe.2023577
    [10] Fadoua El Moustaid, Amina Eladdadi, Lafras Uys . Modeling bacterial attachment to surfaces as an early stage of biofilm development. Mathematical Biosciences and Engineering, 2013, 10(3): 821-842. doi: 10.3934/mbe.2013.10.821
  • In humans particularly in children, adenovirus is one of the most common viruses that cause respiratory illnesses. Knowing how to detect adenovirus proficiently and rapidly will help reinforce surveillance of adenovirus infections, detect epidemic situations in real-time, and understand the trend of virus epidemics, which will allow effective actions to be taken quickly. The rapid detection of antiviral antibodies or viral antigens in clinical samples can be achieved by molecular diagnostic techniques like PCR, Real-Time PCR, LAMP, mPCR-RLB, PCR-ELISA, Tem-PCR, Gene Chip, and so on. Some of the molecular diagnostic methods are relatively economical, exceedingly sensitive and explicit. There are several commercially accessible molecular diagnostic techniques that enable their use in clinical laboratories all over the world. In this review, the principles, characteristics, and applications of molecular biology surveillance methods commonly used in labs and clinics for the detection of human adenoviruses are examined and highlighted.



    In this paper, the finite-difference method (FDM) is applied to solve the following two-dimensional (2D) Keller-Segel model [1,2]:

    {ut+(χuv)=dΔu,vt=Δvv+u,(x,y,t)Ω×(0,T]. (1.1)

    The model (1.1) describes the evolutionary process of cell density u(x,y,t) and a chemical stimulus (chemoattractant) concentration v(x,y,t) over a time t and location (x,y), where Ω={(x,y)|ax,yb}R2 is a convex bounded domain and a and b are constants. and Δ are gradient and Laplacian operators, respectively. χ>0 represents the chemotactic sensitivity constant; d>0 denotes the diffusion rate of cells. In addition, the initial conditions related with (1.1) are given as

    u(x,y,0)=u0(x,y),v(x,y,0)=v0(x,y),(x,y)Ω, (1.2)

    and the boundary conditions are assumed to be homogeneous Neumann boundary conditions, that is,

    un=vn=0,(x,y,t)Ω×(0,T], (1.3)

    where Ω is the boundary of Ω and n is the outward normal of Ω. With this condition (1.3), the total mass

    Massu=Ωu(x,y,t)dxdy=Ωu(x,y,0)dxdy

    is conserved as temporal evolution. Besides the above, the model (1.1) has the following form of free energy:

    E(u,v)(t)=Ω(uln(u)+χ2|v|2+χ2v2χuv)dxdy,t>0. (1.4)

    From mathematical analysis, the following equation can be verified by direct calculation of Eq. (1.4),

    E(u,v)t=Ω(u|(ln(u)χv)|2+χ(vt)2)dxdy0,t>0.

    Many early works show that the free energy (1.4) is decreasing over time and is mainly employed to demonstrate the existence of solutions for the chemotaxis system; see [3,4] and the references therein.

    Chemotaxis refers to the directional movement, which includes toward or away from the higher concentrations of cells or microorganisms that are stimulated by chemical substances in the external environment along the gradient directions of the concentration for stimuli. Chemotaxis phenomena play a crucial role in numerous intricate biological evolutions, such as bacterial aggregation, angiogenesis, pattern formation, embryonic development and so on[5]. From as early the 1950 to 1980, the chemotaxis phenomena have been extensively studied by many applied mathematicians and biologists, and a class of models of partial differential systems closely associated with taxis have been proposed; see [1,2,6,7,8,9,10]. Among them the most classical is the above-named system (1.1) (first proposed by Keller and Segel [1,2] in 1970), which simulated the aggregation phenomenon for Amoebae and Dictyostelium, as well as the traveling wave migration phenomenon for Escherichia coli in a capillary filled with nutrients.

    Since the model was proposed, many researchers have systematically analyzed the properties of its solution, including its global existence, asymptotic profile [11], global boundedness [12,13], finite-time blow-up [14,15,16], etc. Particularly, if the initial mass Ωu(x,y)dxdy of the cells in the 2D case satisfies a critical threshold value, its solution will blow up in finite time [7,8,9,14,15,16]. This blow-up denotes a mathematical concept of the bacterial aggregation arising in real biotic environments [7,17,18]. However, it is arduous if we want to obtain the analytical solutions of the model due to the strong nonlinear characteristics of the model itself. Meanwhile, it is still difficult to better numerically capture the blow-up or spike solutions. Therefore, it is desperately needed to establish a high accuracy and more stable numerical method to investigate the properties of the solutions for the chemotaxis system given by (1.1)–(1.3). At the same time, the numerical investigation of the chemotaxis system also facilitates the theoretical exploration of its dynamic behavior.

    In recent years, some numerical methods have been involved via the investigation of the chemotaxis systems [3,4,16,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35]. For example, Saito and Suzuki [19] proposed a conservative numerical scheme by using the FDM to solve a parabolic-elliptic coupling chemotaxis system. Meanwhile, in order to obtain a positivity-preserving scheme under total mass that is conservative, Saito [20] proposed a conservative scheme by using the FDM for the system given by (1.1)–(1.3) with nonlinear diffusion. Xiao et al. [21] derived a semi-implicit scheme by using a characteristic finite element method (FEM) to simulate the blow-up solutions of the chemotaxis system on surfaces, as well as pattern formulation and aggregation phenomenon of bacteria. And, their method has second-order accuracy in L2-norm and H1-norm errors. Epshteyn and Kurganov [22] introduced chemotaxis concentration gradient variables to rewrite the original Keller-Segel model given by (1.1)–(1.3), and they designed an internal penalty discontinuous Galerkin (DG) method for the rewritten system to border on the solution of the original system. Li et al. [23] used a local DG (i.e., LDG) method to modify that in [22], and they obtained the optimal convergence rates based on a specific finite element space before blow-up for the chemotaxis system. They also deduced a positivity-preserving P1 LDG scheme and proved that it is L1-stable. The LDG method was also used to solve the system given by (1.1)–(1.3) in [3], and the energy dissipation with the LDG discretization was proved. Sulman and Nguyen [24] proposed an adaptive moving mesh method by applying an implicit-explicit FEM to solve the system given by (1.1)–(1.3). The method has second-order accuracy in the spatiotemporal directions, and it is noteworthy that the obtained solutions are positive at all time steps if the initial values of the system (1.1) are positive. Qiu et al. [25] proposed a new scheme by using a interface-corrected direct DG method to solve this model (1.1), and their method satisfies the positivity-preserving requiremet without losing third-order accuracy. Based on the gradient flow structure, Shen and Xu [4] proposed a class of numerical schemes for solving the Keller-Segel model given by (1.1)–(1.3). Among them, the first-order accuracy scheme satisfies the mass conservation, bounded positivity, unique solvability and energy dissipation of the original differential equation, while the second-order accuracy scheme satisfies the first three properties. Chen et al. [16] analyzed the error of the numerical scheme proposed in [4], and they deduced the finite-time blow-up of non-radial numerical solutions under certain assumptions. Based on the generalized smoothed particle hydrodynamics meshless method, Dehghan and Abbaszadeh [26] established a second-order-accuracy numerical scheme for some chemotaxis models with the blow-up phenomena during tumor growth, and they numerically simulated the blow-up problems of the original chemotaxis model. Filbet [27] employed the finite volume method (FVM) to approximate the solution of the system given by (1.1)–(1.3), and they simulated the blow-up problems. Chertock and Kurganov [28] developed a center-upward scheme with second-order accuracy by using the FVM based on their previous work in [29] and the blow-up problems of the system (1.1)–(1.3), the pattern formation of bacteria were also simulated. Epshteyn [30] proposed an upwind-difference potentials scheme to solve the problems in complex geometries for the system given by (1.1)–(1.3). In addition, some other numerical methods [31,32,33,34,35], such as the fractional step (or operator splitting) method [31,32], hybrid finite-volume-finite-difference methods [33], generalized FDM [35], etc., were also used to solve the chemotaxis system given by (1.1)–(1.3).

    Although some numerical methods mentioned above for solving the system given by (1.1)–(1.3) can achieve high accuracy in the spatial direction, such as those in [25,30] (third-order accuracy) and [3,33] (fourth-order accuracy), most numerical methods have low accuracy, especially in the temporal direction; see [4,16,19,20,21,22,23,24,26,27,28,29,31,32,34,35]. Meanwhile, many high-accuracy numerical methods are non-compact in space, and the stability conditions are relatively harsh. In other words, if we want to employ these methods to solve real problems, we must take a small time step length to satisfy their stability conditions, which will expend expensive computational time. However, the high-order compact (HOC) difference scheme has attracted many researchers because of its strong advantages, such as fewer computational nodes, small computational errors, better numerical stability and non-complicated boundaries. Meanwhile, the backward differentiation formula with fourth-order accuracy (BDF-4), which is an A-stabilized method and appears first in [36], has been verified to be a feasible method to obtain high accuracy, and it has a relatively large stability condition range [37,38]. In addition, the major difficulty in solving the system comes from the nonlinearities, such as the chemotaxis term (χuv). One issue is that the coupling form will increase the discrete difficulty when we want to obtain the high accuracy and satisfy mass conservative schemes. Another difficulty is positivity preservation since u and v in the system (1.1)-(1.3) have real biological significance in complex chemotaxis phenomenon and cannot have negative values. To acquire the high-accuracy, positivity-preserving numerical solutions for the system given by (1.1)–(1.3), in this work, our purpose was to derive a compact difference scheme to approximate the solutions of the original chemotaxis system, and the scheme has space-time fourth-order accuracy and is stable, as well as positivity-preserving.

    The remainder of this paper is organized as follows. Some preliminary preparations with basic symbols, definitions and theorems are provided in Section 2. In Section 3, we deduce an HOC scheme for the system given by (1.1)–(1.3), and give the computational strategies for the initial time steps and the nonlinear terms. In Section 4, a time advancement algorithm combined with a multigrid method and a positivity-preserving algorithm are proposed. In Section 5, some numerical examples are employed to verify the accuracy, stability, positivity-preserving property, mass conservation and energy dissipation. The finitetime blow-up problems for the chemotaxis system given by (1.1)–(1.3) are simulated by using the proposed method. Finally the conclusion is provided in the end.

    First, the domain {(x,y,t)|ax,yb,0tT} is divided by uniform meshes N2×M, M,NZ+. Denote h=(ba)/N to represent the spatial step size, and τ=T/M stands for the temporal step length. Let xi=a+ih,yj=a+jh and tn=nτ, 0i,jN, 0nM, and mark the mesh points (xi,yj,tn). Figure 1 shows the 2D spatial mesh point stencil.

    Figure 1.  Spatial mesh point stencil.

    Second, we define Ωh={(xi,yj)|0i,jN}; its discrete boundary Ωh={(0,j),(N,j)|0jN}{(i,0),(i,N)|1iN1}; let Ωτ={tn|0nM} and Ωhτ=Ωh×Ωτ. For any mesh function wWhτ={wni,j|0i,jN,0nM} defined on Ωhτ, we have

    wn12i,j=12(wni,j+wn1i,j),δtwn12i,j=1τ(wni,jwn1i,j),δxwni,j=12h(wni+1,jwni1,j),δywni,j=12h(wni,j+1wni,j1),δ2xwni,j=1h2(wni+1,j2wni,j+wni1,j),δ2ywni,j=1h2(wni,j+12wni,j+wni,j1),Δtwni,j=112τ(25wni,j48wn1i,j+36wn2i,j16wn3i,j+3wn4i,j).

    Next, for simplicity, we let pwςp(x,y,t):≜wςp(x,y,t), pZ+, where ς represents x, y or t. And, denote A=[A1,A2,A3,A4,A5]=[25,48,36,16,3], B=[B1,B2,B3]=[1,10,1], S=[S1,S2,S3,S4]=[1,34,13,116]. To obtain a higher accuracy on the boundary described by Eq.(1.3), we employ the Taylor expansion with the Peano remainder to cope with Eq. (1.3). Thus, the following theorem holds.

    Theorem 2.1. Denote Wh={w(xi,yj)|0i,jN}; for mapping w:ΩhWh, we have the following:

    (1) If w(x,y)C5,0([x0,x4]×[y0,yN]), then wx(x0,yj)=112h5k=1Akw(xk1,yj)+O0,j(h4);

    (2) If w(x,y)C5,0([xN4,xN]×[y0,yN]), then wx(xN,yj)=112hN1k=N5ANkw(xk+1,yj)+ON,j(h4);

    (3) If w(x,y)C0,5([x0,xN]×[y0,y4]), then wy(xi,y0)=112h5k=1Akw(xi,yk1)+Oi,0(h4);

    (4) If w(x,y)C0,5([x0,xN]×[yN4,yN]), then wy(xi,yN)=112hN1k=N5ANkw(xi,yk+1)+Oi,N(h4), where the local truncation errors are

    O0,j(h4)=h464k=110k5Skwx5(x0+ksh,yj)(1s)4ds,0jN,ON,j(h4)=h464k=110k5Skwx5(xNksh,yj)(1s)4ds,0jN,Oi,0(h4)=h464k=110k5Skwy5(xi,y0+ksh)(1s)4ds,0iN,Oi,N(h4)=h464k=110k5Skwy5(xi,yNksh)(1s)4ds,0iN.

    Proof. First, we prove (1). We assume that w(x,y)Ck+1,0([xi1,xi+1]×[y0,yN]), and, according to the Taylor expansion with the Peano remainder, we have

    w(xi±h,yj)=kl=0(1)lhll!wxl(xi,yj)+hk+1k!10wxk+1(xi±sh,yj)(1s)kds, (2.1)

    where 1iN1 and 0jN. We suppose that w(x,y)C5,0([x0,x4]×[y0,yN]) and expand w(x1,yj), w(x2,yj), w(x3,yj) and w(x4,yj) at (x0,yj) by using Eq. (2.1) above, that is, we respectively take k=1,2,3,4, and have

    w(xk,yj)=4l=0(kh)ll!wxl(x0,yj)+(kh)52410wx5(x0+ksh,yj)(1s)4ds. (2.2)

    Then, we perform the following the operation: Eq. (2.2)k=1+α× Eq. (2.2)k=2+β× Eq. (2.2)k=3+γ× Eq. (2.2)k=4 for Eq. (2.2), and denote E=[E1,E2,E3,E4]=[1,α,β,γ], where α,β and γ, are constants; then, we can obtain

    4k=1Ekw(xk,yj)=5m=14k=1(kh)m1(m1)!Ekwxm1(x0,yj)+h5244k=1k5Ek10wx5(x0+ksh,yj)(1s)4ds. (2.3)

    Suppose that the coefficients of wx2(x0,yj), wx3(x0,yj) and wx4(x0,yj) in Eq. (2.3) equal to 0; we obtain

    {4k=1k2Ek=0,4k=1k3Ek=0,4k=1k4Ek=0.{α=34,β=13,γ=116.

    Substituting them into Eq. (2.3), denote A=[25,48,36,16,3] and S=[1,34,13,116], and we have

    wx(x0,yj)=112h5k=1Akw(xk1,yj)+O0,j(h4),0jN,

    where O0,j(h4)=h464k=110k5Skwx5(x0+ksh,yj)(1s)4ds. Similarly, (2) and (3) are also easily obtained. The proof is complete.

    According to Theorem 2.1, and by considering wx(x0,yj)=0, wx(xN,yj)=0, wy(xi,y0)=0 and wy(xi,yN)=0, we can easily derive the following fourth-order-accuracy boundary approximation formulas:

    w(x0,yj)125[48w(x1,yj)36w(x2,yj)+16w(x3,yj)3w(x4,yj)],w(xN,yj)125[48w(xN1,yj)36w(xN2,yj)+16w(xN3,yj)3w(xN4,yj)],w(xi,y0)125[48w(xi,y1)36w(xi,y2)+16w(xi,y3)3w(xi,y4)],w(xi,yN)125[48w(xi,yN1)36w(xi,yN2)+16w(xi,yN3)3w(xi,yN4)].

    In addition, we suppose that w(x,y,t)C6,6,5([xi1,xi+1]×[yj1,yj+1]×[tn,tn1]) for 1i,jN1, 1nM, and we denote I=[1,24,81,64] to derive the following truncation errors based on the Taylor expansion with the Peano remainder, that is,

    (Oxx)ni,j(h4)=h424102k=1[13(1s)315(1s)5]wx6(xi+(1)k1sh,yj,tn)ds,(Ot)ni,j(τ4)=τ4610(1μ)44k=1Ikwt5(xi,yj,tnkμτ)dμ,(Ox)ni,j(h4)=h412102k=1[13(1s)314(1s)4]wx5(xi+(1)k1sh,yj,tn)ds,(Ot)n12i,j(τ2)=τ21610(1μ)22k=1wt3(xi,yj,tn12+(1)k1μτ2)dμ,˜On12i,j(τ2)=τ2410(1μ)2k=1wt2(xi,yj,tn12+(1)k1μτ2)dμ.

    Similarly, we can easily derive the expressions of (Oyy)ni,j(h4), (Oy)ni,j(h4), (Ox)n1i,j(h4), (Oy)n1i,j(h4), (Oxx)n12i,j(h4) and (Oyy)n12i,j(h4). To facilitate mathematical analysis below, the following definitions are given.

    Definition 2.2. For any mesh function {wi,j|0i,jN}, the average operators A and B are defined as

    Awi,j={1254k=1Ak+1wk,j,i=0,0jN,1123k=1Bkwi+k2,j,1i,jN1,125N1k=N4ANk+1wk,j,i=N,0jN,Bwi,j={1254k=1Ak+1wi,k,j=0,0iN,1123k=1Bkwi,j+k2,1i,jN1,125N1k=N4ANk+1wi,k,j=N,0iN.

    Definition 2.3. The maximum norm and 2-norm errors are defined as

    ||||=max0i,jN|wMi,jw(xi,yj,tM)|,||||2=h2Ni,j=0[wMi,jw(xi,yj,tM)]2,

    where w(xi,yj,tM) and wMi,j stand for the exact and numerical solutions at the discrete mesh point (xi,yj,tM), respectively.

    Definition 2.4. Denote wmax(t):≜max(x,y)w(x,y,t); the relatively maximum norm and 2-norm errors in the temporal dimension are defined as

    Rel=wmax(t)wmax(t)||||wmax=max0nM|wmax(tn)wmax(tn)|max0nM|wmax(tn)|,Rel2=wmax(t)wmax(t)||2||wmax2=τMn=0[wmax(tn)wmax(tn)]2τMn=0[wmax(tn)]2,

    where wmax(t) represents the reference solution.

    Definition 2.5. The convergence rate is defined as

    Rate=log[Lν(h1)/Lν(h2)]log(h1/h2),

    where ν stands for or 2 and Lν(h1) and Lν(h2) are corresponding ||||ν norm errors which are closely related to h1 and h2, respectively.

    In this part, an HOC scheme is derived to border on the solution of the Keller-Segel system given by (1.1)–(1.3). To facilitate the numerical analysis, we rewrite an equivalent form for the system given by (1.1)–(1.3), that is,

    {Ut+Fx+Gy=DΔU+R,a<x,y<b,0<tT,(3.1)U(x,y,0)=U0(x,y),ax,yb,(3.2)Ux(a,y,t)=Ux(b,y,t)=0,0<tT,(3.3)Uy(x,a,t)=Uy(x,b,t)=0,0<tT,(3.4)

    where

    U:≜[uv],F:≜[χuvx0],G:≜[χuvy0],D:≜[d1],R:≜[0v+u],U0:≜[u0v0].

    Equation (3.1) is a nonlinear advection-diffusion-reaction equation. Then, we define

    U={[u(xi,yj,tn),v(xi,yj,tn)]:≜Uni,j|0i,jN,0nM},F={[χu(xi,yj,tn)vx(xi,yj,tn),0]:≜Fni,j|0i,jN,0nM},G={[χu(xi,yj,tn)vy(xi,yj,tn),0]:≜Gni,j|0i,jN,0nM},R={[0,u(xi,yj,tn)v(xi,yj,tn)]:≜Rni,j|0i,jN,0nM}

    on Ωhτ, respectively, and employ the FDM to discretize Eqs. (3.1)–(3.4) for the interior points.

    Now, we focus on Eq. (3.1) at (xi,yj,tn) for 1i,jN1 and 0<nM; we have

    Ut(xi,yj,tn)+Fx(xi,yj,tn)+Gy(xi,yj,tn)=D[Uxx(xi,yj,tn)+Uyy(xi,yj,tn)]+R(xi,yj,tn). (3.5)

    First, the following formulas are applied to discretize the second derivative terms on the right-hand side of Eq. (3.5), i.e.,

    Uxx(xi,yj,tn)=A1δ2xUni,j+A1(Oxx)ni,j(h4),Uyy(xi,yj,tn)=B1δ2yUni,j+B1(Oyy)ni,j(h4), (3.6)

    where 1i,jN1, and δ2x, and δ2y are the central difference operators. Then, for the nonlinear advection terms Fx(xi,yj,tn) and Gy(xi,yj,tn) on the left-hand side of Eq. (3.5), the following Padé compact schemes [39] are applied to compute them, that is,

    (1+h26δ2x)Fx(xi,yj,tn)=δxFni,j+(Ox)ni,j(h4),(1+h26δ2y)Gy(xi,yj,tn)=δyGni,j+(Oy)ni,j(h4), (3.7)

    where 1i,jN1 and 0nM. To be consistent with the accuracy of the interior points above, the following fourth-order boundary schemes [40] are employed, i.e.,

    Fx(x0,yj,tn)+1415Fx(x0+h,yj,tn)=1h(AF0)+(Ox)n0,j(h4),0jN,0<nM, (3.8)
    Fx(xN,yj,tn)1415Fx(xNh,yj,tn)=1h(BFN)+(Ox)nN,j(h4),0jN,0<nM, (3.9)
    Gy(xi,y0,tn)+1415Gy(xi,y0+h,tn)=1h(AG0)+(Oy)ni,0(h4),0iN,0<nM, (3.10)
    Gy(xi,yN,tn)1415Gy(xi,yNh,tn)=1h(BGN)+(Oy)ni,N(h4),0iN,0<nM, (3.11)

    where

    A=[18475,703180,8930,6730,7790,41300],B=[5225,1067180,6710,4110,13390,69300],F0=[Fn0,j,Fn1,j,Fn2,j,Fn3,j,Fn4,j,Fn5,j],FN=[FnN,j,FnN1,j,FnN2,j,FnN3,j,FnN4,j,FnN5,j],G0=[Gni,0,Gni,1,Gni,2,Gni,3,Gni,4,Gni,5],GN=[Gni,N,Gni,N1,Gni,N2,Gni,N3,Gni,N4,Gni,N5],(Ox)n0,j(h4)=160h4Fx5(x0,yj,tn)+1691800h5Fx6(x0,yj,tn)+O(h6),(Ox)nN,j(h4)=160h4Fx5(xN,yj,tn)2811800h5Fx6(xN,yj,tn)+O(h6),(Oy)ni,0(h4)=160h4Fy5(xi,y0,tn)+1691800h5Fy6(xi,y0,tn)+O(h6),(Oy)ni,N(h4)=160h4Fy5(xi,yN,tn)2811800h5Fy6(xi,yN,tn)+O(h6).

    Next, we employ the BDF-4 [36] to cope with the temporal derivative on the left-hand side of Eq. (3.5), i.e., for 1i,jN1 and 4nM,

    Ut(xi,yj,tn)=112τ(25Uni,j4k=1Ak+1Unki,j)+(Ot)ni,j(τ4):≜ΔtUni,j+(Ot)ni,j(τ4). (3.12)

    Substituting Eqs. (3.6)–(3.12) into Eq. (3.5), using the definitions above, combining Eqs. (3.2)–(3.4) and applying Theorem 2.1 for 4nM, we have

    {AB[ΔtUni,j+(Fx)ni,j+(Gy)ni,j]=D(Bδ2x+Aδ2y)Uni,j+ABRni,j+Qni,j,1i,jN1,(3.13)U0i,j=U0(xi,yj),                                                                            0i,jN,(3.14)Un0,j=AUn0,j+Qn0,j,          UnN,j=AUnN,j+QnN,j,                             0jN,(3.15)Uni,0=BUni,0+Qni,0,          Uni,N=BUni,N+Qni,N,                               0iN,(3.16)

    where Qni,j=AB(Ot)ni,j(τ4)+{[B(Oxx)ni,j+A(Oyy)ni,j]+AB[(Ox)ni,j+(Ox)n0,j+(Ox)nN,j+(Oy)ni,j+(Oy)ni,0+(Oy)ni,N]}(h4), Qn0,j=On0,j(h4), QnN,j=OnN,j(h4), Qni,0=Oni,0(h4) and Qni,N=Oni,N(h4). Then, there exist positive constants c1, c2, c3, c4 and c5; it holds that

    |Qni,j|{c1h4,i=0,0jN,4nM,c2h4,j=0,0iN,4nM,c3(τ4+h4),1i,jN1,4nM,c4h4,i=N,0jN,4nM,c5h4,j=N,0iN,4nM.

    Omitting Qni,j in Eqs. (3.13)–(3.16) and replacing U, Fx, Gy and R with U, Fx, Gy and R, respectively, the following HOC scheme for solving the chemotaxis system given by (1.1)–(1.3) for 4nM can be obtained:

    {AB[ΔtUni,j+(Fx)ni,j+(Gy)ni,j]=D(Bδ2x+Aδ2y)Uni,j+ABRni,j,   1i,jN1,(3.17)U0i,j=U0(xi,yj),                                                                                0i,jN,(3.18)Un0,j=AUn0,j,               UnN,j=AUnN,j,                                                 0jN,(3.19)Uni,0=BUni,0,               Uni,N=BUni,N,                                                      0iN.(3.20)

    Remark 3.1. For the function values containing the unknown time steps Uni,j in Eqs. (3.17)–(3.20), the following method is used to compute them, where k=1,2,3,, which stands for the nonlinear cyclic iterative parameter, Un,() denotes the approximate convergent solution at the nth time step and σ represents a small amount. That is, for the nth time step, the following is performed:

    Approximate Un,(k) using Un,(k1) and solve vn,(k)x and vn,(k)y using Eqs. (3.3), (3.4) and (3.7);

    Compute the chemotaxis terms Fn,(k)x and Gn,(k)y using Eqs. (3.7)–(3.11);

    Solve the linear system given by (3.17)–(3.20) for Un,(k);

    Set Un,()Un,(k); if Un,(k)Un,(k1)<σ, then proceed to the (n+1)th time step.

    Remark 3.2. Equations (3.17)–(3.20) demonstrate fully implicit compact scheme including five time steps, and its truncation error is O(τ4+h4), i.e., it has space-time fourth-order accuracy. It is noteworthy that the compactness here only refers to the spatial direction, since no more than nine mesh points are required for the 2D spatial mesh subdomain. However, it is non-compact in the temporal direction.

    Remark 3.3. For the HOC scheme given by (3.17)–(3.20), besides the values of U0 being known, the values of U1, U2 and U3 are also needed; then, we can use the scheme given by (3.17)–(3.20) to compute the values of Un(n4) in turn. Thus, in the following part, we discuss the computational strategies for solving U1, U2 and U3.

    Next, we consider Eq. (3.1) at (xi,yj,tn12) for 1i,jN1 and 1n3; we have

    Ut(xi,yj,tn12)+Fx(xi,yj,tn12)+Gy(xi,yj,tn12)=D[Uxx(xi,yj,tn12)+Uyy(xi,yj,tn12)]+R(xi,yj,tn12). (3.21)

    First, we employ the Crank-Nicolson (C-N) method to cope with the temporal derivative term in Eq. (3.21), that is,

    Ut(xi,yj,tn12)=δtUn12i,j+(Ot)n12i,j(τ2),1i,jN1,1n3.

    Then, we employ the following method for the remaining parts, i.e.,

    Θ(xi,yj,tn12)=12[Θ(xi,yj,tn)+Θ(xi,yj,tn1)]+(˜OΘ)n12i,j(τ2),1i,jN1,1n3,

    where Θ could be U,Fx,Gy and R in Eq. (3.21). Then, Eq. (3.21) becomes

    δtUn12i,j+12[Fx(xi,yj,tn)+Fx(xi,yj,tn1)+Gy(xi,yj,tn)+Gy(xi,yj,tn1)]=12D[Uxx(xi,yj,tn)+Uxx(xi,yj,tn1)+Uyy(xi,yj,tn)+Uyy(xi,yj,tn1)]+12[R(xi,yj,tn)+R(xi,yj,tn1)]+[(Ot)n12i,j+(˜Oxx)n12i,j+(˜Oyy)n12i,j+(˜OFx)n12i,j+(˜OGx)n12i,j+(˜OR)n12i,j](τ2),1i,jN1,1n3. (3.22)

    We use Eq. (3.6) to discretize Uxx(xi,yj,tn), Uxx(xi,yj,tn1), Uyy(xi,yj,tn) and Uyy(xi,yj,tn1) in the spatial direction of Eq. (3.22), and we use Eqs. (3.7)–(3.11) to compute Fx(xi,yj,tn), Fx(xi,yj,tn1), Gy(xi,yj,tn) and Gy(xi,yj,tn1) in the spatial direction. According to Eqs. (3.2)–(3.4) and Theorem 2.1, after rearrangement, the following form can be obtained:

    {[1τAB12D(Bδ2x+Aδ2y)]Uni,j+12AB[(Fx)ni,j+(Fx)n1i,j+(Gy)ni,j+(Gy)n1i,j]=[1τAB+12D(Bδ2x+Aδ2y)]Un1i,j+12AB(Rni,j+Rn1i,j)+Pni,j,1i,jN1,1n3,(3.23)U0i,j=U0(xi,yj),0i,jN,(3.24)Un0,j=AUn0,j+Pn0,j,UnN,j=AUnN,j+PnN,j,0jN,1n3,(3.25)Uni,0=BUni,0+Pni,0,Uni,N=BUni,N+Pni,N,0iN,1n3,(3.26)

    where Pni,j=AB[(Ot)n12i,j+(˜Oxx)n12i,j+(˜Oyy)n12i,j+(˜OFx)n12i,j+(˜OGx)n12i,j+(˜OR)n12i,j](τ2)+12{[B((Oxx)ni,j+(Oxx)n1i,j)+A(Oyy)ni,j+A(Oyy)n1i,j]+AB[(Ox)ni,j+(Ox)n0,j+(Ox)nN,j+(Oy)ni,j+(Oy)ni,0+(Oy)ni,N+(Ox)n1i,j+(Ox)n10,j+(Ox)n1N,j+(Oy)n1i,j+(Oy)n1i,0+(Oy)n1i,N]}(h4). Pn0,j=On0,j(h4), PnN,j=OnN,j(h4), Pni,0=Oni,0(h4), and Pni,N=Oni,N(h4). Then, there exist positive constants c6, c7, c8, c9 and c10 such that

    |Pni,j|{c6h4,i=0,0jN,1n3,c7h4),j=0,0iN,1n3,c8(τ2+h4),1i,jN1,1n3,c9h4,i=N,0jN,1n3,c10h4,j=N,0iN,1n3.

    Omitting Pni,j in Eqs. (3.23)–(3.26) and replacing U, Fx, Gy and R with U, Fx, Gy and R, respectively, we can get the following scheme for the initial time steps to solve the chemotaxis system given by (1.1)–(1.3) as follows:

    {[1τAB12D(Bδ2x+Aδ2y)]Uni,j+12AB[(Fx)ni,j+(Fx)n1i,j+(Gy)ni,j+(Gy)n1i,j]=[1τAB+12D(Bδ2x+Aδ2y)]Un1i,j+12AB(Rni,j+Rn1i,j),1i,jN1,1n3,(3.27)U0i,j=U0(xi,yj),0i,jN,(3.28)Un0,j=AUn0,j,UnN,j=AUnN,j,0jN,1n3,(3.29)Uni,0=BUni,0,Uni,N=BUni,N,0iN,1n3.(3.30)

    Remark 3.4. The same method in the previous section can be used to compute (Fx)ni,j, (Fx)n1i,j, (Gy)ni,j, (Gy)n1i,j, Rni,j and Rn1i,j in Eq. (3.27), as described in Remark 3.1.

    Remark 3.5. Equations (3.27)–(3.30) constitute a two-level implicit scheme, and its truncation error is O(τ2+h4).

    Remark 3.6. Since the C-N method is used to discretize the time derivative term, the scheme given by (3.27)–(3.30) for the initial time steps is unconditionally stable. On the other hand, the formula (3.12) is a k-step method for time integration. In this work, only the fourth-order BDF scheme (k=4) is used, which is widely employed to solve stiff differential equations [37]. Therefore, according to [36,37,38], the HOC scheme given by (3.17)–(3.20) is A-stable.

    As described in Subsections 3.1 and 3.2 above, we use the two-level implicit scheme given by (3.27)–(3.30) to compute U1i,j, U2i,j and U3i,j based on a known U0i,j. The Richardson extrapolation technique is applied to enhance the accuracy of the scheme given by (3.27)–(3.30) in time and ensure that it is consistent with that of the scheme given by (3.17)–(3.20); the following extrapolation formula is used, that is,

    ˆUni,j(h,τ)=13[4U2ni,j(h,τ/2)Uni,j(h,τ)],0i,jN,1n3. (3.31)

    Here, ˆUni,j(h,τ) stands for the extrapolated solution at the nth time step. Uni,j(h,τ) represents the computed solution obtained by using the scheme given by (3.27)–(3.30) with the temporal step length τ, whereas U2ni,j(h,τ/2) denotes the computed solution obtained by using the scheme given by (3.27)–(3.30) with τ/2.

    Remark 3.7. Equation (3.31) can extrapolate the accuracy in the temporal direction from second to fourth order.

    In this part, based on the compact difference schemes proposed above, we will establish the corresponding numerical algorithm. On one hand, because the slow convergence speed of the classical iterative method (for instance, Gauss–Seidel, Jacobi or successive over-relaxation iterations) used to solve the algebraic systems arise out of the fully implicit scheme at each time step, we aim to employ a multigrid method to accelerate the convergence speed. On the other hand, we want to establish a positivity-preserving algorithm to obtain the non-negative solutions of u and v without losing the fourth-order accuracy. Finally, we structure a time advancement algorithm based on these algorithms above to solve the system given by (1.1)–(1.3).

    Since the idea of the multigrid method was proposed in the 1930s, until Professor Brandt published his pioneering work [41] in 1977, this method was increasingly applied to solve engineering and technical problems. Up to now, the multigrid method has been widely applied in numerous disciplines and engineering fields, such as computational fluid dynamics and computational biology. The multigrid method has been theoretically proved to be a first-rank numerical computational method for linear elliptic problems [41,42]. Its convergence rate is independent of the mesh scale, and the computational speed of the existing computational program using the classical iterative method can be increased by employing the multigrid method with 1 to 2 orders of magnitude, which is especially suitable for application in large-scale engineering numerical computational problems. Nowadays, the multigrid method is also used more and more to solve non-elliptic problems to speed up the convergence of each time step and satisfy the needs of practical problems [43,44].

    The multigrid method is implemented by using a cycling algorithm. A multigrid cycle V includes three elements: relaxation, restriction and interpolation operators. We use alternating direction line Gauss–Seidel (ALGS) [45] iteration for the relaxation operator, as well as half-weighted and fully weighted restriction operators and a bilinear interpolation operator [42]. As the schemes given by (3.27)–(3.30) and (3.17)–(3.20) are both nonlinear, we employ a multigrid full approximation scheme (FAS) [42,44] to accommodate nonlinearities. For simplicity, we formally denote the algebraic equations arising from Eqs. (3.27)–(3.30) and (3.17)–(3.20) at each time step as

    LhUh=Fh. (4.1)

    The multigrid FAS algorithm is implemented recursively. Here, we only give the following two-level FAS V(ν1,ν2) cycle algorithm based on Eq. (4.1):

    Algorithm 4.1: Two-level FAS V(ν1,ν2) cycle algorithm
     Step 1:  Using ALGS iteration ν1 times at the fine mesh level to solve LhUh = Fh and computing
    its residual, i.e., Resh = FhLhUh;
     Step 2:  Restricting Uh and Resh to the next coarse mesh: ¯U2h I2hhUh, ¯Res2h ¯I2hhResh;
     Step 3: Using ALGS iteration ν1 times at the coarse mesh level to solve L2hU2h = L2h¯U2h+¯Res2h;
     Step 4:  Interpolating the correction errors from the coarse to the fine mesh levels: Δ Uh Ih2h(U2h¯U2h);
     Step 5:  Updating Uh at the fine mesh level: Uh Uh+Δ Uh;
     Step 6:  Using ALGS iteration ν2 times at the fine mesh level to solve LhUh = Fh.

    Remark 4.1. Here, Lh and L2h are the difference operators; I2hh (half-weighted) and ¯I2hh (fully weighted) are the restriction operators, which are used to restrict the approximate value and residual from the fine mesh level h to coarse mesh level 2h; and, Ih2h is the interpolation operator, which is used to transfer the error correction from the coarse to fine mesh levels. Resh is the residual at the fine mesh level, and Res2h is the residual at the adjacent coarse mesh level; ν1 and ν2 represent the pre-smooth and post-smooth numbers, respectively.

    Remark 4.2. In this work, we employ fully multi-level algorithms for all calculations. For instance, if the finest grid is 256×256, an eight-level algorithm is used. Under such conditions, the coarsest grid level only has 2×2 grids.

    When the schemes given by (3.27)–(3.30) and (3.17)–(3.20) are employed to approximate the solution of the chemotaxis system, it may generate negative values where blow-up occurs. To obtain the non-negative values of u and v for all time levels, we define the following average solution for each time step:

    ¯Uni,j=1ΔxiΔyjΩi,jU(x,y,tn)dxdy, (4.2)

    where Δxi=xi+1xi1, Δyj=yj+1yj1, 1i,jN1 and Ωi,j={[xi1,xi+1]×[yj1,yj+1]|1i,jN1}. Similar to the numerical integration formula of the one-variable function (i.e., Simpson formula), we structure the numerical integration of the two-variable function U(x,y) on Ωi,j for each time step, i.e., for 1i,jN1,

    Ωi,jU(x,y,tn)dxdy16ΔxiΔyj3k=1[U(xi+k2,yj,tn)+U(xi,yj+k2,tn)]. (4.3)

    Substituting Eq. (4.3) into Eq. (4.2) at each time step, we obtain the following average solutions ¯U, that is, for 1i,jN1 and 0<nM, we have

    ¯Uni,j163k=1[Uni+k2,j+Uni,j+k2]. (4.4)

    Next, a positivity-preserving limiter (PPL) [25,46] is employed to eliminate the negative values of Uni,j, 0i,jN1, 0<nM, that is,

    ˜Uni,j=min{|¯Uni,j¯Uni,jϖni,j|,1}(Uni,j¯Uni,j)+¯Uni,j,0i,jN,0<nM, (4.5)

    where ϖni,j=min{Uni,j,Uni1,j,Uni+1,j,Uni,j+1,Uni,j1},1i,jN1,0<nM. For 1i,jN1 and 0<nM, we can get the positivity-preserving solution ˜Uni,j0 without losing the fourth-order accuracy obtained by the proposed scheme. Finally, we establish the following algorithm for the nth time step:

    Algorithm 4.2: Positivity-preserving algorithm
     1.  Compute Uni, j using the schemes given by (3.27) - (3.30) and (3.17) - (3.20);
     2.  Compute the averaged solution ¯Uni, j using Eq. (4.4);
     3. Compute the positive solution ˜Uni, j using the PPL (4.5).

    On the basis of the derivation process described in Section 3 and the proposed algorithms in Subsections 4.1 and 4.2, we employ the difference schemes given by (3.27)–(3.30) and (3.17)–(3.20) to approximate the solution of the original system given by (1.1)–(1.3), which mainly includes the following three steps: first, provide the initial value of U0i,j; second, we employ Eqs. (3.7)–(3.11) and (3.27)–(3.31) to compute the values of U1i,j, U2i,j and U3i,j; finally, we employ Eqs. (3.7)–(3.11) and (3.17)–(3.20) to compute the value of Uni,j (n4). As (3.27)–(3.30) and (3.17)–(3.20) involve nonlinear parts, we establish the following algorithm with a nonlinear iteration tactic to solve them.

    Algorithm 4.3: Time advancement algorithm
     Establish the initial values: U0i,j=U0(xi,yj);
     for n=1,2,,M do 
       if 1n3 then 
        R2(n1)i,j=R2(n1),()i,j(h,τ/2),R(n1)i,j=R(n1),()i,j(h,τ);
        Computing the following chemotaxis parts at the (n1)th time step using Eqs. (3.7)(3.11);
        (Fx)2(n1)i,j=(Fx)2(n1),()i,j(h,τ/2),(Fx)(n1)i,j=(Fx)(n1),()i,j(h,τ);
        (Gy)2(n1)i,j=(Gy)2(n1),()i,j(h,τ/2),(Gy)(n1)i,j=(Gy)(n1),()i,j(h,τ);
     for k=1,2,3, do 
        Computing  Un0,j, UnN,j, Uni,0 and Uni,N using Eqs. (3.29)(3.30), i.e.,
        U2n,k0,j=AU2n,k10,j(h,τ/2),U2n,kN,j=AU2n,k1N,j(h,τ/2);Un,k0,j=AUn,k10,j(h,τ),Un,kN,j=AUn,k1N,j(h,τ);
        U2n,ki,0=AU2n,k1i,0(h,τ/2),U2n,ki,N=AU2n,k1i,N(h,τ/2);Un,ki,0=AUn,k1i,0(h,τ),Un,ki,N=AUn,k1i,N(h,τ);
        R2n,(k)i,j=R2n,(k1)i,j(h,τ/2),Rn,(k)i,j=Rn,(k1)i,j(h,τ);
        Computing the following chemotaxis parts at the (n)th time step using Eqs. (3.7)(3.11);
        (Fx)2n,(k)i,j=(Fx)2n,(k1)i,j(h,τ/2),(Fx)n,(k)i,j=(Fx)n,(k1)i,j(h,τ);
        (Gy)2n,(k)i,j=(Gy)2n,(k1)i,j(h,τ/2),(Gy)n,(k)i,j=(Gy)n,(k1)i,j(h,τ);
        Computing  U2n,(k)i,j(h,τ/2) and Un,(k)i,j(h,τ)  using Eqs.  (3.27)(3.30) and  Algorithm 4.1; 
         if ||U2n,(k)i,j(h,τ/2)U2n,(k1)i,j(h,τ/2)||<σand||Un,(k)i,j(h,τ)Un,(k1)i,j(h,τ)||<σ then 
          U2n,()i,j(h,τ/2)=U2n,(k)i,j(h,τ/2),Un,()i,j(h,τ)=Un,(k)i,j(h,τ);
         end if 
       end for 
        Computing the extrapolated solutions ˆUn, ()i, j using Eq.(3.31);
     else 
       for  k=1,2,3,, do 
        Computing  Un0,j, UnN,j, Uni,0 and Uni,N using Eqs. (3.15)(3.16) i.e., 
        Un,k0,j=AUn,k10,j(h,τ),Un,kN,j=AUn,k1N,j(h,τ);Un,ki,0=AUn,k1i,0(h,τ),Un,ki,N=AUn,k1i,N(h,τ);
        Rn,(k)i,j=Rn,(k1)i,j(h,τ);
        Computing the following chemotaxis parts at the  (n)th time step using Eqs. (3.7)(3.11);
        (Fx)n,(k)i,j=(Fx)n,(k1)i,j(h,τ),(Gy)n,(k)i,j=(Gy)n,(k1)i,j(h,τ);
        Computing   Un,(k)i,j(h,τ) using Eqs. (3.17)(3.20) and Algorithm 4.1;
         if ||Un,(k)i,j(h,τ)Un,(k1)i,j(h,τ)||<σ then Un,()i,j(h,τ)=Un,(k)i,j(h,τ);
         end if 
        end for 
       end if 
       Computing the positive solutions ˜Uni,j using Algorithm 4.2.
     end for 

    In this part, we give several numerical experiments to test the various properties of the proposed scheme and algorithm when solving the system given by (1.1)–(1.3). First, in Subsection 5.1 below, the accuracy and stability of the proposed scheme and algorithm in the absence of PPL and in the presence of a PPL when solving the chemotaxis system are tested. Second, in Subsection 5.2 below, we simulate the blow-up phenomena for the system given by (1.1)–(1.3), compare the obtained results from the proposed scheme and algorithm before and after using the positivity-preserving algorithm and verify the mass conservation and energy dissipation of the proposed method.

    We consider the following general type of Eq. (1.1) with source terms to test the accuracy, that is,

    {ut=u(uv)4u2+sin(x+y)+2u2[cos(2x+2y)2sin(x+y)][2+sin(x+y)]2,0x,y2π,t>0,vt=vv+u4v2+sin(x+y),0x,y2π,t>0.

    We use the periodic boundary conditions, and its analytical solution is u(x,y,t)=v(x,y,t)=(2+sin(x+y))e2t. By observing the expression of this analytical solution, we can see that the global solutions of this problem are positive values in the entire physical domain. Therefore, in this computation, we will use the difference scheme proposed in this work to approximate the solution for this problem in the absence of a PPL.

    In Table 1, we take τ=0.1h and test the convergence of the HOC scheme at the final time T=0.1 for Problem 5.1.1. From this table, the fourth-order accuracy is obtained by using the proposed scheme. The obtained results for the LDG method in [3] (the time step size τ=0.001h2 is used to solve this problem) represent third-order accuracy, which shows the advantage of our method, that is, we can obtain higher accuracy with a larger time step. Table 2 displays the convergence of the HOC scheme at T=1 and T=5 for Problem 5.1.1, where τ=0.5h. Table 3 shows the ||||2 errors of u(x,y,t) and v(x,y,t) for different step ratios λ=τ/h for Problem 5.1.1, where T=1 and τ=λh. From these tables, we obtain that the proposed method still converges at a fourth-order rate when computing long times and large step ratios, which indicates that our method has better stability.

    Table 1.  Convergence of the HOC scheme for Problem 5.1.1 at T=0.1.
    N LDG scheme [3] HOC scheme
    ||u|| Rate ||u||2 Rate ||u|| Rate ||u||2 Rate
    8 5.79E-03 2.53E-02 9.415E-03 3.974E-02
    16 8.79E-04 2.69 3.87E-03 2.71 3.359E-04 4.81 1.379E-03 4.85
    32 1.19E-04 2.91 5.18E-04 2.90 1.912E-05 4.13 7.711E-05 4.16
    64 1.49E-05 3.00 6.57E-05 2.98 1.112E-06 4.10 4.431E-06 4.12

     | Show Table
    DownLoad: CSV
    Table 2.  Convergence of the HOC scheme for Problem 5.1.1 with τ=0.5h.
    N T=1 T=5
    |||| Rate ||||2 Rate |||| Rate ||||2 Rate
    Cell density u(x,y,t)
    20 1.472E-03 6.477E-03 5.428E-06 2.415E-05
    40 1.046E-04 3.81 4.621E-04 3.81 2.289E-07 4.57 1.018E-06 4.57
    80 6.404E-06 4.03 2.832E-05 4.03 1.840E-08 3.64 5.259E-08 4.27
    160 3.902E-07 4.04 1.726E-06 4.04 6.967E-10 4.72 3.043E-09 4.11
    Chemoattractant concentration v(x,y,t)
    20 1.459E-03 6.488E-03 5.428E-06 2.415E-05
    40 1.041E-04 3.81 4.629E-04 3.81 2.289E-07 4.57 1.018E-06 4.57
    80 6.382E-06 4.03 2.837E-05 4.03 1.840E-08 3.64 5.259E-08 4.27
    160 3.889E-07 4.04 1.729E-06 4.04 6.967E-10 4.72 3.043E-09 4.11

     | Show Table
    DownLoad: CSV
    Table 3.  ||||2 errors of the HOC scheme for Problem 5.1.1 with T=1,τ=λh.
    N λ ||u||2 ||v||2
    0.4 2.8366512125E-05 2.8485324431E-05
    64 0.8 4.2917900191E-04 4.2928242249E-04
    1.6 5.7178346316E-03 5.7179298620E-03
    0.8 2.8033654005E-05 2.8041090035E-05
    128 1.6 4.1876675539E-04 4.1877318296E-04
    3.2 5.5523910128E-03 5.5523969211E-03
    1.6 2.7710430090E-05 2.7710912582E-05
    256 3.2 4.1352388720E-04 4.1352430255E-04
    6.4 5.4716502752E-03 5.4716506553E-03

     | Show Table
    DownLoad: CSV

    Next, we consider the following general type of Eq. (1.1) with two additional fluxes[24]:

    {ut=u(χuv)+r1(x,y,t),0x,y2π,t>0,vt=vv+u+r2(x,y,t),0x,y2π,t>0,

    where χ=1, and the homogeneous Neumann boundary condition (1.3) is applied.

    Case 1: The additional fluxes are r1(x,y,t)=2exp(2t)(cos2(x)+cos(x)cos(y)+cos2(y)1), r2(x,y,t)=0, with the initial condition u(x,y,0)=v(x,y,0)=cos(x)+cos(y),0x,y2π, and the analytical solution[23,24] u(x,y,t)=v(x,y,t)=exp(t)(cos(x)+cos(y)).

    Case 2: The additional fluxes are r1(x,y,t)=2e4tcos(2x+2y), r2(x,y,t)=0, with the initial condition u(x,y,0)=v(x,y,0)=cos(x+y),0x,y2π, and the analytical solution u(x,y,t)=v(x,y,t)=e2tcos(x+y).

    Table 4 shows the convergence of the HOC scheme when we take the time step size τ=0.001 at T=0.1 when using the proposed method to solve Problem 5.1.2 Case 1: the obtained results are compared with those in [24]. According to Table 4, the results obtained via the proposed method in the absence of a PPL can converge at the fourth-order rate. After using the positivity-preserving algorithm, although the errors obtained via the proposed scheme in the presence of a PPL decreases slightly, it still converges at the fourth-order rate and is better than the results in [24], as it is only second-order accuracy in [24]. In Tables 5 and 6, the results computed by using the proposed method for Problem 5.1.2 are given respectively. And, the convergence of the HOC scheme in the presence and absence of a PPL are compared when τ=0.5h, T=1 and T=5. By observing these tables, it can be found that the computed results for the proposed method can achieve the fourth-order accuracy in both cases, including with and without a PPL. Finally, we take the spatial meshes N=64,128 and 256, respectively, at the final time T=1. By using the proposed method to solve Problem 5.1.2 Case 1, we obtain the ||||2 norm errors computed with and without a PPL for different step ratios λ. The computed results are listed in Table 7. On the basis of this table, with the increase of the step ratio λ, the proposed method can converge well regardless of PPL existence, which also reflects that the proposed HOC difference scheme has better stability.

    Table 4.  Convergence of the HOC scheme for Problem 5.1.2 Case 1 with τ=0.001,T=0.1.
    Ref. [24] HOC scheme
    Without PPL With PPL
    N ||u||2 Rate ||u||2 Rate ||u||2 Rate
    10 1.02E-01 9.275E-03 2.170E-01
    20 2.74E-02 1.90 6.000E-04 3.95 2.870E-02 2.92
    40 6.30E-03 2.12 2.624E-05 4.52 2.469E-03 3.54
    80 1.45E-03 2.12 1.593E-06 4.04 1.476E-04 4.06

     | Show Table
    DownLoad: CSV
    Table 5.  Convergence of the HOC scheme for Problem 5.1.2 Case 1 when τ=0.5h,T=1.
    N Without PPL With PPL
    |||| Rate ||||2 Rate |||| Rate ||||2 Rate
    Cell density u(x,y,t)
    16 4.670E-03 1.022E-02 4.846E-03 1.424E-02
    32 1.009E-04 5.53 2.447E-04 5.38 2.684E-04 4.17 7.271E-04 4.29
    64 1.686E-06 5.90 5.902E-06 5.37 2.048E-05 3.71 4.974E-05 3.87
    128 7.683E-08 4.46 2.329E-07 4.66 1.508E-06 3.76 3.539E-06 3.81
    Chemoattractant concentration v(x,y,t)
    16 3.154E-03 8.025E-03 3.457E-03 1.217E-02
    32 9.167E-05 5.10 2.391E-04 5.07 2.769E-04 3.64 7.280E-04 4.06
    64 2.233E-06 5.36 6.687E-06 5.16 2.129E-05 3.70 5.074E-05 3.84
    128 4.873E-08 5.52 1.986E-07 5.07 1.558E-06 3.77 3.574E-06 3.83

     | Show Table
    DownLoad: CSV
    Table 6.  Convergence of the HOC scheme for Problem 5.1.2 Case 1 when τ=0.5h,T=5.
    N Without PPL With PPL
    ||u|| Rate ||u||2 Rate ||u|| Rate ||u||2 Rate
    16 2.693E-03 1.418E-02 4.424E-03 2.430E-02
    32 6.517E-05 5.37 3.530E-04 5.33 2.106E-04 4.39 1.172E-03 4.37
    64 1.631E-06 5.32 9.106E-06 5.28 1.331E-05 3.98 7.517E-05 3.96
    128 4.703E-08 5.12 2.672E-07 5.09 8.955E-07 3.89 5.098E-06 3.88

     | Show Table
    DownLoad: CSV
    Table 7.  ||||2 errors of the HOC scheme for Problem 5.1.2 Case 1 with T=1,τ=λh.
    N λ Without PPL With PPL
    ||u||2 ||v||2 ||u||2 ||v||2
    0.4 6.4023871303E-06 7.2616742789E-06 6.1464555035E-05 6.2581655559E-05
    64 0.8 1.2105436765E-05 1.1596839477E-05 3.5050251455E-05 3.5272598246E-05
    1.6 1.4457697771E-04 1.4260648214E-04 1.4620731799E-04 1.4406453705E-04
    0.8 9.4507876168E-07 8.4039918076E-07 2.4504429125E-06 2.4129793104E-06
    128 1.6 1.4859666027E-05 1.4234907841E-05 1.4967935575E-05 1.4314292650E-05
    3.2 1.5091392012E-04 1.4863144808E-04 1.5093408941E-04 1.4864245029E-04
    1.6 1.0545853150E-06 9.9391781457E-07 1.0618094946E-06 9.9866861508E-07
    256 3.2 1.5143113254E-05 1.4499368302E-05 1.5145042129E-05 1.4500157168E-05
    6.4 1.5233575478E-04 1.4973374665E-04 1.5233626588E-04 1.4973395499E-04

     | Show Table
    DownLoad: CSV

    Table 8 lists the numerical convergence of HOC scheme Problem 5.1.2 Case 2 in the absence and presence of a PPL, where T=0.2 and τ=0.1h, and it compares the computed results with those obtained via the method in [25]. The results in [25] have third-order accuracy before and after using a PPL, while the results computed via the method in this work converge at the fourth-order rate for both cases, including without and with a PPL, and the ||||, ||||2 errors are better than those in [25]. This also reflects the superiority of the HOC difference scheme proposed in this work.

    Table 8.  Convergence of the HOC scheme for Problem 5.1.2 Case 2 at T=0.2.
    N Ref. [25] HOC scheme
    ||u|| Rate ||u||2 Rate ||u|| Rate ||u||2 Rate
    Without PPL
    20 3.65E-03 3.04E-03 6.189E-04 9.217E-04
    40 4.55E-04 3.00 3.67E-04 3.05 3.554E-05 4.12 7.089E-05 3.70
    80 5.68E-05 3.00 4.58E-05 3.00 2.535E-06 3.81 4.677E-06 3.92
    160 7.02E-06 3.01 5.77E-06 2.99 1.741E-07 3.86 3.077E-07 3.93
    With PPL
    20 3.55E-02 3.62E-02 6.337E-04 1.019E-03
    40 4.37E-03 3.02 4.54E-03 3.00 3.788E-05 4.06 7.260E-05 3.81
    80 5.40E-04 3.02 5.61E-04 3.02 2.562E-06 3.89 4.706E-06 3.95
    160 6.74E-05 3.01 6.92E-05 3.02 1.745E-07 3.88 3.082E-07 3.93

     | Show Table
    DownLoad: CSV

    In the 2D case of the system given by (1.1)–(1.3), when the initial mass satisfies certain critical values, the solutions will blow up in a finite time. To capture the blow-up time of the solutions, referring to [47], the following adaptive technology is used to obtain the optimal time step sizes, i.e.,

    τn=hUn1,n0. (5.1)

    The computation is terminated if ||U||10ν, where ν denotes the maximum order of magnitude of cell density when the actual problem experiences a blow up. And, tn=n1m=0τm is regarded as the approximation of the blow-up time T.

    In a rectangular region [0.5,0.5]×[0.5,0.5], we first take χ=d=1 and test the initial-boundary value problem (IBVP) [23,24,25] for the chemotaxis system given by (1.1)–(1.3) with the homogeneous Neumann boundary conditions given by (1.3), where its initial conditions are as follow:

    {u(x,y,0)=840e84(x2+y2),0.5x,y0.5,t>0,v(x,y,0)=420e42(x2+y2),0.5x,y0.5,t>0.

    The solution of this problem will blow up in finite time at the center of this rectangular region.

    (1) Finite-time blow up. Figure 2 plots the numerical solutions of cell density u(x,y,t) as obtained by using the proposed method in the presence of a PPL, to solve the IBVP 5.2.1 at t=0, t=5×105, t=104 and t=1.2×104, where the initial time step size τ0=5×106 and the fixed spatial mesh numbers are 100×100. In [22], the blow-up time t is approximately equal to 1.21×104, which is verified by the authors. Based on observation, we can expect to see that Problem 5.2.1 displays blow-up phenomena around t1.2×104 at the center of the rectangular region [0.5,0.5]×[0.5,0.5]. It shows that the maximum peak values of the cell density u(x,y,t) gradually increase over time and present a very sharp aiguille structure at the central zone. We find that u(x,y,t) is strictly positive during time evolution by using the proposed method in the presence of a PPL. Figure 3 displays the projections of u(x,y,t) on the xou plane for Problem 5.2.1 at two pre-blow-up times t=104 and t=1.2×104, with τ0=5×106 and N=100. They intuitively show that the proposed method in this work can achieve positivity-preserving capability.

    Figure 2.  Cell density u(x,y,t) with PPL for Problem 5.2.1.
    Figure 3.  Projection for u(x,y,t) on the xou plane for Problem 5.2.1.

    (2) Non-negativity. To highlight the performance of positivity-preserving capability for our method in the presence of a PPL, we employ the proposed scheme in the absence and presence of a PPL to solve the IBVP 5.2.1 in two different space grids N=50 and N=100 at pre-blow-up times t=104 and t=1.2×104, where τ0=5×106. In Figure 4, we can see that u(x,y,t) has negative values in the absence of a PPL, and that all negative values are eliminated after the PPL is used. At the same time, the negative values decrease with the increase of spatial grid numbers, and the solutions become smoother. However, for the same number of grids, the oscillation of the blow-up area increases with time.

    Figure 4.  One-dimensional profile of u(x,y,t) along y=0 for Problem 5.2.1.

    (3) Influence of chemotaxis sensitivity coefficient χ. We employ the proposed method to compute the blow-up solution in a finite time for this IBVP 5.2.1 under different chemotaxis sensitivity coefficients χ=2 and χ=3 at the same pre-blow-up time t=5×105. The results are shown in Figure 5, where τ0=106 and N=100. We find that, for different values of χ, the cell density u(x,y,t) achieves negative values in the absence of a PPL. After using the PPL, the algorithm effectively eliminates the negative values. At the same time, the maximum peak value when χ=3 is one order of magnitude higher than that when χ=2. It can be seen that the peak value is gradually increased over χ. Therefore, the sensitivity intensity χ also affects the value of u(x,y,t) and the blow-up degree of Problem 5.2.1.

    Figure 5.  One-dimensional profile of u(x,y,t) along y=0 for Problem 5.2.1 with t=5×105.

    (4) Chemoattractant concentration. In Figure 6, we plot the numerical solutions of chemoattractant concentration v(x,y,t) for Problem 5.2.1 when t=104, τ0=5×106 and N=100, and we have compared the computed results in the absence and presence of a PPL. Observing Figures 4b and 6, we find that, although we applied the same parameters in the absence of a PPL, u(x,y,t) achieves negative values, but v(x,y,t) does not have negative values, and that the numerical solutions in the presence of a PPL is highly consistent with those in the absence of a PPL. In summary, our method can effectively eliminate these negative values in this computation. It is consistent with the results obtained in the absence of a PPL if the obtained solutions do not have negative values. Meanwhile, Problem 5.2.1 always experiences a finite-time blow-up at the center of the rectangular region under this initial condition.

    Figure 6.  Chemoattractant concentration v(x,y,t) for Problem 5.2.1 with t=104.

    (5) Accuracy and mass conservation. Since the analytical solution for this IBVP 5.2.1 is very difficult to obtain, to verify the accuracy and mass conservation, we take the fine mesh 800×800 as the reference solution. The convergence of the HOC scheme in the presence of a PPL are tested by using the relative errors between the maximum value Umax(t) of u(x,y,t) and the reference solutions Umax(t); the results are shown in Table 9. We find that u(x,y,t) can converge at a fourth-order relative error convergence rate. Meanwhile, we plot umax(t), which is the time evolution of the maximum value of u(x,y,t), on the left of Figure 7, and the numerical energy with time on the right of Figure 7, using the proposed scheme in the presence of a PPL to solve this IBVP 5.2.1 with PPL at τ0=5×106, N=50,80, 100. We find that the energy is strictly positive and monotonically decreasing with the evolution of time, which verifies the energy dissipation nature of the original problem. In addition, in Table 10, we give the discrete masses of u(x,y,t) and v(x,y,t) as obtained by using the proposed scheme in the absence and presence of a PPL to solve this IBVP 5.2.1 at different times T with τ0=5×106,N=100. According to this table, we can find that the proposed scheme in the absence of a PPL well verifies the mass conservation. After using the PPL, it is still conserved before blow up, and the mass conservation is slightly affected near the blow-up time.

    Table 9.  Relative ||||, ||||2 errors and convergence rates of umax(t) for Problem 5.2.1 with τ0=106,t=5×105.
    N Rel Rate Rel2 Rate
    100 3.204E-02 2.104E-02
    200 1.864E-03 4.10 1.197E-03 4.14
    300 3.290E-04 4.28 2.134E-04 4.25
    400 7.872E-05 4.97 5.902E-05 4.47
    500 3.332E-05 3.85 2.154E-05 4.52

     | Show Table
    DownLoad: CSV
    Table 10.  Massu and Massv of the HOC scheme for Problem 5.2.1 at time T with τ0=5×106,N=100.
    T Without PPL With PPL
    Massu Massv Massu Massv
    0 31.4159265323 31.4156968041 31.4159265323 31.4156968041
    1.0×105 31.4159265336 31.4157188735 31.4159265336 31.4157188735
    2.0×105 31.4159265339 31.4157222761 31.4159265339 31.4157222761
    3.0×105 31.4159265341 31.4157255274 31.4159265341 31.4157255274
    4.0×105 31.4159265344 31.4157287900 31.4159265344 31.4157287900
    5.0×105 31.4159265348 31.4157319661 31.4159265348 31.4157319661
    7.0×105 31.4159265354 31.4157381680 31.4292462953 31.4157381680
    1.0×104 31.4159265366 31.4157472806 37.7315206418 31.4157952866
    1.2×104 31.4159265376 31.4157533114 48.2548641463 31.4159976499

     | Show Table
    DownLoad: CSV
    Figure 7.  Left: Plots of umax(t): the time evolution of the maximum value for u(x,y,t); Right: Plots of the free energy with time for Problem 5.2.1 with the PPL at N=50,80 and 100.

    In the section, we consider the following IBVP [24] for the chemotaxis system given by (1.1)–(1.3) with the boundary conditions given by (1.3) for the rectangular region [0.5,0.5]×[0.5,0.5]. Its initial condition is given as follows:

    u(x,y,0)=1000e100[(x0.15)2+(y0.15)2],v(x,y,0)=0.

    In the computation, we take χ=d=1 and employ the proposed scheme in the presence of a PPL to simulate the temporal evolution of u(x,y,t) and v(x,y,t). The solution for this Problem 5.2.2 will blow up in finite time [8,9,24] at the corner of the rectangular region [0.5,0.5]×[0.5,0.5]. Next, we take the initial time step size τ0=103 and fix a mesh of 200×200 to compute them.

    Figures 8 and 9 respectively display the evolutionary processes of u(x,y,t) and v(x,y,t) over time t. It can be seen in these figures that the maximum value of u(x,y,t) gradually moves toward the corner of the rectangular region [0.5,0.5]×[0.5,0.5] in the form of a Gaussian profile. When it is close to the corner of the rectangular region, the value of u(x,y,t) increases rapidly, finally blowing up in a finite time. At the same time, v(x,y,t) also increases rapidly at the corner of the rectangular region.

    Figure 8.  Evolution of u(x,y,t) over time t for Problem 5.2.2.
    Figure 9.  Evolution of v(x,y,t) over time t for Problem 5.2.2.

    In this work, first, a fourth-order compact difference scheme for solving a 2D Keller-Segel model was derived and the computing strategies for the initial time steps and the nonlinear chemotaxis terms has been presented. Then, a multigrid algorithm was established to improve the computational efficiency and the binary function integration method has been used to establish the positivity-preserving algorithm. Third, the accuracy and reliability of the proposed method in the absence and presence of a PPL have been respectively verified by several numerical examples with flux source terms. And, we simulated the blow-up phenomena for the chemotaxis system in the center and corner of the rectangular region, respectively. Finally, the energy dissipation and mass conservation were also verified by numerical experiments. In summary, compared with the classical high-order numerical methods for solving the chemotaxis system in the literature, this proposed method has the following advantages:

    ● The proposed HOC scheme is compact in the spatial directions and fully implicit, as no more than nine mesh points are required for the 2D spatial mesh subdomain.

    ● The truncation error of the proposed HOC scheme is O(τ4+h4), i.e., it has space-time fourth-order accuracy. Thus, we can take large time step sizes to approximate the solution of this chemotaxis system.

    ● The proposed numerical algorithm can effectively filter the negative values in the solution process to ensure that the values of cell density are not negative at all time steps, while maintaining the original fourth-order accuracy without loss.

    ● The computed results show that the proposed method is more accurate than most such schemes reported in the literature.

    Nevertheless, our method also has some disadvantages. First, the new scheme has five layers in time. When it is used to compute the actual problem, three start-up time steps are required, which is slightly more troublesome for programming. Second, the proposed positivity-preserving algorithm would slightly affect its mass conservation when the actual problem blows up. The above deficiencies are also the driving force for our next work. At the same time, it is also possible to extend our approach to solve more complex chemotaxis and haptotaxis systems, which will be incorporated into our next work.

    Support of the study was received from the National Natural Science Foundation of China (12161067, 12001015, 12261067), National Natural Science Foundation of Ningxia, China (2022AAC02023) and National Youth Top-Notch Talent Support Program of Ningxia, China.

    All authors assert no conflict of interest.

    The data that support the findings of this study are available from the authors upon reasonable request.



    [1] Gallardo J, Pérez-Illana M, Martín-González N, et al. (2021) Adenovirus structure: what is new?. In J Mol Sci 22: 5240. https://doi.org/10.3390/ijms22105240
    [2] Ebner K, Pinsker W, Lion T (2005) Comparative sequence analysis of the hexon gene in the entire spectrum of human adenovirus serotypes: phylogenetic, taxonomic, and clinical implications. J Virol 79: 12635-12642. https://doi.org/10.1128/JVI.79.20.12635-12642.2005
    [3] Suraka B, Usman U, Tijjani A (2022) A brief review on the molecular biology of human adenoviruses. Baghdad J Biochem Appl Biol Sci 3: 166-182. https://doi.org/10.47419/bjbabs.v3i03.146
    [4] Rowe WP, Huebner RJ, Gilmore LK, et al. (1953) Isolation of a cytopathogenic agent from human adenoids undergoing spontaneous degeneration in tissue culture. Proc Soc Ex Biol Med 84: 570-573. https://doi.org/10.3181/00379727-84-20714
    [5] Lion T (2014) Adenovirus infections in immunocompetent and immunocompromised patients. Clin Microbiol Rev 27: 441-462. https://doi.org/10.1128/CMR.00116-13
    [6] Hur SJ, Kim DH, Chun SC, et al. (2013) Effect of adenovirus and influenza virus infection on obesity. Life Sci 93: 531-535. https://doi.org/10.1016/j.lfs.2013.08.016
    [7] Mullis K, Faloona F, Scharf S, et al. (1986) Specific enzymatic amplification of DNA in vitro: the polymerase chain reaction. Cold Spring Harb Symp Quant Biol 51: 263-273. https://doi.org/10.1101/sqb.1986.051.01.032
    [8] Wu X, Zhang J, Lan W, et al. (2022) Molecular typing and rapid identification of human adenoviruses associated with respiratory diseases using universal PCR and sequencing primers for the three major capsid genes: penton base, hexon, and fiber. Front Microbiol 13: 911694. https://doi.org/10.3389/fmicb.2022.911694
    [9] Shieh WJ (2022) Human adenovirus infections in pediatric population–an update on clinico-pathologic correlation. Biomed J 45: 38-49. https://doi.org/10.1016/j.bj.2021.08.009to
    [10] Pehler-Harrington K, Khanna M, Waters CR, et al. (2004) Rapid detection and identification of human adenovirus species by adenoplex, a multiplex PCR-enzyme hybridization assay. J Clin Microbiol 42: 4072-4076. https://doi.org/10.1128/JCM.42.9.4072-4076.2004
    [11] Al-Siyabi T, Binkhamis K, Wilcox M, et al. (2013) A cost effective real-time PCR for the detection of adenovirus from viral swabs. Virol J 10: 184. https://doi.org/10.1186/1743-422X-10-184
    [12] Casas I, Avellon A, Mosquera M, et al. (2005) Molecular identification of adenoviruses in clinical samples by analyzing a partial hexon genomic region. J Clin Microbiol 43: 6176-6182. https://doi.org/10.1128/JCM.43.12.6176-6182.2005
    [13] He JW, Jiang S (2005) Quantification of enterococci and human adenoviruses in environmental samples by real-time PCR. Appl Environ Microbio 71: 2250-2255. https://doi.org/10.1128/AEM.71.5.2250-2255.2005
    [14] Dolskiy AA, Grishchenko IV, Yudkin DV (2020) Cell cultures for virology: usability, advantages, and prospects. Int J Mo Sci 21: 7978. https://doi.org/10.3390/ijms21217978
    [15] Shike H, Shimizu C, Kanegaye J, et al. (2005) Quantitation of adenovirus genome during acute infection in normal children. Pediatr Infect Dis J 24: 29-33. https://doi.org/10.1097/01.inf.0000148882.97399.79
    [16] Choi S, Jiang SC (2005) Real-time PCR quantification of human adenoviruses in urban rivers indicates genome prevalence but low infectivity. Appl Enviro Microb 71: 7426-7433. https://doi.org/10.1128/AEM.71.11.7426-7433.2005
    [17] Morozumi M, Shimizu H, Matsushima Y, et al. (2014) Evaluation of new immunochromatographic assay kit for adenovirus detection in throat swab: Comparison with culture and real-time PCR results. J Infect Chemother 20: 303-306. https://doi.org/10.1016/j.jiac.2014.01.005
    [18] Jones MS, Hudson NR, Gibbins C, et al. (2011) Evaluation of type-specific real-time PCR assays using the LightCycler and J.B.A.I.D.S. for detection of adenoviruses in species HAdV-C. PloS One 6: e26862. https://doi.org/10.1371/journal.pone.0026862
    [19] Feeney SA, Armstrong VJ, Mitchell SJ, et al. (2011) Development and clinical validation of multiplex TaqMan® assays for rapid diagnosis of viral gastroenteritis. J Med Virol 83: 1650-1656. https://doi.org/10.1002/jmv.22162
    [20] Li F, Zhao LQ, Deng J, et al. (2013) Detecting human adenoviruses in respiratory samples collected from children with acute respiratory infections by loop-mediated isothermal amplification. Chinese J Pediatr 51: 52-57.
    [21] Metzger-Boddien C, Kehle J (2005) Development and evaluation of a sensitive PCR-ELISA for detection of adenoviruses in feces. Intervirology 48: 297-300. https://doi.org/10.1159/000085098
    [22] Murali S, Langston AA, Nolte FS, et al. (2009) Detection of respiratory viruses with a multiplex polymerase chain reaction assay (multicode-PLx respiratory virus panel) in patients with hematologic malignancies. Leukemia & lymphoma 50: 619-624. https://doi.org/10.1080/10428190902777665
    [23] van Maarseveen NM, Wessels E, de Brouwer CS, et al. (2010) Diagnosis of viral gastroenteritis by simultaneous detection of adenovirus group F, astrovirus, rotavirus group A, norovirus genogroups I and II, and sapovirus in two internally controlled multiplex real-time PCR assays. J Clin Microbiol J Clin Microbiol 49: 205-210. https://doi.org/10.1016/j.jcv.2010.07.019
    [24] Lalli MA, Langmade JS, Chen X, et al. (2021) Rapid and extraction-free detection of SARS-CoV-2 from saliva by colorimetric reverse-transcription loop-mediated isothermal amplification. Clin Chem 67: 415-424. https://doi.org/10.1093/clinchem/hvaa267
    [25] Foo PC, Nurul Najian AB, Muhamad NA, et al. (2020) Loop-mediated isothermal amplification (LAMP) reaction as viable PCR substitute for diagnostic applications: a comparative analysis study of LAMP, conventional PCR, nested PCR (nPCR) and real-time PCR (qPCR) based on Entamoeba histolytica DNA derived from faecal sample. BMC Biotechnol 20–34. https://doi.org/10.1186/s12896-020-00629-8
    [26] Wang Y, Kong F, Yang Y, et al. (2008) A multiplex PCR-based reverse line blot hybridization (mPCR/RLB) assay for detection of bacterial respiratory pathogens in children with pneumonia. Pediatr Pulm 43: 150-159. https://doi.org/10.1002/ppul.20749
    [27] O'Sullivan MV, Zhou F, Sintchenko V, et al. (2011) Multiplex PCR and reverse line blot hybridization assay (mPCR/RLB). J Vis Exp 54: e2781. https://doi.org/10.3791/2781
    [28] Pierce VM, Elkan M, Leet M, et al. (2012) Comparison of the Idaho technology film array system to real-time PCR for detection of respiratory pathogens in children. J Clin Microbiol 50: 364-371. https://doi.org/10.1128/JCM.05996-11
    [29] Zink S, Cirino NM, Egan C (2011) Multiplex PCR product detection and discrimination. Mol Microbiol: Diagn Prin Pract 325–341. https://doi.org/10.1128/9781555816834.ch21
    [30] Silbereisen A Monoclonal antibodies as tools in antigen detection assay and vaccine development: design of a sensitive detection test for Brucella bacteria and profiling of the malaria vaccine candidate antigen reticulocyte-binding homolog 2 (PfRH2), University_of_Basel (2015).
    [31] Puppe W, Weigl JA, Aron G, et al. (2004) Evaluation of a multiplex reverse transcriptase PCR ELISA for the detection of nine respiratory tract pathogens. J Clin Virol 30: 165-174. https://doi.org/10.1016/j.jcv.2003.10.003
    [32] Puppe W, Weigl J, Gröndahl B, et al. (2013) Validation of a multiplex reverse transcriptase PCR ELISA for the detection of 19 respiratory tract pathogens. Infection 41: 77-91. https://doi.org/10.1007/s15010-012-0298-6
    [33] Zulauf BJ Multiplex real-time PCR in the detection and differentiation of bovine respiratory disease pathogens (2007). Available from: https://ir.library.oregonstate.edu/concern/graduate_thesis_or_dissertations/0g354j20f
    [34] Zhang Y, Cao L, Xu Z, et al. (2020) Evaluation of a multiplex PCR assay for detection of respiratory viruses and Mycoplasma pneumoniae in oropharyngeal swab samples from outpatients. J Clin Lab Anal 34: e23032. https://doi.org/10.1002/jcla.23032
    [35] Lee JA, Kim NH, Kim SJ, et al. (2005) Rapid identification of human adenovirus types 3 and 7 from respiratory specimens via multiplex type-specific PCR. J Clin Microbiol 43: 5509-5514. https://doi.org/10.1128/JCM.43.11.5509-5514.2005
    [36] Chin EL, da Silva C, Hegde M (2013) Assessment of clinical analytical sensitivity and specificity of next-generation sequencing for detection of simple and complex mutations. BMC Genet 14: 6. https://doi.org/10.1186/1471-2156-14-6
    [37] Reijans M, Dingemans G, Klaassen CH, et al. (2008) RespiFinder: a new multiparameter test to differentially identify fifteen respiratory viruses. J Clin Microbiol 46: 1232-1240. https://doi.org/10.1128/JCM.02294-07
    [38] Thevendran R, Citartan M (2022) Assays to estimate the binding affinity of aptamers. Talanta 238: 122971. https://doi.org/10.1016/j.talanta.2021.122971
    [39] Shaffer SM, Joshi RP, Chambers BS, et al. (2015) Multiplexed detection of viral infections using rapid in situ RNA analysis on a chip. Lab Chip 15: 3170-3182. https://doi.org/10.1039/c5lc00459d
    [40] Renois F, Talmud D, Huguenin A, et al. (2010) Rapid detection of respiratory tract viral infections and coinfections in patients with influenza-like illnesses by use of reverse transcription-PCR DNA microarray systems. J Clin Microbiol 48: 3836-3842. https://doi.org/10.1128/JCM.00733-10
    [41] Daly GM, Leggett RM, Rowe W, et al. (2015) Host subtraction, filtering and assembly validations for novel viral discovery using next generation sequencing data. PloS one 10: e0129059. https://doi.org/10.1371/journal.pone.0129059
    [42] Fujii H, Kakiuchi S, Tsuji M, et al. (2018) Application of next-generation sequencing to detect acyclovir-resistant herpes simplex virus type 1 variants at low frequency in thymidine kinase gene of the isolates recovered from patients with hematopoietic stem cell transplantation. J Virol Methods 251: 123-128. https://doi.org/10.1016/j.jviromet.2017.10.019
    [43] Zhang W, Ehrhardt A (2017) Getting genetic access to natural adenovirus genomes to explore vector diversity. Virus genes 53: 675-683. https://doi.org/10.1007/s11262-017-1487-2
    [44] Capobianchi MR, Giombini E, Rozera G (2013) Next-generation sequencing technology in clinical virology. Clin Microbiol Infec 19: 15-22. https://doi.org/10.1111/1469-0691.12056
    [45] Reta DH, Tessema TS, Ashenef AS, et al. (2020) Molecular and immunological ddiagnostic ttechniques of medical viruses. Int J Microbiol 2020: 8832728. https://doi.org/10.1155/2020/8832728
    [46] Salez N, Vabret A, Leruez-Ville M, et al. (2015) Evaluation of four commercial multiplex molecular tests for the diagnosis of acute respiratory infections. PloS One 10: e0130378. https://doi.org/10.1371/journal.pone.0130378
  • Reader Comments
  • © 2023 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(4635) PDF downloads(251) Cited by(1)

Figures and Tables

Figures(9)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog