Deep Amplicon Sequencing Research Articles

Monitoring molecular markers associated with antimalarial drug resistance in south-east Senegal from 2021 to 2023.

Since 2006, artemisinin-based combination therapies (ACTs) have been introduced in Senegal in response to chloroquine resistance (CQ-R) and have shown high efficacy against Plasmodium falciparum. However, the detection of the PfKelch13R515K mutation in Kaolack, which confers artemisinin resistance in vitro, highlights the urgency of strengthening antimalarial drug surveillance to achieve malaria elimination by 2030. To assess the proportion of P. falciparum parasites carrying molecular signatures associated with antimalarial resistance (PfKelch13, Pfmdr1, Pfcrt, dhfr and dhps) in isolates collected at Kédougou using multiplex amplicon deep sequencing. Venous blood samples were collected from patients diagnosed with P. falciparum infection over a 3-year period (2021, 2022 and 2023). Parasite DNA was extracted, and multiplex amplicon sequencing was used to investigate gene polymorphisms. Analysis of PfKelch13 did not reveal any non-synonymous mutations. Pfcrt mutations were present in 45% of the samples, mainly K76T (44%) and I356T (36%). The dominant Pfmdr-1 allele was Y184F (62%). The sextuple mutant 51I/59R/108N + 436A/437G/613S dhfr/dhps was observed in 10% of the samples. The absence of PfKelch13 mutants suggests that ACT efficacy remains uncompromised, although clinical outcome studies are required to confirm this. Analysis of Pfcrt and Pfmdr-1 shows that CQ-R alleles, probably from previous CQ use, are slowly decreasing. Likewise, the detection of the dhfr/dhps sextuple mutant highlights the need to monitor sulfadoxine-pyrimethamine resistance and the emergence of 581G. There is therefore a need for continued antimalarial resistance surveillance in Senegal.

Molecular epidemiology and genetic diversity of disappearing Plasmodium vivax in southern Thailand

The evolution of genetic diversity and population structure of Plasmodium vivax as malaria elimination approaches remains unclear. This study analyzed the genetic variation and molecular epidemiology of P. vivax from Yala Province in southern Thailand, an area in the pre-elimination phase. Seventy P. vivax isolates, collected between 2017 and 2020, were genotyped for domain II of pvdbp and the 42-kDa region of pvmsp1 using amplicon deep sequencing. Data from Yala province were compared to published data from Tak province, where transmission was higher. Key analyses included nucleotide diversity (π), haplotype diversity (Hd), natural selection, recombination rates, and complexity of infection (COI). Genetic diversity in Yala was relatively low (π = 0.008dbp and 0.014msp1; Hd = 0.774dbp and 0.407msp1) compared to Tak (π = 0.012dbp and 0.027msp1; Hd = 0.849dbp and 0.962msp1). In Yala, polyclonal infections were found in 53.7% of pvdbpII and 47.8% of pvmsp142 isolates, with average COI of 1.6 and 1.7. Both genes were under balancing selection. Distinct genetic differences were found between Yala and Tak in pvmsp142, providing a local genotypic profile useful for tracing parasite origins.

Photodynamic Therapy with Protoporphyrin IX Precursors Using Artificial Daylight Improves Skin Antisepsis for Orthopedic Surgeries.

Classical preoperative skin antisepsis is insufficient in completely eliminating bacterial skin colonization for arthroplasty. In contrast, photodynamic therapy (PDT) with red light and methyl-aminolevulinate (MAL), combined with skin antisepsis, led to the absence of bacterial growth in healthy participants, though with local skin erythema, posing an obstacle for orthopedic surgery. Therefore, we explored whether artificial daylight PDT (PDT-DL) was superior to red light. Twenty healthy participants were allocated to either 5-aminolevulinic acid-(5-ALA) PDT-DL (n = 10) or MAL-PDT-DL (n = 10) before antisepsis with povidone-iodine/alcohol. Skin swabs from the groin were taken to cultivate bacteria at baseline, after PDT-DL, and after the subsequent antisepsis. Additional swabs were taken on day 4 before and after antisepsis without PDT. The contralateral groin of each participant and of ten additional healthy volunteers served as the control (n = 30). In selected participants, 16S rRNA-based amplicon deep sequencing was performed. All participants showed a baseline bacterial colonization. After a PDT-DL with skin antisepsis, bacterial growth occurred in three (30%) and in one (10%) participants with 5-ALA and MAL, respectively, compared to the sixteen (55%) participants in the control group. On day 4, three (30%) participants per group showed positive cultures post antisepsis. Adverse effects were reported in six (60%) and zero (0%) participants for 5-ALA- and MAL-PDT-DL, respectively. The skin bacteriome changes correlated with the bacterial culture results. The MAL-PDT-DL with skin antisepsis significantly increased bacterial reduction on the skin without adverse effects. This offers an opportunity to prevent infections in arthroplasty patients and reduce antibiotic use, thus contributing to antibiotic stewardship goals emphasized in the One Health approach.

Detection of low frequency artemisinin resistance mutations, C469Y, P553L and A675V, and fixed antifolate resistance mutations in asymptomatic primary school children in Kenya

BackgroundTo understand the emergence and spread of drug-resistant parasites in malaria-endemic areas, accurate assessment and monitoring of antimalarial drug resistance markers is critical. Recent advances in next-generation sequencing (NGS) technologies have enabled the tracking of drug-resistant malaria parasites.MethodsIn this study, we used Targeted Amplicon Deep Sequencing (TADS) to characterise the genetic diversity of the Pfk13, Pfdhfr, Pfdhps, and Pfmdr1 genes among primary school-going children in 15 counties in Kenya (Bungoma, Busia, Homa Bay, Migori, Kakamega, Kilifi, Kirinyaga, Kisii, Kisumu, Kwale, Siaya, Tana River, Turkana, Vihiga and West Pokot). A total of 920 dried blood spot (DBS) samples collected from 121 selected primary schools within the country were used to extract genomic DNA. A nested polymerase chain reaction (PCR) was used to generate amplicons that were sequenced to determine the prevalence of known and novel polymorphisms.ResultsPfk13 mutations associated with artemisinin resistance were present as mixed genotype infections for the C469Y mutation in 23 samples (4%), the A675V mutation in 2 samples (1.7%), and the P553L mutation in 7 samples (1.2%). The A578S mutation, was also identified in mixed infections, appearing in 15.2% of the 87 samples analysed. The Pfdhfr 51I and 108 N pyrimethamine-resistance mutations were at fixation (100% frequency), and the Pfmdr1 Y184F mutation, which is linked to reduced susceptibility to several antimalarial drugs, especially those used in combination therapies for malaria treatment, was detected in 97.5% of the samples as mixed-genotype infections.ConclusionThe genomic surveillance of asymptomatic school children in Kenya provides an early warning signal of at least 1 of the 3 validated artemisinin resistance mutations circulating in all regions in Western Kenya sampled except Homa Bay and Kisii Counties. These signals in asymptomatic and mixed infections would have been missed without deep sequencing.

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

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

Adaptation of bacterial natural single guide RNA (tracr-L) for efficient plant genome editing.

A long tracrRNA (tracr-L),which naturally act as single guide RNA,and its truncated version, Δtracr-L, from S. pyogenes, efficiently induce Cas9-mediated double-strand breaks (DSBs) in plant genomic loci, as demonstrated by in vitro cleavage assay and protoplast transfection. CRISPR-Cas system provides a form of immune memory in prokaryotes and archaea, protecting them against viruses and foreign genetic elements. In Streptococcus pyogenes, this system includes the pre-crRNA along with another non-coding RNA, tracrRNA, which aids in CRISPR-based immunity. In S. pyogenes, two distinct tracrRNAs are produced: a long form (tracr-L) and a short form (tracr-S). The tracr-S regulates crRNA biogenesis and Cas9 cleavage, while tracr-L suppresses CRISPR-Cas expression by targeting the Cas9 promoter to prevent autoimmunity. Deleting 79 nucleotides from tracr-L results in Δtracr-L, which retains similar functionality in gene repression. This study investigates, for thefirst time, the effectiveness of tracr-L, and Δtracr-L in genome editing within plant systems. In vitro cleavage assays using purified Cas9 and synthesized sgRNAs targeting the Cas9 gene, OsPDS, and the OsSWEET11 promoter revealed that across all target sites, tracr-S demonstrated the highest cleavage efficiency compared to tracr-L and Δtracr-L. For in vivo genome editing, we transfected rice protoplasts with tracr-L, Δtracr-L, and tracr-S, targeting three rice genes: OsPDS, OsSPL14, and the promoter of OsSWEET14. Amplicon deep sequencing revealed various types of indels at the target regions across all three tracrRNA versions, indicating comparable levels of efficiency. This study establishes the utility of both the long-form tracrRNA (tracr-L) and its truncated variant (Δtracr-L) in eukaryote genome editing. These two new forms of tracrRNA provide proof of concept and expand the CRISPR-Cas toolkit for plant genome editing applications, and for eukaryotes more broadly.

Genetic diversity of Plasmodium vivax population in northeast Myanmar assessed by amplicon sequencing of PvMSP1 and PvMSP3α

This study aimed to assess the baseline genetic diversity of the Plasmodium vivax population in an endemic area of northeast Myanmar at the onset of the malaria elimination campaign in the Greater Mekong Subregion. We genotyped 125 P vivax clinical samples at two merozoite surface protein loci, PvMSP1 and PvMSP3α, by amplicon deep sequencing. Our study revealed that the parasite population in this region was highly diverse, identifying 60 PvMSP1 and 98 PvMSP3α haplotypes, with haplotype diversity of 0.929 and 0.944, respectively. Remarkably, 97.6 % (122/125) of the patients harbored multiclonal infections with a mean multiplicity of infection of 4.18, indicating a relatively high transmission intensity. Neutrality tests and network analysis suggested a recent parasite population expansion, consistent with the concurrent malaria outbreak in the region. These findings underscore the existence of a highly diverse P. vivax population at the China-Myanmar border, highlighting the need for effective malaria control strategies to achieve the goal of regional malaria elimination.

Deep-amplicon sequencing of the complete beta-tubulin gene in Trichuris trichiura before and after albendazole treatment

Concerns about the emergence of benzimidazole resistance in soil-transmitted helminths (STH) infections, particularly against Trichuris trichiura, have arisen. Previous studies of veterinary nematodes have linked benzimidazole resistance to single-nucleotide polymorphisms (SNPs) at three specific codons in the beta-tubulin gene, but similar associations in STH have not been consistently observed. In this work, we screened the complete beta-tubulin gene previously linked to benzimidazole resistance in T. trichiura by deep-amplicon sequencing to identify genetic variants and associate levels of diversity with drug response to albendazole. We used 99 DNA samples extracted from T. trichiura pooled eggs, previously semi-purified from human stool samples collected in Manhiça district, Mozambique. We obtained a set of 39 amplicons of the complete gene by subjecting the pooled eggs to long-read PCR and subsequently sequencing them. Of those amplicons, 22 and 17 were obtained from stool samples collected before, and 21 days after albendazole treatment, respectively. We observed genetic variation across the whole gene sequence, in both exons and introns; however, none were associated with the previously proposed resistance-associated SNPs, and none were predicted to significantly affect protein function. No significant differences in genetic diversity were observed between pre- and post-treatment samples. Using publicly available genome-wide data, we also analysed a second beta-tubulin isotype in the T. trichiura genome. We focused on detecting the canonical SNPs and assessing for signatures of genetic selection around this second isotype gene. This analysis did not reveal evidence supporting this second isotype's role in anthelmintic resistance. Despite the limitations of our study, such as a small sample size, particularly paired pre- and post-treatment samples (n = 6), or a restricted geographical area, we found no evidence linking either of the two beta-tubulin genes to benzimidazole resistance in T. trichiura, suggesting that genetic markers of drug resistance likely exist outside the beta-tubulin genes.

Application of a new highly multiplexed amplicon sequencing tool to evaluate Plasmodium falciparum antimalarial resistance and relatedness in individual and pooled samples from Dschang, Cameroon.

Resistance to antimalarial drugs remains a major obstacle to malaria elimination. Multiplexed, targeted amplicon sequencing is being adopted for surveilling resistance and dissecting the genetics of complex malaria infections. Moreover, genotyping of parasites and detection of molecular markers drug resistance in resource-limited regions requires open-source protocols for processing samples, using accessible reagents, and rapid methods for processing numerous samples including pooled sequencing. P lasmodium f alciparum Streamlined Multiplex Antimalarial Resistance and Relatedness Testing (Pf-SMARRT) is a PCR-based amplicon panel consisting of 15 amplicons targeting antimalarial resistance mutations and 9 amplicons targeting hypervariable regions. This assay uses oligonucleotide primers in two pools and a non-proprietary library and barcoding approach. We evaluated Pf-SMARRT using control mocked dried blood spots (DBS) at varying levels of parasitemia and a mixture of 3D7 and Dd2 strains at known frequencies, showing the ability to genotype at low parasite density and recall within-sample allele frequencies. We then piloted Pf-SMARRT to genotype 100 parasite isolates collected from uncomplicated malaria cases at three health facilities in Dschang, Western Cameroon. Antimalarial resistance genotyping showed high levels of sulfadoxine-pyrimethamine resistance mutations, including 31% prevalence of the DHPS A613S mutation. No K13 candidate or validated artemisinin partial resistance mutations were detected, but one low-level non-synonymous change was observed. Pf-SMARRT's hypervariable targets, used to assess complexity of infections and parasite diversity and relatedness, showed similar levels and patterns compared to molecular inversion probe (MIP) sequencing. While there was strong concordance of antimalarial resistance mutations between individual samples and pools, low-frequency variants in the pooled samples were often missed. Overall, Pf-SMARRT is a robust tool for assessing parasite relatedness and antimalarial drug resistance markers from both individual and pooled samples. Control samples support that accurate genotyping as low as 1 parasite per microliter is routinely possible.

Genotyping methods to distinguish Plasmodium falciparum recrudescence from new infection for the assessment of antimalarial drug efficacy: an observational, single-centre, comparison study

Distinguishing Plasmodium falciparum recrudescence from new infections is crucial for the assessment of antimalarial drug efficacy against Pfalciparum. We aimed to compare the efficacy of different genotyping methods to assess their effect on drug efficacy estimates, particularly in patients from high-transmission settings with polyclonal infections. In this head-to-head comparison study, we compared five different genotyping methods currently used: fast capillary electrophoresis (F-CE) using msp1, msp2, and glurp; high-resolution capillary electrophoresis (H-CE) using msp1, msp2, and glurp; H-CE using microsatellites; targeted amplicon deep sequencing (TADS) using single nucleotide polymorphism (SNP)-rich markers; and high-resolution melting (HRM) analysis using msp1 and msp2. We assessed their sensitivity in detecting minority clones in polyclonal infections, their reproducibility, and the genetic diversity of the markers used. Our study used four well characterised Pfalciparum laboratory strains mixed in varying ratios, and 20paired samples collected from an in-vivo clinical trial. The experiments were performed at the Swiss Tropical and Public Health Institute in Basel, Switzerland between May 5, 2020, and Aug 23,2021. H-CE using msp1 and msp2 and TADS revealed the highest sensitivity in detecting minority clones (up to ratios of 1:100 for H-CE and 50:1:1:1 for TADS in the FCB1:HB3 and 3D7:K1:HB3:FCB1 laboratory strain mixtures, respectively), highest reproducibility (intra-assay: 99% and 91% for H-CE and TADS, respectively; inter-assay: 98% and 92% for H-CE and TADS, respectively), and highest genetic diversity in the used markers (up to 36and32unique genotypes in 20paired samples for H-CE using msp2 and TADS using cpmp, respectively). Microsatellites assessed by H-CE had a lower genetic diversity compared with msp1, msp2, and glurp assessed by H-CE and the SNP-rich markers assessed by TADS, with a maximum of 13unique genotypes, and some genotypes having allelic frequencies larger than 30%.Markers used by TADS gave the most consistent results in distinguishing recrudescence from new infection across all methods (in 18of 20pairs of samples vs 15of 20pairs for H-CE). WHO currently recommends replacing glurp with microsatellites. However, in this study, the replacement of glurp with microsatellites did not change the genotyping outcome, probably due to the lower genetic diversity of microsatellites. More studies with large sample sizes are required to identify the most suitablemicrosatellites that could replace glurp. Our study indicates that TADS should be considered the gold standardfor genotyping to distinguish recrudescence from new infection, and that it should be used to validate other methods. Swiss Tropical and Public Health Institute.

