Vitamin C facilitates direct cardiac reprogramming by inhibiting reactive oxygen species.

BACKGROUND: After myocardial infarction, the lost myocardium is replaced by fibrotic tissue, eventually progressively leading to myocardial dysfunction. Direct reprogramming of fibroblasts into cardiomyocytes via the forced overexpression of cardiac transcription factors Gata4, Mef2c, and Tbx5 (GMT) offers a promising strategy for cardiac repair. The limited reprogramming efficiency of this approach, however, remains a significant challenge. METHODS: We screened seven factors capable of improving direct cardiac reprogramming of both mice and human fibroblasts by evaluating small molecules known to be involved in cardiomyocyte differentiation or promoting human-induced pluripotent stem cell reprogramming. RESULTS: We found that vitamin C (VitC) significantly increased cardiac reprogramming efficiency when added to GMT-overexpressing fibroblasts from human and mice in 2D and 3D model. We observed a significant increase in reactive oxygen species (ROS) generation in human and mice fibroblasts upon Doxy induction, and ROS generation was subsequently reduced upon VitC treatment, associated with increased reprogramming efficiency. However, upon treatment with dehydroascorbic acid, a structural analog of VitC but lacking antioxidant properties, no difference in reprogramming efficiency was observed, suggesting that the effect of VitC in enhancing cardiac reprogramming is partly dependent of its antioxidant properties. CONCLUSIONS: Our findings demonstrate that VitC supplementation significantly enhances the efficiency of cardiac reprogramming, partially by suppressing ROS production in the presence of GMT.

Genomic analysis of Plasmodium vivax describes patterns of connectivity and putative drivers of adaptation in Ethiopia.

Ethiopia has the greatest burden of Plasmodium vivax in Africa, but little is known about the epidemiological landscape of parasites across the country. We analysed the genomic diversity of 137 P. vivax isolates collected nine Ethiopian districts from 2012 to 2016. Signatures of selection were detected by cross-country comparisons with isolates from Thailand (n = 104) and Indonesia (n = 111), representing regions with low and high chloroquine resistance respectively. 26% (35/137) of Ethiopian infections were polyclonal, and 48.5% (17/35) of these comprised highly related clones (within-host identity-by-descent > 25%), indicating frequent co-transmission and superinfection. Parasite gene flow between districts could not be explained entirely by geographic distance, with economic and cultural factors hypothesised to have an impact on connectivity. Amplification of the duffy binding protein gene (pvdbp1) was prevalent across all districts (16-75%). Cross-population haplotype homozygosity revealed positive selection in a region proximal to the putative chloroquine resistance transporter gene (pvcrt-o). An S25P variant in amino acid transporter 1 (pvaat1), whose homologue has recently been implicated in P. falciparum chloroquine resistance evolution, was prevalent in Ethiopia (96%) but not Thailand or Indonesia (35-53%). The genomic architecture in Ethiopia highlights circulating variants of potential public health concern in an endemic setting with evidence of stable transmission.

Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas.

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

Human Plaque Myofibroblasts to Study Mechanisms of Atherosclerosis.

Background Plaque myofibroblasts are critical players in the initiation and advancement of atherosclerotic disease. They are involved in the production of extracellular matrix, the formation of the fibrous cap, and the underlying lipidic core via modulation processes in response to different environmental cues. Despite clear phenotypic differences between myofibroblast cells and healthy vascular smooth muscle cells, smooth muscle cells are still widely used as a cellular model in atherosclerotic research. Methods and Results Here, we present a conditioned outgrowth method to isolate and culture myofibroblast cells from plaques. We obtained these cells from 27 donors (24 carotid and 3 femoral endarterectomies). We show that they keep their proliferative capacity for 8 passages, are transcriptionally stable, retain donor-specific gene expression programs, and express extracellular matrix proteins (FN1, COL1A1, and DCN) and smooth muscle cell markers (ACTA2, MYH11, and CNN1). Single-cell transcriptomics reveals that the cells in culture closely resemble the plaque myofibroblasts. Chromatin immunoprecipitation sequencing shows the presence of histone H3 lysine 4 dimethylation at the MYH11 promoter, pointing to their smooth muscle cell origin. Finally, we demonstrated that plaque myofibroblasts can be efficiently transduced (>97%) and are capable of taking up oxidized low-density lipoprotein and undergoing calcification. Conclusions In conclusion, we present a method to isolate and culture cells that retain plaque myofibroblast phenotypical and functional capabilities, making them a suitable in vitro model for studying selected mechanisms of atherosclerosis.

Multi-ancestry genome-wide study identifies effector genes and druggable pathways for coronary artery calcification.

Coronary artery calcification (CAC), a measure of subclinical atherosclerosis, predicts future symptomatic coronary artery disease (CAD). Identifying genetic risk factors for CAC may point to new therapeutic avenues for prevention. Currently, there are only four known risk loci for CAC identified from genome-wide association studies (GWAS) in the general population. Here we conducted the largest multi-ancestry GWAS meta-analysis of CAC to date, which comprised 26,909 individuals of European ancestry and 8,867 individuals of African ancestry. We identified 11 independent risk loci, of which eight were new for CAC and five had not been reported for CAD. These new CAC loci are related to bone mineralization, phosphate catabolism and hormone metabolic pathways. Several new loci harbor candidate causal genes supported by multiple lines of functional evidence and are regulators of smooth muscle cell-mediated calcification ex vivo and in vitro. Together, these findings help refine the genetic architecture of CAC and extend our understanding of the biological and potential druggable pathways underlying CAC.

X chromosome inactivation skewing is common in advanced carotid atherosclerotic lesions in females and predicts secondary peripheral artery events.

BACKGROUND AND AIM: Sex differences in atherosclerosis have been described with female plaques being mostly perceived as stable and fibrous. Sex-specific mechanisms such as mosaic loss of the Y chromosome in men have been linked to cardiovascular health. In women, X-linked mechanisms such as X chromosome inactivation (XCI) skewing is common in several tissues. Yet, information on the role of XCI in female atherosclerotic plaques is lacking. Here, we investigated the presence of XCI skewing in advanced atherosclerotic lesions and its association with cardiovascular risk factors, histological plaque data, and clinical data. METHODS: XCI skewing was quantified in 154 atherosclerotic plaque and 55 blood DNA samples of women included in the Athero-Express study. The skewing status was determined performing the HUMARA assay. Then, we studied the relationship of XCI skewing in female plaque and cardiovascular risk factors using regression models. In addition, we studied if plaque XCI predicted plaque composition, and adverse events during 3-years follow-up using Cox proportional hazard models. RESULTS: XCI skewing was detected in 76 of 154 (49.4%) plaques and in 27 of 55 (67%) blood samples. None of the clinical risk factors were associated with plaque skewing. Plaque skewing was more often detected in plaques with a plaque hemorrhage (OR [95% CI]: 1.44 [1.06-1.98], P = 0.02). Moreover, skewed plaques were not associated with a higher incidence of composite and major events but were specifically associated with peripheral artery events during a 3-year follow-up period in a multivariate model (HR [95%CI]: 1.46 [1.09-1.97]; P = 0.007). CONCLUSIONS: XCI skewing is common in carotid plaques of females and is predictive for the occurrence of peripheral artery events within 3 years after carotid endarterectomy.

Sex-differences have been observed in the development of atherosclerosis between men and women. Women tend to have more stable and fibrous plaques compared to men. Sex-specific mechanisms such as mosaic loss of the Y chromosome in men, were associated with cardiovascular health. In women, despite X-linked mechanisms like X chromosome inactivation (XCI) skewing was identified in various tissues. However, its relationship with atherosclerosis has not yet been investigated. In our study, we explored if prevalence of XCI skewing in advanced atherosclerotic lesions related to cardiovascular risk factors, histological plaque data, and clinical information. We found that XCI skewing was present in approximately 50% of human plaques, particularly those with plaque hemorrhage. Interestingly, we did not find any notable relationship between plaque skewing and clinical risk factors. However, we found that XCI was more present in women with peripheral artery events during the 3 years period following carotid endarterectomy. In summary, our findings indicate that XCI skewing is commonly observed in carotid plaques among females and may serve as a predictive factor for the occurrence of peripheral artery events within 3 years after carotid endarterectomy.

