SeekDeep: single-base resolution de novo clustering for amplicon deep sequencing

ONT sequencing identifies a high prevalence of crt sensitive, triple mutant dhfr and single mutant dhps parasites within an ANC population in Nigeria

Background

Malaria in pregnancy is a major public health issue, particularly among vulnerable populations in malaria-endemic sub-Saharan African countries. To mitigate its risks, WHO recommends sulphadoxine-pyrimethamine (SP) for chemoprevention and artemisinin-based combination therapy (ACT) to treat uncomplicated Plasmodium falciparum malaria. These interventions have helped to alleviate the risk associated with malaria in pregnancy; however, in the context of the emergence of SP- and ACT-resistant P. falciparum, maintained efficacy is under threat. Molecular surveillance is a reliable tool to monitor the emergence of resistance where molecular markers are known. Thus, the objective of the study was to use a multiplexed amplicon Oxford Nanopore sequencing approach to assess the molecular markers for antimalarial resistance among pregnant women in Nigeria.

Methods

Dried blood spots (DBS) were collected from pregnant women who received IPTp-SP at the enrollment and follow-up visits. P. falciparum genomic DNA was extracted by the Chelex® method and Pf18S qPCR was used to detect parasite DNA in each sample. With nested PCR assays, fragments of Pfdhps, Pfdhfr, Pfmdr1, Pfcrt, Pfk13 and Pfama1 genes were amplified and multiplexed amplicon-based sequencing was conducted on the minION Oxford Nanopore Technology.

Result

In total, 251 pregnant women were enrolled in the study and 457 DBS samples were collected. P. falciparum genomic DNA was detected in 12% (56/457) of the samples, 31 at baseline and the remaining during the follow-up visits. Pfama1, pfk13, Pfdhps, Pfdhfr, Pfmdr1 and Pfcrt were successfully sequenced in a single run. Notably, k13 artemisinin resistance mutations were absent, the frequencies of Pfdhfr and Pfdhps SP resistance haplotypes, IRN for pyrimethamine resistance and ISGKA/IAGKA associated with sulphadoxine resistance were 82% (36/44) and 64% (27/42), respectively, and the Pfcrt CVIET resistant haplotype was at approximately 22% (7/32).

Conclusion and recommendations

Here a multiplexed amplicon-based ONT assay established that triple mutant Pfdfhr-IRN, double mutant Pfdhps-SG haplotypes and the chloroquine sensitive strain were prevalent among pregnant women in Nigeria.

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

Background

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.

Methods

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.

Results

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.

Conclusion

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.

A retrospective analysis of P. falciparum drug resistance markers detects an early (2016/17) high prevalence of the k13 C469Y mutation in asymptomatic infections in Northern Uganda

The emergence of drug-resistant Plasmodium falciparum parasites in sub-Saharan Africa will substantially challenge malaria control. Here, we evaluated the frequency of common drug resistance markers among adolescents from Northern Uganda with asymptomatic infections. We used an established amplicon deep sequencing strategy to screen dried blood spot samples collected from 2016 to 2017 during a reported malaria epidemic within the districts of Kitgum and Pader in Northern Uganda. We screened single-nucleotide polymorphisms within: kelch13 (Pfk13), dihydropteroate synthase (Pfdhps), multidrug resistance-1 (Pfmdr1), dihydrofolate reductase (Pfdhfr), and apical membrane antigen (Pfama1) genes. Within the study population, the median age was 15 years (14.3-15.0, 95% CI), and 54.9% (78/142) were Plasmodium positive by 18S rRNA qPCR, which were subsequently targeted for sequencing analysis. We observed a high frequency of resistance markers particularly for sulfadoxine-pyrimethamine (SP), with no wild-type-only parasites observed for Pfdhfr (N51I, C59R, and S108N) and Pfdhps (A437G and K540E) mutations. Within Pfmdr1, mixed infections were common for NF/NY (98.5%). While for artemisinin resistance, in kelch13, there was a high frequency of C469Y (34%). Using the pattern for Pfama1, we found a high level of polygenomic infections with all individuals presenting with complexity of infection greater than 2 with a median of 6.9. The high frequency of the quintuple SP drug-resistant parasites and the C469Y artemisinin resistance-associated mutation in asymptomatic individuals suggests an earlier high prevalence than previously reported from symptomatic malaria surveillance studies (in 2016/2017). Our data demonstrate the urgency for routine genomic surveillance programs throughout Africa and the value of deep sequencing.

Sensitive and modular amplicon sequencing of Plasmodium falciparum diversity and resistance for research and public health

Targeted amplicon sequencing is a powerful and efficient tool to interrogate the P. falciparum genome and generate actionable data from infections to complement traditional malaria epidemiology. For maximum impact, genomic tools should be multi-purpose, robust, sensitive and reproducible. We developed, characterized, and implemented MAD4HatTeR, an amplicon sequencing panel based on Multiplex Amplicons for Drug, Diagnostic, Diversity, and Differentiation Haplotypes using Targeted Resequencing, along with a bioinformatic pipeline for data analysis. MAD4HatTeR targets 165 highly diverse loci, focusing on multiallelic microhaplotypes; key markers for drug and diagnostic resistance, including duplications and deletions; and csp and potential vaccine targets. In addition, it can detect non-falciparum Plasmodium species. We used laboratory control and field sample data to demonstrate the high sensitivity and robustness of the panel. The successful implementation of this method in five laboratories, including three in malaria-endemic African countries, showcases its feasibility in generating reproducible data across laboratories. Finally, we introduce an analytical approach to detect gene duplications and deletions from amplicon sequencing data. MAD4HatTeR is thus a powerful research tool and a robust resource for malaria public health surveillance and control.

MIPs: multi-locus intron polymorphisms in species identification and population genomics

The study of species groups in which the presence of interspecific hybridization or introgression phenomena is known or suspected involves analysing shared bi-parentally inherited molecular markers. Current methods are based on different categories of markers among which the classical microsatellites or the more recent genome wide approaches for the analyses of thousands of SNPs or hundreds of microhaplotypes through high throughput sequencing. Our approach utilizes intron-targeted amplicon sequencing to characterise multi-locus intron polymorphisms (MIPs) and assess genetic diversity. These highly variable intron regions, combined with inter-specific transferable loci, serve as powerful multiple-SNP markers potentially suitable for various applications, from species and hybrid identification to population comparisons, without prior species knowledge. We developed the first panel of MIPs highly transferable across fish genomes, effectively distinguishing between species, even those closely related, and populations with different structures. MIPs offer versatile, hypervariable nuclear markers and promise to be especially useful when multiple nuclear loci must be genotyped across different species, such as for the monitoring of interspecific hybridization. Moreover, the relatively long sequences obtained ease the development of single-locus PCR-based diagnostic markers. This method, here demonstrated in teleost fishes, can be readily applied to other taxa, unlocking a new source of genetic variation.

Synergism between remdesivir and ribavirin leads to SARS-CoV-2 extinction in cell culture

BACKGROUND AND PURPOSE

There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis.

EXPERIMENTAL APPROACH

The effect of combinations of remdesivir and ribavirin on the infectivity of SARS-CoV-2 in cell culture has been tested. Viral populations were monitored by ultra-deep sequencing, and the decrease of infectivity as a result of the treatment was measured.

KEY RESULTS

Remdesivir and ribavirin exerted a synergistic inhibitory activity against SARS-CoV-2, quantified both by CompuSyn (Chou-Talalay method) and Synergy Finder (ZIP-score model). In serial passage experiments, virus extinction was readily achieved with remdesivir-ribavirin combinations at concentrations well below their cytotoxic 50 value, but not with the drugs used individually. Deep sequencing of treated viral populations showed that remdesivir, ribavirin, and their combinations evoked significant increases of the number of viral mutations and haplotypes, as well as modification of diversity indices that characterize viral quasi-species.

CONCLUSION AND IMPLICATIONS

SARS-CoV-2 extinction can be achieved by synergistic combination treatments based on lethal mutagenesis. In addition, the results offer prospects of triple drug treatments for effective SARS-CoV-2 suppression.

High-throughput system for the thermostability analysis of proteins

Thermal stability of proteins is a primary metric for evaluating their physical properties. Although researchers attempted to predict it using machine learning frameworks, their performance has been dependent on the quality and quantity of published data. This is due to the technical limitation that thermodynamic characterization of protein denaturation by fluorescence or calorimetry in a high-throughput manner has been challenging. Obtaining a melting curve that derives solely from the target protein requires laborious purification, making it far from practical to prepare a hundred or more samples in a single workflow. Here, we aimed to overcome this throughput limitation by leveraging the high protein secretion efficacy of Brevibacillus and consecutive treatment with plate-scale purification methodologies. By handling the entire process of expression, purification, and analysis on a per-plate basis, we enabled the direct observation of protein denaturation in 384 samples within 4 days. To demonstrate a practical application of the system, we conducted a comprehensive analysis of 186 single mutants of a single-chain variable fragment of nivolumab, harvesting the melting temperature (Tm) ranging from -9.3 up to +10.8°C compared to the wild-type sequence. Our findings will allow for data-driven stabilization in protein design and streamlining the rational approaches.

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.

Effectiveness of an Unsupervised Primaquine Regimen for Preventing Plasmodium vivax Malaria Relapses in Northeast Myanmar: A Single-Arm Nonrandomized Observational Study

BACKGROUND

Plasmodium vivax presents a significant challenge for malaria elimination in the Greater Mekong Subregion. We evaluated the effectiveness of primaquine for reducing relapses of vivax malaria.

METHODS

Patients with uncomplicated P vivax malaria from eastern Myanmar received chloroquine (25-mg base/kg given in 3 days) plus unsupervised PQ (0.25 mg/kg/d for 14 days) without screening for glucose-6-phosphate dehydrogenase deficiency and were followed for a year.

RESULTS

A total of 556 patients were enrolled to receive the chloroquine/primaquine treatment from February 2012 to August 2013. During the follow-up, 38 recurrences were detected, presenting a cumulative recurrence rate of 9.1% (95% CI, 4.1%-14.1%). Genotyping at the pvmsp1 and pvmsp3α loci by amplicon deep sequencing and model prediction indicated that 13 of the 27 recurrences with genotyping data were likely due to relapses. Notably, all confirmed relapses occurred within the first 6 months.

CONCLUSIONS

The unsupervised standard dose of primaquine was highly effective as a radical cure for P vivax malaria in eastern Myanmar. The high presumed effectiveness might have benefited from the health messages delivered during the enrollment and follow-up activities. Six-month follow-ups in the Greater Mekong Subregion are sufficient for detecting most relapses.

Molecular tools are crucial for malaria elimination

The eradication of Plasmodium parasites, responsible for malaria, is a daunting global public health task. It requires a comprehensive approach that addresses symptomatic, asymptomatic, and submicroscopic cases. Overcoming this challenge relies on harnessing the power of molecular diagnostic tools, as traditional methods like microscopy and rapid diagnostic tests fall short in detecting low parasitaemia, contributing to the persistence of malaria transmission. By precisely identifying patients of all types and effectively characterizing malaria parasites, molecular tools may emerge as indispensable allies in the pursuit of malaria elimination. Furthermore, molecular tools can also provide valuable insights into parasite diversity, drug resistance patterns, and transmission dynamics, aiding in the implementation of targeted interventions and surveillance strategies. In this review, we explore the significance of molecular tools in the pursuit of malaria elimination, shedding light on their key contributions and potential impact on public health.

Genotyping Plasmodium falciparum gametocytes using amplicon deep sequencing

BACKGROUND

Understanding the dynamics of gametocyte production in polyclonal Plasmodium falciparum infections requires a genotyping method that detects distinct gametocyte clones and estimates their relative frequencies. Here, a marker was identified and evaluated to genotype P. falciparum mature gametocytes using amplicon deep sequencing.

METHODS

A data set of polymorphic regions of the P. falciparum genome was mined to identify a gametocyte genotyping marker. To assess marker resolution, the number of unique haplotypes in the marker region was estimated from 95 Malawian P. falciparum whole genome sequences. Specificity of the marker for detection of mature gametocytes was evaluated using reverse transcription-polymerase chain reaction of RNA extracted from NF54 mature gametocytes and rings from a non-gametocyte-producing strain of P. falciparum. Amplicon deep sequencing was performed on experimental mixtures of mature gametocytes from two distinct parasite clones, as well as gametocyte-positive P. falciparum field isolates to evaluate the quantitative ability and determine the limit of detection of the genotyping approach.

RESULTS

A 400 bp region of the pfs230 gene was identified as a gametocyte genotyping marker. A larger number of unique haplotypes was observed at the pfs230 marker (34) compared to the sera-2 (18) and ama-1 (14) markers in field isolates from Malawi. RNA and DNA genotyping accurately estimated gametocyte and total parasite clone frequencies when evaluating agreement between expected and observed haplotype frequencies in gametocyte mixtures, with concordance correlation coefficients of 0.97 [95% CI: 0.92-0.99] and 0.92 [95% CI: 0.83-0.97], respectively. The detection limit of the genotyping method for male gametocytes was 0.41 pfmget transcripts/µl [95% CI: 0.28-0.72] and for female gametocytes was 1.98 ccp4 transcripts/µl [95% CI: 1.35-3.68].

CONCLUSIONS

A region of the pfs230 gene was identified as a marker to genotype P. falciparum gametocytes. Amplicon deep sequencing of this marker can be used to estimate the number and relative frequency of parasite clones among mature gametocytes within P. falciparum infections. This gametocyte genotyping marker will be an important tool for studies aimed at understanding dynamics of gametocyte production in polyclonal P. falciparum infections.

Exploiting long read sequencing to detect azole fungicide resistance mutations in Pyrenophora teres using unique molecular identifiers

Resistance to fungicides is a global challenge as target proteins under selection can evolve rapidly, reducing fungicide efficacy. To manage resistance, detection technologies must be fast and flexible enough to cope with a rapidly increasing number of mutations. The most important agricultural fungicides are azoles that target the ergosterol biosynthetic enzyme sterol 14α-demethylase (CYP51). Mutations associated with azole resistance in the Cyp51 promoter and coding sequence can co-occur in the same allele at different positions and codons, increasing the complexity of resistance detection. Resistance mutations arise rapidly and cannot be detected using traditional amplification-based methods if they are not known. To capture the complexity of azole resistance in two net blotch pathogens of barley we used the Oxford Nanopore MinION to sequence the promoter and coding sequence of Cyp51A. This approach detected all currently known mutations from biologically complex samples increasing the simplicity of resistance detection as multiple alleles can be profiled in a single assay. With the mobility and decreasing cost of long read sequencing, we demonstrate this approach is broadly applicable for characterizing resistance within known agrochemical target sites.