Application of multiple mosaic callers improves post-zygotic mutation detection from exome sequencing data

PurposeThe gold standard for identification of post-zygotic variants (PZVs) is droplet digital polymerase chain reaction or high-depth sequencing across multiple tissues types. These approaches are yet to be systematically implemented for monogenic disorders. We developed PZV detection pipelines for correct classification of de novo variants. MethodOur pipelines detect PZV in parents (gonosomal mosaicism [pGoM]) and children (somatic mosaicism, “M3”). We applied them to research exome sequencing (ES) data from the Australian Cerebral Palsy Biobank (n = 145 trios) and Simons Simplex Collection (n = 405 families). Candidate mosaic variants were validated using deep amplicon sequencing or droplet digital polymerase chain reaction. Results69.2% (M3trio), 63.9% (M3single), and 92.7% (pGoM) of detected variants were validated, with 48.6%, 56.7%, and 26.2% of variants, respectively, meeting strict criteria for mosaicism. In the Australian Cerebral Palsy Biobank, 16.6% of probands and 20.7% of parents had at least 1 true-positive somatic or pGoM variant, respectively. A large proportion of PZVs detected in Simons Simplex Collection parents (79.8%) and child (94.5%) were not previously reported. We reclassified 3.7% to 8.0% of germline de novo variants as mosaic. ConclusionMany PZVs were incorrectly classified as germline variants or missed by previous approaches. Systematic application of our pipelines could increase genetic diagnostic rate, improve estimates of recurrence risk in families, and benefit novel disease gene identification.

Allele-specific CRISPR-Cas9 editing inactivates a single nucleotide variant associated with collagen VI muscular dystrophy

The application of allele-specific gene editing tools can expand the therapeutic options for dominant genetic conditions, either via gene correction or via allelic gene inactivation in situations where haploinsufficiency is tolerated. Here, we used allele-targeted CRISPR/Cas9 guide RNAs (gRNAs) to introduce inactivating frameshifting indels at a single nucleotide variant in the COL6A1 gene (c.868G>A; G290R), a variant that acts as dominant negative and that is associated with a severe form of congenital muscular dystrophy. We expressed spCas9 along with allele-targeted gRNAs, without providing a repair template, in primary fibroblasts derived from four patients and one control subject. Amplicon deep-sequencing for two gRNAs tested showed that single nucleotide deletions accounted for the majority of indels introduced. While activity of the two gRNAs was greater at the G290R allele, both gRNAs were also active at the wild-type allele. To enhance allele-selectivity, we introduced deliberate additional mismatches to one gRNA. One of these optimized gRNAs showed minimal activity at the WT allele, while generating productive edits and improving collagen VI matrix in cultured patient fibroblasts. This study strengthens the potential of gene editing to treat dominant-negative disorders, but also underscores the challenges in achieving allele selectivity with gRNAs.

Origins, diversity, and adaptive evolution of DWV in the honey bees of the Azores: the impact of the invasive mite Varroa destructor.

Deformed wing virus (DWV) is a honey bee virus, whose emergence from relative obscurity is driven by the recent host-switch, adaptation, and global dispersal of the ectoparasitic mite Varroa destructor (a highly efficient vector of DWV) to reproduction on honey bees (Apis mellifera). Our study examines how varroa affects the continuing evolution of DWV, using the Azores archipelago, where varroa is present on only three out of the eight Islands, as a natural experimental system for comparing different evolutionary conditions and trajectories. We combined qPCR of 494 honey bee colonies sampled across the archipelago with amplicon deep sequencing to reveal how the DWV genetic landscape is altered by varroa. Two of the varroa-free Islands were also free of DWV, while a further two Islands were intriguingly dominated by the rare DWV-C major variant. The other four Islands, including the three varroa-infested Islands, were dominated by the common DWV-A and DWV-B variants. The varroa-infested Islands had, as expected, an elevated DWV prevalence relative to the uninfested Islands, but not elevated DWV loads, due the relatively high prevalence and loads of DWV-C on the varroa-free Islands. This establishes the Azores as a stable refuge for DWV-C and provides the most convincing evidence to date that at least some major strains of DWV may be capable of not just surviving, but actually thriving in honey bees in the absence of varroa-mediated transmission. We did not detect any change in DWV genetic diversity associated with island varroa status but did find a positive association of DWV diversity with virus load, irrespective of island varroa status.

