Risk of death following chikungunya virus disease in the 100 Million Brazilian Cohort, 2015-18: a matched cohort study and self-controlled case series.

BACKGROUND: Chikungunya virus outbreaks have been associated with excess deaths at the ecological level. Previous studies have assessed the risk factors for severe versus mild chikungunya virus disease. However, the risk of death following chikungunya virus disease compared with the risk of death in individuals without the disease remains unexplored. We aimed to investigate the risk of death in the 2 years following chikungunya virus disease. METHODS: We used a population-based cohort study and a self-controlled case series to estimate mortality risks associated with chikungunya virus disease between Jan 1, 2015, and Dec 31, 2018, in Brazil. The dataset was created by linking national databases for social programmes, notifiable diseases, and mortality. For the matched cohort design, individuals with chikungunya virus disease recorded between Jan 1, 2015, and Dec 31, 2018, were considered as exposed and those who were arbovirus disease-free and alive during the study period were considered as unexposed. For the self-controlled case series, we included all deaths from individuals with a chikungunya virus disease record, and each individual acted as their own control according to different study periods relative to the date of disease. The primary outcome was all-cause natural mortality up to 728 days after onset of chikungunya virus disease symptoms, and secondary outcomes were cause-specific deaths, including ischaemic heart diseases, diabetes, and cerebrovascular diseases. FINDINGS: In the matched cohort study, we included 143 787 individuals with chikungunya virus disease who were matched, at the day of symptom onset, to unexposed individuals using sociodemographic factors. The incidence rate ratio (IRR) of death within 7 days of chikungunya symptom onset was 8·40 (95% CI 4·83-20·09) as compared with the unexposed group and decreased to 2·26 (1·50-3·77) at 57-84 days and 1·05 (0·82-1·35) at 85-168 days, with IRR close to 1 and wide CI in the subsequent periods. For the secondary outcomes, the IRR of deaths within 28 days after disease onset were: 1·80 (0·58-7·00) for cerebrovascular diseases, 3·75 (1·33-17·00) for diabetes, and 3·67 (1·25-14·00) for ischaemic heart disease, and there was no evidence of increased risk in the subsequent periods. For the self-controlled case series study, 1933 individuals died after having had chikungunya virus disease and were included in the analysis. The IRR of all-cause natural death within 7 days of symptom onset of chikungunya virus disease was 8·75 (7·18-10·66) and decreased to 1·59 (1·26-2·00) at 57-84 days and 1·09 (0·92-1·29) at 85-168 days. For the secondary outcomes, the IRRs of deaths within 28 days after disease onset were: 2·73 (1·50-4·96) for cerebrovascular diseases, 8·43 (5·00-14·21) for diabetes, and 2·38 (1·33-4·26) for ischaemic heart disease, and there was no evidence of increased risk at 85-168 days. INTERPRETATION: Chikungunya virus disease is associated with an increased risk of death for up to 84 days after symptom onset, including deaths from cerebrovascular diseases, ischaemic heart diseases, and diabetes. This study highlights the need for equitable access to approved vaccines and effective anti-chikungunya virus therapeutics and reinforces the importance of robust vector-control efforts to reduce viral transmission. FUNDING: Brazilian National Research Council (CNPq), Fundação de Amparo à Pesquisa do Estado da Bahia, Wellcome Trust, and UK Medical Research Council. TRANSLATION: For the Portuguese translation of the abstract see Supplementary Materials section.

Complex network analysis of arboviruses in the same geographic domain: Differences and similarities.

Arbovirus can cause diseases with a broad spectrum from mild to severe and long-lasting symptoms, affecting humans worldwide and therefore considered a public health problem with global and diverse socio-economic impacts. Understanding how they spread within and across different regions is necessary to devise strategies to control and prevent new outbreaks. Complex network approaches have widespread use to get important insights on several phenomena, as the spread of these viruses within a given region. This work uses the motif-synchronization methodology to build time varying complex networks based on data of registered infections caused by Zika, chikungunya, and dengue virus from 2014 to 2020, in 417 cities of the state of Bahia, Brazil. The resulting network sets capture new information on the spread of the diseases that are related to the time delay in the synchronization of the time series among different municipalities. Thus the work adds new and important network-based insights to previous results based on dengue dataset in the period 2001-2016. The most frequent synchronization delay time between time series in different cities, which control the insertion of edges in the networks, ranges 7 to 14 days, a period that is compatible with the time of the individual-mosquito-individual transmission cycle of these diseases. As the used data covers the initial periods of the first Zika and chikungunya outbreaks, our analyses reveal an increasing monotonic dependence between distance among cities and the time delay for synchronization between the corresponding time series. The same behavior was not observed for dengue, first reported in the region back in 1986, either in the previously 2001-2016 based results or in the current work. These results show that, as the number of outbreaks accumulates, different strategies must be adopted to combat the dissemination of arbovirus infections.

Socioeconomic risk markers of congenital Zika syndrome: a nationwide, registry-based study in Brazil.

While it is well known that socioeconomic markers are associated with a higher risk of arbovirus infections, research on the relationship between socioeconomic factors and congenital Zika syndrome (CZS) remains limited. This study investigates the relationship between socioeconomic risk markers and live births with CZS in Brazil. We conducted a population-based study using data from all registered live births in Brazil (Live Births Information System) linked with the Public Health Event Record from 1 January 2015 to 31 December 2018. We used logistic regression models to estimate the OR and 95% CIs of CZS based on a three-level framework. In an analysis of 11 366 686 live births, of which 3353 had CZS, we observed that live births of self-identified black or mixed race/brown mothers (1.72 (95% CI 1.47 to 2.01) and 1.37 (95% CI 1.24 to 1.51)) were associated with a higher odds of CZS. Live births from single women compared with married women and those from women with less than 12 years of education compared with those with more than 12 years of education also had higher odds of CZS. In addition, live births following fewer prenatal care appointments had increased odds of CZS in the nationwide data. However, in the analyses conducted in the Northeast region (where the microcephaly epidemic started before the link with Zika virus was established and before preventive measures were known or disseminated), no statistical association was found between the number of prenatal care appointments and the odds of CZS. This study shows that live births of the most socially vulnerable women in Brazil had the greatest odds of CZS. This disproportionate distribution of risk places an even greater burden on already socioeconomically disadvantaged groups, and the lifelong disabilities caused by this syndrome may reinforce existing social and health inequalities.

Population-based surveillance for congenital zika virus syndrome: a latent class analysis of recorded cases from 2015-2018.

OBJECTIVE: This study aims to describe clinical findings and determine the medium-term survival of congenital zika syndrome (CZS) suspected cases. METHODS: A retrospective cohort study using routine register-based linked data. It included all suspected cases of CZS born in Brazil from January 1, 2015, to December 31, 2018, and followed up from birth until death, 36 months, or December 31, 2018, whichever came first. Latent class analysis was used to cluster unconfirmed cases into classes with similar combinations of anthropometry at birth, imaging findings, maternally reported rash, region, and year of birth. Kaplan-Meier curves were plotted, and Cox proportional hazards models were fitted to determine mortality up to 36 months. RESULTS: We followed 11,850 suspected cases of CZS, of which 28.3% were confirmed, 9.3% inconclusive and 62.4% unconfirmed. Confirmed cases had almost two times higher mortality when compared with unconfirmed cases. Among unconfirmed cases, we identified three distinct clusters with different mortality trajectories. The highest mortality risk was observed in those with abnormal imaging findings compatible with congenital infections (HR = 12.6; IC95%8.8-18.0) and other abnormalities (HR = 11.6; IC95%8.6-15.6) compared with those with normal imaging findings. The risk was high in those with severe microcephaly (HR = 8.2; IC95%6.4-10.6) and macrocephaly (HR = 6.6; IC95%4.5-9.7) compared with normal head size. CONCLUSION: Abnormal imaging and head circumference appear to be the main drivers of the increased mortality among suspected cases of CZS. We suggest identifying children who are more likely to die and have a greater need to optimise interventions and resource allocation regardless of the final diagnoses.

Mortality from Congenital Zika Syndrome - Nationwide Cohort Study in Brazil.

BACKGROUND: Prenatal exposure to Zika virus has potential teratogenic effects, with a wide spectrum of clinical presentation referred to as congenital Zika syndrome. Data on survival among children with congenital Zika syndrome are limited. METHODS: In this population-based cohort study, we used linked, routinely collected data in Brazil, from January 2015 through December 2018, to estimate mortality among live-born children with congenital Zika syndrome as compared with those without the syndrome. Kaplan-Meier curves and survival models were assessed with adjustment for confounding and with stratification according to gestational age, birth weight, and status of being small for gestational age. RESULTS: A total of 11,481,215 live-born children were followed to 36 months of age. The mortality rate was 52.6 deaths (95% confidence interval [CI], 47.6 to 58.0) per 1000 person-years among live-born children with congenital Zika syndrome, as compared with 5.6 deaths (95% CI, 5.6 to 5.7) per 1000 person-years among those without the syndrome. The mortality rate ratio among live-born children with congenital Zika syndrome, as compared with those without the syndrome, was 11.3 (95% CI, 10.2 to 12.4). Among infants born before 32 weeks of gestation or with a birth weight of less than 1500 g, the risks of death were similar regardless of congenital Zika syndrome status. Among infants born at term, those with congenital Zika syndrome were 14.3 times (95% CI, 12.4 to 16.4) as likely to die as those without the syndrome (mortality rate, 38.4 vs. 2.7 deaths per 1000 person-years). Among infants with a birth weight of 2500 g or greater, those with congenital Zika syndrome were 12.9 times (95% CI, 10.9 to 15.3) as likely to die as those without the syndrome (mortality rate, 32.6 vs. 2.5 deaths per 1000 person-years). The burden of congenital anomalies, diseases of the nervous system, and infectious diseases as recorded causes of deaths was higher among live-born children with congenital Zika syndrome than among those without the syndrome. CONCLUSIONS: The risk of death was higher among live-born children with congenital Zika syndrome than among those without the syndrome and persisted throughout the first 3 years of life. (Funded by the Ministry of Health of Brazil and others.).

