Applying a multi-strain dengue model to epidemics data

Dengue disease transmission is a complex vector-borne disease, mainly due to the co-circulation of four serotypes of the virus. Mathematical models have proved to be a useful tool to understand the complexity of this disease. In this work, we extend the model studied by Esteva et al., 2003, originally proposed for two serotypes, to four circulating serotypes. Using epidemic data of dengue fever in Iquitos (Peru) and San Juan (Puerto Rico), we estimate numerically the co-circulation parameter values for selected outbreaks using a bootstrap method, and we also obtained the Basic Reproduction Number, R₀, for each serotype, using both analytical calculations and numerical simulations. Our results indicate that the impact of co-circulation of serotypes in population dynamics of dengue infection is such that there is a reduced effect from DENV-3 to DENV-4 in comparison to no-cross effect for epidemics in Iquitos. Concerning San Juan epidemics, also comparing to no-cross effect, we also observed a reduced effect from the predominant serotype DENV-3 to both DENV-2 and DENV-1 epidemics neglecting the very small number of cases of DENV-4.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

Assessing the nationwide impact of COVID-19 mitigation policies on the transmission rate of SARS-CoV-2 in Brazil

COVID-19 is now identified in almost all countries in the world, with poorer regions being particularly more disadvantaged to efficiently mitigate the impacts of the pandemic. In the absence of efficient therapeutics or large-scale vaccination, control strategies are currently based on non-pharmaceutical interventions, comprising changes in population behavior and governmental interventions, among which the prohibition of mass gatherings, closure of non-essential establishments, quarantine and movement restrictions. In this work we analyzed the effects of 707 governmental interventions published up to May 22, 2020, and population adherence thereof, on the dynamics of COVID-19 cases across all 27 Brazilian states, with emphasis on state capitals and remaining inland cities. A generalized SEIR (Susceptible, Exposed, Infected and Removed) model with a time-varying transmission rate (TR), that considers transmission by asymptomatic individuals, is presented. We analyze the effect of both the extent of enforced measures across Brazilian states and population movement on the changes in the TR and effective reproduction number. The social mobility reduction index, a measure of population movement, together with the stringency index, adapted to incorporate the degree of restrictions imposed by governmental regulations, were used in conjunction to quantify and compare the effects of varying degrees of policy strictness across Brazilian states. Our results show that population adherence to social distance recommendations plays an important role for the effectiveness of interventions and represents a major challenge to the control of COVID-19 in low- and middle-income countries.

Assessing the nationwide impact of COVID-19 mitigation policies on the transmission rate of SARS-CoV-2 in Brazil

COVID-19 is now identified in almost all countries in the world, with poorer regions being particularly more disadvantaged to efficiently mitigate the impacts of the pandemic. In the absence of efficient therapeutics or large-scale vaccination, control strategies are currently based on non-pharmaceutical interventions, comprising changes in population behavior and governmental interventions, among which the prohibition of mass gatherings, closure of non-essential establishments, quarantine and movement restrictions. In this work we analyzed the effects of 707 governmental interventions published up to May 22, 2020, and population adherence thereof, on the dynamics of COVID-19 cases across all 27 Brazilian states, with emphasis on state capitals and remaining inland cities. A generalized SEIR (Susceptible, Exposed, Infected and Removed) model with a time-varying transmission rate (TR), that considers transmission by asymptomatic individuals, is presented. We analyze the effect of both the extent of enforced measures across Brazilian states and population movement on the changes in the TR and effective reproduction number. The social mobility reduction index, a measure of population movement, together with the stringency index, adapted to incorporate the degree of restrictions imposed by governmental regulations, were used in conjunction to quantify and compare the effects of varying degrees of policy strictness across Brazilian states. Our results show that population adherence to social distance recommendations plays an important role for the effectiveness of interventions and represents a major challenge to the control of COVID-19 in low- and middle-income countries.

Mathematical modeling of COVID-19 in 14.8 million individuals in Bahia, Brazil

COVID-19 is affecting healthcare resources worldwide, with lower and middle-income countries being particularly disadvantaged to mitigate the challenges imposed by the disease, including the availability of a sufficient number of infirmary/ICU hospital beds, ventilators, and medical supplies. Here, we use mathematical modelling to study the dynamics of COVID-19 in Bahia, a state in northeastern Brazil, considering the influences of asymptomatic/non-detected cases, hospitalizations, and mortality. The impacts of policies on the transmission rate were also examined. Our results underscore the difficulties in maintaining a fully operational health infrastructure amidst the pandemic. Lowering the transmission rate is paramount to this objective, but current local efforts, leading to a 36% decrease, remain insufficient to prevent systemic collapse at peak demand, which could be accomplished using periodic interventions. Non-detected cases contribute to a ∽55% increase in R0. Finally, we discuss our results in light of epidemiological data that became available after the initial analyses.

Investigating extradomiciliary transmission of tuberculosis: An exploratory approach using social network patterns of TB cases and controls and the genotyping of Mycobacterium tuberculosis

Extradomiciliary contacts have been overlooked in the study of TB transmission due to difficulties in identifying actual contacts in large populations. Complex network analysis provides a framework to model the structure of contacts, specially extradomiciliary ones. We conducted a study of incident sputum-positive TB cases and healthy controls occurring in a moderate TB burden city. Cases and controls were interviewed to obtain data regarding the usual locations of residence, work, study, and leisure. Mycobacterium tuberculosis isolated from sputum was genotyped. The collected data were used to build networks based on a framework of putative social interactions indicating possible TB transmission. A user-friendly open source environment (GraphTube) was setup to extract information from the collected data. Networks based on the likelihood of patient-patient, patient-healthy, and healthy-healthy contacts were setup, depending on a constraint of geographical distance of places attended by the volunteers. Using a threshold for the geographical distance of 300 m, the differences between TB cases and controls are revealed. Several clusters formed by social network nodes with high genotypic similarity were characterized. The developed framework provided consistent results and can be used to support the targeted search of potentially infected individuals and to help to understand the TB transmission.

Heterogeneities in dengue spatial-temporal transmission in Brazilian cities and its influence on the optimal age of vaccination

Background

The development of a safe and effective vaccine is considered crucial for dengue transmission control since vetor control has been failed; some potential candidates are currently in test, and in this context theoretical studies are necessary to evaluate vaccination strategies such as the age groups that should be vaccinated, the percentage of the population at risk, and the target geographic regions to make dengue control feasible and optimal.

Methods

A partial differential model is used to mimics dengue transmission in human population in order to estimate the optimal vaccination age, using data collected from dengue reported cases in ten cities of Brazil from 2001 to 2014. For this purpose, the basic reproduction number of the disease was minimized assuming a single-dose vaccination strategy, equal vaccine efficacy for all circulating serotypes, and no vaccine failure. Numerical methods were used to assess the optimal vaccination age and its confidence age range.

Results

The results reveal complex spatial-temporal patterns associated to the disease transmission, highlighting the heterogeneity in defining the target population for dengue vaccination. However, the values obtained for the optimal age of vaccination, as targeting individuals under 13 years old, are compatible with the ones reported in similar studies in Brazil. The results also show that the optimal age for vaccination in general does not match with the age of the highest number of cases.

Conclusions

The variation of the optimal age for vaccination across the country reflects heterogeneities in dengue spatial-temporal transmission in Brazilian cities, and can be used to define the target population and cities to optimize vaccination strategies in a context of high cost and low quantity of available vaccine.

Comparison of complex networks and tree-based methods of phylogenetic analysis and proposal of a bootstrap method

Complex networks have been successfully applied to the characterization and modeling of complex systems in several distinct areas of Biological Sciences. Nevertheless, their utilization in phylogenetic analysis still needs to be widely tested, using different molecular data sets and taxonomic groups, and, also, by comparing complex networks approach to current methods in phylogenetic analysis. In this work, we compare all the four main methods of phylogenetic analysis (distance, maximum parsimony, maximum likelihood, and Bayesian) with a complex networks method that has been used to provide a phylogenetic classification based on a large number of protein sequences as those related to the chitin metabolic pathway and ATP-synthase subunits. In order to perform a close comparison to these methods, we selected Basidiomycota fungi as the taxonomic group and used a high-quality, manually curated and characterized database of chitin synthase sequences. This enzymatic protein plays a key role in the synthesis of one of the exclusive features of the fungal cell wall: the presence of chitin. The communities (modules) detected by the complex network method corresponded exactly to the groups retrieved by the phylogenetic inference methods. Additionally, we propose a bootstrap method for the complex network approach. The statistical results we have obtained with this method were also close to those obtained using traditional bootstrap methods.

Stochastic dynamics for reinfection by transmitted diseases

The use of stochastic models to study the dynamics of infectious diseases is an important tool to understand the epidemiological process. For several directly transmitted diseases, reinfection is a relevant process, which can be expressed by endogenous reactivation of the pathogen or by exogenous reinfection due to direct contact with an infected individual (with smaller reinfection rate σβ than infection rate β). In this paper, we examine the stochastic susceptible, infected, recovered, infected (SIRI) model simulating the endogenous reactivation by a spontaneous reaction, while exogenous reinfection by a catalytic reaction. Analyzing the mean-field approximations of a site and pairs of sites, and Monte Carlo (MC) simulations for the particular case of exogenous reinfection, we obtained continuous phase transitions involving endemic, epidemic, and no transmission phases for the simple approach; the approach of pairs is better to describe the phase transition from endemic phase (susceptible, infected, susceptible (SIS)-like model) to epidemic phase (susceptible, infected, and removed or recovered (SIR)-like model) considering the comparison with MC results; the reinfection increases the peaks of outbreaks until the system reaches endemic phase. For the particular case of endogenous reactivation, the approach of pairs leads to a continuous phase transition from endemic phase (SIS-like model) to no transmission phase. Finally, there is no phase transition when both effects are taken into account. We hope the results of this study can be generalized for the susceptible, exposed, infected, and removed or recovered (SEIR_{I}^{E}) model, for which the state exposed (infected but not infectious), describing more realistically transmitted diseases such as tuberculosis. In future work, we also intend to investigate the effect of network topology on phase transitions when the SIRI model describes both transmitted diseases (σ<1) and social contagions (σ>1).

What are the Evolutionary Origins of Mitochondria? A Complex Network Approach

Mitochondria originated endosymbiotically from an Alphaproteobacteria-like ancestor. However, it is still uncertain which extant group of Alphaproteobacteria is phylogenetically closer to the mitochondrial ancestor. The proposed groups comprise the order Rickettsiales, the family Rhodospirillaceae, and the genus Rickettsia. In this study, we apply a new complex network approach to investigate the evolutionary origins of mitochondria, analyzing protein sequences modules in a critical network obtained through a critical similarity threshold between the studied sequences. The dataset included three ATP synthase subunits (4, 6, and 9) and its alphaproteobacterial homologs (b, a, and c). In all the subunits, the results gave no support to the hypothesis that Rickettsiales are closely related to the mitochondrial ancestor. Our findings support the hypothesis that mitochondria share a common ancestor with a clade containing all Alphaproteobacteria orders, except Rickettsiales.

Interpreting measures of tuberculosis transmission: a case study on the Portuguese population

Background

Tuberculosis remains a high burden for Human society despite considerable investments in its control. Unique features in the history of infection and transmission dynamics of tuberculosis pose serious limitations on the direct interpretation of surveillance data and call for models that incorporate latent processes and simulate specific interventions.

Methods

A transmission model was adjusted to the dataset of active tuberculosis cases reported in Portugal between 2002 and 2009. We estimated key transmission parameters from the data (i.e. time to diagnosis, treatment length, default proportion, proportion of pulmonary TB cases). Using the adjusted model to the Portuguese case, we estimated the total burden of tuberculosis in Portugal. We further performed sensitivity analysis to heterogeneities in susceptibility to infection and exposure intensity.

Results

We calculated a mean time to diagnose of 2.81 months and treatment length of 8.80 months in Portugal. The proportion defaulting treatment was calculated as 0.04 and the proportion of pulmonary cases as 0.75. Using these values, we estimated a TB burden of 1.6 million infected persons, corresponding to more than 15% of the Portuguese population. We further described the sensitivity of these estimates to heterogeneity.

Conclusions

We showed that the model reproduces well the observed dynamics of the Portuguese data, thus demonstrating its adequacy for devising control strategies for TB and predicting the effects of interventions.

Stochastic dynamics of dengue epidemics

We use a stochastic Markovian dynamics approach to describe the spreading of vector-transmitted diseases, such as dengue, and the threshold of the disease. The coexistence space is composed of two structures representing the human and mosquito populations. The human population follows a susceptible-infected-recovered (SIR) type dynamics and the mosquito population follows a susceptible-infected-susceptible (SIS) type dynamics. The human infection is caused by infected mosquitoes and vice versa, so that the SIS and SIR dynamics are interconnected. We develop a truncation scheme to solve the evolution equations from which we get the threshold of the disease and the reproductive ratio. The threshold of the disease is also obtained by performing numerical simulations. We found that for certain values of the infection rates the spreading of the disease is impossible, for any death rate of infected mosquitoes.

Detecting network communities: an application to phylogenetic analysis

This paper proposes a new method to identify communities in generally weighted complex networks and apply it to phylogenetic analysis. In this case, weights correspond to the similarity indexes among protein sequences, which can be used for network construction so that the network structure can be analyzed to recover phylogenetically useful information from its properties. The analyses discussed here are mainly based on the modular character of protein similarity networks, explored through the Newman-Girvan algorithm, with the help of the neighborhood matrix . The most relevant networks are found when the network topology changes abruptly revealing distinct modules related to the sets of organisms to which the proteins belong. Sound biological information can be retrieved by the computational routines used in the network approach, without using biological assumptions other than those incorporated by BLAST. Usually, all the main bacterial phyla and, in some cases, also some bacterial classes corresponded totally (100%) or to a great extent (>70%) to the modules. We checked for internal consistency in the obtained results, and we scored close to 84% of matches for community pertinence when comparisons between the results were performed. To illustrate how to use the network-based method, we employed data for enzymes involved in the chitin metabolic pathway that are present in more than 100 organisms from an original data set containing 1,695 organisms, downloaded from GenBank on May 19, 2007. A preliminary comparison between the outcomes of the network-based method and the results of methods based on Bayesian, distance, likelihood, and parsimony criteria suggests that the former is as reliable as these commonly used methods. We conclude that the network-based method can be used as a powerful tool for retrieving modularity information from weighted networks, which is useful for phylogenetic analysis.

Complex weighted networks have been applied to uncover organizing principles of complex biological, technological, and social systems. We propose herein a new method to identify communities in such structures and apply it to phylogenetic analysis. Recent studies using this theory in genomics and proteomics contributed to the understanding of the structure and dynamics of cellular complex interaction webs. Three main distinct molecular networks have been investigated based on transcriptional and metabolic activity, and on protein interaction. Here we consider the evolutionary relationship between proteins throughout phylogeny, employing the complex network approach to perform a comparative study of the enzymes related to the chitin metabolic pathway. We show how the similarity index of protein sequences can be used for network construction, and how the underlying structure is analyzed by the computational routines of our method to recover useful and sound information for phylogenetic studies. By focusing on the modular character of protein similarity networks, we were successful in matching the identified networks modules to main bacterial phyla, and even some bacterial classes. The network-based method reported here can be used as a new powerful tool for identifying communities in complex networks, retrieving useful information for phylogenetic studies.

Modelling the dynamics of dengue real epidemics

In this work, we use a mathematical model for dengue transmission with the aim of analysing and comparing two dengue epidemics that occurred in Salvador, Brazil, in 1995-1996 and 2002. Using real data, we obtain the force of infection, Λ, and the basic reproductive number, R(0), for both epidemics. We also obtain the time evolution of the effective reproduction number, R(t), which results in a very suitable measure to compare the patterns of both epidemics. Based on the analysis of the behaviour of R(0) and R(t) in relation to the adult mosquito control parameter of the model, we show that the control applied only to the adult stage of the mosquito population is not sufficient to stop dengue transmission, emphasizing the importance of applying the control to the aquatic phase of the mosquito.

Periodic forcing in a three-level cellular automata model for a vector-transmitted disease

A periodically forced two-dimensional cellular automata model is used to reproduce and analyze the complex spatiotemporal patterns observed in the transmission of vector infectious diseases. The system, which comprises three population levels, is introduced to describe complex features of the dynamics of the vector-transmitted dengue epidemics, known to be very sensitive to seasonal variables. The three coupled levels represent the human, the adult, and immature vector populations. The dynamics includes external seasonality forcing, human and mosquito mobility, and vector control effects. The model parameters, even if bounded to well-defined intervals obtained from reported data, can be selected to reproduce specific epidemic outbursts. In the current study, explicit results are obtained by comparison with actual data retrieved from the time series of dengue epidemics in two cities in Brazil. The results show fluctuations that are not captured by mean-field models. It also reveals the qualitative behavior of the spatiotemporal patterns of the epidemics. In the extreme situation of the absence of external periodic drive, the model predicts a completely distinct long-time evolution. The model is robust in the sense that it is able to reproduce the time series of dengue epidemics of different cities, provided that the forcing term takes into account the local rainfall modulation. Finally, an analysis is provided of the effect of the dependence between epidemics threshold and vector control actions, both in the presence and absence of human mobility factor.

Tsallis scaling and the long-range Ising chain: A transfer matrix approach

A numerically efficient transfer matrix (TM) approach is introduced to investigate the long-range Ising spin chain. Results obtained within this procedure are primarily used to verify the Tsallis scaling hypothesis for long-range systems with an alpha power-law decay of the interaction constants, both in the extensive (alpha>1) and nonextensive (alpha<1) regimes. Results for finite-size systems, taking into account all interactions between spins up to 24 sites apart, show that the conjecture is satisfied with a very good precision (less than 0.004%) for all temperature intervals. This TM procedure is further used to investigate several other thermodynamic and critical properties of this system, and it may also be extended to similar one-dimensional long-range systems.