Bidirectional transfer of human cytomegalovirus strains in donor and recipient seropositive lung transplant patients

AbstractDonor and recipient human cytomegalovirus (HCMV) seropositive (D+R+) lung transplant recipients (LTRs) often harbor multiple strains of HCMV, likely due to transmitted donor (D) strains and reactivated recipient (R) strains. To date, the extent and timely occurrence of each likely source in shaping the post‐transplantation (post‐Tx) strain population is unknown. Here, we deciphered the D and R origin of the post‐Tx HCMV strain composition in blood, bronchoalveolar lavage (BAL), and CD45+ BAL cell subsets. We investigated either D and/or R formalin‐fixed paraffin‐embedded blocks or fresh D lung tissue from four D+R+ LTRs obtained before transplantation. HCMV strains were characterized by short amplicon deep sequencing. In two LTRs, we show that the transplanted lung is reseeded by R strains within the first 6 months after transplantation, likely by infiltrating CD14+ CD163+/− alveolar macrophages. In three LTRs, we demonstrate both rapid D‐strain dissemination and persistence in the transplanted lung for >1 year post‐Tx. Broad inter‐host diversity contrasts with intra‐host genotype sequence stability upon transmission, during follow‐up and across compartments. In D+R+ LTRs, HCMV strains of both, D and R origin can emerge first and dominate long‐term in subsequent episodes of infection, indicating replication of both sources despite pre‐existing immunity.

Deep Amplicon Sequencing Reveals Culture Selection of Mycobacterium Tuberculosis by Clinical Samples.

The commonly-used drug susceptibility testing (DST) relies on bacterial culture and faces shortcomings such as long turnaround time and clone/subclone selection. We developed a targeted deep amplification sequencing (DAS) method directly applied to clinical specimens. In this DAS panel, we examined 941 drug-resistant mutations associated with 20 anti-tuberculosis drugs with an initial amount of 4 pg DNA and reduced clinical testing time from 20 days to two days. A prospective study was conducted using 115 clinical specimens mainly with Xpert® Mycobacterium tuberculosis/rifampicin (Xpert MTB/RIF) assay positive to evaluate drug-resistant mutation detection. DAS was performed on culture-free specimens, while culture-dependent isolates were used for phenotypic DST, DAS, and whole-genome sequencing (WGS). For in silico molecular DST, our result based on DAS panel revealed the similar accuracy to three published reports based on WGS. For 82 isolates, application of DAS showed better sensitivity (93.03% vs. 92.16%), specificity (96.10% vs. 95.02%), and accuracy (91.33% vs. 90.62%) than Mykrobe software using WGS. Compared to culture-dependent WGS, culture-free DAS provides a full picture of sequence variation at population level, exhibiting in detail the gain-and-loss variants caused by bacterial culture. Our study performs a systematic verification of the advantages of DAS in clinical applications and comprehensively illustrates the discrepancy in Mycobacterium tuberculosis before and after culture.

Use of Deep-Amplicon Sequencing (DAS), Real-Time PCR and In Situ Hybridization to Detect H. pylori and Other Pathogenic Helicobacter Species in Feces from Children.