Population-based genomic study of Plasmodium vivax malaria in seven Brazilian states and across South America.

Background: Brazil is a unique and understudied setting for malaria, with complex foci of transmission associated with human and environmental conditions. An understanding of the population genomic diversity of P. vivax parasites across Brazil can support malaria control strategies. Methods: Through whole genome sequencing of P. vivax isolates across 7 Brazilian states, we use population genomic approaches to compare genetic diversity within country (n = 123), continent (6 countries, n = 315) and globally (26 countries, n = 885). Findings: We confirm that South American isolates are distinct, have more ancestral populations than the other global regions, with differentiating mutations in genes under selective pressure linked to antimalarial drugs (pvmdr1, pvdhfr-ts) and mosquito vectors (pvcrmp3, pvP45/48, pvP47). We demonstrate Brazil as a distinct parasite population, with signals of selection including ABC transporter (PvABCI3) and PHIST exported proteins. Interpretation: Brazil has a complex population structure, with evidence of P. simium infections and Amazonian parasites separating into multiple clusters. Overall, our work provides the first Brazil-wide analysis of P. vivax population structure and identifies important mutations, which can inform future research and control measures. Funding: AI is funded by an MRC LiD PhD studentship. TGC is funded by the Medical Research Council (Grant no. MR/M01360X/1, MR/N010469/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1). SC is funded by Medical Research Council UK grants (MR/M01360X/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1) and Bloomsbury SET (ref. CCF17-7779). FN is funded by The Shloklo Malaria Research Unit - part of the Mahidol Oxford Research Unit, supported by the Wellcome Trust (Grant no. 220211). ARSB is funded by São Paulo Research Foundation - FAPESP (Grant no. 2002/09546-1). RLDM is funded by Brazilian National Council for Scientific and Technological Development - CNPq (Grant no. 302353/2003-8 and 471605/2011-5); CRFM is funded by FAPESP (Grant no. 2020/06747-4) and CNPq (Grant no. 302917/2019-5 and 408636/2018-1); JGD is funded by FAPESP fellowships (2016/13465-0 and 2019/12068-5) and CNPq (Grant no. 409216/2018-6).

Female gene networks are expressed in myofibroblast-like smooth muscle cells in vulnerable atherosclerotic plaques.

Women presenting with coronary artery disease (CAD) 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. Here, we show sex-stratified gene regulatory networks (GRNs) from human carotid atherosclerotic tissue. Prioritization of these networks identified two main SMC GRNs in late-stage atherosclerosis. Single-cell RNA-sequencing mapped these GRNs to two SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like GRN was mostly expressed in plaques that were vulnerable in females. Finally, mice orthologs of the female myofibroblast-like genes showed retained expression in advanced plaques from female mice but were downregulated in male mice during atherosclerosis progression. Female atherosclerosis is driven by GRNs that promote a fibrous vulnerable plaque rich in myofibroblast-like SMCs.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes.

BACKGROUND: Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes. METHODS: We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome. RESULTS: We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs. CONCLUSIONS: Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Deep learning-derived cardiovascular age shares a genetic basis with other cardiac phenotypes.

Artificial intelligence (AI)-based approaches can now use electrocardiograms (ECGs) to provide expert-level performance in detecting heart abnormalities and diagnosing disease. Additionally, patient age predicted from ECGs by AI models has shown great potential as a biomarker for cardiovascular age, where recent work has found its deviation from chronological age ("delta age") to be associated with mortality and co-morbidities. However, despite being crucial for understanding underlying individual risk, the genetic underpinning of delta age is unknown. In this work we performed a genome-wide association study using UK Biobank data (n=34,432) and identified eight loci associated with delta age ([Formula: see text]), including genes linked to cardiovascular disease (CVD) (e.g. SCN5A) and (heart) muscle development (e.g. TTN). Our results indicate that the genetic basis of cardiovascular ageing is predominantly determined by genes directly involved with the cardiovascular system rather than those connected to more general mechanisms of ageing. Our insights inform the epidemiology of CVD, with implications for preventative and precision medicine.

De Novo Assembly of Plasmodium knowlesi Genomes From Clinical Samples Explains the Counterintuitive Intrachromosomal Organization of Variant SICAvar and kir Multiple Gene Family Members.

Plasmodium knowlesi, a malaria parasite of Old World macaque monkeys, is used extensively to model Plasmodium biology. Recently, P. knowlesi was found in the human population of Southeast Asia, particularly Malaysia. P. knowlesi causes uncomplicated to severe and fatal malaria in the human host with features in common with the more prevalent and virulent malaria caused by Plasmodium falciparum. As such, P. knowlesi presents a unique opportunity to develop experimental translational model systems for malaria pathophysiology informed by clinical data from same-species human infections. Experimental lines of P. knowlesi represent well-characterized genetically stable parasites, and to maximize their utility as a backdrop for understanding malaria pathophysiology, genetically diverse contemporary clinical isolates, essentially wild-type, require comparable characterization. The Oxford Nanopore PCR-free long-read sequencing platform was used to sequence and de novo assemble P. knowlesi genomes from frozen clinical samples. The sequencing platform and assembly pipelines were designed to facilitate capturing data and describing, for the first time, P. knowlesi schizont-infected cell agglutination (SICA) var and Knowlesi-Interspersed Repeats (kir) multiple gene families in parasites acquired from nature. The SICAvar gene family members code for antigenically variant proteins analogous to the virulence-associated P. falciparum erythrocyte membrane protein (PfEMP1) multiple var gene family. Evidence presented here suggests that the SICAvar family members have arisen through a process of gene duplication, selection pressure, and variation. Highly evolving genes including PfEMP1family members tend to be restricted to relatively unstable sub-telomeric regions that drive change with core genes protected in genetically stable intrachromosomal locations. The comparable SICAvar and kir gene family members are counter-intuitively located across chromosomes. Here, we demonstrate that, in contrast to conserved core genes, SICAvar and kir genes occupy otherwise gene-sparse chromosomal locations that accommodate rapid evolution and change. The novel methods presented here offer the malaria research community not only new tools to generate comprehensive genome sequence data from small clinical samples but also new insight into the complexity of clinically important real-world parasites.

Skin microbiome alters attractiveness to Anopheles mosquitoes.

BACKGROUND: Some people produce specific body odours that make them more attractive than others to mosquitoes, and consequently are at higher risk of contracting vector-borne diseases. The skin microbiome can break down carbohydrates, fatty acids and peptides on the skin into volatiles that mosquitoes can differentiate. RESULTS: Here, we examined how skin microbiome composition of women differs in relation to level of attractiveness to Anopheles coluzzii mosquitoes, to identify volatiles in body odour and metabolic pathways associated with individuals that tend to be poorly-attractive to mosquitoes. We used behavioural assays to measure attractiveness of participants to An. coluzzii mosquitoes, 16S rRNA amplicon sequencing of the bacteria sampled from the skin and gas chromatography of volatiles in body odour. We found differences in skin microbiome composition between the poorly- and highly-attractive groups, particularly eight Amplicon Sequence Variants (ASVs) belonging to the Proteobacteria, Actinobacteria and Firmicutes phyla. Staphylococcus 2 ASVs are four times as abundant in the highly-attractive compared to poorly-attractive group. Associations were found between these ASVs and volatiles known to be attractive to Anopheles mosquitoes. Propanoic pathways are enriched in the poorly-attractive participants compared to those found to be highly-attractive. CONCLUSIONS: Our findings suggest that variation in attractiveness of people to mosquitoes is related to the composition of the skin microbiota, knowledge that could improve odour-baited traps or other next generation vector control tools.

Transcriptomic-based clustering of human atherosclerotic plaques identifies subgroups with different underlying biology and clinical presentation.

Histopathological studies have revealed key processes of atherosclerotic plaque thrombosis. However, the diversity and complexity of lesion types highlight the need for improved sub-phenotyping. Here we analyze the gene expression profiles of 654 advanced human carotid plaques. The unsupervised, transcriptome-driven clustering revealed five dominant plaque types. These plaque phenotypes were associated with clinical presentation and showed differences in cellular compositions. Validation in coronary segments showed that the molecular signature of these plaques was linked to coronary ischemia. One of the plaque types with the most severe clinical symptoms pointed to both inflammatory and fibrotic cell lineages. Further, we did a preliminary analysis of potential circulating biomarkers that mark the different plaques phenotypes. In conclusion, the definition of the plaque at risk for a thrombotic event can be fine-tuned by in-depth transcriptomic-based phenotyping. These differential plaque phenotypes prove clinically relevant for both carotid and coronary artery plaques and point to distinct underlying biology of symptomatic lesions.

Artemether-lumefantrine treatment failure of uncomplicated Plasmodium falciparum malaria in travellers coming from Angola and Mozambique.

The failure of artemisinin combination therapy (ACT) in malaria patients returning from endemic regions may be driven by parasite resistance to this treatment. ACT is used globally as the first-line treatment for Plasmodium falciparum malaria. However, artemisinin-resistant strains of P. falciparum have emerged and spread across Southeast Asia, with the risk of reaching high malaria burden regions in Africa and elsewhere. Here, we report on two malaria imported cases from Africa with possible parasite resistance to the ACT artemether-lumefantrine (AL). Case presentation: Two middle-aged males returning from Angola and Mozambique developed malaria symptoms in Portugal, where they were diagnosed and received treatment with AL as hospital inpatients. After apparent cure and discharge from hospital, these individuals returned to hospital showing signs of late clinical failure. Molecular analysis was performed across a number of drug resistance associated genes. No evidence of pfk13-mediated artemisinin resistance was found. Both subjects had complete parasite clearance after treatment with non-ACT antimalarials. Conclusion: Our case-studies highlights the need for close monitoring of signs of unsatisfactory antimalarial efficacy among AL treated patients and the possible implication of other genes or mutations in the parasite response to ACTs.

Using deep learning to identify recent positive selection in malaria parasite sequence data.

BACKGROUND: Malaria, caused by Plasmodium parasites, is a major global public health problem. To assist an understanding of malaria pathogenesis, including drug resistance, there is a need for the timely detection of underlying genetic mutations and their spread. With the increasing use of whole-genome sequencing (WGS) of Plasmodium DNA, the potential of deep learning models to detect loci under recent positive selection, historically signals of drug resistance, was evaluated. METHODS: A deep learning-based approach (called "DeepSweep") was developed, which can be trained on haplotypic images from genetic regions with known sweeps, to identify loci under positive selection. DeepSweep software is available from https://github.com/WDee/Deepsweep . RESULTS: Using simulated genomic data, DeepSweep could detect recent sweeps with high predictive accuracy (areas under ROC curve > 0.95). DeepSweep was applied to Plasmodium falciparum (n = 1125; genome size 23 Mbp) and Plasmodium vivax (n = 368; genome size 29 Mbp) WGS data, and the genes identified overlapped with two established extended haplotype homozygosity methods (within-population iHS, across-population Rsb) (~ 60-75% overlap of hits at P < 0.0001). DeepSweep hits included regions proximal to known drug resistance loci for both P. falciparum (e.g. pfcrt, pfdhps and pfmdr1) and P. vivax (e.g. pvmrp1). CONCLUSION: The deep learning approach can detect positive selection signatures in malaria parasite WGS data. Further, as the approach is generalizable, it may be trained to detect other types of selection. With the ability to rapidly generate WGS data at low cost, machine learning approaches (e.g. DeepSweep) have the potential to assist parasite genome-based surveillance and inform malaria control decision-making.

Flavivirus integrations in Aedes aegypti are limited and highly conserved across samples from different geographic regions unlike integrations in Aedes albopictus.

Mosquitoes of the genus Aedes are the main vectors of many viruses, e.g. dengue and Zika, which affect millions of people each year and for which there are limited treatment options. Understanding how Aedes mosquitoes tolerate high viral loads may lead to better disease control strategies. Elucidating endogenous viral elements (EVEs) within vector genomes may give exploitable biological insights. Previous studies have reported the presence of a large number of EVEs in Aedes genomes. Here we investigated if flavivirus EVEs are conserved across populations and different Aedes species by using ~ 500 whole genome sequence libraries from Aedes aegypti and Aedes albopictus, sourced from colonies and field mosquitoes across continents. We found that nearly all flavivirus EVEs in the Ae. aegypti reference genome originate from four separate putative viral integration events, and that they are highly conserved across geographically diverse samples. By contrast, flavivirus EVEs in the Ae. albopictus reference genome originate from up to nine distinct integration events and show low levels of conservation, even within samples from narrow geographical ranges. Our analysis suggests that flaviviruses integrated as long sequences and were subsequently fragmented and shuffled by transposable elements. Given that EVEs of Ae. aegypti and Ae. albopictus belong to different phylogenetic clades and have very differing levels of conservation, they may have different evolutionary origins and potentially different functional roles.

An inter-laboratory study to investigate the impact of the bioinformatics component on microbiome analysis using mock communities.

Despite the advent of whole genome metagenomics, targeted approaches (such as 16S rRNA gene amplicon sequencing) continue to be valuable for determining the microbial composition of samples. Amplicon microbiome sequencing can be performed on clinical samples from a normally sterile site to determine the aetiology of an infection (usually single pathogen identification) or samples from more complex niches such as human mucosa or environmental samples where multiple microorganisms need to be identified. The methodologies are frequently applied to determine both presence of micro-organisms and their quantity or relative abundance. There are a number of technical steps required to perform microbial community profiling, many of which may have appreciable precision and bias that impacts final results. In order for these methods to be applied with the greatest accuracy, comparative studies across different laboratories are warranted. In this study we explored the impact of the bioinformatic approaches taken in different laboratories on microbiome assessment using 16S rRNA gene amplicon sequencing results. Data were generated from two mock microbial community samples which were amplified using primer sets spanning five different variable regions of 16S rRNA genes. The PCR-sequencing analysis included three technical repeats of the process to determine the repeatability of their methods. Thirteen laboratories participated in the study, and each analysed the same FASTQ files using their choice of pipeline. This study captured the methods used and the resulting sequence annotation and relative abundance output from bioinformatic analyses. Results were compared to digital PCR assessment of the absolute abundance of each target representing each organism in the mock microbial community samples and also to analyses of shotgun metagenome sequence data. This ring trial demonstrates that the choice of bioinformatic analysis pipeline alone can result in different estimations of the composition of the microbiome when using 16S rRNA gene amplicon sequencing data. The study observed differences in terms of both presence and abundance of organisms and provides a resource for ensuring reproducible pipeline development and application. The observed differences were especially prevalent when using custom databases and applying high stringency operational taxonomic unit (OTU) cut-off limits. In order to apply sequencing approaches with greater accuracy, the impact of different analytical steps needs to be clearly delineated and solutions devised to harmonise microbiome analysis results.

Distinctive genetic structure and selection patterns in Plasmodium vivax from South Asia and East Africa.

Despite the high burden of Plasmodium vivax malaria in South Asian countries, the genetic diversity of circulating parasite populations is not well described. Determinants of antimalarial drug susceptibility for P. vivax in the region have not been characterised. Our genomic analysis of global P. vivax (n = 558) establishes South Asian isolates (n = 92) as a distinct subpopulation, which shares ancestry with some East African and South East Asian parasites. Signals of positive selection are linked to drug resistance-associated loci including pvkelch10, pvmrp1, pvdhfr and pvdhps, and two loci linked to P. vivax invasion of reticulocytes, pvrbp1a and pvrbp1b. Significant identity-by-descent was found in extended chromosome regions common to P. vivax from India and Ethiopia, including the pvdbp gene associated with Duffy blood group binding. Our investigation provides new understanding of global P. vivax population structure and genomic diversity, and genetic evidence of recent directional selection in this important human pathogen.

Failure of rapid diagnostic tests in Plasmodium falciparum malaria cases among travelers to the UK and Ireland: Identification and characterisation of the parasites.

OBJECTIVES: Our objective was to systematically investigate false-negative histidine-rich protein 2 rapid diagnostic tests (HRP2-RDT) in imported Plasmodium falciparum malaria cases from travelers to the UK and the Republic of Ireland (RoI). METHODS: Five imported malaria cases in travellers returning to the UK and RoI from East Africa were reported to the PHE Malaria Reference Laboratory as negative according to histidine-rich protein (HRP2)-RDT. The cases were systematically investigated using microscopic, RDT, molecular, genomic, and in in vitro approaches. RESULTS: In each case, HRP2-RDT was negative, whereas microscopy confirmed the presence of P. falciparum. Further analysis revealed that the genes encoding HRP2 and HRP3 were deleted in three of the five cases. Whole-genome sequencing in one of these isolates confirmed deletions in P. falciparum chromosomes 8 and 13. Our study produced evidence that the fourth case, which had high parasitemia at clinical presentation, was a rare example of antigen saturation ('prozone-like effect'), leading to a false negative in the HRP2-RDT, while the fifth case was due to low parasitemia. CONCLUSIONS: False-negative HRP2-RDT results with P. falciparum are concerning. Our findings emphasise the necessity of supporting the interpretation of RDT results with microscopy, in conjunction with clinical observations, and sets out a systematic approach to identifying parasites carrying pfhrp2 and pfhrp3 deletions.

Drug resistance profile and clonality of Plasmodium falciparum parasites in Cape Verde: the 2017 malaria outbreak.

BACKGROUND: Cape Verde is an archipelago located off the West African coast and is in a pre-elimination phase of malaria control. Since 2010, fewer than 20 Plasmodium falciparum malaria cases have been reported annually, except in 2017, when an outbreak in Praia before the rainy season led to 423 autochthonous cases. It is important to understand the genetic diversity of circulating P. falciparum to inform on drug resistance, potential transmission networks and sources of infection, including parasite importation. METHODS: Enrolled subjects involved malaria patients admitted to Dr Agostinho Neto Hospital at Praia city, Santiago island, Cape Verde, between July and October 2017. Neighbours and family members of enrolled cases were assessed for the presence of anti-P. falciparum antibodies. Sanger sequencing and real-time PCR was used to identify SNPs in genes associated with drug resistance (e.g., pfdhfr, pfdhps, pfmdr1, pfk13, pfcrt), and whole genome sequencing data were generated to investigate the population structure of P. falciparum parasites. RESULTS: The study analysed 190 parasite samples, 187 indigenous and 3 from imported infections. Malaria cases were distributed throughout Praia city. There were no cases of severe malaria and all patients had an adequate clinical and parasitological response after treatment. Anti-P. falciparum antibodies were not detected in the 137 neighbours and family members tested. No mutations were detected in pfdhps. The triple mutation S108N/N51I/C59R in pfdhfr and the chloroquine-resistant CVIET haplotype in the pfcrt gene were detected in almost all samples. Variations in pfk13 were identified in only one sample (R645T, E668K). The haplotype NFD for pfmdr1 was detected in the majority of samples (89.7%). CONCLUSIONS: Polymorphisms in pfk13 associated with artemisinin-based combination therapy (ACT) tolerance in Southeast Asia were not detected, but the majority of the tested samples carried the pfmdr1 haplotype NFD and anti-malarial-associated mutations in the the pfcrt and pfdhfr genes. The first whole genome sequencing (WGS) was performed for Cape Verdean parasites that showed that the samples cluster together, have a very high level of similarity and are close to other parasites populations from West Africa.