An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples [version 1; peer review: 2 approved]

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed. Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination. This manuscript from the MalariaGEN consortium, a data-sharing community of teams working on Plasmodium falciparum genomic epidemiology, presents the new release of curated P. falciparum genomes from isolates collected in 73 locations in Africa, Asia, South America and Oceania. Based on robust and perfectly detailed methods (ranging from the treatment of the blood samples, the DNA extraction, the Illumina and computational platforms developed to produce genome sequencing for variant discovery and genotype calling), they analyzed 7000 P. falciparum genome sequences and provided numerous exciting data. For instance, they found that variations (SNPs and indels) in P. falciparum genome affected about a quarter of the 23 Mb genome (and mostly coding regions), or that duplication genotypes are frequent around mdr1, plasmepsin2/3 and gch1, which are known to be associated with antimalarial drug resistance (including mefloquine, piperaquine and sulfadoxine/pyrimethamine). Moreover, population genetic analyses conducted on this largest available data resource, depict a comprehensive picture of P. falciparum parasite populations globally and sub populations at continental level. In the results, a large section is devoted to the description of the geographic patterns of validated molecular markers (SNPs and CNVs) associated with antimalarial drug resistance. By compiling data on all samples collected from 2002–2015, they present clear profiles of drug resistance by regional sub-populations for the most used antimalarial drugs. Finally, they reveal a global landscape regarding a major challenge for malaria elimination, that are deletions in hrp2 and 3 genes linked The analysis of whole-genome sequences obtained from Plasmodium falciparum is particularly challenging due to the presence of hypervariable regions, highly repetitive sequences, and frequent mixture of parasites due to multiple infections of the host. The authors of this study describe a curated list of over three million high-confidence polymorphisms obtained from the genome sequence analysis of more than 7000 samples of P. falciparum collected by several studies in 73 locations in Africa, Asia, South America and Oceania. This work, reporting a laudable effort to substantially enrich publicly available genome data of P. falciparum worldwide, is of paramount importance for the field. The contribution goes in line with authors' previous consortia publications, extending largely the number of available data that can be analysed via web with powerful data analysis pipelines. By providing open access to a curated list of polymorphisms based on reproducible and high-quality protocols for the sequencing and analysis of P. falciparum genomes this study is likely to decrease the difficulties that have delayed the research on genomic epidemiology and population genomics of P. falciparum . Among other advances, studies in this area are likely to have important implications for a better understanding of the evolution towards drug resistance of the different global parasite populations ultimately contributing for a better control of this devastating disease. The manuscript is very well written and clear. It presents eight genetically distinct populations of parasites each endemic to different word regions, including South America, West Africa, Central Africa, East Africa, South Asia, West Southeast Asia, East Southeast Asia and Oceania. An interesting genetic and geographic characterization of the eight parasite populations is also shown. Of note, the finding of higher within-host diversity in the parasite populations endemic to Africa, the identification of single nucleotide polymorphism with high levels of geographic differentiation, and further characterization of geographic patterns of drug resistance and polymorphisms with potential


Introduction
A major obstacle to malaria elimination is the great capacity of the parasite and vector populations to evolve in response to malaria control interventions. The widespread use of chloroquine and DDT in the 1950's led to high levels of drug and insecticide resistance, and the same pattern has been repeated for other first-line antimalarial drugs and insecticides. Over the past 15 years, mass distribution of pyrethroid-treated bednets in Africa and worldwide use of artemisinin combination therapy (ACT) has led to substantial reductions in malaria prevalence and mortality, but there are rapidly increasing levels of resistance to ACT in Southeast Asian parasites and of pyrethroid resistance in African mosquitoes. A deep understanding of local patterns of resistance and the continually changing nature of the local parasite and vector populations is necessary to manage the use of drugs and insecticides and to deploy public health resources for maximum sustainability and impact.
Current methods for genetic surveillance of the parasite population are largely based on targeted genotyping of specific loci, e.g. known markers of drug resistance. Whole genome sequencing of malaria parasites is currently more expensive and complex, particularly at the stage of data analysis, but it is an important adjunct to targeted genotyping, as it provides a more comprehensive picture of parasite genetic variation. It is particularly important for discovery of new drug resistance markers and for monitoring patterns of gene flow and evolutionary adaptation in the parasite population.
The Plasmodium falciparum Community Project (Pf Community Project) was established with the aim of integrating parasite genome sequencing into clinical and epidemiological studies of malaria (www.malariagen.net/projects). It forms part of the Malaria Genomic Epidemiology Network (MalariaGEN), a global data-sharing network comprising multiple partner studies, each with its own research objectives and led by a local investigator 1 . Genome sequencing was performed centrally, and partner studies were free to analyse and publish the genetic data produced on their own samples, in line with MalariaGEN's guiding principles on equitable data sharing [1][2][3] . A programme of capacity building for research into parasite genetics was developed at multiple sites in Africa alongside the Pf Community Project 4 .
The first phase of the project focused on developing simple methods to obtain purified parasite genome DNA from small blood samples collected in the field 5,6 and on establishing reliable computational methods for variant discovery and genotype calling from short-read sequencing data 7 . This presented a number of analytical challenges due to long tracts of highly repetitive sequence and hypervariable regions within the P. falciparum genome, and also because a single infection can contain a complex mixture of genotypes. Once a reliable analysis pipeline was in place, a process was established for periodic data releases to partners, with continual improvements in data quality as new analytical methods were developed.
Data from the Pf Community Project were initially released through a companion project called Pf3k, whose goal was to bring together leading analysts from multiple institutions to benchmark and standardise methods of variant discovery and genotyping calling. A visual analytics web application was developed 8 for researchers to explore the data. The open dataset was enlarged in 2016 when multiple partner studies contributed to a consortial publication on 3,488 samples from 23 countries 9 .
Data produced by the Pf Community Project have been used to address a broad range of research questions, both by the groups that generated samples and data and by the wider research community, and have generated over 50 previous publications (refs 5-55). These data have become a key resource for the epidemiology and population genetics of antimalarial drug resistance 9-22 and an important platform for the discovery of new genetic markers and mechanisms of resistance through genome-wide association studies 23-27 and combined genometranscriptome analysis 28 . The data have also been used to study gene deletions that cause failure of rapid diagnostic tests 29 ; to characterise genetic variation in malaria vaccine antigens 30,31 ; to screen for new vaccine candidates 32 ; to investigate specific host-parasite interactions 33,34 ; and to describe the evolutionary adaptation and diversification of local parasite populations 7,9,12,35-40 .
The Pf Community Project data also provide an important resource for developing and testing new analytical and computational methods. A key area of methods development is quantification of within-host diversity 7,41-46 , estimation of inbreeding 7,47 , and deconvolution of mixed infections into individual strains 48,49 . The data have also been used to develop and test methods for estimating identity by descent 50,51 , imputation 52 , typing structural variants 53 , designing other SNP genotyping platforms 54 and data visualisation 8,55 . In a companion study we performed whole genome sequencing of experimental genetic crosses of P. falciparum, and this provided a benchmark to test the accuracy of our genotyping methods, and to conduct an in-depth analysis of indels, structural variants and recombination events which are complicated to ascertain in these population genetic samples 56 .

Variant discovery and genotyping
We used the Illumina platform to produce genome sequencing data on all samples and we mapped the sequence reads against the P. falciparum 3D7 v3 reference genome. The median depth of coverage was 73 sequence reads averaged across the whole genome and across all samples. We constructed an analysis pipeline for variant discovery and genotyping, including stringent quality control filters that took into account the unusual features of the P. falciparum genome, incorporating lessons learnt from our previous work 7,56 and the Pf3k project, as outlined in the Methods section.  Countries are grouped into eight geographic regions based  on their geographic and genetic characteristics. For each country, the table reports: the number of distinct  sampling locations; the total number of samples sequenced; the number of high-quality samples included in  the analysis; and the percentage of samples collected between 2012-2015, the most recent sampling period  in the dataset. Eight samples were obtained from travellers returning from an endemic country, but where the  precise site of the infection could not be determined. These were reported from Ghana (3 sequenced samples/2 analysis set samples), Kenya (2/1), Uganda (2/1) and Mozambique (1/1). "Lab samples" contains all sequences obtained from long-term in vitro cultured and adapted isolates, e.g. laboratory strains. The breakdown by site is reported in Supplementary In the first stage of analysis we discovered variation at over six million positions, corresponding to about a quarter of the 23 Mb P. falciparum genome (Supplementary Data; Supplementary Table 3). These included 3,168,721 single nucleotide polymorphisms (SNPs): these were slightly more common in coding than non-coding regions and were mostly biallelic. The remaining 2,882,975 variants were predominantly short indels but also included more complex combinations of SNPs and indels: these were much more abundant in non-coding than coding regions, and mostly had at least three alleles. The predominance of indels in non-coding regions has been previously observed and is most likely a consequence of the extreme AT bias which leads to many short repetitive sequences 56,57 .
For the purpose of this analysis, we excluded all variants in subtelomeric and internal hypervariable regions, mitochondrial and apicoplast genomes, and some other regions of the genome where the mapping of short sequence reads is prone to a high error rate due to extremely high rates of variation 56 . A total of 1,838,733 SNPs (of which 1,626,886 were biallelic) and 1,276,027 indels (or SNP/indel combinations) passed all these filters. The pass rate for SNPs in coding regions (66%) was considerably higher than that for SNPs in non-coding regions (47%), indels in coding regions (37%) and indels in noncoding regions (47%). Finally, we removed samples with a low genotyping success rate or other quality control issues. We also removed replicates and 41 samples with genetic markers of infection by multiple Plasmodium species, leaving 5,970 high-quality samples from 28 countries (Table 1).
We used coverage and read pair analysis to determine duplication genotypes around mdr1, plasmepsin2/3 and gch1, each of which are associated with drug resistance. For each of these three genes we discovered many different sets of breakpoints (29, 10 and 3 pairs of breakpoints for mdr1, gch1, and plasmepsin 2/3, respectively), including complex rearrangements 58 that to the best of our knowledge have not been observed before in Plasmodium species (Supplementary Data; Supplementary Note, Supplementary Tables 4-6). We also used sequence reads coverage to identify large structural variants that appear to delete or disrupt hrp2 and hrp3, an event that can cause rapid diagnostic tests to malfunction.
The population genetic analyses in this paper are based on the filtered dataset of high-quality SNP genotypes in 5,970 samples. These data are openly available, together with annotated genotyping data on 6 million putative variants in all 7,113 samples, plus details of partner studies and sampling locations, at www.malariagen.net/resource/26.

Global population structure
The genetic structure of the global parasite population reflects its geographic regional structure 7,9,10 as illustrated by a neighbourjoining tree and a principal component analysis of all samples based on their SNP genotypes ( Figure 1). Based on these observations we grouped the samples into eight geographic regions: West Africa, Central Africa, East Africa, South Asia, the western part of Southeast Asia, the eastern part of Southeast Asia, Oceania and South America. Each of these can be viewed as a regional sub-population of parasites, which is more or less differentiated from other regional sub-populations depending on rates of gene flow and other factors. The different regions encompass a range of epidemiological and environmental settings, varying in transmission intensity, vector species and history of antimalarial drug usage. Note these regional classifications are intentionally broad, and therefore overlook many interesting aspects of local population structure, e.g. a distinctive Ethiopian sub-population can be identified by more detailed analysis of African samples 12 .
Genetically mixed infections were considerably more common in Africa than other regions, consistent with the high intensity of malaria transmission in Africa (Figure 2a) The average nucleotide diversity across the global sample collection was 0.040% (median=0.028%), i.e. two randomlyselected samples differ by an average of 4 nucleotide positions per 10kb. Levels of nucleotide diversity vary greatly across the genome 56 and also geographically ( Figure 2b). Distributions of values were highest in Africa, followed by Bangladesh, but the scale of regional differences was relatively modest, ranging from an average of 0.030% in Eastern Southeast Asia to 0.040% in West Africa (median=0.019% and 0.028% respectively; Figure 2b). In other words, the nucleotide diversity of each regional parasite population was not much less than that of the global parasite population. This is consistent with the idea that the global P. falciparum population has a common African origin and that historically there must have been significant levels of migration.
All regional sub-populations showed very low levels of linkage disequilibrium relative to human populations, e.g. r 2 decayed to <0.1 within 500 bp ( Figure 2c). As expected, African populations had the highest rates of LD decay, implying the highest levels of haplotype diversity.
Geographic patterns of population differentiation and gene flow Parasite sub-populations in different locations naturally tend to differentiate over time unless there is sufficient gene flow to counterbalance genetic drift. Genome-wide estimates of F ST provide an indicator of this process of genetic differentiation, which is partly determined by geographic distance (Figure 3). For example, we observe much greater genetic differentiation  between South America and South Asia (genome-wide average F ST 0.22) or between Africa and Oceania (0.20) than between sub-regions within Asia (<0.1) or within Africa (<0.02).
These data reveal some interesting exceptions to the general rule that genome-wide F ST is correlated with geographic distance. For example, African parasites are more strongly differentiated from Southeast Asian parasites (genome-wide average F ST 0.20) than they are from parasites in neighbouring Bangladesh (0.11). If this is examined in more detail, there is an unexpectedly steep gradient of genetic differentiation at the geographical boundary between South Asia and Southeast Asia, i.e. parasites sampled in Myanmar and Western Thailand are much more strongly differentiated from parasites sampled in Bangladesh (genome-wide F ST 0.07) than would be expected given that these are neighbouring countries. As discussed later, Southeast Asia is the global epicentre of antimalarial drug resistance, and these observations add to a growing body of evidence that Southeast Asian parasites have acquired a wide range of genomic features that are likely due to natural selection rather than genetic drift 23,40 . It is noteworthy that the level of genetic differentiation between western and eastern parts of Southeast Asia (genomewide F ST 0.05) is greater than between West Africa and East Africa (0.02) although the geographic distances are much greater in Africa. This is likely due to the lower intensity of malaria transmission in Southeast Asia, and in particular the presence of a malaria-free corridor running through Thailand, which act as barriers to gene flow across the region 23,40 .

Genes with high levels of geographic differentiation
The F ST metric can also be calculated for individual variants to identify specific genes that have acquired high levels of geographic differentiation relative to the genome as a whole. This can be done either at the global level (to identify variants that are highly differentiated between different regions of the world) or at the local level (to identify variants that are highly differentiated between different sampling locations within a region).
To identify variants that are strongly differentiated at the global level, we began by estimating F ST for each SNP across all of the eight regional sub-populations. The group of SNPs Each point represents one pairwise comparison between two regional parasite populations. The x-axis reports the geographic separation between the two populations, measured as great-circle distance between the centre of mass of each population and without taking into account natural barriers. The y-axis reports the genetic differentiation between the two populations, measured as average genome-wide F ST . Points are coloured based on the regional populations they represent: between African populations (red); between Asian populations (blue); between Southeast Asia (as a whole) and Oceania, Africa or South America (purple); all the rest (orange).
with the highest global F ST levels were found to be strongly enriched for non-synonymous mutations, suggesting that the process of differentiation is at least in part due to natural selection ( Figure 4). After ranking all SNPs according to their global F ST value, we calculated a global differentiation score for each gene based on the highest-ranking non-synonymous SNP within the gene (see Methods). All genes are ranked according to their global differentiation score in the accompanying data release, and those with the highest score are listed in Supplementary Table 7 (Supplementary Data). The most highly differentiated gene, p47, is known to interact with the mosquito immune system 59 and has two variants (S242L and V247A) that are at fixation in South America but absent in other geographic regions. Also among the five most highly differentiated genes are gig (implicated in gametocytogenesis 60 ), pfs16, (expressed on the surface of gametes 61 ) and ctrp (expressed on the ookinete cell surface and essential for mosquito infection 62 ). Thus, four of the five most highly differentiated parasite genes are involved in the process of transmission by the mosquito vector, raising the possibility that this reflects evolutionary adaptation of the P. falciparum population to the different Anopheles species that transmit malaria in different geographical regions.
It is more difficult to characterise variants that are strongly differentiated at the local level, due to smaller sample sizes and various sources of sampling bias, but a crude estimate can be obtained by analysis of each of the six geographical regions with samples from multiple countries. F ST was estimated for each SNP across different sampling locations within each geographical region, and the results for different regions were combined by a heuristic approach to obtain a local differentiation score for each gene (see Methods

Geographic patterns of drug resistance Classification of samples based on markers of drug resistance.
Antimalarial drug resistance represents a major focus of research for many partner studies within the Pf Community Project, and this dataset therefore contains a significant body of data that have appeared in previous reports on drug resistance. Readers are referred to these publications for more detailed analyses of local patterns of resistance 9-14,16-22 and of resistance to specific drugs including chloroquine 16,21 , sulfadoxine-pyrimethamine 16,19,21 and artemisinin combination therapy 9-11,13-15,17,18,21,22 .
Here we have classified all samples into different types of drug resistance based on published genetic markers and current knowledge of the molecular mechanisms (see www.malariagen.net/resource/26 for details of the heuristic used). Below we summarise the overall profile of drug resistance types in the regional sub-populations: this is intended simply to provide context for users of this dataset, and should not be regarded as a statement of the current epidemiological situation. The Supplementary Notes (Supplementary Data) contain a more detailed description of the geographical distribution of haplotypes, CNV breakpoints, interactions between genes, and variants associated with less commonly used antimalarial drugs.
In the accompanying data release, we also identify samples with mdr1, plasmepsin2/3 and gch1 gene amplifications that can affect drug resistance.
Chloroquine resistance. Samples were classified as chloroquine resistant if they carried the crt 76T allele. As shown in Table 2, this was found in almost all samples from Southeast Asia, South America and Oceania. It was also found across Africa but at lower frequencies, particularly in East Africa where chloroquine resistance is known to have declined since chloroquine was discontinued [63][64][65] . Supplementary Table 11 (Supplementary Data) shows the geographical distribution of different crt haplotypes (based on amino acid positions [72][73][74][75][76] which is consistent with the theory that chloroquine resistance spread from Southeast Asia to Africa with multiple independent origins in South America and Oceania 66,67 . The crt locus is also relevant to other types of drug resistance, e.g. crt variants that are relatively specific to Southeast Asia form the genetic background of artemisinin resistance, and newly emerging crt alleles have been associated with the spread of ACT failure due to piperaquine resistance 13,14,22,68 . Sulfadoxine-pyrimethamine resistance. Clinical resistance to sulfadoxine-pyrimethamine (SP) is determined by multiple

Resistance to artemisinin combination therapy.
We classified samples as artemisinin resistant based on the World Health Organization classification of non-synonymous mutations in the propeller region of the kelch13 gene that have been associated with delayed parasite clearance 71 . By this definition, artemisinin resistance was confined to Southeast Asia but, as previously reported, this dataset contains a substantial number of non-synonymous kelch13 propeller SNPs occurring at <5% frequency in Africa and elsewhere 9 . The most common ACT formulations in Southeast Asia are artesunate-mefloquine (AS-MQ) and dihydroartemisinin-piperaquine (DHA-PPQ). We classified samples as mefloquine resistant if they had mdr1 amplification 72 or as piperaquine resistant if they had plasmepsin 2/3 amplification 25 . Mefloquine resistance was observed throughout Southeast Asia and was most common in the western part. Piperaquine resistance was confined to eastern Southeast Asia with a notable concentration in western Cambodia. Elsewhere 11,13 we describe the kel1/pla1 lineage of artemisinin-and piperaquine-resistant parasites that expanded in western Cambodia during 2008-13, and then spread to other countries during 2013-18, causing high rates of DHA-PPQ treatment failure across eastern Southeast Asia: since the current dataset extends only to 2015 it captures only the first phase of the kel1/pla1 lineage expansion.

HRP2/3 deletions that affect rapid diagnostic tests
Rapid diagnostic tests (RDTs) provide a simple and inexpensive way to test for parasites in the blood of patients who are suspected to have malaria, and have become a vital tool for malaria control 73,74 . The most widely used RDTs are designed to detect P. falciparum histidine-rich protein 2 and cross-react with histidine-rich protein 3, encoded by the hrp2 and hrp3 genes respectively. Parasites with gene deletions of hrp2 and/or hrp3 have emerged as an important cause of RDT failure in a number of locations [75][76][77][78][79] . It is difficult to devise a simple genetic assay to monitor for risk of RDT failure because hrp2 and hrp3 deletions comprise a diverse mixture of large structural variations with multiple independent origins, and both genes are located in subtelomeric regions of the genome with very high levels of natural variation 29,80-83 . In the absence of a wellvalidated algorithmic method, we visually inspected sequence read coverage and identified samples with clear evidence of large structural variants that disrupted or deleted the hrp2 and hrp3 genes. We took a conservative approach: samples that appeared to have a mixture of deleted and non-deleted genotypes were classified as non-deleted.
Deletions were found at relatively high frequency in Peru (8 of 21 samples had hrp2 deletions, 14 had hrp3 deletions and 6 had both) but were not seen in samples from Colombia and were relatively rare outside South America. Oceania was the only other region where we observed hrp2 deletions, but at very low frequency (4%, n=3/80), and also had hrp3 deletions (25%) though no combined deletions were seen. Deletions of hrp3 only were more geographically widespread than hrp2 deletions, being common in Ethiopia (43%, n=9/21) and in Senegal (7%, n=6/84), and at relatively low frequency (<5%) in Kenya, Cambodia, Laos, and Vietnam (Supplementary Data; Supplementary Table 13). Note that these findings might underestimate the true prevalence of hrp2/hrp3 deletions, due to sampling bias (our samples were primarily collected from RDTpositive cases) and also because we focused on large structural variants and did not consider polymorphisms that might also cause RDT failure but would require more sophisticated analytical approaches. There is a need for more reliable diagnostics of hrp2 and hrp3 deletions, and we hope that these open data will accelerate this important area of applied methodological research.

Discussion
This open dataset comprises sequence reads and genotype calls on over 7,000 P. falciparum samples from MalariaGEN partner studies in 28 countries. After excluding variants and samples that failed to meet stringent quality control criteria, the dataset contains high-quality genotype calls for 3 million polymorphisms including SNPs, indels, CNVs and large structural variations, in almost 6,000 samples. The data can be analysed in their entirety or can be filtered to select for specific genes, or geographical locations, or samples with particular genotypes. This is twice the sample size of our previous consortial publication 9 and is the largest available data resource for analysis of P. falciparum population structure, gene flow and evolutionary adaptation. Each sample has been annotated to show its profile of resistance to six major antimalarial drugs and whether it carries structural variations that can cause RDT failure. The classification scheme is heuristic and based on a subset of known genetic markers, so it should not be treated as a failsafe predictor of the phenotype of a particular sample. Our purpose in providing these annotations is to make it easy for users without specialist training in genetics to explore the global dataset and to analyse any subset of samples for key features that are relevant to malaria control.
An important function of this curated dataset is to provide information on the provenance and key features of samples associated with each partner study, thus allowing the findings reported in different publications to be linked and compared.
Data produced by the Pf Community Project have been analysed in more than 50 publications (refs 5-55) and a few examples will serve to illustrate the diverse ways in which the data are being used. An analysis of samples collected across Africa by Amambua-Ngwa, Djimde and colleagues found evidence that parasite population structure overlaps with historical patterns of human migration and that the P. falciparum population in Ethiopia is significantly diverged from other parts of the continent 12 . A series of studies by Amato, Miotto and colleagues have documented the evolution of a multidrugresistant lineage of P. falciparum that originated in Western Cambodia over ten years ago and is now expanding rapidly across Southeast Asia, acquiring additional resistance mutations as it spreads 11,13,14 . McVean and colleagues have developed a computational method for deconvolution of the haplotypic structure of mixed infections, allowing analysis of the pedigree structure of parasites that are cotransmitted by the same mosquito 49 . Bahlo and colleagues have developed a different haplotype-based method to describe the relatedness structure of the parasite population and to identify new genomic loci with evidence of recent positive selection 50 .
A recent report from the World Health Organization highlights the need for improved surveillance systems in sustaining malaria control and achieving the long-term goal of malaria eradication 84 . To be of practical value for national malaria control programmes, genetic data must address well-defined use cases and be readily accessible 85 . Amplicon sequencing technologies provide a powerful new tool for targeted genotyping that could feasibly be implemented locally in malariaendemic countries 86,87 , but there remains a need for the international malaria control community to generate and share whole genome sequencing data, e.g. to monitor for newly emerging forms of drug resistance and to understand regional patterns of parasite migration. The next generation of long-read sequencing technologies will improve the precision of population genomic inference, e.g. by enabling analysis of hypervariable regions of the genome, and of pedigree structures within mixed infections. The accuracy with which the resistance phenotype of a sample can be predicted from genome sequencing data will also improve as we gain better functional understanding of the polygenic determinants of drug resistance.
Thus, the next few years are likely to see major advances in both the scale and information content of parasite genomic data. The practical value for malaria control will be greatly enhanced by the progressive acquisition of longitudinal time-series data, particularly if this is linked to other sources of epidemiological data and translated into reliable, actionable information with sufficient rapidity to allow control programmes to monitor the impact of their interventions on the parasite population in near real time. The Pf Community Project provides proof of concept that systems can be developed for groups in different countries to share data, to analyse it using standardised methods, and to make it readily accessible to other researchers and the malaria control community.

Methods
Here we summarise the bioinformatics methods used to produce and analyse the data; further details are available at www.malariagen.net/resource/26.

Ethical approval
All samples in this study were derived from blood samples obtained from patients with P. falciparum malaria, collected with informed consent from the patient or a parent or guardian. At each location, sample collection was approved by the appropriate local and institutional ethics committees. Standard laboratory protocols were used to determine DNA quantity and proportion of human DNA in each sample as previously described 7,56 .

Data generation and curation
Reads mapping to the human reference genome were discarded before all analyses, and the remaining reads were mapped to the P. falciparum 3D7 v3 reference genome using bwa mem 88 version 0.7.15. "Improved" BAMs were created using the Picard tools CleanSam, FixMateInformation and MarkDuplicates version 2.6.0 and GATK v3 base quality score recalibration. All lanes for each sample were merged to create sample-level BAM files.
We discovered potential SNPs and indels by running GATK's HaplotypeCaller 89 independently across each of the 7,182 sample-level BAM files and genotyped these for each of the 16 reference sequences (14 chromosomes, 1 apicoplast and 1 mitochondria) using GATK's CombineGVCFs and GenotypeGCVFs.
SNPs and indels were filtered using GATK's Variant Quality Score Recalibration (VQSR). Variants with a VQSLOD score ≤ 0 were filtered out. Functional annotations were applied using snpEff 90 version 4.1. Genome regions were annotated using vcftools version 0.1.10 and masked if they were outside the core genome. Unless otherwise specified, we used biallelic SNPs that pass all quality filters for all the analysis.
We removed 69 samples from lab studies to create the release VCF files which contain 7,113 samples. VCF files were converted to ZARR format and subsequent analyses were mainly performed using scikit-allel version 1.1.18 and the ZARR files.
We identified species using nucleotide sequence from reads mapping to six different loci in the mitochondrial genome, using custom java code (available at https://github.com/malariagen/ GeneticReportCard). The loci were located within the cox3 gene (PF3D7_MIT01400), as described in a previously published species detection method 91 . Alleles at various mitochondrial positions within the six loci were genotyped and used for classification as shown in Supplementary Table 14 (Supplementary  Data).
We created a final analysis set of 5,970 samples after removing replicate, low coverage, suspected contaminations or mislabelling and mixed-species samples.

Genotyping of drug resistance markers and samples classification
We used two complementary methods to determine tandem duplication genotypes around mdr1, plasmepsin2/3 and gch1, namely a coverage-based method and a method based on position and orientation of reads near discovered duplication breakpoints. In brief, the outline algorithm is: (1) Determine copy number at each locus using a coverage based hidden Markov model (HMM); (2) Determine breakpoints of identified duplications by manual inspection of reads and face-away read pairs around all sets of breakpoints; (3) for each locus in each sample, initially set copy number to that determined by the HMM if ≤ 10 CNVs discovered in total, else consider undetermined; (4) if faceaway pairs provide self-sufficient evidence for the presence or absence of the amplification, override the HMM call; (5) for each locus in each sample, set the breakpoint to be that with the highest proportion of face-away reads.
We genotyped deletions in hrp2 and hrp3 by manual inspection of sequence read coverage plots.
The procedure used to map genetic markers to inferred resistance status classification is described in detail for each drug in the accompanying data release (https://www.malariagen. net/resource/26).
In brief, we called amino acids at selected loci by first determining the reference amino acids and then, for each sample, applying all variations using the GT field of the VCF file. The amino acid and copy number calls generated were used to classify all samples into different types of drug resistance. Our methods of classification were heuristic and based on the available data and current knowledge of the molecular mechanisms. Each type of resistance was considered to be either present, absent or unknown for a given sample.

Population-level analysis and characterisation
We calculate genetic distance between samples using biallelic SNPs that pass filters using a method previously described 9 . In addition to calculating genetic distance between all pairs of samples from the current data set, we also calculated the genetic distance between each sample and the lab strains 3D7, 7G8, GB4, HB3 and Dd2 from the Pf3k project.
The matrix of genetic distances was used to generate neighbour-joining trees and principal coordinates. Based on these observations we grouped the samples into eight geographic regions: South America, West Africa, Central Africa, East Africa, South Asia, the western part of Southeast Asia, the eastern part of Southeast Asia and Oceania, with samples assigned to region based on the geographic location of the sampling site. Five samples from returning travellers were assigned to region based on the reported country of travel.
F WS was calculated using custom python scripts using the method previously described 7 . Nucleotide diversity (π) was calculated in non-overlapping 25 kbp genomic windows, only considering coding biallelic SNPs to reduce the ascertainment bias caused by poor accessibility of non-coding regions. LD decay (r 2 ) was calculated using the method of Rogers and Huff and biallelic SNPs with low missingness and regional allele frequency >10%. Mean F ST between populations was calculated using Hudson's method.
Allele frequencies stratified by geographic regions and sampling sites were calculated using the genotype calls produced by GATK. We defined the local differentiation score for each gene using the second highest of these six ranks (N), to ensure that the gene was highly ranked in at least two populations, i.e. to minimise the chance of artefactually ranked a gene highly due to a single variant in a single population. The final local differentiation score was normalised to ensure that the range of possible scores was between 0 and 1, local differentiation score was defined as 1 An earlier version of this article can be found on bioRxiv (DOI: https://doi.org/10.1101/824730). • Study information: Details of the 49 contributing partner studies, including description, contact information and key people.

Data availability
• Sample provenance and sequencing metadata: sample information including partner study information, location and year of collection, ENA accession numbers, and QC information for 7,113 samples from 28 countries.
• Measure of complexity of infections: characterisation of within-host diversity (FWS) for 5,970 QC pass samples.
• Inferred resistance status classification: classification of 5,970 QC pass samples into different types of resistance to 10 drugs or combinations of drugs and to RDT detection: chloroquine, pyrimethamine, sulfadoxine, mefloquine, artemisinin, piperaquine, sulfadoxine-pyrimethamine for treatment of uncomplicated malaria, sulfadoxinepyrimethamine for intermittent preventive treatment in pregnancy, artesunate-mefloquine, dihydroartemisininpiperaquine, hrp2 and hrp3 genes deletions.
• Drug resistance markers to inferred resistance status: details of the heuristics utilised to map genetic markers to resistance status classification.
• Gene differentiation: estimates of global and local differentiation for 5,561 genes.
• Short variants genotypes: Genotype calls on 6,051,696 SNPs and short indels in 7,113 samples from 29 countries, available both as VCF and zarr files.

Didier Menard
Malaria Genetics and Resistance Unit, Parasites and Insect Vectors Department, Institut Pasteur, Paris, France This manuscript from the MalariaGEN consortium, a data-sharing community of teams working on Plasmodium falciparum genomic epidemiology, presents the new release of curated P. falciparum genomes from isolates collected in 73 locations in Africa, Asia, South America and Oceania.
Based on robust and perfectly detailed methods (ranging from the treatment of the blood samples, the DNA extraction, the Illumina and computational platforms developed to produce genome sequencing for variant discovery and genotype calling), they analyzed 7000 P. falciparum genome sequences and provided numerous exciting data. For instance, they found that variations (SNPs and indels) in P. falciparum genome affected about a quarter of the 23 Mb genome (and mostly coding regions), or that duplication genotypes are frequent around mdr1, plasmepsin2/3 and gch1, which are known to be associated with antimalarial drug resistance (including mefloquine, piperaquine and sulfadoxine/pyrimethamine). Moreover, population genetic analyses conducted on this largest available data resource, depict a comprehensive picture of P. falciparum parasite populations globally and sub populations at continental level. In the results, a large section is devoted to the description of the geographic patterns of validated molecular markers (SNPs and CNVs) associated with antimalarial drug resistance. By compiling data on all samples collected from 2002-2015, they present clear profiles of drug resistance by regional sub-populations for the most used antimalarial drugs. Finally, they reveal a global landscape regarding a major challenge for malaria elimination, that are deletions in hrp2 and 3 genes linked with false negative results of HRP2-based malaria RDT.
Written in a very clear way, it must be point out that the authors have made huge efforts so that these data are understandable for a general audience, especially for the non-experts in genomics or for policy makers in malaria endemic countries. Their data effectively depict the main challenges currently encountered in the fight against malaria: the monitoring of the strategies deployed by the assessment of the impact on P. falciparum parasite populations, the geographical evolution of antimalarial drug resistances and the effectiveness of diagnostic tools used in malaria endemic areas (i.e. malaria RDT).
Of note, the authors fairly expose the main issues and drawbacks related to the methods used (i.e. the analytical challenges due to long tracts of highly repetitive sequence and hypervariable regions within the P. falciparum genome, and the challenges of studying a complex mixture of genotypes from polyclonal infections), Although, I am impressed by the work done by the consortium, I have several minor comments that could improve the manuscript: Sample collection -P. falciparum samples investigated are not from systematic sampling collections dedicated to this study but rather from multiple studies conducted by groups with different objectives and from heterogeneous populations (patients living in malaria endemic areas, travelers, etc.) . I think this issue should be discussed in the manuscript. Yes impact in rapid diagnostic tests. We do not have major criticisms of the study.

Our minor suggestions for the improval of the manuscript focus on:
Increasing the accessibility of the table listing polymorphisms in supplementary data. The