Detecting Helicobacter pylori in fecal samples is easier and more comfortable than invasive techniques, especially in children. Thus, the objective of the present work was to detect H. pylori in feces from children by molecular methods as an alternative for diagnostic and epidemiological studies. Forty-five fecal samples were taken from pediatric patients who presented symptoms compatible with H. pylori infection. HpSA test, culture, real-time quantitative PCR (qPCR), fluorescence in situ hybridization (FISH), direct viable count associated with FISH (DVC-FISH), and Illumina-based deep-amplicon sequencing (DAS) were applied. No H. pylori colonies were isolated from the samples. qPCR analysis detected H. pylori in the feces of 24.4% of the patients. In comparison, DVC-FISH analysis showed the presence of viable H. pylori cells in 53.3% of the samples, 37% of which carried 23S rRNA mutations that confer resistance to clarithromycin. After DAS, H. pylori-specific 16S rDNA sequences were detected in 26 samples. In addition, DNA from H. hepaticus was identified in 10 samples, and H. pullorum DNA was detected in one sample. The results of this study show the presence of H. pylori, H. hepaticus, and H. pullorum in children's stools, demonstrating the coexistence of more than one Helicobacter species in the same patient. The DVC-FISH method showed the presence of viable, potentially infective H. pylori cells in a high percentage of the children's stools. These results support the idea that fecal-oral transmission is probably a common route for H. pylori and suggest possible fecal-oral transmission of other pathogenic Helicobacter species.

Analytic optimization of Plasmodium falciparum marker gene haplotype recovery from amplicon deep sequencing of complex mixtures.

Molecular epidemiologic studies of malaria parasites and other pathogens commonly employ amplicon deep sequencing (AmpSeq) of marker genes derived from dried blood spots (DBS) to answer public health questions related to topics such as transmission and drug resistance. As these methods are increasingly employed to inform direct public health action, it is important to rigorously evaluate the risk of false positive and false negative haplotypes derived from clinically-relevant sample types. We performed a control experiment evaluating haplotype recovery from AmpSeq of 5 marker genes (ama1, csp, msp7, sera2, and trap) from DBS containing mixtures of DNA from 1 to 10 known P. falciparum reference strains across 3 parasite densities in triplicate (n = 270 samples). While false positive haplotypes were present across all parasite densities and mixtures, we optimized censoring criteria to remove 83% (148/179) of false positives while removing only 8% (67/859) of true positives. Post-censoring, the median pairwise Jaccard distance between replicates was 0.83. We failed to recover 35% (477/1365) of haplotypes expected to be present in the sample. Haplotypes were more likely to be missed in low-density samples with <1.5 genomes/μL (OR: 3.88, CI: 1.82-8.27, vs. high-density samples with ≥75 genomes/μL) and in samples with lower read depth (OR per 10,000 reads: 0.61, CI: 0.54-0.69). Furthermore, minority haplotypes within a sample were more likely to be missed than dominant haplotypes (OR per 0.01 increase in proportion: 0.96, CI: 0.96-0.97). Finally, in clinical samples the percent concordance across markers for multiplicity of infection ranged from 40%-80%. Taken together, our observations indicate that, with sufficient read depth, the majority of haplotypes can be successfully recovered from DBS while limiting the false positive rate.

Worldwide absence of canonical benzimidazole resistance-associated mutations within β-tubulin genes from Ascaris.

The giant roundworm Ascaris is an intestinal nematode, causing ascariasis by infecting humans and pigs worldwide. Recent estimates suggest that Ascaris infects over half a billion people, with chronic infections leading to reduced growth and cognitive ability. Ascariasis affects innumerable pigs worldwide and is known to reduce production yields via decreased growth and condemnation of livers. The predominant anthelminthic drugs used to treat ascariasis are the benzimidazoles. Benzimidazoles interact with β-tubulins and block their function, and several benzimidazole resistance-associated mutations have been described in the β-tubulins of ruminant nematodes. Recent research on ascarids has shown that these canonical benzimidazole resistance-associated mutations are likely not present in the β-tubulins of Ascaris, Ascaridia or Parascaris, even in phenotypically resistant populations. To further determine the putative absence of key β-tubulin polymorphisms, we screened two β-tubulin isotypes of Ascaris, highly expressed in adult worms. Using adult and egg samples of Ascaris obtained from pigs and humans worldwide, we performed deep amplicon sequencing to look for canonical resistance-associated mutations in Ascaris β-tubulins. Subsequently, we examined these data in closer detail to study the population dynamics of Ascaris and genetic diversity within the two isotypes and tested whether genotypes appeared to partition across human and pig hosts. In the 187 isolates, 69 genotypes were found, made up of eight haplotypes of β-tubulin isotype A and 20 haplotypes of isotype B. Single nucleotide polymorphisms were seen at 14 and 37 positions for β-tubulin isotype A and isotype B, respectively. No evidence of any canonical benzimidazole resistance-associated mutations was found in either human- or pig-derived Ascaris isolates. There was, however, a difference in the genetic diversity of each isotype and distribution of β-tubulin genotypes between human- and pig-derived Ascaris. Statistical tests of population differentiation show significant differences (p < 0.001) between pig- and human-derived worms; however, more diversity was seen between worms from different populations than worms from different hosts. Our work suggests an absence of canonical β-tubulin mutations within Ascaris, but alternative modes of anthelminthic resistance may emerge necessitating continued genetic scrutiny alongside monitoring of drug efficacy.

