Research article Special Issues

Cross immunity protection and antibody-dependent enhancement in a distributed delay dynamic model

  • Dengue fever is endemic in tropical and subtropical countries, and certain important features of the spread of dengue fever continue to pose challenges for mathematical modelling. Here we propose a system of integro-differential equations (IDE) to study the disease transmission dynamics that involve multi-serotypes and cross immunity. Our main objective is to incorporate and analyze the effect of a general time delay term describing acquired cross immunity protection and the effect of antibody-dependent enhancement (ADE), both characteristics of Dengue fever. We perform qualitative analysis of the model and obtain results to show the stability of the epidemiologically important steady solutions that are completely determined by the basic reproduction number and the invasion reproduction number. We establish the global dynamics by constructing a suitable Lyapunov functional. We also conduct some numerical experiments to illustrate bifurcation structures, indicating the occurrence of periodic oscillations for a specific range of values of a key parameter representing ADE.

    Citation: Vanessa Steindorf, Sergio Oliva, Jianhong Wu. Cross immunity protection and antibody-dependent enhancement in a distributed delay dynamic model[J]. Mathematical Biosciences and Engineering, 2022, 19(3): 2950-2984. doi: 10.3934/mbe.2022136

    Related Papers:

    [1] Nadwan Majeed Ali, Handri Ammari . Design of a hybrid wind-solar street lighting system to power LED lights on highway poles. AIMS Energy, 2022, 10(2): 177-190. doi: 10.3934/energy.2022010
    [2] David Ludwig, Christian Breyer, A.A. Solomon, Robert Seguin . Evaluation of an onsite integrated hybrid PV-Wind power plant. AIMS Energy, 2020, 8(5): 988-1006. doi: 10.3934/energy.2020.5.988
    [3] Virendra Sharma, Lata Gidwani . Recognition of disturbances in hybrid power system interfaced with battery energy storage system using combined features of Stockwell transform and Hilbert transform. AIMS Energy, 2019, 7(5): 671-687. doi: 10.3934/energy.2019.5.671
    [4] Zeid Al Qaisi, Qais Alsafasfeh, Ahmad Harb . Stability impact of integrated small scale hybrid (PV/Wind) system with electric distribution network. AIMS Energy, 2018, 6(5): 832-845. doi: 10.3934/energy.2018.5.832
    [5] Rashid Al Badwawi, Mohammad Abusara, Tapas Mallick . Speed control of synchronous machine by changing duty cycle of DC/DC buck converter. AIMS Energy, 2015, 3(4): 728-739. doi: 10.3934/energy.2015.4.728
    [6] Aaron St. Leger . Demand response impacts on off-grid hybrid photovoltaic-diesel generator microgrids. AIMS Energy, 2015, 3(3): 360-376. doi: 10.3934/energy.2015.3.360
    [7] Virendra Sharma, Lata Gidwani . Optimistic use of battery energy storage system to mitigate grid disturbances in the hybrid power system. AIMS Energy, 2019, 7(6): 688-709. doi: 10.3934/energy.2019.6.688
    [8] Amer Braik, Asaad Makhalfih, Ag Sufiyan Abd Hamid, Kamaruzzaman Sopian, Adnan Ibrahim . Impact of photovoltaic grid-tied systems on national grid power factor in Palestine. AIMS Energy, 2022, 10(2): 236-253. doi: 10.3934/energy.2022013
    [9] Tilahun Nigussie, Wondwossen Bogale, Feyisa Bekele, Edessa Dribssa . Feasibility study for power generation using off- grid energy system from micro hydro-PV-diesel generator-battery for rural area of Ethiopia: The case of Melkey Hera village, Western Ethiopia. AIMS Energy, 2017, 5(4): 667-690. doi: 10.3934/energy.2017.4.667
    [10] Yassine Charabi, Sabah Abdul-Wahab . The optimal sizing and performance assessment of a hybrid renewable energy system for a mini-gird in an exclave territory. AIMS Energy, 2020, 8(4): 669-685. doi: 10.3934/energy.2020.4.669
  • Dengue fever is endemic in tropical and subtropical countries, and certain important features of the spread of dengue fever continue to pose challenges for mathematical modelling. Here we propose a system of integro-differential equations (IDE) to study the disease transmission dynamics that involve multi-serotypes and cross immunity. Our main objective is to incorporate and analyze the effect of a general time delay term describing acquired cross immunity protection and the effect of antibody-dependent enhancement (ADE), both characteristics of Dengue fever. We perform qualitative analysis of the model and obtain results to show the stability of the epidemiologically important steady solutions that are completely determined by the basic reproduction number and the invasion reproduction number. We establish the global dynamics by constructing a suitable Lyapunov functional. We also conduct some numerical experiments to illustrate bifurcation structures, indicating the occurrence of periodic oscillations for a specific range of values of a key parameter representing ADE.



    Abbreviations:

    PKS: polyketide synthase

    NRPS: nonribosomal polypeptide-synthetase

    PT: phylotype

    HPLC: high pressure liquid chromatography

    MS: molecular mass

    MeCN: acetonitrile

    nt: nucleotide


    1. Introduction

    Members of the Actinomycetales are widely distributed in terrestrial as well as in aquatic and also marine ecosystems, where they are supposed to play an important role in the decomposition and recycling of biomaterials [1]. Maldonado et al. described first obligate marine species of the new genus Salinispora: Salinispora arenicola and Salinispora tropica [10]. Furthermore, two new genera, Salinibacterium and Serinicoccus were described containing solely marine isolates that tolerate high concentrations of sodium chloride [11,12]. They even have been recovered from the world’s deepest ocean trench, the Mariana Trench [2,3,4]. A controversial discussion exists about the question whether they are metabolically active in marine sediments or have been washed into the marine environment where they are just prevailing as spores. Culture experiments using marine salinity showed that many strains grow very well under these conditions and therefore appear well adapted to the marine environment [5,6]. In fact, evidence accumulates that they are active members of the microbial community in marine sediments [6,7,8,9].

    A specific feature of the actinomycetes is their enormous contribution to the production of chemically diverse secondary metabolites and clinically relevant antibiotics. Indeed, the majority of bioactive compounds of microbial origin originates from actinomycetes and in particular from representatives of the genera Streptomyces and Micromonospora [13,14]. Novel biological active compounds have also been found in actinomycetes isolated from the marine realm [1,15,16,17]. These findings are a clear indication that actinomycetes from the marine environment may represent a valuable source for the discovery of novel natural products which has so far not been intensively studied. The genus Micromonospora was described by Orskov in 1923 and is represented by a total of 61 species and 7 subspecies in 2016 [18,19]. The first antibiotic isolated from a Micromonospora strain was micromonosporin [20]. After the discovery of the aminoglycoside antibiotic gentamicin in 1963 an intensive screening for new natural products in Micromonospora strains started [21]. Today more than 400 compounds are known that have been isolated from members of the genus Micromonospora [22]. The great potential of marine Micromonospora isolates to produce bioactive natural products was also highlighted by a study of marine actinomycetes isolated from Pacific sediments [23] which demonstrated strong anti-tumor activity of several culture extracts of Micromonospora strains. Cytotoxic activity against several human tumor cell lines was also observed for the levantilides, 20-membered cytotoxic macrolides produced by a deep-sea Micromonospora strain isolated from a Mediterranean deep-sea sediment sample in a previous study [24].

    We have studied the bacterial diversity in the eastern basin of the Mediterranean Sea, which is one of the most oligotrophic regions of the world’s oceans [28]. The extreme depletion of nutrients, in particular phosphorus, results in low primary production in waters distant from the coast [25]. Based on the poor primary production, minor amounts of organic matter reach the deep-sea floor of the Eastern Mediterranean Sea, in particular in one of the deep basins, the Ierapetra Basin, which is located approximately 30 nautic miles southeast of Crete. Studies on the macrofauna community showed that sediments in the Eastern Mediterranean Sea can be considered as event-driven ecosystems and the biota need to be adapted to few and pulsed nutrient input events [26]. In addition, the Mediterranean deep sea is characterized by a comparably high temperature, which may exceeds 13 °C in deep waters [27]. The unique and extreme conditions of the Eastern Mediterranean deep sea require specific adaptation of its inhabitants and therefore these sediments are an attractive source for the study of actinomycetes. In a previous study on the cultured biodiversity from Mediterranean deep sea sediments, representatives of Bacillus and Actinobacteria were found asmajor groups, the Actinobacteria representing almost a third of all isolated bacteria [28]. This result animated to perform the present study in which a highly diverse group of Actinobacteria were selectively enriched from Mediterranean sediments (4400 m and 2800 m deep). In particular, members of the genus Micromonospora, which are rarely described from deep-sea habitats, were identified and their potential to produce antibiotic compounds was studied.


    2. Materials and Method


    2.1. Sample collection

    All samples were obtained during the Meteor research cruise 71 leg 2 southwards of Crete in January 2007. Sediment samples from two different sites at 4400 m (Ierapetra Basin; 34°30.296N, 26°11.507E) and 2800 m depth (Herodotos Plain; 33°42.989N, 26°20.329E) were collected using a multiple corer. The uppermost 5 cm of each sediment core were aseptically sub-sampled (from the Herodotos Plain also a sub-sample at 25 cm depth) and stored at 4 °C until further treatment.


    2.2. Isolation of Actinobacteria

    15 Micromonospora strains as well as one representative each of the genera Knoellia (A207) and Marmoricola (A214), respectively, were isolated after pre-treatments with pressure incubation and/or heat treatment of the sediment samples (Table 1).

    Table 1. Cultivation conditions and phylogenetic analyses of Micromonospora strains from the Eastern Mediterranean deep sea sediment
    Cluster PT Strain 16S rRNA gene sequence lenght accession [nt] number Next related type strainand accession number Similarity (%) Pressure pretreat-ment Heat treatment Isolationmedium (temperature) Sampleorigin
    A 12 S61 1414 FM992759a) M. pattaloongensis JCM 12883T AB275607 98.8 no no Chitin (25 °C) HP
    B 5 S20 1388 FM992749a) M. auratinigra DSM 44815T AB159779 99.4 no no CFB (25 °C) HP
    B 5 S51 1349 FM992776a) M. auratinigra DSM 44815T AB159779 99.3 no no CFB (25 °C) IB
    B 5 A77 1351 FR714833b) M. auratinigra DSM 44815T AB159779 99.3 no yes XJ4 (25 °C) IB
    C 11 S32a 1361 FM992754a) M. tulbaghiae DSM 45142T EU196562 100.0 no no MW (25 °C) HP
    D 13 S71 1437 FM992746a) M. chaiyaphumensis JCM 12873T AB196710 99.7 no no MWY (25 °C) IB
    E 1 S29 1417 FM992751a) M. purpureochromogenes DSM 43821T X92611 98.6 no no CFB (25 °C) IB
    E 4 A68 1353 KX256200c) M. purpureochromogenes DSM 43821T X92604 98.6 yes yes XJ4 (25 °C) IB
    E 6 A75a 1384 KX256202c) M. purpureochromogenes DSM 43821T X92611 98.8 no yes XJ4 (25 °C) IB
    E 6 A172 1472 KX256204c) M. purpureochromogenes DSM 43821T X92611 98.9 yes yes XJ4 (13.5 °C) IB
    E 6 A203 1443 KX256210c) M. purpureochromogenes DSM 43821T X92611 98.7 no yes XJ4 (13.5 °C) IB
    E 6 A165 1371 KX256203c) M. purpureochromogenes DSM 43821T X92611 98.7 yes yes MWY (25 °C) IB
    E 7 A177 1474 KX256206c) M. purpureochromogenes DSM 43821T X92611 98.3 no yes XJ4 (25 °C) IB
    E 8 A200 1476 KX256207c) M. purpureochromogenes DSM 43821T X92611 98.9 no yes HA (13.5 °C) IB
    E 8 A208 1471 KX256212c) M. purpureochromogenes DSM 43821T X92611 99.1 yes yes HA (13.5 °C) IB
    E 10 A176 1471 KX256205c) M. purpureochromogenes DSM 43821T X92611 98.8 no yes XJ4 (25 °C) IB
    E 14 A201 1097 KX256208c) M. purpureochromogenes DSM 43821T X92611 99.7 no yes XJ4 (13.5 °C) IB
    F 2 D117x 1342 FM992779c) M. saelicesensis DSM 44871T AJ783993 99.0 yes no CFB (25 °C) IB
    G 3 A74 1363 KX256201c) M. saelicesensis DSM 44871T AJ783993 99.5 no yes XJ4 (25 °C) IB
    G 3 A64 1472 KX256199c) M. saelicesensis DSM 44871T AJ783993 99.6 no yes HA (25 °C) HP
    G 9 A202 1437 KX256209c) M. saelicesensis DSM 44871T AJ783993 99.2 no yes XJ4 (13.5 °C) IB
    IB: Ierapetra Basin (4400 m); HP; Herodotos Plain (2800 m); PT: phylotype; a) from [28]; b) from [24]; c) this study;
     | Show Table
    DownLoad: CSV

    Pressure incubation. In order to selectively enrich deep-sea bacteria (not solely Actinomycetes) that are adapted to elevated hydrostatic pressure and spontaneous nutrient input, sediment samples were pre-incubated at in situ pressure and in situ temperature in sea water supplemented with N-acetyl-D-glucosamine prior to plating on agar media. For this purpose 10 ml of sediment were transferred into a plastic bag with addition of 9 ml of sterile-filtered Mediterranean seawater (taken from the multicorer) and 1 ml N-acetyl-D-glucosamine solution with a concentration of 100 µM. N-acetyl-D-glucosamine (the monomer of chitin) is supposed to be an important carbon and nitrogen source for the microbial deep-sea community. The samples were incubated in compression-proof steel tubes at 280 bar (samples from the Herodotos Plain) or 440 bar (samples from the Ierapetra Basin) hydrostatic pressure and at 13.5 °C for 6 d.

    Heat treatment. For selective isolation of Actinobacteria and other spore-forming bacteria, 2 g of the sediment samples were dried in a sterile Petri dish at room temperature for 3 months. Dry sediment was then heated to 120 °C for 1 hour and re-suspended in 19 ml deionized water. 200 µl of this suspension were plated on the isolation media.

    Combination of both pre-treatments. A few bags were dried after the pressure incubation and in addition exposed to the heat treatment.

    The following isolation media were used: a) CFB-medium: 15 g agar, 1 g tryptone, 0.5 g yeast extract, 0.5 g CaCl2·2H2O, 0.5 g MgCl7H2O in 1 L Mediterranean seawater; b) XJ4-medium: 15 g agar, 0.1 g histidine, 1 g raffinose, 0.5 g NaHPO4, 1.7 g KCl, 0.05 g MgSO4·4H2O, 0.01 g FeSO7H2O, 0.02 g CaCO3, 0.5 mg each of thiamine, riboflavine, niacine, pyridoxin, Ca-panthotenate, inositol and p-aminobenzoic acid, 0.25 mg biotin and 50 ppm K2Cr2O7 in 1 L of aquadest with addition of cycloheximide (50µg/ml); c) HA-medium: 1 g humic acid, 0.01 g FeSO4, 0.5 g Na2HPO4, 1.7 g KCl, 0.02 g CaCO3, 0.5 g asparagine and 50 ppm K2Cr2O7 in 1 L of deionized water; and d) MWY- medium: 15 g agar and 0.001 g yeast extract in 1 L Mediterranean seawater added with 50 µg/mL cycloheximide and penicillin G. Agar plates were incubated at 13.5 °C (in situ temperature) and 25 °C and checked regularly for growth. Using a binocular microscope the morphology of all colonies was checked and those appearing morphologically different were transferred to fresh agar medium until pure cultures were obtained. Isolates were stored at −80 °C using the Cryobank system (Mast Diagnostica GmbH, Reinfeld, Germany).

    With the aim to test the capability of all Micromonospora sp. strains to cope with the availability of low nutrient concentrations, they were cultivated on a medium consistent of 1L Mediterranean Sea water supplemented with 15 g agar at 13.5 °C.


    2.3. Phylogenetic analyses

    DNA extraction, 16S rRNA gene amplification and sequencing were performed as described previously [29]. 16S rRNA gene sequences obtained during this study, were submitted to the GenBank database and were assigned accession numbers KX256199 - KX256213.Sequences were compared to NCBI database sequences using BLAST search [30]. Sequence alignment of Micromonospora sp. isolates, next related type strains and further representatives of the genus was performed using the ARB database [31]. Phylogenetic trees were calculated applying maximum likelihood analysis with PhyML software [32]. This analysis was performed assuming the GTR evolution model with an optimized gamma distribution parameter alpha and 100 bootstrap replicates [33]. For systematic consideration, sequences with similarities higher than 99.5 % were grouped as a phylotype together using Mothur [34].


    2.4. Amplification of PKS-I, PKS-II and NRPS genes

    DNA extracts of Micromonospora strainswere analyzed for the presence genes for representative natural product biosynthesis pathways such as the polyketide synthase I and II (PKS-I and PKS-II) and non-ribosomal polypeptide synthetase (NRPS). The sequences of the respective primers are listed in Table 2. Detection of PKS-II genes was performed using the primers and PCR conditions as described by Mets-Ketel et al. [35]. For gene amplification, 50 pMol primer, 3.5 U TaqPolymerase, 10% DMSO and 5 µl template were applied. PKS-I and NRPS PCR assays were modified according to Ayuso-Sacido and Genilloud using PuReTaq® Ready-To-Go PCR Beads (GE Healthcare, Solingen, Germany) and 50 pMol of the primers K1 and M6R for the PKS-I gene [36]. The same system was used for amplification of NRPS genes with primers A3 and A7R and in addition primers A2f and A3r were used for amplification of NRPS according to the protocol of Martens et al., resulting in fragment length of 700 nt (NRPSlarge) and 300 nt, (NRPSshort), respectively [37]. PCR products were purified and sequenced using the BigDye Terminator v1.1 Sequencing Kit (Applied Biosystems, Life Technologies GmbH, Darmstadt, Germany) in a 3730-DNA-Analyzer (Applied Biosystems) as specified by the manufacturer. Sequences were edited by ChromasPro v.1.33, translated to amino acid sequences using BioEdit and compared to the NCBI database using BLAST search [30,38]. Matches of the blastx search were checked for sequence similarities and coverage regardingpreferably genes of Micromonospora sp.

    Table 2. Primers used in this study for the detection of PKS and NRPS gene fragments
    Name Sequence Target Sequence length [nt] Reference
    K1 5x-TSA AGT CSA ACA TCG GBC A-3x PKS-I 1200 - 1400 [36]
    M6R 5x-CGC AGG TTS CSG TAC CAG TA-3x
    PKS-KSαf 5x-TSG CST GCT TCG AYG CSA TC-3x PKS-II ~ 600 [35]
    PKS-KSαr 5x-TGG AAN CCG CCG AAB CCG CT-3x
    A3 5x-GCS TAC SYS ATS TAC ACS TCS GG-3x NRPSlong ~ 700 [36]
    A7R 5x-SAS GTC VCC SGT SCG GTA S-3x
    A2f 5x-AAG GCN GGC GSB GCS TAY STG CC-3x NRPSshort ~ 300 [37]
    A3r 5x-TTG GGB IKB CCG GTS GIN CCS GAG GTG-3x
     | Show Table
    DownLoad: CSV

    2.5. Production of culture extracts

    For HPLC-UV/MS analysis of the metabolite profiles and for antimicrobial test assays, culture extracts of all Micromonospora strains were prepared. Strains were cultivated in 100 ml liquid starch-peptone medium (1 g starch, 0.5 g soy-peptone and 1.5 g Tropic Marine® sea salt in 100 ml deionized water) using 300 ml Erlenmeyer-flasks with one baffle at 120 rpm. After incubation for 8 days at 28 °C 100 ml ethyl acetate (Picograde®, Promochem, Wesel, Germany) were added and the culture broth was homogenized using an Ultra-Turrax T25 basic disperser at 17500 rpm for 30 sec. The supernatants were separated, dried and re-suspended in 1 ml methanol.

    With the aim to study the effect of salt on the production of the compounds of interest, i.e. S20-4, S20-5. S20-6, and S20-7, culture extracts of strain S20 were prepared after growth in starch-peptone medium with Mediterranean Sea water (3.86%) and 30 g/L Tropic Marine® sea salt, respectively, for 8 days at 28 °C.


    2.6. HPLC analysis of culture extracts

    Metabolite profiles of culture extracts were screened by HPLC-UV/MS. For reversed-phase HPLC analysis, dried extracts were resolved in methanol. The measurements were carried out using a Phenomenex Onyx Monolithic C18 column (100 by 3.00 mm) on a VWR Hitachi Elite LaChrom system coupled to an electrospray ionization ion trap detector (Esquire 4000, Bruker Daltonics, Bremen). The solvents, water and acetonitrile (MeCN), were supplemented with 0.1% formic acid. The gradient of MeCN was as follows: 0 min 5 %; 4 min 60 %; 6 min 100 %; flow rate: 2ml/min).


    2.7. Fermentation of strain S20 and isolation of its metabolites

    Micromonospora strain S20 was cultivated in 10 L of liquid starch-peptone medium (10 g starch, 5 g soy-peptone, 15 g Tropic Marine® sea salt and 1 g CaCO3 in 1L deionized water) using 2 L shaken Erlenmeyer-flasks with one baffle and filled with 1 L medium each. After 6 days of incubation at 28 °C and 120 rpm, the culture supernatant was separated from the cells by centrifugation at 8.000 rpm for 5 min (Beckmann_J2-MC). Afterwards the supernatant (9.5 L) was acidified by adding 20 ml HCl (15%) and transferred to an Amberlite XAD16-column. After washing the column with 500 ml 40% ethanol bound metabolites were eluted with 2 L 80% (v/v) ethanol and subsequently with 2 L 100% (v/v) ethanol. Both eluates were collected and freeze-dried. The dried extracts were dissolved in 5 ml methanol each and analyzed by HPLC-UV/MS for the presence of the compounds S20-4, S20-5, S20-6, and S20-7. Since the desired substances were found in both extracts, they were combined for further analyses. The compounds were isolated and purified by preparative reversed phase HPLC using a Phenomenex Gemini-NX column (C18, 5 μ, 110A, 100 × 21.20 mm) and a VWR Hitachi Elite LaChrom system consisting of an L-7150 pump, an L-2450 diode array detector, and an L-2200 autosampler. A water (A) / methanol (B) gradient (50% B to 100% B in 17 min; flow 15 ml/min) with 0.1% formic acid added to both solvents was applied.


    2.8. Antimicrobialtests of culture extracts and metabolites

    Antimicrobial activity of the culture extracts of the Micromonospora strains as well as of the metabolites produced by the strain S20 (S20-4, S20-5, S20-6, S20-7) was tested against Bacillus subtilis (DSM 347), Staphylococcus lentus (DSM 6672), Xanthomonas campestris (DSM 2405), andCandida glabrata (DSM 6425) using a resazurin assay according to Nagel at al. [38]. The concentration of culture extracts in the bioassay was 0.01 mg/200 µl, the compounds were tested at concentrations of 100 µM. A negative (no compound) and a positive control (100 µM of chloramphenicol for bacterial test strains and 100 µM of cycloheximid for the yeast) were measured on the same microtiter plate. All tests were run in triplicates.


    3. Results


    3.1. Isolation of Actinobacteria

    Actinobacteria from Mediterranean deep-sea sediments were isolated after different pre-treatment procedures of the sediment samples, mostly from the uppermost 5 cm (Table 1). Among 23 actinobacterial strains 21 belonged to the genus Micromonospora and one strain each affiliated to Knoellia and Marmoricola.

    Ten Micromonospora strains as well as Marmoricola isolate A214 were obtained after treatment with dry heat on selective agar media. Marmoricola sp. strain A214 was isolated from the Herodotos Plain sediment after cultivation on XJ4-medium at 13.5 °C. It was the only isolate obtained from the 25 cm depth section of the sediment core. Four Micromonospora strains as well as Knoellia strain A207 were isolated after the combined pre-treatment approach, after incubation at in situ pressure with N-acetyl-D-glucosamine and subsequent heat treatment. The Knoellia sp. strain A207 was isolated originated from the Ierapetra Basin on HA-medium at 25 °C. Micromonospora sp. strain D117x was isolated after incubation at in situ pressure with enrichment on N-acetyl-D-glucosamine, without heat treatment. Six Micromonospora strains (S20, S29, S32a, S51, S61, and S71) were previously isolated from untreated sediment samples [28].


    3.2. Identification and phylogenetic relations of the isolates

    All 21 strains affiliating to the genus Micromonospora could be assigned to 14 different phylotypes (PTs) with an exclusion limit of ≥ 99.5% 16S gene sequence similarity (Table 1,Figure 1). The phylogenetic analysis of 16S rDNA gene sequences of the Micromonospora isolates revealed that most deep-sea strains were phylogenetically distinct from known reference organisms (Figure 1). These phylotypes formed at least seven distinct clusters, in part with a single isolate and phylotype representing the cluster (clusters A, C, D, F). Cluster B (strains S20, S51 and A77) affiliates with Micromonospora auratinigra as next relative, cluster G (strains A64, A74, A202) affiliates with M. saelicesensis and cluster E with M. purpureochromogenes (Table 1,Figure 1). Interestingly all strains of cluster G and cluster E (except strain S29) were isolated after heat treatment of the sediments and 4 out of 11 strains of cluster E in addition survived the pressure incubation. Despite their similarity to known reference species of Micromonospora, all isolates of clusters E, G and F (a single strain obtained after pressure treatment) formed separate branches in the phylogenetic tree (Figure 1). In addition to Micromonospora, members of rare actinomycetes were observed in the pre-treated sediment samples. Strain A207 affiliated to Knoelliasubterranea DSM 12372T (AJ294413) with a sequence similarity of 97.9% (1476 nt) and strain A214 belonged to Marmoricola aurantiacus DSM 12652T (Y18629) with a sequence similarity of 99.2% (1439 nt).

    Figure 1. Maximum likelihood phylogenetic subtree of 16S rRNA gene sequences from type strains affiliating to the genus Micromonospora as references and the Micromonospora strains derived from Eastern Mediterranean Sea deep-sea sediments. The tree was constructed with 100 bootstrap replicates using phyML software. Despite the low bootstrap values, multiple calculations of the tree always resulted in identical clusters of the sequences. PT = number of phylotypes

    3.3. Potentialof Micromonospora strainsto produce natural products

    Culture extracts of almost half of the Micromonospora strains showed antimicrobial activity against Gram-positive bacteria but none inhibited growth of Gram-negative and eukaryotic test strain. Ten strains showed antimicrobial activity against B. subtilis and most of these in addition exhibited antimicrobial activity against S. lentus (Table 3). In addition, 10 strains without antibiotic activity contained at least one of the tested biosynthetic genes for natural products of the PKS and/or NRPS type.

    Table 3. Antibiotic activity of deep-sea Micromonospora strains and detection of gene clusters for PKS-I, PKS-II and NRPS
    Strain Cluster Antibiotic activity [%] Detection of biosynthetic gene clusters
    B.subtilis S. lentus PKS-I PKS--II NRPSshort NRPSlarge
    S61 A 99 - -
    S20 B 83 -
    S51 B 100 23 - √√
    A77 B - -
    S32a C - - - - - -
    S71 D 96 82 - - - -
    S29 E - - - - -
    A68 E 100 100 - √√ √√
    A75a E - - - √√
    A172 E - - - √√
    A203 E 100 100 -
    A165 E 96 - - √√
    A177 E - - - √√ √√ -
    A200 E 100 100 √√ √√ √√ -
    A208 E - - - √√
    A176 E - - - √√
    A201 E - - √√ √√
    D117x F 42 52 - -
    A74 G 100 69 - - -
    A64 G - - - √√ √√
    A202 G - - √√
    : positive PCR amplification; √√: positive PCR amplification and confirmation by sequence information; -: negative test result.
     | Show Table
    DownLoad: CSV

    In order to substantiate the potential of producing natural products of known biosynthetic pathways, the presence of gene clusters for polyketide synthases of type I and type II (PKS-I, PKS-II) and for non-ribosomal peptide synthetase (NRPS) was analysed in all Micromonosporastrains (Table 3,Figure 2). The amplification of non-ribosomal peptide synthetase genes by two primer sets and PCR protocols resulted in gene fragments of different length (NRPSshort and NRPSlarge). With a few exception all strains contained NRPS genes, most exhibited both gene fragments and a few strains showed amplicons only in one or the other protocol. While genes for PKS-II were present in all but 5 isolates, PKS-I was detected only in a few strains of clusters B, E and G. Sequence similarity to PKS-I and PKS-II genes was in the range of 72% and 91% and that of NRPS genes to reference sequences was >50% and (Table 4).

    Table 4. NRPS, PKS-I and PKS-II gene products related to amino acid sequences derived from the respective PCR amplificates of Micromonospora strains
    Strain Next related protein (protein, producer, accession number) Amino acid match Similarity (%)
    NRPS
    A64 non-ribosomal peptide synthetase (sanglifehrin A), Streptomyces flaveolus DSM9954, ACY06285.1 45/89 50
    A68 non-ribosomal peptide synthetase, Micromonospora sp. ATCC 39149, ZP_04604018.1 45/85 52
    A177 non-ribosomal peptide synthetase, Micromonospora sp. ATCC 39149, ZP_04604018.1 47/88 54
    A200 non-ribosomal peptide synthetase, Micromonospora sp. ATCC 39149, ZP_04604018.1 46/84 54
    A201 non-ribosomal peptide synthetase, Micromonospora sp. L5, ZP_06399932.1 46/81 56
    PKS-I
    A200 type I polyketide synthase , Streptomycesdiastaticus var. 108, AAR87133.1 27/34 79
    PKS-II
    S51 beta-ketoacyl synthase, Micromonospora sp. ATCC 39149, EEP70048.1 130/146 89
    A64 SaqA (saquayamycin Z), Micromonospora sp. Tu 6368, ACP19353.1 62/86 72
    A68 type II ketosynthase, Micromonospora sp. 8g107, EU868817.1 114/142 80
    A75a beta-ketoacyl synthase, Micromonospora sp. L5, EFC60411.1 79/86 91
    A165 SaqA (saquayamycin Z), Micromonospora sp. Tu 6368, ACP19353.1 133/171 77
    A172 type II ketosynthase, Micromonospora sp. 8g107, ACG49888.1 171/200 85
    A176 beta-ketoacyl synthase, Micromonospora sp. ATCC 39149, EEP70048.1 123/135 91
    A177 ketosynthase, Streptomyces albospinus NBRC 13846, BAF43344 43/52 82
    A200 type II ketosynthase, Micromonospora sp. 8g107, ACG49888.1 129/153 84
    A201 beta-ketoacyl synthase, Micromonospora sp. ATCC 39149, EEP70048.1 117/142 82
    A202 beta-ketoacyl synthase, Micromonospora sp. ATCC 39149, EEP70048.1 138/171 80
    A208 type II ketosynthase, Micromonospora sp. 8g107, ACG49888.1 141/172 81
     | Show Table
    DownLoad: CSV
    Figure 2. Presence of PKS-I, PKS-II and NRPS biosynthesis gene clusters in deep-sea Micromonospora isolates

    For those strains, which showed no inhibitory effects on the growth of the test strains, but were positive for PKS and/or NRPS it is expected that the cultivation conditions used in this study were not appropriate to stimulate the biosynthesis of antibiotic compounds.


    3.4. Chemical metabolite profiles of selected strains

    Culture extracts of Micromonospora strains showed large variation in the secondary metabolite profiles as analyzed by HPLC-MS/UV analysis. Similar metabolite profiles were found in representatives of cluster E (strains A177, A200, A203 and A208). While strains A64 and A74 were phylogenetically closely related by affiliating to cluster G, they showed different metabolite profiles. In contrast, the already known compounds levantilide A and its derivative levantilide B were produced by members of different clusters, i.e. the original producer strain A77 [24] and strain S51, which affiliates both to cluster B, but also strain A74 (cluster G).

    The HPLC-MS/UV profile of strain S20 revealed four unknown metabolites, which were not produced by any of the other strains tested. The main compound S20-4 exhibited a molecular weight of 940 (Figure 3). The production of this substance was initiated at the end of the log-phase and lasts throughout the stationary phase, but was independent from the salt concentration. It was found at salinities up to 3.0% of Tropic Marine® salt and also when seawater from the habitat (salinity of 3.9%) was used for cultivation. Three derivatives of this major compound (S20-5, S20-6, and S20-7) showing the molecular weights of 954, 982, and 940, respectively, were observed. Interestingly, these compounds showed variations in the bioactivity towards B. subtilis. While S20-7 exhibited no antibiotic activity, the compounds S20-4, S20-5 and S20-6 inhibited growth of B. subtilis with 99 %, 100 % and 32 % at 100 µM, respectively.

    Figure 3. HPLC-UV profile of Micromonospora strain S20. Top: The DAD spectra of the chromatogram of the crude extract showed several compounds. The strongest UV signal was exhibited by the unknown compound (m/z = 940) which is marked with an arrow. Bottom: The UV spectrum of thiscompound is characterized by two maxima at 298 nm and 310 nm as well as a shoulder at 330nm

    4. Discussion

    Even thoughreports on Actinobacteria isolated from marine habitats are accumulating, so far only few reports deal with Actinobacteria from deep-sea sediments. In the present study, Actinobacteria were isolated from an extreme oligotrophic marine deep-sea habitat, the majority of which was classified as belonging to the genus Micromonospora. This is the first detailed study on the phylogenetic diversity of deep-sea Micromonospora strains from the oligotrophic deep Mediterranean Sea and on their potential to produce bioactive natural products. Diversity, bioactivity and the potential for secondary metabolite biosynthesis of these bacteria were analyzed in this study.


    4.1. Isolationand phylogenetic analyses of the Actinobacteria

    A number of treatments have been applied to improve the isolation success of Actinobacteria. Dry-heat has been found to be useful for the selective enrichment of non-streptomycete actinobacteria [40] and we have used this pre-treatment with dry heat to selectively enrich actinobacteria and to isolate members of Micromonospora, Knoellia and Marmoricola. The selectivity of this pre-treatment is demonstrated by the observation, that members of Streptomycetes were completely missing after the heat-treatment approach applied in this study. This finding demonstrates that in addition to the media and cultivation conditions the sample pre-treatment is important for the isolation of a diverse array of Actinobacteria. The increase in isolation efficiency of these Actinobacteria is of particular interest, because the diversity of non-streptomycete Actinobacteria is hardly covered in eubacterial environmental gene libraries.

    The high diversity and number of phylotypes of the genus Micromonospora (7 clusters, 14 PTs) recovered from the deep-sea sediments in this study is remarkable. In a previous study on these sediments only two phylotypes of Streptomyces (a total of 8 strains) were recovered [28]. This is in accordance with a study at the Pacific coast off San Francisco in whichrecovery of actinomycetes from shallow versus deep sediments was compared [5]. These authors discovered an increasing diversity of Micromonospora relative to Streptomyces with increasing water depth. In the Trondheim Fjord (Norway), a higher abundance of Micromonospora strains was observed among isolates from deeper sediments (450 m) compared to shallow water sediments and Micromonospora isolates were the major group of actinobacteria in the deeper sediment samples [41].

    Noteworthy, most Micromonospora isolates of the present study phylogenetically clustered apart from already described type strains of this genus and distinct clusters could be detected by phylogenetic analysis of the 16S rDNA gene sequences. Clusters E and F (D117x, A165, A68, A208 and A172) contain isolates obtained after pretreatment at in situ pressure and may be considered as (Table 1). As these strains obviously are adapted to changes in hydrostatic pressure (during sampling and incubation procedures) and were selected after pressure incubation, they might be considered to be well adapted to the conditions prevailing in the Mediterranean deep sea sediments. The survival of long lasting treatments of high hydrostatic pressure at low temperature apparently is common to Micromonospora strains and was already noticed earlier [40]. All Micromonospora isolates grew well under characteristic features of the deep Eastern Mediterranean Sea, at sea water from the habitat, at minor nutrient concentrations (i.e. Mediterranean Sea water supplemented with 15 g agar) and at low temperature (13.5 °C). In conclusion, these strains cope well the conditions prevailing in situ, appear well adapted to the conditions of the Mediterranean deep sea and therefore may be regarded as indigenous inhabitants of the investigated sediments.


    4.2. Natural productsfromdeep-sea Micromonospora strains

    As the deep sea generally is characterized by oligotrophy with spontaneous and punctual nutrient input, inhabitants must be adapted to these features and compete with others when nutrients become available. Therefore, the use of antimicrobial compounds in order to outcompete competitors when nutrient become available appears to be one reasonable strategy and might explain why secondary metabolite production is an important feature of deep-sea bacteria.

    Antimicrobial activity of the Micromonospora strains was specifically targeted towards Gram-positive bacteria and approx. half of the tested strains inhibited growth of B. subtilis and half of these strains inhibited in addition S. lentus. If one considers, that Gram-positive bacteria and in particular members of the genus Bacillus are important players in the marine sediments, they can be expected being primary competitors of Actinobacteria. According to a former study, isolates from the Mediterranean deep-sea sediments were dominated by a highly diverse group of Gram-positive bacteria, in particular by members of the genus Bacillus, all of which cope quite well with minor nutrient concentrations [28]. Since the Bacillus strains and the Micromonospora strains have similar nutrient requirements, the production of antibiotics might be a useful tool of Micromonospora to compete against Bacillus in marine environments. Inhibition of Bacillus might be of advantage for the occurrence, persistence and growth of the antibiotic producing Micromonospora community. The importance of bioactive secondary metabolites for Micromonospora is supported by the detection of multiple biosynthesis genes for secondary metabolite production (PKS-I, PKS-II and NRPS) in the great majority of these isolates (86%). As far as tested, the Micromonospora strains provide a rich source of compounds, which will be identified in further studies.


    5. Conclusion

    Growth characteristics of Micromonospora strains isolated from Mediterranean deep sea sediments and their affiliation to phylogenetic clusters apart from known type strains indicate their potential to adapt to the conditions prevailing in situ. At least some of them can cope with pressure conditions prevailing in these sediments. Therefore, they may be regarded as indigenous bacteria of these deep-sea sediments. They have a broad genetic potential to produce secondary metabolites, which quite likely is significantly exceeding those detected in the culture extracts. As we know, that cultivation conditions have strong impact on the production of secondary metabolites [42], the choice of media and cultivation conditions are of major importance for the production of secondary metabolites. Further experiments are expected to significantly enlarge the number of metabolites produced by the Micromonospora strains isolated from deep sea sediments.


    Acknowledgments

    The authors are grateful to Arlette Wenzel-Storjohann for taking care about the antibiotic assays, to Dr. Markus Schilhabel and his team at the Institute for Clinical Molecular Biology (University Hospital Schleswig-Holstein, Kiel, Germany) for providing Sanger sequencing and to the DFG for support through the Clusters of Excellence “Inflammation at interfaces” and “Future Ocean”.


    Conflict of Interest

    Authors declare no conflict of interest in this paper.




    [1] World Health Organization (WHO), 2018. Available from: https://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue.
    [2] D. J. Gubler, E. E. Ooi, S. Vasudevan, J. Farrar, Dengue and dengue hemorrhagic fever, CABI (2014).
    [3] M. G. Guzman, M. Alvarez, S. B. Halstead, Secondary infection as a risk factor for dengue hemorrhagic fever/dengue shock syndrome: an historical perspective and role of antibody-dependent enhancement of infection, {Arch. Virol.}, 158 (2013), 1445–1459. https://doi.org/10.1007/s00705-013-1645-3. doi: 10.1007/s00705-013-1645-3
    [4] N. G. Reich, S. Shrestha, A. A. King, P. Rohani, J. Lessler, S. Kalayanarooj, et al., Interactions between serotypes of dengue highlight epidemiological impact of cross-immunity, {J. R. Soc. Interface}, 10 (2013), 4–14. https://doi.org/10.1098/rsif.2013.0414. doi: 10.1098/rsif.2013.0414
    [5] B. Adams, M. Boots, Modelling the relationship between antibody-dependent enhancement and immunological distance with application to dengue, {J. Theor. Biol.}, 242 (2006), 337–346. https://doi.org/10.1016/j.jtbi.2006.03.002. doi: 10.1016/j.jtbi.2006.03.002
    [6] S. Bianco, L. B. Shaw, I. B. Schwartz, Epidemics with multistrain interactions: the interplay between cross immunity and antibody-dependent enhancement, {Chaos}, 19 (2009), 043123. https://doi.org/10.1063/1.3270261. doi: 10.1063/1.3270261
    [7] L. Billings, I. B. Schwartz, L. B. Shaw, M. McCrary, D. S. Burke, D. A. T. Cummings, Instabilities in multiserotype disease models with antibody-dependent enhancement, {J. Theor. Biol.}, 246 (2007), 18–27. https://doi.org/10.1016/j.jtbi.2006.12.023. doi: 10.1016/j.jtbi.2006.12.023
    [8] K. Hu, C. Thoens, S. Bianco, S. Edlund, M. Davis, J. Douglas, et al., The effect of antibody-dependent enhancement, cross immunity, and vector population on the dynamics of dengue fever, { J. Theor. Biol.}, 319 (2013), 62–74. https://doi.org/10.1016/j.jtbi.2012.11.021. doi: 10.1016/j.jtbi.2012.11.021
    [9] M. Aguiar, N. Stollenwerk, A new chaotic attractor in a basic multi-strain epidemiological model with temporary cross-immunity, preprint (2007), arXiv: 0704.3174v1.
    [10] M. Aguiar, B. Kooi, N. Stollenwerk, Epidemiology of dengue fever: A model with temporary cross immunity and possibly secondary infection shows bifurcations and chaotic behaviors in wide parameter region, Math. Model. Nat. Phenom., 3 (2008), 48–70. https://doi.org/10.1051/mmnp:2008070. doi: 10.1051/mmnp:2008070
    [11] B. W. Kooi, M. Aguiar, N. Stollenwerk, Analysis of an asymmetric two-strain dengue model, Math. Biosci., 248 (2014), 128–139. https://doi.org/10.1016/j.mbs.2013.12.009. doi: 10.1016/j.mbs.2013.12.009
    [12] P. Van den Driessche, Some epidemiological Models with delays. In: Differential Equations and Application to biology and to Industry, World Sci. (1996), 507–520.
    [13] K. Nah, Y. Nakata, G. Rost, Malaria dynamics with long incubation period in host, Comp. Math. with App., 68 (2014), 915–930. https://doi.org/10.1016/j.camwa.2014.05.001. doi: 10.1016/j.camwa.2014.05.001
    [14] D. Chen, Z. Xu, Global dynamics of a delayed diffusive two-strain disease model, Diff. Eq. App., 8 (2016), 99–122. https://doi.org/10.7153/dea-08-07. doi: 10.7153/dea-08-07
    [15] J. Guan, L. Wu., M. Chen, X. Dong, H. Tang, Z. Chen, The stability and Hopf Bifurcation of the dengue fever model with time delay, It. J. Pure App. Math., 37 (1973), 139–156.
    [16] K. Hattaf, A. A. Lashari, Y. Louartassi, N. Yousfi, A delayed SIR epidemic model with a general incidence rate, Elect. J. Qual. The. Differ. Equ., 3 (2013), 1–9. https://doi.org/10.14232/ejqtde.2013.1.3. doi: 10.14232/ejqtde.2013.1.3
    [17] C. Huang, J. Cao, F. Wen, X. Yang, X., Stability Analysis of SIR Model with Distributed Delay on Complex Networks, PLoS ONE, 11 (2016), e0158813. https://doi.org/10.1371/journal.pone.0158813. doi: 10.1371/journal.pone.0158813
    [18] J. Xu, Y. Geng, Y. Zhou, Global stability of a multi-group model with distributed delay and vaccination, Math. Meth. App. Sci., 40 (2017), 1475–1486. https://doi.org/10.1002/mma.4068. doi: 10.1002/mma.4068
    [19] L. C. Katzelnick, L. Gresh, M. E. Halloran, J. C. Mercado, G. Kuan, A. Gordon, et al., Antibody-dependent enhancement of severe dengue disease in humans, { Science}, 358 (2017), 929–932. https://doi.org/10.1126/science.aan6836. doi: 10.1126/science.aan6836
    [20] A. Rothman, Immunity to dengue virus: a tale of original antigenic sin and tropical cytokine storms, Nat. Rev. Immunol., 11 (2011), 532–543. https://doi.org/10.1038/nri3014. doi: 10.1038/nri3014
    [21] L. Wang, Y. Li, L. Pang, Dynamics Analysis of an Epidemiological Model with media impact and two delays, Math. Prob. Eng., (2016), Article ID 1598932. https://doi.org/10.1155/2016/1598932.
    [22] K. L. Cooke, P. Van den Driessche, Analysis of an SEIRS epidemic model with two delays, J. Math. Biol., 35 (1996), 240–260.
    [23] P. Van den Driessche, L. Wang, X. Zou, Modeling diseases with latency and relapse, Math. Biosci. Eng., 4 (2007), 205–219. https://doi.org/10.3934/mbe.2007.4.205. doi: 10.3934/mbe.2007.4.205
    [24] P. Van den Driessche, J. Watmough, Further Notes on the Basic Reproduction Number, Math. Ep., (2008), 159–178.
    [25] M. Martcheva, Introduction to Mathematical Epidemiology, Springer, New York (2015).
    [26] C. S. VinodKumar, N. K. Kalapannavar, K. G. Basavarajappa, D. Sanjay, C. Gowli, N. G. Nadig, et al., Episode of coexisting infections with multiple dengue virus serotypes in central Karnataka, J. Inf. Public. Health, 6 (2013), 302–306. https://doi.org/10.1016/j.jiph.2013.01.004. doi: 10.1016/j.jiph.2013.01.004
    [27] F. Brauer, Asymptotic stability of a class of integro-differential equation, J. Diff. Eq., 28 (1978), 180–188.
    [28] R. K. Miller, Nonlinear Volterra Integral Equations, W. A. Benjamin, California, (1971).
    [29] H. K. Hethcote, H. W. Stech, P. Van den Driessche, Nonlinear Oscillations in Epidemic Models, J. Appl. Math., 40 (1981).
    [30] R. K. Miller, Asymptotic stability properties of linear Volterra Integro-differential equations, J. Diff. Eq., 10 (1971), 485–506.
    [31] X. Feng, K. Wang, F. Zhang, Z. Teng, Threshold dynamics of a nonlinear multigroup epidemic model with two infinite distributed delays, Math. Meth. Appl. Sci., 40 (2016), 2762–2771. https://doi.org/10.1002/mma.4196. doi: 10.1002/mma.4196
    [32] J. Wang, Y. Takeuchi, A multi-group SVEIR epidemic model with distributed delay and vaccination, Int. J. Bio., 5 (2012), 1260001. https://doi.org/10.1142/S1793524512600017. doi: 10.1142/S1793524512600017
    [33] M. Y. Li, Z. Shuai, C. Wang, Global stability of multi-group epidemic models with distributed delay, J. Math. Anal. Appl., 361 (2010), 38–47. https://doi.org/10.1016/j.nonrwa.2011.11.016. doi: 10.1016/j.nonrwa.2011.11.016
    [34] G. Röst, J. Wu, SEIR epidemiological model with varying infectivity and infinite delay, Math. Bio. Eng, 5 (2008), 389–402. https://doi.org/10.3934/mbe.2008.5.389. doi: 10.3934/mbe.2008.5.389
    [35] S. I. Grossman, R. K. Miller, Nonlinear Volterra Integro-differential systems with L1-Kernels, J. Diff. Eq., 13 (1973), 551–556.
    [36] T. A. Burton, Volterra Integral and Differential Equations, in Mathematics in Science and Engineerin, Elsevier, (2005), 1–355.
    [37] J. P. LaSalle, The Stability of Dynamical Systems, In Regional Conference Series in Applied Mathematics, SIAM, USA (1976).
    [38] IBGE/DPE/Coordenação de População e Indicadores Sociais, Gerência de Estudos e Análises da Dinâmica Demográfica, 2018. Avaiable from: http://tabnet.datasus.gov.br/cgi/idb2012/a11tb.htm.
    [39] World Health Organization (2009) Dengue: Guidelines for Diagnosis, Treatment, Prevention and Control, New edition, A joint publication of the World Health Organization (WHO) and the Special Programme for Research and Training in Tropical Diseases (TDR), 2009. Available from: https://www.who.int/tdr/publications/documents/dengue-diagnosis.pdf.
    [40] N. Ferguson, R. Anderson, S. Gupta, The effect of antibody-dependent enhancement on the transmission dynamics and persistence of multiple strain pathogens, Proc. Natl. Acad. Sci., 96 (1999), 790–794.
    [41] S. B. Maier, X. Huanga, E. Massad, M. Amaku, M. N. Burattini, D. Greenhalgh, Analysis of the optimal vaccination age for dengue in Brazil with a tetravalent dengue vaccine, Math. Biosci., 294 (2017), 15—32. https://doi.org/10.1016/j.mbs.2017.09.004. doi: 10.1016/j.mbs.2017.09.004
    [42] E. Massad, F. A. Coutinho, M. N. Burattini, L. F. Lopez, The risk of yellow fever in a dengue-infested area, {Trans. R. Soc. Trop. Med. Hyg.}, 95 (2001), 370–374. https://doi.org/10.1016/s0035-9203(01)90184-1. doi: 10.1016/s0035-9203(01)90184-1
    [43] R. C. Reiner, S. T. Stoddard, B. M. Forshey, A. A. King, A. M. Ellis, A. L. Lloyd, et al., Time-varying, serotype-specific force of infection of dengue virus, Proc. Natl. Acad. Sci., 111 (2014), E2694–E2702. https://doi.org/10.1073/pnas.1314933111. doi: 10.1073/pnas.1314933111
    [44] N. M. Ferguson, C. A. Donnelly, R. M. Anderson, Transmission dynamics and epidemiology of dengue: Insights from age-stratified sero-prevalence surveys, { Phil. Trans. R. Soc. B Biol. Sci.}, 354 (1999), 757—768.
    [45] L. Edelstein-Keshet, Mathematical Models in Biology, SIAM PA, USA (2005).
    [46] J. D. Murray, Mathematical Biology I. An Introduction, 3rd ed, Springer (2002).
    [47] S. Lynch, Dynamical systems with applications using MATLAB, Springer, (2004).
    [48] G. Huang, A. Liu, A note on global stability for a heroin epidemic model with distributed delay, App. Math. Let., 26 (2013), 687–691. https://doi.org/10.1016/j.aml.2013.01.010. doi: 10.1016/j.aml.2013.01.010
    [49] J. Hale, S. M. V. Lunel, Introduction to Functional Differential Equations, Vol. 9, Applied Mathematical Science, New York, (1993).
  • This article has been cited by:

    1. Yi-Lei Nie, Yun-Dan Wu, Chuan-Xi Wang, Ru Lin, Yang Xie, Dong-Sheng Fang, Hong Jiang, Yun-Yang Lian, Structure elucidation and antitumour activity of a new macrolactam produced by marine-derived actinomycete Micromonospora sp. FIM05328, 2018, 32, 1478-6419, 2133, 10.1080/14786419.2017.1366479
    2. Jutta Wiese, Johannes F. Imhoff, Marine bacteria and fungi as promising source for new antibiotics, 2019, 80, 02724391, 24, 10.1002/ddr.21482
    3. Ming-Hong Chen, Yun-Yang Lian, Dong-Sheng Fang, Li Chen, Jia Jia, Wen-Long Zhang, Ru Lin, Yang Xie, Hong-Kai Bi, Hong Jiang, Identification and antimicrobial properties of a new alkaloid produced by marine-derived Verrucosispora sp. FIM06-0036, 2019, 1478-6419, 1, 10.1080/14786419.2019.1689498
    4. Agnieszka A. Pilarska, Agnieszka Wolna-Maruwka, Alicja Niewiadomska, Jarosław Grządziel, Anna Gałązka, Emil Paluch, Klaudia Borowiak, Krzysztof Pilarski, Quantitative and Qualitative Changes in the Genetic Diversity of Bacterial Communities in Anaerobic Bioreactors with the Diatomaceous Earth/Peat Cell Carrier, 2022, 11, 2073-4409, 2571, 10.3390/cells11162571
    5. Agnieszka A. Pilarska, Agnieszka Wolna-Maruwka, Alicja Niewiadomska, Krzysztof Pilarski, Mariusz Adamski, Aleksandra Grzyb, Jarosław Grządziel, Anna Gałązka, Silica/Lignin Carrier as a Factor Increasing the Process Performance and Genetic Diversity of Microbial Communities in Laboratory-Scale Anaerobic Digesters, 2021, 14, 1996-1073, 4429, 10.3390/en14154429
    6. Uğur ÇİĞDEM, Ayten KUMAŞ, Fadime ÖZDEMİR KOÇAK, Sakarya Nehir Sedimentinden İzole Edilen Aktinobakterilerin Antimikrobiyal ve Bitki Gelişim Teşvik Edici Özelliklerinin Belirlenmesi, 2021, 11, 1309-4726, 239, 10.31466/kfbd.889423
    7. Pankaj Verma, Vikas Pandey, Seyieleno C. Seleyi, Abirami Alagarsamy, Gopal Dharani, Exploring the hidden treasures: Deep-sea bacterial community structure in the Bay of Bengal and their metabolic profile, 2024, 10, 2296-7745, 10.3389/fmars.2023.1308953
  • Reader Comments
  • © 2022 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(3047) PDF downloads(134) Cited by(6)

Article outline

Figures and Tables

Figures(12)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog