Deep Amplicon Sequencing Research Articles

Proton pump inhibitor treatment aggravates bacterial translocation in patients with advanced cirrhosis and portal hypertension.

Long-term prescription of proton pump inhibitors (PPIs) in patients with cirrhosis is common practice. However, in recent years, several observational studies have reported increased complications and negative prognostic effects of PPI treatment in these patients. Judging the significance of these associations is complicated by the fact that a plausible underlying pathomechanism has not been identified so far. In the present study, we address this important issue by investigating the impact of PPI treatment on subclinical bacterial translocation from the gut into the blood stream in patients with advanced cirrhosis and portal hypertension. Indeed, we report significantly aggravated bacterial translocation in cirrhosis patients receiving PPI treatment. This finding is highly relevant, as bacterial translocation is known to promote the development of complications and impair prognosis in patients with cirrhosis. Hence, the present study could establish a plausible link between PPI treatment and adverse effects in cirrhosis.

Efficient detection of somatic UBA1 variants and clinical scoring system predicting patients with variants in VEXAS syndrome.

To efficiently detect somatic UBA1 variants and establish a clinical scoring system predicting patients with pathogenic variants in VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome. Eighty-nine Japanese patients with clinically suspected VEXAS syndrome were recruited [81 males and 8 females; median age of onset 69.3 years (interquartile range 62.1-77.6)]. Peptide nucleic acid-clamping PCR (PNA-PCR), regular PCR targeting exon 3 clustering UBA1 variants and subsequent Sanger sequencing were conducted for variant screening. Partitioning digital PCR or targeted amplicon deep sequencing was also performed to evaluate the variant allele frequency (VAF). We developed our clinical scoring system to predict UBA1 variant-positive and -negative patients and assessed the diagnostic value of our system using receiver operating characteristics (ROC) curve analysis. Forty patients (44.9%) with reported pathogenic UBA1 variants were identified, including a case having a variant with VAF of 1.7%, using a highly sensitive method. Our clinical scoring system considering age >50 years, cutaneous lesions, lung involvement, chondritis and macrocytic anaemia efficiently predicted patients with UBA1 variants (the area under the curve for the scoring total was 0.908). Genetic screening with the combination of regular PCR and PNA-PCR detected somatic UBA1 variants with high sensitivity and specificity. Our scoring system could efficiently predict patients with UBA1 variants.

Molecular Challenges in the Diagnosis of X-Linked Chronic Granulomatous Disease: CNVs, Intronic Variants, Skewed X-Chromosome Inactivation, and Gonosomal Mosaicism.

Chronic granulomatous disease (CGD) is a prototypical inborn error of immunity affecting phagocytes, in which these cells are unable to produce reactive oxygen species. CGD is caused by defects in genes encoding subunits of the NADPH oxidase enzyme complex (CYBA, CYBB, CYBC1, NCF1, NCF2, NCF4); inflammatory responses are dysregulated, and patients are highly susceptible to recurrent severe bacterial and fungal infections. X-linked CGD (XL-CGD), caused by mutations in the CYBB gene, is the most common and severe form of CGD. In this study, we describe the analytical processes undertaken in 3 families affected with XL-CGD to illustrate several molecular challenges in the genetic diagnosis of this condition: in family 1, a girl with a heterozygous deletion of CYBB exon 13 and skewed X-chromosome inactivation (XCI); in family 2, a boy with a hemizygous deletion of CYBB exon 7, defining its consequences at the mRNA level; and in family 3, 2 boys with the same novel intronic variant in CYBB (c.1151 + 6T > A). The variant affected the splicing process, although a small fraction of wild-type mRNA was produced. Their mother was a heterozygous carrier, while their maternal grandmother was a carrier in form of gonosomal mosaicism. In summary, using a variety of techniques, including an NGS-based targeted gene panel and deep amplicon sequencing, copy number variation calling strategies, microarray-based comparative genomic hybridization, and cDNA analysis to define splicing defects and skewed XCI, we show how to face and solve some uncommon genetic mechanisms in the diagnosis of XL-CGD.

A mixed amplicon metabarcoding and sequencing approach for surveillance of drug resistance to levamisole and benzimidazole in Haemonchus spp.

Anthelmintic-resistant parasitic nematodes present a significant threat to sustainable livestock production worldwide. The ability to detect the emergence of anthelmintic resistance at an early stage, and therefore determine which drugs remain most effective, is crucial for minimising production losses. Despite many years of research into the molecular basis of anthelmintic resistance, no molecular-based tools are commercially available for the diagnosis of resistance as it emerges in field settings. We describe a mixed deep amplicon sequencing approach to determine the frequency of the levamisole (LEV)-resistant single nucleotide polymorphism (SNP) within arc-8 exon 4 (S168T) in Haemonchus spp., coupled with benzimidazole (BZ)-resistant SNPs within β-tubulin isotype-1 and the internal transcribed spacer-2 (ITS-2) nemabiome. This constitutes the first known multi-drug and multi-species molecular diagnostic developed for helminths of veterinary importance. Of the ovine, bovine, caprine and camelid Australian field isolates we tested, S168T was detected in the majority of Haemonchus spp. populations from sheep and goats, but rarely at a frequency greater than 16%; an arbitrary threshold we set based on whole genome sequencing (WGS) of LEV-resistant Haemonchus contortus GWBII. Overall, BZ resistance was far more prevalent in Haemonchus spp. than LEV resistance, confirming that LEV is still an effective anthelmintic class for small ruminants in New South Wales, Australia. The mixed amplicon metabarcoding approach described herein paves the way towards the use of large scale sequencing as a surveillance technology in the field, the results of which can be translated into evidence-based recommendations for the livestock sector.

Hi-plex deep amplicon sequencing for identification, high-resolution genotyping and multidrug resistance prediction of Mycobacterium leprae directly from patient biopsies by using Deeplex Myc-Lep.

Expansion of antimicrobial resistance monitoring and epidemiological surveillance are key components of the WHO strategy towards zero leprosy. The inability to grow Mycobacterium leprae invitro precludes routine phenotypic drug susceptibility testing, and only limited molecular tests are available. We evaluated a culture-free targeted deep sequencing assay, for mycobacterial identification, genotyping based on 18 canonical SNPs and 11 core variable-number tandem-repeat (VNTR) markers, and detection of rifampicin, dapsone and fluoroquinolone resistance-associated mutations in rpoB/ctpC/ctpI, folP1, gyrA/gyrB, respectively, and hypermutation-associated mutations in nth. The limit of detection (LOD) was determined using DNA of M.leprae reference strains and from 246 skin biopsies and 74 slit skin smears of leprosy patients, with genome copies quantified by RLEP qPCR. Sequencing results were evaluated versus whole genome sequencing (WGS) data of 14 strains, and versus VNTR-fragment length analysis (FLA) results of 89 clinical specimens. The LOD for sequencing success ranged between 80 and 3000 genome copies, depending on the sample type. The LOD for minority variants was 10%. All SNPs detected in targets by WGS were identified except in a clinical sample where WGS revealed two dapsone resistance-conferring mutations instead of one by Deeplex Myc-Lep, due to partial duplication of the sulfamide-binding domain in folP1. SNPs detected uniquely by Deeplex Myc-Lep were missed by WGS due to insufficient coverage. Concordance with VNTR-FLA results was 99.4% (926/932 alleles). Deeplex Myc-Lep may help improve the diagnosis and surveillance of leprosy. Gene domain duplication is an original putative drug resistance-related genetic adaptation in M.leprae. EDCTP2 programme supported by the European Union (grant number RIA2017NIM-1847 -PEOPLE). EDCTP, R2Stop: Effect:Hope, The Mission To End Leprosy, the Flemish Fonds Wetenschappelijk Onderzoek.

Polymorphism analysis of drug resistance markers in Plasmodium falciparum isolates from Benin

Like most countries in sub-Saharan African countries, Benin continues to bear a heavy malaria burden. In 2014, the National Malaria Control Programme (NMCP) changed its treatment policy, and recommended the use of artemisinin-based combination therapy (ACT) as first-line treatment for uncomplicated Plasmodium falciparum cases. The study presented here was conducted to investigate the impact of current antimalarial drug resistance on the country. Molecular surveillance targeting the Pfcrt, Pfmdr1, Pfkelch13, dhfr, and dhps genes was carried out on samples from patients positive for P. falciparum malaria by microscopy, LAMP and PCR diagnostic test. Molecular analysis was performed using targeted amplicon deep sequencing (TADS). In addition, the frequency of parasites with dual deletion of the histidine-rich protein 2 and 3 genes (pfhrp2 and pfhrp3), known to be responsible of the performance of HRP-based malaria rapid diagnostic tests (HRP-RDT), was estimated. Fifty-three falciparum samples collected at the Saint Jean de Dieu hospital in Tanguiéta, Benin, were tested. No Pfkelch13 validated or candidate artemisinin partial resistant variants were identified. A marked prevalence of Asn51Ile (N51I), Cys59Arg (C59R), and Ser108Asn (S108N) mutant alleles was found in the dhfr gene, representing the most frequent genotype (64%). Five-point mutations were detected in dhps, Ile431Val (I431V), Ser436Ala (S436A), Ala437Gly (A437G), Ala581Gly (A581G), Ala613Ser (A613S) of which the third was the most common (92%). No mutation was identified in dhps Lys540Glu (K540E). The quintuple mutant genotype resulting from the combination of the dhfr triple mutant (51I/59R/108N) with the dhps double mutant 436A/437G was detected at a frequency of 30%. Low levels of mutations in Pfcrt and no mutation at codon 86 in the Pfmdr1 DNA fragment were observed, whereas a high level of Tyr184Phe (Y184F) polymorphism in the Pfmdr1 gene was found. These results could be indicative, over a decade after the implementation of ACT therapy, of the return of chloroquine-sensitive but artemether-lumefantrine resistant falciparum genotypes in Benin. There was no evidence of HRP2 and HRP3 deletions. Data from the present study support the need for routine monitoring of molecular markers of antimalarial drug resistance as part of surveillance activities aimed to make informed treatment policy decisions at the national level.

Deep amplicon sequencing reveals extensive allelic diversity in the erg11/CYP51 promoter and allows multi-population DMI fungicide resistance monitoring in the canola pathogen Leptosphaeria maculans

Continued use of fungicides provides a strong selection pressure towards strains with mutations to render these chemicals less effective. Previous research has shown that resistance to the demethylation inhibitor (DMI) fungicides, which target ergosterol synthesis, in the canola pathogen Leptosphaeria maculans has emerged in Australia and Europe. The change in fungicide sensitivity of individual isolates was found to be due to DNA insertions into the promoter of the erg11/CYP51 DMI target gene. Whether or not these were the only types of mutations and how prevalent they were in Australian populations was explored in the current study. New isolates with reduced DMI sensitivity were obtained from screens on DMI-treated plants, revealing eight independent insertions in the erg11 promoter. A novel deep amplicon sequencing approach applied to populations of ascospores fired from stubble identified an additional undetected insertion allele and quantified the frequencies of all known insertions, suggesting that, at least in the samples processed, the combined frequency of resistant alleles is between 0.0376% and 32.6%. Combined insertion allele frequencies positively correlated with population-level measures of in planta resistance to four different DMI treatments. Additionally, there was no evidence for erg11 coding mutations playing a role in conferring resistance in Australian populations. This research provides a key method for assessing fungicide resistance frequency in stubble-borne populations of plant pathogens and a baseline from which additional surveillance can be conducted in L. maculans. Whether or not the observed resistance allele frequencies are associated with loss of effective disease control in the field remains to be established.

Exploring the mycobiota of bromeliads phytotelmata in Brazilian Campos Rupestres.

The phytotelmata is a water-filled tank on a terrestrial plant, and it plays an important role in bromeliad growth and ecosystem functioning. Even though previous studies have contributed to elucidate the composition of the prokaryotic component of this aquatic ecosystem, its mycobiota (fungal community) is still poorly known. In the present work, ITS2 amplicon deep sequencing was used to examine the fungal communities inhabiting the phytotelmata of two bromeliads species that coexist in a sun-exposed rupestrian field of Southeastern Brazil, namely Aechmea nudicaulis (AN) and Vriesea minarum (VM). Ascomycota was the most abundant phylum in both bromeliads (57.1 and 89.1% in AN and VM respectively, on average), while the others were present in low abundance (< 2%). Mortierellomycota and Glomeromycota were exclusively observed in AN. Beta-diversity analysis showed that samples from each bromeliad significantly clustered together. In conclusion, despite the considerable within-group variation, the results suggested that each bromeliad harbor a distinct fungi community, what could be associated with the physicochemical characteristics of the phytotelmata (mainly total nitrogen, total organic carbon, and total carbon) and plant morphological features.

Genetic variations of Plasmodium falciparum circumsporozoite protein and the impact on interactions with human immunoproteins and malaria vaccine efficacy.

In October 2021, the world's first malaria vaccine RTS,S was endorsed by WHO for broad use in children, despite its low efficacy. This study examined polyclonal infections and the associations of parasite genetic variations with binding affinity to human leukocyte antigen (HLA). Multiplicity of infection was determined by amplicon deep sequencing of PfMSP1. Genetic variations in PfCSP were examined across 88 samples from Ghana and analyzed together with 1655 PfCSP sequences from other African and non-African isolates. Binding interactions of PfCSP peptide variants and HLA were predicted using NetChop and HADDOCK. High polyclonality was detected among infections, with each infection harboring multiple non-3D7 PfCSP variants. Twenty-seven PfCSP haplotypes were detected in the Ghanaian samples, and they broadly represented PfCSP diversity across Africa. The number of genetic differences between 3D7 and non-3D7 PfCSP variants does not influence binding to HLA. However, CSP peptide length after proteolytic degradation significantly affects its molecular weight and binding affinity to HLA. Despite the high diversity of HLA, the majority of the HLAI and II alleles interacted/bound with all Ghana CSP peptides. Multiple non-3D7 strains among P. falciparum infections could impact the effectiveness of RTS,S. Longer peptides of the Th2R/Th3R CSP regions should be considered in future versions of RTS,S.

Amplicon Deep Sequencing Reveals Multiple Genetic Events Lead to Treatment Failure with Atovaquone-Proguanil in Plasmodium falciparum.

Atovaquone-proguanil (AP) is used as treatment for uncomplicated malaria, and as a chemoprophylactic agent against Plasmodium falciparum. Imported malaria remains one of the top causes of fever in Canadian returning travelers. Twelve sequential whole-blood samples before and after AP treatment failure were obtained from a patient diagnosed with P. falciparum malaria upon their return from Uganda and Sudan. Ultradeep sequencing was performed on the cytb, dhfr, and dhps markers of treatment resistance before and during the episode of recrudescence. Haplotyping profiles were generated using three different approaches: msp2-3D7 agarose and capillary electrophoresis, and cpmp using amplicon deep sequencing (ADS). A complexity of infection (COI) analysis was conducted. De novo cytb Y268C mutants strains were observed during an episode of recrudescence 17 days and 16 h after the initial malaria diagnosis and AP treatment initiation. No Y268C mutant reads were observed in any of the samples prior to the recrudescence. SNPs in the dhfr and dhps genes were observed upon initial presentation. The haplotyping profiles suggest multiple clones mutating under AP selection pressure (COI > 3). Significant differences in COI were observed by capillary electrophoresis and ADS compared to the agarose gel results. ADS using cpmp revealed the lowest haplotype variation across the longitudinal analysis. Our findings highlight the value of ultra-deep sequencing methods in the understanding of P. falciparum haplotype infection dynamics. Longitudinal samples should be analyzed in genotyping studies to increase the analytical sensitivity.

Angie-LAMP for diagnosis of human eosinophilic meningitis using dog as proxy: A LAMP assay for Angiostrongylus cantonensis DNA in cerebrospinal fluid

Angiostrongylus cantonensis (rat lungworm) is recognised as the leading cause of human eosinophilic meningitis, a serious condition observed when nematode larvae migrate through the CNS. Canine Neural Angiostrongyliasis (CNA) is the analogous disease in dogs. Both humans and dogs are accidental hosts, and a rapid diagnosis is warranted. A highly sensitive PCR based assay is available but often not readily accessible in many jurisdictions. An alternative DNA amplification assay that would further improve accessibility is needed. This study aimed to assess the diagnostic utility of a newly designed LAMP assay to detect DNA of globally distributed and invasive A. cantonensis and Angiostrongylus mackerrasae, the other neurotropic Angiostrongylus species, which is native to Australia. Cerebrospinal fluid (CSF) from dogs with a presumptive diagnosis of A. cantonensis infection (2020-2022) were received for confirmatory laboratory testing and processed for DNA isolation and ultrasensitive Angiostrongylus qPCR targeting AcanR3390. A newly designed LAMP assay targeting the same gene target was directly compared to the reference ultrasensitive qPCR in a diagnostic laboratory setting to determine the presence of A. cantonensis DNA to diagnose CNA. The LAMP assay (Angie-LAMP) allowed the sensitive detection of A. cantonensis DNA from archived DNA specimens (Kappa = 0.81, 95%CI 0.69-0.92; n = 93) and rapid single-step lysis of archived CSF samples (Kappa = 0.77, 95%CI 0.59-0.94; n = 52). Only A. cantonensis DNA was detected in canine CSF samples, and co-infection with A. mackerrasae using amplicon deep sequencing (ITS-2 rDNA) was not demonstrated. Both SYD.1 and AC13 haplotypes were detected using sequencing of partial cox1. The Angie-LAMP assay is a useful molecular tool for detecting Angiostrongylus DNA in canine CSF and performs comparably to a laboratory Angiostrongylus qPCR. Adaptation of single-step sample lysis improved potential applicability for diagnosis of angiostrongyliasis in a clinical setting for dogs and by extension, to humans.

Single nucleotide variants in Pseudomonas aeruginosa populations from sputum correlate with baseline lung function and predict disease progression in individuals with cystic fibrosis.

The severity and progression of lung disease are highly variable across individuals with cystic fibrosis (CF) and are imperfectly predicted by mutations in the human gene CFTR, lung microbiome variation or other clinical factors. The opportunistic pathogen Pseudomonas aeruginosa (Pa) dominates airway infections in most CF adults. Here we hypothesized that within-host genetic variation of Pa populations would be associated with lung disease severity. To quantify Pa genetic variation within CF sputum samples, we used deep amplicon sequencing (AmpliSeq) of 209 Pa genes previously associated with pathogenesis or adaptation to the CF lung. We trained machine learning models using Pa single nucleotide variants (SNVs), microbiome diversity data and clinical factors to classify lung disease severity at the time of sputum sampling, and to predict lung function decline after 5 years in a cohort of 54 adult CF patients with chronic Pa infection. Models using Pa SNVs alone classified lung disease severity with good sensitivity and specificity (area under the receiver operating characteristic curve: AUROC=0.87). Models were less predictive of lung function decline after 5 years (AUROC=0.74) but still significantly better than random. The addition of clinical data, but not sputum microbiome diversity data, yielded only modest improvements in classifying baseline lung function (AUROC=0.92) and predicting lung function decline (AUROC=0.79), suggesting that Pa AmpliSeq data account for most of the predictive value. Our work provides a proof of principle that Pa genetic variation in sputum tracks lung disease severity, moderately predicts lung function decline and could serve as a disease biomarker among CF patients with chronic Pa infections.

Evolution of Plasmodium falciparum antimalarial drug resistance markers post-adoption of artemisinin-based combination therapies in Yaounde, Cameroon

ObjectivesThis study aimed to investigate the evolution of Plasmodium falciparum antimalarial drug resistance markers by comparing the pre- and post-adoption of artemisinin-based combination therapies (ACTs) in Yaounde, Cameroon. MethodsThe molecular characterization of known antimalarial drug resistance markers (Pfcrt, Pfmdr1, Pfdhfr, Pfdhps, and Pfk13) in P. falciparum-positive samples collected in 2014 and 2019-2020 was achieved using nested polymerase chain reaction, followed by targeted amplicon deep sequencing on the Illumina MiSeq platform. Data derived were compared with those published during the pre-ACT adoption period from 2004 to 2006. ResultsA high prevalence of Pfmdr1 184F, Pfdhfr 51I/59R/108N, and Pfdhps 437G mutant alleles was observed during the post-ACT adoption period. The Pfcrt 76T and Pfmdr1 86Y mutant alleles significantly declined between 2004 and 2020 (P <0.0001). Conversely, the resistance markers to antifolates, Pfdhfr 51I/59R/108N and Pfdhps 437G, significantly increased during the same study period (P <0.0001). We identified nine mutations in the propeller domains of Pfk13; although they were all present in single parasite isolates, none of them are known to confer artemisinin resistance. ConclusionThis study documented a near-complete reversion to sensitive parasites for markers conferring resistance to the 4-aminoquinolines and arylamino alcohols in Yaounde. In contrast, the Pfdhfr mutations associated with pyrimethamine resistance are moving toward saturation.

Application of a universal parasite diagnostic test to biological specimens collected from animals.

A previously described universal parasite diagnostic (nUPDx) based on PCR amplification of the 18S rDNA and deep-amplicon sequencing, can detect human blood parasites with a sensitivity comparable to real-time PCR. To date, the efficacy of this assay has only been assessed on human blood. This study assessed the utility of nUPDx for the detection of parasitic infections in animals using blood, tissues, and other biological sample types from mammals, birds, and reptiles, known to be infected with helminth, apicomplexan, or pentastomid parasites (confirmed by microscopy or PCR), as well as negative samples. nUPDx confirmed apicomplexan and/or nematode infections in 24 of 32 parasite-positive mammals, while also identifying several undetected coinfections. nUPDx detected infections in 6 of 13 positive bird and 1 of 2 positive reptile samples. When applied to 10 whole parasite specimens (worms and arthropods), nUPDx identified all to the genus or family level, and detected one incorrect identification made by morphology. Babesia sp. infections were detected in 5 of the 13 samples that were negative by other diagnostic approaches. While nUPDx did not detect PCR/microscopy-confirmed trichomonads or amoebae in cloacal swabs/tissue from 8 birds and 2 reptiles due to primer template mismatches, 4 previously undetected apicomplexans were detected in these samples. Future efforts to improve the utility of the assay should focus on validation against a larger panel of tissue types and animal species. Overall, nUPDx shows promise for use in both veterinary diagnostics and wildlife surveillance, especially because species-specific PCRs can miss unknown or unexpected pathogens.

Competitive survival of clonal serial Pseudomonas aeruginosa isolates from cystic fibrosis airways in human neutrophils

Chronic airway infections with Pseudomonas aeruginosa are the major co-morbidity in most people with cystic fibrosis (CF) sustained by neutrophils as the major drivers of lung inflammation, damage, and remodeling. Phagocytosis assays were performed with clonal consortia of longitudinal P.aeruginosa airway isolates collected from people with CF since the onset of lung colonization until patient's death or replacement by another clone. The extra- and intracellular abundance of individual strains was assessed by deep amplicon sequencing of strain-specific single nucleotide variants in the bacterial genome. The varied microevolution of the accessory genome of the P.aeruginosa clones during mild and severe courses of infection corresponded with a differential persistence of clonal progeny in the neutrophil phagosome. By simultaneously exposing the ancestor and its progeny to the same habitat, the study recapitulated the time lapse of the temporal change of the fitness of the clone to survive in neutrophils.