Associations between cesarean delivery and child mortality: A national record linkage longitudinal study of 17.8 million births in Brazil.

BACKGROUND: There is an increasing use of cesarean delivery (CD) based on preference rather than on medical indication. However, the extent to which nonmedically indicated CD benefits or harms child survival remains unclear. Our hypothesis was that in groups with a low indication for CD, this procedure would be associated with higher child mortality and in groups with a clear medical indication CD would be associated with improved child survival chances. METHODS AND FINDINGS: We conducted a population-based cohort study in Brazil by linking routine data on live births between January 1, 2012 and December 31, 2018 and assessing mortality up to 5 years of age. Women with a live birth who contributed records during this period were classified into one of 10 Robson groups based on their pregnancy and delivery characteristics. We used propensity scores to match CD with vaginal deliveries (1:1) and prelabor CD with unscheduled CD (1:1) and estimated associations with child mortality using Cox regressions. A total of 17,838,115 live births were analyzed. After propensity score matching (PSM), we found that live births to women in groups with low expected frequencies of CD (Robson groups 1 to 4) had a higher death rate up to age 5 years if they were born via CD compared with vaginal deliveries (HR = 1.25, 95% CI: 1.22 to 1.28; p < 0.001). The relative rate was greatest in the neonatal period (HR = 1.39, 95% CI: 1.34 to 1.45; p < 0.001). There was no difference in mortality rate when comparing offspring born by a prelabor CD to those born by unscheduled CD. For the live births to women with a CD in a prior pregnancy (Robson group 5), the relative rates for child mortality were similar for those born by CD compared with vaginal deliveries (HR = 1.05, 95% CI: 1.00 to 1.10; p = 0.024). In contrast, for live births to women in groups with high expected rates of CD (Robson groups 6 to 10), the child mortality rate was lower for CD than for vaginal deliveries (HR = 0.90, 95% CI: 0.89 to 0.91; p < 0.001), particularly in the neonatal period (HR = 0.84, 95% CI: 0.83 to 0.85; p < 0.001). Our results should be interpreted with caution in clinical practice, since relevant clinical data on CD indication were not available. CONCLUSIONS: In this study, we observed that in Robson groups with low expected frequencies of CD, this procedure was associated with a 25% increase in child mortality. However, in groups with high expected frequencies of CD, the findings suggest that clinically indicated CD is associated with a reduction in child mortality.

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.

Classification algorithm for congenital Zika Syndrome: characterizations, diagnosis and validation.

Zika virus was responsible for the microcephaly epidemic in Brazil which began in October 2015 and brought great challenges to the scientific community and health professionals in terms of diagnosis and classification. Due to the difficulties in correctly identifying Zika cases, it is necessary to develop an automatic procedure to classify the probability of a CZS case from the clinical data. This work presents a machine learning algorithm capable of achieving this from structured and unstructured available data. The proposed algorithm reached 83% accuracy with textual information in medical records and image reports and 76% accuracy in classifying data without textual information. Therefore, the proposed algorithm has the potential to classify CZS cases in order to clarify the real effects of this epidemic, as well as to contribute to health surveillance in monitoring possible future epidemics.

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.

New Insights on the Zika Virus Arrival in the Americas and Spatiotemporal Reconstruction of the Epidemic Dynamics in Brazil.

Zika virus (ZIKV) became a worldwide public health emergency after its introduction in the Americas. Brazil was implicated as central in the ZIKV dispersion, however, a better understanding of the pathways the virus took to arrive in Brazil and the dispersion within the country is needed. An updated genome dataset was assembled with publicly available data. Bayesian phylogeography methods were applied to reconstruct the spatiotemporal history of ZIKV in the Americas and with more detail inside Brazil. Our analyses reconstructed the Brazilian state of Pernambuco as the likely point of introduction of ZIKV in Brazil, possibly during the 2013 Confederations Cup. Pernambuco played an important role in spreading the virus to other Brazilian states. Our results also underscore the long cryptic circulation of ZIKV in all analyzed locations in Brazil. Conclusions: This study brings new insights about the early moments of ZIKV in the Americas, especially regarding the Brazil-Haiti cluster at the base of the American clade and describing for the first time migration patterns within Brazil.

Interdependence between confirmed and discarded cases of dengue, chikungunya and Zika viruses in Brazil: A multivariate time-series analysis.

The co-circulation of different arboviruses in the same time and space poses a significant threat to public health given their rapid geographic dispersion and serious health, social, and economic impact. Therefore, it is crucial to have high quality of case registration to estimate the real impact of each arboviruses in the population. In this work, a Vector Autoregressive (VAR) model was developed to investigate the interrelationships between discarded and confirmed cases of dengue, chikungunya, and Zika in Brazil. We used data from the Brazilian National Notifiable Diseases Information System (SINAN) from 2010 to 2017. There were three peaks in the series of dengue notification in this period occurring in 2013, 2015 and in 2016. The series of reported cases of both Zika and chikungunya reached their peak in late 2015 and early 2016. The VAR model shows that the Zika series have a significant impact on the dengue series and vice versa, suggesting that several discarded and confirmed cases of dengue could actually have been cases of Zika. The model also suggests that the series of confirmed and discarded chikungunya cases are almost independent of the cases of Zika, however, affecting the series of dengue. In conclusion, co-circulation of arboviruses with similar symptoms could have lead to misdiagnosed diseases in the surveillance system. We argue that the routinely use of mathematical and statistical models in association with traditional symptom-surveillance could help to decrease such errors and to provide early indication of possible future outbreaks. These findings address the challenges regarding notification biases and shed new light on how to handle reported cases based only in clinical-epidemiological criteria when multiples arboviruses co-circulate in the same population.

Impact evaluation of Zika epidemic on congenital anomalies registration in Brazil: An interrupted time series analysis.

This study aimed to assess the impact of the Zika epidemic on the registration of birth defects in Brazil. We used an interrupted time series analysis design to identify changes in the trends in the registration of congenital anomalies. We obtained monthly data from Brazilian Live Birth Information System and used two outcome definitions: 1) rate of congenital malformation of the brain and eye (likely to be affected by Zika and its complications) 2) rate of congenital malformation not related to the brain or eye unlikely to be causally affected by Zika. The period between maternal infection with Zika and diagnosis of congenital abnormality attributable to the infection is around six months. We therefore used September 2015 as the interruption point in the time series, six months following March 2015 when cases of Zika started to increase. For the purposes of this analysis, we considered the period from January 2010 to September 2015 to be "pre-Zika event," and the period from just after September 2015 to December 2017 to be "post-Zika event." We found that immediately after the interruption point, there was a great increase in the notification rate of congenital anomalies of 14.9/10,000 live births in the brain and eye group and of 5.2/10,000 live births in the group not related with brain or eye malformations. This increase in reporting was in all regions of the country (except in the South) and especially in the Northeast. In the period "post-Zika event", unlike the brain and eye group which showed a monthly decrease, the group without brain or eye malformations showed a slow but significant increase (relative to the pre-Zika trend) of 0.2/10,000 live births. These findings suggest an overall improvement in the registration of birth malformations, including malformations that were not attributed to Zika, during and after the Zika epidemic.

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.

Application of spatio-temporal scan statistics for the detection of areas with increased risk for American visceral leishmaniasis in the state of Bahia, Brazil.

Mirroring the global increase of registered cases of American visceral leishmaniasis (AVL), this infection has become a growing public health problem in Brazil during the last several years. As the traditional approach to control employed by the governmental health agencies has failed to reduce the incidence and epidemic outbreaks of this illness, we propose a re-evaluation of the national strategy of intervention and monitoring. Our thinking is based on a series of spatio-temporal scan statistics of the west-central region of the state of Bahia covering the 11-year period from 1994 to 2004. By analyzing the situation, spatially and temporally, we show that the disease is a not only a growing focal threat but that it is also appearing in the form of endemic clusters in the cities. The areas where the disease has been found have been classified according to the degree of risk of infection for humans and canines. The overall objective of this study was to identify areas of increased risk of AVL, including its seasonality, and to suggest ways and means to improve the detection of the disease. The findings presented here should not only be of interest for the efforts to control AVL in the study area but also be useful for developing control strategies in other endemic regions of Brazil.