The ecology of the genome and the dynamics of the biological dark matter.

Transposable elements (TEs) are essential components of the eukaryotic genomes. While mostly deleterious, evidence is mounting that TEs provide the host with beneficial adaptations. How 'selfish' or 'parasitic' DNA persists until it helps species evolution is emerging as a major evolutionary puzzle, especially in asexual taxa where the lack of sex strongly impede the spread of TEs. Since occasional but unchecked TE proliferations would ultimately drive host lineages toward extinction, asexual genomes are typically predicted to be free of TEs, which contrasts with their persistence in asexual taxa. We designed innovative 'Eco-genomic' models that account for both host demography and within-host molecular mechanisms of transposition and silencing to analyze their impact on TE dynamics in asexual genome populations. We unraveled that the spread of TEs can be limited to a stable level by density-dependent purifying selection when TE copies are over-dispersed among lineages and the host demographic turn-over is fast. We also showed that TE silencing can protect host populations in two ways; by preventing TEs with weak effects to accumulate or by favoring the elimination of TEs with large effects. Our predictions may explain TE persistence in known asexual taxa that typically show fast demography and where TE copy number variation between lineages is expected. Such TE persistence in asexual taxa potentially has important implications for their evolvability and the preservation of sexual reproduction.

Epidemiological modeling of Trypanosoma cruzi: Low stercorarian transmission and failure of host adaptive immunity explain the frequency of mixed infections in humans.

People living in areas with active vector-borne transmission of Chagas disease have multiple contacts with its causative agent, Trypanosoma cruzi. Reinfections by T. cruzi are possible at least in animal models leading to lower or even hardly detectable parasitaemia. In humans, although reinfections are thought to have major public health implications by increasing the risk of chronic manifestations of the disease, there is little quantitative knowledge about their frequency and the timing of parasite re-inoculation in the course of the disease. Here, we implemented stochastic agent-based models i) to estimate the rate of re-inoculation in humans and ii) to assess how frequent are reinfections during the acute and chronic stages of the disease according to alternative hypotheses on the adaptive immune response following a primary infection. By using a hybrid genetic algorithm, the models were fitted to epidemiological data of Argentinean rural villages where mixed infections by different genotypes of T. cruzi reach 56% in humans. To explain this percentage, the best model predicted 0.032 (0.008-0.042) annual reinfections per individual with 98.4% of them occurring in the chronic phase. In addition, the parasite escapes to the adaptive immune response mounted after the primary infection in at least 20% of the events of re-inoculation. With these low annual rates, the risks of reinfection during the typically long chronic stage of the disease stand around 14% (4%-18%) and 60% (21%-70%) after 5 and 30 years, with most individuals being re-infected 1-3 times overall. These low rates are better explained by the weak efficiency of the stercorarian mode of transmission than a highly efficient adaptive immune response. Those estimates are of particular interest for vaccine development and for our understanding of the higher risk of chronic disease manifestations suffered by infected people living in endemic areas.

Invasion speeds of Triatoma dimidiata, vector of Chagas disease: An application of orthogonal polynomials method.

Demographic processes and spatial dispersal of Triatoma dimidiata, a triatomine species vector of Chagas disease, are modeled by integrodifference equations to estimate invasion capacity of this species under different ecological conditions. The application of the theory of orthogonal polynomials and the steepest descent method applied to these equations, allow a good approximation of the abundance of the adult female population and the invasion speed. We show that: (1) under the same mean conditions of demography and dispersal, periodic spatial dispersal results in an invasion speed 2.5 times larger than the invasion speed when spatial dispersal is continuous; (2) when the invasion speed of periodic spatial dispersal is correlated to adverse demographic conditions, it is 34.7% higher as compared to a periodic dispersal that is correlated to good demographic conditions. From our results we conclude, in terms of triatomine population control, that the invasive success of T. dimidiata may be most sensitive to the probability of transition from juvenile to adult stage. We discuss our main theoretical predictions in the light of observed data in different triatomines species found in the literature.

Intrusive versus domiciliated triatomines and the challenge of adapting vector control practices against Chagas disease.

Chagas disease prevention remains mostly based on triatomine vector control to reduce or eliminate house infestation with these bugs. The level of adaptation of triatomines to human housing is a key part of vector competence and needs to be precisely evaluated to allow for the design of effective vector control strategies. In this review, we examine how the domiciliation/intrusion level of different triatomine species/populations has been defined and measured and discuss how these concepts may be improved for a better understanding of their ecology and evolution, as well as for the design of more effective control strategies against a large variety of triatomine species. We suggest that a major limitation of current criteria for classifying triatomines into sylvatic, intrusive, domiciliary and domestic species is that these are essentially qualitative and do not rely on quantitative variables measuring population sustainability and fitness in their different habitats. However, such assessments may be derived from further analysis and modelling of field data. Such approaches can shed new light on the domiciliation process of triatomines and may represent a key tool for decision-making and the design of vector control interventions.

Assessing the risk of West Nile Virus seasonal outbreaks and its vector control in an urbanizing bird community: An integrative R0-modelling study in the city of Merida, Mexico.

Urbanization is a global trend associated with key socio-economic issues, one of them being to control the transmission of infectious diseases to a urban fraction of the world's population that shall reach 68% in 2050. While urban growth has been shown to favor mosquito species responsible for the transmission of the West Nile Virus (WNV), a major human arbovirosis, the effects of concomitant changes in the host bird communities remain hard to anticipate albeit essential to quantify disease risk and to plan control initiatives. We developed a R0 modelling of WNV transmission in a urban bird community to assess the risk of outbreak in Merida, one of the cities with the highest growth rate in Mexico. The model was parameterized using ecological and epidemiological data collected over the past 15-years on the local vector, Culex quinquefasciatus, and avian community. We identified a 3-weeks summer period during which the vector population strongly amplifies the WNV enzootic transmission and lead to a significant risk of outbreaks in humans. Extensive sensitivity analyses showed that urbanization induced changes in the bird community could lead to an up-to 6-fold increase in the duration of the risk period, while the daily risk could rise by 40%. Interestingly, the increase in Quiscalus mexicanus abundance had 4-5 times larger impact than any other change in the bird community. In such a context, annihilating the current and future risk of WNV outbreaks in Merida requires reducing the mosquito population by 13% and up to 56%, respectively. This study provides an integrative assessment of the current and future risks of WNV outbreak in the fast urbanizing city of Merida, and points toward the implementation of epidemiological monitoring combined with preemptive measures targeting both C. quinquefasciatus and Q. mexicanus populations, as they are expected to have synergistic effects.

Ecology, evolution, and epidemiology of zoonotic and vector-borne infectious diseases in French Guiana: Transdisciplinarity does matter to tackle new emerging threats.

French Guiana is a European ultraperipheric region located on the northern Atlantic coast of South America. It constitutes an important forested region for biological conservation in the Neotropics. Although very sparsely populated, with its inhabitants mainly concentrated on the Atlantic coastal strip and along the two main rivers, it is marked by the presence and development of old and new epidemic disease outbreaks, both research and health priorities. In this review paper, we synthetize 15 years of multidisciplinary and integrative research at the interface between wildlife, ecosystem modification, human activities and sociodemographic development, and human health. This study reveals a complex epidemiological landscape marked by important transitional changes, facilitated by increased interconnections between wildlife, land-use change and human occupation and activity, human and trade transportation, demography with substantial immigration, and identified vector and parasite pharmacological resistance. Among other French Guianese characteristics, we demonstrate herein the existence of more complex multi-host disease life cycles than previously described for several disease systems in Central and South America, which clearly indicates that today the greater promiscuity between wildlife and humans due to demographic and economic pressures may offer novel settings for microbes and their hosts to circulate and spread. French Guiana is a microcosm that crystallizes all the current global environmental, demographic and socioeconomic change conditions, which may favor the development of ancient and future infectious diseases.

Patterns of house infestation dynamics by non-domiciliated Triatoma dimidiata reveal a spatial gradient of infestation in rural villages and potential insect manipulation by Trypanosoma cruzi.

OBJECTIVE: Chagas disease is a major vector-borne parasitic disease in Latin America, primarily transmitted to humans by triatomine vectors. Non-domiciliated triatomine species such as Triatoma dimidiata in the Yucatan peninsula, Mexico, can transiently invade houses and are emerging as a major challenge to control Trypanosoma cruzi transmission to humans. We analyzed the spatio-temporal spreading of house infestation by T. dimidiata in four rural villages. METHODS: Triatomines were collected in four rural villages over a 2 years period, and the spatio-temporal patterns of infestation were analyzed. RESULTS: Triatomines were consistently more abundant at the periphery of villages than in centers, indicating a much higher risk of T. cruzi transmission at the periphery. Male T. dimidiata were found further in the center of the village, while females remained closer to the periphery, suggesting differential dispersal capabilities between sexes, although the timing of dispersal appeared identical. Surprisingly, infected females were consistently collected in houses much further from the surrounding bushes than non-infected females, while the distribution of males was unaffected by their T. cruzi infection status, suggesting an increased dispersal capability in infected females. CONCLUSION: The spatial structure of infestation should be taken into account for the prioritization of vector control activities within villages, and spatially targeted interventions may be explored. A potential vector manipulation by T. cruzi, observed for the first time in triatomines, may favor parasite transmission to new hosts.

