Hybrid selection for sequencing pathogen genomes from clinical samples

Strengthening Bordetella pertussis genomic surveillance by direct sequencing of residual positive specimens

Whole-genome sequencing (WGS) of microbial pathogens recovered from patients with infectious disease facilitates high-resolution strain characterization and molecular epidemiology. However, increasing reliance on culture-independent methods to diagnose infectious diseases has resulted in few isolates available for WGS. Here, we report a novel culture-independent approach to genome characterization of Bordetella pertussis, the causative agent of pertussis and a paradigm for insufficient genomic surveillance due to limited culture of clinical isolates. Sequencing libraries constructed directly from residual pertussis-positive diagnostic nasopharyngeal specimens were hybridized with biotinylated RNA "baits" targeting B. pertussis fragments within complex mixtures that contained high concentrations of host and microbial background DNA. Recovery of B. pertussis genome sequence data was evaluated with mock and pooled negative clinical specimens spiked with reducing concentrations of either purified DNA or inactivated cells. Targeted enrichment increased the yield of B. pertussis sequencing reads up to 90% while simultaneously decreasing host reads to less than 10%. Filtered sequencing reads provided sufficient genome coverage to perform characterization via whole-genome single nucleotide polymorphisms and whole-genome multilocus sequencing typing. Moreover, these data were concordant with sequenced isolates recovered from the same specimens such that phylogenetic reconstructions from either consistently clustered the same putatively linked cases. The optimized protocol is suitable for nasopharyngeal specimens with diagnostic IS481 Ct < 35 and >10 ng DNA. Routine implementation of these methods could strengthen surveillance and study of pertussis resurgence by capturing additional cases with genomic characterization.

A novel hybridization capture method for direct whole genome sequencing of clinical specimens to inform Legionnaires' disease investigations

The unprecedented precision and resolution of whole genome sequencing (WGS) can provide definitive identification of infectious agents for epidemiological outbreak tracking. WGS approaches, however, are frequently impeded by low pathogen DNA recovery from available primary specimens or unculturable samples. A cost-effective hybrid capture assay for Legionella pneumophila WGS analysis directly on primary specimens was developed. DNA from a diverse range of sputum and autopsy specimens PCR-positive for L. pneumophila serogroup 1 (LPSG1) was enriched with this method, and WGS was performed. All tested specimens were determined to be enriched for Legionella reads (up to 209,000-fold), significantly improving the discriminatory power to compare relatedness when no clinical isolate was available. We found the WGS data from some enriched specimens to differ by less than five single-nucleotide polymorphisms (SNPs) when compared to the WGS data of a matched culture isolate. This testing and analysis retrospectively provided previously unconfirmed links to environmental sources for clinical specimens of sputum and autopsy lung tissue. The latter provided the additional information needed to identify the source of these culture-negative cases associated with the South Bronx 2015 Legionnaires' disease (LD) investigation in New York City. This new method provides a proof of concept for future direct clinical specimen hybrid capture enrichment combined with WGS and bioinformatic analysis during outbreak investigations.IMPORTANCELegionnaires' disease (LD) is a severe and potentially fatal type of pneumonia primarily caused by inhalation of Legionella-contaminated aerosols from man-made water or cooling systems. LD remains extremely underdiagnosed as it is an uncommon form of pneumonia and relies on clinicians including it in the differential and requesting specialized testing. Additionally, it is challenging to obtain clinical lower respiratory specimens from cases with LD, and when available, culture requires specialized media and growth conditions, which are not available in all microbiology laboratories. In the current study, a method for Legionella pneumophila using hybrid capture by RNA baiting was developed, which allowed us to generate sufficient genome resolution from L. pneumophila serogroup 1 PCR-positive clinical specimens. This new approach offers an additional tool for surveillance of future LD outbreaks where isolation of Legionella is not possible and may help solve previously unanswered questions from past LD investigations.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

Improved metagenome assemblies through selective enrichment of bacterial genomic DNA from eukaryotic host genomic DNA using ATAC-seq

Genomics can be used to study the complex relationships between hosts and their microbiota. Many bacteria cannot be cultured in the laboratory, making it difficult to obtain adequate amounts of bacterial DNA and to limit host DNA contamination for the construction of metagenome-assembled genomes (MAGs). For example, Wolbachia is a genus of exclusively obligate intracellular bacteria that live in a wide range of arthropods and some nematodes. While Wolbachia endosymbionts are frequently described as facultative reproductive parasites in arthropods, the bacteria are obligate mutualistic endosymbionts of filarial worms. Here, we achieve 50-fold enrichment of bacterial sequences using ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) with Brugia malayi nematodes, containing Wolbachia (wBm). ATAC-seq uses the Tn5 transposase to cut and attach Illumina sequencing adapters to accessible DNA lacking histones, typically thought to be open chromatin. Bacterial and mitochondrial DNA in the lysates are also cut preferentially since they lack histones, leading to the enrichment of these sequences. The benefits of this include minimal tissue input (<1 mg of tissue), a quick protocol (<4 h), low sequencing costs, less bias, correct assembly of lateral gene transfers and no prior sequence knowledge required. We assembled the wBm genome with as few as 1 million Illumina short paired-end reads with >97% coverage of the published genome, compared to only 12% coverage with the standard gDNA libraries. We found significant bacterial sequence enrichment that facilitated genome assembly in previously published ATAC-seq data sets from human cells infected with Mycobacterium tuberculosis and C. elegans contaminated with their food source, the OP50 strain of E. coli. These results demonstrate the feasibility and benefits of using ATAC-seq to easily obtain bacterial genomes to aid in symbiosis, infectious disease, and microbiome research.

Selective whole-genome sequencing of Plasmodium parasites directly from blood samples by nanopore adaptive sampling

IMPORTANCE

Malaria is caused by parasites of the genus Plasmodium, and reached a global disease burden of 247 million cases in 2021. To study drug resistance mutations and parasite population dynamics, whole-genome sequencing of patient blood samples is commonly performed. However, the predominance of human DNA in these samples imposes the need for time-consuming laboratory procedures to enrich Plasmodium DNA. We used the Oxford Nanopore Technologies' adaptive sampling feature to circumvent this problem and enrich Plasmodium reads directly during the sequencing run. We demonstrate that adaptive nanopore sequencing efficiently enriches Plasmodium reads, which simplifies and shortens the timeline from blood collection to parasite sequencing. In addition, we show that the obtained data can be used for monitoring genetic markers, or to generate nearly complete genomes. Finally, owing to its inherent mobility, this technology can be easily applied on-site in endemic areas where patients would benefit the most from genomic surveillance.

Highly multiplexed targeted sequencing strategy for infectious disease surveillance

BACKGROUND

