Skip to main content

Zika virus disease, microcephaly and Guillain-Barré syndrome in Colombia: epidemiological situation during 21 months of the Zika virus outbreak, 2015–2017

Abstract

Background

The Zika virus disease (ZVD) has had a huge impact on public health in Colombia for the numbers of people affected and the presentation of Guillain-Barre syndrome (GBS) and microcephaly cases associated to ZVD.

Methods

A retrospective descriptive study was carried out, we analyze the epidemiological situation of ZVD and its association with microcephaly and GBS during a 21-month period, from October 2015 to June 2017. The variables studied were: (i) ZVD cases, (ii) ZVD cases in pregnant women, (iii) laboratory-confirmed ZVD in pregnant women, (iv) ZVD cases associated with microcephaly, (v) laboratory-confirmed ZVD associated with microcephaly, and (vi) ZVD associated to GBS cases. Average number of cases, attack rates (AR) and proportions were also calculated. The studied variables were plotted by epidemiological weeks and months. The distribution of ZVD cases in Colombia was mapped across the time using Kernel density estimator and QGIS software; we adopted Kernel Ridge Regression (KRR) and the Gaussian Kernel to estimate the number of Guillain Barre cases given the number of ZVD cases.

Results

One hundred eight thousand eighty-seven ZVD cases had been reported in Colombia, including 19,963 (18.5%) in pregnant women, 710 (0.66%) associated with microcephaly (AR, 4.87 cases per 10,000 live births) and 453 (0.42%) ZVD associated to GBS cases (AR, 41.9 GBS cases per 10,000 ZVD cases). It appears the cases of GBS increased in parallel with the cases of ZVD, cases of microcephaly appeared 5 months after recognition of the outbreak. The kernel density map shows that throughout the study period, the states most affected by the Zika outbreak in Colombia were mainly San Andrés and Providencia islands, Casanare, Norte de Santander, Arauca and Huila. The KRR shows that there is no proportional relationship between the number of GBS and ZVD cases. During the cross validation, the RMSE achieved for the second order polynomial kernel, the linear kernel, the sigmoid kernel, and the Gaussian kernel are 9.15, 9.2, 10.7, and 7.2 respectively.

Conclusions

This study updates the epidemiological analysis of the ZVD situation in Colombia describes the geographical distribution of ZVD and shows the functional relationship between ZVD cases and GBS.

Peer Review reports

Background

In recent decades, Colombia has been an exceptional witness to an avalanche of emerging arboviruses. Outbreaks of Dengue, Venezuelan equine encephalitis, Chikungunya virus -associated arthralgia- and most recently Zika virus disease (ZVD) have caused significant health burden in Colombia. Other emerging viruses have been detected with unknown disease burden, such as West Nile virus, St. Louis encephalitis virus and hantavirus [1]. Guillain-Barre Syndrome (GBS), microcephaly and encephalitis due to arboviruses are novel and emerging public health issues in Colombia [2].

By March 10, 2017, autochthonous cases and vector-borne transmission of ZVD have been confirmed in 35 countries and 13 regions of the Americas [3]. In Latin America, the countries with higher ZVD cases rate (i.e., [autochthonous suspected + autochthonous confirmed]/100,000 pop) are Honduras (396), Belize (239.62), Colombia (221.29), Venezuela (197.72), El Salvador (187.44) and Brazil (168.09) [3, 4]. In Colombia and Latin American, during the last months, the infection with ZVD has had a huge impact on public health and its effects continue to be felt [2].

In October 2015, the first nine cases of ZVD in Colombia were reported in Turbaco (Bolivar State), in the Caribbean area [5]. So far, June 24th, 2017 and nearly 21 months later the disease has spread to the entire country where Aedes aegypti is established. Similar to other countries affected by Zika virus, the increases of congenital syndromes including microcephaly cases and GBS appeared abruptly coincident with the ZVD outbreak [6,7,8,9,10,11,12].

In Colombia, the National Health Institute (INS) conducts intensified surveillance for microcephaly and GBS cases associated with ZVD [13]. As a result, the INS has published preliminary reports about ZVD cases [8], GBS [8] and microcephaly cases associated with ZVD [8], other reports about transmission and epidemiology of Zika virus in Colombia have also been published [9,10,11]. In this study, we analyze the epidemiological situation of ZVD and its association with microcephaly and GBS in Colombia during a 21-month period, from October 2015 to June 2017.

Methods

Type of study and population

We conducted a retrospective descriptive epidemiological study of people with symptoms of ZVD reported by INS. The study included 21 months (from epiweek 40 beginning 10th October 2015, until epiweek 24, beginning 25th June 2017) of notifications of ZVD cases in 37 reporting administrative units of Colombia (32 states and 5 districts).

Data collection and study variables

The public and freely available weekly epidemiological datasets from the National Health Institute of Colombia were used [13,14,15]. We imported these datasets into Microsoft Excel, where it was organized by notification areas (States and Districts). Then, from the accumulated cases of each week, we estimated the cases of the previous one. The variables studied were: (i) ZVD cases, (ii) ZVD cases in pregnant women, (iii) laboratory-confirmed ZVD in pregnant women, (iv) ZVD cases associated with microcephaly, (v) laboratory-confirmed ZVD associated with microcephaly, and (vi) ZVD associated to GBS cases.

For the present study, ZVD cases comprise the patients who were notified to the INS with symptoms of ZVD (fever and rash, and one or more of the following symptoms which are not explained by other medical conditions: non-purulent conjunctivitis, arthralgia’s, myalgia, headache or malaise), with and without laboratory confirmation. Laboratory-confirmed ZVD was defined as the presence of clinical symptoms of ZVD and a positive result for Zika virus RNA on RT-PCR assay [2, 9, 13]. A ZVD cases associated with microcephaly is any product of gestation (abortions, voluntary interruptions of pregnancy, stillbirths and live births) up to 12 months of age, which presents with a present a structural or functional alteration of the central nervous system, like microcephaly and others abnormalities of central nervous system, occurring during the gestation of a woman with symptoms of ZVD, with or without laboratory confirmation [2, 8]. A ZVD associated to GBS cases, is a person with laboratory-confirmed ZVD or with symptoms of ZVD, living in an endemic area of residence or travel within 15 days of symptom onset to an area with laboratory-confirmed Zika virus (ZIKV) circulation [7, 14, 15].

Population data and trends of Colombian population 2015–2017 were obtained by DANE (Statistics and data of Colombian government) [16]. The trends of live births ≤1 years was used as estimator of indirect of pregnant women and the data of live births was also taken from DANE [17, 18].

Data analysis

Descriptive analysis of epidemiological situation of ZVD

ZVD transmission was analyzed by plotting the above mentioned variables against epidemiologic week/month in Microsoft Excel. The distribution of ZVD cases per gender/groups of age and the distribution per states were presented as percentage and attack rates (per each 100,000 population). Microcephaly cases were presented as a proportion over the total of live births during the study period, besides, GBS cases were presented as a proportion over the total ZVD cases.

Geographic distribution of ZVD cases

A spatial analysis was conducted by importing the dataset and shapefiles of Colombian maps [19] into QGIS 2.14.3 (Development Team - Open Source Geospatial Foundation Project, 2016). With this analysis the goal was to describe the geographical distribution of ZVD cases (based on notification areas) and show areas with high concentration of cases or high intensity of cases per unit of surface. In order to assess the intensity of incidence per unit area (km2) we applied the density kernel technique, this technique produce a smoothed surface with the representation of disease distribution [19]. We used a bandwidth of 120 km. The maps were visualized for each trimester of each year and a summarized map for the entire period of study; moreover, to see the intensity of Zika incidences in each map, we established seven categories using as criteria the classification per quantile, the two highest categories show the clusters.

Functional relationship between ZVD indicator and GBS cases

Kernel ridge regression (KRR) was used to estimate the number of GBS cases given the number of ZVD cases. We use this machine learning model with the kernel trick to disclose the functional dependence between the target variable, i.e., the number of GBS cases, and the independent observable variable, i.e., the number of ZVD cases.

In Fig. 3a and b, are depicted the number GBS cases (the target variable) vs ZVD cases (the independent variable), as we can see there is no linear relation among both variables, then a linear regression model underfits (Fig. B), It will fail in the prediction with high error. In such case, there are two options: firstly, we directly map the independent variables into a higher dimensional space, thereafter, the linear ridge regression is applied. The second option is adopting the kernel trick, which allows us to solve the regression problem by applying only inner product between independent variables without mapping them, thereby, the kernel trick is computationally less costly than the first option, because the original space of independent variables has less dimensions than the higher dimensional space of mapped variables. Therefore, we adopted KRR as a regression model. Saunders et al. [20], provide a broader explanation of this model.

Moreover, KRR solves the regression problem in an analytic fashion instead of an algorithmic one. The latter implies tuning more parameters than the analytic solution, e.g., the learning rate and the maximum number of iterations in the gradient descent algorithm. Furthermore, we used the R language and the CVST package to implement, tune, and test the KRR model.

On the other hand, in order to set up the parameters and choose the kernel function, we carried out the 5-fold cross validation technique, where the data set was randomly split in 5 folds, four of them for training and one fold for validation. This was done five times, changing the elements contained into each fold. The accuracy of the regression model was measured through the Root Mean Square Error and the coefficient of determination as well.

Results

Descriptive analysis of epidemiological situation of ZVD

From 10th October 2015 to 24th June 2017, 108,087 ZVD cases were notified by INS, Table 1 shows the distribution of gender and age groups. 66.2% of cases occurred in females, the attack rate was higher among women (286 per 100,000 population) than men (149 per 100,000 population). The highest AR was seen in the groups 25–29 years (375 per 100,000 population), 30–34 years (365 per 100,000 population), 20–24 years (326 per 100,000 population) and 35–39 years (305 per 100,000 population), representing 49.0% of the total cases (Table 1).

Table 1 Attack rates of Zika virus disease cases, by gender and age groups, Colombia, October 10th, 2015-June, 24th 2017

Monthly and weekly distribution of ZVD cases (n = 108,087), including 19,963 (18.5%) in pregnant women, 710 (0.66%) associated with microcephaly (4.87 cases per 10,000 live births), and 453 (0.42%) ZVD associated to GBS cases (41.9 GBS cases per 10,000 ZVD cases) are depicted in Fig. 1. During this period of time, 9802 (9.0%) ZVD cases was confirmed by ZIKV-specific RT-PCR, 6365 (5.91%) in pregnant women and 174 (0.16%) in cases of microcephaly (1.2 microcephaly cases per 10,000 live births).

Fig. 1
figure 1

a Colombia, monthly distribution of Zika virus disease (ZVD) cases, microcephaly cases and Guillain-Barré Syndrome, 2015–2017. b Colombia, weekly distribution of ZVD cases, microcephaly cases and Guillain-Barré Syndrome, 2015–2017. In both figures, the left axis represents the ZVD cases, ZVD cases in pregnant women and laboratory-confirmed (LC) ZVD in pregnant women; the right axis represents the ZVD cases associated with microcephaly, laboratory-confirmed ZVD associated with microcephaly and ZVD associated to GBS cases

Outbreak was first reported in October, 2015 (epiweek 40), in 2016 presented a largest number of ZVD cases (87.8%; 94,947 cases), and currently in 2017 (epiweek 1–25, January to June), have a median of 52 cases (range 15–123) (Fig. 1). The peak of the ZVD cases and ZVD cases in pregnant women was in February 2016 (epiweek 5–8), and was preceded by an increase of GBS cases during the period of November 2015 to January 2016 (peak of GBS in January 2016) (Fig. 1). Five months after the outbreak started, the microcephaly cases were increased in association with ZVD cases in the country (Fig. 1),

Table 2 shows population distribution of ZVD and ARs in Colombia, pregnant women diagnosed by laboratory are mainly concentrated in Norte de Santander, Huila, Meta, Casanare, Arauca, Caquetá, Tolima, Santander and Valle del Cauca.

Table 2 Attack rates per 100,000 population distribution of Zika virus disease in the total Colombian population and pregnant women, October 10th, 2015-June, 24th 2017

Geographic distribution of ZVD cases

The Kernel density map (Fig. 2) shows the areas with high concentration of cases during each quarter of each year and throughout the study period (October 2015–June 2017). In 2015, cases were concentrated in two geographic regions of the country; in the Caribbean region, mainly in states of San Andrés and Providencia Islands (not included in the map), Magdalena, Sucre, Atlantico, Bolívar (the state where the first Zika cases appeared in October 2015) and Córdoba. Another affected area was the Andean region, mainly in the state of Norte of Santander, where the intensity was higher, compared to Cundinamarca, Huila and Tolima.

Fig. 2
figure 2

Geographically distribution of Zika cases in Colombia (October 10th, 2015-June, 24th 2017) using Kernel density estimator maps

