Heterogeneity of Trypanosoma cruzi infection rates in vectors and animal reservoirs in Colombia: a systematic review and meta-analysis

Background The heterogeneity of Trypanosoma cruzi infection rates among triatomines insects and animal reservoirs has been studied in independent studies, but little information has been systematised to allow pooled and comparative estimates. Unravelling the main patterns of this heterogeneity could contribute to a further understanding of T. cruzi transmission in Colombia. Methods A systematic search was conducted in PubMed, Medline, LILACS, Embase, Web of Knowledge, Google Scholar and secondary sources with no filters of language or time and until April 2018. Based on selection criteria, all relevant studies reporting T. cruzi infection rates in reservoirs or triatomines were chosen. For pooled analyses, a random effects model for binomial distribution was used. Heterogeneity among studies is reported as I2. Subgroup analyses included: taxonomic classification, ecotope and diagnostic methods. Publication bias and sensitivity analyses were performed. Results Overall, 39 studies reporting infection rates in Colombia were found (22 for potential reservoirs and 28 for triatomine insects) for a total sample of 22,838 potential animals and 11,307 triatomines evaluated for T. cruzi infection. We have found evidence of 38/71 different animal species as potential T. cruzi reservoirs and 14/18 species as triatomine vectors for T. cruzi. Among animals, the species with the highest pooled prevalence were opossum (Didelphis marsupialis) with 48.0% (95% CI: 26–71%; I2 = 88%, τ2 = 0.07, P < 0.01) and domestic dog (Canis lupus familiaris) with 22.0% (95% CI: 4–48%; I2 = 96%, τ2 = 0.01, P < 0.01). Among triatomines, the highest prevalence was found for Triatoma maculata in the peridomestic ecotope (68.0%, 95% CI: 62–74%; I2 = 0%, τ2 = 0, P < 0.0001), followed by Rhodnius prolixus (62.0%, 95% CI: 38–84%; I2 = 95%, τ2 = 0.05, P < 0.01) and Rhodnius pallescens (54.0%, 95% CI: 37–71%; I2 = 86%, τ2 = 0.035, P < 0.01) in the sylvatic ecotope. Conclusions To our knowledge, this is the first systematic and quantitative analyses of triatomine insects and potential animal reservoirs for T. cruzi infection in Colombia. The results highlight a marked heterogeneity between species and provide initial estimates of infection rates heterogeneity. Electronic supplementary material The online version of this article (10.1186/s13071-019-3541-5) contains supplementary material, which is available to authorized users.


Background
Chagas disease is caused by the protozoan Trypanosoma cruzi, a parasite largely distributed across 21 Latin American countries where it has been associated with increased morbidity and mortality and a high economic burden [1,2]. Transmission through triatomine vectors is the most frequent route, followed in frequency by congenital, blood transfusions, transplants and oral mechanisms [3].
Trypanosoma cruzi transmission is highly complex as it is the cause of a zoonotic parasitosis present in more than 150 animal reservoirs and 147 triatomine vectors [4,5]. Additionally, six different discreet typing units (DTUs) of T. cruzi have been described so far [4,6,7]. The role that different animals and triatomine species play in T. cruzi transmission to humans depends on the local conditions that facilitate contact with humans in the different type of cycles: sylvatic, peridomestic and domestic [8]. Importantly, domestic transmission from non-native triatomines is subject to control via insecticide spraying and housing improvement, and there has been a great success controlling these vectors across Latin American countries [9][10][11]. However, sylvatic transmission is much more complex and control measures are less clear [12]; continued population growth and deforestation across sylvatic areas increase the likelihood that humans establish closer contact with both animal reservoirs and triatomine vectors from the sylvatic cycles [8]. Furthermore, there is an increased report of outbreaks arising from oral transmission in both rural areas and major urban populations of Brazil, Colombia and Venezuela, also suggesting an increased contact between sylvatic and domestic ecotopes [13].
In Colombia, 26 triatomine species have been reported as potential vectors for T. cruzi. Among them, Rhodnius prolixus, Triatoma dimidiata, T. maculata and T. venosa are the most frequently found close to humans, within houses or in the surrounding area (peridomicile) [6,19]. However, a large number of sylvatic species have the capacity to occasionally invade homes and to lead to human infection. The importance of these species is indicated by their distribution, density, transmission efficiency, feeding behaviour and the prevalence of infection [20,21]. The main triatomine species in domiciliated habitats in Colombia is R. prolixus. Important efforts to interrupt domestic transmission from this species to humans have been made in the last decades and great progress has been made in some parts of the country [22]. However, Chagas outbreaks resulting from oral transmission have been continuously reported in recent decades, highlighting the importance of gaining a better understanding of the sylvatic cycle [16].
Documenting and quantifying general trends and heterogeneity in infection rates among the different reservoir hosts and triatomine vector species is paramount to understanding the role that these species play in T. cruzi transmission to humans. Here, we conducted a literature review and meta-analysis of the existent evidence for T. cruzi presence and infection rates in vectors and animal reservoirs in Colombia.

Methods
The systematic review and meta-analysis were carried out following the verification checklist PRISMA to ensure standardised methods (Additional file 1).

Search strategy
A search for observational studies of T. cruzi infection among animal reservoirs and triatomine vectors was conducted in PubMed, Medline, Index Medicus, LILACS, Embase and Web of Knowledge without language or time filter until 5th April 2018. Additionally, a secondary search through Google Scholar and Colombia University thesis repositories was conducted. For details of the search algorithm, see Additional file 2: Table S1.

Selection criteria and data extraction
Selection criteria included: (i) study type as cross-sectional studies that included triatomines or potential animal reservoirs; (ii) conduced in Colombian territory; (iii) reporting details of the number of tested individuals and the number of positive individuals; (iv) with details on diagnostic method, ecotope of capture; (v) with a minimum sample size of 10 subjects for pooled calculations.
An Excel spreadsheet was used to report characteristics of the studies, including reference, year of the study, place, ecotope, sample size, diagnostic method and taxonomic classification. For quantitative analyses, both triatomines and potential animal reservoirs were grouped by taxonomic characteristics (order, family, genus and species) as much as possible.
Diagnostic methods for potential reservoirs were classified as: (i) serological methods that do not account for current infection and have high cross-reactivity with other tripanosomatid species; (ii) non-serological methods, including both parasitological (blood spread smear, hemoculture and xenodiagnostic) and molecular methods (PCR), that allow establishing solid evidence as T. cruzi reservoirs; and (iii) mixed methods, which use a combination of the two previous methods but with no clear information to allow disaggregated data.

Quality assessment
A content quality assessment was conducted. The items evaluated included research question, sampling method, study period, diagnostic method, and potential bias. We allocated a score to each item in order to create an overall classification of high, moderate, or low quality [26] (see Additional file 2: Figure S1).

Statistical analysis Descriptive statistics
Descriptive statistics are presented as proportion of infection and 95% confidence intervals (95% CI). Frequency tables and maps of geographical distribution are presented.

Pooled analysis
A meta-analysis of binary outcome with binomial distribution was conducted in R v.3.2.2. specifically using the 'meta' library [27]. In order to stabilise variance a Freeman-Tukey double-arccosine transformation was used [27,28]. Given the fact that the study design used for the multiple species, the sample size, geographical characteristics and methods used could differ between studies, a random effects model was considered more appropriate [29], where tau-squared (τ 2 ) represents the variance between studies [27,29]. For assessing heterogeneity between studies, the I 2 -and Cochranʼs Q tests are presented. I 2 indicates the proportion of the total variation attributed to heterogeneity, i.e. an assessment of the degree of coherence between studies [I 2 = 100% × (Q − df )/Q]. The Cochranʼs Q is calculated as the sum of the squared deviations of each study's estimate that allows weighting each study's contribution [30].

Publication bias and sensitivity analysis
Publication bias was assessed when possible (for metaanalysis with at least five studies included), by using the Eggerʼs test, funnel plots and Trim fill technique evaluating the potential effect of non-included studies on the general estimates [31]. Sensitivity analysis was conducted by re-assessing general effect while sequentially omitting one study at a time [32] and also restricting the analyses to only high quality studies.

General findings
Overall, 960 potentially useful studies were found. After reading all abstracts and based on the selection criteria indicated above, 40 studies containing relevant information of the presence of T. cruzi infection in animals or triatomines were selected, including one study on culture parasite isolates [33] that was relevant for qualitative analyses but not for prevalence calculations. Therefore, for prevalence calculations we used 39 studies; of these, 22 included data for prevalence in potential animal reservoir and 28 for prevalence in triatomine vectors (see Fig. 1 for details). Table 1 shows a general description of the final 39 studies included in the quantitative analyses.
Geographical distribution of the samples used in the studies for potential animal reservoirs and triatomine vectors include 9 and 16 departments (first administrative level), respectively, and also the special territory of the Sierra Nevada de Santa Marta (SNSM). The departments with the highest number of studies were Casanare, Córdoba, Meta and Bolívar for reservoirs and Casanare, Bolívar, César, Guajira, Magdalena and SNSM for triatomines (see Fig. 2) (for details, see Additional file 2: Table S2.

Potential reservoirs
The number of subjects (potential animal reservoirs) analysed in the 22 selected studies was 22,838 of which 1585 (6.94%) were found infected with Trypanosoma spp. Eighteen of these studies reported the capture method; Sherman ® and Tomahawak ® were the most commonly used. All 22 studies reported the ecotope of collection and diagnostic method. The geographical distribution of the studies included classified by taxonomic order are shown in Fig. 3. Overall, 71 animal species of nine mammal orders were studied, out of which 38 were found infected with T. cruzi ( Table 2).
The geographical distribution of potential animal reservoirs and proportion of infection found by location by taxonomic order are presented in Fig. 4. Further details on the general characteristics of the studies and proportion of infection are shown in Additional file 2: Table S3.

Triatomine insects
The total sample size among studies included was 11,307 triatomines of which 1721 (15.22%) presented Trypanosoma spp. infection; of these, 1354 (11.97%) presented infection by T. cruzi. The capture protocols most commonly used were active search and live bait traps. Out of the 28 studies, 23 registered disaggregated data for ecotope of collection. Diagnostic methods used were mainly parasitological (direct microscopic examination of extruded faeces) followed by PCR and mice inoculation (see Additional file 2: Table S4) for further details.
In the selected papers, 18 triatomine species were studied of which four (Belminus ferroae, Eratyrus mucronatus, Psammolestes arthuri and Panstrongylus lignarius) were found with no evidence and 14 with solid evidence of T. cruzi infection (Table 3). Table 3 shows the geographical distribution of all triatomine species studied, except for Triatoma dispar, Rhodnius robustus and R. colombiensis as the studies reporting these species did not have geographical information. The point prevalence of Trypanosoma spp. in triatomine insects and the geographical distribution is presented in Fig. 5.

Pooled analyses Pooled prevalence of T. cruzi infection in animal reservoirs
Out of the nine orders of potential reservoirs reported in this study, seven presented sufficient sample size and additional information for calculating pooled estimates of the prevalence of T. cruzi infection. The aggregated estimates showed that the species with highest prevalence was Didelphis marsupialis, followed by Canis familiaris (Table 4). Further details of the disaggregated results are shown in Additional file 2: Table S5. Two orders, Artiodactyla and Lagomorpha, were excluded from the quantitative analyses due to small sample size. Within the order Artiodactyla (i.e. white-tailed deer, collared peccary and boar), a prevalence ranging between 0-5% was registered. For the order Lagomorpha, only one study was found for the European rabbit Oryctolagus cuniculus with a prevalence of 22% (2/9).

Table 1 General information on the 39 studies included for quantitative analysis of T. cruzi infection in potential animal reservoirs and triatomine vectors
Notes: Letters (a) and (b) refer to different studies by the same first author in a single year. The study by León et al. [33] was included as part of the qualitative analyses (Tables 2 and 3) but not as part of the meta-analyses as this is a study on parasite isolates only The two orders with the highest number of studies were Didelphimorphia and Carnivora, with 11 studies each. The results are presented in Fig. 6. Not all of these studies were included in meta-analyses due to sample size restrictions. Eleven studies reported Didelphimorphia. However, only five of them were included for pooled estimates due to sample size restrictions. The estimated prevalence for common opossum (Didelphis marsupialis) was 48.0% (95% CI: 26-71%; I 2 = 88%, τ 2 = 0.07, P < 0.01). The order Carnivora was the only one with reports for prevalence based on serological tests. The most relevant species, the domestic dog (Canis lupus familiaris), presented a high pooled prevalence of T. cruzi infection measured by serological tests (24.0%, 95% CI: 1-64%; I 2 = 98%, τ 2 = 0.17, P < 0.01). Only one study included the kinkajou (Potos flavus) and presented a point prevalence of T. cruzi of 44.4%. The order Chiroptera (bats) and the order Rodentia (rodents) were found with a pooled prevalence of 15.0% (95% CI: 5-29%; I 2 = 96%, τ 2 = 0.04, P < 0.01) and 6.0% (95% CI: 2-12%, I 2 = 77%, τ 2 = 0.01, P < 0.01), respectively.

Publication bias and sensitivity analysis
Eggerʼs test was used for cases with at least five studies included for either reservoirs or triatomines. Only one animal reservoir family (Canidae) was found with evidence of potential publication bias. After applying the Trimfill method, estimated prevalence changed from 14.0% (95% CI: 1-37%) to 32.0% (95% CI: 10-58%). Publication bias results are presented in Additional file 2: Table S7 and Figures S2-S4.
On the other hand, sensitivity analysis by repeating all pooled estimates omitting one study at a time did not

Quality assessment
Among the 39 studies included in the meta-analysis, 67% were classified as high-quality, 21% as moderatequality and 12% as low-quality. The main problems found in the low-quality studies included unclear information on the study period, lack of details on taxonomic classification, and not enough information on selection methods of the subjects. The details are presented in Additional file 2: Tables S8, S9 and Figure  S7.

Discussion
This systematic review obtained 39 studies that included 71 species of potential reservoirs and 18 species of triatomine insects investigated for the presence of T. cruzi infection in Colombia. Among them, we found solid evidence of infection in 38 reservoirs and 14 triatomine vector species. To the best of our knowledge, this is the first study that systematically identified this evidence and analysed the heterogeneity of infection rates of T. cruzi
Within the order Didelphimorphia, D. marsupialis (common opossum) presented the highest pooled prevalence, 48.0%. The epidemiological importance of this mammal is highlighted for several reasons including its high susceptibility to infection, ability to act as both reservoir and vector, its persistent high infection rates and its highly adaptive behaviour that allows it to live close to both domestic and sylvatic habitats and to transport infection between households [17,34]. Opossums are mammals with night habits and a clear ability to adapt to artificial ecosystems. In both rural and urban areas, these animals frequently visit houses in search for food, establishing close contact with humans [14,17,35]. These animals have been implicated in outbreaks resulting from oral transmission in different areas including in Colombia, suggesting they are able to contaminate human food sources via anal secretions [36][37][38]. The high pooled prevalence found for this species in our study and the already documented tolerance to TcI infection (the same DTU identified in human cases from oral transmission) confirms opossum as the most important species so far identified for transmission to humans in the sylvatic cycle of T. cruzi [4,39,40].
Our results for the order Carnivora show that Canis lupus familiaris presented the highest pooled prevalence (17.0%) in the domestic habitat, with a non-significant difference between serological (24.0%; 95% CI: 1-64%) and non-serologicalal tests (14.0%; 95% CI: 1-31%). Interestingly, diverging results have been found in other countries using parasitological methods [18]. Whereas in the Argentinian Chaco, Canis lupus familiaris has been reported as the most important animal in the domestic cycle of T. cruzi, with both a high prevalence and a high parasitemia [41,42], other studies in Brazil have found the domestic dog only with a high seroprevalence, but with a low parasitaemia [43]. It has been suggested that domestic dogs have a paramount role for linking sylvatic and domestic cycles and suggesting that this family could be used as an indicator of recent local transmission [39]. Another interesting contrast occurs with regards to the family Felidae. Species of this family have been reported infected in studies from Argentina and Brazil [44,45],  12:308 whereas in our systematic review in Colombia only one study for this family was identified, which reported no T. cruzi infection [46].
Another important finding from our review is the estimated pooled prevalence for the orders Chiroptera and Rodentia, 15.0% (95% CI: 5-29%) and 6.0% (95% CI: 2-12%), respectively. The order Chiroptera (to which different species of bats belong) has been reported commonly sharing habitats with triatomines. Their ability to travel long distances and their long lifespan that ranges between 15-50 years [47] make them a very important taxonomic order to keep T. cruzi transmission active [48,49]. Within the Chiroptera, the frugivorous Neotropical bats were the most commonly studied and reported with infected species [50,51]. Interestingly, it has been documented that T. cruzi infection in bats could occur via different routes that include oral and congenital transmission, the latter of which could explain the high T. cruzi prevalence even in non-insectivorous bats [18].
An additional discrete typing unit identified in bats, TcBat, has been found genetically related to TcI [4]. In recent years, TcBat was isolated from a human case, a five-year-old girl in Colombia [52]. The authors of this case study suggest vectorial or oral as potential infection sources [52]. The high prevalence of T. cruzi infection in the order Chiroptera found in our study confirms that species of this order are very relevant reservoirs in sylvatic T. cruzi transmission. Additionally, TcBat detected in humans also reveals the potential relevance of Chiroptera to T. cruzi human infection. TcBat is not usually investigated via molecular diagnostic methods so that the real magnitude of human infection by TcBat still remains unknown [4,52]. Determining the importance of animal reservoirs in the transmission of T. cruzi depends on several factors, such as their ability to be infectious in natural conditions and the concentration of parasites in the blood (parasitemia), which in turn is determined by the species of the animal reservoir, the genotype of T. cruzi, the nutritional status of the reservoir, the route of infection and the concomitant infection that ensures its diffusion and maintenance in nature [34]. With regard to rodents, the importance of these animals as T. cruzi reservoirs includes their ability to maintain the parasite for long periods of time [53]. Additionally, it has been suggested that rodents act as secondary reservoirs in both domestic and sylvatic habitats where they can play a role in the oral transmission route as part of the food chain to larger predators [18,54,55].
Rhodnius prolixus was the species with the highest number of studies, which allowed us to perform subgroup analysis which evidenced an important difference between sylvatic and domestic cycles, 62.0% (95% CI: 38-84%) and 39.0% (95% CI: 4-82%), respectively. While this triatomine species is the most important domiciliated species reported in Venezuela, Colombian and Central America, our results also highlight the importance of sylvatic colonies. Due to the highly anthropophilic behaviour of this species, its ability to become domiciliated and its high transmission efficiency [57], our results add evidence of the relevance of this species in Colombia, both in domestic and sylvatic cycles. Interestingly, the sylvatic presence of R. prolixus has been reported increasingly in both the native palm trees (Attalea butyracea) and introduced agro-industrial palms (Elaeis guineensis) across the plains of Colombia [51,58,59].
The triatomine species with the highest prevalence found in our study was T. maculata but only located in the northeast part of the country [60,61]. In the past, due to its ornithophilic behaviour, it was traditionally considered not very relevant for human transmission [62]. However, recent studies demonstrate their potential importance due to their high prevalence rates and colonisation index [5,60,61,[63][64][65].
Another triatomine species with a very high prevalence was R. pallescens with geographical predominance in the north-west part of the country (Caribbean region). This vector has been also associated with Attalea butyracea palms [63,, has been reported frequently in Costa Rica and Nicaragua [67] and it is considered the main vector in Panama [68,69]. Importantly, this species has been found commonly circulating in areas of oral outbreaks in Colombia [37]. Table 4 Results of point or pooled prevalence estimates (meta-analysis using random effects model) by serological vs non-serological methods for the different orders of potential Trypanosoma cruzi animal reservoirs studied in Colombia Notes: Non-serological: includes parasitological methods (blood smears hemoculture and xenodiagnostic) and molecular methods (PCR) Serological: includes ELISA, IFAT or rapid tests based on the whole parasite or recombined antigens; +: positive  [70]. Eratyrus mucronatus has been reported in two studies, but it remains unclear whether or not they were found with natural infection [46,58]. Psammolestes arthuri has been reported in two studies associated to nests of birds but with no natural infection [65,71]. Finally, P. lignarius has been reported as Fig. 6 Meta-analysis for pooled prevalence of T. cruzi infection using random effects model for the orders Carnivora (a) and Didelphimorphia (b) without disaggregation by diagnostic tests used. Numbers 1 and 2 in parentheses indicate the same study but a different geographical distribution infected by Trypanosoma spp. based on non-speciesspecific diagnostic method (microscopy) [71].
Finally, other triatomine species that include R. robustus, R. colombiensis and T. dispar and the animal species Caluromys lanatus have been found infected with T. cruzi in Colombia and are reported in a study of parasite isolates [33]. However, as that study does not correspond to a prevalence assessment, it was not possible to use those results for the quantitative analyses.

Public health implications of these results
Several types of dynamic models for Chagas disease have been developed [72]. Most of them have been restricted to domiciliated habitats and in some cases included sylvatic cycle with a single vector and single host. The results from our study show that the heterogeneity of species involved in the sylvatic cycle and the heterogeneity of infection rates among them are worth considering when modelling particularly the sylvatic cycle, as some pioneer studies have started to consider by using multihost models [73].
Several cultural aspects related to human contact with some of the animal species described in our study could be potentially relevant. Several animal species with clear evidence as T. cruzi reservoirs are considered common food sources in indigenous and rural communities. Factors associated with the hunting, sacrifice and manipulation of these animal species could represent a potential risk for accidental contamination with T. cruzi in humans. Of particular importance, several species from the order Cingulata (popularly known armadillos and cachicamos) and from the order Chiroptera (bats) are believed to have curative properties in some communities; for this reason their blood is consumed immediately after sacrifice as part of curative rituals, a practice described even from the pre-Columbian era [74,75].
The results from this study have several potential public health implications. First, to our knowledge, this study is a first baseline, accounting for the diversity of reservoirs and triatomine vectors participating in T. cruzi Fig. 7 Meta-analysis for pooled prevalence of T. cruzi infection using random effects model for Rhodnius prolixus by ecotope of collection, and without disaggregation by diagnostic tests used transmission in Colombia, which can help understanding the broader risk to which human populations are exposed to. Secondly, the evidence of T. cruzi infection in animal reservoirs from different orders with a lifespan shorter than humans could be considered as an indicator of active transmission in a specific location. Finally, understanding the cultural relations between humans and several wild species studied here could potentially help inform populations on the potential risks associated with certain cultural practices.

Study limitations
Several relevant limitations arise from this research. First, the lack of sufficient details from several original papers reduces the possibility of more disaggregated analyses for some species. Secondly, our study evidences a high heterogeneity of study designs which makes it more difficult to obtain better pooled estimates [30]. Although publication bias was only detected in one of the species, we are aware that, due to the small number of studies, the test used could be underpowered and the fact that we only found this test as significant for one of the families with the highest number of studies (Canidae) confirms this limitation [31,32]. Another relevant limitation is that no studies were excluded based on their quality assessment, also due to the very limited number of studies. However, when this sensitivity analysis was possible, we did not find difference between using all or high-quality-only Table 5 Results of point or pooled prevalence estimates (meta-analysis using random effects model) for triatomine species by specific and non-specific T. cruzi infection in Colombia Notes: T. cruzi: includes all T. cruzi-specific methods (culture, PCR and mice inoculation); T. spp.: includes non-species-specific diagnostic methods (direct observation based on microscopic examination of extruded faeces); +: positive studies. Although we are aware of this limitation, we are also aware of the limitations in conducting studies with animals and particularly in sylvatic areas. The reduced sample size of many of the eligible studies leads to very broad confidence intervals; although this is a limitation, it is also a more reliable measure of the uncertainly in the prevalence estimates from such heterogeneous types of studies.
Finally, an important limitation is related to diagnostic methods. An important recent study in triatomine vectors has identified that optical microscopy used for routine surveillance missed ~50-75% of infections and wrongly scored as infected ~7% of the bugs when only one slide per vector was used [75]. For animals, we found the use of various tests with high specificity for T. cruzi detection such as xenodiagnostic, haemoculture and coloured smears, and in the case of the family Canidae serological methods were used and in some cases with PCR confirmation [76,77]. For vectors, however, most of the studies were conducted by direct microscopic observation which, as discussed above, is not highly speciesspecific. In the more recent studies PCR has also been included which certainly improves the quality of the measurement. We acknowledge that the use of optic microscopy can lead primarily to an underestimate of the prevalence infection rates in vectors. However, due the lack of data on the operative performance (sensitivity and specificity) of some of the techniques used in the original studies which can vary between years and operators, unfortunately we could not adjust the prevalence estimates in our meta-analysis. We have conducted subgroup analysis by diagnostic techniques in order to address this issue in vectors (see Table 5). This analysis suggests that techniques such as PCR (that can also differentiate species) are significantly more sensitive and able to detect highest values of the prevalence.

Conclusions
To our knowledge, this is the first systematic evaluation and quantitative analysis of T. cruzi infection rates in triatomine insects and potential animal reservoirs in Colombia. The results also confirm solid evidence of the presence of T. cruzi infection in 38 out of 71 animal species and in 14 out of 18 triatomine species studied. Among animals, the species with the highest pooled prevalence were opossums (Didelphis marsupialis) and domestic dogs (Canis lupus familiaris). Among triatomines, the highest prevalence was found in T. maculata from peridomestic ecotopes followed by R. prolixus and R. pallescens, both from sylvatic ecotopes. Further research is needed to better understand the interaction between these species, humans and parasites.

Additional files
Additional file 1. PRISMA checklist.
Additional file 2. Text S1. Additional text on quality assessment. Table S1. Boolean algorithms for literature search as of 5th of April 2018. Table S2. Pooled prevalence estimates of T. cruzi infection in both triatomines and potential animal reservoirs by departments (firs administrative units) in Colombia. Table S3. Detailed data for the 18 studies chosen for potential T. cruzi reservoirs in Colombia. Table S4. Detailed data for the 28 studies chosen for Trypanosoma cruzi triatomine vectors in Colombia. Table S5. Results of point or pooled prevalence estimates (meta-analysis using random effects model) by detailed diagnostic methods for the different orders of potential Trypanosoma cruzi animal reservoirs studied in Colombia. Table S6. Results of point or pooled prevalence estimates (meta-analysis using random effects model) of T. cruzi infection by detailed diagnostic methods for the different orders of triatomines studied in Colombia. Table S7. Results of Eggerʼs test for publication bias assessment for estimates of T. cruzi infection rates in animal reservoirs and triatomine vectors in Colombia. Table S8. Results of point or pooled prevalence estimates (meta-analysis using random effects model) for the different orders of potential Trypanosoma cruzi reservoirs studied in Colombia, using only high-quality reports. Table S9. Results of point or pooled prevalence estimates (meta-analysis using random effects model) of T. cruzi infection for the different orders of triatomines studied in Colombia, using only high-quality reports. Figure S1. Content quality assessment of the 39 studies used in the meta-analysis. Figure S2. Funnel plot of the 95% CI pseudo limits for the estimates of T. cruzi prevalence in animal reservoirs in Colombia.