Global efforts to characterize diseases of poverty are hampered by lack of affordable and comprehensive detection platforms, resulting in suboptimal allocation of health care resources and inefficient disease control. Next generation sequencing (NGS) can provide accurate data and high throughput. However, shotgun and metagenome-based NGS approaches are limited by low concentrations of microbial DNA in clinical samples, requirements for tailored sample and library preparations plus extensive bioinformatics analysis. Here, we adapted molecular inversion probes (MIPs) as a cost-effective target enrichment approach to characterize microbial infections from blood samples using short-read sequencing. We designed a probe panel targeting 2 bacterial genera, 21 bacterial and 6 fungi species and 7 antimicrobial resistance markers (AMRs).

RESULTS

Our approach proved to be highly specific to detect down to 1 in a 1000 pathogen DNA targets contained in host DNA. Additionally, we were able to accurately survey pathogens and AMRs in 20 out of 24 samples previously profiled with routine blood culture for sepsis.

CONCLUSIONS

Overall, our targeted assay identifies microbial pathogens and AMRs with high specificity at high throughput, without the need for extensive sample preparation or bioinformatics analysis, simplifying its application for characterization and surveillance of infectious diseases in medium- to low- resource settings.

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.

Genome Capture Sequencing Selectively Enriches Bacterial DNA and Enables Genome-Wide Measurement of Intrastrain Genetic Diversity in Human Infections

Within-host evolution produces genetic diversity in bacterial strains that cause chronic human infections. However, the lack of facile methods to measure bacterial allelic variation in clinical samples has limited understanding of intrastrain diversity’s effects on disease. Here, we report a new method termed genome capture sequencing (GenCap-Seq) in which users inexpensively make hybridization probes from genomic DNA or PCR amplicons to selectively enrich and sequence targeted bacterial DNA from clinical samples containing abundant human or nontarget bacterial DNA. GenCap-Seq enables accurate measurement of allele frequencies over targeted regions and is scalable from specific genes to entire genomes, including the strain-specific accessory genome. The method is effective with samples in which target DNA is rare and inhibitory and DNA-degrading substances are abundant, including human sputum and feces. In proof-of-principle experiments, we used GenCap-Seq to investigate the responses of diversified Pseudomonas aeruginosa populations chronically infecting the lungs of people with cystic fibrosis to in vivo antibiotic exposure, and we found that treatment consistently reduced intrastrain genomic diversity. In addition, analysis of gene-level allele frequency changes suggested that some genes without conventional resistance functions may be important for bacterial fitness during in vivo antibiotic exposure. GenCap-Seq’s ability to scalably enrich targeted bacterial DNA from complex samples will enable studies on the effects of intrastrain and intraspecies diversity in human infectious disease.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

High-throughput method for the hybridisation-based targeted enrichment of long genomic fragments for PacBio third-generation sequencing

Hybridisation-based targeted enrichment is a widely used and well-established technique in high-throughput second-generation short-read sequencing. Despite the high potential to genetically resolve highly repetitive and variable genomic sequences by, for example PacBio third-generation sequencing, targeted enrichment for long fragments has not yet established the same high-throughput due to currently existing complex workflows and technological dependencies. We here describe a scalable targeted enrichment protocol for fragment sizes of >7 kb. For demonstration purposes we developed a custom blood group panel of challenging loci. Test results achieved > 65% on-target rate, good coverage (142.7×) and sufficient coverage evenness for both non-paralogous and paralogous targets, and sufficient non-duplicate read counts (83.5%) per sample for a highly multiplexed enrichment pool of 16 samples. We genotyped the blood groups of nine patients employing highly accurate phased assemblies at an allelic resolution that match reference blood group allele calls determined by SNP array and NGS genotyping. Seven Genome-in-a-Bottle reference samples achieved high recall (96%) and precision (99%) rates. Mendelian error rates were 0.04% and 0.13% for the included Ashkenazim and Han Chinese trios, respectively. In summary, we provide a protocol and first example for accurate targeted long-read sequencing that can be used in a high-throughput fashion.

Probe Capture Enrichment Methods for HIV and HCV Genome Sequencing and Drug Resistance Genotyping

Human immunodeficiency virus (HIV) infections remain a significant public health concern worldwide. Over the years, sophisticated sequencing technologies such as next-generation sequencing (NGS) have emerged and been utilized to monitor the spread of HIV drug resistance (HIVDR), identify HIV drug resistance mutations, and characterize transmission dynamics. Similar applications also apply to the Hepatitis C virus (HCV), another bloodborne viral pathogen with significant intra-host genetic diversity. Several advantages to using NGS over conventional Sanger sequencing include increased data throughput, scalability, cost-effectiveness when batched sample testing is performed, and sensitivity for quantitative detection of minority resistant variants. However, NGS alone may fail to detect genomes from pathogens present in low copy numbers. As with all sequencing platforms, the primary determinant in achieving quality sequencing data is the quality and quantity of the initial template input. Samples containing degraded RNA/DNA and/or low copy number have been a consistent sequencing challenge. To overcome this limitation probe capture enrichment is a method that has recently been employed to target, enrich, and sequence the genome of a pathogen present in low copies, and for compromised specimens that contain poor quality nucleic acids. It involves the hybridization of sequence-specific DNA or RNA probes to a target sequence, which is followed by an enrichment step via PCR to increase the number of copies of the targeted sequences after which the samples are subjected to NGS procedures. This method has been performed on pathogens such as bacteria, fungus, and viruses and allows for the sequencing of complete genomes, with high coverage. Post NGS, data analysis can be performed through various bioinformatics pipelines which can provide information on genetic diversity, genotype, virulence, and drug resistance. This article reviews how probe capture enrichment helps to increase the likelihood of sequencing HIV and HCV samples that contain low viral loads and/or are compromised.

Trends and Developments in the Detection of Pathogens in Central Nervous System Infections: A Bibliometric Study

Introduction

Rapid, sensitive, and specific laboratory assays are critical for the diagnosis and management of central nervous system (CNS) infections. The purpose of this study is to explore the intellectual landscape of research investigating methods for the detection of pathogens in patients with CNS infections and to identify the development trends and research frontier in this field.

Methods

A bibliometric study is conducted by analyzing literature retrieved from the Web of Science (WoS) Core Collection Database for the years 2000 to 2021. CiteSpace software is used for bibliometric analysis and network visualization, including co-citation analysis of references, co-occurrence analysis of keywords, and cooperation network analysis of authors, institutions, and countries/regions.

Results

A total of 2,282 publications are eventually screened, with an upward trend in the number of publications per year. The majority of papers are attributed to the disciplines of MICROBIOLOGY, INFECTIOUS DISEASES, IMMUNOLOGY, NEUROSCIENCES & NEUROLOGY, and VIROLOGY. The co-citation analysis of references shows that recent research has focused on the largest cluster “metagenomic next-generation sequencing”; the results of the analysis of the highest-cited publications and the citation burst of publications reveal that there is a strong interest stimulated in metagenomic next-generation sequencing. The co-occurrence analysis of keywords indicates that “infection”, “pathogen”, “diagnosis”, “gene”, “virus”, “polymerase chain reaction”, “cerebrospinal fluid”, “epidemiology”, and “metagenomic next-generation sequencing” are the main research priorities in the field of pathogen detection for CNS infections, and the keyword with the highest strength of burst is “metagenomic next-generation sequencing”. Collaborative network analysis reveals that the USA, the Centers for Disease Control and Prevention of USA, and XIN WANG and JENNIFER DIEN BARD are the most influential country, institution, and researchers, respectively.

