Processing math: 100%
Research article Special Issues

Analysis of a mathematical model with nonlinear susceptibles-guided interventions

  • In this paper, we considered a mathematical model describing the nonlinear susceptibles-guided vaccination and isolation strategies, incorporating the continuously saturated treatment. In this strategy, we find that the disease-free periodic solution can always exist, and consequently the control reproduction number can be defined through analyzing the stability of the disease-free periodic solution. Also, we discussed the existence and stability of the positive order-1 periodic solution from two points of view. Initially, we investigated the transcritical and pitchfork bifurcation of the Poincaré map with respect to key parameters, and proved the existence of a stable or an unstable positive order-1 periodic solution near the disease-free periodic solution. For another aspect, by studying the properties of the Poincaré map, we verified the existence of the positive order-1 periodic solution in a large range of the control parameters, especially, we verified the co-existence of finite or infinite countable different positive order-1 periodic solutions. Furthermore, numerical simulations show that the unstable order-1 periodic solution can co-exist with the stable order-1, or order-2, or order-3 periodic solution. The finding implies that the nonlinear susceptibles-triggered feedback control strategy can induce much rich dynamics, which suggests us to carefully choose key parameters to ensure the stability of the disease-free periodic solution, indicating that infectious diseases die out.

    Citation: Qian Li, Yanni Xiao. Analysis of a mathematical model with nonlinear susceptibles-guided interventions[J]. Mathematical Biosciences and Engineering, 2019, 16(5): 5551-5583. doi: 10.3934/mbe.2019276

    Related Papers:

    [1] Théo Lebeaupin, Hafida Sellou, Gyula Timinszky, Sébastien Huet . Chromatin dynamics at DNA breaks: what, how and why?. AIMS Biophysics, 2015, 2(4): 458-475. doi: 10.3934/biophy.2015.4.458
    [2] Bertrand R. Caré, Pierre-Emmanuel Emeriau, Ruggero Cortini, Jean-Marc Victor . Chromatin epigenomic domain folding: size matters. AIMS Biophysics, 2015, 2(4): 517-530. doi: 10.3934/biophy.2015.4.517
    [3] Larisa I. Fedoreyeva, Boris F. Vanyushin, Ekaterina N. Baranova . Peptide AEDL alters chromatin conformation via histone binding. AIMS Biophysics, 2020, 7(1): 1-16. doi: 10.3934/biophy.2020001
    [4] E.N. Baranova, R.M. Sarimov, A.A. Gulevich . Stress induced «railway for pre-ribosome export» structure as a new model for studying eukaryote ribosome biogenesis. AIMS Biophysics, 2019, 6(2): 47-67. doi: 10.3934/biophy.2019.2.47
    [5] Vladimir B. Teif, Andrey G. Cherstvy . Chromatin and epigenetics: current biophysical views. AIMS Biophysics, 2016, 3(1): 88-98. doi: 10.3934/biophy.2016.1.88
    [6] Claudia Tanja Mierke . Physical role of nuclear and cytoskeletal confinements in cell migration mode selection and switching. AIMS Biophysics, 2017, 4(4): 615-658. doi: 10.3934/biophy.2017.4.615
    [7] Michael J. Hamilton, Matthew D. Young, Silvia Sauer, Ernest Martinez . The interplay of long non-coding RNAs and MYC in cancer. AIMS Biophysics, 2015, 2(4): 794-809. doi: 10.3934/biophy.2015.4.794
    [8] Raghvendra P Singh, Guillaume Brysbaert, Marc F Lensink, Fabrizio Cleri, Ralf Blossey . Kinetic proofreading of chromatin remodeling: from gene activation to gene repression and back. AIMS Biophysics, 2015, 2(4): 398-411. doi: 10.3934/biophy.2015.4.398
    [9] Larisa I. Fedoreyeva, Tatiana A. Dilovarova, Boris F. Vanyushin, Inna A. Chaban, Neonila V. Kononenko . Regulation of gene expression in Nicotiana tabacum seedlings by the MKASAA peptide through DNA methylation via the RdDM pathway. AIMS Biophysics, 2022, 9(2): 113-129. doi: 10.3934/biophy.2022011
    [10] Audrey Lejart, Gilles Salbert, Sébastien Huet . Cytosine hydroxymethylation by TET enzymes: From the control of gene expression to the regulation of DNA repair mechanisms, and back. AIMS Biophysics, 2018, 5(3): 182-193. doi: 10.3934/biophy.2018.3.182
  • In this paper, we considered a mathematical model describing the nonlinear susceptibles-guided vaccination and isolation strategies, incorporating the continuously saturated treatment. In this strategy, we find that the disease-free periodic solution can always exist, and consequently the control reproduction number can be defined through analyzing the stability of the disease-free periodic solution. Also, we discussed the existence and stability of the positive order-1 periodic solution from two points of view. Initially, we investigated the transcritical and pitchfork bifurcation of the Poincaré map with respect to key parameters, and proved the existence of a stable or an unstable positive order-1 periodic solution near the disease-free periodic solution. For another aspect, by studying the properties of the Poincaré map, we verified the existence of the positive order-1 periodic solution in a large range of the control parameters, especially, we verified the co-existence of finite or infinite countable different positive order-1 periodic solutions. Furthermore, numerical simulations show that the unstable order-1 periodic solution can co-exist with the stable order-1, or order-2, or order-3 periodic solution. The finding implies that the nonlinear susceptibles-triggered feedback control strategy can induce much rich dynamics, which suggests us to carefully choose key parameters to ensure the stability of the disease-free periodic solution, indicating that infectious diseases die out.


    Abbreviations

    HP1Heterochromatin Protein 1
    PcG proteinsPolycomb group of proteins
    PRC1Polycomb Repressive Complex 1
    PRC2Polycomb Repressive Complex 2
    Pax7Paired box protein Pax-7
    3CChromosome Conformantion Capture
    4CChromosome Conformantion Capture on Chip
    Hi-CHigh resolution Capture
    PREPolycomb Response elements
    TADTopological Associated Domain
    HOXDHomeobox protein Hox D
    MET-2Histone-lysine N-methyltransferase met-2
    SET-25SET-25
    NBNuclear Body
    FRAPFluorescence Recovery After Photobleaching
    HOXHomeobox protein Hox
    ESCEmbryonic Stem Cell
    LADLamina Associated Domain
    ASCAdipose Stem Cell
    Klf4Kruppel-like factor 4
    Oct4Octamer-binding transcription factor 4
    NPCNeuronal Precursor Cell
    HGPSHutchinson-Gilford Progeria Syndrome
    NPCNuclear Pore Complex
    AFAtrial Fibrillation
    NupNucleoporin
    NADNucleolar Associated Domain
    LIN28AProtein lin-28 homolog A
    BERBase Excision Repair
    NERNucleotide Excision Repair
    MMRMismatch Repair
    DSBDouble Strand Break
    NHEJNon Homologous End Joining
    HRHomologous Recombination
    Ku70 (XRCC6)X-ray repair cross-complementing protein 6
    Ku80 (XRCC5)X-ray repair cross-complementing protein 5
    DNA-PKcsDNA-dependent Protein Kinase
    ATMAtaxia Telangiectasia Mutated
    ATRAtaxia Telangiectasia and Rad3-related protein
    XRCC4X-ray repair cross-complementing protein 4
    RPAReplication Protein A
    Rad 51DNA repair protein RAD51 homolog 1
    Rad 52DNA repair protein RAD52 homolog
    DDRDNA Damage Response
    53BP1Tumor suppressor p53-binding protein 1
    MRE11Meiotic Recombination 11
    NSB1Nijmegen breakage syndrome protein 1
    ROSReactive Oxygen Species
    NADPHNicotinamide Adenine Dinucleotide Phosphate
    RbRetinoblastoma protein
    p53Cellular tumor antigen p53
    NF-kBNuclear factor NF-kappa-B
    FoxOForkhead box protein O
    HIV-1Human Immunodeficiency Virus 1

    1. The First Level of the Epigenome Complexity: Histone Modifications

    1.1.The histone structure

    In eukaryotic cells, DNA is packaged inside the nucleus through its ordered and repetitive aggregation with histones and is bound by several proteins and RNA molecules: the overall complex of nucleic acids and their associated proteins is known as “chromatin”. The first level of chromatin organization is the nucleosome: four histone proteins; H2A, H2B, H3, H4 surrounded by DNA [1]. Each nucleosome is separated from the next one by a DNA linker, variable in length, which is associated to histone H1 and or its variants [2,3]. These DNA-histones complexes fold in 30 nm fibers. This process is driven in part by the histone H1, which regulates the intranuclear electrostatic balance and by nucleosome repeats which ensure the stabilization of chromatin fibers [4,5,6]. In interphase the chromatin shows different levels of compaction, clearly visible under electron microscopy, that allowed scientists to classify it into euchromatin and heterochromatin, observed for the first time by Emil Heizt (1928) [7].The more decondensed euchromatin is generally enriched with expressed genes (reviewed in [8,9]). In fact, being more accessible, euchromatin favours the binding of transcription factors and RNA polymerase [10] and is correlated to high transcriptional activity [11]. On the other hand, heterochromatin is more condensed, gene-poor and is further divided into “facultative” and “constitutive”. The constitutive heterochromatin includes pericentric and telomeric regions [12,13,14], both critically important for the higher order structure folding and the maintenance of intact chromosomes [8]. The facultative heterochromatin is a repressive environment able to rapidly shift between activated and repressed state through the crosstalk with epigenetic regulators that lead to a chromatin reorganization [8].

    1.2. Histone signatures

    Histones are characterized by amino and carboxy-protruding tails that are the target of several modifications, critically important for regulating the balance between activated or repressed genes. These modifications mainly include acetylation, methylation, phosphorylation, ubiquitination and sumoylation (reviewed in [15]). Multiple histone modifications can also coexist on the same tail, dictating specific biological readouts [16]. Euchromatin and heterochromatin are distinguishable biochemically by different covalent histone modifications [17] (Figure 1). Euchromatin is enriched with H3 and H4 acetylation and H3K4 methylation [18,19] while facultative heterochromatin is enriched with hypoacetylated histones, H4K20me1, H3K9me2, H3K27me3, H2AK119ub1 (reviewed in [20]). On the other hand, constitutive heterochromatin is marked with H3K9me3 catalyzed by Suv39H [21], which in turn stimulates the recruitment of Heterochromatin Protein 1(HP1) [22], so contributing to form a closed chromatin structure normally associated with a repressive transcriptional state. However, despite this classification, histone marks should not be considered exclusively restricted to specific chromatin compartments because they were also found in other chromatin states. For instance, although H3K9me3 and HP1 are considered key features of heterochromatin, they have also been found in euchromatic region in flies [23,24] and in mammals [25], thus also suggesting their role in reinforcing gene repression in euchromatin.

    Figure 1. Travelling into the depths of the nucleus: chromatin architecture. The scheme represents the epigenome structure. The nuclear membrane, defines the border of nucleus. The transcriptionally activated chromatin, euchromatin (green), appears decondensed, tightly associated to nuclear pore and enriched by histone modifications H3ac, H4ac and H3K4me. The heterochromatin (red) is a repressive environment for gene expression and appears at nuclear periphery, often associated to nuclear lamina, condensed and enriched of H3K27me3, H3K9me3, H2A119ub. Inside the nucleus, euchromatin and heterochromatin give rise to several grades of higher order structures: chromosome loops, Topological Associated Domains (TADs), Lamin Associated Domains (LADs) and chromosomal territories. Also the nucleolus, the “assembly-chain” of ribosomes, associates with specific DNA regions: the Nucleolar Associated Domains (NADs), that surround the highly transcribed region of nucleolus, giving rise to another grade of chromatin organization.

    2. The Second Level of the Epigenome Complexity: Chromatin Modifiers

    At the chromatin level, specific epigenetic factors are able to post-translationally modify the histone proteins. Then, specific set of histone modifications generate a chromatin environment that is recognized by specific binding proteins; these, in turn, influence gene expression and other chromatin functions and represent an additional layer of gene regulation complexity. One of the most studied families of epigenetic regulators that include both histone modifiers and chromatin binding proteins are the Polycomb group of proteins (PcG proteins) which act as a repressor of gene expression (reviewed in [26]). PcG proteins aggregate into different complexes. The most studied are Polycomb Repressive Complex 1 and 2 (PRC1 and PRC2) that deposit H2AK119ub and H3K27me3 marks, respectively. The catalytic subunit responsible for the histone signature placement are E3 ubiquitin-protein ligase RING1 (Ring 1A/B) for PRC1 and Enhancer of zeste 2 (Ezh2) for PRC2 (reviewed in [27]). A classical, hierarchical model suggests that the PRC2-dependent H3K27me3 signature recruits the PRC1 complex which, through the ubiquitination of H2AK119, inhibits RNA polymerase II-dependent transcriptional elongation [28,29,30]. In recent years, several studies have challenged this model (reviewed in [31]) by showingindependent functions and recruitment mechanisms for PRC1 and PRC2. Thus, it has been demonstrated that PRC1 and PRC2 complexes share only a subset of binding sites [28,32,33] and that there are several PRC1 complexes containing distinct subunits that differ in genomic binding [34]. Moreover, PRC1 can recruit PRC2 [32,35,36] while recruitment of PRC1 and PRC2 can also be mediated by noncoding RNA (reviewed in [37]), specific factors or structures [38,39,40,41] and by their own histone marks [42,43]. These findings suggest that PRC recruitment relies on several alternative mechanisms and that the combination of mechanisms tethering PcG to chromatin is probably locus-specific. The plasticity of the facultative heterochromatin to switch on or off the gene expression through histone modifications is extremely important for lineage commitment, development and cell differentiation, when a proper timing of gene expression is needed [28,44,45]. Most of differentiation genes in mammals are under the control of Polycomb group of proteins [46]. Upon differentiation stimuli, PcG proteins leave lineage specific muscle promoters and bind genes important for stemness maintenance [46]. This dynamics has been extensively shown in muscle differentiation: in muscle stem cells (satellite cells) muscle specific genes are maintained repressed by PcG proteins; at the onset of muscle differentiation PcG proteins are displaced from muscle gene and are relocalized at stemness genes, such as Pax7 (marker of satellite cells), ensuring the correct timing of the muscle differentiation [47,48,49,50,51].

    3. The Third Level of the Epigenome Complexity: Chromatin Higher Order Structures

    3.1. DNA looping

    In addition to histone modifications, the chromatin can fold in specific higher order structures that favour the activation or the repression of genes (reviewed in [52]). In particular, structure-mediated gene expression control regards genes that are regulated by a large set of enhancers, insulators or repressors, which are often located at a considerable distance from the target gene. Thus, the establishment of topologically distinct chromatin domains could play a fundamental role in the modulation of gene expression (Figure 1). PcG proteins also exert gene repression by mediating chromatin higher order structures (reviewed in [26,53]). The development of Chromosome Conformation Capture (3C) technology [54] and its derivative technologies (reviewed in [55]), such as Chromosome Conformation Capture on chip (4C) [56] and High resolution Capture (Hi-C) [57] shed light on the chromatin contacts occurring in the nucleus allowing the high-through-put mapping of the genome conformation. The use of these technologies has allowed important advances in understanding PcG functions, demonstrating that the coordinated action of PcG proteins is required to mediate the formation of multi-looped structures where all the major PcG targets are gathered together, by cis and trans interactions [58,59,60,61]. This PcG-mediated DNA conformation is cell cycle regulated and fundamental for the maintenance of gene repression [58,59,62,63]. Notably, although PcG recruitment is not conserved from fly to mammals, being mediated in Drosophila by specific DNA sequences called Polycomb Response Elements (PRE) [64], the higher order structures of PcG-regulated genes are similarly organized in multi-loops among different species, thus suggesting a functional role for DNA higher order structures in the control of PcG mediated transcriptional repression.

    3.2. Topological Associated Domains (TADs)

    At the genome wide level, the combination of all DNA-DNA interactions determines the formation of specific macro-domains of active and inactive chromatin, named Topological Associated Domains (TADs) observed in human, in mouse and in Drosophila melanogaster genomes [65,66,67] (Figure 1). TADs are characterized by a central region of high chromatin interactions where reside the tissue-specific genes, surrounded by less-interacting boundary regions enriched of housekeeping genes [52]. Although TADs share a common organization in different cell-types and are partially conserved from human to mouse, suggesting a very stable conformation, they can be subjected to important structural rearrangements. During the mouse hind limb development, the HOXD genes are expressed following specific timing from early genes to later genes; this transition is guided by the reorganization of TADs which allow the switch of HOXD locus interactions from one compartment to another one [68]. This dynamics is required for the formation of an intermediate area that will develop into the wrist.

    3.3. Chromosome territories

    Another level of chromatin organization is represented by the chromosome position. Several studies have shown that chromosomes are not randomly assembled and packaged in the nucleus, but occupy specific regions called chromosome territories (reviewed in [69]) (Figure 1). The positioning of these “territories” is cell specific and conserved between human and other primates, suggesting a functional role of specific chromosome organization inside the nucleus [52,70,71]. In fact, chromosome territories have an internal structure thought to modulate the accessibility of proteins and consequently regulates the gene expression through its remodeling (reviewed in [72,73]). In particular, it has been shown that within the territories the gene-rich regions are separated from heterochromatin and generally are positioned at the periphery of chromosome territories while gene-poor regions are at the interior of chromosome territories [74,75]. How chromosome territories form and occupy specific regions inside the nucleus is still unclear, but researchers have suggested two different hypothesis. The first suggests that preferential positioning of chromosome territories inside the nucleus is self-organized and guided in part by the interchromatin compartment, which creates channels for the interconnection between the chromosome territories and in part by the content of gene-rich or gene-poor areas of each specific chromosome [72,76,77]. According to this hypothesis chromosome folding within territories should allow that the gene-rich zones (euchromatin) will assume a position compatible with expression, whereas not commonly expressed genes, tissue specific or developmental regulated genes (heterochromatin) will be buried inside the chromosome territories and expressed only after a conformational change [72]. However, some exceptions to this role have been described, showing ubiquitously and tissue-specific genes located within the chromosome territories after their activation, suggesting that the inner region of chromosome territories is not an insurmountable barrier to transcription machinery and transcription factors [74,78]. The second emerging hypothesis supports the involvement of nuclear proteins and nuclear envelope in the chromosome positioning (reviewed in ([76,79]). This hypothesis does not conflict with the first one and could also explain how the chromosome can be oriented in the nuclear space. Thus, even though some progress has been made on the functional role of chromosomal territories in controlling gene expression, there is still a long way to go before discovering the mechanism behind their positioning and assembling inside the nucleus.

    4. The Fourth Level of the Epigenome Complexity: The Nuclear Environment

    4.1. Nuclear periphery

    Inside the nucleus the coexistence of different environments ensure an additional level of gene transcription regulation (reviewed in [52]). Indeed, chromatin is non-homogenously distributed inside the nucleus, as shown by electron microscopy. In particular, most of the darkly stained chromatin (heterochromatin) is found around the nucleoli or at the nuclear periphery interrupted with euchromatin at nuclear pores (Figure 1). One of the most studied microenvironments is the nuclear periphery directly linked to a repressive transcriptional state (reviewed in [80,81]). Indeed, artificial tethering of a reporter gene to the nuclear periphery favours the establishment of silent chromatin, by mediating the access to local high concentration of repressors and chromatin modifiers, as observed in yeast and human cells (reviewed in [80] and [82,83]). Although endogenous molecular mechanisms driving the periphery localization are still under investigation, a genetic screen in worms revealed that histone H3 lysine 9 (H3K9) methylation marks target repetitive heterochromatin to the nuclear envelope and that the combined elimination of the genes encoding two histone H3K9 methyl transferase, MET-2 and SET-25, determines the loss of peripheral localization [84]. Interestingly, loss of SET-25 alone de-represses the transgene array, but does not lead to array delocalization, suggesting that H3K9me1 or me2, catalyzed by SET-25, are sufficient to anchor chromatin to the nuclear envelope and that further modification, generating H3K9me3, is necessary for transcriptional repression [84]. In line with and further corroborating these observations, in budding yeast, several highly regulated classes of genes were described as translocating to the nuclear periphery concurrent with their activation (reviewed in [85]). Interestingly, it has been shown that when an inducible active gene is repressed again, its persistent presence at the nuclear periphery allows a rapid reactivation upon stimuli. This process has been described as “transcriptional memory, ” since it allows the cell to “remember” recently activated loci (reviewed in [85]), so suggesting that the nuclear periphery environment can establish an epigenetic memory. Thus, in the nuclear periphery there is an alternation on permissive and non-permissive environment that can influence a given gene’s transcriptional state.

    4.2. The nuclear bodies (NBs)

    Within the nuclear space the non-uniform spatial distribution of macromolecules gives rise to a formation of a variety of functionally distinct nuclear bodies (NBs), highly dynamic structures formed by protein-protein interactions and having a large impact on gene regulation (reviewed in [86,87]). NBs are nuclear niches, in which machineries responsible for specific nuclear functions, such as transcription (transcriptional factories), replication (replication foci) or repression (repressive bodies) are situated. The presence of NBs implies a non-diffusive transport of factors involved in chromatin regulation. This compartmentalization of nuclear processes avoids the stochastic non-directional transport of nuclear regulators and speeds up the dynamics in the nucleus [87]. The most prominent NB in the nucleus is the nucleolus, responsible of ribosomal RNA transcription (Figure 1). Interestingly the nucleolus has a specific structural organization in three distinct regions that mirror nucleolar functions (reviewed in [86]). An example of epigenetic regulators forming aggregates in the nucleus are again the PcG proteins, forming structures called PcG bodies, often localized close to pericentric heterochromatin [88] (Figure 2). The formation of PcG bodies is mediated by intrinsic [89]

    Figure 2. The Polycomb bodies. The immunofluorescence depicts in green the Polycomb bodies, stained with the PRC1 component: Bmi-1 for murine C2C12 myoblasts and human HeLa cancer cells and PC for Drosophila embryonic Schneider cells (S2). PcG proteins aggregate in these distinct domains and create a repressive environment for gene expression, clustering the repressed target genes.

    and extrinsic [90] protein-protein interactions. Studies on the kinetics of PcG proteins exchange with the body have revealed that these protein aggregates are very stable and that distinct PcG factors exhibit large variability in their mobilities [91]. In terms of kinetics, three distinguishable components form the foci: fast, slow and an immobile on the time scales of measurements [91]. The slow binding component amplitudes directly correlate with the intensities and sizes of PcG foci, this suggesting that the PcG bodies contain multiple PcG targets. On the other hand, FRAP analysis has revealed that the size and the kinetics inside the foci is mostly slower than PcG complexes localized outside the foci [89,91], suggesting a coordinated overlap between genome higher order structures and PcG proteic compartments. In fact, previous studies revealed that Drosophila Hox clusters co-localize in PcG bodies in the tissue where they are repressed, this nuclear localization being required for efficient silencing [58]. In line with these findings, PcG bodies are clearly distinguishable at the end of early embryogenesis; they progressively increase in size and number during stages 5-11 of late embryogenesis when the transcriptional programs need to be maintained through cell division [92]. Similarly, the mobility of PcG bodies increases at the onset of ESC (Embryonic Stem Cell) differentiation and decreases as differentiation progresses [93]. Further studies of the composition of different PcG foci compartments could shed light on the coordinated functions of PcG complexes.

    In spite of all classifications of epigenetic signatures, the epigenome is highly dynamic and can rapidly respond to the environmental stimuli. The occasional presence of repressive marks on the heterochromatin, as in the case of H3K9me3 and HP1 [11], the heterochromatin’s ability to switch from a repressive to an active state (or vice-versa) and recent findings showing that epigenetic repressors can be bound to active transcribed chromatin [94,95] all underline the thin line separating euchromatin from heterochromatin. All these evidence also demonstrate that chromatin cannot be considered a static entity rigorously classified in an active or repressive state, but should be considered dynamic, and that this plasticity is the key for precise gene expression.

    5. Nuclear Architecture and Gene Expression Regulation: The Nuclear Lamina

    5.1. Lamin proteins

    The nuclear lamina is a proteinaceous layer tightly associated with the inner nuclear membrane. The main components of the nuclear lamina are the lamin proteins [96], classified as intermediate filament V and further divided in vertebrates into A-type lamins (lamin A and lamin C) and the B-type lamins (lamin B1 and lamin B2) [97,98]. Whereas B-type lamins are ubiquitously expressed and firmly anchored at the nuclear membrane, A-type lamins are mostly absent from early embryos [99,100], developmentally regulated, expressed in lineage-committed progenitor cells and differentiated cells [101,102]. Moreover, A-type lamins are present in both peri-nuclear and nucleoplasmatic compartments [103]. A and B-type lamins, form distinct homopolimers that interweave each other at the inner membrane to create an interconnected plot [104]. One of the lamina’s roles is to maintain nuclear shape and integrity. Indeed, knocking down lamin B1 results in an enlargement of Lamin A and Lamin B2 network and in the formation of nuclear blebs. These alterations are accompanied by a perturbation of epigenetic signatures, as a reduction of H3K9me2 and an increase of H3 acetylation, suggesting that the impairment of the nuclear envelope affects the proper organization of chromatin [104]. Other cellular functions are attributed to the lamins such as transcription, cell cycle regulation, development and differentiation (reviewed in [105]).

    5.2. Lamina Associated Domains (LADs)

    The importance of lamins and nuclear lamina as epigenetic regulators is reinforced by the fact that the nuclear lamina interacts with the genome at specific DNA sequences called Lamina Associated Domains (LADs) (Figure 1). These domains show variable lengths from 0.1 to 10 Mb and generally create an environment that maintains genes repressed [106,107] and marked by H3K9me2, H3K9me3 or H3K27me3 [108,109,110,111]. Interestingly, the borders of LADs include promoters that often contain transcription-associated histone modification indicating the transcriptional start site of the activated genes [112]. This characteristic may suggest that LADs define genomic regions in which genes are repressed, neighboured by activated genes. Some LADs are cell-specific, suggesting a possible role in cell identity and differentiation [106,113]. Recently, it has been shown in worms that in a lamin-defective background muscle promoters do not re-localize from the nuclear periphery to a more internal location in differentiating muscle cells [114]; this further indicates that the genome architecture dynamics mediated by Lamin A are necessary for a correct differentiation. This was confirmed by studies in mammals, showing that during the differentiation of adipose stem cells (ASCs) into adipocytes, genes that regulate a specific cell lineage, are released from the lamina while genes that maintain the cell undifferentiated are retained near the nuclear lamina [110]. Similarly, during the differentiation of Embryonic stem cells (ESCs) into neural precursor cells (NPCs), the so named “stemness genes” such as Nanog, Klf4 and Oct4, exhibit significantly increased interactions with the nuclear lamina in NPCs compared to ESCs [115]. In parallel, other LADs, seem to be constitutively associated to the lamina and conserved between species, suggesting a role in the organization of chromosome architecture [106,107]. One possibility is that the presence of these constitutive LADs could help the correct DNA organization (and gene expression) in daughter cells after mitosis, acting as a “bookmarks” of the genome architecture [116]. However, in contrast with this hypothesis, recent evidence has shown that LADs-nuclear lamina interactions are not rigid and fixed domains, but on the contrary, are dynamic and organized in a random manner after cell division; this stochastic organization is related to H3K9me2 and guided in part by Lamin A [108,117]. Moreover, in spite of all the efforts until now to characterize the LADs and the intriguing hypothesis concerning the role of LADs during genome organization, a recent work has questioned some aspects of LADs, asserting that nuclear lamins are not essential for the formation of LADs [118]. Thus it is clear that the research is still far from clearly defining the role of these domains in genome organization and regulation of gene expression.

    5.3. Lamins and mechanotrasduction

    A-type Lamin proteins, being situated in a strategic position at the interface between the nuclear membrane and chromatin, cover an important role also in mechanotrasduction, the process by which the external mechanical stimuli are converted into biochemical signals inside cells (reviewed in [119]). An intriguing hypothesis suggests that mechanical signals from the extracellular environment can be transmitted physically by the contractile cytoskeleton to the nucleus through the nuclear membrane and the nuclear lamina [120] and then translated into gene expression response through Lamin-mediated distortion of chromatin higher order structure and changes in transcription factor accessibility [121]. Recent findings support this hypothesis, showing that stress, matrix stiffness and the environment act on gene expression by perturbing the Lamin A expression and consequently the lamin A:B ratio. The altered Lamin A content, in turn, influences specific cell lineage determination and differentiation [122].

    5.4. The laminopathies

    Mutations in lamin proteins can cause a wide set of diseases, named laminopathies, affecting specific tissues or acting upon multiple types of tissue creating overlapping or systemic phenotypes. Only few mutations in Lamin B have been described as being compatible with life [123,124] so confirming its fundamental role in sustaining the nuclear structure [125,126]. On the other hand, pathology with mutations related to Lamin A are more frequent and give rise to muscular dystrophy, lipodystrophy, neuropathy, cardiomyopathy and HGPS syndrome (Hutchinson-Gilford progeria syndrome) (reviewed in [126]). Many of the laminopathies are characterized by nuclear abnormalities [127,128]. Thus it is reasonable to assume that in these diseases, mutations in lamins or alteration of lamin A:B ratio act on gene expression in different ways: on the one hand by increasing nuclear stress sensitive and improper signalling with an aberrant mechanotrasduction [129] and on the other hand by altering the chromatin organization and modulating gene expression epigenetically. Both these aspects could be involved in the pathogenesis and disease progression [122,129]. Hutchinson-Gilford Progeria Syndrome (HGPS) is one of the most known and studied diseases belonging to laminopathies and is characterized by a premature aging. In 90% of cases, the disease is caused by a point mutation 1824C > T that activates a cryptic site of splicing creating an immature form of Lamin A/C (called progerin) that accumulates at the nuclear membrane, giving rise to abnormal nuclear shapes (reviewed in [130]). Several epigenetic marks have been found altered in HGPS, such as a reduction of H3K27me3, H3K9me3 and lamina-heterochromatin associations [131,132]. Although it is not clear if these alterations precede the senescent phenotype, it has been proposed that epigenetic dysfunctions and aberrant chromatin higher order structures could play a role in HGPS pathogenesis and progression [131]. This suggests a direct interplay between the alteration of the nuclear membrane and epigenome misregulation. Further confirming the idea that Lamin A/C plays a key role in epigenome regulation, studies on laminopathies have revealed that distinct mutations in Lamin A gene can give rise to the same clinical condition while a single mutation can result in different phenotypes and disease severity [133,134,135], suggesting an involvement of the individual epigenetic background to lamin dependent diseases.In line with this hypothesis, recently we found that Lamin A/C maintains transcriptional repression at muscle-specific genes by sustaining the intranuclear positioning of PcG bodies and the activity of PcG proteins [136]. Lamin A/C depletion does indeed lead to an earlier muscle differentiation onset due to a premature PcG displacement from their targets.

    6. Nuclear Architecture and Gene Expression Regulation: The Nuclear Pore

    The nuclear pore complexes (NPCs) are multi-component protein complexes that form selectively permeable channels through the nuclear envelope, mediating nucleo-cytoplasmic transport of molecules between the nucleus and the cytoplasm (reviewed in [137]). The NPC is composed of multiple copies of nucleoporins. The overall structure of the NPC is highly conserved and displays nucleoporins organized into an eightfold symmetrical ring structures lining the circumference of the pore where they interact with the membrane. Various studies have implicated the NPCs as the binding platform for several highly transcribed housekeeping genes and genes strongly induced by changes in environmental conditions [137] (Figure 1). The presence of several expressed loci in the close proximity of the nuclear pore inspired the gene gating hypothesis, in which active genes are found near nuclear pore complexes to facilitate efficient export of mRNAs from the nucleus [138,139,140]. 3C analysis confirmed this hypothesis showing that some inducible genes are arranged in chromatin loops the ends of which are anchored at the NPC [141]. However, NPCs have also been associated with transcriptional repression [142,143,144,145], suggesting distinctive roles for individual components of the NPC in facilitating both gene activation and repression through its crosstalk with the epigenome. The NPC control of gene transcription is further confirmed by their central role during cell differentiation [38,146,147,148,149,150,151]. Intriguingly, NPC have been also found to be involved in the maintenance of the pluripotency via controlling levels of the pluripotency factors in the nucleus [152,153], underlying functional roles of the NPC in regulating developmental states and transitions. As described above for lamins, also mutations in key components of the NPC have been described to be associated with specific diseases. Thus, inherited cases of a cardiac disorder, the atrial fibrillation (AF), have been mapped to a missense mutation in the human nucleoporin Nup155 [154]. Studies in mice and human have shown that the impairment of Nup155 nuclear localization affects mRNA and protein transport leading to cardiovascular disease. Other studies have shown that mutations in the nucleoporin Nup62 are responsible for the familial form of infantile bilateral striatal necrosis, determining the degeneration of the basal ganglia in humans [155].

    7. Nuclear Architecture and Gene Expression Regulation: The Nucleolus

    Although heterochromatic satellite regions were described near the nucleolus in the 1930s the role of the nucleolus as chromatin organizer is a recent hypothesis. Two independent high-throughput sequencing studies have isolated and sequenced the genomic DNA associated with purified nucleoli [156,157], finding the nucleolar-associated domains (NADs) (Figure 1). These genomic regions are gene poor and enriched for satellite DNA repeats. Interestingly, NAD peaks overlap with previously published LADs [156] and live imaging-based experiments have revealed an exchange between nuclear lamin and nucleolus border after mitosis. These data were confirmed by an independent study done by tracking LADs during the cell cycle [108]. Intriguingly, while Lamin B tends to interact only with LADs at the nuclear periphery, the subfraction of Lamin A that localize in the nucleoplasm preferentially accumulates around the nucleoli [117]. This suggests an interplay between the nucleolus and the nuclear membrane. Taken together these findings suggest that the nucleolus, as well as the nuclear envelope, could be involved in various nuclear functions, including gene transcription and the maintenance of chromatin structure. In line with this hypothesis, two seminal papers have provided insights into independent nucleolar mechanisms regulating pluripotency in embryonic stem cells (ESCs) [158,159]. In the first paper the authors found that the long noncoding RNA, pRNA, determines the progressive condensation of euchromatin into heterochromatin at nucleolus reducing epigenetic plasticity during differentiation [159]. In the second work authors showed that the methylation of the key pluripotency factor LIN28A leads to nucleolar localization and sequestration, preventing LIN28A processing and maintaining stemness programs [158]. These studies support the view that nuclear components are an important scaffold for spatial genome organization, which bears direct relevance to the establishment of cell type specific gene expression programs. Further studies are necessary to elucidate how the nuclear structures work together in a coordinated manner.

    8. When the Nuclear Architecture Works Badly

    In recent years, what is emerging is that, besides the plasticity of the chromatin fundamental for fine-regulated process, the nuclear architecture can also influence important cellular processes and is a hallmark of the healthy cell. Indeed an increasing number of tissue-specific pathologies in humans and tissue-specific phenotypes in model organisms have been described for mutations in a variety of the nuclear architecture components (reviewed in [160]) (Table 1).

    Table 1. Diseases associated with mutations in nuclear architecture components. The table reports some of the most importantpathologies related to mutations in lamins and nuclear pore complex (NPC), two main component of the nuclear envelope.
    NUCLEAR ARCHITECTURE COMPONENTSMUTATIONSASSOCIATED DISEASEREFERENCE
    Lamin A/C1824 C Cryptic splice variant at LMNA gene < THutchinson-Gilford Progeria Syndrome (HGPS)[130]
    Q6X, R25G/P, R50S/P, R133P, S143F, H222P, R249Q, Y267C, M371K, R377H/L, R453W, W498R, R527P/H, L530P, R644CEmery Dreifuss Muscular dystrophy (EDMD)[134]
    S143P, E161K, R190W/Q, N195K, R377H/L, R644CDilated Cardiomyopathy (DCM)[134]
    R482W/Q, K486N, R60G, G62G, R133LFamilial Partial Lipodystrophy 2 (FPLD2)[134]
    Lamin B1Partial duplication of LMNB1 geneLeukodystrophy (ADLD)[123]
    Lamin B2R215Q, A407TAcquired Partial Lipodystrophy (APL)[124]
    Nup155R391HAtrial fibrillation (AF)[154]
    Nup62Q391PInfatile bilateral striatal necrosis[155]
     | Show Table
    DownLoad: CSV

    During cancerogenesis the nuclear shape is used as parameter to identify the grade of malignancy (reviewed in [161,162]); in spite of all molecular mechanisms found to be involved in cancer, cyto-pathologists still make their official diagnosis by basing it on nuclear morphological features. Thus, it is possible that a number of different mechanisms involved in cancer progression finally determine the same nuclear abnormalities which in turn have functional consequence. The other view is that the nuclear architecture can drive the chromosomal instability and the aberrant gene expression underlying cancer pathogenesis and progression. In line with this hypothesis, proteomics screening methods for the identification of diagnostic and prognostic biomarkers in prostate cancer showed an increased expression of Lamin A in higher risk tumours and indicated a link to the progression of the disease and tumor aggressiveness [163]. Moreover, expression of Lamin A in pre-metastatic colon adenocarcinoma cell line determines an increase of spindle-like morphology and motility, two parameters normally associated with metastatic behaviour [164]. Live-imaging studies on a prostatic cancer cell line have confirmed these results, showing unexpected dynamics of the nuclear envelope, with cells showing irregular nuclei and cells showing dynamic interphase deflections of the lamina [165,166]. Intriguingly, these nuclear alterations were not dependent on actin microfilaments or microtubules. Thus the authors suggest that an intranuclear force, possibly chromatin-based, actively deforms the interphase nuclear envelope in a subset of chromosomally unstable cells. On the other hand, it is now well known that cancer cells are strongly influenced by their biomechanical environment (reviewed in [167]). Correspondingly, the nuclear envelope, in particular Lamin A, can influence nuclear mechanical properties [168,169], being at the crossroad between nucleoskeleton and cytoskeleton. Taken together such evidence leads to the hypothesis that the nuclear envelope can directly mediate the crosstalk between the tumor environment and the epigenome. Further studies will elucidate if the unique microenvironment of a tumor can influence the epigenome through remodeling of the nuclear architecture and if is possible to counteract tumor progression by reverting nuclear structure abnormalities.

    9. The Role of the Nuclear Architecture in DNA Repair

    9.1. An overview on DNA repair

    When DNA damage happens, distinct pathways are activated to repair depending on the type of damage [170]. DNA damage involving only one strand is repaired through Base Excision Repair (BER) or Nucleotide excision repair (NER), using undamaged strand as a template for repair [170]. Alternatively, cells activate Mismatch Repair (MMR) that repairs mismatch, deletion or insertion which all occur during recombination or replication [170]. The DNA could also be subjected to the formation of adducts such as alkylation that is toxic for the cell. These lesions are repaired “directly” through DNA repair by reversal of DNA damage [171]. The double strand break (DSB), in mammalian cells are repaired with non-homologous end-joining (NHEJ) and homologous recombination (HR) (reviewed in [172]) (Figure 3). NHEJ doesnot require a homologous template and performing the religation of the DNA ends can be error prone, making errors caused by the end processing [172]. Indeed, during NHEJ, the ends of damaged DNA are recognized and bound by the heterodimer Ku70/Ku80, which in turn recruits DNA-PKcs, stimulating its kinase activity. DNA-PKcs autophosporylates, so phosphorylating histone H2AX(γH2AX). This can be also phosphorylated directly by ATM and or ATR [173], two key sensors for DNA damage, which can activate the checkpoint of DSB repair in distinct manner and timing [174]. These events stimulate the recruitment of other downstream components of the NHEJ pathway including DNA ligase IV that religate the ends with the help of XRCC4 [175]. On the other hand, the homologous recombination (HR) is more accurate than NHEJ and requires a sister chromatid used as a template [170,172] (Figure 3). During homologous recombination, one filament of DNA is processed by several proteins to produce a 3’ overhang that is recognized and bound by RPA. Subsequently, Rad51 and Rad52 bind the single strand of DNA and search the homologous sequence. Then, the DNA/protein complex invades the duplex template and starts the branch migration, forming a displacement loop (D-loop) between the invading 3' overhang strand and the homologous chromosome. A DNA polymerase extends the 3’ end of the invading strand; this changes the D-loop to a Holliday junction. Finally a DNA-Ligase and Resolvase complete the process [176].

    Figure 3. DNA repair mechanisms. Schematic representation of the two main double strand break pathways: non-homologous end-joining (NHEJ) and the homologous recombination (HR). During NHEJ, the ends of damaged DNA are recognized and bound by the heterodimer Ku70/Ku80 that in turn recruits DNA-PKcs, that phosphorylates histone H2AX. DNA ends are finally joined by DNA ligase IV. During the HR, 3’ DNA overhangs are produced at DSB. ssDNA is recognized by RPA that recruits several factors such as Rad51 and Rad52. The DNA/proteins complex invades the duplex template and starts the branch migration, forming a displacement loop (D-loop) between the invading 3' overhang strand and the homologous chromosome. A DNA polymerase extends the 3’ end of the invading strand forming a Holliday junction and finally a DNA-Ligase with a Resolvase complete the process.

    9.2. DNA repair and the nuclear environment

    Emerging findings have revealed that the DNA repair is another cellular process that can be influenced by the nuclear position. Tethering a DNA sequence containing a double strand break (DSB) at nuclear lamina determines a delay in DNA Damage Response (DDR) and a repression of the homologous recombination (HR) repair pathway, whereas the non-homologous end-joining (NHEJ) is not influenced [177]. On the other hand, the forced localization of the same DNA sequence at nuclear pores did not give the same results, confirming again the difference between the two peripheral compartments. At the genome wide level, gene-rich regions and euchromatin marked by H3K36me3 are preferentially repaired by HR [177,178] while, when damaged, the heterochromatin associated with LADs, characterized by enrichment in A/T repeats undergoes the NHEJ rather than HR [107]. This choice probably is to avoid the genomic instability caused by the recombination of A/T repeats mediated by HR. These observations suggest that the grade of organization, and the state of chromatin in term of condensation and transcribed gene content, could partly regulate the choice of repair pathway after DSB. Parallel studies in yeast have confirmed these findings, further showing that different perinuclear anchorage sites define distinct outcomes on repair, arguing that the spatial segregation of damage at the nuclear envelope participates selectively in the choice of survival pathways [179]. The mobility of the chromatin after the DSB could be an important determinant of different DNA repair choice (reviewed in [180]). In fact through single-particle tracking it has been shown that DSBs move within a larger radius than does its undamaged counterpart [181]. This depends on chromatin remodelers [181,182] and proteins involved in HR pathway [181,183]. However, the chromatin conformational changes and epigenetic signatures involved in these processes need to be elucidated. Intriguingly, Lamin A seems to limit the DSB movement in mammalian cells [184]. This confirms that nuclear lamina, closely associated to heterochromatin, favours the NHEJ and further suggesting that Lamin A and the nuclear lamina can influence the DNA repair choice through regulation of chromatin dynamics and structure. Recently, a role for chromatin structure in DSB repair has been described also at the nucleolus level [185]. When DSB occurs at ribosomal gene (rDNA) buried inside nucleolus a nucleolar re-organization was observed, this allowing the damaged DNA to be relocated at the nucleolar periphery. Reorganization renders rDNA accessible to repair factors normally excluded from nucleoli. Importantly, this chromatin re-organization is thought to predominantly allow the HR repair rather than NHEJ. The exact factors driving this choice remain unknown. One possibility is that the suppression of NHEJ in nucleolar damage is dependent on the acetylation of histone H4 on K16 which inhibits the recognition of H4K20me2 by 53BP1, an important protein that promotes NHEJ [186,187]. These findings further support the view that the “chromatin remodelers” acting on histone modifications and chromatin organization are able to dictate the DNA repair choice, enlarging the roles of chromatin modifiers also in this field [188]. On the other hand, the nuclear components, such as lamins or nucleoporins, probably through their influence on the epigenome, can also determine the DNA repair choice. In line with this idea, Hutchinson-Gilford Progeria Syndrome (HGPS) patient-derived cells with nuclear shape abnormalities as described above, show DNA repair defects. In particular, HGPS fibroblasts accumulate DNA damage, showing increased levels of the marker of DNA damage such as the phosphorylated histone variant γH2AX and an impaired recruitment on DNA damage of 53BP1, Rad51 [172,173,174], MRE11 and phospho NBS1 [189], key components of HR repair pathway (Figure 4). These findings suggest that both the HR and the NHEJ pathway could be affected in HGPS. It has been proposed that the accumulation of DNA damage in progeria could depend on the accumulation of reactive oxygen species (ROS), produced by mitochondria and NADPH oxidases. This leads to oxidative stress possibly due to chromatin alterations [190]. In line with this hypothesis it has been shown that while DNA damage induced by the etoposide could be repaired in progeroid fibroblasts, ROS-induced DSB are unrepairable; this suggests that the accumulation of DNA damage could be mainly due to higher ROS levels in progeroid fibroblasts [191]. In response to oxidative stress, the gene expression is physiologically modulated in order to activate various redox-transcription factors, such as Rb, p53, NF-kB and FoxO (reviewed in [192]), fundamental for the resistance to oxidative stress. Accumulation of farnesylated prelamin A in progeric laminopathies results in increased production of ROS, leading respectively to senescence and apoptosis and to a reduction of antioxidant enzymes [193,194,195](Figure 4); this indicates that, under normal condition, lamins could act as recruiters and/or stabilizers also for antioxidant enzyme [190]. In line with this hypothesis, in Atypical Werner Syndrome, another diseases associated with Lamin A/C mutations [196,197] the level of the antioxidant enzymes is lower than normal [198]. Such evidence underlines the role of lamins as a shield against the ROS and shows that aberrant lamins lead to a lowering of this defence so adding another possible cause to the aging phenotype of laminopathies.

    Figure 4. DNA repair in Hutchinson-Gilford progeria syndrome. The figure compares DNA damage response in normal (A) and progeria fibroblasts (B). A typical marker of DNA damage, γH2AX (yellow dots) is indicated. In wt cells double strand breaks (DSB) are preferentially processed by NHEJ in the nuclear periphery or by HR in the nuclear interior (A). In progeric fibroblasts the accumulation of progerin determines an altered nuclear envelope leading to an increase of genome instability. This determines an accumulation of DNA damage and γH2AX due to the displacement of NHEJ and HR pathway components. On the other hand, the accumulation of prelamin A in progeric cells also determines an increase in reactive oxygen species (ROS) and a decrease in antioxidant enzymes, so favouring the formation of additional DSBs.

    10. Conclusions

    As summarized in this review, the organization of the nucleus and the compartmentalization of chromatin are fundamental for gene expression regulation and, as recent studies show, can both influence specific cellular processes such as DNA repair. Recently, DNA integration has been indicated as another mechanism that could be dependent on the nuclear architecture. Previous studies on Human immunodeficiency virus 1 (HIV-1) have already suggested a preference for an integration within actively expressed chromosomal locations [199]. However it was unclear whether transcription directly or indirectly helps integration (reviewed in [200]). Recently, it has been shown that HIV-1 preferentially integrates in genes localized at the periphery, close to the nuclear pores and characterized by active transcription chromatin marks before viral infection [201]. After integration, the viral DNA needs the activity of the nucleoporin Nup 153 for transcriptional activation while, for replication purposes, it needs the activity of Tpr, also located at the nuclear basket site [201,202]. Overall such evidence promotes a new concept supporting the role of nuclear architecture in HIV integration. Thus it is tempting to speculate that other cellular processes (such as DNA replication, cellular migration or the stress response) can also be dependent on the nuclear architecture and on its influence on epigenome regulation.

    Acknowledgments

    We are grateful to Prof. Martin Bennett for “proofreading work” of the manuscript. This work was supported by grants from the Italian Ministry of Research and University (Futuro in Ricerca RBFR106S1Z_001) and the flagship CNR project (Epigen). Authors declare that there are no conflicts of interest.

    Conflict of Interest

    All authors declare no conflicts of interest in this paper.



    [1] R. M. Anderson and R. M. May, Infectious Diseases of Humans, Dynamics and Control, Oxford University, Oxford,1991.
    [2] V. Capasso, Mathematical Structures of Epidemic Systems, Lecture Notes in Biomathematics, vol. 97, Springer-Verlag, Berlin, 1993.
    [3] H. W. Hethcote, The mathematics of infectious diseases, SIAM Rev., 42 (2000).
    [4] O. Diekmann and J. A. P. Heesterbeek, Mathematical Epidemiology of Infectious Diseases. Model Building, Analysis and Interpretation, Wiley Ser. Math. Comput. Biol., John Wiley and Sons, Chichester, England, 2000.
    [5] F. Brauer and P. van den Driessche, Models for the transmission of disease with immigration of infectives, Math. Biosci., 171 (2001), 143–154.
    [6] M. A. Nowak and A. R. McLean, A mathematical model of vaccination against HIV to prevent the development of AIDS, Proc. Biol. Sci., 246 (1991), 141–146.
    [7] N. M. Ferguson, D. A. T. Cummings, S. Cauchemez, et al., Strategies for containing an emerging in fuenza pandemic in Southeast Asia, Nature, 437 (2005), 209–214.
    [8] Y. N. Xiao, T. T. Zhao and S. Y. Tang, Dynamics of an infectious disease with media/ psychology induced non-sooth incidence, Math. Biosci. Eng., 10 (2013), 445–461.
    [9] A. L. Wang and Y. N. Xiao, A Filippov system describing media effects on the spread of infectious diseases, Nonlinear Anal. Hybri., 11 (2014), 84–97.
    [10] J. A. Cui, X. X. Mu and H. Wan, Saturation recovery leads to multiple endemic equilibria and backward bifurcation, J. Theor. Biol., 254 (2008), 275–283.
    [11] J. L. Wang, S. Q. Liu, B. W. Zheng, et al., Qualitative and bifurcation analysis using an SIR model with a saturated treatment function, Math. Comput. Model., 55 (2012), 710–722.
    [12] X. Zhang and X. N. Liu, Backward bifurcation of an epidemic model with saturated treatment function, J. Math. Anal. Appl., 348 (2008), 433–443.
    [13] L. H. Zhou and M. Fan, Dynamics of an SIR epidemic model with limited medical resources revisited, Nonlinear Anal. Real., 13 (2012), 312–324.
    [14] H. Wan and J. A. Cui, Rich dynamics of an epidemic model with saturation recovery, J. Appl. Math., 314958 (2013).
    [15] Z. H. Zhang and Y. H. Suo, Qualitative analysis of a SIR epidemic model with saturated treatment rate, J. Appl. Math. Comput., 34 (2010), 177–194.
    [16] S. Y. Tang and R. A. Cheke, Stage-dependent impulsive models of integrated pest management (IPM) strategies and their dynamic consequences, J. Math. Biol., 50 (2005), 257–292.
    [17] S. Y. Tang, Y. N. Xiao and R. A. Cheke, Dynamical analysis of plant disease models with cultural control strategies and economic thresholds, Math. Comput. Simul., 80 (2010), 894–921.
    [18] A. B. Sabin, Measles, killer of millions in developing countries: strategies of elimination and continuing control, Eur. J. Epidemiol., 7 (1991), 1–22.
    [19] B. Shulgin, L. Stone and Z. Agur, Pulse vaccination strategy in the SIR epidemic model, Bull. Math. Biol., 60 (1998), 1123–1148.
    [20] L. Stone, B. Shulgin and Z. Agur, Theoretical examination of the pulse vaccination policy in the SIR epidemic model, Math. Comput. Model., 31 (2000), 207–215.
    [21] A. D'Onofrio, On pulse vaccination strategy in the SIR epidemic model with vertical transmission, Appl. Math. Lett., 18 (2005), 729–732.
    [22] F. L. Black, Measles endemicity in insular populations: Critical community size and its evolutionary implication, J. Theor. Biol., 11 (1966), 207–211.
    [23] W. J. Moss and P. Strebel, Biological feasibility of measles eradication, J. Infect. Dis., 204 (2011), 47–53.
    [24] Q. Q. Zhang, B. Tang and S. Y. Tang, Vaccination threshold size and backward bifurcation of SIR model with state-dependent pulse control, J. Theor. Biol., 455 (2018), 75–85.
    [25] Q. Li and Y. N. Xiao, Dynamical behaviour and bifurcation analysis of the SIR model with continuous treatment and state-dependent impulsive control, Int. J. Bifurcat. Chaos, (2019), Accepted.
    [26] A. L. Wang, Y. N. Xiao and R. Smith?, Using non-smooth models to determine thresholds for microbial pest management, J. Math. Biol., (2019).
    [27] A. L. Wang, Y. N. Xiao and R. Smith?, Multiple equilibria in a non-smooth epidemic model with medical-resource constraints, Bull. Math. Biol., 81 (2019), 963–994.
    [28] W. J. Qin, S. Y. Tang and R. A. Cheke, Nonlinear pulse vaccination in an SIR epidemic model with resource limitation, Abstr. Appl. Anal., 670263 (2013).
    [29] J. Yang and S. Y. Tang, Holling type II predator–prey model with nonlinear pulse as state-dependent feedback control, J. Comput. Appl. Math., 291 (2016), 225–241.
    [30] S. Y. Tang, B. Tang, A. L. Wang, et al., Holling II predator-prey impulsive semi-dynamic model with complex Poincaré map, Nonlinear Dyn., 81 (2015), 1575–1596.
    [31] S. Y. Tang, W. H. Pang, R. A. Cheke, et al., Global dynamics of a state-dependent feedback control system, Adv. Diff. Equ., 322 (2015), DOI: 10.1186/s13662-015-0661-x.
    [32] J. M. Grandmont, Nonlinear difference equations, bifurcations and chaos: an introduction, Research in Economics, 62 (2008), 122–177.
    [33] Y. N. Xiao, X. X. Xu and S. Y. Tang, Sliding mode control of outbreaks of emerging infectious diseases, Bull. Math. Biol., 74 (2012), 2403–2422.
    [34] Y. P. Yang, Y. N. Xiao and J. H. Wu, Pulse HIV vaccination: feasibility for virus eradication and optimal vaccination schedule, Bull. Math. Biol., 75 (2013), 725–751.
    [35] Y. Tian, K. B. Sun, A. Kasperski, et al., Nonlinear modelling and qualitative analysis of a real chemostat with pulse feeding, Discrete Dyn. Nat. Soc., 640594 (2011), 1–18.
    [36] S. Y. Tang and W. H. Pang, On the continuity of the function describing the times of meeting impulsive set and its application, Math. Biosci. Eng., 14 (2017), 1399–1406.
    [37] X. N. Liu, Y. Takeuchi and S. Iwami, SVIR epidemic model with vaccination strategy, J. Theor. Biol., 253 (2008), 1–11.
  • This article has been cited by:

    1. Vladimir B. Teif, Andrey G. Cherstvy, Chromatin and epigenetics: current biophysical views, 2016, 3, 2377-9098, 88, 10.3934/biophy.2016.1.88
    2. F. Marullo, E. Cesarini, L. Antonelli, F. Gregoretti, G. Oliva, C. Lanzuolo, Nucleoplasmic Lamin A/C and Polycomb group of proteins: An evolutionarily conserved interplay, 2016, 7, 1949-1034, 103, 10.1080/19491034.2016.1157675
    3. Federica Marasca, Fabrizia Marullo, Chiara Lanzuolo, 2016, Chapter 15, 978-1-4939-6378-2, 167, 10.1007/978-1-4939-6380-5_15
    4. Andrea Bianchi, Pierluigi Giuseppe Manti, Federica Lucini, Chiara Lanzuolo, Mechanotransduction, nuclear architecture and epigenetics in Emery Dreifuss Muscular Dystrophy: tous pour un, un pour tous, 2018, 9, 1949-1034, 321, 10.1080/19491034.2018.1460044
    5. Endre Sebestyén, Fabrizia Marullo, Federica Lucini, Cristiano Petrini, Andrea Bianchi, Sara Valsoni, Ilaria Olivieri, Laura Antonelli, Francesco Gregoretti, Gennaro Oliva, Francesco Ferrari, Chiara Lanzuolo, SAMMY-seq reveals early alteration of heterochromatin and deregulation of bivalent genes in Hutchinson-Gilford Progeria Syndrome, 2020, 11, 2041-1723, 10.1038/s41467-020-20048-9
    6. Austė Kanapeckaitė, Neringa Burokienė, Asta Mažeikienė, Graeme S. Cottrell, Darius Widera, Biophysics is reshaping our perception of the epigenome: from DNA-level to high-throughput studies, 2021, 1, 26670747, 100028, 10.1016/j.bpr.2021.100028
    7. Philina Santarelli, Valentina Rosti, Maria Vivo, Chiara Lanzuolo, 2024, 158, 9780128234358, 375, 10.1016/bs.ctdb.2024.01.014
    8. Maria Vivo, Valentina Rosti, Sara Cervone, Chiara Lanzuolo, Chromatin plasticity in mechanotransduction, 2024, 88, 09550674, 102376, 10.1016/j.ceb.2024.102376
  • Reader Comments
  • © 2019 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(5309) PDF downloads(518) Cited by(3)

Article outline

Figures and Tables

Figures(8)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog