Plasmodium vivax genomic surveillance in the Peruvian Amazon with Pv AmpliSeq assay.

BACKGROUND: Plasmodium vivax is the most predominant malaria species in Latin America, constituting 71.5% of malaria cases in 2021. With several countries aiming for malaria elimination, it is crucial to prioritize effectiveness of national control programs by optimizing the utilization of available resources and strategically implementing necessary changes. To support this, there is a need for innovative approaches such as genomic surveillance tools that can investigate changes in transmission intensity, imported cases and sources of reintroduction, and can detect molecular markers associated with drug resistance. METHODOLOGY/PRINCIPAL FINDINGS: Here, we apply a modified highly-multiplexed deep sequencing assay: Pv AmpliSeq v2 Peru. The tool targets a newly developed 41-SNP Peru barcode for parasite population analysis within Peru, the 33-SNP vivaxGEN-geo panel for country-level classification, and 11 putative drug resistance genes. It was applied to 230 samples from the Peruvian Amazon (2007-2020), generating baseline surveillance data. We observed a heterogenous P. vivax population with high diversity and gene flow in peri-urban areas of Maynas province (Loreto region) with a temporal drift using all SNPs detected by the assay (nSNP = 2909). In comparison, in an indigenous isolated area, the parasite population was genetically differentiated (FST = 0.07-0.09) with moderate diversity and high relatedness between isolates in the community. In a remote border community, a clonal P. vivax cluster was identified, with distinct haplotypes in drug resistant genes and ama1, more similar to Brazilian isolates, likely representing an introduction of P. vivax from Brazil at that time. To test its applicability for Latin America, we evaluated the SNP Peru barcode in P. vivax genomes from the region and demonstrated the capacity to capture local population clustering at within-country level. CONCLUSIONS/SIGNIFICANCE: Together this data shows that P. vivax transmission is heterogeneous in different settings within the Peruvian Amazon. Genetic analysis is a key component for regional malaria control, offering valuable insights that should be incorporated into routine surveillance.

Genomic surveillance of malaria parasites in an indigenous community in the Peruvian Amazon.

Hard-to-reach communities represent Peru's main challenge for malaria elimination, but information about transmission in these areas is scarce. Here, we assessed Plasmodium vivax (Pv) and P. falciparum (Pf) transmission dynamics, resistance markers, and Pf hrp2/3 deletions in Nueva Jerusalén (NJ), a remote, indigenous community in the Peruvian Amazon with high population mobility. We collected samples from November 2019 to May 2020 by active (ACD) and passive case detection (PCD) in NJ. Parasites were identified with microscopy and PCR. Then, we analyzed a representative set of positive-PCR samples (Pv = 68, Pf = 58) using highly-multiplexed deep sequencing assays (AmpliSeq) and compared NJ parasites with ones from other remote Peruvian areas using population genetics indexes. The ACD intervention did not reduce malaria cases in the short term, and persistent malaria transmission was observed (at least one Pv infection was detected in 96% of the study days). In Nueva Jerusalen, the Pv population had modest genetic diversity (He = 0.27). Pf population had lower diversity (He = 0.08) and presented temporal clustering, one of these clusters linked to an outbreak in February 2020. Moreover, Pv and Pf parasites from NJ exhibited variable levels of differentiation (Pv Fst = 0.07-0.52 and Pf Fst = 0.11-0.58) with parasites from other remote areas. No artemisin resistance mutations but chloroquine (57%) and sulfadoxine-pyrimethamine (35-67%) were detected in NJ's Pf parasites. Moreover, pfhrp2/3 gene deletions were common (32-50% of parasites with one or both genes deleted). The persistent Pv transmission and the detection of a Pf outbreak with parasites genetically distinct from the local ones highlight the need for tailored interventions focusing on mobility patterns and imported infections in remote areas to eliminate malaria in the Peruvian Amazon.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana.

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. In low transmission settings, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. We performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n = 1,409) through estimation of identity-by-descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 multi-isolate clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally observed 61 nonsynonymous substitutions that increased markedly in frequency over the study period as well as a novel pfk13 mutation (G718S). However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions, given that clones carrying drug resistance polymorphisms do not demonstrate enhanced persistence or higher abundance than clones carrying polymorphisms of comparable frequency that are unrelated to resistance. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Molecular Surveillance of Artemisinin-Resistant Plasmodium falciparum Parasites in Mining Areas of the Roraima Indigenous Territory in Brazil.