In 2016, all geographic regions of the country showed at least one state with high concentration of cases, the Andean region and Orinoquia were the most affected during the four quarters (Fig. 2). In the first quarter of 2016, the cases were concentrated in four states in the Caribbean region (cluster in San Andres and Providencia, Atlantico, Magdalena and Córdoba), four in the Andean region (cluster in Norte de Santander, Huila and Tolima, Santander), 3 in the Orinoquia region (cluster in Casanare and Arauca, Meta), 2 in the Amazonia region (Cluster in Amazonas and Caqueta) and 1 from the Pacifica region (Valle del Cauca) (Fig. 2). In the second quarter of 2016, the number of clusters is decreased in comparison with the previous quarter, so clusters were located in Casanare, Arauca, Valle del Cauca and Santander. In the same quarter the number of ZVD cases also decreased in several areas, however, in Meta, Guaviare, Vichada, Quindio and Cauca an increased was seen. Between the third quarter of 2016 and the second one of 2017, the intensity of cases decreased in the whole country.

In general, throughout the study period (October 2015–June 2017), the states San Andrés and Providencia, Casanare, Norte de Santander, Arauca, Huila, Valle del Cauca, Tolima, Santander, Amazonas, Meta, Atlantico, Magdalena, Caquetá, Cundinamarca, Córdoba, Sucre, Guaviare, Cesar and Putumayo were most affected by ZVD in Colombia (Fig. 2).

Functional relationship between ZVD indicator and GBS cases

Figure 1 suggests that there is a relationship between GBS and ZVD cases in the country due to GBS cases started happening more frequently when ZVD cases increased, however, Fig. 3a and b shows that there is no proportional relationship between the number of GBS and ZVD cases, e.g., in the third week of 2016 were reported 3878 and 41 of ZVD and GBS cases, respectively, and in week #7, the number of reported GBS cases drop down to 14 cases, although the number of ZVD cases increased to 5695. Instances, like this one, evidence that the increasing of ZVD cases does not imply more GBS cases. Other times the more ZVD cases take place, the more GBS cases happen as well, e.g., in the second week of 2016 2888 and 22 of ZVD and GBS cases respectively were reported, and both increased in the third week until 3878 and 41 of ZVD and GBS cases, respectively.

Fig. 3
figure 3

a Prediction of the number of Guillain-Barré Syndrome (GBS) cases given the number of Zika virus disease cases (ZVD) using kernel ridge regression and Gaussian Kernel, Colombia, October 10th, 2015-June, 24th 2017. b Linear regression for predicting the number of GBS cases given the number of ZVD cases, Colombia, October 10th, 2015-June, 24th 2017

From the geometrical point of view, a straight line does not fit the data points (Fig. 3a and b). Figure 3a shows the estimation of the number of Guillain-Barre cases given the number of ZVD cases. This prediction has been estimated through the kernel ridge regression and the Gaussian Kernel. Its accuracy, measured in terms of Root Mean Square Error (RMSE) and the coefficient of determination are 1.77 and 80.7%, respectively. In Fig. 3a, both the red line and crosses correspond with the estimation (or prediction) while the circles correspond with the actual cases. The red crosses (and the line) are very near to the circles, even some crosses overlap circles. This means the prediction has a low error, i.e., the prediction model is very accurate.

On the other hand, Fig. 3b. shows the linear-regression-based estimation of the number of GBS cases given the number of ZVD cases. The linear regression accuracy in terms of RMSE and the coefficient of determination are 4.87 and 58.53%, respectively. The straight line and crosses shown in Fig. 3b are farer from the circles than in Fig. 3a. This means the linear model has underfitting (is less accurate).

During the cross validation, the RMSE achieved for the second order polynomial kernel, the linear kernel, the sigmoid kernel, and the Gaussian kernel are 9.15, 9.2, 10.7, and 7.2 respectively. Finally, the best setting for the Gaussian kernel is where the regularization parameter and radius are equal to 5 × 10–3 and 10–4, respectively, where its coefficient of determination and RMSE over the whole dataset are 0.95 and 0.96, respectively.

Discussion

This study provides the first general analysis of ZVD in Colombia in the period from October 2015 to June 2017. The paper shows that Zika virus disease has had a huge impact on public health for the numbers of people affected and the presentation of Guillain-Barre syndrome and microcephaly cases associated to Zika virus disease. This study updates the epidemiological analysis of the Zika virus disease situation, it describes the geographical distribution of Zika virus disease and shows the numerical relationship between Zika virus disease with microcephaly cases and Guillain-Barre syndrome.