SARS-CoV-2 mutant spectra as variant of concern nurseries: endless variation?

Introduction

SARS-CoV-2 isolates of a given clade may contain low frequency genomes that encode amino acids or deletions which are typical of a different clade.

Methods

Here we use high resolution ultra-deep sequencing to analyze SARS-CoV-2 mutant spectra.

Results

In 6 out of 11 SARS-CoV-2 isolates from COVID-19 patients, the mutant spectrum of the spike (S)-coding region included two or more amino acids or deletions, that correspond to discordant viral clades. A similar observation is reported for laboratory populations of SARS-CoV-2 USA-WA1/2020, following a cell culture infection in the presence of remdesivir, ribavirin or their combinations. Moreover, some of the clade-discordant genome residues are found in the same haplotype within an amplicon.

Discussion

We evaluate possible interpretations of these findings, and reviewed precedents for rapid selection of genomes with multiple mutations in RNA viruses. These considerations suggest that intra-host evolution may be sufficient to generate minority sequences which are closely related to sequences typical of other clades. The results provide a model for the origin of variants of concern during epidemic spread─in particular Omicron lineages─that does not require prolonged infection, involvement of immunocompromised individuals, or participation of intermediate, non-human hosts.

Targeted amplicon deep sequencing of ama1 and mdr1 to track within-host P. falciparum diversity throughout treatment in a clinical drug trial

Introduction

Antimalarial therapeutic efficacy studies are routinely conducted in malaria-endemic countries to assess the effectiveness of antimalarial treatment strategies. Targeted amplicon sequencing (AmpSeq) uniquely identifies and quantifies genetically distinct parasites within an infection. In this study, AmpSeq of Plasmodium falciparum apical membrane antigen 1 ( ama1), and multidrug resistance gene 1 ( mdr1), were used to characterise the complexity of infection (COI) and drug-resistance genotypes, respectively.

Methods

P. falciparum-positive samples were obtained from a triple artemisinin combination therapy clinical trial conducted in 30 children under 13 years of age between 2018 and 2019 in Kilifi, Kenya. Nine of the 30 participants presented with recurrent parasitemia from day 26 (624h) onwards. The ama1 and mdr1 genes were amplified and sequenced, while msp1, msp2 and glurp data were obtained from the original clinical study.

Results

The COI was comparable between ama1 and msp1, msp2 and glurp; overall, ama1 detected more microhaplotypes. Based on ama1, a stable number of microhaplotypes were detected throughout treatment until day 3. Additionally, a recrudescent infection was identified with an ama1 microhaplotype initially observed at 30h and later in an unscheduled follow-up visit. Using the relative frequencies of ama1 microhaplotypes and parasitemia, we identified a fast (<1h) and slow (>5h) clearing microhaplotype. As expected, only two mdr1 microhaplotypes (NF and NY) were identified based on the combination of amino acid polymorphisms at codons 86 and 184.

Conclusions

This study highlights AmpSeq as a tool for highly-resolution tracking of parasite microhaplotypes throughout treatment and can detect variation in microhaplotype clearance estimates. AmpSeq can also identify slow-clearing microhaplotypes, a potential early sign of selection during treatment. Consequently, AmpSeq has the capability of improving the discriminatory power to distinguish recrudescences from reinfections accurately.

Incipient functional SARS-CoV-2 diversification identified through neural network haplotype maps

Since its introduction in the human population, SARS-CoV-2 has evolved into multiple clades, but the events in its intrahost diversification are not well understood. Here, we compare three-dimensional (3D) self-organized neural haplotype maps (SOMs) of SARS-CoV-2 from thirty individual nasopharyngeal diagnostic samples obtained within a 19-day interval in Madrid (Spain), at the time of transition between clades 19 and 20. SOMs have been trained with the haplotype repertoire present in the mutant spectra of the nsp12- and spike (S)-coding regions. Each SOM consisted of a dominant neuron (displaying the maximum frequency), surrounded by a low-frequency neuron cloud. The sequence of the master (dominant) neuron was either identical to that of the reference Wuhan-Hu-1 genome or differed from it at one nucleotide position. Six different deviant haplotype sequences were identified among the master neurons. Some of the substitutions in the neural clouds affected critical sites of the nsp12-nsp8-nsp7 polymerase complex and resulted in altered kinetics of RNA synthesis in an in vitro primer extension assay. Thus, the analysis has identified mutations that are relevant to modification of viral RNA synthesis, present in the mutant clouds of SARS-CoV-2 quasispecies. These mutations most likely occurred during intrahost diversification in several COVID-19 patients, during an initial stage of the pandemic, and within a brief time period.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

Over recent years, progress in molecular markers for genotyping malaria parasites has enabled informative studies of epidemiology and transmission dynamics. Results have highlighted the value of these tools for surveillance to support malaria control and elimination strategies. There are many different types and panels of markers available for malaria parasite genotyping, and for end users, the nuances of these markers with respect to 'use case', resolution, and accuracy, are not well defined. This review clarifies issues surrounding different molecular markers and their application to malaria control and elimination. We describe available marker panels, use cases, implications for different transmission settings, limitations, access, cost, and data accuracy. The information provided can be used as a guide for molecular epidemiology and surveillance of malaria.

Replicative fitness and pathogenicity of primate lentiviruses in lymphoid tissue, primary human and chimpanzee cells: relation to possible jumps to humans

BACKGROUND

Simian immunodeficiency viruses (SIV) have been jumping between non-human primates in West/Central Africa for thousands of years and yet, the HIV-1 epidemic only originated from a primate lentivirus over 100 years ago.

METHODS

This study examined the replicative fitness, transmission, restriction, and cytopathogenicity of 22 primate lentiviruses in primary human lymphoid tissue and both primary human and chimpanzee peripheral blood mononuclear cells.

FINDINGS

Pairwise competitions revealed that SIV from chimpanzees (cpz) had the highest replicative fitness in human or chimpanzee peripheral blood mononuclear cells, even higher fitness than HIV-1 group M strains responsible for worldwide epidemic. The SIV strains belonging to the "HIV-2 lineage" (including SIVsmm, SIVmac, SIVagm) had the lowest replicative fitness. SIVcpz strains were less inhibited by human restriction factors than the "HIV-2 lineage" strains. SIVcpz efficiently replicated in human tonsillar tissue but did not deplete CD4+ T-cells, consistent with the slow or nonpathogenic disease observed in most chimpanzees. In contrast, HIV-1 isolates and SIV of the HIV-2 lineage were pathogenic to the human tonsillar tissue, almost independent of the level of virus replication.

INTERPRETATION

Of all primate lentiviruses, SIV from chimpanzees appears most capable of infecting and replicating in humans, establishing HIV-1. SIV from other Old World monkeys, e.g. the progenitor of HIV-2, replicate slowly in humans due in part to restriction factors. Nonetheless, many of these SIV strains were more pathogenic than SIVcpz. Either SIVcpz evolved into a more pathogenic virus while in humans or a rare SIVcpz, possibly extinct in chimpanzees, was pathogenic immediately following the jump into human.

FUNDING

Support for this study to E.J.A. was provided by the NIH/NIAID R01 AI49170 and CIHR project grant 385787. Infrastructure support was provided by the NIH CFAR AI36219 and Canadian CFI/Ontario ORF 36287. Efforts of J.A.B. and N.J.H. was provided by NIH AI099473 and for D.H.C., by VA and NIH AI AI080313.

CHOIR improves significance-based detection of cell types and states from single-cell data

Clustering is a critical step in the analysis of single-cell data, as it enables the discovery and characterization of putative cell types and states. However, most popular clustering tools do not subject clustering results to statistical inference testing, leading to risks of overclustering or underclustering data and often resulting in ineffective identification of cell types with widely differing prevalence. To address these challenges, we present CHOIR (clustering hierarchy optimization by iterative random forests), which applies a framework of random forest classifiers and permutation tests across a hierarchical clustering tree to statistically determine which clusters represent distinct populations. We demonstrate the enhanced performance of CHOIR through extensive benchmarking against 14 existing clustering methods across 100 simulated and 4 real single-cell RNA-seq, ATAC-seq, spatial transcriptomic, and multi-omic datasets. CHOIR can be applied to any single-cell data type and provides a flexible, scalable, and robust solution to the important challenge of identifying biologically relevant cell groupings within heterogeneous single-cell data.

Plasmodium falciparum population structure inferred by msp1 amplicon sequencing of parasites collected from febrile patients in Kenya

BACKGROUND

Multiplicity of infection (MOI) is an important measure of Plasmodium falciparum diversity, usually derived from the highly polymorphic genes, such as msp1, msp2 and glurp as well as microsatellites. Conventional methods of deriving MOI lack fine resolution needed to discriminate minor clones. This study used amplicon sequencing (AmpliSeq) of P. falciparum msp1 (Pfmsp1) to measure spatial and temporal genetic diversity of P. falciparum.

METHODS

264 P. falciparum positive blood samples collected from areas of differing malaria endemicities between 2010 and 2019 were used. Pfmsp1 gene was amplified and amplicon libraries sequenced on Illumina MiSeq. Sequences were aligned against a reference sequence (NC_004330.2) and clustered to detect fragment length polymorphism and amino acid variations.

RESULTS

Children < 5 years had higher parasitaemia (median = 23.5 ± 5 SD, p = 0.03) than the > 5-14 (= 25.3 ± 5 SD), and those > 15 (= 25.1 ± 6 SD). Of the alleles detected, 553 (54.5%) were K1, 250 (24.7%) MAD20 and 211 (20.8%) RO33 that grouped into 19 K1 allelic families (108-270 bp), 14 MAD20 (108-216 bp) and one RO33 (153 bp). AmpliSeq revealed nucleotide polymorphisms in alleles that had similar sizes, thus increasing the K1 to 104, 58 for MAD20 and 14 for RO33. By AmpliSeq, the mean MOI was 4.8 (± 0.78, 95% CI) for the malaria endemic Lake Victoria region, 4.4 (± 1.03, 95% CI) for the epidemic prone Kisii Highland and 3.4 (± 0.62, 95% CI) for the seasonal malaria Semi-Arid region. MOI decreased with age: 4.5 (± 0.76, 95% CI) for children < 5 years, compared to 3.9 (± 0.70, 95% CI) for ages 5 to 14 and 2.7 (± 0.90, 95% CI) for those > 15. Females' MOI (4.2 ± 0.66, 95% CI) was not different from males 4.0 (± 0.61, 95% CI). In all regions, the number of alleles were high in the 2014-2015 period, more so in the Lake Victoria and the seasonal transmission arid regions.

CONCLUSION

These findings highlight the added advantages of AmpliSeq in haplotype discrimination and the associated improvement in unravelling complexity of P. falciparum population structure.

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

Molecular epidemiologic studies of malaria parasites 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, haplotypes can be successfully recovered from DBS while limiting the false positive rate.

Distinct life history strategies underpin clear patterns of succession in microparasite communities infecting a wild mammalian host

Individual animals in natural populations tend to host diverse parasite species concurrently over their lifetimes. In free-living ecological communities, organismal life histories shape interactions with their environment, which ultimately forms the basis of ecological succession. However, the structure and dynamics of mammalian parasite communities have not been contextualized in terms of primary ecological succession, in part because few datasets track occupancy and abundance of multiple parasites in wild hosts starting at birth. Here, we studied community dynamics of 12 subtypes of protozoan microparasites (Theileria spp.) in a herd of African buffalo. We show that Theileria communities followed predictable patterns of succession underpinned by four different parasite life history strategies. However, in contrast to many free-living communities, network complexity decreased with host age. Examining parasite communities through the lens of succession may better inform the effect of complex within host eco-evolutionary dynamics on infection outcomes, including parasite co-existence through the lifetime of the host.

Neural stem and progenitor cells support and protect adult hippocampal function via vascular endothelial growth factor secretion

Adult neural stem and progenitor cells (NSPCs) reside in the dentate gyrus (DG) of the hippocampus throughout the lifespan of most mammalian species. In addition to generating new neurons, NSPCs may alter their niche via secretion of growth factors and cytokines. We recently showed that adult DG NSPCs secrete vascular endothelial growth factor (VEGF), which is critical for maintaining adult neurogenesis. Here, we asked whether NSPC-derived VEGF alters hippocampal function independent of adult neurogenesis. We found that loss of NSPC-derived VEGF acutely impaired hippocampal memory, caused neuronal hyperexcitability and exacerbated excitotoxic injury. We also found that NSPCs generate substantial proportions of total DG VEGF and VEGF disperses broadly throughout the DG, both of which help explain how this anatomically-restricted cell population could modulate function broadly. These findings suggest that NSPCs actively support and protect DG function via secreted VEGF, thereby providing a non-neurogenic functional dimension to endogenous NSPCs.

Examining the Early Distribution of the Artemisinin-Resistant Plasmodium falciparum kelch13 R561H Mutation in Areas of Higher Transmission in Rwanda

Background

Artemisinin resistance mutations in Plasmodium falciparum kelch13 (Pfk13) have begun to emerge in Africa, with Pfk13-R561H being the first reported in Rwanda in 2014, but limited sampling left questions about its early distribution and origin.

Methods

We genotyped P. falciparum positive dried blood spot (DBS) samples from a nationally representative 2014-2015 Rwanda Demographic Health Surveys (DHS) HIV study. DBS were subsampled from DHS sampling clusters with >15% P. falciparum prevalence, as determined by rapid testing or microscopy done during the DHS study (n clusters = 67, n samples = 1873).

Results

We detected 476 parasitemias among 1873 residual blood spots from a 2014-2015 Rwanda Demographic Health Survey. We sequenced 351 samples: 341/351 were wild-type (97.03% weighted), and 4 samples (1.34% weighted) harbored R561H that were significantly spatially clustered. Other nonsynonymous mutations found were V555A (3), C532W (1), and G533A (1).

Conclusions

Our study better defines the early distribution of R561H in Rwanda. Previous studies only observed the mutation in Masaka as of 2014, but our study indicates its presence in higher-transmission regions in the southeast of the country at that time.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Atypical Mutational Spectrum of SARS-CoV-2 Replicating in the Presence of Ribavirin

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

Accurate estimation of molecular counts from amplicon sequence data with unique molecular identifiers

MOTIVATION