Multidrug- and artemisinin-resistant (ART-R) Plasmodium falciparum (Pf) parasites represent a challenge for malaria elimination worldwide. Molecular monitoring in the Kelch domain region (pfk13) gene allows tracking mutations in parasite resistance to artemisinin. The increase in illegal miners in the Roraima Yanomami indigenous land (YIL) could favor ART-R parasites. Thus, this study aimed to investigate ART-R in patients from illegal gold mining areas in the YIL of Roraima, Brazil. A questionnaire was conducted, and blood was collected from 48 patients diagnosed with P. falciparum or mixed malaria (Pf + P. vivax). The DNA was extracted and the pfk13 gene was amplified by PCR. The amplicons were subjected to DNA-Sanger-sequencing and the entire amplified fragment was analyzed. Among the patients, 96% (46) were from illegal mining areas of the YIL. All parasite samples carried the wild-type genotypes/ART-sensitive phenotypes. These data reinforce the continued use of artemisinin-based combination therapies (ACTs) in Roraima, as well as the maintenance of systematic monitoring for early detection of parasite populations resistant to ART, mainly in regions with an intense flow of individuals from mining areas, such as the YIL. This is especially true when the achievement of falciparum malaria elimination in Brazil is planned and expected by 2030.

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.

Molecular characterization of the Plasmodium falciparum k13 gene helix domain in samples from native communities of Condorcanqui, Amazonas, Perú / Caracterización molecular del dominio de la hélice del gen k13 de Plasmodium falciparum en muestras de comunidades nativas de Condorcanqui, Amazonas, Perú.

Introduction. Resistance of Plasmodium falciparum to different antimalarial drugs is an obstacle to disease elimination. The artemisinin-resistant genotype of P. falciparum can be assessed by examining polymorphisms in the helix domain of the Pfk13 gene. The World Health Organization recommends these mutations as molecular markers to detect artemisinin-resistant in countries where P. falciparum malaria is endemic. Objective. To identify artemisinin resistance-related mutations present in the helix domain of the P. falciparum k13 gene. Materials and methods. We collected a total of 51 samples through passive case detection, positive for Plasmodium by microscopy, from six communities in the district of Río Santiago in Condorcanqui, Amazonas. Molecular species confirmation was performed by real-time PCR and Pfk13 helix domain was amplified and sequenced by capillary electrophoresis. The obtained sequences were compared with the wild type 3D7 reference strain. Results. A total of 51 positive samples were confirmed for P. falciparum from the communities of Ayambis, Chapiza, Palometa, Muchinguis, Alianza Progreso and Caterpiza. DNA sequences alignment showed the absence of resistance-associated mutations in the k13 gene of the collected samples. Discussion. The obtained results are consistent with similar studies conducted in other South American countries, including Perú, so these data provide a baseline for artemisinin-resistance molecular surveillance in the Amazon region and reinforce the efficacy of artemisinin-based combination therapy in this area.

Introducción. La resistencia de Plasmodium falciparum a diferentes fármacos antipalúdicos es un obstáculo para eliminar la enfermedad. El genotipo resistente de P. falciparum a la artemisinina puede evaluarse examinando los polimorfismos en el dominio de la hélice del gen Pfk13. La Organización Mundial de la Salud recomienda utilizar estas mutaciones como marcadores moleculares para detectar la resistencia a la artemisinina en países donde la malaria por P. falciparum es endémica. Objetivo. Identificar mutaciones relacionadas con la resistencia a artemisinina presentes en el dominio de la hélice del gen k13 de P. falciparum. Materiales y métodos. Mediante la detección pasiva de casos, se recolectaron 51 muestras positivas por microscopía para Plasmodium, provenientes de seis comunidades del distrito de Río Santiago en Condorcanqui, Amazonas. Se realizó la confirmación molecular de la especie mediante PCR en tiempo real y el dominio de la hélice del gen Pfk13 se amplificó y secuenció por electroforesis capilar. Las secuencias obtenidas se compararon con la cepa de referencia 3D7 de fenotipo silvestre. Resultados. Se confirmó un total de 51 muestras positivas para P. falciparum, provenientes de las comunidades de Ayambis, Chapiza, Palometa, Muchinguis, Alianza Progreso y Caterpiza. Después del alineamiento de las secuencias de ADN, se determinó que las muestras no presentaron mutaciones asociadas con resistencia en el gen K13. Discusión. Los resultados obtenidos son coherentes con estudios similares realizados en otros países de Sudamérica, incluyendo Perú. Estos datos proporcionan una línea base para la vigilancia molecular de resistencia a artemisinina en la región Amazonas y refuerzan la eficacia de la terapia combinada con artemisinina en esta área.

Targeted and whole-genome sequencing reveal a north-south divide in P. falciparum drug resistance markers and genetic structure in Mozambique

Mozambique is one of the four African countries which account for over half of all malaria deaths worldwide, yet little is known about the parasite genetic structure in that country. We performed P. falciparum amplicon and whole genome sequencing on 2251 malaria-infected blood samples collected in 2015 and 2018 in seven provinces of Mozambique to genotype antimalarial resistance markers and interrogate parasite population structure using genome-wide microhaplotyes. Here we show that the only resistance-associated markers observed at frequencies above 5% were pfmdr1-184F (59%), pfdhfr-51I/59 R/108 N (99%) and pfdhps-437G/540E (89%). The frequency of pfdhfr/pfdhps quintuple mutants associated with sulfadoxine-pyrimethamine resistance increased from 80% in 2015 to 89% in 2018 (p < 0.001), with a lower expected heterozygosity and higher relatedness of microhaplotypes surrounding pfdhps mutants than wild-type parasites suggestive of recent selection. pfdhfr/pfdhps quintuple mutants also increased from 72% in the north to 95% in the south (2018; p < 0.001). This resistance gradient was accompanied by a concentration of mutations at pfdhps-436 (17%) in the north, a south-to-north increase in the genetic complexity of P. falciparum infections (p = 0.001) and a microhaplotype signature of regional differentiation. The parasite population structure identified here offers insights to guide antimalarial interventions and epidemiological surveys.

Fixation of pfcrt chloroquine resistance alleles in Plasmodium falciparum clinical isolates collected from unrest tribal agencies of Pakistan / Fixação de alelos de resistência à cloroquina pfcrt em isolados clínicos de Plasmodium falciparum coletados de agitações de agências tribais do Paquistão

Abstract Plasmodium falciparum resistance to Chloroquine (CQ) is a significant cause of mortality and morbidity worldwide. There is a paucity of documented data on the prevalence of CQ-resistant mutant haplotypes of Pfcrt and Pfmdr1 genes from malaria-endemic war effected Federally Administered Tribal Areas of Pakistan. The objective of this study was to investigate the prevalence of P. falciparum CQ-resistance in this area. Clinical isolates were collected between May 2017 and May 2018 from North Waziristan and South Waziristan agencies of Federally Administrated Trial Area. Subsequently, Giemsa-stained blood smears were examined to detect Plasmodium falciparum. Extraction of malarial DNA was done from microscopy positive P. falciparum samples, and P. falciparum infections were confirmed by nested PCR (targeting Plasmodium small subunit ribosomal ribonucleic acid (ssrRNA) genes). All PCR confirmed P. falciparum samples were sequenced by pyrosequencing to find out mutation in Pfcrt gene at codon K76T and in pfmdr1 at codons N86Y, Y184F, N1042D, and D1246Y. Out of 121 microscopies positive P. falciparum cases, 109 samples were positive for P. falciparum by nested PCR. Pfcrt K76T mutation was found in 96% of isolates, Pfmdr1 N86Y mutation was observed in 20%, and 11% harboured Y184F mutation. All samples were wild type for Pfmdr1 codon N1042D and D1246Y. In the FATA, Pakistan, the frequency of resistant allele 76T remained high despite the removal of CQ. However, current findings of the study suggest complete fixation of P. falciparum CQ-resistant genotype in the study area.

Resumo A resistência do Plasmodium falciparum à cloroquina (CQ) é uma causa significativa de mortalidade e morbidade em todo o mundo. Há uma escassez de dados documentados sobre a prevalência de haplótipos mutantes CQ-resistentes dos genes Pfcrt e Pfmdr1 da guerra endêmica da malária em áreas tribais administradas pelo governo federal do Paquistão. O objetivo deste estudo foi investigar a prevalência de resistência a CQ de P. falciparum nesta área. Isolados clínicos foram coletados entre maio de 2017 e maio de 2018 nas agências do Waziristão do Norte e do Waziristão do Sul da Área de Ensaio Administrada Federalmente. Posteriormente, esfregaços de sangue corados com Giemsa foram examinados para detectar Plasmodium falciparum. A extração do DNA da malária foi feita a partir de amostras de P. falciparum positivas para microscopia, e as infecções por P. falciparum foram confirmadas por nested PCR (visando genes de ácido ribonucleico ribossômico de subunidade pequena de Plasmodium (ssrRNA)). Todas as amostras de P. falciparum confirmadas por PCR foram sequenciadas por pirosequenciamento para descobrir a mutação no gene Pfcrt no códon K76T e em pfmdr1 nos códons N86Y, Y184F, N1042D e D1246Y. De 121 microscopias de casos positivos de P. falciparum, 109 amostras foram positivas para P. falciparum por nested PCR. A mutação Pfcrt K76T foi encontrada em 96% dos isolados, a mutação Pfmdr1 N86Y foi observada em 20% e 11% abrigou a mutação Y184F. Todas as amostras eram do tipo selvagem para o códon N1042D e D1246Y de Pfmdr1. No FATA, Paquistão, a frequência do alelo resistente 76T permaneceu alta apesar da remoção de CQ. No entanto, as descobertas atuais do estudo sugerem a fixação completa do genótipo resistente a CQ de P. falciparum na área de estudo.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis.

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

Polymorphism in the pvcrt-o and pvmdr-1 genes of Plasmodium vivax associated with a better prognosis for malaria.

The early molecular identification of strains of Plasmodium vivax that have a worse prognosis is important to stratify the risk of complications and choice of conduct made by medical teams. Thus, the aim of the present study was to associate the presence of polymorphisms in the pvmdr-1 and pvcrt-o resistance genes of P. vivax in patients with better or worse prognosis. This cross-sectional epidemiological study was conducted based on data obtained from the records of 120 patients diagnosed with malaria in the Brazilian Amazon. The T958M and F1076L mutations of the pvmdr-1 gene had a frequency of 3.3 and 4.2%, respectively, and primo-infected patients had a 17 times greater chance of being infected with protozoa with the T958M mutation compared to patients with previous episodes. Regarding pvcrt-o, the C393T and T786C polymorphisms had a frequency of 14.2 and 3.3%, respectively, and self-declared white patients had a 3.1 times greater chance of being infected with protozoa with the C393T polymorphism. In addition, patients with this pvcrt-o polymorphism had lower concentrations of C-reactive protein, indicating a better prognosis. These data present clues of genetic indicators useful for assessing the virulence of the parasite and the prognosis of patients with vivax malaria.

Drug resistance and population structure of Plasmodium falciparum and Plasmodium vivax in the Peruvian Amazon.