Virus-host coexistence in phytoplankton through the genomic lens.

Virus-microbe interactions in the ocean are commonly described by "boom and bust" dynamics, whereby a numerically dominant microorganism is lysed and replaced by a virus-resistant one. Here, we isolated a microalga strain and its infective dsDNA virus whose dynamics are characterized instead by parallel growth of both the microalga and the virus. Experimental evolution of clonal lines revealed that this viral production originates from the lysis of a minority of virus-susceptible cells, which are regenerated from resistant cells. Whole-genome sequencing demonstrated that this resistant-susceptible switch involved a large deletion on one chromosome. Mathematical modeling explained how the switch maintains stable microalga-virus population dynamics consistent with their observed growth pattern. Comparative genomics confirmed an ancient origin of this "accordion" chromosome despite a lack of sequence conservation. Together, our results show how dynamic genomic rearrangements may account for a previously overlooked coexistence mechanism in microalgae-virus interactions.

Clues about the genetic basis of adaptation emerge from comparing the proteomes of two Ostreococcus ecotypes (Chlorophyta, Prasinophyceae).

We compared the proteomes of two picoplanktonic Ostreococcus unicellular green algal ecotypes to analyze the genetic basis of their adaptation with their ecological niches. We first investigated the function of the species-specific genes using Gene Ontology databases and similarity searches. Although most species-specific genes had no known function, we identified several species-specific functions involved in various cellular processes, which could be critical for environmental adaptations. Additionally, we investigated the rate of evolution of orthologous genes and its distribution across chromosomes. We show that faster evolving genes encode significantly more membrane or excreted proteins, consistent with the notion that selection acts on cell surface modifications that is driven by selection for resistance to viruses and grazers, keystone actors of phytoplankton evolution. The relationship between GC content and chromosome length also suggests that both strains have experienced recombination since their divergence and that lack of recombination on the two outlier chromosomes could explain part of their peculiar genomic features, including higher rates of evolution.

Unravelling cis-regulatory elements in the genome of the smallest photosynthetic eukaryote: phylogenetic footprinting in Ostreococcus.

We used a phylogenetic footprinting approach, adapted to high levels of divergence, to estimate the level of constraint in intergenic regions of the extremely gene dense Ostreococcus algae genomes (Chlorophyta, Prasinophyceae). We first benchmarked our method against the Saccharomyces sensu stricto genome data and found that the proportion of conserved non-coding sites was consistent with those obtained with methods using calibration by the neutral substitution rate. We then applied our method to the complete genomes of Ostreococcus tauri and O. lucimarinus, which are the most divergent species from the same genus sequenced so far. We found that 77% of intergenic regions in Ostreococcus still contain some phylogenetic footprints, as compared to 88% for Saccharomyces, corresponding to an average rate of constraint on intergenic region of 17% and 30%, respectively. A comparison with some known functional cis-regulatory elements enabled us to investigate whether some transcriptional regulatory pathways were conserved throughout the green lineage. Strikingly, the size of the phylogenetic footprints depends on gene orientation of neighboring genes, and appears to be genus-specific. In Ostreococcus, 5' intergenic regions contain four times more conserved sites than 3' intergenic regions, whereas in yeast a higher frequency of constrained sites in intergenic regions between genes on the same DNA strand suggests a higher frequency of bidirectional regulatory elements. The phylogenetic footprinting approach can be used despite high levels of divergence in the ultrasmall Ostreococcus algae, to decipher structure of constrained regulatory motifs, and identify putative regulatory pathways conserved within the green lineage.

Trypanosoma cruzi transmission dynamics in a synanthropic and domesticated host community.

Trypanosoma cruzi is the causative agent of Chagas disease, a Neglected Tropical Disease affecting 8 million people in the Americas. Triatomine hematophagous vectors feed on a high diversity of vertebrate species that can be reservoirs or dead-end hosts, such as avian species refractory to T. cruzi. To understand its transmission dynamics in synanthropic and domesticated species living within villages is essential to quantify disease risk and assess the potential of zooprophylaxis. We developed a SI model of T. cruzi transmission in a multi-host community where vector reproduction and parasite transmission depend on a triatomine blood-feeding rate accounting for vector host preferences and interference while feeding. The model was parameterized to describe T. cruzi transmission in villages of the Yucatan peninsula, Mexico, using the information about Triatoma dimidiata vectors and host populations accumulated over the past 15 years. Extensive analyses of the model showed that dogs are key reservoirs and contributors to human infection, as compared to synanthropic rodents and cats, while chickens or other domesticated avian hosts dilute T. cruzi transmission despite increasing vector abundance. In this context, reducing the number of dogs or increasing avian hosts abundance decreases incidence in humans by up to 56% and 39%, respectively, while combining such changes reduces incidence by 71%. Although such effects are only reached over >10-years periods, they represent important considerations to be included in the design of cost-effective Integrated Vector Management. The concomitant reduction in T. cruzi vector prevalence estimated by simulating these zooprophylactic interventions could indeed complement the removal of colonies from the peridomiciles or the use of insect screens that lower vector indoor abundance by ~60% and ~80%. These new findings reinforce the idea that education and community empowerment to reduce basic risk factors is a cornerstone to reach and sustain the key objective of interrupting Chagas disease intra-domiciliary transmission.

Demographic and dispersal constraints for domestic infestation by non-domicilated chagas disease vectors in the Yucatan Peninsula, Mexico.

Chagas disease is one of the most important diseases in Latin America. Insecticides have been sprayed to control domiciliated vectors. However, some triatomine species are not strictly domiciliated, and the transmission risk posed by immigrants is identified as a major challenge. The design of new control strategies requires disentangling the importance of demography and immigration in vector occurrence inside houses. Using a population dynamics model, we confirmed that dispersal can explain satisfactorily the domestic abundance of Triatoma dimidiata in Yucatan, Mexico. A surprisingly low fecundity was also required (no more than one to two female offspring per female per trimester). A wide range of survival probabilities was possible, although the best fit was obtained for a very low immature survival (< or = 0.01/trimester). Our model predicted that domestic populations are not sustainable, and up to 90% of the individuals found in houses are immigrants. We discuss the potential of different strategies to control the transmission of Chagas disease by non-domiciliated vectors.

Differential effects of caffeine on the antidepressant-like effect of amitriptyline in female rat subpopulations with low and high immobility in the forced swimming test.

The interaction of caffeine (1 mg/kg) and amitriptyline (15 mg/kg) on the immobility time (IT) during Porsolt's forced swimming test (FST) was investigated in female Wistar rats. Akaike's Information Criterion indicated that the ITs recorded from 142 rats during the first day of the FST followed a bimodal distribution. Hence, the median (125.5 s) was used to classify the animals in subpopulations with low (<125.5 s, LI-rats) or high (>125.5 s, HI-rats) immobility. The paired t-test was used to compare the change of ITs between the first and second swimming sessions. Vehicle-treated animals had a significant increase of ITs during the second day of the test, either in LI-rats (77+/-12 s vs. 196+/-8 s, P<0.0001, n=6) or HI-rats (150+/-8 s vs. 201+/-10 s, P<0.02, n=6). In LI-rats amitriptyline only prevented the increase of ITs during the second session (74+/-27 s vs. 97+/-42 s, n=12), whereas in HI-rats the antidepressant produced a significant decrease of ITs during the second session (161+/-22 s vs. 118+/-32 s, n=7, P<0.02). While caffeine alone prevented the increase of ITs in both groups, the methylxanthine abolished the effect of amitriptyline in HI-rats (165+/-23 s vs. 165+/-46 s, n=9), leaving the antidepressant action unaffected in LI-rats (87+/-23 s vs. 96+/-58 s, n=9). These results suggest that the anti-immobility effect of amitriptyline in HI-rats is mediated in part by endogenous adenosine.

To delay once or twice: the effect of hypobiosis and free-living stages on the stability of host-parasite interactions.

