Sex differences in the genetic and molecular mechanisms of coronary artery disease.

Sex differences in coronary artery disease (CAD) presentation, risk factors and prognosis have been widely studied. Similarly, studies on atherosclerosis have shown prominent sex differences in plaque biology. Our understanding of the underlying genetic and molecular mechanisms that drive these differences remains fragmented and largely understudied. Through reviewing genetic and epigenetic studies, we identified more than 40 sex-differential candidate genes (13 within known CAD loci) that may explain, at least in part, sex differences in vascular remodeling, lipid metabolism and endothelial dysfunction. Studies with transcriptomic and single-cell RNA sequencing data from atherosclerotic plaques highlight potential sex differences in smooth muscle cell and endothelial cell biology. Especially, phenotypic switching of smooth muscle cells seems to play a crucial role in female atherosclerosis. This matches the known sex differences in atherosclerotic phenotypes, with men being more prone to lipid-rich plaques, while women are more likely to develop fibrous plaques with endothelial dysfunction. To unravel the complex mechanisms that drive sex differences in CAD, increased statistical power and adjustments to study designs and analysis strategies are required. This entails increasing inclusion rates of women, performing well-defined sex-stratified analyses and the integration of multi-omics data.

Female Gene Networks Are Expressed in Myofibroblast-Like Smooth Muscle Cells in Vulnerable Atherosclerotic Plaques.

BACKGROUND: Women presenting with coronary artery disease more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. METHODS: Gene regulatory networks were created using RNAseq gene expression data from human carotid atherosclerotic plaques. The networks were prioritized based on sex bias, relevance for smooth muscle biology, and coronary artery disease genetic enrichment. Network expression was linked to histologically determined plaque phenotypes. In addition, their expression in plaque cell types was studied at single-cell resolution using single-cell RNAseq. Finally, their relevance for disease progression was studied in female and male Apoe-/- mice fed a Western diet for 18 and 30 weeks. RESULTS: Here, we identify multiple sex-stratified gene regulatory networks from human carotid atherosclerotic plaques. Prioritization of the female networks identified 2 main SMC gene regulatory networks in late-stage atherosclerosis. Single-cell RNA sequencing mapped these female networks to 2 SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like network was mostly expressed in plaques that were vulnerable in women. Finally, the mice ortholog of key driver gene MFGE8 (milk fat globule EGF and factor V/VIII domain containing) showed retained expression in advanced plaques from female mice but was downregulated in male mice during atherosclerosis progression. CONCLUSIONS: Female atherosclerosis is characterized by gene regulatory networks that are active in fibrous vulnerable plaques rich in myofibroblast-like SMCs.

Plasma Proteomic Patterns Show Sex Differences in Early Concentric Left Ventricular Remodeling.

BACKGROUND: Concentric remodeling (cRM) can precede heart failure with preserved ejection fraction (HFpEF), a condition prevalent in women. METHODS: Patients (n=60 593, 54.2% women) visiting outpatient clinics of Cardiology Centers of the Netherlands were analyzed for cRM, HFpEF development, and mortality risk. We studied risk factors for relative wall thickness both sex-stratified and in women and men combined. Biomarker profiling was performed (4534 plasma proteins) in a substudy involving 557 patients (65.4% women) to identify pathways involved in cRM. RESULTS: cRM was present in 23.5% of women and 27.6% of men and associated with developing HFpEF (HR, 2.15 [95% CI, 1.51-2.99]) and mortality risk (HR, 1.09 [95% CI, 1.00-1.19]) in both sexes. Age, heart rate, and hypertension were statistically significantly stronger risk factors for relative wall thickness in women than men. Higher circulating levels of IFNA5 (interferon alpha-5) were associated with higher relative wall thickness in women only. Pathway analysis revealed differential pathway activation by sex and increased expression of inflammatory pathways in women. CONCLUSIONS: cRM is prevalent in approximately 1 in 4 women and men visiting outpatient cardiology clinics and associated with HFpEF development and mortality risk in both sexes. Known risk factors for cRM were more strongly associated in women than men. Proteomic analysis revealed inflammatory pathway activation in women, with a central role for IFNA5. Differential biologic pathway activation by sex in cRM may contribute to the female predominance of HFpEF and holds promise for identification of new therapeutic avenues for prevention and treatment of HFpEF. REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier: NCT001747.

Deep neural networks reveal novel sex-specific electrocardiographic features relevant for mortality risk.

Aims: Incorporation of sex in study design can lead to discoveries in medical research. Deep neural networks (DNNs) accurately predict sex based on the electrocardiogram (ECG) and we hypothesized that misclassification of sex is an important predictor for mortality. Therefore, we first developed and validated a DNN that classified sex based on the ECG and investigated the outcome. Second, we studied ECG drivers of DNN-classified sex and mortality. Methods and results: A DNN was trained to classify sex based on 131 673 normal ECGs. The algorithm was validated on internal (68 500 ECGs) and external data sets (3303 and 4457 ECGs). The survival of sex (mis)classified groups was investigated using time-to-event analysis and sex-stratified mediation analysis of ECG features. The DNN successfully distinguished female from male ECGs {internal validation: area under the curve (AUC) 0.96 [95% confidence interval (CI): 0.96, 0.97]; external validations: AUC 0.89 (95% CI: 0.88, 0.90), 0.94 (95% CI: 0.93, 0.94)}. Sex-misclassified individuals (11%) had a 1.4 times higher mortality risk compared with correctly classified peers. The ventricular rate was the strongest mediating ECG variable (41%, 95% CI: 31%, 56%) in males, while the maximum amplitude of the ST segment was strongest in females (18%, 95% CI: 11%, 39%). Short QRS duration was associated with higher mortality risk. Conclusion: Deep neural networks accurately classify sex based on ECGs. While the proportion of ECG-based sex misclassifications is low, it is an interesting biomarker. Investigation of the causal pathway between misclassification and mortality uncovered new ECG features that might be associated with mortality. Increased emphasis on sex as a biological variable in artificial intelligence is warranted.

Plasma Testosterone Levels and Atherosclerotic Plaque Gene Expression in Men With Advanced Atherosclerosis.

Aims: Low plasma testosterone levels have been shown to predict worse outcome in men with severe atherosclerotic disease. We hypothesized that a low plasma testosterone level affects disease risk through changes in gene expression in atherosclerotic plaques. Therefore, we studied plasma testosterone levels in relation to gene expression levels in atherosclerotic plaque tissue of men with advanced atherosclerotic disease. Methods: Plasma testosterone levels were measured in 203 men undergoing carotid endarterectomy. The corresponding atherosclerotic plaque tissue was used for RNA sequencing. First, we assessed how often the androgen receptor gene was expressed in the plaque. Second, correlations between plasma testosterone levels and pre-selected testosterone-sensitive genes were assessed. Finally, differences within the RNA expression profile of the plaque as a whole, characterized into gene regulatory networks and at individual gene level were assessed in relation to testosterone levels. Results: Testosterone plasma levels were low with a median of 11.6 nmol/L (IQR: 8.6-13.8). RNA-seq of the plaque resulted in reliable expression data for 18,850 genes to be analyzed. Within the RNA seq data, the androgen-receptor gene was expressed in 189 out of 203 (93%) atherosclerotic plaques of men undergoing carotid endarterectomy. The androgen receptor gene expression was not associated with testosterone plasma levels. There were no significant differences in gene expression of atherosclerotic plaques between different endogenous testosterone levels. This remained true for known testosterone-sensitive genes, the complete transcriptomic profile, male-specific gene co-expression modules as well as for individual genes. Conclusion: In men with severe atherosclerotic disease the androgen receptor is highly expressed in plaque tissue. However, plasma testosterone levels were neither associated with pre-selected testosterone sensitive genes, gene expression profiles nor gene regulatory networks in late-stage atherosclerotic plaques. The effect of testosterone on gene expression of the late-stage atherosclerotic plaque appears limited, suggesting that alternate mechanisms explain its effect on clinical outcomes.

Genetic diversity of the Plasmodium falciparum GTP-cyclohydrolase 1, dihydrofolate reductase and dihydropteroate synthetase genes reveals new insights into sulfadoxine-pyrimethamine antimalarial drug resistance.

Plasmodium falciparum parasites resistant to antimalarial treatments have hindered malaria disease control. Sulfadoxine-pyrimethamine (SP) was used globally as a first-line treatment for malaria after wide-spread resistance to chloroquine emerged and, although replaced by artemisinin combinations, is currently used as intermittent preventive treatment of malaria in pregnancy and in young children as part of seasonal malaria chemoprophylaxis in sub-Saharan Africa. The emergence of SP-resistant parasites has been predominantly driven by cumulative build-up of mutations in the dihydrofolate reductase (pfdhfr) and dihydropteroate synthetase (pfdhps) genes, but additional amplifications in the folate pathway rate-limiting pfgch1 gene and promoter, have recently been described. However, the genetic make-up and prevalence of those amplifications is not fully understood. We analyse the whole genome sequence data of 4,134 P. falciparum isolates across 29 malaria endemic countries, and reveal that the pfgch1 gene and promoter amplifications have at least ten different forms, occurring collectively in 23% and 34% in Southeast Asian and African isolates, respectively. Amplifications are more likely to be present in isolates with a greater accumulation of pfdhfr and pfdhps substitutions (median of 1 additional mutations; P<0.00001), and there was evidence that the frequency of pfgch1 variants may be increasing in some African populations, presumably under the pressure of SP for chemoprophylaxis and anti-folate containing antibiotics used for the treatment of bacterial infections. The selection of P. falciparum with pfgch1 amplifications may enhance the fitness of parasites with pfdhfr and pfdhps substitutions, potentially threatening the efficacy of this regimen for prevention of malaria in vulnerable groups. Our work describes new pfgch1 amplifications that can be used to inform the surveillance of SP drug resistance, its prophylactic use, and future experimental work to understand functional mechanisms.

Selective whole genome amplification of Plasmodium malariae DNA from clinical samples reveals insights into population structure.

The genomic diversity of Plasmodium malariae malaria parasites is understudied, partly because infected individuals tend to present with low parasite densities, leading to difficulties in obtaining sufficient parasite DNA for genome analysis. Selective whole genome amplification (SWGA) increases the relative levels of pathogen DNA in a clinical sample, but has not been adapted for P. malariae parasites. Here we design customized SWGA primers which successfully amplify P. malariae DNA extracted directly from unprocessed clinical blood samples obtained from patients with P. malariae-mono-infections from six countries, and further test the efficacy of SWGA on mixed infections with other Plasmodium spp. SWGA enables the successful whole genome sequencing of samples with low parasite density (i.e. one sample with a parasitaemia of 0.0064% resulted in 44% of the genome covered by ≥ 5 reads), leading to an average 14-fold increase in genome coverage when compared to unamplified samples. We identify a total of 868,476 genome-wide SNPs, of which 194,709 are unique across 18 high-quality isolates. After exclusion of the hypervariable subtelomeric regions, a high-quality core subset of 29,899 unique SNPs is defined. Population genetic analysis suggests that P. malariae parasites display clear geographical separation by continent. Further, SWGA successfully amplifies genetic regions of interest such as orthologs of P. falciparum drug resistance-associated loci (Pfdhfr, Pfdhps, Pfcrt, Pfk13 and Pfmdr1), and several non-synonymous SNPs were detected in these genes. In conclusion, we have established a robust SWGA approach that can assist whole genome sequencing of P. malariae, and thereby facilitate the implementation of much-needed large-scale multi-population genomic studies of this neglected malaria parasite. As demonstrated in other Plasmodia, such genetic diversity studies can provide insights into the biology underlying the disease and inform malaria surveillance and control measures.

Discordant bioinformatic predictions of antimicrobial resistance from whole-genome sequencing data of bacterial isolates: an inter-laboratory study.

Antimicrobial resistance (AMR) poses a threat to public health. Clinical microbiology laboratories typically rely on culturing bacteria for antimicrobial-susceptibility testing (AST). As the implementation costs and technical barriers fall, whole-genome sequencing (WGS) has emerged as a 'one-stop' test for epidemiological and predictive AST results. Few published comparisons exist for the myriad analytical pipelines used for predicting AMR. To address this, we performed an inter-laboratory study providing sets of participating researchers with identical short-read WGS data from clinical isolates, allowing us to assess the reproducibility of the bioinformatic prediction of AMR between participants, and identify problem cases and factors that lead to discordant results. We produced ten WGS datasets of varying quality from cultured carbapenem-resistant organisms obtained from clinical samples sequenced on either an Illumina NextSeq or HiSeq instrument. Nine participating teams ('participants') were provided these sequence data without any other contextual information. Each participant used their choice of pipeline to determine the species, the presence of resistance-associated genes, and to predict susceptibility or resistance to amikacin, gentamicin, ciprofloxacin and cefotaxime. We found participants predicted different numbers of AMR-associated genes and different gene variants from the same clinical samples. The quality of the sequence data, choice of bioinformatic pipeline and interpretation of the results all contributed to discordance between participants. Although much of the inaccurate gene variant annotation did not affect genotypic resistance predictions, we observed low specificity when compared to phenotypic AST results, but this improved in samples with higher read depths. Had the results been used to predict AST and guide treatment, a different antibiotic would have been recommended for each isolate by at least one participant. These challenges, at the final analytical stage of using WGS to predict AMR, suggest the need for refinements when using this technology in clinical settings. Comprehensive public resistance sequence databases, full recommendations on sequence data quality and standardization in the comparisons between genotype and resistance phenotypes will all play a fundamental role in the successful implementation of AST prediction using WGS in clinical microbiology laboratories.

A molecular barcode to inform the geographical origin and transmission dynamics of Plasmodium vivax malaria.

Although Plasmodium vivax parasites are the predominant cause of malaria outside of sub-Saharan Africa, they not always prioritised by elimination programmes. P. vivax is resilient and poses challenges through its ability to re-emerge from dormancy in the human liver. With observed growing drug-resistance and the increasing reports of life-threatening infections, new tools to inform elimination efforts are needed. In order to halt transmission, we need to better understand the dynamics of transmission, the movement of parasites, and the reservoirs of infection in order to design targeted interventions. The use of molecular genetics and epidemiology for tracking and studying malaria parasite populations has been applied successfully in P. falciparum species and here we sought to develop a molecular genetic tool for P. vivax. By assembling the largest set of P. vivax whole genome sequences (n = 433) spanning 17 countries, and applying a machine learning approach, we created a 71 SNP barcode with high predictive ability to identify geographic origin (91.4%). Further, due to the inclusion of markers for within population variability, the barcode may also distinguish local transmission networks. By using P. vivax data from a low-transmission setting in Malaysia, we demonstrate the potential ability to infer outbreak events. By characterising the barcoding SNP genotypes in P. vivax DNA sourced from UK travellers (n = 132) to ten malaria endemic countries predominantly not used in the barcode construction, we correctly predicted the geographic region of infection origin. Overall, the 71 SNP barcode outperforms previously published genotyping methods and when rolled-out within new portable platforms, is likely to be an invaluable tool for informing targeted interventions towards elimination of this resilient human malaria.

Rapid and iterative genome editing in the malaria parasite Plasmodium knowlesi provides new tools for P. vivax research.

Tackling relapsing Plasmodium vivax and zoonotic Plasmodium knowlesi infections is critical to reducing malaria incidence and mortality worldwide. Understanding the biology of these important and related parasites was previously constrained by the lack of robust molecular and genetic approaches. Here, we establish CRISPR-Cas9 genome editing in a culture-adapted P. knowlesi strain and define parameters for optimal homology-driven repair. We establish a scalable protocol for the production of repair templates by PCR and demonstrate the flexibility of the system by tagging proteins with distinct cellular localisations. Using iterative rounds of genome-editing we generate a transgenic line expressing P. vivax Duffy binding protein (PvDBP), a lead vaccine candidate. We demonstrate that PvDBP plays no role in reticulocyte restriction but can alter the macaque/human host cell tropism of P. knowlesi. Critically, antibodies raised against the P. vivax antigen potently inhibit proliferation of this strain, providing an invaluable tool to support vaccine development.

Artemisinin resistance-associated markers in Plasmodium falciparum parasites from the China-Myanmar border: predicted structural stability of K13 propeller variants detected in a low-prevalence area.

BACKGROUND: Malaria reduction and future elimination in China is made more difficult by the importation of cases from neighboring endemic countries, particularly Myanmar, Laos, and Vietnam, and increased travel to Africa by Chinese nationals. The increasing prevalence of artemisinin resistant parasites across Southeast Asia highlights the importance of monitoring the parasite importation into China. Artemisinin resistance in the Mekong region is associated with variants of genes encoding the K13 kelch domain protein (pf13k), found in specific genetic backgrounds, including certain alleles of genes encoding the chloroquine resistance transporter (pfcrt) and multidrug resistance transporter PgH1 (pfmdr1). METHODS: In this study we investigated the prevalence of drug resistance markers in 72 P. falciparum samples from uncomplicated malaria infections in Tengchong and Yingjiang, counties on the Yunnan-Myanmar border. Variants of pf13k, pfcrt and pfmdr1 are described. RESULTS: Almost all parasites harboured chloroquine-resistant alleles of pfcrt, whereas pfmdr1 was more diverse. Major mutations in the K13 propeller domain associated with artemisinin resistance in the Mekong region (C580Y, R539T and Y493H) were absent, but F446I and two previously undescribed mutations (V603E and V454I) were identified. Protein structural modelling was carried out in silico on each of these K13 variants, based on recently published crystal structures for the K13 propeller domain. Whereas F446I was predicted to elicit a moderate destabilisation of the propeller structure, the V603E substitution is likely to lead to relatively high protein instability. We plotted these stability estimates, and those for all previously described variants, against published values for in vivo parasitaemia half-life, and found that quadratic regression generates a useful predictive algorithm. CONCLUSION: This study provides a baseline of P. falciparum resistance-associated mutations prevalent at the China-Myanmar border. We also show that protein modelling can be used to generate testable predictions as to the impact of pfk13 mutations on in vivo (and potentially in vitro) artemisinin susceptibility.

Know Your Heart: Rationale, design and conduct of a cross-sectional study of cardiovascular structure, function and risk factors in 4500 men and women aged 35-69 years from two Russian cities, 2015-18.

Russia has one of the highest rates of cardiovascular disease in the world. The International Project on Cardiovascular Disease in Russia (IPCDR) was set up to understand the reasons for this. A substantial component of this study was the Know Your Heart Study devoted to characterising the nature and causes of cardiovascular disease in Russia by conducting large cross-sectional surveys in two Russian cities Novosibirsk and Arkhangelsk. The study population was 4542 men and women aged 35-69 years recruited from the general population. Fieldwork took place between 2015-18. There were two study components: 1) a baseline interview to collect information on socio-demographic characteristics and cardiovascular risk factors, usually conducted at home, and 2) a comprehensive health check at a primary care clinic which included detailed examination of the cardiovascular system. In this paper we describe in detail the rationale for, design and conduct of these studies.

Global analysis of Plasmodium falciparum histidine-rich protein-2 (pfhrp2) and pfhrp3 gene deletions using whole-genome sequencing data and meta-analysis.

Many rapid diagnostic tests (RDT) used on suspected malaria cases are based on the detection of the protein encoded by the Plasmodium falciparum histidine-rich protein-2 (pfhrp2) gene, which shares a high sequence homology with pfhrp3 in the 3D7 reference genome. Parasite isolates showing pfhrp2 and pfhrp3 gene deletions have been emerging over the years, but a comprehensive genetic analysis of these variants is still lacking. With this purpose, genomic data from experimental P. falciparum genetic crosses between different laboratory lines (3D7, HB3, DD2, 7G8 and GB4) were first analysed (n = 98). The frequency of pfhrp2 deletions was consistent with a Mendelian prediction in HB3 × DD2 (56.7%; 95%CI = (39.5%-72.9%)). Moreover, the pfhrp2 and pfhrp3 deletions segregated independently of each other in the same genetic cross. Analysis of 3D7 × HB3 and 7G8 × GB4 estimated the probability of spontaneously generating a pfhrp2 deletion during sexual recombination to be up to 6.2%. Next, whole genome sequence data from 1970 P. falciparum isolates collected globally were analysed. Nine samples displayed depth of coverage consistent with pfhrp2 deletions (0.5%), but the corresponding split-read analysis could not confirm deletions in seven of these samples. Twenty-eight isolates had evidence of pfhrp3 deletions (1.4%), which are widespread in Southeast Asia. Finally, a meta-analysis of published data revealed a positive mean association between the frequencies of pfhrp2 and pfhrp3 deletions in Africa and South America. This result suggested a shared selective pressure acting on these genetic variants. In conclusion, evidence of genetic selection on both pfhrp2 and pfhrp3 deletions was presented, but experimental crosses do not provide evidence of a fitness cost of these variants. Further work is urgently needed to accurately determine the prevalence and the degree of association between these genetic variants, and the respective impact on diagnostic accuracy of many in-use RDT.

Enrichment of Clinically Relevant Organisms in Spontaneous Preterm-Delivered Placentas and Reagent Contamination across All Clinical Groups in a Large Pregnancy Cohort in the United Kingdom.

In this study, differences in the placental microbiota from term and preterm deliveries in a large pregnancy cohort in the United Kingdom were studied by using 16S-targeted amplicon sequencing. The impacts of contamination from DNA extraction, PCR reagents, and the delivery itself were also examined. A total of 400 placental samples from 256 singleton pregnancies were analyzed, and differences between spontaneous preterm-, nonspontaneous preterm-, and term-delivered placentas were investigated. DNA from recently delivered placentas was extracted, and screening for bacterial DNA was carried out by using targeted sequencing of the 16S rRNA gene on the Illumina MiSeq platform. Sequenced reads were analyzed for the presence of contaminating operational taxonomic units (OTUs) identified via sequencing of negative extraction and PCR-blank samples. Differential abundances and between-sample (beta) diversity metrics were then compared. A large proportion of the reads sequenced from the extracted placental samples mapped to OTUs that were also found for negative extractions. Striking differences in the compositions of samples were also observed, according to whether the placenta was delivered abdominally or vaginally, providing strong circumstantial evidence for delivery contamination as an important contributor to observed microbial profiles. When OTU- and genus-level abundances were compared between the groups of interest, a number of organisms were enriched in the spontaneous preterm-delivery cohort, including organisms that have been associated previously with adverse pregnancy outcomes, specifically Mycoplasma spp. and Ureaplasma spp. However, analyses of the overall community structure did not reveal convincing evidence for the existence of a reproducible "preterm placental microbiome."IMPORTANCE Preterm birth is associated with both psychological and physical disabilities and is the leading cause of infant morbidity and mortality worldwide. Infection is known to be an important cause of spontaneous preterm birth, and recent research has implicated variation in the "placental microbiome" in the risk of preterm birth. Consistent with data from previous studies, the abundances of certain clinically relevant species differed between spontaneous preterm- and nonspontaneous preterm- or term-delivered placentas. These results support the view that a proportion of spontaneous preterm births have an intrauterine-infection component. However, an additional observation from this study was that a substantial proportion of sequenced reads were contaminating reads rather than DNA from endogenous, clinically relevant species. This observation warrants caution in the interpretation of sequencing outputs from low-biomass samples such as the placenta.

Analysis of nuclear and organellar genomes of Plasmodium knowlesi in humans reveals ancient population structure and recent recombination among host-specific subpopulations.

The macaque parasite Plasmodium knowlesi is a significant concern in Malaysia where cases of human infection are increasing. Parasites infecting humans originate from genetically distinct subpopulations associated with the long-tailed (Macaca fascicularis (Mf)) or pig-tailed macaques (Macaca nemestrina (Mn)). We used a new high-quality reference genome to re-evaluate previously described subpopulations among human and macaque isolates from Malaysian-Borneo and Peninsular-Malaysia. Nuclear genomes were dimorphic, as expected, but new evidence of chromosomal-segment exchanges between subpopulations was found. A large segment on chromosome 8 originating from the Mn subpopulation and containing genes encoding proteins expressed in mosquito-borne parasite stages, was found in Mf genotypes. By contrast, non-recombining organelle genomes partitioned into 3 deeply branched lineages, unlinked with nuclear genomic dimorphism. Subpopulations which diverged in isolation have re-connected, possibly due to deforestation and disruption of wild macaque habitats. The resulting genomic mosaics reveal traits selected by host-vector-parasite interactions in a setting of ecological transition.

Modest heterologous protection after Plasmodium falciparum sporozoite immunization: a double-blind randomized controlled clinical trial.

BACKGROUND: A highly efficacious vaccine is needed for malaria control and eradication. Immunization with Plasmodium falciparum NF54 parasites under chemoprophylaxis (chemoprophylaxis and sporozoite (CPS)-immunization) induces the most efficient long-lasting protection against a homologous parasite. However, parasite genetic diversity is a major hurdle for protection against heterologous strains. METHODS: We conducted a double-blind, randomized controlled trial in 39 healthy participants of NF54-CPS immunization by bites of 45 NF54-infected (n = 24 volunteers) or uninfected mosquitoes (placebo; n = 15 volunteers) against a controlled human malaria infection with the homologous NF54 or the genetically distinct NF135.C10 and NF166.C8 clones. Cellular and humoral immune assays were performed as well as genetic characterization of the parasite clones. RESULTS: NF54-CPS immunization induced complete protection in 5/5 volunteers against NF54 challenge infection at 14 weeks post-immunization, but sterilely protected only 2/10 and 1/9 volunteers against NF135.C10 and NF166.C8 challenge infection, respectively. Post-immunization plasma showed a significantly lower capacity to block heterologous parasite development in primary human hepatocytes compared to NF54. Whole genome sequencing showed that NF135.C10 and NF166.C8 have amino acid changes in multiple antigens targeted by CPS-induced antibodies. Volunteers protected against heterologous challenge were among the stronger immune responders to in vitro parasite stimulation. CONCLUSIONS: Although highly protective against homologous parasites, NF54-CPS-induced immunity is less effective against heterologous parasite clones both in vivo and in vitro. Our data indicate that whole sporozoite-based vaccine approaches require more potent immune responses for heterologous protection. TRIAL REGISTRATION: This trial is registered in clinicaltrials.gov, under identifier NCT02098590 .

Genomic variation in Plasmodium vivax malaria reveals regions under selective pressure.

BACKGROUND: Although Plasmodium vivax contributes to almost half of all malaria cases outside Africa, it has been relatively neglected compared to the more deadly P. falciparum. It is known that P. vivax populations possess high genetic diversity, differing geographically potentially due to different vector species, host genetics and environmental factors. RESULTS: We analysed the high-quality genomic data for 46 P. vivax isolates spanning 10 countries across 4 continents. Using population genetic methods we identified hotspots of selection pressure, including the previously reported MRP1 and DHPS genes, both putative drug resistance loci. Extra copies and deletions in the promoter region of another drug resistance candidate, MDR1 gene, and duplications in the Duffy binding protein gene (PvDBP) potentially involved in erythrocyte invasion, were also identified. For surveillance applications, continental-informative markers were found in putative drug resistance loci, and we show that organellar polymorphisms could classify P. vivax populations across continents and differentiate between Plasmodia spp. CONCLUSIONS: This study has shown that genomic diversity that lies within and between P. vivax populations can be used to elucidate potential drug resistance and invasion mechanisms, as well as facilitate the molecular barcoding of the parasite for surveillance applications.