Conclusions

Exploring more advanced laboratory assays to improve the diagnostic accuracy of pathogens is essential for CNS infection research. Metagenomic next-generation sequencing is emerging as a novel useful unbiased approach for diagnosing infectious diseases of the CNS.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Challenges for Cryptosporidium Population Studies

Cryptosporidiosis is ranked sixth in the list of the most important food-borne parasites globally, and it is an important contributor to mortality in infants and the immunosuppressed. Recently, the number of genome sequences available for this parasite has increased drastically. The majority of the sequences are derived from population studies of Cryptosporidium parvum and Cryptosporidium hominis, the most important species causing disease in humans. Work with this parasite is challenging since it lacks an optimal, prolonged, in vitro culture system, which accurately reproduces the in vivo life cycle. This obstacle makes the cloning of isolates nearly impossible. Thus, patient isolates that are sequenced represent a population or, at times, mixed infections. Oocysts, the lifecycle stage currently used for sequencing, must be considered a population even if the sequence is derived from single-cell sequencing of a single oocyst because each oocyst contains four haploid meiotic progeny (sporozoites). Additionally, the community does not yet have a set of universal markers for strain typing that are distributed across all chromosomes. These variables pose challenges for population studies and require careful analyses to avoid biased interpretation. This review presents an overview of existing population studies, challenges, and potential solutions to facilitate future population analyses.

Single-cell sequencing of the small and AT-skewed genome of malaria parasites

Single-cell genomics is a rapidly advancing field; however, most techniques are designed for mammalian cells. We present a single-cell sequencing pipeline for an intracellular parasite, Plasmodium falciparum, with a small genome of extreme base content. Through optimization of a quasi-linear amplification method, we target the parasite genome over contaminants and generate coverage levels allowing detection of minor genetic variants. This work, as well as efforts that build on these findings, will enable detection of parasite heterogeneity contributing to P. falciparum adaptation. Furthermore, this study provides a framework for optimizing single-cell amplification and variant analysis in challenging genomes.

The molecular basis of antimalarial drug resistance in Plasmodium vivax

Plasmodium vivax is the most geographically widespread cause of human malaria and is responsible for the majority of cases outside of the African continent. While great progress has been made towards eliminating human malaria, drug resistant parasite strains pose a threat towards continued progress. Resistance has arisen to multiple antimalarials in P. vivax, including to chloroquine, which is currently the first line therapy for P. vivax in most regions. Despite its importance, an understanding of the molecular mechanisms of drug resistance in this species remains elusive, in large part due to the complex biology of P. vivax and the lack of in vitro culture. In this review, we will cover the extent and challenges of measuring clinical and in vitro drug resistance in P. vivax. We will consider the roles of candidate drug resistance genes. We will highlight the development of molecular approaches for studying P. vivax biology that provide the opportunity to validate the role of putative drug resistance mutations as well as identify novel mechanisms of drug resistance in this understudied parasite. Validated molecular determinants and markers of drug resistance are essential for the rapid and cost-effective monitoring of drug resistance in P. vivax, and will be useful for optimizing drug regimens and for informing drug policy in control and elimination settings.

•

Drug resistance is emerging in Plasmodium vivax, an important cause of malaria.

•

The complex biology of P. vivax and the limited range of research tools make it difficult to identify drug resistance.

•

The molecular mechanisms of drug resistance in P. vivax remain elusive.

•

This review highlights the extent of drug resistance, the putative mechanisms of resistance and new technologies for the study of P. vivax drug resistance.

Population genomic evidence of Plasmodium vivax Southeast Asian origin

Plasmodium vivax is the most common and widespread human malaria parasite. It was recently proposed that P. vivax originates from sub-Saharan Africa based on the circulation of its closest genetic relatives (P. vivax-like) among African great apes. However, the limited number of genetic markers and samples investigated questions the robustness of this hypothesis. Here, we extensively characterized the genomic variations of 447 human P. vivax strains and 19 ape P. vivax-like strains collected worldwide. Phylogenetic relationships between human and ape Plasmodium strains revealed that P. vivax is a sister clade of P. vivax-like, not included within the radiation of P. vivax-like By investigating various aspects of P. vivax genetic variation, we identified several notable geographical patterns in summary statistics in function of the increasing geographic distance from Southeast Asia, suggesting that P. vivax may have derived from a single area in Asia through serial founder effects.

Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples

Background

Whole-genome sequencing (WGS) is becoming increasingly useful to study the biology, epidemiology, and ecology of malaria parasites. Despite ease of sampling, DNA extracted from dried blood spots (DBS) has a high ratio of human DNA compared to parasite DNA, which poses a challenge for downstream genetic analyses. The effects of multiple methods for DNA extraction, digestion of methylated DNA, and amplification were evaluated on the quality and fidelity of WGS data recovered from DBS.

Methods

Low parasite density mock DBS samples were created, extracted either with Tween-Chelex or QIAamp, treated with or without McrBC, and amplified with one of three different amplification techniques (two sWGA primer sets and one rWGA). Extraction conditions were evaluated on performance of sequencing depth, percentiles of coverage, and expected SNP concordance.

Results

At 100 parasites/μL, Chelex-Tween-McrBC samples had higher coverage (5 × depth = 93% genome) than QIAamp extracted samples (5 × depth = 76% genome). The two evaluated sWGA primer sets showed minor differences in overall genome coverage and SNP concordance, with a newly proposed combination of 20 primers showing a modest improvement in coverage over those previously published.

Conclusions

Overall, Tween-Chelex extracted samples that were treated with McrBC digestion and are amplified using 6A10AD sWGA conditions had minimal dropout rate, higher percentages of coverage at higher depth, and more accurate SNP concordance than QiaAMP extracted samples. These findings extend the results of previously reported methods, making whole genome sequencing accessible to a larger number of low density samples that are commonly encountered in cross-sectional surveys.

Next-generation diagnostics: virus capture facilitates a sensitive viral diagnosis for epizootic and zoonotic pathogens including SARS-CoV-2

BACKGROUND

The detection of pathogens in clinical and environmental samples using high-throughput sequencing (HTS) is often hampered by large amounts of background information, which is especially true for viruses with small genomes. Enormous sequencing depth can be necessary to compile sufficient information for identification of a certain pathogen. Generic HTS combining with in-solution capture enrichment can markedly increase the sensitivity for virus detection in complex diagnostic samples.

METHODS

A virus panel based on the principle of biotinylated RNA baits was developed for specific capture enrichment of epizootic and zoonotic viruses (VirBaits). The VirBaits set was supplemented by a SARS-CoV-2 predesigned bait set for testing recent SARS-CoV-2-positive samples. Libraries generated from complex samples were sequenced via generic HTS (without enrichment) and afterwards enriched with the VirBaits set. For validation, an internal proficiency test for emerging epizootic and zoonotic viruses (African swine fever virus, Ebolavirus, Marburgvirus, Nipah henipavirus, Rift Valley fever virus) was conducted.

RESULTS

The VirBaits set consists of 177,471 RNA baits (80-mer) based on about 18,800 complete viral genomes targeting 35 epizootic and zoonotic viruses. In all tested samples, viruses with both DNA and RNA genomes were clearly enriched ranging from about 10-fold to 10,000-fold for viruses including distantly related viruses with at least 72% overall identity to viruses represented in the bait set. Viruses showing a lower overall identity (38% and 46%) to them were not enriched but could nonetheless be detected based on capturing conserved genome regions. The internal proficiency test supports the improved virus detection using the combination of HTS plus targeted enrichment but also points to the risk of cross-contamination between samples.

CONCLUSIONS

The VirBaits approach showed a high diagnostic performance, also for distantly related viruses. The bait set is modular and expandable according to the favored diagnostics, health sector, or research question. The risk of cross-contamination needs to be taken into consideration. The application of the RNA-baits principle turned out to be user friendly, and even non-experts can easily use the VirBaits workflow. The rapid extension of the established VirBaits set adapted to actual outbreak events is possible as shown for SARS-CoV-2. Video abstract.

Next-generation diagnostics: virus capture facilitates a sensitive viral diagnosis for epizootic and zoonotic pathogens including SARS-CoV-2

BACKGROUND

The detection of pathogens in clinical and environmental samples using high-throughput sequencing (HTS) is often hampered by large amounts of background information, which is especially true for viruses with small genomes. Enormous sequencing depth can be necessary to compile sufficient information for identification of a certain pathogen. Generic HTS combining with in-solution capture enrichment can markedly increase the sensitivity for virus detection in complex diagnostic samples.

METHODS

A virus panel based on the principle of biotinylated RNA baits was developed for specific capture enrichment of epizootic and zoonotic viruses (VirBaits). The VirBaits set was supplemented by a SARS-CoV-2 predesigned bait set for testing recent SARS-CoV-2-positive samples. Libraries generated from complex samples were sequenced via generic HTS (without enrichment) and afterwards enriched with the VirBaits set. For validation, an internal proficiency test for emerging epizootic and zoonotic viruses (African swine fever virus, Ebolavirus, Marburgvirus, Nipah henipavirus, Rift Valley fever virus) was conducted.

RESULTS

The VirBaits set consists of 177,471 RNA baits (80-mer) based on about 18,800 complete viral genomes targeting 35 epizootic and zoonotic viruses. In all tested samples, viruses with both DNA and RNA genomes were clearly enriched ranging from about 10-fold to 10,000-fold for viruses including distantly related viruses with at least 72% overall identity to viruses represented in the bait set. Viruses showing a lower overall identity (38% and 46%) to them were not enriched but could nonetheless be detected based on capturing conserved genome regions. The internal proficiency test supports the improved virus detection using the combination of HTS plus targeted enrichment but also points to the risk of cross-contamination between samples.

CONCLUSIONS

The VirBaits approach showed a high diagnostic performance, also for distantly related viruses. The bait set is modular and expandable according to the favored diagnostics, health sector, or research question. The risk of cross-contamination needs to be taken into consideration. The application of the RNA-baits principle turned out to be user friendly, and even non-experts can easily use the VirBaits workflow. The rapid extension of the established VirBaits set adapted to actual outbreak events is possible as shown for SARS-CoV-2. Video abstract.

Whole genome sequencing of Theileria parva using target capture

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0-33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva.

Culture-free genome-wide locus sequence typing (GLST) provides new perspectives on Trypanosoma cruzi dispersal and infection complexity

Analysis of genetic polymorphism is a powerful tool for epidemiological surveillance and research. Powerful inference from pathogen genetic variation, however, is often restrained by limited access to representative target DNA, especially in the study of obligate parasitic species for which ex vivo culture is resource-intensive or bias-prone. Modern sequence capture methods enable pathogen genetic variation to be analyzed directly from host/vector material but are often too complex and expensive for resource-poor settings where infectious diseases prevail. This study proposes a simple, cost-effective 'genome-wide locus sequence typing' (GLST) tool based on massive parallel amplification of information hotspots throughout the target pathogen genome. The multiplexed polymerase chain reaction amplifies hundreds of different, user-defined genetic targets in a single reaction tube, and subsequent agarose gel-based clean-up and barcoding completes library preparation at under 4 USD per sample. Our study generates a flexible GLST primer panel design workflow for Trypanosoma cruzi, the parasitic agent of Chagas disease. We successfully apply our 203-target GLST panel to direct, culture-free metagenomic extracts from triatomine vectors containing a minimum of 3.69 pg/μl T. cruzi DNA and further elaborate on method performance by sequencing GLST libraries from T. cruzi reference clones representing discrete typing units (DTUs) TcI, TcIII, TcIV, TcV and TcVI. The 780 SNP sites we identify in the sample set repeatably distinguish parasites infecting sympatric vectors and detect correlations between genetic and geographic distances at regional (< 150 km) as well as continental scales. The markers also clearly separate TcI, TcIII, TcIV and TcV + TcVI and appear to distinguish multiclonal infections within TcI. We discuss the advantages, limitations and prospects of our method across a spectrum of epidemiological research.

Capture-based enrichment of Theileria parva DNA enables full genome assembly of first buffalo-derived strain and reveals exceptional intra-specific genetic diversity

Theileria parva is an economically important, intracellular, tick-transmitted parasite of cattle. A live vaccine against the parasite is effective against challenge from cattle-transmissible T. parva but not against genotypes originating from the African Cape buffalo, a major wildlife reservoir, prompting the need to characterize genome-wide variation within and between cattle- and buffalo-associated T. parva populations. Here, we describe a capture-based target enrichment approach that enables, for the first time, de novo assembly of nearly complete T. parva genomes derived from infected host cell lines. This approach has exceptionally high specificity and sensitivity and is successful for both cattle- and buffalo-derived T. parva parasites. De novo genome assemblies generated for cattle genotypes differ from the reference by ~54K single nucleotide polymorphisms (SNPs) throughout the 8.31 Mb genome, an average of 6.5 SNPs/kb. We report the first buffalo-derived T. parva genome, which is ~20 kb larger than the genome from the reference, cattle-derived, Muguga strain, and contains 25 new potential genes. The average non-synonymous nucleotide diversity (π_N) per gene, between buffalo-derived T. parva and the Muguga strain, was 1.3%. This remarkably high level of genetic divergence is supported by an average Wright’s fixation index (F_ST), genome-wide, of 0.44, reflecting a degree of genetic differentiation between cattle- and buffalo-derived T. parva parasites more commonly seen between, rather than within, species. These findings present clear implications for vaccine development, further demonstrated by the ability to assemble nearly all known antigens in the buffalo-derived strain, which will be critical in design of next generation vaccines. The DNA capture approach used provides a clear advantage in specificity over alternative T. parva DNA enrichment methods used previously, such as those that utilize schizont purification, is less labor intensive, and enables in-depth comparative genomics in this apicomplexan parasite.

An estimated 50 million cattle in sub-Saharan Africa are at risk of the deadly livestock disease East coast fever (ECF), caused by the parasite Theileria parva, which imposes tremendous economic hardship on smallholder farmers. An existing ECF vaccine protects against strains circulating among cattle, but not against T. parva derived from African Cape buffalo, its main wildlife carrier. Understanding this difference in protective efficacy requires characterization of the genetic diversity in T. parva strains associated with each mammalian host, a goal that has been hindered by the proliferation of T. parva in nucleated host cells, with much larger genomes. Here we adapted a sequence capture approach to target the whole parasite genome, enabling enrichment of parasite DNA over that of the host. Choices in protocol development resulted in nearly 100% parasite genome specificity and sensitivity, making this approach the most successful yet to generate T. parva genome sequence data in a high-throughput manner. The analyses uncovered a degree of genetic differentiation between cattle- and buffalo-derived genotypes that is akin to levels more commonly seen between species. This approach, which will enable an in-depth T. parva population genomics study from cattle and buffalo in the endemic regions, can easily be adapted to other intracellular pathogens.

Ultraconserved element bait set for trypanosomatida target enrichment and phylogenetics

Lack of knowledge of taxonomic biodiversity and reliable genetic markers in Trypanosomatidae limit our understanding of their phylogenetic relationships. Ultraconserved elements (UCEs) have improved phylogenetic analyses and inferences in many vertebrate and invertebrate taxa. However, it is unknown whether protozoans have these markers, their abundance, and if these could be reliably used for phylogenetics. In this study I design a target enrichment bait set for UCE loci for this group. In silico testing showed good loci recovery rates across 63 taxa and produced consistent, highly supported phylogenetic trees. This bait set adds a new resource of useful genetic markers for Trypanosomatidae phylogenetics.

Human Cytomegalovirus Genomes Sequenced Directly From Clinical Material: Variation, Multiple-Strain Infection, Recombination, and Gene Loss

The genomic characteristics of human cytomegalovirus (HCMV) strains sequenced directly from clinical pathology samples were investigated, focusing on variation, multiple-strain infection, recombination, and gene loss. A total of 207 datasets generated in this and previous studies using target enrichment and high-throughput sequencing were analyzed, in the process enabling the determination of genome sequences for 91 strains. Key findings were that (i) it is important to monitor the quality of sequencing libraries in investigating variation; (ii) many recombinant strains have been transmitted during HCMV evolution, and some have apparently survived for thousands of years without further recombination; (iii) mutants with nonfunctional genes (pseudogenes) have been circulating and recombining for long periods and can cause congenital infection and resulting clinical sequelae; and (iv) intrahost variation in single-strain infections is much less than that in multiple-strain infections. Future population-based studies are likely to continue illuminating the evolution, epidemiology, and pathogenesis of HCMV.

Optimization of parasite DNA enrichment approaches to generate whole genome sequencing data for Plasmodium falciparum from low parasitaemia samples

Background

Owing to the large amount of host DNA in clinical samples, generation of high-quality Plasmodium falciparum whole genome sequencing (WGS) data requires enrichment for parasite DNA. Enrichment is often achieved by leukocyte depletion of infected blood prior to storage. However, leukocyte depletion is difficult in low-resource settings and limits analysis to prospectively-collected samples. As a result, approaches such as selective whole genome amplification (sWGA) are being used to enrich for parasite DNA. However, sWGA has had limited success in generating reliable sequencing data from low parasitaemia samples. In this study, enzymatic digestion with MspJI prior to sWGA and whole genome sequencing was evaluated to determine whether this approach improved genome coverage compared to sWGA alone. The potential of sWGA to cause amplification bias in polyclonal infections was also examined.

Methods

DNA extracted from laboratory-created dried blood spots was treated with a modification-dependent restriction endonuclease, MspJI, and filtered via vacuum filtration. Samples were then selectively amplified using a previously reported sWGA protocol and subjected to WGS. Genome coverage statistics were compared between the optimized sWGA approach and the previously reported sWGA approach performed in parallel. Differential amplification by sWGA was assessed by comparing WGS data generated from lab-created mixtures of parasite isolates, from the same geographical region, generated with or without sWGA.

Results

MspJI digestion did not enrich for parasite DNA. Samples that underwent vacuum filtration (without MspJI digestion) prior to sWGA had the highest parasite DNA concentration and displayed greater genome coverage compared to MspJI + sWGA and sWGA alone, particularly for low parasitaemia samples. The optimized sWGA (filtration + sWGA) approach was successfully used to generate WGS data from 218 non-leukocyte depleted field samples from Malawi. Sequences from lab-created mixtures of parasites did not show evidence of differential amplification of parasite strains compared to directly sequenced samples.

Conclusion

This optimized sWGA approach is a reliable method to obtain WGS data from non-leukocyte depleted, low parasitaemia samples. The absence of amplification bias in data generated from mixtures of isolates from the same geographic region suggests that this approach can be appropriately used for molecular epidemiological studies.

Spatially distinct physiology of Bacteroides fragilis within the proximal colon of gnotobiotic mice

A complex microbiota inhabits various microenvironments of the gut, with some symbiotic bacteria having evolved traits to invade the epithelial mucus layer and reside deep within intestinal tissue of animals. Whether these distinct bacterial communities across gut biogeographies exhibit divergent behaviors remains largely unknown. Global transcriptomic analysis to investigate microbial physiology in specific mucosal niches has been hampered technically by overabundance of host RNA. Herein, we employed hybrid selection RNA sequencing (hsRNA-Seq) to enable detailed spatial transcriptomic profiling of a prominent human commensal as it colonizes the colonic lumen, mucus or epithelial tissue of mice. Compared to conventional RNA-Seq, hsRNA-Seq increased reads mapping to the Bacteroides fragilis genome by 48- and 154-fold in mucus and tissue, respectively, allowing for high fidelity comparisons across biogeographic sites. Near the epithelium, B. fragilis up-regulated numerous genes involved in protein synthesis, indicating that bacteria inhabiting the mucosal niche are metabolically active. Further, a specific sulfatase (BF3086) and glycosyl hydrolase (BF3134) were highly induced in mucus and tissue compared to bacteria in the lumen. In-frame deletion of these genes impaired in vitro growth on mucus as a carbon source, as well as mucosal colonization of mice. Mutants in either B. fragilis gene displayed a fitness defect in competing for colonization against bacterial challenge, revealing the importance of site-specific gene expression for robust host-microbial symbiosis. As a versatile tool, hsRNA-Seq can be deployed to explore the in vivo spatial physiology of numerous bacterial pathogens or commensals.

