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.

Causes of death in children with congenital Zika syndrome in Brazil, 2015 to 2018: A nationwide record linkage study.

BACKGROUND: Children with congenital Zika syndrome (CZS) have severe damage to the peripheral and central nervous system (CNS), greatly increasing the risk of death. However, there is no information on the sequence of the underlying, intermediate, immediate, and contributing causes of deaths among these children. The aims of this study are describe the sequence of events leading to death of children with CZS up to 36 months of age and their probability of dying from a given cause, 2015 to 2018. METHODS AND FINDINGS: In a population-based study, we linked administrative data on live births, deaths, and cases of children with CZS from the SINASC (Live Birth Information System), the SIM (Mortality Information System), and the RESP (Public Health Event Records), respectively. Confirmed and probable cases of CZS were those that met the criteria established by the Brazilian Ministry of Health. The information on causes of death was collected from death certificates (DCs) using the World Health Organization (WHO) DC template. We estimated proportional mortality (PM%) among children with CZS and among children with non-Zika CNS congenital anomalies (CA) by 36 months of age and proportional mortality ratio by cause (PMRc). A total of 403 children with confirmed and probable CZS who died up to 36 months of age were included in the study; 81.9% were younger than 12 months of age. Multiple congenital malformations not classified elsewhere, and septicemia unspecified, with 18 (PM = 4.5%) and 17 (PM = 4.2%) deaths, respectively, were the most attested underlying causes of death. Unspecified septicemia (29 deaths and PM = 11.2%) and newborn respiratory failure (40 deaths and PM = 12.1%) were, respectively, the predominant intermediate and immediate causes of death. Fetuses and newborns affected by the mother's infectious and parasitic diseases, unspecified cerebral palsy, and unspecified severe protein-caloric malnutrition were the underlying causes with the greatest probability of death in children with CZS (PMRc from 10.0 to 17.0) when compared to the group born with non-Zika CNS anomalies. Among the intermediate and immediate causes of death, pneumonitis due to food or vomiting and unspecified seizures (PMRc = 9.5, each) and unspecified bronchopneumonia (PMRc = 5.0) were notable. As contributing causes, fetus and newborn affected by the mother's infectious and parasitic diseases (PMRc = 7.3), unspecified cerebral palsy, and newborn seizures (PMRc = 4.5, each) were more likely to lead to death in children with CZS than in the comparison group. The main limitations of this study were the use of a secondary database without additional clinical information and potential misclassification of cases and controls. CONCLUSION: The sequence of causes and circumstances involved in the deaths of the children with CZS highlights the greater vulnerability of these children to infectious and respiratory conditions compared to children with abnormalities of the CNS not related to Zika.

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.

Dengue, Zika, and Chikungunya viral circulation and hospitalization rates in Brazil from 2014 to 2019: An ecological study.

BACKGROUND: In addition to their direct pathogenic effects, arthropod-borne (arboviruses) have been hypothesized to indirectly contribute to hospitalizations and death through decompensation of pre-existing comorbidities. Using nationwide data routinely collected from 1 January 2014 to 31 December 2019 in Brazil, we investigated whether local increases in arbovirus notifications were associated with excess hospitalization. METHODS: We estimated the relative risks for the association between municipality- and state-level increases in arboviral case notifications and age-standardized hospitalization rates (i.e., classified as direct or indirect based on ICD-10 codes) using Bayesian multilevel models with random effects accounting for temporal and geographic correlations. For municipality-level analyses, we excluded municipalities with <200 notifications of a given arbovirus and further adjusted the models for the local Gini Index, Human Development Index, and Family Healthcare Strategy (Estratégia de Saúde da Família) coverage. Models for dengue, Zika, and chikungunya were performed separately. RESULTS: From 2014 to 2019, Brazil registered 7,566,330 confirmed dengue cases, 159,029 confirmed ZIKV cases, and 433,887 confirmed CHIKV cases. Dengue notifications have an endemic and seasonal pattern, with cases present in 5334 of the 5570 (95.8%) Brazilian municipalities and most (69.5%) registered between February and May. Chikungunya notifications followed a similar seasonal pattern to DENV but with a smaller incidence and were restricted to 4390 (78.8%) municipalities. ZIKV was only notified in 2581 (46.3%) municipalities. Increases in dengue and chikungunya notifications were associated with small increases in age-standardized arbovirus-related hospitalizations, but no consistent association was found with all-cause or other specific indirect causes of hospitalization. Zika was associated to increases in hospitalizations by neurological diseases. CONCLUSIONS: Although we found no clear association between increased incidence of the three arboviruses and excess risks of all-cause or indirect hospitalizations at the municipality- and state-levels, follow-up investigations at the individual-level are warranted to define any potential role of acute arbovirus infection in exacerbating risks of hospitalization from underlying conditions.

Use of soil moisture active passive satellite data and WorldClim 2.0 data to predict the potential distribution of visceral leishmaniasis and its vector Lutzomyia longipalpis in Sao Paulo and Bahia states, Brazil.

Visceral leishmaniasis (VL) is a neglected tropical disease transmitted by Lutzomyia longipalpis, a sand fly widely distributed in Brazil. Despite efforts to strengthen national control programs reduction in incidence and geographical distribution of VL in Brazil has not yet been successful; VL is in fact expanding its range in newly urbanized areas. Ecological niche models (ENM) for use in surveillance and response systems may enable more effective operational VL control by mapping risk areas and elucidation of eco-epidemiologic risk factors. ENMs for VL and Lu. longipalpis were generated using monthly WorldClim 2.0 data (30-year climate normal, 1-km spatial resolution) and monthly soil moisture active passive (SMAP) satellite L4 soil moisture data. SMAP L4 Global 3-hourly 9-km EASE-Grid Surface and Root Zone Soil Moisture Geophysical Data V004 were obtained for the first image of day 1 and day 15 (0:00-3:00 hour) of each month. ENM were developed using MaxEnt software to generate risk maps based on an algorithm for maximum entropy. The jack-knife procedure was used to identify the contribution of each variable to model performance. The three most meaningful components were used to generate ENM distribution maps by ArcGIS 10.6. Similar patterns of VL and vector distribution were observed using SMAP as compared to WorldClim 2.0 models based on temperature and precipitation data or water budget. Results indicate that direct Earth-observing satellite measurement of soil moisture by SMAP can be used in lieu of models calculated from classical temperature and precipitation climate station data to assess VL risk.

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.

Case Fatality Rate Related to Microcephaly Congenital Zika Syndrome and Associated Factors: A Nationwide Retrospective Study in Brazil †.

BACKGROUND: The clinical manifestations of microcephaly/congenital Zika syndrome (microcephaly/CZS) have harmful consequences on the child's health, increasing vulnerability to childhood morbidity and mortality. This study analyzes the case fatality rate and child-maternal characteristics of cases and deaths related to microcephaly/CZS in Brazil, 2015-2017. METHODS: Population-based study developed by linkage of three information systems. We estimate frequencies of cases, deaths, case fatality rate related to microcephaly/CZS according to child and maternal characteristics and causes of death. Multivariate logistic regression models were applied. RESULTS: The microcephaly/CZS case fatality rate was 10% (95% CI 9.2-10.7). Death related to microcephaly/CZS was associated to moderate (OR = 2.15; 95% CI 1.63-2.83), and very low birth weight (OR = 3.77; 95% CI 2.20-6.46); late preterm births (OR = 1.65; 95% CI 1.21-2.23), Apgar < 7 at 1st (OR = 5.98; 95% CI 4.46-8.02) and 5th minutes (OR = 4.13; 95% CI 2.78-6.13), among others. CONCLUSIONS: A high microcephaly/CZS case fatality rate and important factors associated with deaths related to this syndrome were observed. These results can alert health teams to these problems and increase awareness about the factors that may be associated with worse outcomes.

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.

Environmental risk factors of leishmaniasis in Bahia State, Brazil using NASA Earth observation satellites / Fatores de risco ambientais para leishmanioses no estado da Bahia, Brasil, utilizando satélites da NASA de observação da Terra

NASA's Earth Observing Satellites (EOS) were used to calculate three vegetation indices, extract precipitation and elevation data, and then evaluate their applicability for assessing risk of visceral leishmaniasis (VL) and cutaneous leishmaniasis (CL) in Bahia State, Brazil. Regression models showed that either form of leishmaniasis can be predicted by NDVI, NDMI, NDWI data products and TRMM) precipitation data (R2= 0.370; p<0.001). Elevation was not significantly associated with the distribution of either VL or CL. In areas of high annual precipitation, CL was 3.6 times more likely to occur than VL. For vegetative moisture (NDMI), CL was 2.11 times more likely to occur than VL. Odds of CL occurrence increased to 5.5 times when vegetation (NDVI) and 13.5 times when liquid water content of vegetation canopies (NDWI) was considered. Areas at risk of CL and VL were mapped based on the selected explanatory variables. Accuracy of models were assessed using area under the receiver operating characteristic curve (AUC=0.72). We propose that statewide scale risk models based on use of EOS products will be a useful tool at 1 km2 spatial resolution to enable health workers to identify and target high risk areas to prevent transmission of leishmaniasis.(AU)

Os satélites de observação da Terra (SOT) da NASA foram usados para calcular três índices de vegetação, extrair dados de precipitação e elevação e avaliar sua aplicabilidade para identificar o risco para leishmaniose visceral (LV) e leishmaniose tegumentar (LT) no Estado da Bahia, Brasil. Modelos de regressão mostraram que ambas as formas de leishmaniose podem ser preditas pelos NDVI, NDMI, NDWI e precipitação TRMM (R2 = 0,370; p<0,001). A elevação não foi significativamente associada à distribuição de LV ou LT. Em áreas de alta precipitação anual, a LT foi 3,6 vezes mais provável de ocorrer do que a LV. Para a umidade vegetativa (NDMI), a LT apresentou 2,11 maior probabilidade de ocorrer do que a LV. As chances de ocorrência de LT aumentaram para 5,5 vezes em relação com a vegetação (NDVI) e 13,5 vezes quando o conteúdo de água líquida dos dosséis da vegetação (NDWI) foi considerado. Áreas em risco de LT e LV foram mapeadas com base nas variáveis explicativas selecionadas. A precisão dos modelos foi avaliada usando a área sob curva característica de operação do receptor (Curva COR=0,72). Propusemos que os modelos de risco em escala estadual baseados no uso de produtos SOT são uma ferramenta útil na resolução espacial de 1 km2 por permitir que profissionais de saúde identifiquem e direcionem áreas de alto risco para evitar a transmissão da leishmaniose. (AU)