According to the retrieved data Zika’s attack rate among the total population of the country was 222 ZVD cases per 100,000 populations. This attack rate is similar to that observed in other Latin American countries such as Brazil, Venezuela, Costa Rica and Salvador [20]. Among the pregnant population, higher attack rates compared with total population were observed, but this situation has already been described previously in Colombia [21] and may be associated with an enhanced surveillance in this population sub-group due to congenital malformations that produces the diseases, such as microcephaly and other congenital malformations of the central nervous system.

Similar to other countries (e.g., Brazil) [22, 23], the most affected population in Colombia were women (AR, 286 per 100,000 population) and people between 20 and 39 years old (AR, 305 to 375 per 100,000 population). Pregnant women (18.5%), and people older than 65 years (3.18%) are considered as groups at risk by INS [14], because their post-infection condition might be more complicated than other groups. On the other hand, the rest of the affected population includes the economically active population, and thus it might be costly for the nation [24]. We recommend to carry out another study including the disability adjusted life years (DALYs) needs to be calculated in the near future.

The increasing number of both GBS and microcephaly cases due to ZKD has been described [7, 8]. By 24th June 2017, 0.22% (AR, 222 ZVD cases per 100,000 population) of Colombian population has been affected by ZVD [16], including almost 19,963 cases of infected pregnant women. This is an alarming situation because congenital defects increase during ZVD transmission, as has been observed in Colombia, Brazil and other regions of America [21]. The number of cases confirmed in laboratories is lower (9.10%) than reported in other countries [3]. Most cases confirmed in laboratories correspond to those notified in pregnant women (64.94%). During pregnancy the infection with ZIKV is associated to microcephaly and other abnormalities of central nervous system [25].

In March 2016, five months after the Zika virus outbreak, microcephaly cases were 1156 increased, which means that cases of microcephaly may correspond to congenital 1157 infections of women infected by ZIKV. Currently (June 2017, epiweek 25) Colombia has 710 ZVD cases associated with microcephaly; 174 were confirmed by laboratory. The increase of microcephaly cases during ZIKV transmission represents new information for clinicians because this active surveillance just started to be carried out after the first confirmed cases in Brazil. Prior to 2016, In Colombia, surveillance of microcephaly cases and other abnormalities of central nervous system, included evaluation for pathogens of the STORCH complex (syphilis, toxoplasmosis, other agents, rubella, cytomegalovirus, and herpes) [26]. ZIKV antibody testing was not done in Colombia due to high cross-reactivity with other endemic flavivirus [26].

The cases of ZVD in Colombia are decreasing, nevertheless, it is possible that new cases of microcephaly associated to ZVD will be present in the next few months when currently infected pregnant women deliver their new babies. Currently, CDC and INS are collaborating to understand long-term effects of Zika virus infection during pregnancy in Colombia, this surveillance will track for 5 years 5.000 pregnant women, as well as their spouses and children [27]. On the other hand, Zika virus also affects the hearing, visual and cognitive capacities of new-borns [28] and in Colombia this level of impact remains unknown.

Currently with the emergence of Chikungunya, Dengue and Zika a new challenge of differential diagnosis must be faced in the country, however, this emergence will present a high cost for the public health system of Colombia [1]. Geographically, the most affected states correspond to the Caribbean, Andean and Orinoquía regions of Colombia.

The Zika virus seems to have spread from the Caribbean area towards the Andean region in Colombia. It is interesting to note that there are hotspots in Norte de Santander, reporting area that are close to the Venezuelan border, where commercial activity and border crossing is intense. We believe that the high impact of ZIKV in Colombia was influenced by people moving across borders, because it is possible that people moving throughout the borders increased the mobility of the virus. We should remember that Colombia after Honduras y Belize is the country with the third highest density of ZIKV-infected people [3]. With the exception of areas of the country at 2000 m above sea level (i.e. Bogota, DC), all the 32 Colombian states were affected, the spread of the virus affected small, medium and large provinces; some cases were located at or more than 2000 m above sea level, despite the absence of known vectors for ZIKV at high elevation, these high altitude infections may be explained by travel to endemic areas or sexual transmission.