Amplicon sequencing is widely applied to explore heterogeneity and rare variants in genetic populations. Resolving true biological variants and quantifying their abundance is crucial for downstream analyses, but measured abundances are distorted by stochasticity and bias in amplification, plus errors during polymerase chain reaction (PCR) and sequencing. One solution attaches unique molecular identifiers (UMIs) to sample sequences before amplification. Counting UMIs instead of sequences provides unbiased estimates of abundance. While modern methods improve over naïve counting by UMI identity, most do not account for UMI reuse or collision, and they do not adequately model PCR and sequencing errors in the UMIs and sample sequences.

RESULTS

We introduce Deduplication and Abundance estimation with UMIs (DAUMI), a probabilistic framework to detect true biological amplicon sequences and accurately estimate their deduplicated abundance. DAUMI recognizes UMI collision, even on highly similar sequences, and detects and corrects most PCR and sequencing errors in the UMI and sampled sequences. DAUMI performs better on simulated and real data compared to other UMI-aware clustering methods.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/DormanLab/AmpliCI.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Irrigation-Induced Environmental Changes Sustain Malaria Transmission and Compromise Intervention Effectiveness

BACKGROUND

Irrigated agriculture enhances food security, but it potentially promotes mosquito-borne disease transmission and affects vector intervention effectiveness. This study was conducted in the irrigated and nonirrigated areas of rural Homa Bay and Kisumu Counties, Kenya.

METHODS

We performed cross-sectional and longitudinal surveys to determine Plasmodium infection prevalence, clinical malaria incidence, molecular force of infection (molFOI), and multiplicity of infection. We examined the impact of irrigation on the effectiveness of the new interventions.

RESULTS

We found that irrigation was associated with >2-fold higher Plasmodium infection prevalence and 3-fold higher clinical malaria incidence compared to the nonirrigated area. Residents in the irrigated area experienced persistent, low-density parasite infections and higher molFOI. Addition of indoor residual spraying was effective in reducing malaria burden, but the reduction was more pronounced in the nonirrigated area than in the irrigated area.

CONCLUSIONS

Our findings collectively suggest that irrigation may sustain and enhance Plasmodium transmission and affects intervention effectiveness.

Classification of disease recurrence using transition likelihoods with expectation-maximization algorithm

When an infectious disease recurs, it may be due to treatment failure or a new infection. Being able to distinguish and classify these two different outcomes is critical in effective disease control. A multi-state model based on Markov processes is a typical approach to estimating the transition probability between the disease states. However, it can perform poorly when the disease state is unknown. This article aims to demonstrate that the transition likelihoods of baseline covariates can distinguish one cause from another with high accuracy in infectious diseases such as malaria. A more general model for disease progression can be constructed to allow for additional disease outcomes. We start from a multinomial logit model to estimate the disease transition probabilities and then utilize the baseline covariate's transition information to provide a more accurate classification result. We apply the expectation-maximization (EM) algorithm to estimate unknown parameters, including the marginal probabilities of disease outcomes. A simulation study comparing our classifier to the existing two-stage method shows that our classifier has better accuracy, especially when the sample size is small. The proposed method is applied to determining relapse vs reinfection outcomes in two Plasmodium vivax treatment studies from Cambodia that used different genotyping approaches to demonstrate its practical use.

High-resolution species assignment of Anopheles mosquitoes using k-mer distances on targeted sequences

The ANOSPP amplicon panel is a genus-wide targeted sequencing panel to facilitate large-scale monitoring of Anopheles species diversity. Combining information from the 62 nuclear amplicons present in the ANOSPP panel allows for a more senstive and specific species assignment than single gene (e.g. COI) barcoding, which is desirable in the light of permeable species boundaries. Here, we present NNoVAE, a method using Nearest Neighbours (NN) and Variational Autoencoders (VAE), which we apply to k-mers resulting from the ANOSPP amplicon sequences in order to hierarchically assign species identity. The NN step assigns a sample to a species-group by comparing the k-mers arising from each haplotype's amplicon sequence to a reference database. The VAE step is required to distinguish between closely related species, and also has sufficient resolution to reveal population structure within species. In tests on independent samples with over 80% amplicon coverage, NNoVAE correctly classifies to species level 98% of samples within the An. gambiae complex and 89% of samples outside the complex. We apply NNoVAE to over two thousand new samples from Burkina Faso and Gabon, identifying unexpected species in Gabon. NNoVAE presents an approach that may be of value to other targeted sequencing panels, and is a method that will be used to survey Anopheles species diversity and Plasmodium transmission patterns through space and time on a large scale, with plans to analyse half a million mosquitoes in the next five years.

Hypnozoite depletion in successive Plasmodium vivax relapses

Genotyping Plasmodium vivax relapses can provide insights into hypnozoite biology. We performed targeted amplicon sequencing of 127 relapses occurring in Indonesian soldiers returning to malaria-free Java after yearlong deployment in malarious Eastern Indonesia. Hepatic carriage of multiple hypnozoite clones was evident in three-quarters of soldiers with two successive relapses, yet the majority of relapse episodes only displayed one clonal population. The number of clones detected in relapse episodes decreased over time and through successive relapses, especially in individuals who received hypnozoiticidal therapy. Interrogating the multiplicity of infection in this P. vivax relapse cohort reveals evidence of independent activation and slow depletion of hypnozoites over many months by multiple possible mechanisms, including parasite senescence and host immunity.

Investigating relapse patterns in infections of Plasmodium vivax, a parasite that causes malaria, is challenging due to concurrent reinfection events alongside true relapses in most clinical cohorts. We performed sequencing on P. vivax samples from a cohort of Indonesian soldiers who were exposed to the parasite while deployed in a malaria-endemic region and then experienced relapses after their return to a region with no current malaria transmission. From these true relapses, we show that most infected individuals harbor multiple lineages of hypnozoites (latent liver stage parasites that reactivate to cause relapse) but individual relapses are largely driven by a single hypnozoite lineage or clone. Additionally, the average number of parasite clones detected in each relapse decreases over time. These findings suggest that P. vivax hypnozoites activate independently from each other and that their population in the liver decreases over time after the initial infection, possibly due to immune clearance or loss of parasite viability.

Characterization of triatomine bloodmeal sources using direct Sanger sequencing and amplicon deep sequencing methods

Knowledge of host associations of blood-feeding vectors may afford insights into managing disease systems and protecting public health. However, the ability of methods to distinguish bloodmeal sources varies widely. We used two methods—Sanger sequencing and amplicon deep sequencing—to target a 228 bp region of the vertebrate Cytochrome b gene and determine hosts fed upon by triatomines (n = 115) collected primarily in Texas, USA. Direct Sanger sequencing of PCR amplicons was successful for 36 samples (31%). Sanger sequencing revealed 15 distinct host species, which included humans, domestic animals (Canis lupus familiaris, Ovis aries, Gallus gallus, Bos taurus, Felis catus, and Capra hircus), wildlife (Rattus rattus, Incilius nebulifer, Sciurus carolinensis, Sciurus niger, and Odocoileus virginianus), and captive animals (Panthera tigris, Colobus spp., and Chelonoidis carbonaria). Samples sequenced by the Sanger method were also subjected to Illumina MiSeq amplicon deep sequencing. The amplicon deep sequencing results (average of 302,080 usable reads per sample) replicated the host community revealed using Sanger sequencing, and detected additional hosts in five triatomines (13.9%), including two additional blood sources (Procyon lotor and Bassariscus astutus). Up to four bloodmeal sources were detected in a single triatomine (I. nebulifer, Homo sapiens, C. lupus familiaris, and S. carolinensis). Enhanced understanding of vector-host-parasite networks may allow for integrated vector management programs focusing on highly-utilized and highly-infected host species.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Low Complexity of Infection Is Associated With Molecular Persistence of Plasmodium falciparum in Kenya and Tanzania