Reduction of Human DNA Contamination in Clinical Cerebrospinal Fluid Specimens Improves the Sensitivity of Metagenomic Next-Generation Sequencing

Metagenomics next-generation sequencing (mNGS) is increasingly available for the detection of obscure infectious diseases of the central nervous system. However, human DNA contamination from elevated white cells, one of the characteristic cerebrospinal fluid (CSF) features in meningitis patients, greatly reduces the sensitivity of mNGS in the pathogen detection. Currently, effective approaches to selectively reduce host DNA contamination from clinical CSF samples are still lacking. In this study, a total of 20 meningitis patients were enrolled, including 10 definitively diagnosed tuberculous meningitis (TBM) and 10 definite cryptococcal meningitis (CM) cases. To evaluate the effect of reduced human DNA in the sensitivity of mNGS detection, three specimen-processing protocols were performed: (i) To remove human DNA, saponin, a nonionic surfactant, was used to selectively lyse white cells in CSF followed by DNase treatment prior to the extraction of DNA; (ii) to reduce host DNA, CSF was centrifuged to remove human cells, and the supernatant was collected for DNA extraction; and (iii) DNA extraction from the unprocessed specimens was set as the control. We found that saponin processing significantly elevated the NGS unique reads for Cryptococcus (P < 0.01) compared with the control but had no effects for Mycobacterium tuberculosis (P > 0.05). However, detection of centrifuged supernatants improved the NGS unique reads for both TBM and CM compared with controls (P < 0.01). Our results demonstrate that the use of mNGS of centrifuged supernatants from clinical CSF samples in patients with TBM and CM is a simple and effective method to improve the sensitivity of pathogen detection.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

Plasmodium vivax chloroquine resistance links to pvcrt transcription in a genetic cross

Mainstay treatment for Plasmodium vivax malaria has long relied on chloroquine (CQ) against blood-stage parasites plus primaquine against dormant liver-stage forms (hypnozoites), however drug resistance confronts this regimen and threatens malaria control programs. Understanding the basis of P. vivax chloroquine resistance (CQR) will inform drug discovery and malaria control. Here we investigate the genetics of P. vivax CQR by a cross of parasites differing in drug response. Gametocytogenesis, mosquito infection, and progeny production are performed with mixed parasite populations in nonhuman primates, as methods for P. vivax cloning and in vitro cultivation remain unavailable. Linkage mapping of progeny surviving >15 mg/kg CQ identifies a 76 kb region in chromosome 1 including pvcrt, an ortholog of the Plasmodium falciparum CQR transporter gene. Transcriptional analysis supports upregulated pvcrt expression as a mechanism of CQR.

Whole-Genome Single-Nucleotide Polymorphism (SNP) Analysis Applied Directly to Stool for Genotyping Shiga Toxin-Producing Escherichia coli: an Advanced Molecular Detection Method for Foodborne Disease Surveillance and Outbreak Tracking

Whole-genome sequencing (WGS) of pathogens from pure culture provides unparalleled accuracy and comprehensive results at a cost that is advantageous compared with traditional diagnostic methods. Sequencing pathogens directly from a primary clinical specimen would help circumvent the need for culture and, in the process, substantially shorten the time to diagnosis and public health reporting. Unfortunately, this approach poses significant challenges because of the mixture of multiple sequences from a complex fecal biomass. The aim of this project was to develop a proof of concept protocol for the sequencing and genotyping of Shiga toxin-producing Escherichia coli (STEC) directly from stool specimens. We have developed an enrichment protocol that reliably achieves a substantially higher DNA yield belonging to E. coli, which provides adequate next-generation sequencing (NGS) data for downstream bioinformatics analysis. A custom bioinformatics pipeline was created to optimize and remove non-E. coli reads, assess the STEC versus commensal E. coli population in the samples, and build consensus sequences based on population allele frequency distributions. Side-by-side analysis of WGS from paired STEC isolates and matched primary stool specimens reveal that this method can reliably be implemented for many clinical specimens to directly genotype STEC and accurately identify clusters of disease outbreak when no STEC isolate is available for testing.

Capturing sequence diversity in metagenomes with comprehensive and scalable probe design

Metagenomic sequencing has the potential to transform microbial detection and characterization, but new tools are needed to improve its sensitivity. Here we present CATCH, a computational method to enhance nucleic acid capture for enrichment of diverse microbial taxa. CATCH designs optimal probe sets, with a specified number of oligonucleotides, that achieve full coverage of, and scale well with, known sequence diversity. We focus on applying CATCH to capture viral genomes in complex metagenomic samples. We design, synthesize, and validate multiple probe sets, including one that targets the whole genomes of the 356 viral species known to infect humans. Capture with these probe sets enriches unique viral content on average 18-fold, allowing us to assemble genomes that could not be recovered without enrichment, and accurately preserves within-sample diversity. We also use these probe sets to recover genomes from the 2018 Lassa fever outbreak in Nigeria and to improve detection of uncharacterized viral infections in human and mosquito samples. The results demonstrate that CATCH enables more sensitive and cost-effective metagenomic sequencing.

Hybrid Capture-Based Next Generation Sequencing and Its Application to Human Infectious Diseases

This review describes target-enrichment approaches followed by next generation sequencing and their recent application to the research and diagnostic field of modern and past infectious human diseases caused by viruses, bacteria, parasites and fungi.

RNA-Seq Profile Reveals Th-1 and Th-17-Type of Immune Responses in Mice Infected Systemically with Aspergillus fumigatus

With the increasing numbers of immunocompromised hosts, Aspergillus fumigatus emerges as a lethal opportunistic fungal pathogen. Understanding innate and acquired immunity responses of the host is important for a better therapeutic strategy to deal with aspergillosis patients. To determine the transcriptome in the kidneys in aspergillosis, we employed RNA-Seq to obtain single 76-base reads of whole-genome transcripts of murine kidneys on a temporal basis (days 0; uninfected, 1, 2, 3 and 8) during invasive aspergillosis. A total of 6284 transcripts were downregulated, and 5602 were upregulated compared to baseline expression. Gene ontology enrichment analysis identified genes involved in innate and adaptive immune response, as well as iron binding and homeostasis, among others. Our results showed activation of pathogen recognition receptors, e.g., β-defensins, C-type lectins (e.g., dectin-1), Toll-like receptors (TLR-2, TLR-3, TLR-8, TLR-9 and TLR-13), as well as Ptx-3 and C-reactive protein among the soluble receptors. Upregulated transcripts encoding various differentiating cytokines and effector proinflammatory cytokines, as well as those encoding for chemokines and chemokine receptors, revealed Th-1 and Th-17-type immune responses. These studies form a basic dataset for experimental prioritization, including other target organs, to determine the global response of the host against Aspergillus infection.