Several studies report that Guillain-Barre is associated with several etiologies, including ZIKV [1]. In this study, we have adopted kernel ridge regression to predict the number of Guillain-Barre cases given the number of ZVD cases. This model is appropriate where functional dependence is not linear (i.e., it is not proportional), in fact, Fig. 3a and b, shows there is no proportional relationship between the number of GBS and ZVD cases. Instances, like this one, evidence that the increasing of ZVD cases does not imply more GBS cases. It is important to emphasize that the linear kernel correspond to classical linear regression model, and there are two kernel which correspond with the more accurate models than the linear one. These outcomes are coherent with the distribution of the data points which cannot be aligned in a straight line in Fig. 3, which also means that there is no proportional relationship between ZVD and GBS cases, as it was above mentioned. As a consequence, the Gaussian kernel is the best choice for the prediction because with this kernel, KRR reported the smallest RMSE during cross validation, i.e., this is the most accurate kernel for this domain. We have performed 5-fold cross validation for tuning parameters (i.e. the radius and regularization parameters for the Gaussian kernel) and choosing the most appropriated kernel between sigmoid, polynomial, and Gaussian kernel. The cross validation results reveal that the regression model with the latter kernel outperforms the others in terms of accuracy.

The above mentioned regression model accurately predicted the number of Guillain-Barre cases. Nevertheless, there is an over fitting risk due to small sample size compared with those available in other domains. However, we lack another independent dataset to carry out the test, so the best setting achieved through cross validation was tested over the only dataset available to us. So in this study we provide evidence which might confirm the link between both diseases, notwithstanding that there are several reports where GBS has multiple associated etiologies.

As this incidence of ZVD decreases, GBS cases are expected to decrease as well. However, because of multiple etiologies, we expect GBS to persist even when ZVD disappears. In fact, we showed that the functional relationship between both diseases is not proportional, so it is not accurate to state the more Zika cases appear, the more Guillain Barre cases shall take place. Instead, as it is shown in Fig. 3a and b, the more Zika cases happen, the Guillain Barre cases will either increase or decrease, otherwise the linear kernel (i.e., the classical linear regression model) would have outperformed other kernels. Nonetheless, the prediction model adopted in this work is able to capture the regular patterns accurately, in spite of there are other known and latent (hidden) factors which explain the relationship between both illnesses. The drawback of KRR is that this model does not allows us interpret the prediction, however, it is the most accurate in this context.

Conclusions

This work contributes with an update on epidemiological analysis of the situation of ZVD in Colombia. The study describes the geographical distribution of ZIKV infection in the country and shows the association with microcephaly cases and GBS. A kernel ridge regression may be important to predict the number of Guillain-Barre cases given the number of ZVD cases. This study provides evidence that confirms the relationship between both diseases. Zika virus may have some influence of the new presentation of GBS in Colombia. There is a public health concern about the future of microcephaly cases in Colombia. Data essential in ZIKA control such as disaggregated data by age-group and by altitude are currently missing.