Background

Plasmodium falciparum resistance to artemisinin-based combination therapies (ACTs) is a threat to malaria elimination. ACT-resistance in Asia raises concerns for emergence of resistance in Africa. While most data show high efficacy of ACT regimens in Africa, there have been reports describing declining efficacy, as measured by both clinical failure and prolonged parasite clearance times.

Methods

Three hundred children aged 2-10 years with uncomplicated P. falciparum infection were enrolled in Kenya and Tanzania after receiving treatment with artemether-lumefantrine. Blood samples were taken at 0, 24, 48, and 72 h, and weekly thereafter until 28 days post-treatment. Parasite and host genetics were assessed, as well as clinical, behavioral, and environmental characteristics, and host anti-malarial serologic response.

Results

While there was a broad range of clearance rates at both sites, 85% and 96% of Kenyan and Tanzanian samples, respectively, were qPCR-positive but microscopy-negative at 72 h post-treatment. A greater complexity of infection (COI) was negatively associated with qPCR-detectable parasitemia at 72 h (OR: 0.70, 95% CI: 0.53-0.94), and a greater baseline parasitemia was marginally associated with qPCR-detectable parasitemia (1,000 parasites/uL change, OR: 1.02, 95% CI: 1.01-1.03). Demographic, serological, and host genotyping characteristics showed no association with qPCR-detectable parasitemia at 72 h. Parasite haplotype-specific clearance slopes were grouped around the mean with no association detected between specific haplotypes and slower clearance rates.

Conclusions

Identifying risk factors for slow clearing P. falciparum infections, such as COI, are essential for ongoing surveillance of ACT treatment failure in Kenya, Tanzania, and more broadly in sub-Saharan Africa.

Vaccine breakthrough infections with SARS-CoV-2 Alpha mirror mutations in Delta Plus, Iota, and Omicron

Replication of SARS-CoV-2 in the human population is defined by distributions of mutants that are present at different frequencies within the infected host and can be detected by ultra-deep sequencing techniques. In this study, we examined the SARS-CoV-2 mutant spectra of amplicons from the spike-coding (S-coding) region of 5 nasopharyngeal isolates derived from patients with vaccine breakthrough. Interestingly, all patients became infected with the Alpha variant, but amino acid substitutions that correspond to the Delta Plus, Iota, and Omicron variants were present in the mutant spectra of the resident virus. Deep sequencing analysis of SARS-CoV-2 from patients with vaccine breakthrough revealed a rich reservoir of mutant types and may also identify tolerated substitutions that can be represented in epidemiologically dominant variants.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Amplicon sequencing as a potential surveillance tool for complexity of infection and drug resistance markers in Plasmodium falciparum asymptomatic infections

Background

Genotyping Plasmodium falciparum subpopulations in malaria infections is an important aspect of malaria molecular epidemiology to understand within-host diversity and the frequency of drug resistance markers.

Methods

We characterized P. falciparum genetic diversity in asymptomatic infections and subsequent first febrile infections using amplicon sequencing (AmpSeq) of ama1 in Coastal Kenya. We also examined temporal changes in haplotype frequencies of mdr1, a drug-resistant marker.

Results

We found >60% of the infections were polyclonal (complexity of infection [COI] >1) and there was a reduction in COI over time. Asymptomatic infections had a significantly higher mean COI than febrile infections based on ama1 sequences (2.7 [95% confidence interval {CI}, 2.65–2.77] vs 2.22 [95% CI, 2.17–2.29], respectively). Moreover, an analysis of 30 paired asymptomatic and first febrile infections revealed that many first febrile infections (91%) were due to the presence of new ama1 haplotypes. The mdr1-YY haplotype, associated with chloroquine and amodiaquine resistance, decreased over time, while the NY (wild type) and the NF (modulates response to lumefantrine) haplotypes increased.

Conclusions

This study emphasizes the utility of AmpSeq in characterizing parasite diversity as it can determine relative proportions of clones and detect minority clones. The usefulness of AmpSeq in antimalarial drug resistance surveillance is also highlighted.

Targeted Amplicon deep sequencing of ama1 and mdr1 to track within-host P. falciparum diversity throughout treatment in a clinical drug trial

Antimalarial therapeutic efficacy studies are routinely conducted in malaria-endemic countries to assess the effectiveness of antimalarial treatment strategies. Targeted amplicon deep sequencing (TADS) uniquely identifies and quantifies genetically distinct parasites within an infection. In this study, TADS Plasmodium falciparum apical membrane antigen 1 (ama1), and multidrug resistance gene 1 (mdr1), were used to characterize the complexity of infection (COI) and drug-resistance genotypes, respectively. P. falciparum positive samples were obtained from a triple artemisinin combination therapy clinical trial conducted in 30 children under 13 years of age between 2018 and 2019 in Kilifi, Kenya. Of the 30 participants, 9 presented with recurrent parasitemia from day 26 (624h) onwards. The ama1 and mdr1 genes were amplified and sequenced, while msp1, msp2 and glurp data were obtained from the original clinical study. The COI was comparable between ama1 and msp1, msp2 and glurp, however, overall ama1 detected more haplotypes. Based on ama1, a stable number of haplotypes were detected throughout treatment up until day 3. Additionally, a recrudescent infection was identified with an ama1 haplotype initially observed at 30h and later in an unscheduled follow-up visit. Using the relative frequencies of ama1 haplotypes and parasitaemia, we identified a fast (<1h) and slow (>5h) clearing haplotype. As expected, only two mdr1 haplotypes (NF and NY) were identified based on the combination of amino acid polymorphisms at codons 86 and 184. This study highlights TADS as a sensitive tool for tracking parasite haplotypes throughout treatment and can detect variation in haplotype clearance estimates. TADS can also identify slow clearing haplotypes, a potential early sign of selection during treatment. Consequently, TADS has the capability of improving the discriminatory power to accurately distinguish recrudescences from reinfections.

Targeted Amplicon Deep Sequencing for Monitoring Antimalarial Resistance Markers in Western Kenya

Molecular surveillance of Plasmodium falciparum parasites is important to track emerging and new mutations and trends in established mutations and should serve as an early warning system for antimalarial resistance. Dried blood spots were obtained from a Plasmodium falciparum malaria survey in school children conducted across eight counties in western Kenya in 2019. Real-time PCR identified 500 P. falciparum-positive samples that were amplified at five drug resistance loci for targeted amplicon deep sequencing (TADS). The absence of important kelch 13 mutations was similar to previous findings in Kenya pre-2019, and low-frequency mutations were observed in codons 569 and 578. The chloroquine resistance transporter gene codons 76 and 145 were wild type, indicating that the parasites were chloroquine and piperaquine sensitive, respectively. The multidrug resistance gene 1 haplotypes based on codons 86, 184, and 199 were predominantly present in mixed infections with haplotypes NYT and NFT, driven by the absence of chloroquine pressure and the use of lumefantrine, respectively. The sulfadoxine-pyrimethamine resistance profile was a “superresistant” combination of triple mutations in both Pfdhfr (51I 59R 108N) and Pfdhps (436H 437G 540E), rendering sulfadoxine-pyrimethamine ineffective. TADS highlighted the low-frequency variants, allowing the early identification of new mutations, Pfmdr1 codon 199S and Pfdhfr codon 85I and emerging 164L mutations. The added value of TADS is its accuracy in identifying mixed-genotype infections and for high-throughput monitoring of antimalarial resistance markers.