AmpSeqR: an R package for amplicon deep sequencing data analysis

Amplicon sequencing (AmpSeq) is a methodology that targets specific genomic regions of interest for polymerase chain reaction (PCR) amplification so that they can be sequenced to a high depth of coverage. Amplicons are typically chosen to be highly polymorphic, usually with several highly informative, high frequency single nucleotide polymorphisms (SNPs) segregating in an amplicon of 100–200 base pair (bp). This allows high sensitivity detection and quantification of the frequency of each sequence within each sample making it suitable for applications such as low frequency somatic mosaicism detection or minor clone detection in mixed samples. AmpSeq is being increasingly applied to both biological and medical studies, in applications such as cancer, infectious diseases and brain mosaicism studies. Current bioinformatics pipelines for AmpSeq data processing lack downstream analysis, have difficulty distinguishing between true sequences and PCR sequencing errors and artifacts, and often require bioinformatic expertise. We present a new R package: AmpSeqR, designed for the processing of deep short-read amplicon sequencing data, with a focus on infectious diseases. The pipeline integrates several existing R packages combining them with newly developed functions to perform optimal filtering of reads to remove noise and improve the accuracy of the detected sequences data, permitting detection of very low frequency clones in mixed samples. The package provides useful functions including data pre-processing, amplicon sequence variants (ASVs) estimation, data post-processing, data visualization, and automatically generates a comprehensive Rmarkdown report that contains all essential results facilitating easy inclusion into reports and publications. AmpSeqR is publicly available at https://github.com/bahlolab/AmpSeqR.

Molecular evidence of widespread benzimidazole drug resistance in Ancylostoma caninum from domestic dogs throughout the USA and discovery of a novel β-tubulin benzimidazole resistance mutation.

Ancylostoma caninum is an important zoonotic gastrointestinal nematode of dogs worldwide and a close relative of human hookworms. We recently reported that racing greyhound dogs in the USA are infected with A. caninum that are commonly resistant to multiple anthelmintics. Benzimidazole resistance in A. caninum in greyhounds was associated with a high frequency of the canonical F167Y(TTC>TAC) isotype-1 β-tubulin mutation. In this work, we show that benzimidazole resistance is remarkably widespread in A. caninum from domestic dogs across the USA. First, we identified and showed the functional significance of a novel benzimidazole isotype-1 β-tubulin resistance mutation, Q134H(CAA>CAT). Several benzimidazole resistant A. caninum isolates from greyhounds with a low frequency of the F167Y(TTC>TAC) mutation had a high frequency of a Q134H(CAA>CAT) mutation not previously reported from any eukaryotic pathogen in the field. Structural modeling predicted that the Q134 residue is directly involved in benzimidazole drug binding and that the 134H substitution would significantly reduce binding affinity. Introduction of the Q134H substitution into the C. elegans β-tubulin gene ben-1, by CRISPR-Cas9 editing, conferred similar levels of resistance as a ben-1 null allele. Deep amplicon sequencing on A. caninum eggs from 685 hookworm positive pet dog fecal samples revealed that both mutations were widespread across the USA, with prevalences of 49.7% (overall mean frequency 54.0%) and 31.1% (overall mean frequency 16.4%) for F167Y(TTC>TAC) and Q134H(CAA>CAT), respectively. Canonical codon 198 and 200 benzimidazole resistance mutations were absent. The F167Y(TTC>TAC) mutation had a significantly higher prevalence and frequency in Western USA than in other regions, which we hypothesize is due to differences in refugia. This work has important implications for companion animal parasite control and the potential emergence of drug resistance in human hookworms.

Exploring how space, time, and sampling impact our ability to measure genetic structure across Plasmodium falciparum populations.

