Inferred relatedness and heritability in malaria parasites

Abstract

Malaria parasites vary in phenotypic traits of biomedical or biological interest such as growth rate, virulence, sex ratio and drug resistance, and there is considerable interest in identifying the genes that underlie this variation. An important first step is to determine trait heritability (H²). We evaluate two approaches to measuring H² in natural parasite populations using relatedness inferred from genetic marker data. We collected single-clone Plasmodium falciparum infections from 185 patients from the Thailand–Burma border, monitored parasite clearance following treatment with artemisinin combination therapy (ACT), measured resistance to six antimalarial drugs and genotyped parasites using 335 microsatellites. We found strong relatedness structure. There were 27 groups of two to eight clonally identical (CI) parasites, and 74 per cent of parasites showed significant relatedness to one or more other parasites. Initially, we used matrices of allele sharing and variance components (VC) methods to estimate H². Inhibitory concentrations (IC₅₀) for six drugs showed significant H² (0.24 to 0.79, p = 0.06 to 2.85 × 10⁻⁹), demonstrating that this study design has adequate power. However, a phenotype of current interest—parasite clearance following ACT—showed no detectable heritability (H²= 0–0.09, ns) in this population. The existence of CI parasites allows the use of a simple ANOVA approach for quantifying H², analogous to that used in human twin studies. This gave similar results to the VC method and requires considerably less genotyping information. We conclude (i) that H² can be effectively measured in malaria parasite populations using minimal genotype data, allowing rational design of genome-wide association studies; and (ii) while drug response (IC₅₀) shows significant H², parasite clearance following ACT was not heritable in the population studied.