Exploring DNA quality of single cells for genome analysis with simultaneous whole-genome amplification

Single cell genome analysis methods are powerful tools to define features of single cells and to identify differences between them. Since the DNA amount of a single cell is very limited, cellular DNA usually needs to be amplified by whole-genome amplification before being subjected to further analysis. A single nucleus only contains two haploid genomes. Thus, any DNA damage that prevents amplification results in loss of damaged DNA sites and induces an amplification bias. Therefore, the assessment of single cell DNA quality is urgently required. As of today, there is no simple method to determine the quality of a single cell DNA in a manner that will still retain the entire cellular DNA for amplification and downstream analysis. Here, we describe a method for whole-genome amplification with simultaneous quality control of single cell DNA by using a competitive spike-in DNA template.

In Vitro Culture, Drug Sensitivity, and Transcriptome of Plasmodium Vivax Hypnozoites

The unique relapsing nature of Plasmodium vivax infection is a major barrier to malaria eradication. Upon infection, dormant liver-stage forms, hypnozoites, linger for weeks to months and then relapse to cause recurrent blood-stage infection. Very little is known about hypnozoite biology; definitive biomarkers are lacking and in vitro platforms that support phenotypic studies are needed. Here, we recapitulate the entire liver stage of P. vivax in vitro, using a multiwell format that incorporates micropatterned primary human hepatocyte co-cultures (MPCCs). MPCCs feature key aspects of P. vivax biology, including establishment of persistent small forms and growing schizonts, merosome release, and subsequent infection of reticulocytes. We find that the small forms exhibit previously described hallmarks of hypnozoites, and we pilot MPCCs as a tool for testing candidate anti-hypnozoite drugs. Finally, we employ a hybrid capture strategy and RNA sequencing to describe the hypnozoite transcriptome and gain insight into its biology.

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

Methylation-based enrichment facilitates low-cost, noninvasive genomic scale sequencing of populations from feces

Obtaining high-quality samples from wild animals is a major obstacle for genomic studies of many taxa, particularly at the population level, as collection methods for such samples are typically invasive. DNA from feces is easy to obtain noninvasively, but is dominated by bacterial and other non-host DNA. The high proportion of non-host DNA drastically reduces the efficiency of high-throughput sequencing for host animal genomics. To address this issue, we developed an inexpensive capture method for enriching host DNA from noninvasive fecal samples. Our method exploits natural differences in CpG-methylation density between vertebrate and bacterial genomes to preferentially bind and isolate host DNA from majority-bacterial samples. We demonstrate that the enrichment is robust, efficient, and compatible with downstream library preparation methods useful for population studies (e.g., RADseq). Compared to other enrichment strategies, our method is quick and inexpensive, adding only a negligible cost to sample preparation. In combination with downstream methods such as RADseq, our approach allows for cost-effective and customizable genomic-scale genotyping that was previously feasible in practice only with invasive samples. Because feces are widely available and convenient to collect, our method empowers researchers to explore genomic-scale population-level questions in organisms for which invasive sampling is challenging or undesirable.

Culture-independent approaches to chlamydial genomics

The expanding field of bacterial genomics has revolutionized our understanding of microbial diversity, biology and phylogeny. For most species, DNA extracted from culture material is used as the template for genome sequencing; however, the majority of microbes are actually uncultivable, and others, such as obligate intracellular bacteria, require laborious tissue culture to yield sufficient genomic material for sequencing. Chlamydiae are one such group of obligate intracellular microbes whose characterization has been hampered by this requirement. To circumvent these challenges, researchers have developed culture-independent sample preparation methods that can be applied to the sample directly or to genomic material extracted from the sample. These methods, which encompass both targeted [immunomagnetic separation-multiple displacement amplification (IMS-MDA) and sequence capture] and non-targeted approaches (host methylated DNA depletion-microbial DNA enrichment and cell-sorting-MDA), have been applied to a range of clinical and environmental samples to generate whole genomes of novel chlamydial species and strains. This review aims to provide an overview of the application, advantages and limitations of these targeted and non-targeted approaches in the chlamydial context. The methods discussed also have broad application to other obligate intracellular bacteria or clinical and environmental samples.

Early Epstein-Barr Virus Genomic Diversity and Convergence toward the B95.8 Genome in Primary Infection

Over 90% of the world's population is persistently infected with Epstein-Barr virus. While EBV does not cause disease in most individuals, it is the common cause of acute infectious mononucleosis (AIM) and has been associated with several cancers and autoimmune diseases, highlighting a need for a preventive vaccine. At present, very few primary, circulating EBV genomes have been sequenced directly from infected individuals. While low levels of diversity and low viral evolution rates have been predicted for double-stranded DNA (dsDNA) viruses, recent studies have demonstrated appreciable diversity in common dsDNA pathogens (e.g., cytomegalovirus). Here, we report 40 full-length EBV genome sequences obtained from matched oral wash and B cell fractions from a cohort of 10 AIM patients. Both intra- and interpatient diversity were observed across the length of the entire viral genome. Diversity was most pronounced in viral genes required for establishing latent infection and persistence, with appreciable levels of diversity also detected in structural genes, including envelope glycoproteins. Interestingly, intrapatient diversity declined significantly over time (P < 0.01), and this was particularly evident on comparison of viral genomes sequenced from B cell fractions in early primary infection and convalescence (P < 0.001). B cell-associated viral genomes were observed to converge, becoming nearly identical to the B95.8 reference genome over time (Spearman rank-order correlation test; r = -0.5589, P = 0.0264). The reduction in diversity was most marked in the EBV latency genes. In summary, our data suggest independent convergence of diverse viral genome sequences toward a reference-like strain within a relatively short period following primary EBV infection.IMPORTANCE Identification of viral proteins with low variability and high immunogenicity is important for the development of a protective vaccine. Knowledge of genome diversity within circulating viral populations is a key step in this process, as is the expansion of intrahost genomic variation during infection. We report full-length EBV genomes sequenced from the blood and oral wash of 10 individuals early in primary infection and during convalescence. Our data demonstrate considerable diversity within the pool of circulating EBV strains, as well as within individual patients. Overall viral diversity decreased from early to persistent infection, particularly in latently infected B cells, which serve as the viral reservoir. Reduction in B cell-associated viral genome diversity coincided with a convergence toward a reference-like EBV genotype. Greater convergence positively correlated with time after infection, suggesting that the reference-like genome is the result of selection.

Malaria Genomics in the Era of Eradication

The first reference genome assembly for the Plasmodium falciparum malaria parasite was completed over a decade ago, and the impact of this and other genomic resources on malaria research has been significant. Genomic resources for other malaria parasites are being established, even as P. falciparum continues to be the focus of development of new genomic methods and applications. Here we review the impact and applications of genomic data on malaria research, and discuss future needs and directions as genomic data generation becomes less expensive and more decentralized. Specifically, we focus on how population genomic strategies can be utilized to advance the malaria eradication agenda.

Gut Microbiota, Inflammation, and Colorectal Cancer

Colorectal cancer is the second-leading cause of cancer-related deaths in the United States and fourth-leading cause of cancer-related deaths worldwide. While cancer is largely considered to be a disease of genetic and environmental factors, increasing evidence has demonstrated a role for the microbiota (the microorganisms associated with the human body) in shaping inflammatory environments and promoting tumor growth and spread. Herein, we discuss both human data from meta'omics analyses and data from mechanistic studies in cell culture and animal models that support specific bacterial agents as potentiators of tumorigenesis-including Fusobacterium nucleatum, enterotoxigenic Bacteroides fragilis, and colibactin-producing Escherichia coli. Further, we consider how microbes can be used in diagnosing colorectal cancer and manipulating the tumor environment to encourage better patient outcomes in response to immunotherapy treatments.

Characterization of Human Cytomegalovirus Genome Diversity in Immunocompromised Hosts by Whole-Genome Sequencing Directly From Clinical Specimens

Background

Advances in next-generation sequencing (NGS) technologies allow comprehensive studies of genetic diversity over the entire genome of human cytomegalovirus (HCMV), a significant pathogen for immunocompromised individuals.

Methods

Next-generation sequencing was performed on target enriched sequence libraries prepared directly from a variety of clinical specimens (blood, urine, breast milk, respiratory samples, biopsies, and vitreous humor) obtained longitudinally or from different anatomical compartments from 20 HCMV-infected patients (renal transplant recipients, stem cell transplant recipients, and congenitally infected children).

Results

De novo-assembled HCMV genome sequences were obtained for 57 of 68 sequenced samples. Analysis of longitudinal or compartmental HCMV diversity revealed various patterns: no major differences were detected among longitudinal, intraindividual blood samples from 9 of 15 patients and in most of the patients with compartmental samples, whereas a switch of the major HCMV population was observed in 6 individuals with sequential blood samples and upon compartmental analysis of 1 patient with HCMV retinitis. Variant analysis revealed additional aspects of minor virus population dynamics and antiviral-resistance mutations.

Conclusions

In immunosuppressed patients, HCMV can remain relatively stable or undergo drastic genomic changes that are suggestive of the emergence of minor resident strains or de novo infection.

Culture-Independent Diagnostics for Health Security

The past decade has seen considerable development in the diagnostic application of nonculture methods, including nucleic acid amplification-based methods and mass spectrometry, for the diagnosis of infectious diseases. The implications of these new culture-independent diagnostic tests (CIDTs) include bypassing the need to culture organisms, thus potentially affecting public health surveillance systems, which continue to use isolates as the basis of their surveillance programs and to assess phenotypic resistance to antimicrobial agents. CIDTs may also affect the way public health practitioners detect and respond to a bioterrorism event. In response to a request from the Department of Homeland Security, Los Alamos National Laboratory and the Centers for Disease Control and Prevention cosponsored a workshop to review the impact of CIDTs on the rapid detection and identification of biothreat agents. Four panel discussions were held that covered nucleic acid amplification-based diagnostics, mass spectrometry, antibody-based diagnostics, and next-generation sequencing. Exploiting the extensive expertise available at this workshop, we identified the key features, benefits, and limitations of the various CIDT methods for providing rapid pathogen identification that are critical to the response and mitigation of a bioterrorism event. After the workshop we conducted a thorough review of the literature, investigating the current state of these 4 culture-independent diagnostic methods. This article combines information from the literature review and the insights obtained at the workshop.

Longitudinal genomic surveillance of Plasmodium falciparum malaria parasites reveals complex genomic architecture of emerging artemisinin resistance

Background

Artemisinin-based combination therapies are the first line of treatment for Plasmodium falciparum infections worldwide, but artemisinin resistance has risen rapidly in Southeast Asia over the past decade. Mutations in the kelch13 gene have been implicated in this resistance. We used longitudinal genomic surveillance to detect signals in kelch13 and other loci that contribute to artemisinin or partner drug resistance. We retrospectively sequenced the genomes of 194 P. falciparum isolates from five sites in Northwest Thailand, over the period of a rapid increase in the emergence of artemisinin resistance (2001–2014).

Results

We evaluate statistical metrics for temporal change in the frequency of individual SNPs, assuming that SNPs associated with resistance increase in frequency over this period. After Kelch13-C580Y, the strongest temporal change is seen at a SNP in phosphatidylinositol 4-kinase, which is involved in a pathway recently implicated in artemisinin resistance. Furthermore, other loci exhibit strong temporal signatures which warrant further investigation for involvement in artemisinin resistance evolution. Through genome-wide association analysis we identify a variant in a kelch domain-containing gene on chromosome 10 that may epistatically modulate artemisinin resistance.

Conclusions

This analysis demonstrates the potential of a longitudinal genomic surveillance approach to detect resistance-associated gene loci to improve our mechanistic understanding of how resistance develops. Evidence for additional genomic regions outside of the kelch13 locus associated with artemisinin-resistant parasites may yield new molecular markers for resistance surveillance, which may be useful in efforts to reduce the emergence or spread of artemisinin resistance in African parasite populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1204-4) contains supplementary material, which is available to authorized users.

Whole genome sequencing of Plasmodium falciparum from dried blood spots using selective whole genome amplification

Background

Translating genomic technologies into healthcare applications for the malaria parasite Plasmodium falciparum has been limited by the technical and logistical difficulties of obtaining high quality clinical samples from the field. Sampling by dried blood spot (DBS) finger-pricks can be performed safely and efficiently with minimal resource and storage requirements compared with venous blood (VB). Here, the use of selective whole genome amplification (sWGA) to sequence the P. falciparum genome from clinical DBS samples was evaluated, and the results compared with current methods that use leucodepleted VB.

Methods

Parasite DNA with high (>95%) human DNA contamination was selectively amplified by Phi29 polymerase using short oligonucleotide probes of 8–12 mers as primers. These primers were selected on the basis of their differential frequency of binding the desired (P. falciparum DNA) and contaminating (human) genomes.

Results

Using sWGA method, clinical samples from 156 malaria patients, including 120 paired samples for head-to-head comparison of DBS and leucodepleted VB were sequenced. Greater than 18-fold enrichment of P. falciparum DNA was achieved from DBS extracts. The parasitaemia threshold to achieve >5× coverage for 50% of the genome was 0.03% (40 parasites per 200 white blood cells). Over 99% SNP concordance between VB and DBS samples was achieved after excluding missing calls.

Conclusion

The sWGA methods described here provide a reliable and scalable way of generating P. falciparum genome sequence data from DBS samples. The current data indicate that it will be possible to get good quality sequence on most if not all drug resistance loci from the majority of symptomatic malaria patients. This technique overcomes a major limiting factor in P. falciparum genome sequencing from field samples, and paves the way for large-scale epidemiological applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1641-7) contains supplementary material, which is available to authorized users.

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.