The life cycle of many endoparasites can be delayed by free-living infective stages and a developmental arrestment in the host referred to as hypobiosis. We investigated the effects of hypobiosis and its interaction with delay in the free-living stages on host-parasite population dynamics by expanding a previous attempt by Dobson & Hudson. When the parasite life cycle does not include free-living stages, hypobiosis destabilizes the host-parasite interactions, irrespective of the assumptions about the regulation of the host population dynamics. Interestingly, the destabilizing effect varies in a nonlinear way with the duration of hypobiosis, the maximal effect being expected for three to five months delay. When the parasite life cycle involves free-living stages, hypobiosis of short or intermediate duration increases the destabilizing effect of the first time delay. However, hypobiosis of a duration of five months or more can stabilize interactions, irrespective of the regulation of the host population dynamics. Overall, we confirmed that hypobiosis is an unusual time delay as it can stabilize a two-way interaction. Contrary to the previous conclusions, such an atypical effect does not require self-regulation of the host population, but instead depends on the existence of free-living stages.

Evolutionary ecology of Chagas disease; what do we know and what do we need?

The aetiological agent of Chagas disease, Trypanosoma cruzi, is a key human pathogen afflicting most populations of Latin America. This vectorborne parasite is transmitted by haematophageous triatomines, whose control by large-scale insecticide spraying has been the main strategy to limit the impact of the disease for over 25 years. While those international initiatives have been successful in highly endemic areas, this systematic approach is now challenged by the emergence of insecticide resistance and by its low efficacy in controlling species that are only partially adapted to human habitat. In this contribution, we review evidences that Chagas disease control shall now be entering a second stage that will rely on a better understanding of triatomines adaptive potential, which requires promoting microevolutionary studies and -omic approaches. Concomitantly, we show that our knowledge of the determinants of the evolution of T. cruzi high diversity and low virulence remains too limiting to design evolution-proof strategies, while such attributes may be part of the future of Chagas disease control after the 2020 WHO's target of regional elimination of intradomiciliary transmission has been reached. We should then aim at developing a theory of T. cruzi virulence evolution that we anticipate to provide an interesting enrichment of the general theory according to the specificities of transmission of this very generalist stercorarian trypanosome. We stress that many ecological data required to better understand selective pressures acting on vector and parasite populations are already available as they have been meticulously accumulated in the last century of field research. Although more specific information will surely be needed, an effective research strategy would be to integrate data into the conceptual and theoretical framework of evolutionary ecology and life-history evolution that provide the quantitative backgrounds necessary to understand and possibly anticipate adaptive responses to public health interventions.

Detailed ecological associations of triatomines revealed by metabarcoding and next-generation sequencing: implications for triatomine behavior and Trypanosoma cruzi transmission cycles.

Trypanosoma cruzi is the agent of Chagas disease, transmitted by hematophagous triatomine vectors. Establishing transmission cycles is key to understand the epidemiology of the disease, but integrative assessments of ecological interactions shaping parasite transmission are still limited. Current approaches also lack sensitivity to assess the full extent of this ecological diversity. Here we developed a metabarcoding approach based on next-generation sequencing to identify triatomine gut microbiome, vertebrate feeding hosts, and parasite diversity and their potential interactions. We detected a dynamic microbiome in Triatoma dimidiata, including 23 bacterial orders, which differed according to blood sources. Fourteen vertebrate species served as blood sources, corresponding to domestic, synantropic and sylvatic species, although four (human, dog, cow and mice) accounted for over 50% of blood sources. Importantly, bugs fed on multiple hosts, with up to 11 hosts identified per bug, indicating very frequent host-switching. A high clonal diversity of T. cruzi was detected, with up to 20 haplotypes per bug. This analysis provided much greater sensitivity to detect multiple blood meals and multiclonal infections with T. cruzi, which should be taken into account to develop transmission networks, and characterize the risk for human infection, eventually leading to a better control of disease transmission.

Non-randomized controlled trial of the long-term efficacy of an Ecohealth intervention against Chagas disease in Yucatan, Mexico.

Non-domiciliated intrusive triatomine vectors are responsible for a low but significant transmission of Trypanosoma cruzi to humans. Their control is a challenge as insecticide spraying is of limited usefulness, and alternative strategies need to be developed for a sustainable control. We performed a non-randomized controlled trial of an Ecohealth intervention based on window insect screens and community participation to reduce house infestation by Triatoma dimidiata in two rural villages in Yucatan, Mexico. Efficacy of the intervention was measured over a three years follow-up period and entomological indicators showed that the proportion of triatomines found inside houses was significantly reduced in houses with insect screens, which effectively kept more bugs on the outside of houses. Using a previously developed model linking entomological data to the prevalence of infection in human, we predicted that the intervention would lead to a 32% reduction in yearly incidence and in the prevalence of T. cruzi infection. The cost for the coverage of all the windows of a house was of comparable magnitude to what families currently spend on various domestic insecticide, and most screens were still in good conditions after three years. In conclusion, the Ecohealth approach proposed here is effective for the long-term and sustainable control of intrusive T. dimidiata vectors in the Yucatan peninsula, Mexico. This strategy may also be easily adapted to other intrusive triatomine species as well as other regions/countries with comparable eco-epidemiological settings, and would be an excellent component of a larger integrated program for the control of a variety of other vector-borne diseases, bringing additional benefits to the communities. Our results should encourage a further scaling-up of our implementation strategy in additional villages in the region.

