Frédéric Lardeux
email: frederic.lardeux@ird.fr

Keywords: competition and coexistence amongst wild and domestic Bolivian populations of Chagas vector Triatoma infestans (Hemiptera, Reduviidae, Triatominae) Niche invasion, competition and coexistence amongst wild and domestic Bolivian populations of Chagas vector Triatoma infestans (Hemiptera, Reduviidae, Triatominae) Invasion de niche, compétition et coexistence entre populations Triatominae, invasion, competition, coexistence, domiciliation, Chagas Triatominae, invasion, compétition, coexistence, domiciliation, Chagas p.4

diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Chagas disease is an endemic parasitic anthropozoonosis from the American continent caused by Trypanosoma cruzi, a flagellate protozoan which can infect a wide range of wild and domestic mammals as well as humans [START_REF] Jansen | Domestic and wild mammalian reservoirs[END_REF]. Human infections can occur through blood transfusion, congenital transmission, oral route, organ transplant etc., but the most important route by far is through insect vectors (>80% of the cases) [START_REF] Schofield | Triatominae. Biology and control[END_REF]. Insect vectors are hematophagous bugs (Hemiptera, Reduviidae, Triatominae) [START_REF] Lent | Revision of the Triatominae (Hemiptera, Reduviidae), and their significance as vector of Chagas Disease, Bull[END_REF] which acquire the infection by sucking blood from an infected mammal and may retransmit the parasite during the next blood-meal, not through the bite itself but by defecating while sucking blood. The trypanosomes are in the faeces and therefore can penetrate the skin using any mucus membrane or any wound. In the infected host, the disease tends quickly towards a chronic presentation, the parasite causing not immediately apparent but irreversible lesions in the cardiac, digestive and nervous system, frequently leading to the death of the host.

Most species of Triatominae occupy predominantly wild habitats that offer shelter from climatic extremes along with easy access to a blood source: bird nests or rodent burrows are common examples. However, some species can also invade and colonise peridomestic habitats such as hen houses or goat corrals, and some have successfully colonised human dwellings. This tendency to colonize houses is important in terms of disease transmission to humans. Transmission is dependent on various factors and for triatomines the most significant are, among others, the blood trophic preference, the rapidity to which the insect can defecate shortly after blood-feeding and therefore leaving T. cruzi close to possible body entries [START_REF] Zeledon | Epidemiology, modes of transmission and reservoir hosts of Chagas disease[END_REF],

and even the density of the insects in a house [START_REF] Schofield | Density-dependent perception of triatominae bugs bites[END_REF]. However, the most important parameter by far is the tendency of the species to live close to human dwellings (i.e., its degree of domesticity) and therefore to enter frequently in contact with him [START_REF] Dujardin | Les vecteurs de la maladie de Chagas. Recherches taxonomiques, biologiques et génétiques[END_REF]. Indeed, major vectors of Chagas disease such as Triatoma infestans, Rhodnius prolixus, T. dimidiata or T. brasiliensis are amongst the species that have successfully colonized human dwellings in which they now accomplish their life cycle and develop colonies [START_REF] Dujardin | Les vecteurs de la maladie de Chagas. Recherches taxonomiques, biologiques et génétiques[END_REF]. The "domiciliation process" of Triatominae [START_REF] Schofield | Triatominae. Biology and control[END_REF][START_REF] Schofield | The process of domestication in Triatominae[END_REF] is of major importance in transmission dynamics.

Interrelationships between wild (i.e., "sylvatic"), peridomestic and domestic populations of Triatominae (from the same or from different species) should therefore be carefully studied to propose efficient vector control strategies [9 among others]. p. [START_REF] Zeledon | Epidemiology, modes of transmission and reservoir hosts of Chagas disease[END_REF] Unfortunately, the existence of interspecific interactions between triatomine species capable of domiciliation is not yet well documented. A particular species has two ways of being successful in colonizing habitats: it may be very good at colonizing open sites, or it may be very good at displacing other species from a given site. Between these two extremes, species may coexist if none of them have too superior a combination of both traits [START_REF] Hastings | Disturbance, coexistence, history and competition for space[END_REF]. For triatomines, colonizing open sites, and in particular human dwellings, may occur when humans enter new territories. Such a process is observed in Bolivia for R. stali [START_REF] Matias | Domiciliation process of Rhodnius stali (Hemiptera: Reduviidae) in Alto Beni, La Paz, Bolivia[END_REF], T. sordida, Eratyrus mucronatus [START_REF] Noireau | Sylvatic triatomines (Hemiptera: Reduviidae) in Bolivia : trends toward domesticity and possible infection with Trypanosoma cruzi (Kinetoplastida: Trypanosomatidae)[END_REF] or Panstrongylus rufotuberculatus [START_REF] Noireau | Apparent trend to domesticity observed in Panstrongylus rufotuberculatus Champion, 1899 (Hemiptera: Reduviidae), Bolivia Res[END_REF], and in Argentina for T. guasayana [START_REF] Dujardin | Les vecteurs de la maladie de Chagas. Recherches taxonomiques, biologiques et génétiques[END_REF]. Another scenario would be when a desinsectization campaign empties such an ecological niche. Then some new species may colonize the empty niche from the wild environment [START_REF] Gamboa | Dispersion de Rhodnius prolixus en Venezuela[END_REF][START_REF] Bar | Presencia de Triatoma sordida Stal (1859) en ecotopos urbanos de la ciudad de Corrientes, Argentina[END_REF] or the same species may re-colonize the niche from the peridomestic environment that is generally less well controlled by pesticides [16 -20]. Despite these examples of the struggle for domiciliation, there is no clear record of coexistence of colonies of two different species in the same house. Apart of very infrequent observations [START_REF] Espinola | Colonization and competition of Rhodnius prolixus, Rhodnius pictipes, and Triatoma maculata (Hemiptera: Reduviidae, Triatominae) in experimental chicken houses in Venezuela[END_REF], only one single species is generally present in one type of ecotope [22 among others]. It seems there exist competitive exclusion, which specifies that two different species cannot share the same ecological niche: one of the two competitors will always win over the other. The loser will either become extinct or will shift towards a different ecological niche by means of evolution or behaviour [START_REF] Harding | The competitive exclusion principle[END_REF]. This principle might also explain a third scenario of domiciliation which consists of species displacements, such as those observed with T. infestans which appear to be a strong competitor to other species [7, 16, 24 -26].

In Bolivia, T. infestans is well domiciliated and is the main vector of T. cruzi. Domiciliated populations are numerous and widely distributed. There also exist space-limited wild foci which have been identified not only in Bolivia [27 -31], but also in Argentina [32 -35],

Paraguay [START_REF] Velasquez | Aspectos de la enfermedad de Chagas en Paraguay[END_REF], Brazil [START_REF] Barretto | Estudos sôbre reservatorios e vetores silvestres do Trypanosoma cruzi. IV.Infecçao natural do Triatoma infestans encontrado em ecotopos silvestres por tripanossomo semelhante ao T. Cruzi[END_REF] and Chile [START_REF] Bacigalupo | First finding of Chagas disease vectors associated with wild bushes in the Metropolitan Region of Chile[END_REF][START_REF] Bacigalupo | Sylvatic foci of the Chagas disease vector Triatoma infestans in Chile: description of a new focus and challenges for control programs[END_REF]. In these countries, wild foci of T. infestans were observed in the Chaco region and the Andine cordillera which represent only a small part of the geographic distribution of the species in Latin America. Wild and domestic populations belong to the same species and probably wild populations have had domestic ones as ancestors [29, 40 -42]. However, unlike for example, for T. dimidiata [START_REF] Ramirez | Genetic structure of sylvatic, peridomestic and domestic populations of Triatoma dimidiata (Hemiptera: Reduviidae) from an endemic zone of Boyaca, Colombia[END_REF], the two entities do not mix. They behave as two different "species" without or with limited gene flow [29, 42, 44 -46] and can be morphologically distinguished [START_REF] Dujardin | Metric differences between sylvatic and domestic Triatoma infestans (Hemiptera: Reduviidae) in Bolivia[END_REF][START_REF] Dujardin | Genetic analysis of Triatoma infestans following insecticidal control interventions in central Bolivia[END_REF]. However, although re-invasion of the domestic environment after insecticide treatment seems to be the fact of residual domestic T.

p.6

infestans populations [47 -50 among others], the possibility that wild populations may reinvade empty human dwellings cannot be discarded [START_REF] Noireau | Can wild Triatoma infestans foci in Bolivia jeopardize Chagas disease control efforts?[END_REF]. This would be a complication in terms of vector control strategy which therefore would have to take into account the role of wild populations as a potential reservoir for human dwelling reinvasion.

Therefore, it is crucial to better understand (1) how different species (or different populations that do not interbreed) can compete for colonization, (2) if they can manage it, which one is likely to succeed (i.e., what are the key factors for a successful establishment of a species in human dwellings, in particular if another species is already there?) and (3) if two species can coexist in the same ecological niche.

Two hypotheses can be suggested: (1) the two species (or separated "groups") can really share the same ecological niche and in that case invasion, competition and coexistence are part of the natural ecological processes that the two groups may undergo, or (2) the two groups cannot share the same niche because they don't have exactly the same ecological needs. In that case, species cannot interact and there is no possibility of competition or coexistence.

The second hypothesis deals, for example, with characteristics such as trophic preferences, microclimatic preferences, behavioural characteristics or any other genetically driven characteristic that would impede one species from colonizing the niche of the other species. This hypothesis will be discussed in a subsequent article for wild and domestic populations of

T. infestans.

In the present paper, two initial hypothesis will therefore be taken into account: (1) as stated above, the two groups of T. infestans may share the same ecological niche, and (2) the two populations (wild and domestic) are supposed not to interbreed and can thus be considered as two different species. These hypotheses enable competition and therefore, a more general theoretical framework to explore the population dynamics of two competing species of triatomines colonizing human dwellings. In this framework, a new modelling approach for these disease vectors can be proposed: Triatomines can be assimilated to "macroparasites" and human dwellings to "hosts" in a similar manner as for human parasitic diseases modelling [START_REF] Anderson | Infectious Diseases of Humans: Dynamics and Control[END_REF]. Therefore the mathematical model developed for macroparasites colonizing human populations by Bottomley et al. [START_REF] Bottomley | Population biology of multispecies helminth infection: Competition and coexistence[END_REF] can be adapted to "triatomines colonizing houses". The model enables to point out the importance of some life traits of the triatomine species for successful invasion, competition or coexistence. The model results although general in the framework of Triatomine dynamics, are discussed in the light of the problem of wild/domestic p.7 populations of T. infestans in Bolivia and the possible establishment of wild populations in place of domestic ones.

Material and Methods

The model initial hypothesis is that the two populations, i.e. wild and domestic ones, can live in the same biotopes and do not interbreed. Both populations may colonize a domestic environment (inside houses) or a wild environment as described in [START_REF] Noireau | Wild Triatoma infestans, a potential threat that needs to be monitored[END_REF]: rupiculous or arboreal environments consisting in marsupial or rodent burrows, rocky outcrops, hollow trees etc. Model parameters appear in Table 1 and are explained thereafter.

Model for domiciliation of a single species

A colony of triatomines will be defined as "all the insects from the same mother living in one house". From all the triatomines living in an area, the population with the potential to colonize a house and start a colony are fertilized females. This population is of size F. These females may belong to a colony already established in a house and fertilized there, or may come from another house (by flight for example) and may be already fertilized or will be fertilized in this new house. Thereafter, these insects will be referred as "fertilized females" for simplicity. The number of houses that can be colonized is n. In one house i, there are C i colonies of triatomines of one single species. F and C i are random variables.

The rate of production of new females by a colony is r. If the probability that a new-born female develop to the status of "fertilized female" is P, then the colonies in the n houses produce fertilized females at a rate δ 

 n i i C 1
, where δ = r.P.

In one house, "fertilized females" in colonies are lost either through natural mortality μ F or through emigration to another house (intending to start a new colony). The rate at which fertilized females are lost through house colonization can be computed as

F n i i    1 .
where the product Γ i .F represents the contact rate of one house i with fertilized females.

Houses may harbour large colonies of T. infestans, but generally less than the apparent carrying capacity (i.e. number of adequate places to breed) because populations undergo p.8 strong density-dependence effects which limit the number of triatomines in one house. These density-dependence relationships are related to the availability of hosts for blood-feeding [6, 54 -58] and hiding places availability [START_REF] Cecere | Effects of refuge availability on the population dynamics of Triatoma infestans in central Argentina[END_REF]. Therefore, the model should take into account a kind of density-dependent establishment, and it can be assumed that the probability of a "fertilized female" surviving and starting a new colony in a house i is lowered by the number of yet established colonies C i . This probability can therefore be formulated as (1 -ηC i ), where η is the "per colony reduction in probability of starting a new colony". In the model, whenever (1 -ηC i ) becomes negative, it is set to zero.

For model completion, some hypothesis should be formulated on the dynamics of both triatomine colonies and houses. A colony can disappear at a per capita rate μ C and a house at a per capita rate μ H . However, it is assumed that the number of houses is maintained at size n (i.e., each time a house disappears; it is replaced by a new one that can be infected by triatomines). If a house disappears, the colonies in this house disappear as well.

The model can then be formulated in terms of increasing or decreasing number of fertilized

i i C F   
-The number of colonies in one house decreases by one at a rate μ C . C i -The number of colonies in a house is set to 0 at a rate μ H

Model when house exposure to infestation is identical

In a first approach it is assumed that each house has the same contact rate with "fertilized females", i.e., Γ i = Γ for all i. Following the analysis of [START_REF] Bottomley | Population biology of multispecies helminth infection: Competition and coexistence[END_REF], the random variable F (size of the population of "fertilized females") is replaced by its mean m F for simplification and it is therefore possible to derive differential equations for the mean number of "fertilized females" m F (t), the mean number of colonies m C (t) and its variance σ 2 C (t). It is reasonable to consider that the life duration of the "fertilized female" stage (which corresponds to the house-seeking p.9 phase) is shorter than the duration of its corresponding colony. Therefore dm F /dt can be set to zero and the model simplifies to the two following equations for the mean number of colonies and its variance (see [START_REF] Bottomley | Population biology of multispecies helminth infection: Competition and coexistence[END_REF] for details):

dm C /dt = δ'. m C (1 -η.m C ) -((μ C + μ H ) . m C (1) dσ 2 C /dt = δ'. m c (1-η.(2σ 2 c + m c )) + μ c m c + μ H .m 2 c -(2 μ c + μ H ). σ 2 c ( 2
)
Where δ' = δ. Γ.n / (μ L + Γ.n) represents the product of δ, the rate at which a colony produces new "fertilized females", and the probability Γ.n / (μ L + Γ.n) of survival of a "fertilized female" to start a colony.

In a mathematical sense m C (t) and σ 2 C (t)are limiting values of the mean number of colonies and variance when the number of houses tends to infinity. However, these terms also well describe large but finite populations (i.e., >100 houses) as observed in field situations.