A targeted amplicon sequencing panel to simultaneously identify mosquito species and Plasmodium presence across the entire Anopheles genus

Anopheles is a diverse genus of mosquitoes comprising over 500 described species, including all known human malaria vectors. While a limited number of key vector species have been studied in detail, the goal of malaria elimination calls for surveillance of all potential vector species. Here, we develop a multilocus amplicon sequencing approach that targets 62 highly variable loci in the Anopheles genome and two conserved loci in the Plasmodium mitochondrion, simultaneously revealing both the mosquito species and whether that mosquito carries malaria parasites. We also develop a cheap, nondestructive, and high-throughput DNA extraction workflow that provides template DNA from single mosquitoes for the multiplex PCR, which means specimens producing unexpected results can be returned to for morphological examination. Over 1000 individual mosquitoes can be sequenced in a single MiSeq run, and we demonstrate the panel's power to assign species identity using sequencing data for 40 species from Africa, Southeast Asia, and South America. We also show that the approach can be used to resolve geographic population structure within An. gambiae and An. coluzzii populations, as the population structure determined based on these 62 loci from over 1000 mosquitoes closely mirrors that revealed through whole genome sequencing. The end-to-end approach is quick, inexpensive, robust, and accurate, which makes it a promising technique for very large-scale mosquito genetic surveillance and vector control.

Dynamic Classification of Plasmodium vivax Malaria Recurrence: An Application of Classifying Unknown Cause of Failure in Competing Risks

A standard competing risks set-up requires both time to event and cause of failure to be fully observable for all subjects. However, in application, the cause of failure may not always be observable, thus impeding the risk assessment. In some extreme cases, none of the causes of failure is observable. In the case of a recurrent episode of Plasmodium vivax malaria following treatment, the patient may have suffered a relapse from a previous infection or acquired a new infection from a mosquito bite. In this case, the time to relapse cannot be modeled when a competing risk, a new infection, is present. The efficacy of a treatment for preventing relapse from a previous infection may be underestimated when the true cause of infection cannot be classified. In this paper, we developed a novel method for classifying the latent cause of failure under a competing risks set-up, which uses not only time to event information but also transition likelihoods between covariates at the baseline and at the time of event occurrence. Our classifier shows superior performance under various scenarios in simulation experiments. The method was applied to Plasmodium vivax infection data to classify recurrent infections of malaria.

Immediate pools of malaria infections at diagnosis combined with targeted deep sequencing accurately quantifies frequency of drug resistance mutations

Antimalarial resistance surveillance in sub-Saharan Africa is often constrained by logistical and financial challenges limiting its breadth and frequency. At two sites in Ghana, we have piloted a streamlined sample pooling process created immediately by sequential addition of positive malaria cases at the time of diagnostic testing. This streamlined process involving a single tube minimized clinical and laboratory work and provided accurate frequencies of all known drug resistance mutations after high-throughput targeted sequencing using molecular inversion probes. Our study validates this method as a cost-efficient, accurate and highly-scalable approach for drug resistance mutation monitoring that can potentially be applied to other infectious diseases such as tuberculosis.

Unraveling the Complexity of Imported Malaria Infections by Amplicon Deep Sequencing

Imported malaria and recurrent infections are becoming an emerging issue in many malaria non-endemic countries. This study aimed to determine the molecular patterns of the imported malaria infections and recurrence. Blood samples were collected from patients with imported malaria infections during 2016-2018 in Guangxi Zhuang Autonomous Region, China. Next-generation amplicon deep-sequencing approaches were used to assess parasite genetic diversity, multiplexity of infection, relapse, recrudescence, and antimalarial drug resistance. A total of 44 imported malaria cases were examined during the study, of which 35 (79.5%) had recurrent malaria infections within 1 year. The majority (91.4%) had one recurrent malaria episode, whereas two patients had two recurrences and one patient had three recurrences. A total of 19 recurrence patterns (the species responsible for primary and successive clinical episodes) were found in patients returning from malaria epidemic countries. Four parasite species were detected with a higher than usual proportion (46.2%) of non-falciparum infections or mixed-species infections. An increasing trend of recurrence infections and reduced drug treatment efficacy were observed among the cases of imported malaria. The high recurrence rate and complex patterns of imported malaria from Africa to non-endemic countries have the potential to initiate local transmission, thereby undermining efforts to eliminate locally acquired malaria. Our findings highlight the power of amplicon deep-sequencing applications in molecular epidemiological studies of the imported malaria recurrences.

Monoclonal antibodies block transmission of genetically diverse Plasmodium falciparum strains to mosquitoes

Malaria parasite transmission to mosquitoes relies on the uptake of sexual stage parasites during a blood meal and subsequent formation of oocysts on the mosquito midgut wall. Transmission-blocking vaccines (TBVs) and monoclonal antibodies (mAbs) target sexual stage antigens to interrupt human-to-mosquito transmission and may form important tools for malaria elimination. Although most epitopes of these antigens are considered highly conserved, little is known about the impact of natural genetic diversity on the functional activity of transmission-blocking antibodies. Here we measured the efficacy of three mAbs against leading TBV candidates (Pfs48/45, Pfs25 and Pfs230) in transmission assays with parasites from naturally infected donors compared to their efficacy against the strain they were raised against (NF54). Transmission-reducing activity (TRA) was measured as reduction in mean oocyst intensity. mAb 45.1 (α-Pfs48/45) and mAb 4B7 (α-Pfs25) reduced transmission of field parasites from almost all donors with IC80 values similar to NF54. Sequencing of oocysts that survived high mAb concentrations did not suggest enrichment of escape genotypes. mAb 2A2 (α-Pfs230) only reduced transmission of parasites from a minority of the donors, suggesting that it targets a non-conserved epitope. Using six laboratory-adapted strains, we revealed that mutations in one Pfs230 domain correlate with mAb gamete surface binding and functional TRA. Our findings demonstrate that, despite the conserved nature of sexual stage antigens, minor sequence variation can significantly impact the efficacy of transmission-blocking mAbs. Since mAb 45.1 shows high potency against genetically diverse strains, our findings support its further clinical development and may inform Pfs48/45 vaccine design.

AmpliCI: a high-resolution model-based approach for denoising Illumina amplicon data

Motivation

Next-generation amplicon sequencing is a powerful tool for investigating microbial communities. A main challenge is to distinguish true biological variants from errors caused by amplification and sequencing. In traditional analyses, such errors are eliminated by clustering reads within a sequence similarity threshold, usually 97%, and constructing operational taxonomic units, but the arbitrary threshold leads to low resolution and high false-positive rates. Recently developed ‘denoising’ methods have proven able to resolve single-nucleotide amplicon variants, but they still miss low-frequency sequences, especially those near more frequent sequences, because they ignore the sequencing quality information.

Results

