#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast


Authors: Aimee R. Taylor aff001;  Diego F. Echeverry aff003;  Timothy J. C. Anderson aff006;  Daniel E. Neafsey aff002;  Caroline O. Buckee aff001
Authors place of work: Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, Massachusetts, USA aff001;  Infectious Disease and Microbiome Program, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA aff002;  Centro Internacional de Entrenamiento e Investigaciones Médicas (CIDEIM), Cali, Colombia aff003;  Universidad Icesi, Calle 18 No. 122-135, Cali, Colombia aff004;  Departamento de Microbiologia, Facultad de Salud, Universidad del Valle, Cali, Colombia aff005;  Disease Intervention and Prevention Program, Texas Biomedical Research Institute, San Antonio, Texas, USA aff006;  Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, Massachusetts, USA aff007
Published in the journal: Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast. PLoS Genet 16(11): e1009101. doi:10.1371/journal.pgen.1009101
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1009101

Summary

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies.

Keywords:

Cities – Colombia – DNA recombination – Malaria – Malarial parasites – Parasitic diseases – Plasmodium – Single nucleotide polymorphisms

Introduction

In many research fields genetic data are used to help characterise connectivity between geographically distinct biological populations, with numerous applications in conservation, agriculture, and public health. Patterns of genetic similarity between pathogen populations help us understand how the disease spreads. Patterns of relatedness (a measure of genetic similarity) between malaria parasites sampled from different human populations, for instance, help characterise the connectivity between different malaria parasite populations, thus guide the design of targeted public health interventions [1].

Several methods are employed to measure genetic similarity and thus characterise connectivity. Phylogenetic methods, in which genetic distances between individuals are measured in units of mutation [2], are most applicable to rapidly mutating organisms that do not recombine (e.g RNA viruses) [3]. Studies of relatedness, in which relatedness is a measure of probability of inter-individual identity-by-descent (IBD), are applicable to organisms that do recombine (e.g. malaria parasites). Population genetic parameters of allelic variation (e.g. FST) are applicable to all organisms (those that do and do not recombine), but do not generate measures of genetic distance or similarity on an inter-individual level, thus provide less granularity. Moreover, among recombining organisms, inter-population allelic variation tends to accumulate more slowly than inter-individual variation in IBD [4]. As such, analyses of relatedness sometimes recover evidence of nearby and recent connectivity where analyses of FST do not [5].

Malaria parasites are protozoan parasites that undergo an obligate stage of sexual recombination in the mosquito midgut. Like many organisms (e.g. many plants [6, 7]), malaria parasites have a mixed mating system that encompasses both inbreeding and outcrossing. The extent to which malaria parasites outcross depends on transmission intensity and is not fully understood [8]. For outcrossing to occur a mosquito must ingest genetically distinct gametocytes. Humans can be infected by multiple genetically distinct parasite clones that are either co-transmitted via inoculation from a single mosquito, in which case they are likely recombinants so inter-related (unless they derive from different blood meals), or transmitted independently by multiple mosquitoes (a mechanisms coined superinfection by George MacDonald, 1950 [9, 10]), in which case the parasite clones are likely unrelated [11, 12]. The latter can occur in a setting where the entomological inoculation rate is high; recent work suggests co-transmission is important in both low and high transmission settings [12].

Malaria genomic epidemiology studies of connectivity are increasingly common, especially in the context of public health and using genotype (versus whole genome sequence) data [5, 1316]. Using IBD-based relatedness but not FST, evidence of isolation-by-distance among P. falciparum populations along a 100 km stretch of the Thailand-Myanmar border was found [5]. This study was based, in part, on analyses of monoclonal P. falciparum samples genotyped at 93 single nucleotide polymorphisms (SNPs). Based on FST estimated using P. falciparum samples genotyped at 250 SNPs, a different study found evidence of departure from isolation-by-distance among P. falciparum populations along a 500 km stretch of the Colombian-Pacific coast where transmission is mixed (low but high in some regions) and outcrossing limited [13, 17]. In the current study, we re-explore this departure from isolation-by-distance with more granularity using IBD-based relatedness. For the first time in malaria epidemiology (and, to our knowledge, for the first time in ecological and epidemiological studies more generally), we account for uncertainty in relatedness estimates; we also use threshold-free methods to compare parasite populations and identify clonal components. The original study [13] is described in more detail below.

Malaria epidemiology in Colombia is associated with a multitude of ecological, evolutionary and social factors, including human migration due to deforestation, illegal crops, gold mining [1822], and the mass emigration of people fleeing the humanitarian crisis in Venezuela [2326]. Understanding the interplay between e.g. human migration, parasite population connectivity and the spread of antimalarial resistance is critical [18, 20]. For example, if resistance is driven by spread (versus de-novo mutation), targeted efforts to eliminate hotspots of transmission (e.g. in eastern Myanmar [27, 28]) may help to prolong the longevity of compromised antimalarial therapies. To ensure adequate isolation, thereby prevent re-population, units of targeted intervention need to account for parasite population connectivity, which relates to human migration [2931]. In preparation for studies of resistance, Echeverry et al. genotyped P. falciparum samples from four provinces on the Colombian-Pacific coast [13]. Clonality, population structure and linkage disequilibrium (LD) were characterised using a suite of population genetic analyses. The results were highly informative: the vast majority of successfully genotyped P. falciparum samples were deemed monoclonal (325 of 400) with a strong association between incidence and clonality. Among the 325 monoclonal samples, 136 unique haploid multilocus genotypes (MLGs) were identified using relatedness based on identity-by-state (IBS), which is a correlate of IBD [32] (and has been used elsewhere to characterise connectivity between nearby malaria parasite populations [1416]). Of the 136 MLGs, 44 infected two or more patients (max. 28 patients), 45 persisted for two or more days (max. 8 years), and 7 of the 15 most common MLGs were sampled in two or more provinces (max. all four provinces). Panmixia was rejected based on evidence of four sympatric but geographically structured subpopulations; and, overall, LD decayed at a rate that was faster than expected for South American P. falciparum populations (compare with e.g. [33]). Echeverry et al. concluded that evidence of low genetic diversity, persistent MLGs and population structure is consistent with low transmission and limited outcrossing, while evidence of a relatively fast rate of LD decay and of shared MLGs across provinces is consistent with extensive human movement connecting P. falciparum populations.

Although the study by Echeverry et al. features analyses of IBS-based relatedness (i.e. MLGs), evidence of departure from isolation-by-distance was based on FST alone. To explore isolation-by-distance in more granularity while accounting for uncertainty, we compute IBD-based relatedness estimates and confidence intervals for all pairs of 325 monoclonal parasite samples. Akin to previous studies (e.g. [5]), highly-related parasites were classified using a threshold; however, confidence intervals allow uncertainty to be accounted for in this study. For example, in [5] a parasite pair was considered highly-related if its relatedness estimate exceeded 0.5, whereas here a parasite pair is considered highly-related if the lower end-point of the 95% confidence interval around its relatedness estimate exceeds some stated value, which is 0.25 in the main text and 0.5 in sensitivity analyses. This is important because uncertainty can overwhelm relatedness estimated using limited genotype data [32]. Our approach includes two additional contributions. First, we complement our analysis of highly-related parasites with a threshold-free approach that uses a metric called the 1-Wasserstein distance, which can be interpreted as the cost of transporting a distribution of parasite samples from one city to another [18, 34]. Second, we identify groups of statistically indistinguishable parasites, which we call clonal components, using the simple concept of components from graph theory and confidence intervals. Confidence intervals circumvent reliance on an arbitrary clonal threshold (i.e. some number of differences tolerated between parasites samples considered clonal). Graph components circumvent reliance on unsupervised clustering methods that are sensitive to both the definition of genetic similarity and algorithmic specification [35, 36]. Overall, our approach could be adapted to viruses and bacteria that show recombination or reshuffling of segments as well as clonal propagation [3740], to other protozoans (e.g. Toxoplasma, Cryprosporidium [4143]), and to the many fungal pathogens [44], plants [6, 7], and animals with mixed mating systems. Due to our treatment of uncertainty, it is especially relevant for a growing number of studies that plan to estimate IBD-based relatedness using genotype (versus sequence) data.

Results

Relatedness estimates between P. falciparum sample pairs

For all 52650 pairwise comparisons of 325 previously published monoclonal P. falciparum samples with data on 250 biallelic single nucleotide polymorphisms (SNPs) [13], relatedness was estimated using the hidden Markov model (HMM) described in [32]. Relatedness is thus defined as the probability that, at any SNP, the two alleles drawn from the paired monoclonal P. falciparum samples are IBD.

The parasite samples were collected between 1993 and 2007 from symptomatic patients participating in studies at five cities on the Colombian-Pacific coast (S1 Table). Despite considerable uncertainty, all estimates are informative (Fig 1). That is to say, there are no relatedness estimates whose 95% confidence intervals span entirely from zero to one. The vast majority of relatedness estimates were classified unrelated.

Fig. 1. Estimates of relatedness with 95% confidence intervals.
Estimates of relatedness with 95% confidence intervals.
Estimates and confidence intervals are shown for all 325 choose two (52650) P. falciparum sample pairs and are ordered by increasing relatedness estimate. Confidence intervals are coloured according to classifications based on lower and upper confidence interval end-points, where τ is an arbitrary threshold used to classify highly-related pairs. For example, a pair is considered highly-related with τ = 0.25 if the lower end-point of the confidence interval around its relatedness estimate exceeds 0.25. Otherwise stated, if its relatedness estimate is statistically distinguishable from 0.25.

Highly-related P. falciparum sample pair fractions partitioned in space and time

In our main analysis (Fig 2), highly-related parasite samples were classified using an arbitrary threshold of 0.25 (Table 1), which corresponds to the expected relatedness between parasites separated by two outcrossed generations, but is hard to interpret in the context of frequent clonal propagation. Despite few highly-related P. falciparum sample pairs overall, there are three notable observations regarding their fraction partitioned in space and time. First, there is a greater fraction of highly-related sample pairs among those collected closer together in time (Fig 2(A)). Second, the fraction of highly-related sample pairs is generally greater within cities than between, with Guapi having the largest fraction of highly-related pairs and Buenaventura having the lowest (Fig 2(B)). However, third, the fraction shared between Buenaventura and Tumaco is greater than expected given inter-city distance (Fig 2(B)). These observations are largely robust to different high-relatedness thresholds (S1 Fig). Spatial trends evaluated using a threshold-free approach are also consistent: they show a general increase in 1-Wasserstein distance with inter-city distance besides Buenaventura and Tumaco (Fig 3). The 1-Wasserstein distance can be interpreted as the total cost required to transport a distribution of parasite samples from one city to another [18, 34], where the cost of transporting a single parasite to another is equal to one minus relatedness. The small 1-Wasserstein distance between Buenaventura and Tumaco is thus consistent with elevated gene flow between P. falciparum populations sampled from these cities.

Tab. 1. Classification of parasite sample pairs.
Classification of parasite sample pairs.
Classification is based on the lower and upper end-points (LCI and UCI, respectively) of the 95% confidence interval around each relatedness estimate, r ^, where ϵ is an arbitrarily small number to identify LCI ≈ 0 and UCI ≈ 1 given that LCI and UCI ∈ (0, 1) not [0, 1]; and τ is an arbitrary threshold used to classify highly-related pairs. We use ϵ = 0.01 throughout, τ = 0.25 (main analysis) and τ ∈ {0.25, 0.50} (sensitivity analysis).
Fig. 2. Fractions of highly-related sample pairs partitioned in time and space.
Fractions of highly-related sample pairs partitioned in time and space.
(A) Partitioned by time between collection dates. (B) Partitioned by collection city, where the inter-city great-circle distance is the distance in kilometres (km) between city pairs on the Earth’s surface.
Fig. 3. P. falciparum population connectivity assessed using a threshold-free approach.
<i>P. falciparum</i> population connectivity assessed using a threshold-free approach.
1-Wasserstein distance between parasite populations from different cities versus inter-city great-circle distance in kilometres (km).

Fig 4 shows the inter-city P. falciparum population connectivity of Fig 2(B) projected onto a map of the Colombian-Pacific coast. Buenaventura and Tumaco are the two largest official ports on the Colombian-Pacific coast (Buenaventura is the largest) and are connected by frequent marine traffic (www.marinetraffic.com). Although Tumaco is connected to Buenaventura via the Pan-American highway, which connects all sites but Guapi, primary access to Tumaco is via the port due to difficult and unsafe country roads in Nariño. Guapi, which is effectively unreachable by road and not an official port, is connected by marine traffic but with less frequency (www.marinetraffic.com). Consistent with its isolation, the fraction of highly-related parasite pairs is relatively large within Guapi (Fig 2(B)), and very small between Guapi and the two inland cities, Quibdó and Tadó (Figs 2(B) and 4). Moreover and importantly regarding the elevated fraction of highly-related samples pairs within both Guapi and Tadó (Fig 2(B)), all samples from Guapi and Tadó were collected within a single year (S1 Table). The low fraction of highly-related parasite sample pairs within Buenaventura (Fig 2(B)) is in part consistent with it having contributed samples over many years (S1 Table) and with it being the most important port on the Pacific coast (www.marinetraffic.com), i.e. a hub through which human traffic and thus potential parasite mixing is high [13].

Fig. 4. P. falciparum population connectivity based on fractions of highly-related sample pairs.
<i>P. falciparum</i> population connectivity based on fractions of highly-related sample pairs.
The width of each inter-city edge is proportional to the fraction of highly-related sample pairs across cities plotted in Fig 2(B). Note that the edges between Guapi and Quibdó and Guapi and Tadó are plotted but too thin to visually discern.

The apparent association between P. falciparum population connectivity and the frequency of marine traffic raises questions about the latter’s role in malaria transmission. However, other scenarios could lead to these relationships, for example high connectivity could result from a single travel event between Buenaventura and Tumaco, followed by expansion of highly-related and clonal parasites. To further explore the genetic signal that supports this association we next consider clonal components.

Clonal components

We define clonal components as groups of statistically indistinguishable parasite samples identified under a graph theoretic framework: consider a graph whose vertices are parasite samples and whose edges are weighted by relatedness estimates, a clonal component is a sub-graph whose vertices are all connected to one another via edges whose weights are statistically indistinguishable from one (i.e. clonally related, Table 1). In total, 46 distinct clonal components were detected, ranging in size from 2 to 28 statistically indistinguishable parasite samples (Fig 5). They are spatially clustered. Ten of the 46 contain parasite samples collected from two or more cities. Each clonal component besides one (clonal component four) is on average related to at least one other (Fig 5). The unrelated clonal component is almost certainly an artefactual contaminant: it accords with MLG 036 reported in [13], where contamination during in vitro adaptation or DNA manipulation was suspected (MLG 036 contained “two culture-adapted samples from Quibdó and Tadó that were indistinguishable from the Dd2 reference strain from Southeast Asia”—the Dd2 reference strain was included as a control when the data were originally generated [13]).

Fig. 5. Clonal components and the average relatedness between them.
Clonal components and the average relatedness between them.
Vertices depict clonal components, which are groups of two or more statistically indistinguishable parasite samples. CC vertices are plotted using the Fruchterman-Reingold layout algorithm [45], thereby clustering inter-related CCs. The size of each CC vertex is proportional to the number of parasite samples per CC, ranging from 2 to 28 statistically indistinguishable parasite samples. CCs are named in order of the collection date of the earliest parasite sample per CC (S2 Table). CCs with parasite samples collected from two or more cities are depicted as pie charts. Colour denotes the city of parasite sample collection. Edge transparency and weight is proportional to average relatedness, ranging from 0.003 to 0.840. Relatedness estimates that are indistinguishable from zero were set to zero. Edges whose average relatedness is zero are not plotted. Each CC besides CC4 is related to at least one other. CC4 contains two samples (one from Tadó, another from Quibdó). It is likely a contaminant; see main text.

Clonal parasite samples detected in both Buenaventura and Tumaco belong to five distinct clonal components (1, 12, 14, 20 and 40, Fig 5). We thus dismiss a single travel event connecting Buenaventura and Tumaco involving a single parasite clone. We cannot dismiss a single travel event involving multiple parasite clones, however. Based on the proportions of multiclonal infections in the original data from Buenaventura and Tumaco (14% and 19% respectively) [13], the probability that these five clones could be distributed across infections in four or fewer individuals is approximately 0.6. Indeed, three of the five clonal components are inter-related on average (S3 Table). As such, they could derive from co-transmitted recombinant parasites transported in a single individual with a multiclonal infection. On the contrary, the remaining two clonal components have relatedness estimates that are not statistically distinguishable from zero. As such, they could derive from a single superinfected individual, or from different individuals with independent monoclonal infections. Unfortunately, the data required to further evaluate these scenarios (data on the multiplicity of multiclonal infections, and on relatedness within multiclonal infections, e.g. [12]) are not available. Given dates and cities of first detection (S2 Table), it is tempting to suggest some clonal components predate others and originate in specific locations. For example, it is possible that parasite samples from clonal components 1 and 20 in Buenaventura and Tumaco emanated from Guapi, creating a spurious link between Buenaventura and Tumaco. However, because these data are from sparsely sampled symptomatic cases in a setting where clonal propagation is frequent, sample collection chronology is not necessarily representative of the chronology of transmission chain events (S2 Fig).

Regarding transmission chain events, we note that clonal component 20 relates to the three inter-related clonal components (1, 12 and 14) via an intermediate clonal component detected in Tumaco only (clonal component 15) as well as an intermediate parasite sample from Quibdó that does not belong to a clonal component (S3 Fig). These intermediates likely derive from recombination between parasites related to the clonal components they connect. Several connections consistent with recombinants can be found among the relatedness graphs (Fig 5 and S3 Fig). As such, it seems it may be at least theoretically possible to construct approximate P. falciparum transmission chains given more dense sampling of malaria infections on the Colombian-Pacific coast.

Discussion

Here we show that estimates of IBD-based relatedness and their associated uncertainty can be used to uncover evidence of epidemiologically meaningful connectivity between P. falciparum populations on a relatively local spatial scale: along the Colombian-Pacific coast where clonal propagation is frequent [13], extending southward to Ecuador [46, 47]. While our approach largely confirms a previous report based on FST [13], estimates of relatedness provide more granularity while their confidence intervals account for uncertainty thus provide more statistical rigor, e.g. when highly-related parasite sample pairs are classified. Our approach includes two additional contributions: 1-Wasserstein distances are used to compare parasite populations in an entirely threshold-free manner; and clonal components are identified using graph components and confidence intervals, thereby circumventing reliance on an arbitrary clonal threshold. Threshold-free methods are especially important in analyses where thresholds do not have clear interpretations (e.g. 0.5 may correspond to the expected relatedness of siblings in an outcrossed population, but its interpretation is unclear in a population where inbreeding and clonal propagation is common) and thus undermine the cross-comparison of studies. Standardisation will accelerate the maturation of malaria genomic epidemiology and facilitate the translation of research into actionable insight for policy makers [1]. Our overall approach could also be adapted for analyses of other recombining organisms with mixed mating systems.

IBD-based relatedness estimates recovered 1) a large fraction of highly-related parasite sample pairs within Guapi, a city on the Colombian-Pacific coast that is relatively isolated besides infrequent marine traffic; 2) a low fraction of highly-related parasite sample pairs within Buenaventura, the most important port on the Colombian-Pacific coast and thus the least isolated city in this study; and 3) a disproportionally large fraction of highly-related parasite pairs between Buenaventura and Tumaco (departure from isolation-by-distance), where Tumaco is the second largest port on the Colombian-Pacific coast. These observations accord with several published previously: 1) elevated LD in a P. falciparum subpopulation (identified using STRUCTURE [35, 48]) predominant in Guapi; 2) rapid LD decay in a P. falciparum subpopulation predominant in Buenaventura; and 3) lowest genetic differentiation (based on FST estimates) between provinces Valle (Buenaventura) and Nariño (Tumaco) [13]. LD, STRUCTURE and FST analyses all rely on allelic variation. The concordance between results based on relatedness and allelic variation suggests that P. falciparum outbreeding on the Colombian-Pacific coast is infrequent enough that both types of analyses generate insight on approximately the same time scale.

The aforementioned results generate hypotheses around the frequency of marine traffic and malaria transmission on the Colombian-Pacific coast. Notwithstanding long-range windborne dispersal, which may be critical for malaria transmission in Africa [49], anopheline flight range is generally small (around 3.5 km [50]). As such, long-range malaria parasite dispersal on the Colombian-Pacific coast is almost certainly human-mediated. A recent study of P. vivax proposed that human movement across a “malaria corridor” stretching from the northwest to the south of the Colombian-Pacific Coast likely promotes P. vivax gene flow, and that mining activities may provide transmission “contact zones” [51], similarly proposed for P. falciparum [22]. P. falciparum population connectivity is consistent with the human “malaria corridor” hypothesis, especially since it correlates with accessibility, not isolation-by-distance. Both infected humans and mosquitoes are compatible with this hypothesis, i.e. checks for infected Anopheles spp. on boats may be merited [52, 53]. However, relatively high differentiation between populations of An. albimanus (one of the three primary vectors of malaria in Colombia [54]) from Buenaventura and Tumaco [55] points towards human carriage.

The Colombian-Pacific coast has long been associated with international trade, but until recently human migration in the region was largely domestic. The flow of Venezuelan migrants infected with Plasmodium spp. has increased in recent years: of 965, 1774 and 2288 non-domestic malaria cases reported in Colombia in 2017, 2018 and 2019, respectively, 882 (91.4%), 1684 (94.9%), and 2190 (95.7%) were from Venezuela [5658]. Other non-domestic sources of malaria in Colombia include countries elsewhere in South America (e.g. Peru, Panama, French Guyana, Ecuador, Brazil) and several African countries (e.g. Uganda, Republic of the Congo, Nigeria, Ivory Coast, Cameroon, Angola) [58]. Some of the infected Venezuelan nationals are migrating southward to Ecuador and Peru [24]. Other non-domestic cases may be associated with the traffic of people who arrive at Colombian ports with a view towards northward travel e.g. to the USA via Central America and Panama [59]. Genetic surveillance of “international parasites” may help malaria control efforts in Colombia.

The evidence we find of connectivity between P. falciparum populations may be unique to the period of time over which the data were collected (1993-2007). This was a period of historically high malaria case counts in Colombia [17], as well as social instability in the South Pacific region. Contemporary data on more densely sampled cases and on mosquito and human movement are required to characterise extant connectivity, its reach beyond Colombia (see e.g. [47]), and to rule out alternative hypotheses. Regarding alternative hypotheses, heterogeneous vectorial capacity and antimalarial drug pressure could selectively enhance parasite survival in such a way that generates apparent connectivity between Buenaventura and Tumaco, e.g. if parasites are adapted to local vectors whose distributions are more similar between Buenaventura and Tumaco than elsewhere. Although adult An. albimanus B and An. neivai s.l. have been detected in the vicinities of both cities [55, 60], the species distributions in the vicinities of Buenaventura and Tumaco differ more than those in the vicinities of Tumaco and Guapi [60]. As such, heterogeneous vectorial capacity seems an unlikely alternative hypothesis. Similarly, relatedness may be greater among parasites with comparable antimalarial resistance: a recent study of South East Asian P. falciparum parasites found greater relatedness in the recent past among parasites with artemisinin resistance mutations versus without [61]. This study used size-stratified IBD segments to date relatedness [61]. On the Colombian-Pacific coast, IBD segment size inference could help identify some recently related parasites. However, it requires whole genome sequence data and is hard (if not presently impossible) to interpret in the face of frequent clonal propagation [32]. The development of an ancestral recombination model that incorporates transmission-dependent selfing is a research priority in malaria genomic epidemiology and would aid research on other organisms that show both outbreeding and clonal propagation.

Materials and methods

Data

This study relies entirely on previously published data that are publicly available [13, 32]. In the original study by Echeverry et al., finger-prick blood spot samples were obtained from patients with symptomatic uncomplicated malaria [13]. Samples were collected between 1993 and 2007 from five cities in four provinces: Tadó and Quibdó in Chocó, Buenaventura in Valle, Guapi in Cauca and Tumaco in Nariño (S1 Table) [13]. Informed consent was obtained from all the subjects enrolled, as approved by the CIDEIM Institutional Review Board (IRB) [13]. The Colombian-pacific coast is one of the rainiest regions of the world [55, 62]. At that time, Colombia had approximately 100,000 malaria cases per year [13, 17]. Collectively Chocó, Valle, Cauca and Nariño accounted for up to 75% of the P. falciparum cases reported, with relatively high transmission in Chocó and relatively low transmission in Valle and Cauca [13].

The data that feature in this descriptive study also feature in a recent methodological study concerning data requirements for relatedness inference [32]. As in [32], we did not post-process the data in any way besides mapping SNP positions to the P. falciparum 3d7 v3 reference genome and recoding heteroallelic calls as missing (since all samples with fewer than 10 heteroallelic SNP calls were classified monoclonal previously [13]). The monoclonal data include 325 P. falciparum samples with data on 250 biallelic SNPs whose minor allele frequency estimates (the minor allele sample count divided by 325) range from 0.006 to 0.495 (S4 Fig).

Relatedness inference and classification of parasite sample pairs and groups

For each pairwise parasite sample comparison, we generated a relatedness estimate and 95% confidence interval using the HMM and parametric bootstrap described in [32]. Sample pairs were classified as unrelated, related, highly-related and clonal using confidence interval end-points as follows and summarised in Table 1. A pair was classified unrelated if its relatedness estimate, r ^, was statistically indistinguishable from zero with lower confidence interval end-point (LCI) less than ϵ, an arbitrarily small number to identify LCI ≈ 0 and UCI ≈ 1 given that LCI and UCI ∈ (0, 1) not [0, 1]. A pair was classified related if its relatedness estimate, r ^, was statistically distinguishable from zero with LCI > ϵ. A pair was considered highly-related if its relatedness estimate, r ^, was statistically distinguishable from some specified threshold, τ, with LCI > τ. A pair was considered clonal if its relatedness estimate, r ^, was statistically indistinguishable from one with upper confidence interval end-point (UCI) > 1 − ϵ. Note that these classifications are possible because all estimates are informative, i.e. no confidence intervals span the entire zero to one range (Fig 1). These classifications are neither necessarily exclusive nor conversely true: a clonal parasite pair is related, but a related parasite pair is not necessarily clonal. Throughout, ϵ = 0.01. In the main analysis (Fig 2) τ = 0.25, in the sensitivity analysis (S1 Fig) τ ∈ {0.25, 0.50}.

In addition to classifying parasite sample pairs, we classify groups of statistically indistinguishable parasite samples, which we call clonal components because they are defined using the simple concept of components from graph theory. First, we construct a super-graph whose vertices are parasite samples connected by edges that are weighted by relatedness estimates. Within the super-graph, a clonal component is a sub-graph within which all parasite samples are connected to one another (directly or not) via edges whose weights are statistically indistinguishable from one, while being connected to parasites samples outside the sub-graph via edges whose weights are not statistically indistinguishable from one. Clonal components tend to be fully connected (i.e. all parasite samples within the clonal component are directly connected to one another by edges whose weights are statistically indistinguishable from one). The igraph package [63] in R [64] was used to identify clonal components and to visualise them using the Fruchterman-Reingold layout algorithm [45].

Spatiotemporal trends in P. falciparum population connectivity

Spatiotemporal trends in population connectivity were explored visually by partitioning parasite sample pairs by their collection cities and dates, then plotting the per-partition fraction of highly-related pairs. Inter-city great-circle distance was calculated using the Haversine formula, which assumes the earth is spherical. Error bars were constructed by re-sampling per-partition parasite sample pairs 100 times with replacement and taking the 2.5th and 97.5th percentiles of the fraction of highly-related pairs as the lower and upper limits, respectively. Sensitivity to τ = 0.25 (high relatedness threshold used in Fig 2) was explored using an alternative τ = 0.50 (S1 Fig) and also by using a threshold-free approach (Fig 3) as follows.

To explore population connectivity using a threshold-free approach, we calculated 1-Wasserstein distances between groups of parasite samples from different cities using the transport [65] package in R [64]. Specifically, for a pair of cities a and b, we construct a na × nb genetic distance matrix, G, of 1 - r ^ i j (where na and nb are the parasite sample counts from cities a and b, respectively, i = 1, …, na and j = 1, …, nb) and two vectors wa=(1na,…,1na) and wb=(1nb,…,1nb) of length na and nb, respectively. We then calculate the 1-Wasserstein distance, which minimises the total cost of transporting wa to wb, where 1 - r ^ i j is the cost of transporting a single unit, using transport::transport(wa, wb, costm = G, method = "shortsimplex"). This amounts to treating parasite samples from different cities as draws from different distributions, where the 1-Wasserstein distance can be interpreted as the cost required to transport a distribution of parasite samples from one city to another [18, 34]. Since per-city parasite sample sizes differ, transportation requires the expansion (or contraction) of parasite mass in addition to the transportation of individual units. City pairs with smaller 1-Wasserstein distances are interpreted as having greater connectivity between the P. falciparum populations collected from them. Error bars were constructed by re-sampling parasite sample pairs per inter-city partition 100 times with replacement and taking the 2.5th and 97.5th percentiles of the distribution of 1-Wasserstein distances based on the re-sampled sample pairs as the lower and upper limits, respectively.

Supporting information

S1 Table [pdf]
Yearly monoclonal sample counts per city.

S2 Table [pdf]
A summary of all clonal components.

S3 Table [ccs]
Average relatedness between select clonal components.

S1 Fig [lci]
Fractions of highly-related sample pairs partitioned in time and space: Sensitivity to the high-relatedness threshold.

S2 Fig [tif]
Sample collection chronology does not reflect transmission chain event chronology.

S3 Fig [ccs]
Clonal components and singletons and the average relatedness between them.

S4 Fig [tif]
Minor allele frequency estimates.


Zdroje

1. Dalmat R, Naughton B, Kwan-Gett TS, Slyker J, Stuckey EM. Use cases for genetic epidemiology in malaria elimination. Malaria journal. 2019;18(1):163. doi: 10.1186/s12936-019-2784-0

2. Holder M, Lewis PO. Phylogeny estimation: traditional and Bayesian approaches. Nature reviews genetics. 2003;4(4):275–284. doi: 10.1038/nrg1044

3. Biek R, Pybus OG, Lloyd-Smith JO, Didelot X. Measurably evolving pathogens in the genomic era. Trends in ecology & evolution. 2015;30(6):306–313. doi: 10.1016/j.tree.2015.03.009

4. Thompson EA. Identity by descent: variation in meiosis, across genomes, and in populations. Genetics. 2013;194(2):301–326. doi: 10.1534/genetics.112.148825

5. Taylor AR, Schaffner SF, Cerqueira GC, Nkhoma SC, Anderson TJ, Sriprawat K, et al. Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent. PLoS genetics. 2017;13(10):e1007065. doi: 10.1371/journal.pgen.1007065 29077712

6. Grant AG, Kalisz S. Do selfing species have greater niche breadth? Support from ecological niche modeling. Evolution. 2020;74(1):73–88. doi: 10.1111/evo.13870

7. Mattila TM, Laenen B, Slotte T. Population genomics of transitions to selfing in Brassicaceae model systems. In: Statistical Population Genomics. Springer; 2020. p. 269–287.

8. Siegel SV, Rayner JC. Single cell sequencing shines a light on malaria parasite relatedness in complex infections. Trends in Parasitology. 2020;36(2):83–85. doi: 10.1016/j.pt.2019.12.007

9. Macdonald G, et al. The analysis of infection rates in diseases in which super infection occurs. Tropical diseases bulletin. 1950;47:907–915. 14798656

10. Nåsell I. On superinfection in malaria. Mathematical Medicine and Biology: A Journal of the IMA. 1986;3(3):211–227. doi: 10.1093/imammb/3.3.211

11. Nkhoma SC, Nair S, Cheeseman IH, Rohr-Allegrini C, Singlam S, Nosten F, et al. Close kinship within multiple-genotype malaria parasite infections. Proceedings of the Royal Society B: Biological Sciences. 2012;279(1738):2589–2598. doi: 10.1098/rspb.2012.0113 22398165

12. Nkhoma SC, Trevino SG, Gorena KM, Nair S, Khoswe S, Jett C, et al. Co-transmission of Related Malaria Parasite Lineages Shapes Within-Host Parasite Diversity. Cell Host & Microbe. 2020;27(1):93–103. doi: 10.1016/j.chom.2019.12.001 31901523

13. Echeverry DF, Nair S, Osorio L, Menon S, Murillo C, Anderson TJ. Long term persistence of clonal malaria parasite Plasmodium falciparum lineages in the Colombian Pacific region. BMC genetics. 2013;14(1):2. doi: 10.1186/1471-2156-14-2

14. Omedo I, Mogeni P, Rockett K, Kamau A, Hubbart C, Jeffreys A, et al. Geographic-genetic analysis of Plasmodium falciparum parasite populations from surveys of primary school children in Western Kenya. Wellcome open research. 2017;2:29. doi: 10.12688/wellcomeopenres.11228.2 28944299

15. Omedo I, Mogeni P, Bousema T, Rockett K, Amambua-Ngwa A, Oyier I, et al. Micro-epidemiological structuring of Plasmodium falciparum parasite populations in regions with varying transmission intensities in Africa. Wellcome open research. 2017;2:10. doi: 10.12688/wellcomeopenres.10784.1

16. Tessema S, Wesolowski A, Chen A, Murphy M, Wilheim J, Mupiri AR, et al. Using parasite genetic and human mobility data to infer local and cross-border malaria connectivity in Southern Africa. Elife. 2019;8:e43510. doi: 10.7554/eLife.43510

17. Rodríguez JCP, Uribe GÁ, Araújo RM, Narváez PC, Valencia SH. Epidemiology and control of malaria in Colombia. Memórias do Instituto Oswaldo Cruz. 2011;106:114–122. doi: 10.1590/S0074-02762011000900015

18. Feged-Rivadeneira A, Ángel A, González-Casabianca F, Rivera C. Malaria intensity in Colombia by regions and populations. PloS One. 2018;13(9). doi: 10.1371/journal.pone.0203673 30208075

19. Castellanos A, Chaparro-Narváez P, Morales-Plaza CD, Alzate A, Padilla J, Arévalo M, et al. Malaria in gold-mining areas in Colombia. Memorias do Instituto Oswaldo Cruz. 2016;111(1):59–66. doi: 10.1590/0074-02760150382 26814645

