swga: a primer design toolkit for selective whole genome amplification

Distinguishing recrudescence from reinfection in lymphatic filariasis

BACKGROUND

The Global Program to Eliminate Lymphatic Filariasis (GPELF) is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimise interventions, it is critical to differentiate between recrudescence and new infections. Since adult filariae are inaccessible in humans, deriving a method that relies on the offspring microfilariae (mf) is necessary.

METHODS

We developed a genome amplification and kinship analysis-based approach using Brugia malayi samples from gerbils, and applied it to analyse Wuchereria bancrofti mf from humans in Côte d'Ivoire. We examined the pre-treatment genetic diversity in 269 mf collected from 18 participants, and further analysed 1-year post-treatment samples of 74 mf from 4 participants. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference. To enrich parasite DNA from samples contaminated with host DNA, a whole-exome capture panel was created for W. bancrofti.

FINDINGS

By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one participant and recrudescence in three participants. The estimated number of reproductively active adult females ranged between 3 and 11 in the studied participants. Population structure analysis revealed genetically distinct parasites in Côte d'Ivoire compared to samples from other countries. Exome capture identified protein-coding variants with ∼95% genotype concordance rate.

INTERPRETATION

We have generated resources to facilitate the development of molecular genetic tools that can estimate adult worm burdens and monitor parasite populations, thus providing essential information for the successful implementation of GPELF.

FUNDING

