Research article Special Issues

Viral infection dynamics with mitosis, intracellular delays and immune response


  • In this paper, we propose a delayed viral infection model with mitosis of uninfected target cells, two infection modes (virus-to-cell transmission and cell-to-cell transmission), and immune response. The model involves intracellular delays during the processes of viral infection, viral production, and CTLs recruitment. We verify that the threshold dynamics are determined by the basic reproduction number R0 for infection and the basic reproduction number RIM for immune response. The model dynamics become very rich when RIM>1. In this case, we use the CTLs recruitment delay τ3 as the bifurcation parameter to obtain stability switches on the positive equilibrium and global Hopf bifurcation diagrams for the model system. This allows us to show that τ3 can lead to multiple stability switches, the coexistence of multiple stable periodic solutions, and even chaos. A brief simulation of two-parameter bifurcation analysis indicates that both the CTLs recruitment delay τ3 and the mitosis rate r have a strong impact on the viral dynamics, but they do behave differently.

    Citation: Jiawei Deng, Ping Jiang, Hongying Shu. Viral infection dynamics with mitosis, intracellular delays and immune response[J]. Mathematical Biosciences and Engineering, 2023, 20(2): 2937-2963. doi: 10.3934/mbe.2023139

    Related Papers:

    [1] Li Cai, Yu Hao, Pengfei Ma, Guangyu Zhu, Xiaoyu Luo, Hao Gao . Fluid-structure interaction simulation of calcified aortic valve stenosis. Mathematical Biosciences and Engineering, 2022, 19(12): 13172-13192. doi: 10.3934/mbe.2022616
    [2] Rongxing Qin, Lijuan Huang, Wei Xu, Qingchun Qin, Xiaojun Liang, Xinyu Lai, Xiaoying Huang, Minshan Xie, Li Chen . Identification of disulfidptosis-related genes and analysis of immune infiltration characteristics in ischemic strokes. Mathematical Biosciences and Engineering, 2023, 20(10): 18939-18959. doi: 10.3934/mbe.2023838
    [3] Yuedan Wang, Jinke Huang, Jiaqi Zhang, Fengyun Wang, Xudong Tang . Identifying biomarkers associated with the diagnosis of ulcerative colitis via bioinformatics and machine learning. Mathematical Biosciences and Engineering, 2023, 20(6): 10741-10756. doi: 10.3934/mbe.2023476
    [4] Xuesi Chen, Qijun Zhang, Qin Zhang . Predicting potential biomarkers and immune infiltration characteristics in heart failure. Mathematical Biosciences and Engineering, 2022, 19(9): 8671-8688. doi: 10.3934/mbe.2022402
    [5] Jun Wang, Mingzhi Gong, Zhenggang Xiong, Yangyang Zhao, Deguo Xing . Immune-related prognostic genes signatures in the tumor microenvironment of sarcoma. Mathematical Biosciences and Engineering, 2021, 18(3): 2243-2257. doi: 10.3934/mbe.2021113
    [6] Shengjue Xiao, Yufei Zhou, Ailin Liu, Qi Wu, Yue Hu, Jie Liu, Hong Zhu, Ting Yin, Defeng Pan . Uncovering potential novel biomarkers and immune infiltration characteristics in persistent atrial fibrillation using integrated bioinformatics analysis. Mathematical Biosciences and Engineering, 2021, 18(4): 4696-4712. doi: 10.3934/mbe.2021238
    [7] Jie Wang, Md. Nazim Uddin, Rehana Akter, Yun Wu . Contribution of endothelial cell-derived transcriptomes to the colon cancer based on bioinformatics analysis. Mathematical Biosciences and Engineering, 2021, 18(6): 7280-7300. doi: 10.3934/mbe.2021360
    [8] Xiangyue Zhang, Wen Hu, Zixian Lei, Hongjuan Wang, Xiaojing Kang . Identification of key genes and evaluation of immune cell infiltration in vitiligo. Mathematical Biosciences and Engineering, 2021, 18(2): 1051-1062. doi: 10.3934/mbe.2021057
    [9] Gaozhong Sun, Tongwei Zhao . Lung adenocarcinoma pathology stages related gene identification. Mathematical Biosciences and Engineering, 2020, 17(1): 737-746. doi: 10.3934/mbe.2020038
    [10] Sufang Wu, Hua He, Jingjing Huang, Shiyao Jiang, Xiyun Deng, Jun Huang, Yuanbing Chen, Yiqun Jiang . FMR1 is identified as an immune-related novel prognostic biomarker for renal clear cell carcinoma: A bioinformatics analysis of TAZ/YAP. Mathematical Biosciences and Engineering, 2022, 19(9): 9295-9320. doi: 10.3934/mbe.2022432
  • In this paper, we propose a delayed viral infection model with mitosis of uninfected target cells, two infection modes (virus-to-cell transmission and cell-to-cell transmission), and immune response. The model involves intracellular delays during the processes of viral infection, viral production, and CTLs recruitment. We verify that the threshold dynamics are determined by the basic reproduction number R0 for infection and the basic reproduction number RIM for immune response. The model dynamics become very rich when RIM>1. In this case, we use the CTLs recruitment delay τ3 as the bifurcation parameter to obtain stability switches on the positive equilibrium and global Hopf bifurcation diagrams for the model system. This allows us to show that τ3 can lead to multiple stability switches, the coexistence of multiple stable periodic solutions, and even chaos. A brief simulation of two-parameter bifurcation analysis indicates that both the CTLs recruitment delay τ3 and the mitosis rate r have a strong impact on the viral dynamics, but they do behave differently.



    Calcific aortic valve stenosis (CAVS) is a continuous global progressive disease that causes stenosis and contraction of the left ventricular outflow tract in the later stage of the disease, causing destructive damage to the heart that affects hemodynamics [1]. Several epidemiological studies have shown that 2.8% of the elderly (over 75 years old) have varying degrees of CAVS, and as many as 25% of the community population over 65 years old have risk factors for valve sclerosis. Older men, high triglyceride levels, smoking time limit, and early aortic valve replacement have been determined to be associated with the progression of CAVS [2]. A calcified aortic valve often leads to aortic valve stenosis. Inflammatory cell infiltration, lipid accumulation, and tissue fibrosis play a leading role in the initial stage of the mechanism of CAVS [3]. Therefore, exploring the pathophysiological process of CAVS is essential for the diagnosis and treatment of this complex disease with a poor prognosis.

    The immune and inflammatory response is a key link in the pathological process of CAVS [4]. A variety of inflammatory markers such as Toll-like receptor (TLR), interleukin-37, interleukin-6, transforming growth factor-β1 are closely related to aortic valve stenosis caused by calcified aortic valve Related [5,6]. Amyloid P Component Serum (APCS), Heat Shock Protein 90 (HSP90), Protein Disulfide Isomerase Family A Member 3 (PDIA3), Annexin A2 (ANXA2), Toll Like Receptor 7 (TLR7) and other immune-related genes (IRGs) also suggest that they have therapeutic effects in the process of CAVS fibrosis [7,8]. Immune cell infiltration is also closely related to CAVS. CD8 T lymphocytes, macrophages, and regulatory T lymphocytes (Tregs) appear in the pathophysiological process of CAVS [9,10,11]. In addition, statin anti-inflammatory therapy affects CAVS and maybe a drug target for the prevention of related diseases [12,13].

    This study conducted bioinformatics and machine learning analysis of IRGs in CAVS and explored potential regulatory methods and functional differences related to IRGs. Compared with traditional models, machine learning models show superior performance in disease classification and prediction [14]. The use of machine learning models is a novel method for disease diagnosis and prediction [15,16,17]. We use Single-sample GSEA (ssGSEA) and CIBERSORT to explore the relationship between differential immunity-related genes (DEIRGs) and the level of immune infiltration of cell subsets. In addition, to better understand the immune mechanism of CAVS, the potential connection between key immune signals and immune cell subsets was studied. The flow chart of this research analysis is shown in Supplementary Figure 1.

    Figure 1.  GSVA analysis results and differential genes. (A) GSVA analysis of gene sets. (B) DEGs and DEIRGs expression volcano graphs of CAVS and normal controls.

    Download the micro data set from Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/). The following are the screening criteria: 1) Select tissue samples including CAVS patients and normal controls; 2) Exclude samples from mitral and tricuspid valves; 3) No other organic diseases. The data set GSE12644 and GSE51472 of the GPL570 platform ([HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array) were selected as the training data set, and the data set GSE83453 of the GPL10558 platform (Illumina HumanHT-12 V4.0 expression beadchip) was selected as the verification. GSE12644 includes 10 calcified aortic valve samples and 10 normal controls (10 CAVS vs 10 NC). GSE51472 includes 5 calcified aortic valve tissues and 5 normal controls (5 CAVS vs 5 NC). GSE83453 includes 9 aortic valve tissues with stenosis and calcification and 8 normal controls (9 CAVS vs 8 NC). All tissue sample data undergoes background correction and de-batch effect.

    Use GSVA to analyze the data set to find pathways with significant differences between samples, analyzing genes more biologically meaningful [18]. The Benjamini & Hochberg method was used for multiple testing calibrations. The score value > 0.5 and adjust P-value < 0.05 is the cutoff value of pathway enrichment. The Limma package [19] (Version 3.44.3) in the R software (Version 4.0.2; https://www.r-project.org/) is used to screen differentially expressed genes (DEGs). The screening criteria for significant differences are P-value < 0.05, | log2 (Fold Change) | ≥ 1. A total of 2483 IRGs were obtained from the ImmPort (https://www.immport.org/) database, and DEIRGs were identified by matching IRGs and DEGs. The ggplot2 package [20] (Version 3.3.5) draws bar graphs of GSVA and volcano graphs of DEIRGs.

    Enrichment analysis of pathways and functions of DEIRGs to discover possible immune pathways and functions. Gene Ontology (GO) and Disease Ontology (DO) enrichment analysis are sorted by adjust P-value < 0.05 and Count value. We use the ggplot2 package (Version 3.3.5) of R software to visualize the plot. The Molecular Signatures Database (MSigDB) library is a collection of annotated gene sets. We can perform a series of analyses on the predefined gene set in the MSigDB library. Gene Set Enrichment Analysis (GSEA) uses "c2.cp.kegg.v7.2.symbols.gmt" under MSigDB (https://www.gsea-msigdb.org/gsea/msigdb/) [21] as the reference gene set, and the cutoff value is set to adjust P-value < 0.05.

    The obtained DEIRGs were further screened using machine learning methods to obtain key immune signals. Two machine learning methods—Least Absolute Shrinkage Selection Operator (LASSO) [22] and Support Vector Machine Recursive Feature Elimination (SVM-RFE) [23] perform feature screening for differences in gene expression values to obtain more accurate screening predictions. LASSO is an analysis method that can perform feature selection on research models. After screening, it aims to enhance the prediction accuracy and reliability of the model. SVM-RFE adopts the risk minimization principle and the experience error minimization principle. It can be used to improve learning performance to filter models. We use the glmnet package [24] (Version 4.1) and e1071 package [25] (Version 1.7) of the R software to execute the LASSO and SVM-RFE algorithms. Afterward, the external verification data set GSE83453 was verified against the selected key immune signals.

    ssGSEA [26] and CIBERSORT [27] are two tools for analyzing immune infiltrating cell subtypes. CIBERSORT analyzes the infiltration of immune cells between CAVS and normal controls. Subsequently, the relationship between the DEIRGs of CAVS and the subtypes of immune infiltrating cells was established. At the same time, the key immune signals obtained by machine learning screening are correlated with the immune infiltrating cells obtained by two analysis methods. Finally, in the MSigDB library "c5.all.v7.4.symbols.gmt" uses "IMMUNITY" as the keyword to find the immune pathway of interest, calculate the pathway enrichment score and analyze the key immune signals found and their correlation. Spearman correlation analysis was used for the correlation. We used the corrplot package [28] (Version 0.9) of R software to make a related heat map. The pheatmap [29] package (Version 1.0.12) constructs a heat map of immune cells. The vioplot package [30] (Version 0.3.7) is used to compare the levels of immune cells between the two groups.

    The Connectivity Map (CMap) (https://www.broadinstitute.org/connectivity-map-cmap) [31] is a database that analyzes the relationship between genes and their possible targeted drugs. Use CLUE (https://clue.io/) to predict the targeted drugs for CAVS immunotherapy on the 150 most significant up-and down-regulated immune genes in the data set. At the same time, use the Touchstone module to analyze the mechanism of actions (MoA) of the drug of interest and explore potential modes of action.

    All statistical analysis uses R software (Version 4.0.2). Use Student's t test for normally distributed variables and Mann-Whitney U test for abnormally distributed variables to compare the differences between the two groups.

    Use non-parametric unsupervised GSVA for gene sets to find the difference between the gene set in CAVS and normal control in the enrichment pathway, suggesting that it has a greater correlation with multiple immune pathways (Figure 1A). After the integration of GSE12644 and GSE51472 microarray matrix standardization and removal of batch effects, 266 DEGs (164 up-regulated genes, 102 down-regulated genes) were obtained, and 71 significant DEIRGs were obtained through integration with immune genes (60 up-regulated genes and 11 down-regulated genes). The volcano map shows the differences in genes. The DEIRGs related to immunity are marked by black circles. The names of some DEIRGs of interest are also marked (Figure 1B).

    GO enrichment found that the biological processes (BP) mainly focused on reactive oxygen species, oxidative stress, cellular response to oxidative stress, and decreased oxygen levels. Cell component (CC) is mainly related to various components of the synaptic membrane. Molecular function (MF) is mainly related to receptor activity and enzyme binding (Figure 2A). DO found that cardiovascular diseases such as arteriosclerotic cardiovascular disease, arteriosclerosis, atherosclerosis, coronary artery disease, and myocardial infarction were significantly enriched. In addition, it includes lung disease, obstructive lung disease, kidney disease, and so on (Figure 2B). GSEA found that two immune-related pathways, Cytokine-cytokine receptor interaction and Chemokine signaling pathway were significantly enriched. It may be related to the immune-related pathways of CAVS (Figure 2C).

    Figure 2.  Pathway and function enrichment analysis results. (A) Enrichment results of BP, CC, and MF. (B) Results of DO analysis. (C) Two immune-related pathways of CAVS.

    According to the LASSO method, the optimal lambda.min is set as 0.007915132 based on the amount of gene expression. The 6 immune signals screened out are Angiotensin II Receptor Type 1(AGTR1), C-X-C Motif Chemokine Ligand 16 (CXCL16), Leptin Receptor (LEPR), Phospholipid Transfer Protein (PLTP), Secretogranin II (SCG2), Secretory Leukocyte Peptidase Inhibitor (SLPI) (Figure 3A). After screening using SVM-RFE, the first 50 variables were screened by 5x cross-check, and the first 4 immune signals were Angiotensin II Receptor Type 1(AGTR1), Phospholipid Transfer Protein (PLTP), Secretogranin II (SCG2), and Tenascin C (TNC) (Figure 3B). Integrating the two results, AGTR1, PLTP, and SCG2 are considered by us to be the key immune signals of CAVS (Figure 3C). The external validation data set GSE83453 verified it and found that the three key immune signals distinguished well between CAVS and normal control (Figure 3D). The areas under the ROC curve were AGTR1 (AUC = 0.917), PLTP (AUC = 0.875), and SCG2 (AUC = 0.917), with high diagnostic value (Figure 3E).

    Figure 3.  Key immune signals and verification. (A) LASSO. (B) SVM-RFE. (C) 3 key immune signals integrated. (D) Verify the expression of key immune signals in the data set. (E) ROC curve.

    Two methods of CIBERSORT and ssGSEA were used to analyze the immune infiltration of CAVS and normal controls. The relative percentages of the 22 immune cells evaluated by CIBERSORT are displayed in a bar graph (Figure 4A). Correlation analysis between immune cells found that NK cells resting and T cells regulatory (Tregs) were positively correlated, and the correlation reached 0.81. The highest negative correlation between Mast cells resting and Mast cells activated, T cells gamma delta and T cells regulatory (Tregs) reached 0.77(Figure 4B). The violin chart shows that B cells naïve (P < 0.01), Macrophages M0 (P < 0.01), Macrophages M2 (P < 0.01) have a higher degree of discrimination between CAVS and normal controls. B cells naïve and Macrophages M2 are less in CAVS, while Macrophages M0 is more in CAVS (Figure 4C). The PCA chart also shows that CAVS is well distinguished from normal controls (Figure 4D). The correlation results between immune-related differential genes (DEIRGs) and immune infiltrating cells produced by the two methods are represented by heat maps (Figure 5A, B). ssGSEA found that most DEIRGs are positively correlated with more immune cells, while CIBERSORT shows that B cells memory, Macrophages.M2, and NK cells activated are negatively correlated with DEIRGs.

    Figure 4.  CIBERSORT immune infiltration analysis. (A) Percentage of immune cells. (B) Correlation between immune cells. (C) Immune cell difference between CAVS and normal control. (D) PCA chart of CAVS and normal control.
    Figure 5.  Correlation heat map of DEIRGs with CIBERSORT and ssGSEA immune cell subtypes. (A) Correlation of immune cells was obtained by DEIRGs and CIBERSORT analysis. (B) Correlation of immune cells obtained by DEIRGs and ssGSEA analysis.

    Among the three key immune signals of CAVS, AGTR1 is negatively correlated with the immune pathway of interest in MSigDB, while PLTP and SCG2 are positively correlated (Figure 6AC). There is a high degree of positive correlation between most immune pathways. Among the 28 immune cells analyzed by ssGSEA, AGTR1 was negatively correlated with most, while PLTP and SCG2 were positively correlated with most (Figure 6DF). The correlation among 28 kinds of immune cells is shown in the correlation heat map. Among the 22 immune cells analyzed by CIBERSORT, the most significant is that AGTR1 is positively correlated with NK cells activated (r = 0.83, P = 1.32E-08), and NK cells resting is negatively correlated (r = -0.63, P = 0.000194) (Figure 6G). PLTP had the most significant positive correlation with Neutrophils (r = 0.69, P = 3.89E-05), and the most significant negative correlation with B cells naive (r = -0.36, P = 0.045151) (Figure 6H). SCG2 had the most significant positive correlation with Macrophages.M0 (r = 0.68, P = 2.63E-05), and the most significant negative correlation with Macrophages.M2 (r = -0.56, P = 0.001529) (Figure 6I).

    Figure 6.  Key immune signals and specific immune pathways, immune infiltration subcellular correlation. (A) The correlation between AGTR1 and immune pathways. (B) The correlation between PLTP and immune pathways. (C) The correlation between SCG2 and immune pathways. (D) The correlation between AGTR1 and immune cell infiltration of ssGSEA. (E) The correlation between PLTP and immune cell infiltration of ssGSEA. (F) The correlation between SCG2 and immune cell infiltration of ssGSEA. (G) The correlation between AGTR1 and CIBERSORT immune cell infiltration. (H) The correlation between PLTP and CIBERSORT immune cell infiltration. (I) The correlation between SCG2 and CIBERSORT immune cell infiltration.

    CMap analysis found that isoliquiritigenin, parthenolide, pyrrolidine-dithiocarbamate, radicicol, RITA, roscovitine, securinine, midazolam, mitomycin-c, and colforsin are the top ten targeted drugs related to CAVS immunity (Figure 7A). The scores of isoliquiritigenin, parthenolide, and pyrrolidine-dithiocarbamate were 92.40, 88.23, and 84.24. The first three drugs isoliquiritigenin, parthenolide, and pyrrolidine-dithiocarbamate were analyzed by the drug MoA and found to be related to the NFkB pathway inhibitor, Immunostimulant, etc (Figure 7BD).

    Figure 7.  Drug CMap and MoA analysis. (A) Top ten drugs and scores. (B) MoA analysis of isoliquiritigenin. (C) MoA analysis of parthenolide. (D) MoA analysis of pyrrolidine-dithiocarbamate.

    Calcific aortic valve stenosis is the most prevalent valve disease in the world. It exists in large numbers in the elderly, and the disease has caused severe damage to the patients [32]. Once the clinical symptoms of severe CAVS appear, the prognosis is poor without intervention. Although there is still a lack of precise molecular insights into the pathophysiological process of CAVS, early intervention of the disease has become more realistic. More evidence shows that aortic valve calcification is closely related to immune inflammation [33,34]. Therefore, we tried to find DEIRGs and explore the possible role of immune cell infiltration in CAVS. 71 DEIRGs (60 up-regulated genes, 11 down-regulated genes) were identified as biomarkers of CAVS, and the potential enrichment function of DEIRGs was further studied. GO enrichment revealed that DEIRGs are associated with the immune inflammatory response. DO enrichment analysis found that DEIRGs are mainly closely related to cardiovascular diseases. GSEA enrichment revealed that two immune-related pathways, Cytokine-cytokine receptor interaction and Chemokine signaling pathway were significantly enriched. Two machine learning methods, LASSO and SVM-RFE, are used to confirm that AGTR1, PLTP, and SCG2 are the key immune signals of CAVS.

    The protein encoded by AGTR1 is part of the renin-angiotensin system and is used to regulate the balance of blood and body fluids. This gene may play a role in the production of arrhythmia after reperfusion after ischemia or infarcted myocardial blood flow is restored. A study showed that the AGTR1 gene has a moderate level of evidence that may be related to the risk of CAVS [35]. At the same time, the pathogenesis of some cardiovascular diseases is also related to AGTR1 [36,37]. PLTP is one of the lipid transfer proteins, which binds to Apolipoprotein A1(ApoA) [38]. Several studies have shown that ApoA is related to calcification and stenosis of the aortic valve. It is mostly located in its fibrous tissue and co-localizes with the calcified area [39,40,41]. In addition, studies in mice have shown that PLTP deficiency can reduce plasma total cholesterol and triglycerides and prevent the progression of arterial calcification [42]. SCG2 is a member of the granulin family. Studies have found that SCG2 is present in mouse myocardium [43]. In addition, a study also found that SCG2 plays a key role in the development of aortic valve calcification [44]. However, the specific mechanism of SCG2 and CAVS remains to be studied.

    CIBERSORT and ssGSEA analyzed the subtypes of CAVS immune infiltrating cells and found that B cells naïve (P < 0.01) and Macrophages M2 (P < 0.01) were less in CAVS, while Macrophages M0 (P < 0.01) were more in CAVS. ssGSEA and CIBERSORT found that most DEIRGs are positively correlated with immune cells, and DEIRGs are more negatively correlated with B cells memory, Macrophages.M2, and NK cells activated. Among the three key immune signals of CAVS, AGTR1 is negatively correlated with the immune pathway of interest in MSigDB, while PLTP and SCG2 are positively correlated. Analysis of the key immune signals and immune cell subtypes of the three CAVS showed that among the 28 immune cells obtained from ssGSEA, AGTR1 was negatively correlated with the majority, while PLTP and SCG2 were positively correlated with the majority. Among the 22 immune cells analyzed by CIBERSORT, the most significant is that AGTR1 is positively correlated with NK cells activated, and NK cells resting is negatively correlated. PLTP is positively correlated with Neutrophils and negatively correlated with B cells naive. SCG2 is positively correlated with Macrophages.M0 and negatively correlated with Macrophages.M2.

    Finally, we performed CMap analysis on immune genes in the dataset to find possible related drugs. Isoliquiritigenin, parthenolide, pyrrolidine-dithiocarbamate, radicicol, RITA, roscovitine, securinine, midazolam, mitomycin-c, and colforsin are the top ten targeted drugs related to the immune mechanism of CAVS. At the same time, it was found that the top three drugs isoliquiritigenin, parthenolide, and pyrrolidine-dithiocarbamate are related to NFkB pathway inhibitor and Immunostimulant. Based on these results, AGTR1, PLTP, and SCG2 seem to play a key role in CSVA by regulating immune infiltration.

    In this study, we obtained 266 DEGs in CAVS and normal controls, and 71 DEIRGs were obtained through integration with immune genes. Enrichment analysis found that DEIRGs are related to oxidative stress, synaptic membrane components, receptor activity, and a variety of cardiovascular diseases and immune pathways. Two machine learning algorithms identified AGTR1, PLTP, and SCG2 as the key immune signals of CAVS. Immune infiltration found that B cells naïve and Macrophages M2 are less in CAVS, while Macrophages M0 is more in CAVS. At the same time, AGTR1, PLTP, SCG2 are highly correlated with a variety of immune cell subtypes. CMap analysis found that isoliquiritigenin, parthenolide, and pyrrolidine-dithiocarbamate are the top three targeted drugs related to CAVS immunity. Our findings will help improve the understanding of CAVS disease and explain new molecular mechanisms and potential targets.

    We acknowledge GEO database for providing their platforms and contributors for uploading their meaningful datasets.

    The authors declare that they have no conflict of interest.



    [1] M. A. Nowak, S. Bonhoeffer, G. M. Shaw, R. M. May, Anti-viral drug treatment: dynamics of resistance in free virus and infected cell populations, J. Theor. Biol., 184 (1997), 203–217. https://doi.org/10.1006/jtbi.1996.0307 doi: 10.1006/jtbi.1996.0307
    [2] A. S. Perelson, P. W. Nelson, Mathematical analysis of HIV-1 dynamics in vivo, SIAM Rev., 41 (1999), 3–44. https://doi.org/10.1137/S0036144598335107 doi: 10.1137/S0036144598335107
    [3] R. V. Culshaw, S. Ruan, G. Webb, A mathematical model of cell-to-cell spread of HIV-1 that includes a time delay, J. Math. Biol., 46 (2003), 425–444. https://doi.org/10.1007/s00285-002-0191-5 doi: 10.1007/s00285-002-0191-5
    [4] Y. Wang, Y. Zhou, J. Wu, J. Heffernan, Oscillatory viral dynamics in a delayed HIV pathogenesis model, Math. Biosci., 219 (2009), 104–112. https://doi.org/10.1016/j.mbs.2009.03.003 doi: 10.1016/j.mbs.2009.03.003
    [5] H. Shu, L. Wang, J. Watmough, Global stability of a nonlinear viral infection model with infinitely distributed intracellular delays and CTL immune responses, SIAM J. Appl. Math., 73 (2013), 1280–1302. https://doi.org/10.1137/120896463 doi: 10.1137/120896463
    [6] Y. Wang, J. Liu, J. M. Heffernan, Viral dynamics of an HTLV-I infection model with intracellular delay and CTL immune response delay, J. Math. Anal. Appl., 459 (2018), 506–527. https://doi.org/10.1016/j.jmaa.2017.10.027 doi: 10.1016/j.jmaa.2017.10.027
    [7] J. Ren, R. Xu, L. Li, Global stability of an HIV infection model with saturated CTL immune response and intracellular delay, Math. Biosci. Eng., 18 (2021), 57–68. https://doi.org/10.3934/mbe.2021003 doi: 10.3934/mbe.2021003
    [8] A. Sigal, J. T. Kim, A. B. Balazs, E. Dekel, A. Mayo, R. Milo, et al., Cell-to-cell spread of HIV permits ongoing replication despite antiretroviral therapy, Nature, 477 (2011), 95–98. https://doi.org/10.1038/nature10347 doi: 10.1038/nature10347
    [9] S. Iwami, J. S. Takeuchi, S. Nakaoka, F. Mammano, F. Clavel, H. Inaba, et al., Cell-to-cell infection by HIV contributes over half of virus infection, eLife, 4 (2015), e08150. https://doi.org/10.7554/eLife.08150 doi: 10.7554/eLife.08150
    [10] N. L. K. Galloway, G. Doitsh, K. M. Monroe, Z. Yang, I. Muñoz-Arias, D. N. Levy, et al., Cell-to-cell transmission of HIV-1 is required to trigger pyroptotic death of lymphoid-tissue-derived CD4 T cells, Cell Rep., 12 (2015), 1555–1563. https://doi.org/10.1016/j.celrep.2015.08.011 doi: 10.1016/j.celrep.2015.08.011
    [11] W. Hübner, G. P. McNerney, P. Chen, B. M. Dale, R. E. Gordan, F. Y. S. Chuang, et al., Quantitative 3D video microscopy of HIV transfer across T cell virological synapses, Science, 323 (2009), 1743–1747. https://doi.org/10.1126/science.1167525 doi: 10.1126/science.1167525
    [12] P. De Leenheer, H. L. Smith, Virus dynamics: a global analysis, SIAM J. Appl. Math., 63 (2003), 1313–1327. https://doi.org/10.1137/S0036139902406905 doi: 10.1137/S0036139902406905
    [13] L. Wang, M. Y. Li, Mathematical analysis of the global dynamics of a model for HIV infection of CD4+ T cells, Math. Biosci., 200 (2006), 44–57. https://doi.org/10.1016/j.mbs.2005.12.026 doi: 10.1016/j.mbs.2005.12.026
    [14] M. Tsiang, J. F. Rooney, J. J. Toole, C. S. Gibbs, Biphasic clearance kinetics of hepatitis B virus from patients during adefovir dipivoxil therapy, Hepatology, 29 (1999), 1863–1869. https://doi.org/10.1002/hep.510290626 doi: 10.1002/hep.510290626
    [15] X. Lai, X. Zou, Modeling cell-to-cell spread of HIV-1 with logistic target cell growth, J. Math. Anal. Appl., 426 (2015), 563–584. https://doi.org/10.1016/j.jmaa.2014.10.086 doi: 10.1016/j.jmaa.2014.10.086
    [16] Y. Yang, L. Zou, S. Ruan, Global dynamics of a delayed within-host viral infection model with both virus-to-cell and cell-to-cell transmissions, Math. Biosci., 270 (2015), 183–191. https://doi.org/10.1016/j.mbs.2015.05.001 doi: 10.1016/j.mbs.2015.05.001
    [17] J. E. Mittler, B. Sulzer, A. U. Neumann, A. S. Perelson, Influence of delayed viral production on viral dynamics in HIV-1 infected patients, Math. Biosci., 152 (1998), 143–163. https://doi.org/10.1016/S0025-5564(98)10027-5 doi: 10.1016/S0025-5564(98)10027-5
    [18] P. W. Nelson, A. S. Perelson, Mathematical analysis of delay differential equation models of HIV-1 infection, Math. Biosci., 179 (2002), 73–94. https://doi.org/10.1016/S0025-5564(02)00099-8 doi: 10.1016/S0025-5564(02)00099-8
    [19] K. Wang, W. Wang, H. Pang, X. Liu, Complex dynamic behavior in a viral model with delayed immune response, Physica D, 226 (2007), 197–208. https://doi.org/10.1016/j.physd.2006.12.001 doi: 10.1016/j.physd.2006.12.001
    [20] S. S. Chen, C. Y. Cheng, Y. Takeuchi, Stability analysis in delayed within-host viral dynamics with both viral and cellular infections, J. Math. Anal. Appl., 442 (2016), 642–672. https://doi.org/10.1016/j.jmaa.2016.05.003 doi: 10.1016/j.jmaa.2016.05.003
    [21] H. Shu, Y. Chen, L. Wang, Impacts of the cell-free and cell-to-cell infection modes on viral dynamics, J. Dyn. Differ. Equations, 30 (2018), 1817–1836. https://doi.org/10.1007/s10884-017-9622-2 doi: 10.1007/s10884-017-9622-2
    [22] T. Wang, Z. Hu, F. Liao, Stability and Hopf bifurcation for a virus infection model with delayed humoral immunity response, J. Math. Anal. Appl., 411 (2014), 63–74. https://doi.org/10.1016/j.jmaa.2013.09.035 doi: 10.1016/j.jmaa.2013.09.035
    [23] Y. Yang, T. Zhang, Y. Xu, J. Zhou, A delayed virus infection model with cell-to-cell transmission and CTL immune response, Int. J. Bifurcation Chaos, 27 (2017), 1750150. https://doi.org/10.1142/S0218127417501504 doi: 10.1142/S0218127417501504
    [24] J. K. Hale, S. M. V. Lunel, Introduction to Functional Differential Equations, Springer-Verlag, New York, 1993.
    [25] H. Shu, M. Y. Li, Joint effects of mitosis and intracellular delay on viral dynamics: two-parameter bifurcation analysis, J. Math. Biol., 64 (2012), 1005–1020. https://doi.org/10.1007/s00285-011-0436-2 doi: 10.1007/s00285-011-0436-2
    [26] H. R. Thieme, Spectral bound and reproduction number for infinite-dimensional population structure and time heterogeneity, SIAM J. Appl. Math., 70 (2009), 188–211. https://doi.org/10.1137/080732870 doi: 10.1137/080732870
    [27] J. K. Hale, P. Waltman, Persistence in infinite-dimensional systems, SIAM J. Math. Anal., 20 (1989), 388–395. https://doi.org/10.1137/0520025 doi: 10.1137/0520025
    [28] H. Shu, X. Hu, L. Wang, J. Watmough, Delay induced stability switch, multitype bistability and chaos in an intraguild predation model, J. Math. Biol., 71 (2015), 1269–1298. https://doi.org/10.1007/s00285-015-0857-4 doi: 10.1007/s00285-015-0857-4
    [29] E. Beretta, Y. Kuang, Geometric stability switch criteria in delay differential systems with delay dependent parameters, SIAM J. Math. Anal., 33 (2002), 1144–1165. https://doi.org/10.1137/S0036141000376086 doi: 10.1137/S0036141000376086
    [30] J. Xu, Y. Zhou, Bifurcation analysis of HIV-1 infection model with cell-to-cell transmission and immune response delay, Math. Biosci. Eng., 13 (2016), 343–367. https://doi.org/10.3934/mbe.2015006 doi: 10.3934/mbe.2015006
    [31] J. Wu, Symmetric functional differential equations and neural networks with memory, Trans. Am. Math. Soc., 350 (1998), 4799–4838. https://doi.org/10.1090/S0002-9947-98-02083-2 doi: 10.1090/S0002-9947-98-02083-2
    [32] H. Shu, W. Xu, X. S. Wang, J. Wu, Complex dynamics in a delay differential equation with two delays in tick growth with diapause, J. Differ. Equations, 269 (2020), 10937–10963. https://doi.org/10.1016/j.jde.2020.07.029 doi: 10.1016/j.jde.2020.07.029
  • This article has been cited by:

    1. Yanli Du, Kun Wang, Xiannian Zi, Xiao Wang, Meiquan Li, Bo Zhang, Jinshan Ran, Wei Huang, Jing Wang, Cuilian Dong, Hanyi Xiang, Li Lei, Changrong Ge, Yong Liu, Combined Transcriptome and Metabolome Analysis of Stable Knockdown and Overexpression of the CD8A Gene in Chicken T Lymphocytes, 2024, 00325791, 104686, 10.1016/j.psj.2024.104686
    2. Ursula Houessou, Pardis Zamani, Hasanga D. Manikpurage, Zhonglin Li, Nathalie Gaudreault, Manel Dahmene, François Dagenais, Christian Couture, Marie-Annick Clavel, Philippe Pibarot, Patrick Mathieu, Yohan Bossé, Sébastien Thériault, Transcriptomic Signatures of Calcific Aortic Valve Stenosis Severity in Human Tricuspid and Bicuspid Aortic Valves, 2025, 2452302X, 10.1016/j.jacbts.2025.01.018
  • 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(2398) PDF downloads(164) Cited by(1)

Figures and Tables

Figures(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog