Document Detail

The emergence and maintenance of sickle cell hotspots in the Mediterranean.
Jump to Full Text
MedLine Citation:
PMID:  22704979     Owner:  NLM     Status:  MEDLINE    
Genetic disorders of haemoglobin (haemoglobinopathies), including the thalassaemias and sickle cell anaemia, abound in historically malarious regions, due to the protection they provide against death from severe malaria. Despite the overall spatial correlation between malaria and these disorders, inter-population differences exist in the precise combinations of haemoglobinopathies observed. Greece and Italy present a particularly interesting case study: their high frequencies of beta thalassaemia speak to a history of intense malaria selection, yet they possess very little of the strongly malaria protective mutation responsible for sickle cell anaemia, despite historical migrational links with Africa where high frequencies of sickle cell occur. Twentieth century surveys of beta thalassaemia and sickle cell in Greece, Sicily and Sardinia have revealed striking sickle cell 'hotspots' - places where the frequency of sickle cell approaches that seen in Africa while neighbouring populations remain relatively sickle cell free. It remains unclear how these hotspots have been maintained over time without sickle cell spreading throughout the region. Here we use a metapopulation model to show that (i) epistasis between the alpha and beta forms of thalassaemia can restrict the spread of sickle cell through a network of linked subpopulations and (ii) the emergence of sickle cell hotspots requires relatively low levels of gene flow, but the aforementioned epistasis increases the chances of hotspots forming.
Bridget S Penman; Sunetra Gupta; Caroline O Buckee
Related Documents :
22659329 - Cisplatin resistance: a cellular self-defense mechanism resulting from multiple epigene...
21539399 - Polydopamine/liposome coatings and their interaction with myoblast cells.
3779849 - Multiple forms of acid phosphatase and their differential secretion during stalk format...
7493629 - Motility of vinculin-deficient f9 embryonic carcinoma cells analyzed by video, laser co...
6738529 - Overproduction of discoidin i by a temperature-sensitive motility mutant of dictyosteli...
7417709 - Polymorphonuclear motility: measurement by computer-linked image analysis.
7720629 - Ontogeny of growth hormone (gh)-secreting cells during chicken embryonic development: i...
16387649 - Coherence and timing of cell cycle start examined at single-cell resolution.
10545239 - Spontaneous calcium transients are required for neuronal differentiation of murine neur...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2012-06-13
Journal Detail:
Title:  Infection, genetics and evolution : journal of molecular epidemiology and evolutionary genetics in infectious diseases     Volume:  12     ISSN:  1567-7257     ISO Abbreviation:  Infect. Genet. Evol.     Publication Date:  2012 Oct 
Date Detail:
Created Date:  2012-08-27     Completed Date:  2012-12-26     Revised Date:  2013-07-12    
Medline Journal Info:
Nlm Unique ID:  101084138     Medline TA:  Infect Genet Evol     Country:  Netherlands    
Other Details:
Languages:  eng     Pagination:  1543-50     Citation Subset:  IM    
Copyright Information:
Copyright © 2012 Elsevier B.V. All rights reserved.
Department of Zoology, University of Oxford, Oxford OX1 3PS, UK.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Computer Simulation
Endemic Diseases
Epistasis, Genetic
Gene Flow
Gene Frequency
Greece / epidemiology
Italy / epidemiology
Mediterranean Islands / epidemiology
Mediterranean Region / epidemiology
Models, Genetic
alpha-Thalassemia / epidemiology,  genetics*
beta-Thalassemia / epidemiology,  genetics*

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Infect Genet Evol
Journal ID (iso-abbrev): Infect. Genet. Evol
ISSN: 1567-1348
ISSN: 1567-7257
Publisher: Elsevier Science
Article Information
© 2012 Elsevier B.V.
Received Day: 20 Month: 3 Year: 2012
Revision Received Day: 1 Month: 6 Year: 2012
Accepted Day: 5 Month: 6 Year: 2012
pmc-release publication date: Month: 10 Year: 2012
Print publication date: Month: 10 Year: 2012
Volume: 12 Issue: 7
First Page: 1543 Last Page: 1550
ID: 3438445
PubMed Id: 22704979
Publisher Id: MEEGID1314
DOI: 10.1016/j.meegid.2012.06.001

The emergence and maintenance of sickle cell hotspots in the Mediterranean
Bridget S. Penmana Email:
Sunetra Guptaa
Caroline O. Buckeeb
aDepartment of Zoology, University of Oxford, Oxford OX1 3PS, UK
bCenter for Communicable Disease Dynamics, Department of Epidemiology, Harvard School of Public Health, Boston, MA 02115, USA
Corresponding author. Tel.: +44 01865 281309.


The malaria parasite has imposed arguably the strongest evolutionary pressure of any pathogen on human populations. In 1949, JBS Haldane suggested that this selection pressure may have resulted in a high frequency of genetic blood disorders in malaria-endemic populations (Haldane, 1949). This ‘malaria hypothesis’ is now widely accepted to apply to a variety of mutations affecting the alpha or beta globin genes that encode the subunits of haemoglobin (haemoglobinopathies), including those responsible for sickle cell anaemia (Allison, 1954; Williams et al., 2005b); haemoglobin C (Modiano et al., 2001); haemoglobin E (Chotivanich et al., 2002) and the alpha and beta thalassaemias (Hill et al., 1988; Flint et al., 1986). A recent review of the genetics of malaria resistance is provided by Hedrick (2011); Taylor et al. (2012) have also carried out a meta-analysis of the evidence for malaria protection provided by each globin mutation.

Almost every old-world malarious region hosts a range of haemoglobinopathies, but – as reviewed in Flint et al. (1998) – the suite of haemoglobinopathies observed varies between populations. Extensive molecular and epidemiological studies, motivated by the enormous public health significance of these disorders, are starting to elucidate how mutations in the alpha and beta globin genes may interact with one another to affect clinical phenotypes (Kan and Nathan, 1970; Wainscoat et al., 1983a,b; Galanello et al., 1989; Williams et al., 2005a). The haemoglobinopathies thus represent a unique opportunity to compare the relative importance of natural selection, migration and gene-gene interactions (epistasis) in generating spatial genetic variation in human populations.

It has been suggested that heterogeneity in the geographical distribution of haemoglobinopathies results from the relatively recent emergence of strong malaria selection at some point within the last 8000 years, following the increase in human population density associated with the advent of agriculture (Flint et al., 1998). Under this scenario, there simply has not been enough time for population mixing to produce a uniform pattern of mutations. However, genetic spatial heterogeneity can also arise and be maintained in populations through a variety of other evolutionary and ecological processes, (Buckee et al., 2004; Hedrick, 2006; Buckee et al., 2007; Lawton-Rauh, 2008; Cushman et al., 2011; Neto et al., 2011; Lappalainen et al., 2011). We have previously shown that epistatic interactions among the haemoglobinopathies may account for large-scale geographic differences in the relative frequencies of haemoglobinopathies in African and Mediterranean populations (Penman et al., 2009). In Africa, high frequencies of sickle cell (βS) co-exist with frequencies of alpha thalassaemia that do not exceed 50%, whereas in the Mediterranean, diverse alpha and beta thalassaemic mutations are present, and βS is largely absent. We proposed that these contrasting suites of alleles could be maintained by two well-documented genetic interactions: (i) positive epistasis, where coinheriting alpha and beta thalassaemia leads to a less severe blood disorder than that caused by beta thalassaemia alone (Kan and Nathan, 1970; Wainscoat et al., 1983a,b; Galanello et al., 1989; Weatherall and Clegg, 2001a) and (ii) negative epistasis, where coinheriting alpha thalassaemia and sickle cell trait leads to a loss of malaria protection (Williams et al., 2005a; May et al., 2007). The former can allow relatively low frequencies of the thalassaemias to exclude βS from a population (Penman et al., 2009); the latter limits the frequency of alpha thalassaemia when βS is present (Williams et al., 2005a), or could even allow very high levels of alpha thalassaemia to exclude βS (Penman et al., 2011).

Although this framework provides a parsimonious explanation for the broad qualitative patterns observed across the two continents, striking local heterogeneities occur in the distribution of blood disorders within populations, which cannot be accounted for within a deterministic framework. Fig. 1 collates the results of five geographical surveys that cover Greece, Sicily and Sardinia (Barnicot et al., 1963; Stamatoyannopoulos and Fessas, 1964; Siniscalco et al., 1966; Schiliro et al., 1986; Cao et al., 2008), and reveals several interesting features: beta thalassaemia is ubiquitous and βS is limited in its spatial distribution (as previously observed), yet in a number of ‘hotspots’ the βS frequency observed in Greek populations is extremely high (>0.1) and comparable to that seen in Africa. It is unclear how these hotspots have been maintained, and why βS has not spread further across the Mediterranean given (i) the history of strong malaria selection in Greece and southern Italy; (ii) the extremely high degree of malaria protection βS provides, and (iii) the close migrational ties between Africa, Greece, and Italy. Here, we explore these local patterns of genetic variation in the Mediterranean, and develop a metapopulation model to investigate how processes such as migration and selection impact the local heterogeneity of blood disorders. This follows the successful use of linear metapopulation (‘stepping stone’) population genetic models by Livingstone (Livingstone 1969a, 1976, 1989) to explore the distribution of malaria protective beta globin variants, both in West Africa and across the old world. We confirm that epistatic interactions between the thalassaemias can help keep βS frequencies low throughout most of a population. We also show that hotspots of sickle such as those observed in current Mediterranean populations require low levels of gene flow between the hotspot region and the rest of the population if they are to be maintained, but that epistasis increases the range of gene flows at which hotspots are possible. It is widely believed that βS came to the Mediterranean from Africa. Our results suggest that βS could have first been introduced to the Mediterranean >2000 years before the present.


We limit our attention to the three most common haemoglobinopathies in the Mediterranean region (Weatherall and Clegg, 2001b): alpha thalassaemia caused by a single deletion of an alpha globin gene (represented as α+); beta thalassaemia caused by a mutation which significantly reduces but does not completely eliminate beta globin production (β+), and the point mutation responsible for sickle cell anaemia (βS). This scheme captures the most important features of Mediterranean haemoglobinopathies, but necessarily simplifies the situation at the beta globin locus by representing multiple beta thalassaemic mutations with a single variant (see Supplementary material). We consider a metapopulation consisting of 40 subpopulations (demes) linked via migration.

Migration between demes was simulated every generation by letting a certain proportion of each deme (m) be assigned allele frequencies equal to the average allele frequencies across all the demes connected to it. Connections between demes followed a predefined network structure (see Supplementary Methods and Watts and Strogatz, 1998). m represents ‘gene flow’ in the sense of the proportion of alleles in a population which originated outside of that population, and the chosen value of m determines whether the demes represent a small villages or settlements, with strong ties to other villages in the same region, or larger collections of settlements in different regions, with weaker migrational ties.

The frequency of each gene of interest within each deme was calculated using the following standard population genetic equation (see Mandel, 1958) relating the frequency of allele i (pi) in generation t + 1 to the frequency of the allele i in generation t:

[Formula ID: e0005]
where i and j represent the frequencies of alleles at the α-globin locus and take values 1 or 2, and k and l represent the frequencies of allele at the β-globin locus and take values 1,2, or 3. aijkl represents the fitness of genotype ‘ijkl’. We have chosen to ignore the effects of genetic drift, since all the alleles we are concerned with are under strong selection from malaria. α-globin and β-globin are considered completely unlinked, since the alpha and beta globin clusters occur on different chromosomes in humans.

In order to investigate the effects of epistasis, each genotype was allocated two mortality rates: one associated with the severity of its blood disorder, and one associated with the level of malaria protection it experiences (see Table 1). We converted these mortality rates into a measure of relative fitness (see Supplementary material), and considered two different scenarios: one where positive epistasis between alpha and beta thalassaemia was included, and one where it was not. Negative epistasis between alpha thalassaemia and sickle cell trait was always included, since we have strong evidence that it occurs (Williams et al., 2005a; May et al., 2007), and we wanted to assess the impact of positive epistasis explicitly. Whenever we refer to ‘epistasis’ in the results and discussion, we are referring to positive epistasis between alpha and beta thalassaemia.

At the start of each simulation, all demes were assumed to contain a set starting frequency (T) of α+, and the same starting frequency of β+. We ran two different versions of the model – in the first, we assumed that βS was first introduced 100 generations ago (assuming a generation time of 20 years, this is approximately 2000 years ago), at a frequency of 0.001, into a population containing a fixed proportion of the thalassaemias. For the second scenario, we considered 200 generations (4000 years) of malaria selection, and investigated what happened when βS was first introduced at various time points. After the first introduction of βS in any simulation, we allowed it to arrive at a frequency of 0.001 in a randomly chosen subset of subsequent generations. The first introduction of βS was always into deme 20, and subsequent introductions were always into either 20 or 40, which represent trading posts, or other entry points into the region.

Reported beta thalassaemia frequencies from Mediterranean communities range between 0.014 and 0.19 (Fig. 1). Alpha thalassaemia frequencies are less well-known (see further discussion in the Supplementary material). The studies which made up Fig. 1 have been conducted at different scales: villages, sets of villages, or screening centres serving whole regions – so we do not claim to have a clear idea of the true frequency at a particular spatial scale, and nor do we seek to exactly recapitulate the history of the Mediterranean in our simulations. Nevertheless, we have applied some constraints in an attempt to keep our simulations realistic: we aimed to keep ‘modern day’ beta thalassaemic frequencies below 0.19 and ‘modern day’ alpha thalassaemia frequencies below 0.3 (see Supplementary material). This limited the ranges of starting thalassaemia frequencies we tested.

For the sake of simplicity, we have only included one element of stochasticity in this model: the randomly chosen generations in which sickle cell re-challenges the population. Supplementary Fig. S7 illustrates how results vary over 30 repeats, with different sets of parameters. We considered that this minimum of 30 repeats gave a reliable sense of the emerging patterns. 100 repeats were used to generate each cell of a heatmap.

3.1  Epistasis can restrict the spread of βS , despite extensive gene flow

Figs. 2a and b illustrate how gene flow (x axis) and the frequency of thalassaemia present when βS arrives (y axis) affect the spread of βS through the population. So long as the thalassaemias are above a certain threshold when βS arrives (0.08 in this case), the inclusion of positive epistasis between alpha and beta thalassaemia can limit its distribution over the next 100 generations. This result echoes our 2009 paper (Penman et al., 2009), but extends our general observation – that a threshold level of epistatically interacting thalassaemias can restrict βS – to a network of linked subpopulations, some of which are repeatedly challenged by βS.

In our 2009 paper, complete competitive exclusion of βS by beta thalassaemia was possible, and the ‘winning’ allele was determined by: (i) the relative fitness of beta thalassaemia (which is inherently bound up with its blood disorder severity, malaria protectiveness, and whether or not epistasis with alpha thalassaemia is present) and (ii) the frequency of beta thalassaemia when sickle cell arrives (see Supplementary Fig. S4 and Penman et al., 2009). Within the metapopulation presented here, complete exclusion of βS is not possible during the 100 generations we have simulated (the equilibrium point described in Penman et al., 2009 is not reached). However, the average frequency of βS in the population after 100 generations is lower in the presence of epistasis than in its absence, and epistasis renders the entire population less susceptible to the effects of increasing the rate of gene flow (Fig. 2 and Supplementary Fig. S5).

Contrary to our expectations, the precise network structure underlying patterns of migration between demes was not a primary determinant of the spread of sickle cell (Supplementary Fig. S7). Whilst in certain circumstances an effect of network structure is discernible (discussed in the supplementary material), its impact was minor in comparison to the starting frequency of thalassaemia, or the level of gene flow between the demes.

How does a 100 generation simulation play out with the starting thalassaemia allele frequencies of Fig. 2? As shown in Supplementary Fig. S6, alpha thalassaemia frequencies increase over the 100 generations, reaching a ‘present day’ frequency of ∼0.28. Beta thalassaemia frequencies on the other hand, actually decline from their starting frequency of 0.08 (although as can be seen in Supplementary Fig. S4, under epistasis they would eventually recover). We can conclude that (i) a thalassaemia frequency of 0.08 100 generations ago need not be associated with unrealistic beta thalassaemia frequencies today, but also that (ii) for βS to be kept out of a population, some earlier process (such as long term malaria selection, or simply the rapid spread of thalassaemic alleles in small populations) must have led to high thalassaemia levels before βS arrived.

3.2  The likelihood of “sickle hotspots” is increased by epistasis but only under low levels of gene flow

As discussed in the introduction, the distribution of βS in the Mediterranean is highly discontinuous. Based on the observations in Fig. 1, where the Chalkidhiki peninsular in Greece represents a ‘spike’ of βS surrounded by an environment of lower frequencies, we defined a ‘βS hotspot’ as a deme with a βS frequency >0.06, where the median frequency of βS in demes with which it has a direct connection via migration is <0.02. We chose this definition so as to capture the idea of βS being contained within a particular region.

Fig. 3 illustrates the likelihood of obtaining such hotspots with and without epistasis, for different thalassaemia starting frequencies; levels of gene flow and levels of malaria selection. Two results are apparent: (i) hotspot demes are only possible at relatively low levels of interdeme gene flow, but (ii) the inclusion of epistasis and a threshold starting frequency of thalassaemia increases the parameter space where hotspots are possible. Supplementary Figs. S7 and S8 provide a fuller exploration of the effects of network structure and sickle cell challenge on the formation of hotspots. Once again, network structure appears to be one of the least important factors.

The likelihood of hotspots in the scenario captured in Fig. 3 is maximised when epistasis is present and the starting frequency of alpha and beta thalassaemia is 0.08 (Fig. 3b). This result can be understood when viewed in conjunction with Fig. 2b: at lower thalassaemia starting frequencies, βS will occupy a large proportion of the entire network, but at higher thalassaemia starting frequencies, the thalassaemias are so effective at keeping βS out that not even a small hotspot can form.

To explore the formation of hotspots further, we developed the pared-down network presented in Supplementary Fig. S9. Two communities with high internal inter-deme gene flow are linked by a connection of much lower gene flow. βS is introduced at a high frequency into one deme in community one (simulating the arrival, perhaps, of a population transplanted from Africa). By considering the frequencies of βS in the two communities, Supplementary Fig. S9 lets us move away from our somewhat arbitrary definition of a hotspot above, and see that for any chosen ratio of βS frequencies in the two communities (e.g. a ratio of 10:1), the inclusion of epistasis allows that ratio to occur at a higher level of gene flow between the communities.

3.3  Epistasis allows for a longer history of βS challenge in the Mediterranean

We next explored the temporal scales underlying sickle cell hotspot formation and maintenance in the region. As described in the Section 2, we considered a 200 generation long scenario, in which we tested the outcomes when βS was first introduced at different time points. Fig. 4a matches the time of entry of βS to approximate present day Sardinian and Greek/Sicilian patterns, for one such timeline. Details of the criteria for each present day pattern are provided in the figure legend. The inclusion of epistasis increases the range of times for the first introduction of βS that are consistent with a present day Sardinian scenario. The Greek/Sicilian pattern is associated with a relatively narrow window for the first introduction of βS, but the presence of epistasis allows this window to occur earlier – implying a longer history of βS challenge.

The small window of potential first entry times for βS associated with the Greek/Sicilian pattern reflects the transient nature of the hotspots themselves: enough time must pass so that βS attains hotspot frequencies, but not so much time that βS can spread too much into the surrounding demes.

Figs. 4b and c offer snapshots of different stages of the simulation, for a scenario where the first entry of sickle cell was 110 generations ago. In the presence of epistasis, present day high frequencies of βS are limited to its points of entry and their immediate contacts; without epistasis, by contrast βS has spread further into the secondary contacts of the original points of entry. If we allowed these populations to continue to evolve over another 100 generations of malaria selection, we arrive at the scenarios shown at the right hand side. Without epistasis, βS dominates the population, and will eventually take over completely: genetic heterogeneity is thus a short-lived phenomenon. In the presence of epistasis, however, hotspots can be maintained indefinitely. Gene-gene interactions therefore provide a mechanism for generating genetic heterogeneity in populations that is stable in time and space.


Our model illustrates that a threshold level of epistatically interacting thalassaemias can (i) keep βS frequencies very low in most demes in a population, and (ii) extend the circumstances under which βS hotspots are possible, within the overall constraint that hotspots can only emerge in relatively self-contained communities with low levels of gene flow to the rest of the population. We propose that a combination of these two effects could have allowed βS to challenge the Mediterranean for over two millennia, yet still be restricted to a few isolated regions today. These results confirm our previous large-scale model (Penman et al., 2009), and provide a spatially explicit framework incorporating migration, local heterogeneity and plausible time frames.

Published data on the distribution of haemoglobinopathies in the Mediterranean can be resolved into two main patterns: Greece and Sicily, where βS is present and attains pockets of high frequency, and Sardinia, where βS is almost completely absent. Explaining the Sardinian pattern is straightforward: Fig. 2 illustrates that, in the presence of positive epistasis between alpha and beta thalassaemia, a threshold level of the thalassaemias can prevent sickle cell from taking over a population – regardless of the level of gene flow between demes within that population. Explaining the hotspots of the Greek/Sicilian population is more challenging: hotspot demes only form within a relatively constrained region of parameter space. We saw in Fig. 4 that the window for the first introduction of βS associated with the Greek/Sicilian scenario is small both with and without epistasis, but that epistasis causes it to occur earlier. The interaction between alpha and beta thalassaemia slows the rise in frequency of βS, delaying both the time taken to reach hotspot levels and its rate of spread out of the introduction deme.

We have illustrated the behaviour of the model at various rates of gene flow, but what is a reasonable degree of mixing to assume between different demes within a metapopulation? Adams and Kasakoff (1976) explored the relationship between endogamy and population size. They found that sub populations in most societies could be described in terms of ‘80% endogamous groups’: a tribe in a pastoralist society; a valley in the highlands of New Guinea; a set of small villages in a peasant society. Levels of endogamy higher than 80% were typically associated with whole ethnic groups. When the gene flow level in our model (m) is set at 0.1, our demes are consistent with 80% endogamous groups (assuming 20% of marriages bring in a partner from outside the group, contributing 10% of the genetic material to the next generation). We found that gene flows of 1% or lower were necessary to generate hotspot demes. This implies that the type of deme which can exist as a sickle cell hotspot cannot be a single 80% endogamous local population: instead it may be a relatively closed collection of such 80% groups, with geographical or social barriers that limit genetic exchange with other populations.

Metapopulation approaches are widely used in ecology to model the dynamics of populations in fragmented or spatially heterogeneous habitats (reviewed in Hanski and Gilpin, 1991), and were first applied to the haemoglobinopathies by Livingstone in the form of a ‘stepping stone’ population genetic model. In such a model, demes are ranged along a line, most migration occurs between immediate or secondary neighbours, but there is also the possibility of random long-distance migration (Livingstone 1969a, 1976, 1989). Livingstone considered the generation and maintenance of clines in the frequencies of βS in West Africa and beta thalassaemia in Sardinia (Livingstone, 1969a); the overlapping clines of βS, beta thalassaemia and haemoglobin C in West Africa (Livingstone, 1976), and the interaction between βS, beta thalassaemia, haemoglobin C and haemoglobin E on an inter-continental scale (Livingstone, 1989). This series of papers demonstrated that the global distribution of the malaria protective haemoglobinopathies is a product of human migration and gene flow as well as malaria selection. In response to Livingstone’s work, Fix (see review: Fix, 2003 and chapter 4 of Fix, 1999) considered other factors which could lead to the rapid accumulation and spread of malaria protective haemoglobin variants, including kin-structured migration (where families were assumed to migrate together, concentrating the import of particular alleles) and a group selection model in which demes with insufficient frequencies of malaria protective alleles risked becoming extinct.

Both Fix’s and Livingstone’s work aimed to account for present day allele frequencies and clines within realistic historical timescales. Enhancing the spread of malaria protective variants (by assuming long distance gene flow or the other processes mentioned above) was necessary in order that present day patterns could be achieved within the hundreds of generations that malaria is thought to have exacted significant mortality on human populations. The work we have presented here is concerned with the opposite problem: how could a high frequency of βS have accumulated in a population under malaria selection, yet not spread over the course of ∼100 generations? Our model differs from those of Livingstone and Fix in that we considered alleles at both the alpha and beta globin loci, not beta globin alone, and we sought to understand the role of epistasis in creating present day patterns. Epistasis affects the rate of spread of sickle cell, but it can only help to generate hotspots if gene flow in and out of the hotspot population is already low.

Unlike Livingstone’s models, Fix developed models considering whole numbers of individuals. This modification made it easier to consider both kin structured migration and genetic drift. Livingstone also used a model of a finite yet growing population to demonstrate that a founder effect coupled with rapid population expansion could retain the βS gene even in the absence of malaria selection (Livingstone, 1969b). In future efforts to understand the micro-heterogeneity of the haemoglobinopathies in the Mediterranean, spatially explicit individual based models with greater stochasticity and more demographic detail than the framework presented here are likely to prove extremely informative.

Fig. 4c shows that the presence of epistasis can contain the spread of sickle cell indefinitely. However, the population in the third panel of 4c has alpha thalassaemia frequencies of ∼0.57. Since such frequencies are higher than those reported in the twentieth century Mediterranean, we did not consider it a realistic present day outcome. The relatively low frequencies of alpha thalassaemia reported in the Mediterranean are discussed further in the supplementary material; they are a significant constraint on the outcomes possible within the model, and it is unfortunate that we do not have estimates of alpha thalassaemia frequencies at a comparable geographical scale to those for beta thalassaemia. Three of the studies in Fig. 1 (Barnicot et al., 1963; Stamatoyannopoulos and Fessas, 1964; Siniscalco et al., 1966) commented on the patterns of βS and beta thalassaemia in the Mediterranean, and suggested that these two beta globin variants may act to exclude one another from populations. Our model is consistent with the conclusion that βS in the Mediterranean has been restricted by evolutionary interactions between genes – however, we have shown that a combination of both alpha and beta thalassaemia can be more effective at excluding βS than beta thalassaemia alone. Whether or not a beta thalassaemic allele can out-compete βS on its own depends on how much malaria protection it offers and how severe a blood disorder it causes. In Supplementary Fig. S4 we illustrate the range of hypothetical beta thalassaemic alleles that can outcompete βS, and show how this range is extended by the inclusion of epistasis with alpha thalassaemia. It remains a theoretical possibility that the beta thalassaemia found in the Mediterranean has the necessary properties to outcompete βS alone, but we have not considered this possibility here, since we wished to focus on the potential population genetic impact of epistasis.

A recent survey of βS and beta thalassaemia within the Chalkidhiki peninsular (Kalleas et al., 2012) provides more detail about the modern day pattern of beta globin variants in one of the Greek βS hotspots. It is striking that the distribution of βS is highly patchy within the peninsular (see Fig. 1 of Kalleas et al.), whilst thalassaemia is present even in (historically non-malarious) mountainous regions – presumably a signature of gene flow within the peninsular. The most intriguing observation in this study relates to the beta thalassaemic mutations present in the Chalkidhiki peninsular alongside βS. Throughout most of the Greek mainland, the IVS-I-110 beta thalassaemic mutation is more common than the codon 39 beta thalassaemic mutation, but in Chalkidhiki the pattern is reversed. Kalleas et al. note that there is a second region (in Central Greece) where codon 39 exceeds IVS-I-110, and where βS is also present. However, codon 39 also predominates in Sardinia, which is βS free. Kalleas et al. point out that the relative geographical isolation of all of these regions may have contributed to their particular genetic patterns. Nevertheless, it is also tempting to speculate that the suite of beta thalassaemic mutations predominating in a particular region may have made particular populations more or less susceptible to sickle cell invasion. For this paper, we represented all beta thalassaemia in the Mediterranean with a single, relatively severe β+ mutation (see Section 2 and Supplementary material). In reality, codon 39 (which, as we just noted, is very common in Sardinia) is a more severe mutation than IVS-I-110 (see Supplementary material and Weatherall and Clegg, 2001a). Allowing for the diversity of beta thalassaemia in the Mediterranean in future modeling work may shed more light on the βS question.

Sickle cell anaemia in Greece is caused by the Benin haplotype (Boussiou et al., 1991), one of several different chromosomal backgrounds on which the βS mutation is found. This haplotype is generally believed to have originated in central West Africa, but the timing of its evolution and spread is unclear. Pagnier et al. (1984) speculate that the Benin haplotype spread to North Africa along trans-Saharan trade routes via camel caravans from the third century CE onwards. Depending on when and where βS emerged in Africa, however, it could have been imported into Greece and Sicily via Carthaginian slave traders from 650 BCE; Ancient Greek and Roman settlements in North Africa, or the later Arab conquest of Sicily by the ninth Century CE (dates from the Penguin Atlas of World History by Kinder and Hilgemann – see especially pages 39 and 119). The earliest evidence of βS from a population known to be in contact with the Ancient Mediterranean comes from a skeleton containing fossilized sickle erythrocytes, excavated from an island in the Persian Gulf that was colonised by Alexander the Great. This skeleton was dated from 130+/−80 years BCE (Salares et al., 2004).

Are there any specific historical reasons for the contrasting patterns observed in Sardinia, Greece and Sicily? Siniscalco et al. (1966) note the relative isolation of the Sardinian interior: “Ancient Greek colonisation was, in fact, limited to Olbia (the northern portion of the island), while the Romans and Carthaginians only exploited the coastal regions for grain”. Such isolation is in keeping with the first introduction of βS to Sardinia occurring relatively late – although, importantly, we do not require that Sardinia was never exposed to βS. In the case of Sicily, Schiliro et al. (1986) point out that most βS on the island occurs in coastal regions that were the sites of Ancient Greek colonisation, suggesting that βS came to Sicily via Greece. The founders of the ancient Sicilian city states of Syracruse and Gela came from Corinth, Rhodes and Crete – none of which are present day Greek sickle cell hotspots – but it is entirely possible that the Greeks precipitated the spread of sickle cell in their Sicilian settlements by bringing in African slaves. As for Greece itself, Stamatoyannopoulos et al. do not posit any specific historical events as the reason for the Greek pattern of βS, but they do note that βS seems only to be found in regions of high malaria selection. However, given that malaria selection must have been high enough in other areas to elevate the frequency of thalassaemia, we consider it unlikely that the βS hotspot pattern is due solely to insufficient malaria pressure in non-hotspot regions.

The models presented here highlight the importance of a low level of gene flow for the formation of βS hotspots. A pertinent historical question, then, is whether the various βS foci have had any particular reason for isolation. In the Greek cases, the geography of Greece itself (valley communities separated by mountain ranges) could be a factor. The Cholomondas mountains in the central part of Chalkidhiki may have limited gene flow from the coastal communities at the tip of the peninsular, where βS is at its highest.

Overall, we have demonstrated that epistasis enhances the ability of thalassaemias to keep βS from taking over a population. The formation of βS hotspots requires generally low levels of gene flow, but epistasis extends the range of gene flow values where hotspots are possible. Given the potential >2000 year window of βS challenge in the Mediterranean and the high degree of malaria protection βS offers, the circumscription of its present day distribution suggests something has acted to curtail its spread. We posit that epistasis is an important contributor to this unknown force.