References

  1. Mattar S, Tique V, Miranda J, Montes JE, Garzon D. Undifferentiated tropical febrile illness in Córdoba, Colombia Not everything is dengue. J Infect Public Health. 2017;10:507–12. http://0-dx-doi-org.brum.beds.ac.uk/10.1016/j.jiph.2016.09.014.

  2. Mattar S, Ojeda C, Arboleda J, Arrieta G, Bosch I, Botia I, et al. Case report: microcephaly associated with Zika virus infection. Colombia BMC Infectious Diseases. 2017;17:423. https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-017-2522-6.

  3. Pan American Health Organization / World Health Organization. Regional Zika Epidemiological Update (Americas) March 10, 2017. Washington, D.C. D.C. 20037, United States of America. 2017. http://www2.paho.org/hq/index.php?option=com_content&view=article&id=11599&Itemid=41691&lang=en. Accessed 23 Apr 2017.

  4. Pan American Health Organization / World Health Organization. Zika suspected and confirmed cases reported by countries and territories in the Americas Cumulative cases, 2015–2017. Updated as of 20 April 2017. Washington, D.C.: 2017. http://www2.paho.org/hq/index.php?option=com_docman&task=doc_view&Itemid=270&gid=39552&lang=es.

  5. Instituto Nacional de Salud (Colombia). Boletin Epidemiologico Semanal. Semana epidemiológica número 40 de 2015. Bogota D.C. 2017. http://www.ins.gov.co/boletin-epidemiologico/Boletn%20Epidemiolgico/2015%20Boletin%20epidemiologico%20semana%2040.pdf. Accessed 3 Mar 2017.

  6. Instituto Nacional de Salud (Colombia). Semana epidemiológica número 52 de 2016. Bogotá, D.C.: 2017. http://www.ins.gov.co/boletin-epidemiologico/Boletn%20Epidemiolgico/2016%20Bolet%C3%ADn%20epidemiol%C3%B3gico%20semana%2052%20-.pdf. Accessed 17 Feb 2017.

  7. Pacheco O, Beltrán M, Nelson CA, Valencia D, Tolosa N, Farr SL, Padilla AV, Tong VT, Cuevas EL, Espinosa-Bode A, et al. Zika virus disease in Colombia — preliminary report. N Engl J Med. 2016. doi: 10.1056/NEJMoa1604037.

  8. Parra B, Lizarazo J, Jiménez-Arango JA, Zea-Vera AF, González-Manrique G, Vargas J, Angarita JA, Zuñiga G, Lopez-Gonzalez R, Beltran CL, et al. Guillain–Barré syndrome associated with Zika virus infection in Colombia. N Engl J Med. 2016;375:1513–23.

    Article  PubMed  Google Scholar 

  9. Cuevas EL, Tong VT, Rozo N, Valencia D, Pacheco O, Gilboa SM, et al. Preliminary report of microcephaly potentially associated with Zika virus infection during pregnancy — Colombia, January–November 2016. MMWR Morb Mortal Wkly Rep. 2016;65:1409–13.

    Article  PubMed  Google Scholar 

  10. Camacho E, Paternina-Gomez M, Blanco PJ, Osorio JE, Aliota MT. Detection of autochthonous Zika virus transmission in Sincelejo, Colombia. Emerg Infect Dis. 2016;22:927–9.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Rojas DP, Dean NE, Yang Y, Kenah E, Quintero J, Tomasi S, Ramirez EL, Kelly Y, Castro C, Carrasquilla G, et al. The epidemiology and transmissibility of Zika virus in Girardot and San Andres island, Colombia, September 2015 to January 2016. Euro Surveill. 2016;21. http://0-dx-doi-org.brum.beds.ac.uk/10.2807/1560-7917.ES.2016.21.28.30283.

  12. Towers S, Brauer F, Castillo-Chavez C, Falconar AK, Mubayi A, Romero-Vivas CM. Estimate of the reproduction number of the 2015 Zika virus outbreak in Barranquilla, Colombia, and estimation of the relative role of sexual transmission. Epidemics. 2016;17:50–5.

    Article  PubMed  Google Scholar 

  13. Instituto Nacional de Salud (Colombia). Boletín Epidemiológico. Bogotá, D.C.: 2017. http://www.ins.gov.co/boletin-epidemiologico/Paginas/default.aspx. Accessed 12 Feb 2017.

  14. Ministerio de Salud y Protección Social - Instituto Nacional de Salud (Colombia). Circular conjunta externa N° 00000061 de 2015. Vigilancia de la fiebre por virus Zika (ZIKV) en su fase II epidémica y fortalecimiento de la prevencion de la fiebre por virus Zika en grupos de riesgo. Bogotá, D.C.: 2017. http://www.ins.gov.co/Noticias/ZIKA/Circular%20Conj%20061%202015%20Fiebre%20zika.pdf. Accessed 12 Feb 2017.

  15. Ministerio de Salud y Protección Social - Instituto Nacional de Salud (Colombia). Protocolo de vigilancia en salud pública: enfermedad por virus Zika. Bogotá, D.C; 2017. http://bvs.minsa.gob.pe/local/MINSA/3449.pdf. Accesed 14 Jul 2017

  16. Departamento Administrativo Nacional de Estadística (DANE, Colombia). Catálogo de Objetos Geográficos del Marco Geoestadístico Nacional. Dirección de Geoestadística, 2014. Bogotá, D.C.: 2017. http://www.dane.gov.co/. Accessed 19 Sept 2017.

  17. Departamento Administrativo Nacional de Estadística (DANE, Colombia). Estadísticas vitales nacimientos y defunciones. Bogotá, D.C.: 2017. https://www.dane.gov.co/index.php/estadisticas-por-tema/demografia-y-poblacion/nacimientos-y-defunciones. Accessed 18 Sept 2017.

  18. Departamento Administrativo Nacional de Estadística (DANE, Colombia). Catálogo de Objetos Geográficos del Marco Geoestadístico Nacional. Dirección de Geoestadística, 2014. Bogotá, D.C.: 2017. http://www.dane.gov.co/. Accessed 18 Sept 2017.

  19. Kloog I, Haim A, Portnov BA. Using kernel density function as an urban analysis tool: investigating the association between nightlight exposure and the incidence of breast cancer in Haifa, Israel. Comput Environ Urban Syst. 2009;33:55–63.

    Article  Google Scholar 

  20. Saunders C, Gammerman A, Vovk V. Ridge regression learning algorithm in dual variables. En (ICML-1998) Proceedings of the 15th International Conference on Machine Learning. Morgan Kaufmann. 1998:515–21.

  21. Pan American Health Organization / World Health Organization. Zika suspected and confirmed cases reported by countries and territories in the Americas Cumulative cases, 2015–2017. Updated as of 13 July 2017. Washington, D.C.: 2017. http://www.paho.org/hq/index.php?option=com_docman&task=doc_view&Itemid=270&gid=40858&lang=es. Accessed 14 Jul 2017.

  22. Rojas DP, Dean NE, Yang Y, Kenah E, Quintero J, Tomasi S, et al. The epidemiology and transmissibility of Zika virus in Girardot and San Andres island, Colombia, September 2015 to January 2016. Euro Surveill. 2016;21(28). doi: 10.2807/1560-7917.ES.2016.21.28.30283.

  23. Jimenez Corona ME, De la Garza Barroso AL, Rodriguez Martínez JC, et al. Clinical and Epidemiological Characterization of Laboratory-Confirmed Autochthonous Cases of Zika Virus Disease in Mexico. PLoS Currents. 2016;8:ecurrents.outbreaks.a2fe1b3d6d71e24ad2b5afe982824053. doi:https://0-doi-org.brum.beds.ac.uk/10.1371/currents.outbreaks.a2fe1b3d6d71e24ad2b5afe982824053.

  24. Oviedo-Pastrana M, Méndez N, Mattar S, Arrieta G, Gomezcaceres L. Epidemic outbreak of chikungunya in two neighboring towns in the Colombian Caribbean: a survival analysis. Archives of Public Health. 2017;75:1. doi: 10.1186/s13690-016-0169-1.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Rasmussen SA, Jamieson DJ, Honein MA, Petersen LR. Zika virus and birth defects — reviewing the evidence for causality. N Engl J Med. 2016;374:1981–7.

    Article  CAS  PubMed  Google Scholar 

  26. Ministerio de Salud y Protección Social (Colombia). Lineamientos para la detección y manejo clínico integral de anomalías congénitas en fetos expuestos al virus del Zika durante la gestación de pacientes en Colombia. Bogotá, D.C.: 2017. https://www.minsalud.gov.co/sites/rid/Lists/BibliotecaDigital/RIDE/VS/PP/ET/linea-deteccion-manejo-clinico-anomalia-congenitas-fotos-zika.pdf. Accessed 3 Mar 2017.

  27. Center for Disease Control and Prevention. CDC and the Instituto Nacional de Salud of Colombia collaborate to understand long-term effects of Zika virus infection during pregnancy. Atlanta.: 2017. https://www.cdc.gov/media/releases/2016/p0902-cdc-colombia-collaboration.html. Accessed 18 Sept 2017.

  28. Center for Disease Control and Prevention. Microcephaly & Other Birth Defects. Atlanta.: 2017. https://www.cdc.gov/zika/healtheffects/birth_defects.html. Accessed 18 Sept 2017.

Download references

Acknowledgements

The authors would like to thank to Dr. Nicholas Komar, Centers for Disease Control and Prevention, Division of Vector-Borne Disease and to Dr. Alvaro Carrascal, Public Health, American Cancer Society Eastern Division for their critique and reviews of this paper. To University of Cordoba, Investigation office, project sustainability of research groups FMV-01-16.

Funding

No funds were obtained to conduct this epidemiological research.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article.

Author information

Authors and Affiliations

Authors

Contributions

SM, MOP, NM, ICC and GA contributed with the conception and design, acquisition of data, analysis and interpretation of data; all authors joined in the preliminaries and final drafts of the paper. All the authors read and approved the final manuscript.

Corresponding author

Correspondence to Salim Mattar.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Méndez, N., Oviedo-Pastrana, M., Mattar, S. et al. Zika virus disease, microcephaly and Guillain-Barré syndrome in Colombia: epidemiological situation during 21 months of the Zika virus outbreak, 2015–2017. Arch Public Health 75, 65 (2017). https://0-doi-org.brum.beds.ac.uk/10.1186/s13690-017-0233-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13690-017-0233-5

Keywords