A Targeted Deep Sequencing Method to Quantify Endogenous Retrovirus Gag Sequence Variants and Open Reading Frames Expressed in Nonobese Diabetic Mice.

Endogenous retroviruses (ERVs) are involved in autoimmune diseases such as type 1 diabetes (T1D). ERV gene products homologous to murine leukemia retroviruses are expressed in the pancreatic islets of NOD mice, a model of T1D. One ERV gene, Gag, with partial or complete open reading frames (ORFs), is detected in the islets, and it contains many sequence variants. An amplicon deep sequencing analysis was established by targeting a conserved region within the Gag gene to compare NOD with T1D-resistant mice or different ages of prediabetic NOD mice. We observed that the numbers of different Gag variants and ORFs are linked to T1D susceptibility. More importantly, these numbers change during the course of diabetes development and can be quantified to calculate the levels of disease progression. Sequence alignment analysis led to identification of additional markers, including nucleotide mismatching and amino acid consensus at specific positions that can distinguish the early and late stages, before diabetes onset. Therefore, the expression of sequence variants and ORFs of ERV genes, particularly Gag, can be quantified as biomarkers to estimate T1D susceptibility and disease progression.

Validation of deep amplicon sequencing of Dicrocoelium in small ruminants from Northern regions of Pakistan.

Dicrocoelium lancet flukes cause significant production loss in ruminant livestock. Although co-infection with multiple Dicrocoelium species within a host is common, techniques for studying the composition of these complex parasite communities are lacking. The pathogenicity, epidemiology, and therapeutic susceptibility of different helminth species vary, and little is known about the interactions that take place between co-infecting species and their hosts. Here, we describe the first applicationof metabarcoding deep amplicon sequencing method to studythe Dicrocoelium species in sheep and goats. First, rDNA ITS-2 sequences of four Dicrocoelium species (Dicrocoelium dendriticum, Dicrocoelium hospes, Dicrocoelium orientalis, and Dicrocoelium chinensis) were extracted from the NCBI public database. Phylogenetic analysis revealed separate clades of Dicrocoelium species; hence, molecular differentiation between each species is possible in co-infections. Second, 202 flukes belonging to seventeen host populations (morphologically verified as belonging to the Dicrocoelium genus) were evaluated to determine the deep amplicon sequencing read threshold of an individual fluke for each of the four species. The accuracy of the method in proportional quantification of samples collected from single hosts was further assessed. Overall, 198 (98.01%) flukes were confirmed as D. dendriticum and 1.98% produced no reads. The comparison of genetic distances between rDNA ITS-2 revealed 86% to 98% identity between the Dicrocoelium species. Phylogenetic analysis demonstrated a distinct clustering of species, apart from D. orientalis and D. chinensis, which sit very close to each other in a single large clade whereas D. hospes and D. dendriticum are separated into their own clade. In conclusion each sample was identified as D. dendriticum based on the proportion of MiSeq reads and validated the presence of this group of parasites in the Gilgit Baltistan and Khyber Pakhtunkhwa provinces of Pakistan. The metabarcoding deep amplicon sequencing technology and bioinformatics pathway have several potential applications, including species interactions during co-infections, identifying the host and geographical distribution of Dicrocoelium in livestock, drug therapy response evaluation and understanding of the emergence and spread of drug resistance.