Malaria is a major health problem in Peru despite substantial progress achieved by the ongoing malaria elimination program. This study explored the population genetics of 63 Plasmodium falciparum and 170 P. vivax cases collected in the Peruvian Amazon Basin between 2015 and 2019. Microscopy and PCR were used for malaria detection and positive samples were genotyped at neutral and drug resistance-associated regions. The P. falciparum population exhibited a low nucleotide diversity (π = 0.02) whereas the P. vivax population presented a higher genetic diversity (π = 0.34). All P. falciparum samples (n = 63) carried chloroquine (CQ) resistant mutations on Pfcrt. Most P. falciparum samples (53 out of 54) carried sulfadoxine (SD) resistant mutations on Pfdhfr and Pfdhps. No evidence was found of artemisinin resistance mutations on kelch13. Population structure showed that a single cluster accounted for 93.4% of the P. falciparum samples whereas three clusters were found for P. vivax. Our study shows a low genetic diversity for both species with significant differences in genetic sub-structuring. The high prevalence of CQ-resistance mutations could be a result of indirect selection pressures driven by the P. vivax treatment scheme. These results could be useful for public health authorities to safeguard the progress that Peru has achieved towards malaria elimination.

Polymorphisms in exon 11 of the mptl-1 gene and monepantel resistance in Haemonchus contortus.

Chemical treatments are the main strategy to control gastrointestinal nematodes in sheep, and the emergence of anthelmintic resistance, as consequence, results in control failures and leads to economic losses. Thus, molecular tests may constitute an excellent tool for the early detection of anthelmintic resistance-related mutations. Thus, a polymerase chain reaction (PCR)-based genotyping assay followed by polyacrylamide gel electrophoresis (PAGE) was developed to detect polymorphisms in exon 11 of the acetylcholine receptor monepantel-1 gene (mptl-1) that were previously associated with monepantel resistance through a genome-wide study in Haemonchus contortus. DNA samples recovered from individual and pooled third-stage larvae from two susceptible field-derived isolates and five (three in vivo-derived and two field-derived) resistant populations were used. New polymorphisms, including a 6-bp deletion and a 3-bp insertion, were detected in resistant individuals. These indels, confirmed using sequencing of cloned PCR products, are predicted to result in amino acid changes in transmembrane domain 2 (TMD2) of the MPTL-1 protein. The two susceptible isolates showed only the presence of the wild-type allele (100%), whereas lower frequencies of the wild-type allele were detected in monepantel-resistant populations (11.1 to 66.7%). These findings report new polymorphisms in the mptl-1 gene, validate the results obtained through genomic mapping for monepantel resistance, and provide a PCR-based assay to genotype indels located in exon 11 of mptl-1 in H. contortus.

Molecular surveillance for anti-malarial drug resistance and genetic diversity of Plasmodium falciparum after chloroquine and sulfadoxine-pyrimethamine withdrawal in Quibdo, Colombia, 2018.

BACKGROUND: Resistance to anti-malarial drugs is associated with polymorphisms in target genes and surveillance for these molecular markers is important to detect the emergence of mutations associated with drug resistance and signal recovering sensitivity to anti-malarials previously used. METHODS: The presence of polymorphisms in genes associated with Plasmodium falciparum resistance to chloroquine and sulfadoxine-pyrimethamine was evaluated by Sanger sequencing, in 85 P. falciparum day of enrollment samples from a therapeutic efficacy study of artemether-lumefantrine conducted in 2018-2019 in Quibdo, Colombia. Samples were genotyped to assess mutations in pfcrt (codons 72-76), pfdhfr (codons 51, 59, 108, and 164), and pfdhps genes (codons 436, 437, 540, and 581). Further, the genetic diversity of infections using seven neutral microsatellites (NMSs) (C2M34, C3M69, Poly α, TA1, TA109, 2490, and PfPK2) was assessed. RESULTS: All isolates carried mutant alleles for pfcrt (K76T and N75E), and for pfdhfr (N51I and S108N), while for pfdhps, mutations were observed only for codon A437G (32/73, 43.8%). Fifty samples (58.8%) showed a complete neutral microsatellites (NMS) profile. The low mean number of alleles (2 ± 0.57) per locus and mean expected heterozygosity (0.17 ± 0.03) showed a reduced genetic diversity. NMS multilocus genotypes (MMG) were built and nine MMG were identified. CONCLUSIONS: Overall, these findings confirm the fixation of chloroquine and pyrimethamine-resistant alleles already described in the literature, implying that these drugs are not currently appropriate for use in Colombia. In contrast, mutations in the pfdhps gene were only observed at codon 437, an indication that full resistance to sulfadoxine has not been achieved in Choco. MMGs found matched the clonal lineage E variant 1 previously reported in northwestern Colombia.

Population replacement of benzimidazole-resistant Haemonchus contortus with susceptible strains: evidence of changes in the resistance status.

The spread of anthelmintic resistance (AR) in nematode populations threatens the viability of sheep production systems worldwide, and warrants the adoption of sensitive, practical, and standardized tests to detect AR. The aim of this study was to characterize the replacement of an Haemonchus contortus population resistant to benzimidazoles (BZDs) by a susceptible one, by means of both phenotypic and genotypic techniques. Phenotypic methods to assess BZD resistance included in vivo tests, such as the fecal egg count reduction test (FECRT), and in vitro tests, such as the egg hatch assay (EHA). Additionally, genotypification of polymorphisms associated with BZD resistance by sequencing a fragment of the isotype 1 ß-tubulin gene was carried out. The initial, BZD-resistant population (initial Balcarce population) exhibited an egg count reduction (ECR) of 59.3%. Following refugium replacement, the final population (final Balcarce population) exhibited an ECR of 95.2%. For the initial Balcarce population, the median effective dose (ED50) for the EHA was 0.607 µg thiabendazole (TBZ)/mL, with a rate of eclosion at a discriminating dose (EDD) of 0.1 µg TBZ/mL of 76.73%. For the final Balcarce population, ED50 was 0.02 µg TBZ/mL, and EDD was 1.97%. In the initial population, 93% of the analyzed individuals exhibited genotypic combinations associated with BZD resistance (53% Phe/Phe167-Tyr/Tyr200, 37% Phe/Tyr167-Phe/Tyr200, and 3% Phe/Tyr167-Glu/Leu198). Conversely, no combination associated with resistance was found in individuals from the final population. All of the tests were useful for detecting AR to BZDs. The results from the genetic and phenotypical studies were consistent, and the resulting information greatly aided in interpreting the outcomes of the population replacement and the potential impact of this strategy on management of AR.

Moving beyond the state of the art of understanding resistance mechanisms in hookworms: confirming old and suggesting new associated SNPs.

Hookworms represent a serious problem for human and animal health in different parts of the world. One of the suggested control strategies for parasitosis caused by members of the Ancylostomatidae family is mass drug aministration with benzimidazole compounds. This strategy has been proven to lead to the establishment of resistant strains in several nematodes related to SNPs at codons 167, 198 and 200 of the beta-tubulin isotype-1 gene. Through bioassay and in vivo analysis, we successfully isolated an albendazole-resistant A. ceylanicum strain by drug selective pressure. We observed a strong correlation between the presence of SNPs at codon 198 and drug resistance. We also described for the first time, in hookworms, the presence of SNP A200L, already described at low frequencies in ruminant nematodes. The results presented here are important for updating the current knowledge about anthelmintic resistance in hookworms. The answers and the new questions raised may provide a basis for the establishment of more effective control strategies.

Phytochemical modulation of P-Glycoprotein and its gene expression in an ivermectin-resistant Haemonchus contortus isolate in vitro.

Haemonchus contortus is the most important gastrointestinal nematode in small ruminant systems worldwide and has developed resistance to several drugs, including ivermectin (IVM). IVM is not only a veterinary drug but also a safe, broad-spectrum, antiparasitic drug used in humans. One of the main IVM-resistance mechanisms in H. contortus involves P-glycoprotein (PgP), a trans-membrane transport protein that rids worm cells from toxic molecules. This study aimed to evaluate the anthelmintic activity of IVM, alone or combined with main terpenes of essential oils (alpha-terpinene, beta-citronellol, beta-pinene, citronellal, limonene, menthol, and terpinolene) and with phenolic compounds (epicatechin, epigallocatechin, gallocatechin, pentagalloylglucose, procyanidin, and quercetin). All compounds were tested, alone or combined with IVM, against susceptible (HcS) and resistant (HcR) isolates of H. contortus through the larval development test (LDT) and the adult motility assay (AMT) using verapamil (VP), a known PgP modulator, as a control. Results for the LDT determined that the lethal concentration required to kill 50% of nematodes (LC50) with IVM was 10 times greater (0.01 µg/mL) for HcR than for HcS (0.001 µg/mL). The combination IVM + VP inhibited the activity of PgP in HcR resulting in a LC50 = 0.002 ug.mL-1. Although limonene was the least effective and alpha-terpinene the most effective terpene when tested alone against HcR, the best combinations were IVM + limonene and IVM + quercetin both produced LC50 = 0.002 µg/mL (similar to IVM+VP) which were chosen for subsequent tests. Because adult parasites are the final target for anthelmintics, IVM was evaluated in HcS (LC50 = 0.067 µg/mL) and HcR (LC50 =164.94 µg/mL) through the AMT. Results obtained with IVM + VP (LC50 = 0.020 µg/mL) in HcR were similar to IVM + limonene (LC50 = 0.028 µg/mL) and outperformed IVM + quercetin (LC50 = 1.39 µg/mL). RNA extracts from HcR adult worms exposed to IVM, IVM+VP, and IVM + limonene were evaluated for PgP expression by RT-PCR. For most concentrations, PgP-9 was significantly more expressed in worms treated with IVM alone than in worms treated with IVM + VP or IVM + limonene. Our results suggest that limonene is involved in the modulation of the PgP-9 gene and that it can restore the activity of IVM in the HcR isolate down to levels seen in HcS. Limonene is one of the main compounds found in citrus peel and has the potential to be both safe and affordable if used in combination with IVM to restore its anthelmintic effects against multi-drug-resistant H. contortus isolates. Our results also suggest that we may be more successful by combining natural products with failing commercial anthelmintics than trying to find natural substitutes for them.

Evaluation of a parasite-density based pooled targeted amplicon deep sequencing (TADS) method for molecular surveillance of Plasmodium falciparum drug resistance genes in Haiti.

Sequencing large numbers of individual samples is often needed for countrywide antimalarial drug resistance surveillance. Pooling DNA from several individual samples is an alternative cost and time saving approach for providing allele frequency (AF) estimates at a population level. Using 100 individual patient DNA samples of dried blood spots from a 2017 nationwide drug resistance surveillance study in Haiti, we compared codon coverage of drug resistance-conferring mutations in four Plasmodium falciparum genes (crt, dhps, dhfr, and mdr1), for the same deep sequenced samples run individually and pooled. Samples with similar real-time PCR cycle threshold (Ct) values (+/- 1.0 Ct value) were combined with ten samples per pool. The sequencing success for samples in pools were higher at a lower parasite density than the individual samples sequence method. The median codon coverage for drug resistance-associated mutations in all four genes were greater than 3-fold higher in the pooled samples than in individual samples. The overall codon coverage distribution for pooled samples was wider than the individual samples. The sample pools with < 40 parasites/µL blood showed more discordance in AF calls for dhfr and mdr1 between the individual and pooled samples. This discordance in AF estimation may be due to low amounts of parasite DNA, which could lead to variable PCR amplification efficiencies. Grouping samples with an estimated ≥ 40 parasites/µL blood prior to pooling and deep sequencing yielded the expected population level AF. Pooling DNA samples based on estimates of > 40 parasites/µL prior to deep sequencing can be used for rapid genotyping of a large number of samples for these four genes and possibly other drug resistant markers in population-based studies. As Haiti is a low malaria transmission country with very few mixed infections and continued chloroquine sensitivity, the pooled sequencing approach can be used for routine national molecular surveillance of resistant parasites.

ABC-transporter gene expression in ivermectin-susceptible and resistant Haemonchus contortus isolates.