20. Recht J, Siqueira AM, Monteiro WM, Herrera SM, Herrera S, Lacerda MV. Malaria in Brazil, Colombia, Peru and Venezuela: current challenges in malaria control and elimination. Malaria journal. 2017;16(1):273. doi: 10.1186/s12936-017-1925-6

21. Daniels JP. Increasing malaria in Venezuela threatens regional progress. The Lancet Infectious Diseases. 2018;18(3):257. doi: 10.1016/S1473-3099(18)30086-0

22. Knudson A, González-Casabianca F, Feged-Rivadeneira A, Pedreros MF, Aponte S, Olaya A, et al. Spatio-temporal dynamics of Plasmodium falciparum transmission within a spatial unit on the Colombian Pacific Coast. Scientific Reports. 2020;10(1):1–16. doi: 10.1038/s41598-020-60676-1 32111872

23. Grillet ME, Villegas L, Oletta JF, Tami A, Conn JE. Malaria in Venezuela requires response. Science. 2018;359(6375):528–528.

24. Jaramillo-Ochoa R, Sippy R, Farrell DF, Cueva-Aponte C, Beltrán-Ayala E, Gonzaga JL, et al. Effects of political instability in Venezuela on malaria resurgence at Ecuador–Peru border, 2018. Emerging infectious diseases. 2019;25(4):834. doi: 10.3201/eid2504.181355 30698522

25. Daniels JP. Venezuela in crisis. The Lancet Infectious Diseases. 2019;19(1):28. doi: 10.1016/S1473-3099(18)30745-X

26. Rodríguez-Morales AJ, Suárez JA, Risquez A, Villamil-Gómez WE, Paniz-Mondolfi A. Consequences of Venezuela’s massive migration crisis on imported malaria in Colombia, 2016-2018. Travel Medicine and Infectious Disease. 2019;28:98–99. doi: 10.1016/j.tmaid.2019.02.004

27. Parker DM, Landier J, Thu AM, Lwin KM, Delmas G, Nosten FH, et al. Scale up of a Plasmodium falciparum elimination program and surveillance system in Kayin State, Myanmar. Wellcome open research. 2017;2. doi: 10.12688/wellcomeopenres.12741.2 29384151

28. Landier J, Parker DM, Thu AM, Lwin KM, Delmas G, Nosten FH, et al. Effect of generalised access to early diagnosis and treatment and targeted mass drug administration on Plasmodium falciparum malaria in Eastern Myanmar: an observational study of a regional elimination programme. The Lancet. 2018;391(10133):1916–1926. doi: 10.1016/S0140-6736(18)30792-X 29703425

29. Blanton RE. Population genetics and molecular epidemiology of eukaryotes. Microbiology spectrum. 2018;6(6). doi: 10.1128/microbiolspec.AME-0002-2018 30387414

30. Wesolowski A, Taylor AR, Chang HH, Verity R, Tessema S, Bailey JA, et al. Mapping malaria by combining parasite genomic and epidemiologic data. BMC medicine. 2018;16(1):1–8. doi: 10.1186/s12916-018-1232-2

31. Gao B, Saralamba S, Lubell Y, White LJ, Dondorp AM, Aguas R. Determinants of MDA impact and designing MDAs towards malaria elimination. Elife. 2020;9:e51773. doi: 10.7554/eLife.51773

32. Taylor AR, Jacob PE, Neafsey DE, Buckee CO. Estimating relatedness between malaria parasites. Genetics. 2019;212(4):1337–1351. doi: 10.1534/genetics.119.302120

33. Neafsey DE, Schaffner SF, Volkman SK, Park D, Montgomery P, Milner DA, et al. Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence. Genome biology. 2008;9(12):R171. doi: 10.1186/gb-2008-9-12-r171 19077304

34. Peyré G, Cuturi M, et al. Computational optimal transport. Foundations and Trends® in Machine Learning. 2019;11(5-6):355–607. doi: 10.1561/2200000073

35. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–959.

36. Watson JA, Taylor AR, Ashley EA, Dondorp AM, Buckee CO, White NJ, et al. Pre-print: A cautionary note on the use of machine learning algorithms to characterise malaria parasite population structure from genetic distance matrices. bioRxiv. 2020; p. 1–18.

37. Wille M, Holmes EC. The Ecology and Evolution of Influenza Viruses. Cold Spring Harbor Perspectives in Medicine. 2019; p. a038489.

38. Katz EM, Esona MD, Betrapally NS, Lucia A, Neira YR, Rey GJ, et al. Whole-gene analysis of inter-genogroup reassortant rotaviruses from the Dominican Republic: Emergence of equine-like G3 strains and evidence of their reassortment with locally-circulating strains. Virology. 2019;534:114–131. doi: 10.1016/j.virol.2019.06.007 31228725

39. Caugant DA, Brynildsrud OB. Neisseria meningitidis: using genomics to understand diversity, evolution and pathogenesis. Nature Reviews Microbiology. 2019; p. 1–13.

40. Smith JM, Feil EJ, Smith NH. Population structure and evolutionary dynamics of pathogenic bacteria. Bioessays. 2000;22(12):1115–1122. doi: 10.1002/1521-1878(200012)22:12%3C1115::AID-BIES9%3E3.0.CO;2-R

41. Tibayrenc M, Ayala FJ. The clonal theory of parasitic protozoa: 12 years on. Trends in parasitology. 2002;18(9):405–410. doi: 10.1016/S1471-4922(02)02357-7

42. Rajendran C, Su C, Dubey JP. Molecular genotyping of Toxoplasma gondii from Central and South America revealed high diversity within and between populations. Infection, Genetics and Evolution. 2012;12(2):359–368. doi: 10.1016/j.meegid.2011.12.010

43. Nader JL, Mathers TC, Ward BJ, Pachebat JA, Swain MT, Robinson G, et al. Evolutionary genomics of anthroponosis in Cryptosporidium. Nature microbiology. 2019;4(5):826–836. doi: 10.1038/s41564-019-0377-x 30833731

44. Nieuwenhuis BP, James TY. The frequency of sex in fungi. Philosophical Transactions of the Royal Society B: Biological Sciences. 2016;371(1706):20150540. doi: 10.1098/rstb.2015.0540

45. Fruchterman TM, Reingold EM. Graph drawing by force-directed placement. Software: Practice and experience. 1991;21(11):1129–1164.

46. Sáenz FE, Morton LC, Okoth SA, Valenzuela G, Vera-Arias CA, Vélez-Álvarez E, et al. Clonal population expansion in an outbreak of Plasmodium falciparum on the northwest coast of Ecuador. Malaria journal. 2015;14(1):497. doi: 10.1186/s12936-015-1019-2

47. Vera-Arias CA, Castro LE, Gómez-Obando J, Sáenz FE. Diverse origin of Plasmodium falciparum in northwest Ecuador. Malaria journal. 2019;18(1):251. doi: 10.1186/s12936-019-2891-y

48. Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164(4):1567–1587.

49. Huestis DL, Dao A, Diallo M, Sanogo ZL, Samake D, Yaro AS, et al. Windborne long-distance migration of malaria mosquitoes in the Sahel. Nature. 2019;574(7778):404–408. doi: 10.1038/s41586-019-1622-4 31578527

50. Verdonschot PF, Besse-Lototskaya AA. Flight distance of mosquitoes (Culicidae): a metadata analysis to support the management of barrier zones around rewetted and newly constructed wetlands. Limnologica-Ecology and Management of Inland Waters. 2014;45:69–79. doi: 10.1016/j.limno.2013.11.002

51. Pacheco MA, Schneider KA, Céspedes N, Herrera S, Arévalo-Herrera M, Escalante AA. Limited differentiation among Plasmodium vivax populations from the northwest and to the south Pacific Coast of Colombia: A malaria corridor? PLoS neglected tropical diseases. 2019;13(3):e0007310. doi: 10.1371/journal.pntd.0007310

52. Guagliardo SA, Morrison AC, Barboza JL, Requena E, Astete H, Vazquez-Prokopec G, et al. River boats contribute to the regional spread of the dengue vector Aedes aegypti in the Peruvian Amazon. PLoS neglected tropical diseases. 2015;9(4):e0003648. doi: 10.1371/journal.pntd.0003648 25860352

53. Lounibos LP. Invasions by insect vectors of human disease. Annual review of entomology. 2002;47(1):233–266. doi: 10.1146/annurev.ento.47.091201.145206

54. Montoya-Lerma J, Solarte YA, Giraldo-Calderón GI, Quiñones ML, Ruiz-López F, Wilkerson RC, et al. Malaria vector species in Colombia: a review. Memórias do Instituto Oswaldo Cruz. 2011;106:223–238. doi: 10.1590/S0074-02762011000900028 21881778

55. Gutiérrez LA, Naranjo NJ, Cienfuegos AV, Muskus CE, Luckhart S, Conn JE, et al. Population structure analyses and demographic history of the malaria vector Anopheles albimanus from the Caribbean and the Pacific regions of Colombia. Malaria journal. 2009;8(1):259. doi: 10.1186/1475-2875-8-259 19922672

56. Instituto Nacional de Salud Colombia, Dirección de Vigilancia y Analisis del Riesgo en Salud Pública. Boletín Epidemiológico Semanal: semana epidemiológica 52; 2017. Available from: https://www.ins.gov.co/buscador-eventos/Paginas/Vista-Boletin-Epidemilogico.aspx.

57. Instituto Nacional de Salud Colombia, Dirección de Vigilancia y Analisis del Riesgo en Salud Pública. Boletín Epidemiológico Semanal: semana epidemiológica 52; 2018. Available from: https://www.ins.gov.co/buscador-eventos/Paginas/Vista-Boletin-Epidemilogico.aspx.

58. Instituto Nacional de Salud Colombia, Dirección de Vigilancia y Analisis del Riesgo en Salud Pública. Boletín Epidemiológico Semanal: semana epidemiológica 52; 2019. Available from: https://www.ins.gov.co/buscador-eventos/Paginas/Vista-Boletin-Epidemilogico.aspx.

59. Wabgou M, Vargas D, Carabalí JA. Las migraciones internacionales en Colombia. Investigación & Desarrollo. 2012;20(1):142–167.

60. Ahumada ML, Orjuela LI, Pareja PX, Conde M, Cabarcas DM, Cubillos EFG, et al. Spatial distributions of Anopheles species in relation to malaria incidence at 70 localities in the highly endemic Northwest and South Pacific coast regions of Colombia. Malaria Journal. 2016;15(407):1–16. 27515166

61. Shetty AC, Jacob CG, Huang F, Li Y, Agrawal S, Saunders DL, et al. Genomic structure and diversity of Plasmodium falciparum in Southeast Asia reveal recent parasite migration patterns. Nature communications. 2019;10(1):2665. doi: 10.1038/s41467-019-10121-3 31209259

62. Naranjo-Díaz N, Altamiranda M, Luckhart S, Conn JE, Correa MM. Malaria vectors in ecologically heterogeneous localities of the Colombian Pacific region. PLoS One. 2014;9(8):e103769. doi: 10.1371/journal.pone.0103769

63. Csardi G, Nepusz T. The igraph software package for complex network research. InterJournal, complex systems. 2006;1695(5):1–9.

64. R Core Team. R: A Language and Environment for Statistical Computing; 2018. Available from: https://www.R-project.org/.

65. Schuhmacher D, Bähre B, Gottschlich C, Hartmann V, Heinemann F, Schmitzer B. transport: Computation of Optimal Transport Plans and Wasserstein Distances; 2019. Available from: https://cran.r-project.org/package=transport.


Článek vyšel v časopise

PLOS Genetics


2020 Číslo 11
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#