Epidemiologists define the basic reproduction number (R 0 ) of a disease as the expected number of secondary infections arising from a single individual during his entire infectious period, in a population of susceptible individuals where there is no immunity and in the absence of interventions to control the infection [START_REF] Heesterbeek | A brief history of R 0 and a recipe for its calculation[END_REF]. The R 0 concept has also been widely used in ecology where it measures individual reproductive success under ideal conditions (i.e., mean number of offspring produced over the lifetime of an individual) [START_REF] Smith | Mathematical Demography: Selected Papers, Biomathematics[END_REF]. From these definitions, it emerges that when R 0 < 1 each individual (or infected individual if a disease is concerned) produces, on average, less than one new (infected) individual, and therefore the population (or the disease) will not grow. If R 0 > 1 the population will grow (or, from the disease's point of view, the "pathogen" is able to invade the susceptible population of hosts).

This threshold notion is one of the most important and useful characteristic of the R 0 concept. This concept can be applied to the development of triatomine colonies (= the "pathogen") in the population of houses (= the "hosts"). In the present triatomine model, the basic reproduction number for the mean number of colonies can be computed from equation (1) as:

R 0 = δ'/ (μ C + μ H ) = (δ. Γ.n (μ L + Γ.n)) / (μ C + μ H ) (3)
p.10

When R 0 < 1, the triatomine species is unable to establish itself in the houses and therefore, m C = 0. This value is a stable point. When R 0 > 1, the point m C = 0 is unstable and the number of colonies can increase, up to an equilibrium point where m C = (1/η) (1-R 0 -1 ). The mean number of triatomine colonies is therefore dependent on the basic reproduction number R 0 and on the strength of density-dependence η (i.e., the "per colony reduction in probability of starting a new colony when other colonies are already established"). At this equilibrium, assuming that the rate of disappearance of houses μ H is much smaller than the other rates, the variance to mean ratio becomes R 0 -1 . Indeed, when equating equations ( 1) and ( 2) to zero at equilibrium, and equating

R 0 = δ'/ (μ C + μ H ) ≈. δ'/ μ C , straightforward calculations give σ 2 C / m C = R 0 -1 .
Then, if R 0 > 1 (which enables triatomine colonies to grow in number), then σ 2 C / m C < 1 which implies that the distribution of the number of triatomine colonies across the houses is distributed more evenly than at random. However, field observations prove that triatomine bugs are generally aggregated amongst houses, and therefore some houses harbour more colonies than others [START_REF] Schofield | Population dynamics and control of Triatoma infestans[END_REF].

Model when there is heterogeneity in house exposure to infestation

The model should therefore be modified to take into account the aggregated distribution of triatomines and therefore to enable σ 2 C / m C > 1. Several mechanisms can generate such heterogeneity. Amongst them, ecological factors that make some habitats less suitable than others. For example triatomine establishment can be limited if triatomine refuge quantity and quality diminish [START_REF] Schofield | Population dynamics and control of Triatoma infestans[END_REF] or when house walls are improved (this can be assimilated to a kind of "immunity" of the houses). Another likely mechanism is heterogeneous exposure of houses to colonization which can also mimic the above mechanisms. Among various factors that can influence heterogeneous exposure, one can think for example to close favourable infested peridomestic habitats that can generate higher domestic exposure in the vicinity [START_REF] Cecere | Spatio-temporal analysis of reinfestation by Triatoma infestans (Hemiptera: Reduviidae) following insecticide spraying in a rural community in northern Argentina[END_REF][START_REF] Cecere | The role of the peridomiciliary area in the elimination of Triatoma infestans from rural Argentina communities[END_REF]. In the model, heterogeneity can be achieved by enabling a random contact rate of houses with "fertilized females", i.e. replacing in the model the contact rate Γ which was fixed, by independent and identically distributed random variables Γ i (i= 1 to n) with mean m Γ and variance σ 2 Γ (and therefore a coefficient of variation of the exposure rate of

v Γ = σ Γ 2 / m Γ p.11
which measures the degree of heterogeneity). In that case, model development gives the following equation for the mean number of colonies at stable equilibrium:

) 4 ) 1 ( 2 1 ( ) 1 ( 2 1 2 0 2 2 2 1 0 2 2              R R m C      , if v Γ > 1 (4a) and ) 1 ( 2 1 1 0    R m C  , if v Γ =1 (4b)
Where

R 0 = δ'/ (μ C + μ H ) and δ' = n.δ.m Γ / (μ F + n.m Γ )
Again, when R 0 < 1 the triatomine species cannot invade the environment and when R 0 > 1, colonies tend to a non-zero equilibrium given by equation ( 4a) and (4b).

These equations indicate that at equilibrium, the mean number of triatomine colonies is, as before, dependent on (1) the basic reproduction number R 0 of the colonies, (2) the strength of density-dependence η, and now also on (3) the degree of heterogeneity v Γ through which the aggregative distribution of triatomine colonies within houses can be mimicked.

More precisely, equation (4a) indicates that the mean number of colonies increases as the degree of heterogeneity σ Γ decreases. In fact no house can harbour more than 1/η colonies and even houses with very high rates of exposure cannot have a corresponding high number of triatomine colonies. On the contrary, houses with rates of exposure close to zero will harbour a very low number of colonies. Thus, the mean number of colonies decreases as heterogeneity in exposure increases. The maximum number of colonies

m Cmax = (1/η) (1-R 0 -1 )
is reached when σ Γ =0, and m C ≈0 when σ Γ →∞.

Model for two competing species or populations

Two species (or two independent populations as it seems to be the case for wild and domestic populations of T. infestans) can compete for domiciliation only if they are both able to colonize the same ecotope. Clearly, there exist species of triatomines that cannot invade p.12 houses because they cannot encounter their ecological preferences in human dwellings. On the contrary, there may exist situations where continuous exchanges between wild and domestic populations take place, as it is the case, for example with Rhodnius prolixus [START_REF] Fitzpatrick | Molecular genetics reveal that sylvatic Rhodnius prolixus do colonize rural houses[END_REF]. In the present model it is assumed that the two "species" can evenly colonize both ecotopes.

Therefore, if the two species can mix, they may compete for resources such as blood sources [START_REF] Schofield | Density-dependent perception of triatominae bugs bites[END_REF] or even space [START_REF] Cecere | Spatio-temporal analysis of reinfestation by Triatoma infestans (Hemiptera: Reduviidae) following insecticide spraying in a rural community in northern Argentina[END_REF][START_REF] Gomez-Núñez | Desarrollo de un Nuevo método para evaluar la infestación intradomiciliaria por Rhodnius prolixus[END_REF]. The latter can be indirectly the consequence of house improvement for example, which limit shelters for T. infestans populations.

The model for one single species can be extended to two competing species considering that the probability that a "fertilized female" starts a colony in one house depends, as before, on the number of colonies of its own species (as in the single species model), but also now on the number of colonies of the other species already established. Therefore, one house i may then be colonized at time t by C 1i colonies of species 1 and C 2i colonies of species 2 respectively.

The "one species model" can then be modified by replacing δ, C i , Γ i , μ F , F, μ C in the transition rates of the Markov process by δ 1 , C 1i , Γ 1i , μ F1 , F 1 , μ C1 for species 1 and δ 2 , C 2i , Γ 2i , μ F2 , F 2 , μ C2 for species 2. Only the transition rate for increasing the number of colonies by one is modified to take into account species competition, modifying the parameter  for each species (per colony reduction in probability of stating a new colony when other colonies are already established), by the parameter  jk (j, k = 1, 2) which represents the effect of species j on the probability of establishment of species k.

Therefore, the Markov process for species 1 (and similarly for species 2) is:

-The number of fertilized females (species 1) increases by one at a rate  

n i i C 1 1 1



-The number of fertilized females (species 1) decreases by one at a rate

    n i i F F 1 1 1 1
). ( -The number of colonies of species 1 in one house increase by one at a rate ) . . 1 .( .

2 21 1 11 1 1 i i i C C F     
. Whenever this rate becomes <0, it is set to zero. The term

) . . 1 ( 2 21 1 11 i i C C    
represents the probability of establishment of an "infective female" of species 1 in a house where C 1 colonies of species 1 and C 2 colonies of species 2 are already present.

p.13 -The number of colonies (species 1) in one house decreases by one at a rate μ C1 . C1 i -The number of colonies in a house is set to 0 at a rate μ H

Model when house exposure to infestation by each species is identical

Following the same approach as for the one species model, and first considering that exposure to each species is the same for all houses (i.e., Γ 1i = Γ 1 and Γ 2i = Γ 2 for i= 1 to n) and setting dm F1 /dt = dm F2 /dt = 0 (i.e., life expectancy of "fertilized females" is much smaller than that of a colony), the mean number of colonies for species 1 and species 2 can be modelled as:

dm C1 /dt = δ 1 '. m C1 (1 -η 11 .m C1 -η 21 .m C2 ) -((μ C1 + μ H ) . m C1 (5) dm C2 /dt = δ 2 '. m C2 (1 -η 22 .m C2 -η 12 .m C1 ) -((μ C2 + μ H ) . m C2 (6)
Where δ' i = δ i . Γ i .n / (μ Fi + n.Γ i ) i = 1, 2 is the product of the rate at which potentially "infective females" are produced by a colony and the probability that an "infective female" survives to colonize a house.

Theses equations are symmetrical in regard to species 1 or 2, and are generalizations of equation [START_REF] Bottomley | Population biology of multispecies helminth infection: Competition and coexistence[END_REF].

Once again, the behaviour of the system is dependent on the basic reproduction numbers of the two species colonies and the inter-and intra-specific interactions parameters (i.e., the η ij ).

For colonies of species i, the basic reproduction number is R oi = δ' i / (μ Ci +μ H ). As in the "one species model", the mean number of colonies of species i will grow if the reproductive number R oi > 1. If reproductive numbers for both species are >1 then two scenarios are possible: (1) the exclusion of one species by the other and therefore the system approaches a single species equilibrium, or (2) coexistence of both species in a mixed equilibrium.

As a first approach, the probability of establishment of both species can be considered as dependent of the availability of resources in the house (blood access and space). In other words, the effect of species i (i= 1, 2) on the establishment of species i or species j is the same. Therefore in the model,  11 =  12 and  22 =  21. . With this assumption, competitive exclusion and coexistence can be understood by analysis of equations ( 5) and ( 6) in the plane p.14 (m C1 , m C2 ) (Fig. 1). In this plane, points P1, P2, P3 and P4 are computed as: In the phase plane analysis, coexistence is achieved when P 1 > P 2 and P 4 > P 3 .

P 1 = (1/η 21 )(1- R 01 -

Model when there is heterogeneity in house exposure to infestation

However, as for the one species model, heterogeneity in house exposure to infestation should be taken into account. Following [START_REF] Bottomley | Population biology of multispecies helminth infection: Competition and coexistence[END_REF], the model can be slightly modified and therefore, the variation with time of the mean number of colonies m ci (i=1,2) is:

1 1 2 1 2 1 21 1 1 1 1 11 1 1 1 . ' 1 . ' / C c C c C c C C m m m m m m dt dm                                         (7) p.15 1 1 1 1 2 1 2 2 1 21 1 1 1 2 1 1 11 1 1 1 2 1 1 1 1 . ' ) ( ) ( . ' / C C C C c C c C m m m m m m dt d                                     (8) 2 1 2 2 1 2 1 1 1 12 2 2 2 1 2 1 22 1 2 2 1 2 2 2 1 . ' ) ( ) ( . ' / C C C C c C c C m m m m m m dt d                                        (9) With     1 1 . / ' 1 1 1    n m n m F     and H Ci Ci      '
.(i= 1, 2). Similar equations can be derived for m c2 , σ Γ2C2 and σ Γ2C1 of species 2.

As in the single species model, the equation ( 7) describing the colony burden of species 1 in one house (and similarly for species 2) now also depends on the covariance between exposure to house infestation and the burden of colonies in houses (σ ΓiCj ).

Coexistence appears when mutual invadability is possible, i.e., when each species can invade an equilibrium where only the other species is present [START_REF] Hutson | Permanence and the dynamics of biological systems[END_REF]. Conditions of coexistence can therefore be derived from the model as follow: considering equilibrium e 1 and e 2 of each species, mutual invadability implies that e 1 and e 2 must be unstable to permit invasion of one species by the other. At e2, if few « infective females » of species 1 start (few) colonies and as such disturb the equilibrium, invasion will be successful if after some time the mean number of colonies of species 1 (m C1 ) is growing. In a mathematical sense, dm C1 /dt (equation ( 7)) >0. At the beginning of the process, the number of colonies of the invading species 1 is small as compared to species 2 and therefore, from equation [START_REF] Dujardin | Les vecteurs de la maladie de Chagas. Recherches taxonomiques, biologiques et génétiques[END_REF], the rate of increase of species 1 is:

1 2 1 2 1 21 1 1 ' 1 ' c c c m m r                          
When r 1 >0 the equilibrium 2 is unstable, enabling the establishment of species 1. In the same way r 2 determines the stability of equilibrium 1. The effective reproductive numbers of species 1 and 2 at equilibrium eq2 and eq1 respectively are:

R 1 = R 01 .r 1 and R 2 = R 02 .r 2 .
where [START_REF] Jansen | Domestic and wild mammalian reservoirs[END_REF]. Therefore, the condition for coexistence is R 1 >1 and R 2 >1. The two effective reproductive numbers R 1 and R 2 can be reformulated taking into account that at e i (i= 1, 2), the maximum possible value of the mean number of colonies for species i is attained when v Гi (the coefficient of variation for the rate of exposure of species i) is zero.

R 0i = (δ' i / μ' c ) (i=1,

p.16

Then, this maximum value is

mc imax = (1/η ii ) (1-R 0i -1 ). Now, introducing ρ = σ Γ1 Γ2 / σ Γ1. σ Γ2
(the correlation between house susceptibility to species 1 and house susceptibility to species 2) and developing r 1 and r 2 , the effective reproductive numbers R 1 and R 2 at equilibrium eq2 and eq1 respectively can be expressed as:

                        2 2 max 2 2 1 21 01 1 1 c c c m m m R R     (11)                         1 1 max 1 1 2 12 02 2 1 c c c m m m R R     (12)
These effective reproductive numbers appear to be the product of the basic reproductive number (R 0 ) and a term representing the probability of establishment. This probability depends on the burden of resident triatomines in a house (m cmax -m c ) and also on the covariance between the burden of the resident species and the susceptibility of the house to the invading species.

These equations indicate that at each respective equilibrium, the effective reproductive numbers are a decreasing function of ρ (the correlation between house susceptibility to species 1 and house susceptibility to species 2). Therefore, intuitively, decreasing the correlation between exposure rates will encourage coexistence.

When heterogeneity is the same for both species (i.e., v 1i = v 2 = v Г ), then the R i are increasing functions of v Г for ρ <1 and therefore, increasing the heterogeneity in house exposure to both species of triatome will facilitate coexistence. From these equations, it may also be the case that if v 1i >> v 2 (i.e., heterogeneity in house exposure to species 1 is high compared to heterogeneity in house exposure to the other species), then R 1 will be small because ν Г1 / ν Г2 >>1 and it will be very difficult for species 1 to establish itself. Therefore, coexistence is unlikely to succeed when the degree of heterogeneity in house exposure greatly differs between the two species of triatomine.

Results and discussion

p.17

A model has been developed to simulate the invasion and colonization of houses by triatomines. The model is a general framework which can describe all situations where one species of triatomine attempts to invade a particular ecotope, or where two species compete to colonize a particular ecotope. The model has been developed in the framework of T. infestans dynamics where individuals may colonize houses, or where two groups (domestic vs wild) may compete or coexist. In that case it was assumed that the two groups behaved as two different species, with no interbreeding.

Single species invasion and mean number of triatomine colonies per house

Importance of R 0 of triatomine colonies

The model captures the observation that the mean number of triatomine colonies in a house is dependent on the basic reproduction number R 0 of colonies, the strength of densitydependence processes (which govern, for example, access to blood sources and/or shelter for triatomes), and the degree of heterogeneity (i.e., the degree of exposure of houses to infestation by triatomines), through which the aggregative distribution of triatomine colonies within houses can be mimicked. Introducing mechanisms that can generate aggregation is of crucial importance in T. infestans dynamics because field observations indicate clearly the existence of such phenomenon [START_REF] Schofield | Triatominae. Biology and control[END_REF][START_REF] Schofield | Density-dependent perception of triatominae bugs bites[END_REF][START_REF] Schofield | The role of blood intake in density regulation of populations of Triatoma infestans (Klug) (Hemiptera: Reduviidae)[END_REF][START_REF] Smith | Mathematical Demography: Selected Papers, Biomathematics[END_REF][START_REF] Schofield | Density regulation of domestic populations of Triatoma infestans in Brazil[END_REF].

Density dependence processes and heterogeneity of colonization exposure

The model reveals that if density dependence processes are strong, it will be difficult for a new colony to start. Triatomine dynamics are dependent on the R 0 value which is (1) directly proportionate to the rate at which a colony produces new « fertilized females » and the probability of survival of « fertilized females » of starting a colony, and (2) inversely proportionate to the per capita rate of disappearance of a colony. Thus, demographic parameters are of importance, in particular, as one might expect intuitively, fecundity and survival [START_REF] Castañera | A stage-structured stochastic model of the population dynamics of Triatoma infestans, the main vector of Chagas disease[END_REF][START_REF] Oscherov | Interactions between Triatoma infestans and Triatoma sordida (Hemiptera: Reduviidae) in artificial ecotopes: population growth and age structure[END_REF]. However, these parameters are offset by density dependence processes. In the model, the mean number of colonies decreases as heterogeneity in exposure increases.

Heterogeneity of exposure can be understood in terms of an infested peridomicile or an p.18 infested nearby house or structure which can be a source of infestation for a house. Indeed, as heterogeneity increases, there are more houses with very high and very low rates of exposure.

However, since the number of colonies in a house is limited (as assumed in the model), houses with very high rates of exposure cannot have a corresponding high number of colonies. Thus, in total, the mean number of colonies decreases as heterogeneity increases. In the field, T. infestans exhibit such an aggregative distribution with highly infested houses and low infested ones, but always with a limited maximum number of colonies per house.

Moreover, when the peridomestic environment is infested (i.e., high exposure rate for the house), and if the "immunity" of the house is low (i.e. if triatomines can find shelter and accessible blood sources inside) [START_REF] Schofield | Population dynamics and control of Triatoma infestans[END_REF], the house also hosts domestic populations of triatomines.

A partial concluding remark would be that for a species to invade and colonize an ecotope, the model points out the predominant role of intuitive parameters such as the survival of "infective females", the survival of a colony, a good contact rate (exposure of that ecotope to triatomines) and a good fecundity of colonies producing "infective females".

Competition and coexistence of two species

The model has been extended to two "species" of triatomines (i.e., true species or distinct groups of a same species not interbreeding). The basic hypothesis in the model is that the two "species" have the same ecological needs and therefore, can theoretically share the same niche.

Competitive exclusion

In the model, competitive exclusion or coexistence can be described. They depend on the R 0 values of both species colonies and the values of the γ ii and γ ij , , i.e., the "effects" of each species on the establishment of other individuals of its own species and the establishment of individuals of the other species. Generally, the species with the highest basic reproduction number (R 0 for colonies) will exclude the other, but coexistence may also occur in particular conditions (see below). Therefore, because R 0 's are directly dependent on fecundity and mortality rates, it would be interesting to compare these life traits between the two species. p.19

Coexistence

Coexistence may occur if one species affects the establishment of its own species more than the other species. Indeed the phase plane analysis of fig. 1b [START_REF] Canale | Peridomestic distribution of Triatoma garciabesi and Triatoma guasayana in north-west Argentina[END_REF].

The two species have slightly different ecotopes in the peridomicile and therefore, coexistence situations have been described. At a larger geographic scale, it has also been described for Rhodnius neglectus and R. nasutus [START_REF] Batista | Ecological niche modelling and differentiation between Rhodnius neglectus Lent, 1954 and Rhodnius nasutus Stål, 1859 (Hemiptera: Reduviidae: Triatominae) in Brazil[END_REF] and for Triatoma sordida and T. garciabesi [START_REF] Gurgel-Gonçalves | Geometric morphometrics and ecological niche modelling for delimitation of near-sibling triatomine species[END_REF]. As for T. infestans however, it does not seem that R 0 values are >>1, and there is no evidence that one group (domestic for example) may affect its proper establishment more than that of the other group if this latter enters the same site (house, in this case). Even if coexistence is possible, a species may still exclude another one if the R 0 of the competitively inferior species is sufficiently close to 1. In that case, for various situations where the R 0 of the triatomine species are low and species are in low densities, the competitive advantage of one species or the other will vary as R 0 values change between situations. Then, it might appear that in one site, one competitor will be present exclusively while in another site, it will be the other one.

The domestic / wild T. infestans populations.

The model clearly shows that the colonization of houses by wild T. infestans populations is unlikely. Personal field observations by mark-release-recapture technique indicate that T.

infestans individuals do not move easily between wild and domestic ecotopes, even if they are close (50 m for example). Therefore, mutual exposure (i.e., exposure of houses to wild triatomines and exposure of sylvatic environment to domestic triatomines) seems to be close to zero while exposure of each type of ecotope by each respective "species" is high and homogenous. If so, the model indicates that (see equations 11 and 12) it is hard for a wild population of T. infestans for which host heterogeneity in exposure is high (i.e., a low value of exposure for houses and a high value for a wild environment), to invade the already p.20 established domestic population where host heterogeneity is relatively low. For example, for the wild species invading a house, equation 11 indicates that its reproductive number will be small because ν Γ1 /ν Γ2 >>1. Invasion will therefore be very hard. Coexistence of wild and domestic populations of T. infestans is possible, but very unlikely to occur. Indeed, the model demonstrates that coexistence may occur if one species affects the establishment of its proper species more than the other species when both share the same site. For wild / domestic populations, this is unlikely as each group seems to have slightly different ecological needs and thus occupies different ecological niches. Genetics studies tend to distinguish two groups which do not mix [START_REF] Piccinali | Molecular population genetics and phylogeography of the Chagas disease vector Triatoma infestans in South America[END_REF]. However, a recent study indicated that gene flows can exist between the two groups [START_REF] Ceballos | Hidden sylvatic foci of the main vector of Chagas disease Triatoma infestans: Threats to the vector elimination campaign?[END_REF] but as stated, may represent invasion of sylvatic habitats by domestic populations finding here refuges after insecticide sprayings of houses. In that case, there is no competition, but invasion of an empty niche followed by intra specific movements of dispersion as one might expect. Inversely, if domestic individuals have survived vector control actions and still occupy their houses, the model tends to minimize the possible role of wild population in invading such a "not totally empty" niche.

Wild populations of T. infestans are mostly found close to the domestic/peridomestic environment (Brenière, comm. pers.) and the definition of "wild" is indeed difficult. Blood meal sources may help in defining a "wild" population which then feed on wild animals.

However, at present, the wild foci of T. infestans discovered are mostly under human influence, close to human settlings, except maybe in the Chaco region of Bolivia. For most of these foci, then, it seems that the domestic-peridomestic populations of T. infestans have invaded the wild environment and our model can account for such behaviour. If so, a new category of triatomine environment could be proposed: the "PARAdomestic" environment, in addition to the "domestic", "peri-domestic" and "sylvatic" terms yet in use. The paradomestic environment is then under the influence of the domestic -peridomestic environments, with some different characteristics (such as blood sources for triatomines that may be mainly of wild origin) and may permit population to undergo some level of differentiation.

When two groups are considered for coexistence or competition, our model assumption is that they do not interbreed (i.e., are real species). Below this working hypothesis, our model indicates that coexistence of wild and domestic populations of T. infestans is still possible, but very hard to achieve. If so, the colonization of domestic habitats by wild populations is p.21 unlikely. However, if further studies indicated that gene flows amongst wild and domestic T.

infestans are significant (i.e., that they interbreed and mix) other models should be used to describe their interactions. Models taking into account source-sink dynamics [START_REF] Pulliam | Sources, sinks, and population regulation[END_REF] or better, ecological traps [74 -75] should then be of great value and would likely explain the present heterogeneities observed amongst wild and close domestic T. infestans populations, such as genes and phenotypic differences or insecticide resistance. If gene flows exist amongst the two groups, then the population is single and the main question would be to identify likely source-sink dynamics where wild groups cannot persist without close domestic ones (and then without threats for vector control) or to identify mechanisms that enable the persistence of wild populations jeopardizing vector control. If so, distant wild populations without any anthropogenic influence should also be found. Mean (and variance) for contact rate of houses with triatomines

Figure captions

v Γ v Γ = σ Γ / m Γ coefficient

1 --

 1 females and triatomine colonies following a n+1 dimensional Markov process:-The number of fertilized females increases by one at a rate  The number of fertilized females decreases by one at a rate The number of colonies in one house increases by one at a rate )

Fig. 1 .

 1 Fig. 1. Phase plane analysis based on eq. (7) for two competing species. The point P 1 = (1/η 21 )(1-R 01 -1 ) correspond to the mean number of colonies of species 2 (C 2 ) when the isoclines dm C2 /dt=0 crosses the Y-axis. Equivalently, P 2 = (1/η 22 )(1-R 02 -1 ); P 3 = (1/η 11 )(1-R 01 -1 ); and P 4 =(1/η 12 )(1-R 02 -1 ) are defined as for P 1 .

Fig

  Fig.1a. corresponds to competition between the two species for a limiting resource, with

ρ 2 Fig. 1a

 2 Fig.1a

Figure

  Fig.1b

  C2 ). P 2 corresponds to the mean number of colonies of species 2 (i.e., m C2 ) when the isocline dm C2 /dt=0 crosses the y-axis (m C2 ). P 3 corresponds to the mean number of colonies of species 1 (i.e., m C1 ) when the isocline dm C1 /dt=0 crosses the x-axis (m C1 ). P 4 corresponds to the mean number of colonies of species 1 (i.e., m C1 ) when the isocline dm C2 /dt=0 crosses the x-axis (m C1 ). Figure 1a describes competition of the two species under the assumption  11 =  12 and  22 =

	1 ); P 2 = (1/η 22 )(1-R 02	-1 ); P 3 = (1/η 11 )(1-R 01	-1 ); and P 4 =(1/η 12 )(1-R 02	-1 ).
	P 1 corresponds to the mean number of colonies of species 2 (i.e., m C2 ) when the isocline
	dm C1 /dt=0 crosses the y-axis (m

 21 (equality of effects of species 1 (2) on the establishment of its own species or of species 2 (1)) . . The two isoclines dm C2 /dt=0 and dm C1 /dt=0 are parallel and therefore a mixed equilibrium cannot exist. The species with the highest basic reproduction number will exclude the other. In the figure, species 2 will exclude species 1 if R 02 > R 01 since this implies that P 2 > P 1 . This dynamics of competitive exclusion, based on the relative importance of the R 0 has been well described in other biological systems [65, 66 for example].

Figure

1b

describes the second scenario, i.e., coexistence. This scenario is possible if species i affects the establishment of its own species more than the other species j. In the model, this implies that  11 >  12 and  22 >  21.

  The developed formulas indicate that for R 01 >>1 and R 02 >>1, coexistence will exist when interspecific effects are weaker than intraspecific ones. Interspecific effects might be weaker than intraspecific ones if, for example, there is site segregation between the two species within human dwellings. Such situations have been observed for the peridomestic distribution of Triatoma garciabesi and Triatoma guasayana in north-west Argentina

			demonstrated that the two species
	will coexist if P1 > P2 and P4 > P3 [i.e., γ 11 (1-R 02	-1 ) > γ 12 (1-R 01	-1 ) and γ 22 (1-R 01	-1 ) > γ 21 (1-
	R 02	-1 )].

  .1a. corresponds to competition between the two species for a limiting resource, with parameters  11 =  12 and  22 =  21 (equality of effects of species 1 (2) on the establishment of its own species or of species 2 (1)) . Species 2 will exclude species 1 if R 02 > R 01 since this implies that P 2 > P 1 . Fig.1b. The two species will coexist if P 1 >P 2 and P 3 >P 4 .

Table captions Table 1 .

 captions1 Model parameters and variables

Table 1 .

 1 Model parameters and variables

	Parameter	Description
	n	Number of human dwellings considered
	F	Size of the population of triatomines with the potential of colonizing a house (i.e.,
		"fertilized females" which are seeking for a human dwelling)
	m F	Mean number of "fertilized females"
	μ F	Per capita rate of mortality of "fertilized females"
	C i	Number of triatomines colonies (of one species) in the house i (burden of
		triatomine colonies for house i)
	m C	Mean number of triatomine colonies in a house
	μ C	Per capita rate of disappearance of triatomine colonies

i  Contact rate of a house i with "fertilized females"

m Γ (σ Г 2 )

Acknowledgements

The author acknowledges the support of (1) the French "Agence Nationale de la Recherche" (ANR), under grant ALTERNATIVE (ANR-2010 CESA 018 01) ''What alternatives to the chemical control of Triatoma infestans ?" and (2) OMS-TDR / IDRC under grant EBS-LAC n°A90281 "An eco-bio-social approach to implement integrated techniques to control T. infestans in Bolivian poor communities" and (3) Sam Featherston for correcting the English spelling.

Disclosure of interest

The author declares that he has no conflict of interest concerning this article.