Amazonian Triatomine Biodiversity and the Transmission of Chagas Disease in French Guiana: In Medio Stat Sanitas.

The effects of biodiversity on the transmission of infectious diseases now stand as a cornerstone of many public health policies. The upper Amazonia and Guyana shield are hot-spots of biodiversity that offer genuine opportunities to explore the relationship between the risk of transmission of Chagas disease and the diversity of its triatomine vectors. Over 730 triatomines were light-trapped in four geomorphological landscapes shaping French-Guiana, and we determined their taxonomic status and infection by Trypanosoma cruzi. We used a model selection approach to unravel the spatial and temporal variations in species abundance, diversity and infection. The vector community in French-Guiana is typically made of one key species (Panstrongylus geniculatus) that is more abundant than three secondary species combined (Rhodnius pictipes, Panstrongylus lignarius and Eratyrus mucronatus), and four other species that complete the assemblage. Although the overall abundance of adult triatomines does not vary across French-Guiana, their diversity increases along a coastal-inland gradient. These variations unravelled a non-monotonic relationship between vector biodiversity and the risk of transmission of Chagas disease, so that intermediate biodiversity levels are associated with the lowest risks. We also observed biannual variations in triatomine abundance, representing the first report of a biannual pattern in the risk of Chagas disease transmission. Those variations were highly and negatively correlated with the average monthly rainfall. We discuss the implications of these patterns for the transmission of T. cruzi by assemblages of triatomine species, and for the dual challenge of controlling Amazonian vector communities that are made of both highly diverse and mostly intrusive species.

Fundamental factors determining the nature of parasite aggregation in hosts.

The distribution of parasites in hosts is typically aggregated: a few hosts harbour many parasites, while the remainder of hosts are virtually parasite free. The origin of this almost universal pattern is central to our understanding of host-parasite interactions; it affects many facets of their ecology and evolution. Despite this, the standard statistical framework used to characterize parasite aggregation does not describe the processes generating such a pattern. In this work, we have developed a mathematical framework for the distribution of parasites in hosts, starting from a simple statistical description in terms of two fundamental processes: the exposure of hosts to parasites and the infection success of parasites. This description allows the level of aggregation of parasites in hosts to be related to the random variation in these two processes and to true host heterogeneity. We show that random variation can generate an aggregated distribution and that the common view, that encounters and success are two equivalent filters, applies to the average parasite burden under neutral assumptions but it does not apply to the variance of the parasite burden, and it is not true when heterogeneity between hosts is incorporated in the model. We find that aggregation decreases linearly with the number of encounters, but it depends non-linearly on parasite success. We also find additional terms in the variance of the parasite burden which contribute to the actual level of aggregation in specific biological systems. We have derived the formal expressions of these contributions, and these provide new opportunities to analyse empirical data and tackle the complexity of the origin of aggregation in various host-parasite associations.

Ecology, evolution and control of Chagas disease: a century of neglected modelling and a promising future.

More than 100 years after its formal description, Chagas disease remains a major public health concern in Latin America with a yearly burden of 430,000 Disability-Adjusted Life Years (DALYs). The aetiological agent, a protozoan named Trypanosoma cruzi, is mainly transmitted to mammalian hosts by triatomine vectors. Multiple species of mammals and triatomines can harbour and transmit the parasite, and the feeding range of triatomine species typically includes many noncompetent hosts. Furthermore, the transmission of the pathogen can occur via several routes including the typical vector's faeces, but also oral, congenital and blood transfusion routes. These ecological and epidemiological complexities of the disease have hindered many control initiatives. In such a context, mathematical models provide invaluable tools to explore and understand the dynamics of T. cruzi transmission, and to design, optimize and monitor the efficacy of control interventions. We intend here to provide the first review of the mathematical models of Chagas disease, focussing on how they have contributed to our understanding of (1) the population dynamics and control of triatomine vectors, and (2) the epidemiology of T. cruzi infections. We also aim at suggesting promising lines of modelling that could further improve our understanding of the ecology, evolution, and control of the disease.