This work was financially supported by the Bill and Melinda Gates Foundation (https://www.gatesfoundation.org) under grant OPP1201530 (Co-PIs PUF & Gary J. Weil). B. malayi parasite material was generated with support of the Foundation for Barnes Jewish Hospital (PUF). In addition, the development of computational methods was supported by the National Institutes of Health under grants AI144161 (MM) and AI146353 (MM). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Advances in Host Depletion and Pathogen Enrichment Methods for Rapid Sequencing-Based Diagnosis of Bloodstream Infection

Bloodstream infection remains a major cause of morbidity and death worldwide. Timely and appropriate treatment can reduce mortality among critically ill patients. Current diagnostic methods are too slow to inform precise antibiotic choice, leading to the prescription of empirical antibiotics, which may fail to cover the resistance profile of the pathogen, risking poor patient outcomes. Additionally, overuse of broad-spectrum antibiotics may lead to more resistant organisms, putting further pressure on the dwindling pipeline of antibiotics, and risk transmission of these resistant organisms in the health care environment. Therefore, rapid diagnostics are urgently required to better inform antibiotic choice early in the course of treatment. Sequencing offers great promise in reducing time to microbiological diagnosis; however, the amount of host DNA compared with the pathogen in patient samples presents a significant obstacle. To address this, various host-depletion and bacterial-enrichment strategies have been used in samples, such as saliva, urine, or tissue. However, these methods have yet to be collectively integrated and/or extensively explored for rapid bloodstream infection diagnosis. Although most of these workflows possess individual strengths, their lack of analytical/clinical sensitivity and/or comprehensiveness demands additional improvements or synergistic application. Therefore, this review provides a distinctive classification system for these methods based on their working principles to guide future research, discusses their strengths and limitations, and explores potential avenues for improvement.

Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.

New reference genomes to distinguish the sympatric malaria parasites, Plasmodium ovale curtisi and Plasmodium ovale wallikeri

Despite Plasmodium ovale curtisi (Poc) and wallikeri (Pow) being important human-infecting malaria parasites that are widespread across Africa and Asia, little is known about their genome diversity. Morphologically identical, Poc and Pow are indistinguishable and commonly misidentified. Recent rises in the incidence of Poc/Pow infections have renewed efforts to address fundamental knowledge gaps in their biology, and to develop diagnostic tools to understand their epidemiological dynamics and malaria burden. A major roadblock has been the incompleteness of available reference assemblies (PocGH01, PowCR01; ~ 33.5 Mbp). Here, we applied multiple sequencing platforms and advanced bioinformatics tools to generate new reference genomes, Poc221 (South Sudan; 36.0 Mbp) and Pow222 (Nigeria; 34.3 Mbp), with improved nuclear genome contiguity (> 4.2 Mbp), annotation and completeness (> 99% Plasmodium spp., single copy orthologs). Subsequent sequencing of 6 Poc and 15 Pow isolates from Africa revealed a total of 22,517 and 43,855 high-quality core genome SNPs, respectively. Genome-wide levels of nucleotide diversity were determined to be 2.98 × 10-4 (Poc) and 3.43 × 10-4 (Pow), comparable to estimates for other Plasmodium species. Overall, the new reference genomes provide a robust foundation for dissecting the biology of Poc/Pow, their population structure and evolution, and will contribute to uncovering the recombination barrier separating these species.

Distinguishing recrudescence from reinfection in lymphatic filariasis: a genomics-based approach for monitoring worm burden

Background

The Global Program to Eliminate Lymphatic Filariasis is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimize interventions, it is essential to differentiate between recrudescence (poor drug response and persistent infection) and new infections (ongoing transmission). Since adult filariae are inaccessible in humans, an approach that relies on genotyping the offspring microfilariae (mf) is required.

Methods

We utilized Brugia malayi adults and mf obtained from gerbils with a known pedigree to develop and validate our whole-genome amplification and kinship analysis approach. We then sequenced the genomes of Wuchereria bancrofti mf from infected humans from Côte d'Ivoire (CDI), characterized the population genetic diversity, and made inferences about the adult breeders. We developed a whole-exome capture panel for W. bancrofti to enrich parasite nuclear DNA from lower-quality samples contaminated with host DNA.

Results

We established a robust analysis pipeline using B. malayi adult and mf. We estimated the pre-treatment genetic diversity in W. bancrofti from 269 mf collected from 18 individuals, and further analyzed 1-year post-treatment samples of 74 mf from 4 individuals. By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one subject and recrudescence in three subjects. Estimated reproductively active adult females ranged between 3 and 9 in the studied subjects. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference, even when the genetic diversity was low. Population structure analysis revealed genetically distinct parasites among our CDI samples. Sequence composition and variant analysis of whole-exome libraries showed that the hybridization capture approach can effectively enrich parasite nuclear DNA and identify protein-coding variants with ∼95% genotype concordance rate.

Conclusions

We have generated resources to facilitate development of field-deployable genotyping tools that can estimate worm burdens and monitor parasite populations. These tools are essential for the success of lymphatic filariasis MDA programs. With further expansion of the databases to include geographically diverse samples, we will be able to spatially track parasite movement associated with host/vector migration.

MultiPrime: A reliable and efficient tool for targeted next-generation sequencing

We present multiPrime, a novel tool that automatically designs minimal primer sets for targeted next-generation sequencing, tailored to specific microbiomes or genes. MultiPrime enhances primer coverage by designing primers with mismatch tolerance and ensures both high compatibility and specificity. We evaluated the performance of multiPrime using a data set of 43,016 sequences from eight viruses. Our results demonstrated that multiPrime outperformed conventional tools, and the primer set designed by multiPrime successfully amplified the target amplicons. Furthermore, we expanded the application of multiPrime to 30 types of viruses and validated the work efficacy of multiPrime-designed primers in 80 clinical specimens. The subsequent sequencing outcomes from these primers indicated a sensitivity of 94% and a specificity of 89%.

Drug resistance profiling of asymptomatic and low-density Plasmodium falciparum malaria infections on Ngodhe island, Kenya, using custom dual-indexing next-generation sequencing

Malaria control initiatives require rapid and reliable methods for the detection and monitoring of molecular markers associated with antimalarial drug resistance in Plasmodium falciparum parasites. Ngodhe island, Kenya, presents a unique malaria profile, with lower P. falciparum incidence rates than the surrounding region, and a high proportion of sub-microscopic and low-density infections. Here, using custom dual-indexing and Illumina next generation sequencing, we generate resistance profiles on seventy asymptomatic and low-density P. falciparum infections from a mass drug administration program implemented on Ngodhe island between 2015 and 2016. Our assay encompasses established molecular markers on the Pfcrt, Pfmdr1, Pfdhps, Pfdhfr, and Pfk13 genes. Resistance markers for sulfadoxine-pyrimethamine were identified at high frequencies, including a quintuple mutant haplotype (Pfdhfr/Pfdhps: N51I, C59R, S108N/A437G, K540E) identified in 62.2% of isolates. The Pfdhps K540E biomarker, used to inform decision making for intermittent preventative treatment in pregnancy, was identified in 79.2% of isolates. Several variants on Pfmdr1, associated with reduced susceptibility to quinolones and lumefantrine, were also identified (Y184F 47.1%; D1246Y 16.0%; N86 98%). Overall, we have presented a low-cost and extendable approach that can provide timely genetic profiles to inform clinical and surveillance activities, especially in settings with abundant low-density infections, seeking malaria elimination.

Genome-wide genetic variation and molecular surveillance of drug resistance in Plasmodium falciparum isolates from asymptomatic individuals in Ouélessébougou, Mali

Sequence analysis of Plasmodium falciparum parasites is informative in ensuring sustained success of malaria control programmes. Whole-genome sequencing technologies provide insights into the epidemiology and genome-wide variation of P. falciparum populations and can characterise geographical as well as temporal changes. This is particularly important to monitor the emergence and spread of drug resistant P. falciparum parasites which is threatening malaria control programmes world-wide. Here, we provide a detailed characterisation of genome-wide genetic variation and drug resistance profiles in asymptomatic individuals in South-Western Mali, where malaria transmission is intense and seasonal, and case numbers have recently increased. Samples collected from Ouélessébougou, Mali (2019-2020; n = 87) were sequenced and placed in the context of older Malian (2007-2017; n = 876) and African-wide (n = 711) P. falciparum isolates. Our analysis revealed high multiclonality and low relatedness between isolates, in addition to increased frequencies of molecular markers for sulfadoxine-pyrimethamine and lumefantrine resistance, compared to older Malian isolates. Furthermore, 21 genes under selective pressure were identified, including a transmission-blocking vaccine candidate (pfCelTOS) and an erythrocyte invasion locus (pfdblmsp2). Overall, our work provides the most recent assessment of P. falciparum genetic diversity in Mali, a country with the second highest burden of malaria in West Africa, thereby informing malaria control activities.

A fast machine-learning-guided primer design pipeline for selective whole genome amplification

Addressing many of the major outstanding questions in the fields of microbial evolution and pathogenesis will require analyses of populations of microbial genomes. Although population genomic studies provide the analytical resolution to investigate evolutionary and mechanistic processes at fine spatial and temporal scales-precisely the scales at which these processes occur-microbial population genomic research is currently hindered by the practicalities of obtaining sufficient quantities of the relatively pure microbial genomic DNA necessary for next-generation sequencing. Here we present swga2.0, an optimized and parallelized pipeline to design selective whole genome amplification (SWGA) primer sets. Unlike previous methods, swga2.0 incorporates active and machine learning methods to evaluate the amplification efficacy of individual primers and primer sets. Additionally, swga2.0 optimizes primer set search and evaluation strategies, including parallelization at each stage of the pipeline, to dramatically decrease program runtime. Here we describe the swga2.0 pipeline, including the empirical data used to identify primer and primer set characteristics, that improve amplification performance. Additionally, we evaluate the novel swga2.0 pipeline by designing primer sets that successfully amplify Prevotella melaninogenica, an important component of the lung microbiome in cystic fibrosis patients, from samples dominated by human DNA.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Development and Optimization of a Selective Whole-Genome Amplification To Study Plasmodium ovale Spp

Since 2010, the human-infecting malaria parasite Plasmodium ovale spp. has been divided into two genetically distinct species, P. ovale wallikeri and P. ovale curtisi. In recent years, application of whole-genome sequencing (WGS) to P. ovale spp. allowed to get a better understanding of its evolutionary history and discover some specific genetic patterns. Nevertheless, WGS data from P. ovale spp. are still scarce due to several drawbacks, including a high level of human DNA contamination in blood samples, infections with commonly low parasite density, and the lack of robust in vitro culture. Here, we developed two selective whole-genome amplification (sWGA) protocols that were tested on six P. ovale wallikeri and five P. ovale curtisi mono-infection clinical samples. Blood leukodepletion by a cellulose-based filtration was used as the gold standard for intraspecies comparative genomics with sWGA. We also demonstrated the importance of genomic DNA preincubation with the endonuclease McrBC to optimize P. ovale spp. sWGA. We obtained high-quality WGS data with more than 80% of the genome covered by ≥5 reads for each sample and identified more than 5,000 unique single-nucleotide polymorphisms (SNPs) per species. We also identified some amino acid changes in pocdhfr and powdhfr for which similar mutations in P. falciparum and P. vivax are associated with pyrimethamine or cycloguanil resistance. In conclusion, we developed two sWGA protocols for P. ovale spp. WGS that will help to design much-needed large-scale P. ovale spp. population studies. IMPORTANCE Plasmodium ovale spp. has the ability to cause relapse, defined as recurring asexual parasitemia originating from liver-dormant forms. Whole-genome sequencing (WGS) data are of importance to identify putative molecular markers associated with relapse or other virulence mechanisms. Due to low parasitemia encountered in P. ovale spp. infections and no in vitro culture available, WGS of P. ovale spp. is challenging. Blood leukodepletion by filtration has been used, but no technique exists yet to increase the quantity of parasite DNA over human DNA when starting from genomic DNA extracted from whole blood. Here, we demonstrated that selective whole-genome amplification (sWGA) is an easy-to-use protocol to obtain high-quality WGS data for both P. ovale spp. species from unprocessed blood samples. The new method will facilitate P. ovale spp. population genomic studies.

Genomic sequence capture of Plasmodium relictum in experimentally infected birds

Background

Sequencing parasite genomes in the presence of host DNA is challenging. Sequence capture can overcome this problem by using RNA probes that hybridize with the parasite DNA and then are removed from solution, thus isolating the parasite DNA for efficient sequencing.

Methods

Here we describe a set of sequence capture probes designed to target 1035 genes (c. 2.5 Mbp) of the globally distributed avian haemosporidian parasite, Plasmodium relictum. Previous sequence capture studies of avian haemosporidians from the genus Haemoproteus have shown that sequencing success depends on parasitemia, with low-intensity, chronic infections (typical of most infected birds in the wild) often being difficult to sequence. We evaluate the relationship between parasitemia and sequencing success using birds experimentally infected with P. relictum and kept under laboratory conditions.

Results

We confirm the dependence of sequencing success on parasitemia. Sequencing success was low for birds with low levels of parasitemia (< 1% infected red blood cells) and high for birds with higher levels of parasitemia. Plasmodium relictum is composed of multiple lineages defined by their mitochondrial DNA haplotype including three that are widespread (SGS1, GRW11, and GRW4); the probes successfully isolated DNA from all three. Furthermore, we used data from 25 genes to describe both among- and within-lineage genetic variation. For example, two samples of SGS1 isolated from different host species differed by 11 substitutions across those 25 genes.

Conclusions

The sequence capture approach we describe will allow for the generation of genomic data that will contribute to our understanding of the population genetic structure and evolutionary history of P. relictum, an extreme host generalist and widespread parasite.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05373-w.

Selective Whole-Genome Amplification as a Tool to Enrich Specimens with Low Treponema pallidum Genomic DNA Copies for Whole-Genome Sequencing

Syphilis is a sexually transmitted, disseminated acute and chronic infection caused by the bacterial pathogen Treponema pallidum subspecies pallidum . Primary syphilis typically presents as single or multiple mucocutaneous lesions and, if left untreated, can progress through multiple stages with various clinical manifestations.

Redondovirus Diversity and Evolution on Global, Individual, and Molecular Scales

Redondoviridae is a newly established family of circular Rep-encoding single-stranded (CRESS) DNA viruses found in the human ororespiratory tract. Redondoviruses were previously found in ∼15% of respiratory specimens from U.S. urban subjects; levels were elevated in individuals with periodontitis or critical illness. Here, we report higher redondovirus prevalence in saliva samples: four rural African populations showed 61 to 82% prevalence, and an urban U.S. population showed 32% prevalence. Longitudinal, limiting-dilution single-genome sequencing revealed diverse strains of both redondovirus species (Brisavirus and Vientovirus) in single individuals, persistence over time, and evidence of intergenomic recombination. Computational analysis of viral genomes identified a recombination hot spot associated with a conserved potential DNA stem-loop structure. To assess the possible role of this site in recombination, we carried out in vitro studies which showed that this potential stem-loop was cleaved by the virus-encoded Rep protein. In addition, in reconstructed reactions, a Rep-DNA covalent intermediate was shown to mediate DNA strand transfer at this site. Thus, redondoviruses are highly prevalent in humans, found in individuals on multiple continents, heterogeneous even within individuals and encode a Rep protein implicated in facilitating recombination. IMPORTANCE Redondoviridae is a recently established family of DNA viruses predominantly found in the human respiratory tract and associated with multiple clinical conditions. In this study, we found high redondovirus prevalence in saliva from urban North American individuals and nonindustrialized African populations in Botswana, Cameroon, Ethiopia, and Tanzania. Individuals on both continents harbored both known redondovirus species. Global prevalence of both species suggests that redondoviruses have long been associated with humans but have remained undetected until recently due to their divergent genomes. By sequencing single redondovirus genomes in longitudinally sampled humans, we found that redondoviruses persisted over time within subjects and likely evolve by recombination. The Rep protein encoded by redondoviruses catalyzes multiple reactions in vitro, consistent with a role in mediating DNA replication and recombination. In summary, we identify high redondovirus prevalence in humans across multiple continents, longitudinal heterogeneity and persistence, and potential mechanisms of redondovirus evolution by recombination.

Characterizing the genomic variation and population dynamics of Plasmodium falciparum malaria parasites in and around Lake Victoria, Kenya

Characterising the genomic variation and population dynamics of Plasmodium falciparum parasites in high transmission regions of Sub-Saharan Africa is crucial to the long-term efficacy of regional malaria elimination campaigns and eradication. Whole-genome sequencing (WGS) technologies can contribute towards understanding the epidemiology and structural variation landscape of P. falciparum populations, including those within the Lake Victoria basin, a region of intense transmission. Here we provide a baseline assessment of the genomic diversity of P. falciparum isolates in the Lake region of Kenya, which has sparse genetic data. Lake region isolates are placed within the context of African-wide populations using Illumina WGS data and population genomic analyses. Our analysis revealed that P. falciparum isolates from Lake Victoria form a cluster within the East African parasite population. These isolates also appear to have distinct ancestral origins, containing genome-wide signatures from both Central and East African lineages. Known drug resistance biomarkers were observed at similar frequencies to those of East African parasite populations, including the S160N/T mutation in the pfap2mu gene, which has been associated with delayed clearance by artemisinin-based combination therapy. Overall, our work provides a first assessment of P. falciparum genetic diversity within the Lake Victoria basin, a region targeting malaria elimination.

Full Molecular Typing of Neisseria meningitidis Directly from Clinical Specimens for Outbreak Investigation

Neisseria meningitidis is a leading cause of bacterial meningitis and sepsis worldwide and an occasional cause of meningococcal urethritis. When isolates are unavailable for surveillance or outbreak investigations, molecular characterization of pathogens needs to be performed directly from clinical specimens, such as cerebrospinal fluid (CSF), blood, or urine. However, genome sequencing of specimens is challenging because of low bacterial and high human DNA abundances. We developed selective whole-genome amplification (SWGA), an isothermal multiple-displacement amplification-based method, to efficiently enrich, sequence, and de novo assemble N. meningitidis DNA from clinical specimens with low bacterial loads. SWGA was validated with 12 CSF specimens from invasive meningococcal disease cases and 12 urine specimens from meningococcal urethritis cases. SWGA increased the mean proportion of N. meningitidis reads by 2 to 3 orders of magnitude, enabling identification of at least 90% of the 1,605 N. meningitidis core genome loci for 50% of the specimens. The validated method was used to investigate two meningitis outbreaks recently reported in Togo and Burkina Faso. Twenty-seven specimens with low bacterial loads were processed by SWGA before sequencing, and 12 of 27 were successfully assembled to obtain the full molecular typing and vaccine antigen profile of the N. meningitidis pathogen, thus enabling thorough characterization of outbreaks. This method is particularly important for enhancing molecular surveillance in regions with low culture rates. SWGA produces enough reads for phylogenetic and allelic analysis at a low cost. More importantly, the procedure can be extended to enrich other important human bacterial pathogens.

Selective whole genome amplification of Plasmodium malariae DNA from clinical samples reveals insights into population structure

The genomic diversity of Plasmodium malariae malaria parasites is understudied, partly because infected individuals tend to present with low parasite densities, leading to difficulties in obtaining sufficient parasite DNA for genome analysis. Selective whole genome amplification (SWGA) increases the relative levels of pathogen DNA in a clinical sample, but has not been adapted for P. malariae parasites. Here we design customized SWGA primers which successfully amplify P. malariae DNA extracted directly from unprocessed clinical blood samples obtained from patients with P. malariae-mono-infections from six countries, and further test the efficacy of SWGA on mixed infections with other Plasmodium spp. SWGA enables the successful whole genome sequencing of samples with low parasite density (i.e. one sample with a parasitaemia of 0.0064% resulted in 44% of the genome covered by ≥ 5 reads), leading to an average 14-fold increase in genome coverage when compared to unamplified samples. We identify a total of 868,476 genome-wide SNPs, of which 194,709 are unique across 18 high-quality isolates. After exclusion of the hypervariable subtelomeric regions, a high-quality core subset of 29,899 unique SNPs is defined. Population genetic analysis suggests that P. malariae parasites display clear geographical separation by continent. Further, SWGA successfully amplifies genetic regions of interest such as orthologs of P. falciparum drug resistance-associated loci (Pfdhfr, Pfdhps, Pfcrt, Pfk13 and Pfmdr1), and several non-synonymous SNPs were detected in these genes. In conclusion, we have established a robust SWGA approach that can assist whole genome sequencing of P. malariae, and thereby facilitate the implementation of much-needed large-scale multi-population genomic studies of this neglected malaria parasite. As demonstrated in other Plasmodia, such genetic diversity studies can provide insights into the biology underlying the disease and inform malaria surveillance and control measures.

Next-Generation Molecular Surveillance of TriTryp Diseases

Elimination programs targeting TriTryp diseases (Leishmaniasis, Chagas' disease, human African trypanosomiasis) significantly reduced the number of cases. Continued surveillance is crucial to sustain this progress, but parasite molecular surveillance by genotyping is currently lacking. We explain here which epidemiological questions of public health and clinical relevance could be answered by means of molecular surveillance. Whole-genome sequencing (WGS) for molecular surveillance will be an important added value, where we advocate that preference should be given to direct sequencing of the parasite's genome in host tissues instead of analysis of cultivated isolates. The main challenges here, and recent technological advances, are discussed. We conclude with a series of recommendations for implementing whole-genome sequencing for molecular surveillance.

Falciparum malaria from coastal Tanzania and Zanzibar remains highly connected despite effective control efforts on the archipelago

BACKGROUND

Tanzania's Zanzibar archipelago has made significant gains in malaria control over the last decade and is a target for malaria elimination. Despite consistent implementation of effective tools since 2002, elimination has not been achieved. Importation of parasites from outside of the archipelago is thought to be an important cause of malaria's persistence, but this paradigm has not been studied using modern genetic tools.

METHODS

Whole-genome sequencing (WGS) was used to investigate the impact of importation, employing population genetic analyses of Plasmodium falciparum isolates from both the archipelago and mainland Tanzania. Ancestry, levels of genetic diversity and differentiation, patterns of relatedness, and patterns of selection between these two populations were assessed by leveraging recent advances in deconvolution of genomes from polyclonal malaria infections.

RESULTS

Significant decreases in the effective population sizes were inferred in both populations that coincide with a period of decreasing malaria transmission in Tanzania. Identity by descent analysis showed that parasites in the two populations shared long segments of their genomes, on the order of 5 cM, suggesting shared ancestry within the last 10 generations. Even with limited sampling, two of isolates between the mainland and Zanzibar were identified that are related at the expected level of half-siblings, consistent with recent importation.

CONCLUSIONS

These findings suggest that importation plays an important role for malaria incidence on Zanzibar and demonstrate the value of genomic approaches for identifying corridors of parasite movement to the island.

Redondoviridae, a Family of Small, Circular DNA Viruses of the Human Oro-Respiratory Tract Associated with Periodontitis and Critical Illness

The global virome is largely uncharacterized but is now being unveiled by metagenomic DNA sequencing. Exploring the human respiratory virome, in particular, can provide insights into oro-respiratory diseases. Here, we use metagenomics to identify a family of small circular DNA viruses-named Redondoviridae-associated with human diseases. We first identified two redondovirus genomes from bronchoalveolar lavage samples from human lung donors. We then queried thousands of metagenomic samples and recovered 17 additional complete redondovirus genomes. Detections were exclusively in human samples and mostly from respiratory tract and oro-pharyngeal sites, where Redondoviridae was the second most prevalent eukaryotic DNA virus family. Redondovirus sequences were associated with periodontal disease, and abundances decreased with treatment. Some critically ill patients in a medical intensive care unit were found to harbor high levels of redondoviruses in respiratory samples. These results suggest that redondoviruses colonize human oro-respiratory sites and can bloom in several human disorders.

Selective whole genome amplification and sequencing of Coxiella burnetii directly from environmental samples

Whole genome sequencing (WGS) is a widely available, inexpensive means of providing a wealth of information about an organism's diversity and evolution. However, WGS for many pathogenic bacteria remain limited because they are difficult, slow and/or dangerous to culture. To avoid culturing, metagenomic sequencing can be performed directly on samples, but the sequencing effort required to characterize low frequency organisms can be expensive. Recently developed methods for selective whole genome amplification (SWGA) can enrich target DNA to provide efficient sequencing. We amplified Coxiella burnetii (a bacterial select agent and human/livestock pathogen) from 3 three environmental samples that were overwhelmed with host DNA. The 68- to 147-fold enrichment of the bacterial sequences provided enough genome coverage for SNP analyses and phylogenetic placement. SWGA is a valuable tool for the study of difficult-to-culture organisms and has the potential to facilitate high-throughput population characterizations as well as targeted epidemiological or forensic investigations.

Human Migration and the Spread of the Nematode Parasite Wuchereria bancrofti

The human disease lymphatic filariasis causes the debilitating effects of elephantiasis and hydrocele. Lymphatic filariasis currently affects the lives of 90 million people in 52 countries. There are three nematodes that cause lymphatic filariasis, Brugia malayi, Brugia timori, and Wuchereria bancrofti, but 90% of all cases of lymphatic filariasis are caused solely by W. bancrofti (Wb). Here we use population genomics to reconstruct the probable route and timing of migration of Wb strains that currently infect Africa, Haiti, and Papua New Guinea (PNG). We used selective whole genome amplification to sequence 42 whole genomes of single Wb worms from populations in Haiti, Mali, Kenya, and PNG. Our results are consistent with a hypothesis of an Island Southeast Asia or East Asian origin of Wb. Our demographic models support divergence times that correlate with the migration of human populations. We hypothesize that PNG was infected at two separate times, first by the Melanesians and later by the migrating Austronesians. The migrating Austronesians also likely introduced Wb to Madagascar where later migrations spread it to continental Africa. From Africa, Wb spread to the New World during the transatlantic slave trade. Genome scans identified 17 genes that were highly differentiated among Wb populations. Among these are genes associated with human immune suppression, insecticide sensitivity, and proposed drug targets. Identifying the distribution of genetic diversity in Wb populations and selection forces acting on the genome will build a foundation to test future hypotheses and help predict response to current eradication efforts.

Molecular epidemiology of giardiasis from a veterinary perspective

A total of eight Giardia species are accepted. These include: Giardia duodenalis (syn. Giardia intestinalis and Giardia lamblia), which infects humans and animals, Giardia agilis, Giardia ardeae, Giardia psittaci, Giardia muris, Giardia microti, Giardia peramelis and G. cricetidarum, which infect non-human hosts including amphibians, birds, rodents and marsupials. Giardia duodenalis is a species complex consisting of eight assemblages (A-H), with assemblages A and B the dominant assemblages in humans. Molecular studies to date on the zoonotic potential of Giardia in animals are problematic and are hampered by lack of concordance between loci. Livestock (cattle, sheep, goats and pigs) are predominantly infected with G. duodenalis assemblage E, which has recently been shown to be zoonotic, followed by assemblage A. In cats and dogs, assemblages A, B, C, D and F are commonly reported but relatively few studies have conducted molecular typing of humans and their pets and the results are contradictory with some studies support zoonotic transmission but the majority of studies suggesting separate transmission cycles. Giardia also infects a broad range of wildlife hosts and although much less well studied, host-adapted species as well as G. duodenalis assemblages (A-H) have been identified. Fish and other aquatic wildlife represent a source of infection for humans with Giardia via water contamination and/or consumption of undercooked fish and interestingly, assemblage B and A predominated in the two molecular studies conducted to date. Our current knowledge of the transmission dynamics of Giardia is still poor and the development of more discriminatory typing tools such as whole genome sequencing (WGS) of Giardia isolates is therefore essential.

Whole genome sequencing of amplified Plasmodium knowlesi DNA from unprocessed blood reveals genetic exchange events between Malaysian Peninsular and Borneo subpopulations

The zoonotic Plasmodium knowlesi parasite is the most common cause of human malaria in Malaysia. Genetic analysis has shown that the parasites are divided into three subpopulations according to their geographic origin (Peninsular or Borneo) and, in Borneo, their macaque host (Macaca fascicularis or M. nemestrina). Whilst evidence suggests that genetic exchange events have occurred between the two Borneo subpopulations, the picture is unclear in less studied Peninsular strains. One difficulty is that P. knowlesi infected individuals tend to present with low parasitaemia leading to samples with insufficient DNA for whole genome sequencing. Here, using a parasite selective whole genome amplification approach on unprocessed blood samples, we were able to analyse recent genomes sourced from both Peninsular Malaysia and Borneo. The analysis provides evidence that recombination events are present in the Peninsular Malaysia parasite subpopulation, which have acquired fragments of the M. nemestrina associated subpopulation genotype, including the DBPβ and NBPXa erythrocyte invasion genes. The NBPXb invasion gene has also been exchanged within the macaque host-associated subpopulations of Malaysian Borneo. Our work provides strong evidence that exchange events are far more ubiquitous than expected and should be taken into consideration when studying the highly complex P. knowlesi population structure.