Environmental risk factors of leishmaniasis in Bahia State, Brazil using NASA Earth observation satellites / Fatores de risco ambientais para leishmanioses no estado da Bahia, Brasil, utilizando satélites da NASA de observação da Terra

NASA’s Earth Observing Satellites (EOS) were used to calculate three vegetation indices, extract precipitation and elevation data, and then evaluate their applicability for assessing risk of visceral leishmaniasis (VL) and cutaneous leishmaniasis (CL) in Bahia State, Brazil. Regression models showed that either form of leishmaniasis can be predicted by NDVI, NDMI, NDWI data products and TRMM) precipitation data (R2 = 0.370; p<0.001). Elevation was not significantly associated with the distribution of either VL or CL. In areas of high annual precipitation, CL was 3.6 times more likely to occur than VL. For vegetative moisture (NDMI), CL was 2.11 times more likely to occur than VL. Odds of CL occurrence increased to 5.5 times when vegetation (NDVI) and 13.5 times when liquid water content of vegetation canopies (NDWI) was considered. Areas at risk of CL and VL were mapped based on the selected explanatory variables. Accuracy of models were assessed using area under the receiver operating characteristic curve (AUC=0.72). We propose that statewide scale risk models based on use of EOS products will be a useful tool at 1 km2 spatial resolution to enable health workers to identify and target high risk areas to prevent transmission of leishmaniasis.

Os satélites de observação da Terra (SOT) da NASA foram usados para calcular três índices de vegetação, extrair dados de precipitação e elevação e avaliar sua aplicabilidade para identificar o risco para leishmaniose visceral (LV) e leishmaniose tegumentar (LT) no Estado da Bahia, Brasil. Modelos de regressão mostraram que ambas as formas de leishmaniose podem ser preditas pelos NDVI, NDMI, NDWI e precipitação TRMM (R2 = 0,370; p<0,001). A elevação não foi significativamente associada à distribuição de LV ou LT. Em áreas de alta precipitação anual, a LT foi 3,6 vezes mais provável de ocorrer do que a LV. Para a umidade vegetativa (NDMI), a LT apresentou 2,11 maior probabilidade de ocorrer do que a LV. As chances de ocorrência de LT aumentaram para 5,5 vezes em relação com a vegetação (NDVI) e 13,5 vezes quando o conteúdo de água líquida dos dosséis da vegetação (NDWI) foi considerado. Áreas em risco de LT e LV foram mapeadas com base nas variáveis explicativas selecionadas. A precisão dos modelos foi avaliada usando a área sob curva característica de operação do receptor (Curva COR=0,72). Propusemos que os modelos de risco em escala estadual baseados no uso de produtos SOT são uma ferramenta útil na resolução espacial de 1 km2 por permitir que profissionais de saúde identifiquem e direcionem áreas de alto risco para evitar a transmissão da leishmaniose.

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.

Ecological niche models for sand fly species and predicted distribution of Lutzomyia longipalpis (Diptera: Psychodidae) and visceral leishmaniasis in Bahia state, Brazil.

Visceral leishmaniasis is a public health problem in Brazil. This disease is endemic in most of Bahia state, with increasing reports of cases in new areas. Ecological niche models (ENM) can be used as a tool for predicting potential distribution for disease, vectors, and to identify risk factors associated with their distribution. In this study, ecological niche models (ENMs) were developed for visceral leishmaniasis (VL) cases and 12 sand fly species captured in Bahia state. Sand fly data was collected monthly by CDC light traps from July 2009 to December 2012. MODIS satellite imagery was used to calculate NDVI, NDMI, and NDWI vegetation indices, MODIS day and night land surface temperature (LST), enhanced vegetation index (EVI), and 19 Bioclim variables were used to develop the ENM using the maximum entropy approach (Maxent). Mean diurnal range was the variable that most contributed to all the models for sand flies, followed by precipitation in wettest month. For Lutzomyia longipalpis (L. longipalpis), annual precipitation, precipitation in wettest quarter, precipitation in wettest month, and NDVI were the most contributing variables. For the VL model, the variables that contributed most were precipitation in wettest month, annual precipitation, LST day, and temperature seasonality. L. longipalpis was the species with the widest potential distribution in the state. The identification of risk areas and factors associated with this distribution is fundamental to prioritize resource allocation and to improve the efficacy of the state's program for surveillance and control of VL.