A primary use of malaria parasite genomics is identifying highly related infections to quantify epidemiological, spatial, or temporal factors associated with patterns of transmission. For example, spatial clustering of highly related parasites can indicate foci of transmission and temporal differences in relatedness can serve as evidence for changes in transmission over time. However, for infections in settings of moderate to high endemicity, understanding patterns of relatedness is compromised by complex infections, overall high forces of infection, and a highly diverse parasite population. It is not clear how much these factors limit the utility of using genomic data to better understand transmission in these settings. In particular, further investigation is required to determine which patterns of relatedness we expect to see with high quality, densely sampled genomic data in a high transmission setting and how these observations change under different study designs, missingness, and biases in sample collection. Here we investigate two identity-by-state measures of relatedness and apply them to amplicon deep sequencing data collected as part of a longitudinal cohort in Western Kenya that has previously been analysed to identify individual-factors associated with sharing parasites with infected mosquitoes. With these data we use permutation tests, to evaluate several hypotheses about spatiotemporal patterns of relatedness compared to a null distribution. We observe evidence of temporal structure, but not of fine-scale spatial structure in the cohort data. To explore factors associated with the lack of spatial structure in these data, we construct a series of simplified simulation scenarios using an agent based model calibrated to entomological, epidemiological and genomic data from this cohort study to investigate whether the lack of spatial structure observed in the cohort could be due to inherent power limitations of this analytical method. We further investigate how our hypothesis testing behaves under different sampling schemes, levels of completely random and systematic missingness, and different transmission intensities.

Baseline prevalence of molecular marker of sulfadoxine/pyrimethamine resistance in Ebonyi and Osun states, Nigeria: amplicon deep sequencing of dhps-540.

Chemoprevention plays an important role in malaria control strategy. Perennial malaria chemoprevention (PMC) using sulfadoxine/pyrimethamine (SP) is a WHO-approved strategy to combat malaria in young children and may lead to drug pressure. Introducing SP-PMC may therefore be compromised due to the emergence of Plasmodium falciparum resistant to SP, particularly mutation at K540E of the dihydropteroate synthase (dhps) gene. Molecular surveillance of resistance markers can support assessment of antimalarial efficacy and effectiveness. High prevalence of 540E is associated with reduced effectiveness of SP, and areas with more than 50% prevalence are considered unsuitable for intermittent preventative treatment in pregnancy (IPTp) implementation. Assessing 540E prevalence is an important undertaking before implementation of SP-PMC. We conducted a rapid surveillance of dhps-540E to assess the suitability of SP as PMC in field studies from Ebonyi and Osun states in Nigeria. We used an in-house developed amplicon deep-sequencing method targeting part of the dhps gene. Our data reveal that 18.56% of individuals evaluated carried the 540E mutation mixed with the WT K540. Mutant variant 540E alone was not found, and 80% of isolates harboured only WT (K540). Clonal analysis of the sequencing data shows a very low proportion of 540E circulating in both states. Our data show that both states are suitable for SP-PMC implementation and, based on this finding, SP-PMC was implemented in Osun in 2022. Continuous monitoring of 540E will be required to ensure the chemoprevention effectiveness of SP in Nigeria.

Insertion of ten amino acids into 13S globulin zero-repeat subunit improves trypsin digestibility in common buckwheat (Fagopyrum esculentum Moench) seeds

The 13S globulin zero-repeat subunit is resistant to trypsin and may have higher allergenicity than the 1–6 tandem repeat subunits in common buckwheat (Fagopyrum esculentum Moench). To explore alleles useful for lowering allergenicity, amplicon deep sequencing targeting the zero-repeat subunit gene was conducted in bulked genomic DNA from eight cultivars and landraces. The analysis identified a unique allele encoding a zero-repeat subunit with 10 amino acid insertion (10aa) at a position equivalent to the tandem repeat insertion. Prediction of its 3-D structure suggested that 10aa changes the β-hairpin structure in the non-10aa (native) subunit to a random coil, which is also found in 1- and 3- repeat subunits. Homozygotes of the 10aa allele were developed and showed that the 10aa subunit was more digestible than the native subunit. However, the 10aa subunit was still less digestible than the 1–6 repeat subunits, suggesting needs to explore unfunctional alleles.