We introduce AmpliCI, a reference-free, model-based method for rapidly resolving the number, abundance and identity of error-free sequences in massive Illumina amplicon datasets. AmpliCI considers the quality information and allows the data, not an arbitrary threshold or an external database, to drive conclusions. AmpliCI estimates a finite mixture model, using a greedy strategy to gradually select error-free sequences and approximately maximize the likelihood. AmpliCI has better performance than three popular denoising methods, with acceptable computation time and memory usage.

Availability and implementation

Source code is available at https://github.com/DormanLab/AmpliCI.

Supplementary information

Supplementary material are available at Bioinformatics online.

Should deep-sequenced amplicons become the new gold-standard for analysing malaria drug clinical trials?

Regulatory clinical trials are required to ensure the continued supply and deployment of effective antimalarial drugs. Patient follow-up in such trials typically lasts several weeks, as the drugs have long half-lives and new infections often occur during this period. “Molecular correction” is therefore used to distinguish drug failures from new infections. The current WHO-recommended method for molecular correction uses length-polymorphic alleles at highly diverse loci but is inherently poor at detecting low-density clones in polyclonal infections. This likely leads to substantial underestimates of failure rates, delaying the replacement of failing drugs with potentially lethal consequences. Deep-sequenced amplicons (AmpSeq) substantially increase the detectability of low-density clones and may offer a new “gold standard” for molecular correction. Pharmacological simulation of clinical trials was used to evaluate the suitability of AmpSeq for molecular correction. We investigated the impact of factors such as the number of amplicon loci analyzed, the informatics criteria used to distinguish genotyping “noise” from real low-density signals, the local epidemiology of malaria transmission, and the potential impact of genetic signals from gametocytes. AmpSeq greatly improved molecular correction and provided accurate drug failure rate estimates. The use of 3 to 5 amplicons was sufficient, and simple, nonstatistical criteria could be used to classify recurrent infections as drug failures or new infections. These results suggest AmpSeq is strongly placed to become the new standard for molecular correction in regulatory trials, with potential extension into routine surveillance once the requisite technical support becomes established.

Deep sequencing to detect diversity of Trypanosoma cruzi infection in patients co-infected with HIV and Chagas disease

Chagas disease, caused by Trypanosoma cruzi, can reactivate and cause severe acute disease in immunocompromised patients such as those infected with HIV. We conducted amplicon deep sequencing of a 327-base pair fragment of the tcscd5 gene using an Ion Torrent PGM directly from clinical samples from HIV patients with high parasitemia. We describe the within host diversity, both characterizing the discrete typing unit (DTUs) of the infections and confirming the presence of multi-strain infections, directly from clinical samples. This method can rapidly provide information on the genetic diversity of T. cruzi infection, which can have direct impacts on clinical disease.

A novel CRISPR-based malaria diagnostic capable of Plasmodium detection, species differentiation, and drug-resistance genotyping

BACKGROUND

CRISPR-based diagnostics are a new class of highly sensitive and specific assays with multiple applications in infectious disease diagnosis. SHERLOCK, or Specific High-Sensitivity Enzymatic Reporter UnLOCKing, is one such CRISPR-based diagnostic that combines recombinase polymerase pre-amplification, CRISPR-RNA base-pairing, and LwCas13a activity for nucleic acid detection.

METHODS

We developed SHERLOCK assays capable of detecting all Plasmodium species known to cause human malaria and species-specific detection of P. vivax and P. falciparum, the species responsible for the majority of malaria cases worldwide. We further tested these assays using a diverse panel of clinical samples from the Democratic Republic of the Congo, Uganda, and Thailand and pools of Anopheles mosquitoes from Thailand. In addition, we developed a prototype SHERLOCK assay capable of detecting the dihydropteroate synthetase (dhps) single nucleotide variant A581G associated with P. falciparum sulfadoxine resistance.

FINDINGS

The suite of Plasmodium assays achieved analytical sensitivities ranging from 2•5-18•8 parasites per reaction when tested against laboratory strain genomic DNA. When compared to real-time PCR, the P. falciparum assay achieved 94% sensitivity and 94% specificity during testing of 123 clinical samples. Compared to amplicon-based deep sequencing, the dhps SHERLOCK assay achieved 73% sensitivity and 100% specificity when applied to a panel of 43 clinical samples, with false-negative calls only at lower parasite densities.

INTERPRETATION

These novel SHERLOCK assays demonstrate the versatility of CRISPR-based diagnostics and their potential as a new generation of molecular tools for malaria diagnosis and surveillance.

FUNDING

National Institutes of Health (T32GM007092, R21AI148579, K24AI134990, R01AI121558, UL1TR002489, P30CA016086).

Under the Radar: Epidemiology of Plasmodium ovale in the Democratic Republic of the Congo

BACKGROUND

Plasmodium ovale is an understudied malaria species prevalent throughout much of sub-Saharan Africa. Little is known about the distribution of ovale malaria and risk factors for infection in areas of high malaria endemicity.

METHODS

Using the 2013 Democratic Republic of the Congo (DRC) Demographic and Health Survey, we conducted a risk factor analysis for P. ovale infections. We evaluated geographic clustering of infections and speciated to P. ovale curtisi and P. ovale wallikeri through deep sequencing.

RESULTS

Of 18 149 adults tested, we detected 143 prevalent P. ovale infections (prevalence estimate 0.8%; 95% confidence interval [CI], .59%-.98%). Prevalence ratios (PR) for significant risk factors were: male sex PR = 2.12 (95% CI, 1.38-3.26), coprevalent P. falciparum PR = 3.52 (95% CI, 2.06-5.99), and rural residence PR = 2.19 (95% CI, 1.31-3.66). P. ovale was broadly distributed throughout the DRC; an elevated cluster of infections was detected in the south-central region. Speciation revealed P. ovale curtisi and P. ovale wallikeri circulating throughout the country.

CONCLUSIONS

P. ovale persists broadly in the DRC, a high malaria burden country. For successful elimination of all malaria species, P. ovale needs to be on the radar of malaria control programs.

Comparison of Methods for Picking the Operational Taxonomic Units From Amplicon Sequences

With the advent of next-generation sequencing technology, it has become convenient and cost efficient to thoroughly characterize the microbial diversity and taxonomic composition in various environmental samples. Millions of sequencing data can be generated, and how to utilize this enormous sequence resource has become a critical concern for microbial ecologists. One particular challenge is the OTUs (operational taxonomic units) picking in 16S rRNA sequence analysis. Lucky, this challenge can be directly addressed by sequence clustering that attempts to group similar sequences. Therefore, numerous clustering methods have been proposed to help to cluster 16S rRNA sequences into OTUs. However, each method has its clustering mechanism, and different methods produce diverse outputs. Even a slight parameter change for the same method can also generate distinct results, and how to choose an appropriate method has become a challenge for inexperienced users. A lot of time and resources can be wasted in selecting clustering tools and analyzing the clustering results. In this study, we introduced the recent advance of clustering methods for OTUs picking, which mainly focus on three aspects: (i) the principles of existing clustering algorithms, (ii) benchmark dataset construction for OTU picking and evaluation metrics, and (iii) the performance of different methods with various distance thresholds on benchmark datasets. This paper aims to assist biological researchers to select the reasonable clustering methods for analyzing their collected sequences and help algorithm developers to design more efficient sequences clustering methods.