Mammalian efflux transporters of the ATP-binding cassette (ABC) regulate cellular levels of endo- and xenobiotics by transporting molecules across cell membranes and are involved in diverse biological processes. Over-expression of these ABC transporters has been involved in macrocyclic lactone resistance. The main goal of this work was to compare the gene expression of the whole ABC-transporter superfamily in isolates of the sheep nematode Haemonchus contortus with different degrees of susceptibility to ivermectin (IVM). Additionally, the effects of in vivo IVM treatment were evaluated in the resistant isolates. Parasite-free Corriedale lambs were artificially infected with either IVM-susceptible or IVM-resistant H. contortus isolates. The differential expression of ABC transcripts in H. contortus female worms with differential susceptibility to IVM were assessed by RNA-seq. Additionally, the transcription levels of ABC-transporter genes in IVM-resistant adult worms recovered from treated sheep at 12 and 24 h after IVM administration were compared to those of IVM-R worms collected from untreated sheep. The comparative analysis of the ABC-transcripts revealed some minor differences in the expression levels of HCON_00042800 (pgp-3), HCON_00020200.mod (ced-7c), HCON_00085890 (abt-4), HCON_00063000 (pmp-5) and HCON_00116670 (wht-8), indicating that, at transcriptional level, these ABC-genes alone cannot explain resistance in H. contortus. HCON_00130060 (pgp-9.2) was highly differentially expressed in resistant isolates compared to susceptible ones, which agrees with previous reports suggesting that pgp-9 may be one of the most relevant candidates contributing to the multi-genic nature of the IVM resistance trait.

Assessment of Plasmodium falciparum anti-malarial drug resistance markers in pfcrt and pfmdr1 genes in isolates from Honduras and Nicaragua, 2018-2021.

BACKGROUND: Central America and the island of Hispaniola have set out to eliminate malaria by 2030. However, since 2014 a notable upturn in the number of cases has been reported in the Mosquitia region shared by Nicaragua and Honduras. In addition, the proportion of Plasmodium falciparum malaria cases has increased significantly relative to vivax malaria. Chloroquine continues to be the first-line drug to treat uncomplicated malaria in the region. The objective of this study was to evaluate the emergence of chloroquine resistant strains of P. falciparum using a genetic approach. Plasmodium vivax populations are not analysed in this study. METHODS: 205 blood samples from patients infected with P. falciparum between 2018 and 2021 were analysed. The pfcrt gene fragment encompassing codons 72-76 was analysed. Likewise, three fragments of the pfmdr1 gene were analysed in 51 samples by nested PCR and sequencing. RESULTS: All samples revealed the CVMNK wild phenotype for the pfcrt gene and the N86, Y184F, S1034C, N1042D, D1246 phenotype for the pfmdr1 gene. CONCLUSIONS: The increase in falciparum malaria cases in Nicaragua and Honduras cannot be attributed to the emergence of chloroquine-resistant mutants. Other possibilities should be investigated further. This is the first study to report the genotype of pfmdr1 for five loci of interest in Central America.

Population genomics and evidence of clonal replacement of Plasmodium falciparum in the Peruvian Amazon.

Previous studies have shown that P. falciparum parasites in South America have undergone population bottlenecks resulting in clonal lineages that are differentially distributed and that have been responsible for several outbreaks different endemic regions. In this study, we explored the genomic profile of 18 P. falciparum samples collected in the Peruvian Amazon Basin (Loreto) and 6 from the Peruvian North Coast (Tumbes). Our results showed the presence of three subpopulations that matched previously typed lineages in Peru: Bv1 (n = 17), Clonet D (n = 4) and Acre-Loreto type (n = 3). Gene coverage analysis showed that none of the Bv1 samples presented coverage for pfhrp2 and pfhrp3. Genotyping of drug resistance markers showed a high prevalence of Chloroquine resistance mutations S1034C/N1042D/D1246Y in pfmdr1 (62.5%) and K45T in pfcrt (87.5%). Mutations associated with sulfadoxine and pyrimethamine treatment failure were found on 88.8% of the Bv1 samples which were triple mutants for pfdhfr (50R/51I/108N) and pfdhps (437G/540E/581G). Analysis of the pfS47 gene that allows P. falciparum to evade mosquito immune responses showed that the Bv1 lineage presented one pfS47 haplotype exclusive to Loreto and another haplotype that was present in both Loreto and Tumbes. Furthermore, a possible expansion of Bv1 was detected since 2011 in Loreto. This replacement could be a result of the high prevalence of CQ resistance polymorphisms in Bv1, which could have provided a selective advantage to the indirect selection pressures driven by the use of CQ for P. vivax treatment.