- Open Access
Association of sub-microscopic malaria parasite carriage with transmission intensity in north-eastern Tanzania
Malaria Journalvolume 10, Article number: 370 (2011)
In malaria endemic areas, individuals are frequently asymptomatic and may be undetected by conventional microscopy or newer, rapid diagnostic tests. Molecular techniques allow a more accurate assessment of this asymptomatic parasite burden, the extent of which is important for malaria control. This study examines the relative prevalence of sub-microscopic level parasite carriage and clonal complexity of infections (multiplicity of infection) over a range of endemicities in a region of north-eastern Tanzania where altitude is an established proxy of malaria transmission. The PCR prevalence was then compared against other measures of transmission intensity collected in the same area.
This study used 1,121 blood samples collected from a previously conducted cross-sectional malario-metric survey during the short rainy season in 2001 from 13 villages (three at < 600 m, four at 600-1,200 m and six at > 1,200 m in altitude above sea level). Samples were analysed by PCR for carriage of parasites and multiplicity of infection. These data were compared with other measures of transmission intensity collected from the same area.
Parasite prevalence was 34.7% by PCR and 13.6% by microscopy; a 2.5-fold difference in line with other recent observations. This fold difference was relatively consistent at the different altitude bands despite a marked decrease in parasite prevalence with altitude: < 600 m 70.9 vs 28.6, 600-1,200 m 35.5 vs 9.9, > 1,200 m 15.8 vs 5.9. The difference between parasite prevalence by PCR was 3.2 in individuals aged between 15 and 45 years (34.5 vs 10.9) compared with 2.5 in those aged 1-5 (34.0 vs 13.5) though this was not statistically significant. Multiplicity of infection (MOI) ranged from 1.2 to 3.7 and was positively associated with parasite prevalence assessed by both PCR and microscopy. There was no association of MOI and age.
Village level PCR parasite prevalence was strongly correlated with altitude, sero-conversion rate and predicted entomological inoculation rate.
Asymptomatic, low density, multi-clone malaria infection was common in this study area. These infections are important as potential contributors to the infectious reservoir of parasites and need to be identified by control programmes especially in this era where malaria elimination is a focus. High throughput standardized PCR approaches are needed to identify individuals who are malaria carriers.
Microscopy is routinely used for malaria diagnosis and epidemiological studies [1, 2]. However, it has limitations due to the subjective nature and sensitivity of slide reading and its time-consuming nature when carrying out studies involving a large number of individuals [1, 3–5]. There is an increase in the use of rapid diagnostic tests (RDTs) that are based on detection of parasite antigens, although they too have limitations related to sensitivity and discriminating current from recent infections . Molecular diagnostic tools, such as those based on PCR, have 10-100 greater sensitivity compared to microscopy and have been used increasingly for assessing infection [4–6].
In malaria endemic areas, asymptomatic malaria parasite carriers especially adults are not uncommon and, as potential gametocyte carriers, represent an important reservoir for malaria transmission . Many of these asymptomatic infections are present at densities below the limit for microscopic detection and, therefore, use of microscopy is likely to lead to underestimation of the malaria burden. Indeed, in a meta-analysis of community-based studies that use PCR detection of parasites, Okell et al showed that microscopy only detected 50% of the infections identified by PCR .
This current study was conducted in order to examine the relationship between PCR and microscopy in a defined geographical area where altitude has been shown to be a proxy for malaria transmission intensity . Risk factors for PCR carriage were examined at the individual and village level and the correlations between PCR prevalence and other measures of transmission intensity were examined.
Study site and sample collection
The study area, study design, and sampling have been described elsewhere with parasitological, haematological, serological and entomological measures all demonstrating malaria transmission intensity decreases with altitude [9–11]. Briefly, the study area is based in north-eastern Tanzania and runs from Kilimanjaro, through the Eastern Arc Pare and Usambara Mountains to Tanga on the coastal plain. Malaria prevalence has been shown to decrease with altitude and with distance from the coast, which is linked to average annual precipitation. In 2001, cross-sectional surveys were conducted during the short rainy season (October to November 2001) and finger-prick blood samples were collected into EDTA-coated tubes from an age-stratified sample of approximately 250 people from each of 13 villages (Figure 1) from three different age groups: 0-4 (n = 80), 5-14 (n = 80) and 15-45(n = 90) years old. The villages were situated at differing altitudes and in transects selected to be as similar as possible in terms of ethnic group, socio-economic status and accessibility to health care. Full details of demographic and clinical procedures can be found in Drakeley et al .
For PCR analysis, samples were chosen at random from the 13 villages to include a minimum of 60 samples from each village, with 20 from each of three age groups: 0-4, 5-14 and 15-45 years old. A total of 1,121 samples were analysed for the carriage of Plasmodium falciparum using nested PCR parasite. Prevalence by microscopy had been determined previously and results published in Drakeley et al .
DNA extraction and polymerase chain reaction
DNA was extracted from the archived, pelleted samples, using Nucleon kit protocol for DNA extraction according to manufacturer instructions http://www.tepnel.com. A nested PCR was used to amplify species-specific sequences of the small sub-unit ribosomal ribonucleic acid (18S SSU rRNA) genes of P. falciparum as described [12, 13]. DNA from the culture of P. falciparum (3D7 strain) as a positive control, DNA from blood samples of an individual never exposed to malaria as a negative control, and no template as a negative control, were included in each set of PCR for quality control. The PCR was carried out in a tetrad thermo-cycler (PTC-0240, The DNA engine Tetrad2® Thermal Cycler, Bio-Rad, Hercules, California, USA), followed by gel electrophoresis to assign individuals as either parasite positive or negative.
Nested PCR was carried out to determine the numbers of msp2 (FC27 and IC/3D7) clones. To amplify DNA, a primer pair corresponding to the outer conserved region of the polymorphic repetitive block 3 of msp2 was used. Then the PCR product was re-amplified in a pair of reactions using primers specific for FC27 and IC/3D7 allelic types of msp2 as described by Snounou et al . Number of products, corresponding to number of infecting FC27 and IC/3D7 clones, was counted after visualization on ethidium bromide (EtBr) stained 2.5 and 1.5% agarose gel, respectively.
The difference in PCR-based parasite prevalence in different villages was assessed by using chi-square and logistic regression to determine the association of parasite prevalence by PCR with different predictors as independent variables, such as age, distance from the coast (a crude proxy for rainfall) and altitude. PCR prevalence and MOI were examined for correlation with altitude, sero-conversion rate (SCR) for age related acquisition of antibodies to apical membrane antigen-1(AMA-1) and measured and extrapolated entomological inoculation rate (EIR)  that had been collected previously. Regression analysis was conducted in Stata using the 'svy' function to allow for stratified analysis with village as the primary sampling unit.
Ethical clearance was obtained from the London School of Hygiene & Tropical Medicine (LSHTM), Kilimanjaro Christian Medical Centre, Tanzania (KCMC) and the Tanzanian National Medical Research Institute (NIMR).
Microscopy and PCR results
The comparison between the 1,116 samples with both PCR and microscopy results is shown in Table 1. Of these, 62.8% were negative and 9.2% positive by both methods, however, PCR identified more parasite-positive individuals than microscopy: 34.4 vs 12.0% respectively. Thirty-one of 1,116 (2.7%) samples were positive by microscopy but negative by PCR; the majority of these samples (22) had parasite densities ≤ 80 parasites/μl. The samples that were microscopy positive but P. falciparum negative by PCR were analysed for other species. Two of the samples were Plasmodium malariae positive.
Using PCR as a gold-standard technique for malaria parasite detection, the sensitivity and specificity of microscopy were 26.8% (103/384) and 95.8% (701/732), respectively. The positive and negative predictive values were 78.4% (105/134) and 71.2% (699/982).
Comparison between microscopic and PCR parasite prevalence in the study area
The 2.5-fold difference in the prevalence of parasitaemia by PCR compared to microscopy is in line with other recent observations. The difference was observed in all villages at different altitudes (Table 2 and Figure 2) and was relatively consistent at the different altitude bands despite a marked decrease in parasite prevalence with altitude: < 600 m 70.9 vs 28.6, 600-1,200 m 35.5 vs 9.9, > 1,200 m 15.8 vs 5.9. The difference between parasite prevalence by PCR was 3.2 in individuals aged between 15 and 45 years (34.5 vs 10.9) compared with 2.5 in those aged 1-5 (34.0 vs 13.5), though this was not statistically significant. The highest prevalence was observed in the middle (5-14 years) age group and was lowest in adults (15-45 years) (Table 3).
Mixed parasite species infections
This was assessed in samples that were microscopy negative but P. falciparum PCR positive, and microscopy positive but PCR P. falciparum negative. A total of 117 samples were analysed for mixed infections by PCR, 12% of individuals had co-infections of P. falciparum and either P. malariae or P. ovale, and all were from lowlands. None of the samples analysed was Plasmodium vivax. The majority of samples (85.7%) with mixed infection were of P. falciparum and P. malariae. Six samples had Plasmodium ovale infection and out of these, four samples had mixed infection of P. falciparum, P. malariae and P. ovale.
Factors related to carriage of sub-microscopic parasites in the study area
After excluding individuals who were microscopically parasite positive, a total of 982 individuals were further involved in the analysis using logistic regression. In uni-variate logistic regression (Table 4), residents in the Tanga region were five-fold more likely to carry sub-microscopic parasites (OR, 5.00 [95%CI, 3.72-6.73]; p < 0.0001). This was also apparent within the Tanga region with both West Usambara and Tanga coast (Mgome village) transects, which were associated with increased carriage of sub-microscopic parasites; the magnitude being approximately four-fold (OR, 4.46 [95% CI, 0.88-22.49]; p = 0.067) and more than ten-fold (OR, 23.86 [95% CI, 7.00-81.29]; p < 0.0001) respectively. Analysis based on altitude band, regardless of region, indicated that individuals living in areas where the altitude was 600 m-1,200 m (OR, 0.28 [95%CI, 0.08-1.01]; p = 0.052) and above 1,200 m (OR, 0.10 [95%CI, 0.02-0.43]; p = 0.005) were less likely to carry sub-microscopic parasites as compared to individuals living in the area of < 600 m. Use of bednets, sex, fever and treatment were not associated with sub-microscopic parasite carriage in individuals (p > 0.05).
In multi-variate analysis, after adjustment for transect of residence, age, use of bed net, sex, fever and treatment, altitude band was still significantly related to carriage of the sub-microscopic P. falciparum, indicating that individuals living in highlands were less likely to carry sub-microscopic parasites.
Individuals with sub-microscopic parasites had significantly lower mean haemoglobin levels (11.4 vs 12.1 g/dl respectively, F = 19.2 p < 0.0001), suggestive of longer term or chronic infections. However, an interaction with of other factors (such as iron deficiency anaemia, helminth infections, HIV, red blood cell polymorphisms) that are associated with decreased levels of haemoglobin cannot be ruled out. No significant difference in prevalence of sub-microscopic parasites was seen between the sexes or in those who reported fever or recent anti-malarial treatment.
Multiplicity of infection (MOI)
The results showed that the majority of PCR positive subjects carried multiple strains; 76.2% (189/248). The MOI in PCR positive samples was significantly higher in high transmission/low altitude area (P = 0.025). There was no significant association between age and MOI (p > 0.05).
The distribution of msp2 allelic families with more than one clone, differed between the two allelic families; 46.8% (116/248) of parasites from 3D7/IC, 37.9% (94/248) of parasites from FC27 and 55.6% (138/248) of clones from both allelic families.
The use of different measures to correlate malaria and transmission intensity
Comparison of the different measures of transmission intensity (Table 5, Figure 2) indicates that all measures were strongly correlated. Relevant to this study, highly significant correlation co-efficients were observed between altitude and parasite rate by PCR and by microscopy, predicted EIR and sero-conversion rate.
Accurate quantification of the malaria parasite reservoir is important for control given the renewed focus on transmission reduction leading to elimination. Assessment of this reservoir by microscopy has limitations, and molecular techniques, largely based on the PCR, have been increasingly widely used to measure parasite prevalence. In this study, samples were analysed from previously conducted malario-metric surveys in an area where there is a very broad range of transmission intensity, i.e. microscopic parasite rates ranging from 2 to 50% and estimated entomological inoculation rates from < 1 to > 100 infectious bites per person per year. This variation occurs in a relatively small geographical area within similar ethnic groups, such that a number of confounders are minimized. Overall, a 2.5-fold higher parasite prevalence was observed with PCR compared to microscopy and this fold difference was consistent across the different transmission levels. The scale of the difference described here compares well with the two-fold difference documented in the meta-analysis by Okell and colleagues  and also from individual studies in Cambodia  and Solomon Islands .
The consistency of these observations on sub-microscopic parasitaemia suggests that this is a robust finding and one that has important implications for control measures. A study by Schneider et al  found that the contribution of sub-microscopic parasites to malaria transmission in individuals was similar to those individuals having microscopic gametocytes. Though gametocyte carriage was not assessed directly in this study, several studies document a high prevalence of concomitant gametocytes in individuals with asexual parasites (e.g. ). Although, gametocytes are likely to be present at lower densities  than asexual parasites, one interpretation of these observations is that using PCR for asexual parasites could be taken to equate to potential infectiousness.
Other studies [21, 22] have reported that sub-microscopic carriage differs in relation to transmission intensity and age. Strong correlations with malaria transmission intensity, as characterized by altitude and transect (a proxy for rainfall which increases with increasing proximity to the Indian Ocean)  were observed with carriage of sub-microscopic parasites. On a micro-epidemiological level, differences were observed between villages that were similar in altitude and distance from the coast (e.g. Kadando and Kileo) showing that local factors that favour transmission will influence both microscopic and sub-microscopic parasite carriage. Similar findings have been encountered in the highlands in Kenya [23, 24] indicating that small differences in altitude and other factors can influence the ecological environments of mosquitoes.
In this study, older age was not associated with increased carriage of sub-microscopic parasites. The reasons for this are not as obvious as one might expect that individuals with higher levels of immunity (i.e. older ages) would maintain parasites at lower densities [5, 16]. At the lowest transmission intensity (highest altitude band) individuals older than 15 were 3.0 times more likely to have sub-microscopic parasites than younger individuals (chi sq 4.01, p = 0.0425) but this was not significant in multi-variant analysis. Other individual level factors such as sex, history of fever and recent treatment were not associated with parasite carriage. Use of bed nets was found to be inconsistently protective against sub-microscopic carriage; at higher transmission intensities where 50% net users were PCR positive compared to 75% of non net users at the lowest altitude, with the corresponding figures of 14 and 40% in the middle altitude band. However, the effect was not significant in stratified, multi-variate analysis.
The MOI results reported in this study provide a rough estimate of the genetic diversity of MSP2, a polymorphic marker, and suggest diversity was pronounced. Multiple parasites genotypes in individuals living in endemic areas are frequently observed, as was the case in other studies in Senegal [21, 22] and Ghana . A study in Ghana indicated that MOI increases the gaining of immunity and the ability to control parasitaemia and protection from subsequent clinical episodes . MOI did not correlate to age; similar findings observed in a more restricted age range of 2-10 years in Senegal .
Sub-microscopic parasite carriage was strongly correlated with MOI and other previously collected measures of transmission intensity such as parasite prevalence by microscopy, extrapolated EIR and SCR. The relationship presented in Table 5 is for the antigen AMA-1 and, the correlation between PCR prevalence and MOI for Merozoite Surface Protein-1 is similarly strong (0.8525, p = 0.0004 & 0.7716, p = 0.017 respectively). This is not surprising as microscopy and PCR parasite prevalence were highly correlated (coefficient r2 = 0.9, p < 0.001) and the relationship between microscopy and the other measures has been demonstrated previously . However, as PCR methodology becomes increasingly more widespread, it will be important to calibrate PCR prevalence across the full range of transmission intensity. The observations made here and in other papers suggest that PCR prevalence will have a greater dynamic range than microscopy and thus will be an important parameter for assessing progress in low endemicity and pre-elimination settings.
This study shows that a large number of individuals carried sub-microscopic parasites in a malaria endemic area of Tanzania. These individuals are a potential reservoir of infection in the population. There is, therefore, a need to further evaluate the burden of sub-microscopic parasites when considering elimination and eradication of malaria. Employing the multiplex PCR for the different species of Plasmodium and also a more sensitive method, such as stevor PCR suggested elsewhere , may be needed to resolve the discordant results of the conventional PCR and microscopy.
Menge DM, Ernst KC, Vulule JM, Zimmerman PA, Guo H, John CC: Microscopy underestimates the frequency of Plasmodium falciparum infection in symptomatic individuals in a low transmission highland area. Am J Trop Med Hyg. 2008, 79: 173-177.
Msellem MI, Màrtensson A, Rotllant G, Bhattarai A, Strömberg J, Kahigwa E, Garcia M, Petzold M, Olumese P, Ali A, Björkman A: Influence of rapid malaria diagnostic tests on treatment and health outcome in fever patients, Zanzibar-a crossover validation study. PLoS Med. 2009, 6: e1000070-10.1371/journal.pmed.1000070.
Shekalaghe SA, Bousema JT, Kunei KK, Lushino P, Masokoto A, Wolters LR, Mwakalinga S, Mosha FW, Sauerwein RW, Drakeley CJ: Submicroscopic Plasmodium falciparum gametocyte carriage is common in an area of low and seasonal transmission in Tanzania. Trop Med Int Health. 2007, 12: 547-553. 10.1111/j.1365-3156.2007.01821.x.
Greenwood B: The molecular epidemiology of malaria. Trop Med Int Health. 2002, 7: 1012-1021. 10.1046/j.1365-3156.2002.00980.x.
Dal-Bianco MP, Köster KB, Kombila UD, Kun JF, Grobusch MP, Ngoma GM, Matsiegui PB, Supan C, Salazar CL, Missinou MA, Issifou S, Lell B, Kremsner P: High prevalence of asymptomatic Plasmodium falciparum infection in Gabonese adults. Am J Trop Med Hyg. 2007, 77: 939-942.
Babiker HA, Schneider P: Application of molecular methods for monitoring transmission stages of malaria parasites. Biomed Mater. 2008, 3: 034007-10.1088/1748-6041/3/3/034007.
Alves FP, Gil LH, Marrelli MT, Ribolla PE, Camargo EP, Da Silva LH: Asymptomatic carriers of Plasmodium spp. as infection source for malaria vector mosquitoes in the Brazilian Amazon. J Med Entomol. 2005, 42: 777-779. 10.1603/0022-2585(2005)042[0777:ACOPSA]2.0.CO;2.
Okell LC, Ghani AC, Lyons E, Drakeley CJ: Submicroscopic infection in Plasmodium falciparum-endemic populations: a systematic review and meta-analysis. J Infect Dis. 2009, 200: 1509-1517. 10.1086/644781.
Drakeley CJ, Carneiro I, Reyburn H, Malima R, Lusingu JP, Cox J, Theander TG, Nkya WM, Lemnge MM, Riley EM: Altitude-dependent and -independent variations in Plasmodium falciparum prevalence in northeastern Tanzania. J Infect Dis. 2005, 191: 1589-1598. 10.1086/429669.
Drakeley CJ, Corran PH, Coleman PG, Tongren JE, McDonald SL, Carneiro I, Malima R, Lusingu JP, Manjurano A, Nkya WM, Lemnge MM, Cox J, Reyburn H, Riley EM: Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proc Natl Acad Sci USA. 2005, 102: 5108-5113. 10.1073/pnas.0408725102.
Carneiro IA, Drakeley CJ, Owusu-Agyei S, Mmbando B, Chandramohan D: Haemoglobin and haematocrit: is the threefold conversion valid for assessing anaemia in malaria endemic settings?. Malar J. 2007, 6: 67-10.1186/1475-2875-6-67.
Snounou G, Viriyakosol S, Zhu XP, Jarra W, Pinheiro L, do Rosario VE, Thaithong S, Brown KN: High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993, 61: 315-320. 10.1016/0166-6851(93)90077-B.
Singh B, Bobogare A, Cox-Singh J, Snounou G, Abdullah SM, Rahman AH: A genus- and species -specific nested polymerase chain reaction malaria detection assay for epidemiologic studies. Am J Trop Med Hyg. 1999, 60: 687-692.
Snounou G, Zhu X, Siripon N, Jarra W, Thaithong S, Brown KN, Viriyakosol S: Biased distribution of msp1 and msp2 allelic variants in Plasmodium falciparum populations in Thailand. Trans R Soc Trop Med Hyg. 1999, 93: 369-374. 10.1016/S0035-9203(99)90120-7.
Corran P, Coleman P, Riley E, Drakeley C: Serology: a robust indicator of malaria transmission intensity?. Trends Parasitol. 2007, 23: 575-582. 10.1016/j.pt.2007.08.023.
Steenkeste N, Rogers WO, Okell L, Jeanne I, Incardona S, Duval L, Chy S, Hewitt S, Chou M, Socheat D, Babin FX, Ariey F, Rogier C: Sub-microscopic malaria cases and mixed malaria infection in a remote area of high malaria endemicity in Rattanakiri province, Cambodia: implication for malaria elimination. Malar J. 2010, 9: 108-10.1186/1475-2875-9-108.
Harris I, Sharrock WW, Bain LM, Gray KA, Bobogare A, Boaz L, Lilley K, Krause D, Vallely A, Johnson ML, Gatton ML, Shanks GD, Cheng Q: A large proportion of asymptomatic Plasmodium infections with low and sub-microscopic parasite densities in the low transmission setting of Temotu Province, Solomon Islands: challenges for malaria diagnostics in an elimination setting. Malar J. 2010, 9: 254-10.1186/1475-2875-9-254.
Schneider P, Bousema JT, Gouagna LC, Otieno S, Van de Vegte-Bolmer M, Omar SA, Sauerwein RW: Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007, 76: 470-474.
Clyde DF: Mass administration of an antimalarial drug combining 4-aminoquinoline and 8-aminoquinoline in Tanganyika. Bull World Health Organ. 1962, 27: 203-212.
Bousema T, Drakeley C: Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011, 24: 377-410. 10.1128/CMR.00051-10.
Ntoumi F, Contamin H, Rogier C, Bonnefoy S, Trape JF, Mercereau-Puijalon O: Age-dependent carriage of multiple Plasmodium falciparum merozoite surface antigen-2 alleles in asymptomatic malaria infections. Am J Trop Med Hyg. 1995, 52: 81-88.
Vafa M, Troye-Blomberg M, Anchang J, Garcia A, Migot-Nabias F: Multiplicity of Plasmodium falciparum infection in asymptomatic children in Senegal: relation to transmission, age and erythrocyte variants. Malar J. 2008, 7: 17-10.1186/1475-2875-7-17.
Ernst KC, Adoka SO, Kowuor DO, Wilson ML, John CC: Malaria hotspot areas in a highland Kenya site are consistent in epidemic and non-epidemic years and are associated with ecological factors. Malar J. 2006, 5: 78-10.1186/1475-2875-5-78.
Githeko AK, Ayisi JM, Odada PK, Atieli FK, Ndenga BA, Githure JI, Yan G: Topograph and malaria transmission heterogeneity in western Kenya highlands: prospects for focal vector control. Malar J. 2006, 5: 107-10.1186/1475-2875-5-107.
Buchholz U, Kobbe R, Danquah I, Zanger P, Reither K, Abruquah HH, Grobusch MP, Ziniel P, May J, Mockenhaupt FP: Multiplicity of Plasmodium falciparum infection following intermittent preventive treatment in infants. Malar J. 2010, 9: 244-10.1186/1475-2875-9-244.
This study was conducted as part of the Joint Malaria Programme (JMP), a collaboration between the National Institute for Medical Research (NIMR), Kilimanjaro Christian Medical College (KCMC), the London School of Hygiene & Tropical Medicine (LSHTM), and the Centre for Medical Parasitology, University of Copenhagen (CMP). We thank members of the JMP community studies team and laboratory staff at KCMC and NIMR Amani. We thank Federrica Verra for her comments on this manuscript. Finally, we thank Colin Sutherland (LSHTM) for providing the nested PCR protocol and malaria parasite positive and negative controls for the sub-patent parasites carriage genotyping.
This study was funded by the UK Medical Research Council and the LSHTM Gates Malaria Partnership (GMP; Bill and Melinda Gates Foundation). AM was a MalariaGEN-funded doctoral student.
The authors declare that they have no competing interests.
AM, ER, and CD designed the study. HR, CD and RO were responsible for study participants recruitment, and clinical and parasitological examinations. AM, TL and SJ did the genotyping, and AM, LO, TC and CD the statistical analyses. AM, CD, TC and CR wrote the paper with major contributions of the other authors. All authors read and approved the final manuscript.