Adams J.W.,Kasakoff A.B.. Factors Underlying Endogamous Group SizeSmith C.Regional AnalysisYear: 1976Academic PressLondon149
Allison A.C.. Protection afforded by sickle-cell trait against subtertian malareal infectionBr. Med. J.1Year: 195429029413115700
Barnicot N.A.,Allison A.C.,Blumberg B.S.,Deliyannis G.,Krimbas C.,Ballas A.. Haemoglobin types in Greek populationsAnn. Hum. Genet.26Year: 196322923613969412
Bastian, M., Heymann, S., Jacomy, M., 2009. Gephi: An Open Source Software for Exploring and Manipulating Networks.
Boussiou M.,Loukopoulos D.,Christakis J.,Fessas P.. The origin of the sickle mutation in Greece; evidence from β(S) globin gene cluster polymorphismsHemoglobin15Year: 19914594671687685
Buckee C.O.,Koelle K.,Mustard M.J.,Gupta S.. The effects of host contact network structure on pathogen diversity and strain structurePNAS10129Year: 2004108391084415247422
Buckee C.O.,Danon L.,Gupta S.. Host community structure and the maintenance of pathogen diversityProc. R. Soc. B: Biol. Sci.2741619Year: 200717151721
Cao A.,Congiu R.,Sollaino M.C.,Desogus M.F.,Demartis F.R.,Loi D.,Cau M.,Galanello R.. Thalassaemia and glucose-6-phosphate dehydrogenase screening in 13- to 14-year-old students of the Sardinian population: preliminary findingsCommunity Genet.11Year: 200812112818376107
Chotivanich K.,Udomsangpetch R.,Pattanapanyasat K.,Chierakul W.,Simpson J.,Looareesuwan S.,White N.. Hemoglobin E: a balanced polymorphism protective against high parasitemias and thus severe P. falciparum malariaBlood100Year: 20021172117612149194
Cushman S.A.,Shirk A.,Landguth E.L.. Separating the effects of habitat area, fragmentation and matrix resistance on genetic differentiation in complex landscapesLandscape Ecol.27Year: 2011369380
Fix A.. Migration and Colonization in Human MicroevolutionYear: 1999Cambridge University Press
Fix A.. Simulating hemoglobin historyHum. Biol.75Year: 200360761814655880
Flint J.,Harding R.M.,Boyce A.J.,Clegg J.B.. The population genetics of the haemoglobinopathiesBaillieres Clin. Heamatol.11Year: 1998151
Flint J.,Hill A.V.S.,Bowden D.K.. High frequencies of α-thalassaemia are the result of natural selection by malariaNature321Year: 19867447503713863
Galanello R.,Dessi E.,Melis M.A.,Addis M.,Sanna M.A.,Rosatelli C.,Argiolu F.,Giagu N.,Turco M.P.,Carace E.,Pirastu M.,Cao A.. Molecular analysis of β(o)-thalassemia intermedia in SardiniaBlood74Year: 19898238272473806
Haldane J.B.S.. Disease and EvolutionRic. Sci. Suppl: A19Year: 194968
Hanski I.,Gilpin M.. Metapopulation dynamics: brief history and conceptual domainBiol. J. Linn. Soc.42Year: 1991316
Hedrick P.W.. Genetic polymorphism in heterogeneous environments: the age of genomicsAnnu. Rev. Ecol. Evol. Syst.37Year: 20066793
Hedrick P.W.. Population genetics of malaria resistance in humansHeredity107Year: 201128330421427751
Hill A.,Bowden D.,O’Shaughnessy D.,Weatherall D.,Clegg J.. Beta thalassemia in Melanesia: association with malaria and characterization of a common variant (IVS-1 nt 5 G––C)Blood72Year: 19889142898955
Kan Y.W.,Nathan D.G.. Mild thalassemia: the result of interactions of alpha and beta thalassemia genesJ. Clin. Invest.49Year: 19706356425443168
Kalleas C.,Anagnostopoulos K.,Sinopoulou K.,Delaki E.,Margaritis D.,Bourikas G.,Tsatalas C.,Kortsaris A.,Tentes I.. Phenotype and genotype frequency of ß-thalassemia and sickle cell disease carriers in Halkidiki, Northern GreeceHemoglobin36Year: 2012647222188117
Kinder, H., Hilgemann, W., 2003 edition. The Penguin Atlas of World History, Vol. 1, Penguin Books Ltd, London.
Lappalainen T.,Montgomery S.B.,Nica A.C.,Dermitzakis E.T.. Epistatic selection between coding and regulatory variation in human evolution and diseaseAm. J. Hum. Genet.893Year: 201145946321907014
Lawton-Rauh A.. Demographic processes shaping genetic variationCurr. Opin. Plant Biol.11Year: 200810310918353707
Livingstone F.B.. Gene frequency clines of the beta hemoglobin locus in various human populations and their simulation by models involving differential selectionHum. Biol.41Year: 19692232365808784
Livingstone F.B.. The founder effect and deleterious genesAm. J. Phys. Anthropol.30Year: 196955605780379
Livingstone F.B.. Hemoglobin history in West AfricaHum. Biol.48Year: 1976487500789227
Livingstone F.B.. Simulation of the diffusion of the beta globin variants in the old worldHum. Biol.61Year: 19892973092680886
May J.,Evans J.A.,Timmann C.,Ehmen C.,Busch W.,Thye T.,Agbenyega T.,Horstmann R.D.. Hemoglobin variants and disease manifestations in severe falciparum malariaJAMA297Year: 20072220222617519411
Mandel S.P.H.. The stability of a multiple allelic systemHeredity13Year: 1958289302
Modiano D.,Luoni G.,Sirima B.S.,Simpore´ J.,Verra F.,Konate´ A.,Rastrelli E.,Olivieri A.,Calissano C.,Paganotti G.M.,D’Urbano L.,Sanou I.,Sawadogo A.,Modiano G.,Coluzzi M.. Haemoglobin C protects against clinical Plasmodium falciparum malariaNature414Year: 200130530811713529
Neto E.D.C.,de Oliviera D.M.,Rosas A.,Campos P.R.A.. The effect of spatially correlated environments on genetic diversity-area relationshipsJ. Theor. Biol.288Year: 2011576521872606
Pagnier J.,Mears J.G.,Dunda Belkhodja O.. Evidence for the multicentric origin of the sickle cell hemoglobin gene in AfricaProc. Natl. Acad. Sci. USA81Year: 1984177117736584911
Penman B.S.,Pybus O.G.,Weatherall D.J.,Gupta S.. Epistatic interactions between genetic disorders of hemoglobin can explain why the sickle-cell gene is uncommon in the MediterraneanPNAS106Year: 2009212422124619955437
Penman B.S.,Habib S.,Kanchan K.,Gupta S.. Negative epistasis between α+ thalassaemia and sickle cell trait can explain interpopulation variation in South AsiaEvolution65Year: 20113625363222133230
Salares R.,Bouwman A.,Anderung C.. The spread of malaria to Southern Europe in Antiquity: New approaches to old problemsMed. Hist.48Year: 200431132816021928
Schiliro G.,Di Gregorio F.,Romeo M.A.. Incidence of hemoglobin S carriers in SicilyHemoglobin10Year: 198695993957693
Siniscalco M.,Bernini L.,Fillipi G.,Latte B.,Meera Kahn P.,Piomelli S.,Rattazzi M.. Population genetics of haemoglobin variants, thalassaemia and glucose-6-phosphate dehydrogenase deficiency with particular reference to the malaria hypothesisBull. World Health Organ.34Year: 19663793935296398
Stamatoyannopoulos G.,Fessas P.H.. Thalassaemia, glucose-6-phosphate dehydrogenase deficiency, sickling and malarial endemicity in Greece. a study of five areasBr. Med. J.1Year: 196487587914102772
Taylor S.M.,Parobek C.M.,Fairhurst R.M.. Haemoglobinopathies and the clinical epidemiology of malaria: a systematic review and meta-analysisLancet Infect. Dis.12Year: 201245746822445352
Wainscoat J.S.,Kanavakis E.,Wood W.G.. Thalassaemia intermedia in cyprus: the interaction of α and β thalassaemiaBr. J. Haematol.53Year: 19834114166297530
Wainscoat J.S.,Old J.M.,Weatherall D.J.,Orkin S.H.. The molecular basis for the clinical diversity of β-thalassaemia in cypriotsLancet1Year: 1983123512376134037
Watts D.J.,Strogatz S.H.. Collective dynamics of ‘small-world networksNature393Year: 19984404429623998
Weatherall D.J.,Clegg J.B.. The Thalassaemia Syndromesfourth ed.Year: 2001Blackwell ScienceOxford, United Kingdom
Weatherall D.J.,Clegg J.B.. Inherited haemoglobin disorders: an increasing global health problemBull. World Health Organ.79Year: 200170471211545326
Williams T.N.,Mwangi T.W.,Wambua S.,Peto T.E.A.,Weatherall D.J.,Gupta S.,Recker M.,Penman B.S.,Uyoga S.,Macharia A.,Mwacharo J.K.,Snow R.W.,Marsh K.. Negative epistasis between the malaria-protective effects of alpha+ thalassemia and the sickle cell traitNat. Genet.37Year: 20051253125716227994
Williams T.N.,Mwangi T.W.,Wambua S.,Alexander N.D.,Kortok M.,Snow R.W.,Marsh K.. Sickle cell trait and the risk of Plasmodium falciparum malaria and other childhood diseasesJ. Infect. Dis.192Year: 200517818615942909
Appendix A  Supplementary data

Supplementary data 1
Click here for additional data file (mmc1.docx)


BSP is a Sir Henry Wellcome Postdoctoral fellow and a Junior Research Fellow at Merton College, Oxford. COB was funded for this study by a Sir Henry Wellcome Postdoctoral Fellowship. SG is a Royal Society Wolfson Research Fellow and an ERC Advanced Investigator. This work was supported by an ERC Advanced Grant (Diversity). We are also grateful to the Christopher Welch Trust for financial support. Finally, we would like to thank our two reviewers and editor for their very helpful comments on this manuscript.

Article Categories:
  • Article

Keywords: Keywords Sickle cell, Thalassaemia, Malaria, Mediterranean, Epistasis, Metapopulation.

Previous Document:  Intravenous high-dose enzyme replacement therapy with recombinant palmitoyl-protein thioesterase red...
Next Document:  La Salpêtrière Hospital before Charcot: a visit described by Pedro Gonzalez Velasco.