An amendment to this paper has been published and can be accessed via the original article.
Background Natural disasters and infectious diseases result in widespread disruption to human health and livelihood. At the scale of a global pandemic, the co-occurrence of natural disasters is inevitable. However, the impact of natural disasters on the spread of COVID-19 has not been extensively evaluated through epidemiological modelling. Methods We create an agent-based epidemiology model based on COVID-19 clinical, epidemiological, and geographic data. We first model 35 scenarios with varying natural disaster timing and duration for a COVID-19 outbreak in a theoretical region. We then evaluate the potential effect of an eruption of Vesuvius volcano on the spread of COVID-19 in Campania, Italy. Results In a majority of cases, the occurrence of a natural disaster increases the number of disease related fatalities. For a natural disaster fifty days after infection onset, the median increase in fatalities is 2, 59, and 180% for a 2, 14, and 31-day long natural disaster respectively, when compared to the no natural disaster scenario. For the Campania case, the median increase in fatalities is 1.1 and 2.4 additional fatalities per 100,000 for eruptions on day 1 and 100 respectively, and 60.0 additional fatalities per 100,000 for an eruption close to the peak in infections (day 50). Conclusion Our results show that the occurrence of a natural disaster in most cases leads to an increase in infection related fatalities, with wide variance in possible outcomes depending on the timing of the natural disaster relative to the peak in infections and the duration of the natural disaster.
BACKGROUND:The COVID-19 pandemic poses a serious threat to global health, and pathogenic mutations are a major challenge to disease control. We developed a statistical framework to explore the association between molecular-level mutation activity of SARS-CoV-2 and population-level disease transmissibility of COVID-19.METHODS:We estimated the instantaneous transmissibility of COVID-19 by using the time-varying reproduction number (Rt). The mutation activity in SARS-CoV-2 is quantified empirically depending on (i) the prevalence of emerged amino acid substitutions and (ii) the frequency of these substitutions in the whole sequence. Using the likelihood-based approach, a statistical framework is developed to examine the association between mutation activity and Rt. We adopted the COVID-19 surveillance data in California as an example for demonstration.RESULTS:We found a significant positive association between population-level COVID-19 transmissibility and the D614G substitution on the SARS-CoV-2 spike protein. We estimate that a per 0.01 increase in the prevalence of glycine (G) on codon 614 is positively associated with a 0.49% (95% CI: 0.39 to 0.59) increase in Rt, which explains 61% of the Rt variation after accounting for the control measures. We remark that the modeling framework can be extended to study other infectious pathogens.CONCLUSIONS:Our findings show a link between the molecular-level mutation activity of SARS-CoV-2 and population-level transmission of COVID-19 to provide further evidence for a positive association between the D614G substitution and Rt. Future studies exploring the mechanism between SARS-CoV-2 mutations and COVID-19 infectivity are warranted.
Successful clinical drug development requires rational design of combination treatments based on preclinical data. Anti-hepatitis C virus (HCV) drugs exhibit significant diversity in antiviral effect. Dose-response assessments can be used to determine parameters profiling the diverse antiviral effect during combination treatment. In the current study, a combined experimental and mathematical approaches were used to compare and score different combinations of anti-HCV treatments. A “required concentration index” was generated and used to rank the antiviral profile of possible double- and triple-drug combinations against HCV genotype 1b and 2a. Rankings varied based on target HCV genotype. Interestingly, multidrug (double and triple) treatment not only augmented antiviral activity, but also reduced genotype-specific efficacy, suggesting another advantage of multidrug treatment. The current study provides a quantitative method for profiling drug combinations against viral genotypes, to better inform clinical drug development.
BACKGROUND:A hepatitis A epidemic occurred among men who have sex with men (MSM) in Japan in 2017-2018. In this study, we employ a parsimonious mathematical model to epidemiologically investigate the dynamics of infection, aiming to evaluate the effectiveness of campaign-based interventions among MSM to raise awareness of the situation.METHODS:A mathematical model describing a mixture of human-to-human transmission and environmental transmission was fitted to surveillance data. Taking seasonally varying environmental transmission into account, we estimated the reproduction number of hepatitis A virus during the course of epidemic, and, especially, the abrupt decline in this reproduction number following campaign-based interventions.RESULTS:The reproduction number prior to the countermeasures ranged from 2.6 to 3.1 and then began to decrease following campaign-based interventions. After the first countermeasure, the reproduction number decreased, but the epidemic remained supercritical (i.e., Rt > 1). The value of Rt dropped well below one following the second countermeasure, which used web articles to widely disseminate information about the epidemic risk.CONCLUSIONS:Although the effective reproduction number, Rt, changes because of both intrinsic and extrinsic factors, the timing of the examined countermeasures against hepatitis A in the MSM population was consistent with the abrupt declines observed in Rt. Even without vaccination, the epidemic was brought under control, and risky behaviors may have been changed by the increase in situation awareness reached through web articles.
Sexually transmitted infections (STIs) have remained a worldwide public health threat. It is difficult to control the spread of STIs, not only because of heterogeneous sexual transmission between men and women but also because of the complicated effects of sexual transmission among men who have sex with men (MSM) and mother-to-child transmission. Many studies point to the existence of a 'bisexual bridge', where STIs spread from the MSM network via bisexual connections. However, it is unclear how the MSM network affects heterosexual networks as well as mother-to-child transmission. To analyse the effect of MSM on the spread of STIs, we divided the population into four subpopulations: (i) women, (ii) men who have sex with women only (MSW), (iii) men who have sex with both men and women (MSMW), (iv) men who have sex with men exclusively (MSME). We calculated the type-reproduction numbers of these four subpopulations, and our analysis determined what preventive measures may be effective. Our analysis shows the impact of bisexual bridge on the spread of STIs does not outweigh their population size. Since MSM and mother-to-child transmission rates do not have a strong synergistic effect when combined, complementary prevention measures are needed. The methodologies and findings we have provided here will contribute greatly to the future development of public health.
Travelers play a role in triggering epidemics of imported dengue fever because they can carry the virus to other countries during the incubation period. If a traveler carrying dengue virus visits open green space and is bitten by mosquitoes, a local outbreak can ensue. In the present study, we aimed to understand the movement patterns of international travelers in Tokyo using mobile phone data, with the goal of identifying geographical foci of dengue transmission. We analyzed datasets based on mobile phone access to WiFi systems and measured the spatial distribution of international visitors in Tokyo on two specific dates (one weekday in July 2017 and another weekday in August 2017). Mobile phone users were classified by nationality into three groups according to risk of dengue transmission. Sixteen national parks were selected based on their involvement in a 2014 dengue outbreak and abundance of Aedes mosquitoes. We found that not all national parks were visited by international travelers and that visits to cemeteries were very infrequent. We also found that travelers from countries with high dengue prevalence were less likely to visit national parks compared with travelers from dengue-free countries. Travelers from countries with sporadic dengue cases and countries with regional transmission tended to visit common destinations. By contrast, the travel footprints of visitors from countries with continuous dengue transmission were focused on non-green spaces. Entomological surveillance in Tokyo has been restricted to national parks since the 2014 dengue outbreak. However, our results indicate that areas subject to surveillance should include both public and private green spaces near tourist sites.
Background Nondepolarizing neuromuscular blocking drugs (NDNBs) are clinically used to produce muscle relaxation during general anesthesia. To better understand clinical properties of NDNBs, comparative in vitro pharmacologic studies have been performed. In these studies, a receptor binding model, which relies on the assumption that the inhibition, i.e., the effect of an NDNB, is proportional to the receptor occupancy by the drug, has been effectively used to describe obtained experimental data. However, it has not been studied in literature under which conditions the above assumption can be justified nor the assumption still holds in vivo. The purpose of this study is to explore the in vivo relationship between the inhibition and the receptor occupancy by an NDNB and to draw implications on how in vitro experimental results can be used to discuss the in vivo properties of NDNBs. Methods An ordinary differential equation model is employed to simulate physiologic processes of the activation of receptors by acetylcholine (ACh) as well as inhibition by an NDNB. With this model, the degree of inhibition is quantified by the fractional amount of receptors that are not activated by ACh due to the presence of an NDNB. The results are visualized by plotting the fractional amounts of the activated receptors as a function of the receptor occupancy. Results Numerical investigations reflecting in vivo conditions show that the degree of inhibition is not proportional to the receptor occupancy, i.e., there is a nonlinear relationship between the inhibition and the receptor occupancy. However, under a setting of high concentration of ACh reflecting a typical situation of in vitro experiments, the relationship between the inhibition and the receptor occupancy becomes linear, suggesting the validity of the receptor binding model. Also, it is found that the extent of nonlinearity depends on the selectivity of NDNBs for the two binding sites of the receptors. Conclusions While the receptor binding model may be effective for estimating affinity of an NDNB through in vitro experiments, these models do not directly describe in vivo properties of NDNBs, because the nonlinearity between the inhibition and the receptor occupancy causes the modulation of the resultant concentration-effect relationships of NDNBs.
Background This study aimed to jointly model HIV disease progression patterns based on viral load (VL) among adult ART patients adjusting for the time-varying “incremental transients states” variable, and the CD4 cell counts orthogonal variable in a single 5-stage time-homogenous multistate Markov model. We further jointly mapped the relative risks of HIV disease progression outcomes (detectable VL (VL ≥ 50copies/uL) and immune deterioration (CD4 < 350cells/uL) at the last observed visit) conditional not to have died or become loss to follow-up (LTFU). Methods Secondary data analysis of individual-level patients on ART was performed. Adjusted transition intensities, hazard ratios (HR) and regression coefficients were estimated from the joint multistate model of VL and CD4 cell counts. The mortality and LTFU transition rates defined the extent of patients’ retention in care. Joint mapping of HIV disease progression outcomes after ART initiation was done using the Bayesian intrinsic Multivariate Conditional Autoregressive prior model. Results The viral rebound from the undetectable state was 1.78times more likely compared to viral suppression among patients with VL ranging from 50-1000copies/uL. Patients with CD4 cell counts lower than expected had a higher risk of viral increase above 1000copies/uL and death if their VL was above 1000copies/uL (state 2 to 3 ( λ 23 ): HR = 1.83 and ( λ 34 ): HR = 1.42 respectively). Regarding the time-varying effects of CD4 cell counts on the VL transition rates, as the VL increased, ( λ 12 and λ 23 ) the transition rates increased with a decrease in the CD4 cell counts over time. Regardless of the individual’s VL, the transition rates to become LTFU decreased with a decrease in CD4 cell counts. We observed a strong shared geographical pattern of 66% spatial correlation between the relative risks of detectable VL and immune deterioration after ART initiation, mainly in Matabeleland North. Conclusion With high rates of viral rebound, interventions which encourage ART adherence and continual educational support on the barriers to ART uptake are crucial to achieve and sustain viral suppression to undetectable levels. Area-specific interventions which focus on early ART screening through self-testing, behavioural change campaigns and social support strategies should be strengthened in heavily burdened regions to sustain the undetectable VL. Sustaining undetectable VL lowers HIV transmission in the general population and this is a step towards achieving zero HIV incidences by 2030.
Background At the moment we have more than 177 million cases and 3.8 million deaths (as of June 2021) around the world and vaccination represents the only hope to control the pandemic. Imperfections in planning vaccine acquisition and difficulties in implementing distribution among the population, however, have hampered the control of the virus so far. Methods We propose a new mathematical model to estimate the impact of vaccination delay against the 2019 coronavirus disease (COVID-19) on the number of cases and deaths due to the disease in Brazil. We apply the model to Brazil as a whole and to the State of Sao Paulo, the most affected by COVID-19 in Brazil. We simulated the model for the populations of the State of Sao Paulo and Brazil as a whole, varying the scenarios related to vaccine efficacy and compliance from the populations. Results The model projects that, in the absence of vaccination, almost 170 thousand deaths and more than 350 thousand deaths will occur by the end of 2021 for Sao Paulo and Brazil, respectively. If in contrast, Sao Paulo and Brazil had enough vaccine supply and so started a vaccination campaign in January with the maximum vaccination rate, compliance and efficacy, they could have averted more than 112 thousand deaths and 127 thousand deaths, respectively. In addition, for each month of delay the number of deaths increases monotonically in a logarithmic fashion, for both the State of Sao Paulo and Brazil as a whole. Conclusions Our model shows that the current delay in the vaccination schedules that is observed in many countries has serious consequences in terms of mortality by the disease and should serve as an alert to health authorities to speed the process up such that the highest number of people to be immunized is reached in the shortest period of time.
Background In Japan, a part of confirmed patients’ samples have been screened for the variant of concern (VOC), including the variant alpha with N501Y mutation. The present study aimed to estimate the actual number of cases with variant alpha and reconstruct the epidemiological dynamics. Methods The number of cases with variant alpha out of all PCR confirmed cases was estimated, employing a hypergeometric distribution. An exponential growth model was fitted to the growth data of variant alpha cases over fourteen weeks in Tokyo. Results The weekly incidence with variant alpha from 18–24 January 2021 was estimated at 4.2 (95% confidence interval (CI): 0.7, 44.0) cases. The expected incidence in early May ranged from 420–1120 cases per week, and the reproduction number of variant alpha was on the order of 1.5 even under the restriction of contact from January-March, 2021, Tokyo. Conclusions The variant alpha was predicted to swiftly dominate COVID-19 cases in Tokyo, and this has actually occurred by May 2021. Devising the proposed method, any country or location can interpret the virological sampling data.
Background Individuals with asymptomatic severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection can propagate the virus unknowingly and thus have been a focus of public health attentions since the early stages of the pandemic. Understanding viral transmissibility among asymptomatic individuals is critical for successful control of coronavirus disease 2019 (COVID-19). The present study aimed to understand SARS-CoV-2 transmissibility among young asymptomatic individuals and to assess whether symptomatology was associated with transmission of symptomatic vs. asymptomatic infections. Methods We analyzed one of the first-identified clusters of SARS-CoV-2 infections with multiple chains of transmission that occurred among university students in March 2020 in Kyoto prefecture, Japan, using discrete and two-type branching process models. Assuming that the number of secondary cases resulting from either primary symptomatic or asymptomatic cases independently followed negative binomial distributions, we estimated the relative reproduction numbers of an asymptomatic case compared with a symptomatic case. To explore the potential association between symptomatology and transmission of symptomatic vs. asymptomatic incident infections, we also estimated the proportion of secondary symptomatic cases produced by primary symptomatic and asymptomatic cases. Results The reproduction number for a symptomatic primary case was estimated at 1.14 (95% confidence interval [CI]: 0.61–2.09). The relative reproduction number for asymptomatic cases was estimated at 0.19 (95% CI: 0.03–0.66), indicating that asymptomatic primary cases did not result in sufficient numbers of secondary infections to maintain chains of transmission. There was no apparent tendency for symptomatic primary cases to preferentially produce symptomatic secondary cases. Conclusions Using data from a transmission network during the early epidemic in Japan, we successfully estimated the relative transmissibility of asymptomatic cases of SARS-CoV-2 infection at 0.22. These results suggest that contract tracing focusing on symptomatic index cases may be justified given limited testing capacity.
Background : Cancer is one of the leading death causes globally with about 8.2 million deaths per year and an increase in numbers in recent years. About 90% of cancer deaths do not occur due to primary tumors but due to metastases, of which most are not clinically identifiable because of their relatively small size at primary diagnosis and limited technical possibilities. However, therapeutic decisions are formed depending on the existence of metastases and their properties. Therefore non-identified metastases might have huge influence in the treatment outcome. The quantification of clinically visible and invisible metastases is important for the choice of an optimal treatment of the individual patient as it could clarify the burden of non identifiable tumors as well as the future behavior of the cancerous disease. Results : The mathematical model presented in this study gives insights in how this could be achieved, taking into account different treatment possibilities and therefore being able to compare therapy schedules for individual patients with different clinical parameters. The framework was tested on three patients with non-small cell lung cancer, one of the deadliest types of cancer worldwide, and clinical history including platinum-based chemotherapy and PD-L1-targeted immunotherapy. Results yield promising insights into the framework to establish methods to quantify effects of different therapy methods and prognostic features for individual patients already at stage of primary diagnosis.
BACKGROUND:The superchiasmatic nucleus (SCN) serves as the primary circadian (24hr) clock in mammals and is known to control important physiological functions such as the sleep-wake cycle, hormonal rhythms, and neurotransmitter regulation. Experimental results suggest that some of these functions reciprocally influence circadian rhythms, creating a highly complex network. Among the clock's downstream products, orphan nuclear receptors REV-ERB and ROR are particularly interesting because they coordinately modulate the core clock circuitry. Recent experimental evidence shows that REV-ERB and ROR are not only crucial for lipid metabolism but are also involved in dopamine (DA) synthesis and degradation, which could have meaningful clinical implications for conditions such as Parkinson's disease and mood disorders.METHODS:We create a mathematical model consisting of differential equations that express how the circadian variables are influenced by light, how REV-ERB and ROR feedback to the clock, and how REV-ERB, ROR, and BMAL1-CLOCK affect the dopaminergic system. The structure of the model is based on the findings of experimentalists.RESULTS:We compare our model predictions to experimental data on clock components in different light-dark conditions and in the presence of genetic perturbations. Our model results are consistent with experimental results on REV-ERB and ROR and allow us to predict the circadian variations in tyrosine hydroxylase and monoamine oxidase seen in experiments. By connecting our model to an extant model of dopamine synthesis, release, and reuptake, we are able to predict circadian oscillations in extracellular DA and homovanillic acid that correspond well with experimental observations.CONCLUSIONS:The predictions of the mathematical model are consistent with a wide variety of experimental observations. Our calculations show that the mechanisms proposed by experimentalists by which REV-ERB, ROR, and BMAL1-CLOCK influence the DA system are sufficient to explain the circadian oscillations observed in dopaminergic variables. Our mathematical model can be used for further investigations of the effects of the mammalian circadian clock on the dopaminergic system. The model can also be used to predict how perturbations in the circadian clock disrupt the dopaminergic system and could potentially be used to find drug targets that ameliorate these disruptions.
Background Influenza is a public health issue that needs to be addressed strategically. The assessment of detailed infectious profiles is an important part of this effort. Household transmission data play a key role in estimating such profiles. We used diagnostic and questionnaire-based data on influenza patients at a Japanese clinic to estimate the detailed infectious period (as well as incubation period, symptomatic and infectious periods, and extended infectious period after recovery) and the secondary attack ratio (SAR) of influenza for households of various sizes based on a modified Cauchemez-type model. Results The data were from enrolled patients with confirmed influenza who were treated at the Hirotsu Clinic (Kawasaki, Japan) with a neuraminidase inhibitor (NAI) during six northern hemisphere influenza seasons between 2010 and 2016. A total of 2,342 outpatients, representing 1,807 households, were included. For influenza type A, the average incubation period was 1.43 days (95% probability interval: 0.03-5.32 days). The estimated average symptomatic and infective period was 1.76 days (0.33-4.62 days); the extended infective period after recovery was 0.25 days. The estimated SAR rose from 20% to 32% as household size increased from 3 to 5. For influenza type B, the average incubation period, average symptomatic and infective period, and extended infective period were estimated as 1.66 days (0.21-4.61), 2.62 days (0.54-5.75) and 1.00 days, respectively. The SAR increased from 12% to 21% as household size increased from 3 to 5. Conclusion All estimated periods of influenza type B were longer than the corresponding periods for type A. However, the SAR for type B was less than that for type A. These results may reflect Japanese demographics and treatment policy. Understanding the infectious profiles of influenza is necessary for assessing public health measures.
Background The effects of arterial wall compliance on blood flow have been revealed using fluid-structure interaction in last decades. However, microcirculation is not considered in previous researches. In fact, microcirculation plays a key role in regulating blood flow. Therefore, it is very necessary to involve microcirculation in arterial hemodynamics. Objective The main purpose of the present study is to investigate how wall compliance affects the flow characteristics and to establish the comparisons of these flow variables with rigid wall when microcirculation is considered. Methods We present numerical modeling in arterial hemodynamics incorporating fluid-structure interaction and microcirculation. A novel outlet boundary condition is employed to prescribe microcirculation in an idealised model. Results The novel finding in this work is that wall compliance under the consideration of microcirculation leads to the increase of wall shear stress in contrast to rigid wall, contrary to the traditional result that wall compliance makes wall shear stress decrease when a constant or time dependent pressure is specified at an outlet. Conclusions This work provides the valuable study of hemodynamics under physiological and realistic boundary conditions and proves that wall compliance may have a positive impact on wall shear stress based on this model. This methodology in this paper could be used in real model simulations.
Background Understanding the dynamical behavior of dengue transmission is essential in designing control strategies. Mathematical models have become an important tool in describing the dynamics of a vector borne disease. Classical compartmental models are well–known method used to identify the dynamical behavior of spread of a vector borne disease. Due to use of fixed model parameters, the results of classical compartmental models do not match realistic nature. The aim of this study is to introduce time in varying model parameters, modify the classical compartmental model by improving its predictability power. Results In this study, per–capita vector density has been chosen as the time in varying model parameter. The dengue incidences, rainfall and temperature data in urban Colombo are analyzed using Fourier mathematical analysis tool. Further, periodic pattern of the reported dengue incidences and meteorological data and correlation of dengue incidences with meteorological data are identified to determine climate data–driven per–capita vector density parameter function. By considering that the vector dynamics occurs in faster time scale compares to host dynamics, a two dimensional data–driven compartmental model is derived with aid of classical compartmental models. Moreover, a function for per–capita vector density is introduced to capture the seasonal pattern of the disease according to the effect of climate factors in urban Colombo. Conclusions The two dimensional data–driven compartmental model can be used to predict weekly dengue incidences upto 4 weeks. Accuracy of the model is evaluated using relative error function and the model can be used to predict more than 75% accurate data.
The functions of the nasal cavity are very important for maintaining the internal environment of the lungs since the inner walls of the nasal cavity control the temperature and saturation of the inhaled air with water vapor until the nasopharynx is reached. In this paper, three-dimensional computational studies of airflow transport in the models of the nasal cavity were carried out for the usual inspiratory velocity in various environmental conditions. Three-dimensional numerical results are compared with experimental data and calculations of other authors. Numerical results show that during normal breathing, the human nose copes with heat and relative moisture metabolism in order to balance the intra-alveolar conditions. It is also shown in this paper that a normal nose can maintain balance even in extreme conditions, for example, in cold and hot weather. The nasal cavity accelerates heat transfer by narrowing the air passages and swirls from the nasal concha walls of the inner cavity.
Alzheimer’s disease (AD) is a degenerative brain disease. A novel agent-based modelling framework was developed in NetLogo 3D to provide fundamental insights into the potential mechanisms by which a microbe (eg. Chlamydia pneumoniae) may play a role in late-onset AD. The objective of our initial model is to simulate one possible spatial and temporal pathway of bacterial propagation via the olfactory system, which may then lead to AD symptoms. The model maps the bacteria infecting cells from the nasal cavity and the olfactory epithelium, through the olfactory bulb and into the olfactory cortex and hippocampus regions of the brain. Based on the set of biological rules, simulated randomized infection by the microbe led to the formation of beta-amyloid (Aβ) plaque and neurofibrillary (NF) tangles as well as caused immune responses. Our initial simulations demonstrated that breathing in C. pneumoniae can result in infection propagation and significant buildup of Aβ plaque and NF tangles in the olfactory cortex and hippocampus. Our model also indicated how mucosal and neural immunity can play a significant role in the pathway considered. Lower immunities, correlated with elderly individuals, had quicker and more Aβ plaque and NF tangle formation counts. In contrast, higher immunities, correlated with younger individuals, demonstrated little to no such formation. The modelling framework provides an organized visual representation of how AD progression may occur via the olfactory system to better understand disease pathogenesis. The model confirms current conclusions in available research but can be easily adjusted to match future evidence and be used by researchers for their own individual purposes. The goal of our initial model is to ultimately guide further hypothesis refinement and experimental testing to better understand the dynamic system interactions present in the etiology and pathogenesis of AD.
