Development of a capture sequencing assay for enhanced detection and genotyping of tick-borne pathogens

Empirical assessment of the enrichment-based metagenomic methods in identifying diverse respiratory pathogens

Probe-based nucleic acid enrichment represents an effective route to enhance the detection capacity of next-generation sequencing (NGS) in a set of clinically diverse and relevant microbial species. In this study, we assessed the effect of the enrichment-based sequencing on identifying respiratory infections using tiling RNA probes targeting 76 respiratory pathogens and sequenced using both Illumina and Oxford Nanopore platforms. Forty respiratory swab samples pre-tested for a panel of respiratory pathogens by qPCR were used to benchmark the sequencing data. We observed a general improvement in sensitivity after enrichment. The overall detection rate increased from 73 to 85% after probe capture detected by Illumina. Moreover, enrichment with probe sets boosted the frequency of unique pathogen reads by 34.6 and 37.8-fold for Illumina DNA and cDNA sequencing, respectively. This also resulted in significant improvements on genome coverage especially in viruses. Despite these advantages, we found that library pooling may cause reads mis-assignment, probably due to crosstalk issues arise from post-capture PCR and from pooled sequencing, thus increasing the risk of bleed-through signal. Taken together, an overall improvement in the breadth and depth of pathogen coverage is achieved using enrichment-based sequencing method. For future applications, automated library processing and pooling-free sequencing could enhance the precision and timeliness of probe enrichment-based clinical metagenomics.

Turning the needle into the haystack: Culture-independent amplification of complex microbial genomes directly from their native environment

High-throughput sequencing (HTS) has revolutionized microbiology, but many microbes exist at low abundance in their natural environment and/or are difficult, if not impossible, to culture in the laboratory. This makes it challenging to use HTS to study the genomes of many important microbes and pathogens. In this review, we discuss the development and application of selective whole genome amplification (SWGA) to allow whole or partial genomes to be sequenced for low abundance microbes directly from complex biological samples. We highlight ways in which genomic data generated by SWGA have been used to elucidate the population dynamics of important human pathogens and monitor development of antimicrobial resistance and the emergence of potential outbreaks. We also describe the limitations of this method and propose some potential innovations that could be used to improve the quality of SWGA and lower the barriers to using this method across a wider range of infectious pathogens.

Human Babesia odocoilei and Bartonella spp. co-infections in the Americas

BACKGROUND

In recent years, Babesia and Bartonella species co-infections in patients with chronic, nonspecific illnesses have continued to challenge and change the collective medical understanding of "individual pathogen" vector-borne infectious disease dynamics, pathogenesis and epidemiology. The objective of this case series is to provide additional molecular documentation of Babesia odocoilei infection in humans in the Americas and to emphasize the potential for co-infection with a Bartonella species.

METHODS

The development of improved and more sensitive molecular diagnostic techniques, as confirmatory methods to assess active infection, has provided increasing clarity to the healthcare community.

RESULTS

Using a combination of different molecular diagnostic approaches, infection with Babesia odocoilei was confirmed in seven people suffering chronic non-specific symptoms, of whom six were co-infected with one or more Bartonella species.

CONCLUSIONS

We conclude that infection with Babesia odocoilei is more frequent than previously documented and can occur in association with co-infection with Bartonella spp.

Retrospective Analysis of Potential Lyme Disease Clinical Cases in Argentina

Lyme disease, a multisystemic infectious disorder caused by pathogenic spirochetes of the genus Borrelia transmitted by the bite of ticks, typically from the family Ixodidae, pose a significant public health issue worldwide. The Borrelia burgdorferi sensu lato (s.l.) group encompasses the Borrelia Lyme Group (LG), Borrelia Echidna-Reptile Group (REPG), and Borrelia Relapsing Fever Group (RFG), with some species remaining unclassified due to culturing challenges. Research into B. burgdorferi s.l. infection (Lyme Group) has intensified, focusing on its epidemiology, diagnosis, and treatment. Originally identified in North America and Europe, Lyme disease has now become a global concern, with Latin American countries reporting the microorganism, the disease, and/or its vectors. In Argentina, the presence of B. burgdorferi and Lyme disease has sparked significant scientific and medical debate. Ecological changes due to climate and habitat shifts have expanded the geographical distribution of these ticks. Argentina, with its diverse geography and climate, hosts various tick species that could potentially act as Lyme disease vectors, raising important public health questions. The confirmed presence of B. burgdorferi s.l. and Lyme disease in Argentina remains contentious but relevant, necessitating thorough scientific and medical examination. This work aims to enhance understanding and discussion of Lyme disease in Argentina by presenting clinical cases and their laboratory analyses, highlighting the disease's presence and implications in the country. Through documenting suspected clinical cases and analyzing available data on B. burgdorferi and Lyme disease in Argentina, this study seeks to contribute to the understanding of the disease's current status and inform future research, prevention, and control strategies in the region. The goal is to provide a basis for addressing Lyme disease's public health impact in Argentina and promote further investigation into this evolving issue.

Update on tick-borne pathogens detection methods within ticks

Ticks are known vectors of various pathogens, including bacteria, parasites and viruses, that impact both animal and human health. Improving knowledge of the distribution of tick-borne pathogens, combined with their early detection in ticks, are essential steps to fight against tick-borne diseases and mitigate their impacts. Here we give an overview of what are the common methods of pathogen detection in tick samples, including recent developments concerning how to handle tick samples, get access to tick-borne pathogens by chemical or physical disruption of the ticks, and methods used for the RNA/DNA extraction steps. Furthermore, we discuss promising tools that are developed for other sample types such as serum or blood to detect tick-borne pathogens, and those that could be used in the future for tick samples.

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.

RT-PCR genotyping assays to identify SARS-CoV-2 variants in England in 2021: a design and retrospective evaluation study

BACKGROUND

Whole-genome sequencing (WGS) is the gold standard diagnostic tool to identify and genetically characterise emerging pathogen mutations (variants), but cost, capacity, and timeliness limit its use when large populations need rapidly assessing. We assessed the potential of genotyping assays to provide accurate and timely variant information at scale by retrospectively examining surveillance for SARS-CoV-2 variants in England between March and September, 2021, when genotyping assays were used widely for variant detection.

METHODS

We chose a panel of four RT-PCR genotyping assays to detect circulating variants of SARS-COV-2 in England and developed a decision algorithm to assign a probable SARS-CoV-2 variant to samples using the assay results. We extracted surveillance data from the UK Health Security Agency databases for 115 934 SARS-CoV-2-positive samples (March 1-Sept 6, 2021) when variant information was available from both genotyping and WGS. By comparing the genotyping and WGS variant result, we calculated accuracy metrics (ie, sensitivity, specificity, and positive predictive value [PPV]) and the time difference between the sample collection date and the availability of variant information. We assessed the number of samples with a variant assigned from genotyping or WGS, or both, over time.

FINDINGS

Genotyping and an initial decision algorithm (April 10-May 11, 2021 data) were accurate for key variant assignment: sensitivities and PPVs were 0·99 (95% CI 0·99-0·99) for the alpha, 1·00 (1·00-1·00) for the beta, and 0·91 (0·80-1·00) for the gamma variants; specificities were 0·97 (0·96-0·98), 1·00 (1·00-1·00), and 1·00 (1·00-1·00), respectively. A subsequent decision algorithm over a longer time period (May 27-Sept 6, 2021 data) remained accurate for key variant assignment: sensitivities were 0·91 (95% CI 0·74-1·00) for the beta, 0·98 (0·98-0·99) for the delta, and 0·93 (0·81-1·00) for the gamma variants; specificities were 1·00 (1·00-1·00), 0·96 (0·96-0·97), and 1·00 (1·00-1·00), respectively; and PPVs were 0·83 (0·62-1·00), 1·00 (1·00-1·00), and 0·78 (0·59-0·97), respectively. Genotyping produced variant information a median of 3 days (IQR 2-4) after the sample collection date, which was faster than with WGS (9 days [8-11]). The flexibility of genotyping enabled a nine-times increase in the quantity of samples tested for variants by this method (from 5000 to 45 000).

INTERPRETATION

RT-PCR genotyping assays are suitable for high-throughput variant surveillance and could complement WGS, enabling larger scale testing for known variants and timelier results, with important implications for effective public health responses and disease control globally, especially in settings with low WGS capacity. However, the choice of panels of RT-PCR assays is highly dependent on database information on circulating variants generated by WGS, which could limit the use of genotyping assays when new variants are emerging and spreading rapidly.

FUNDING

UK Health Security Agency and National Institute for Health Research Health Protection Research Unit in Emergency Preparedness and Response.

Tick-Borne Co-Infections: Challenges in Molecular and Serologic Diagnoses

Co-infections are a poorly understood aspect of tick-borne diseases. In the United States alone, nineteen different tick-borne pathogens have been identified. The majority of these agents are transmitted by only two tick species, Ixodes scapularis and Amblyomma americanum. Surveillance studies have demonstrated the presence of multiple pathogens in individual ticks suggesting a risk of polymicrobial transmission to humans. However, relatively few studies have explored this relationship and its impact on human disease. One of the key factors for this deficiency are the intrinsic limitations associated with molecular and serologic assays employed for the diagnosis of tick-borne diseases. Limitations in the sensitivity, specificity and most importantly, the capacity for inclusion of multiple agents within a single assay represent the primary challenges for the accurate detection of polymicrobial tick-borne infections. This review will focus on outlining these limitations and discuss potential solutions for the enhanced diagnosis of tick-borne co-infections.

Multiple piroplasm parasites (Apicomplexa: Piroplasmida) in northeastern populations of the invasive Asian longhorned tick, Haemaphysalis longicornis Neumann (Ixodida: Ixodidae), in the United States

Piroplasms, which include the agents of cattle fever and human and dog babesiosis, are a diverse group of blood parasites of significant veterinary and medical importance. The invasive Asian longhorned tick, Haemaphysalis longicornis, is a known vector of piroplasms in its native range in East Asia and invasive range in Australasia. In the USA, H. longicornis has been associated with Theileria orientalis Ikeda outbreaks that caused cattle mortality. To survey invasive populations of H. longicornis for a broad range of piroplasms, 667 questing H. longicornis collected in 2021 from 3 sites in New Jersey, USA, were tested with generalist piroplasm primers targeting the 18S small subunit rRNA (395–515 bp, depending on species) and the cytochrome b oxidase loci (1009 bp). Sequences matching Theileria cervi type F (1 adult, 5 nymphs), an unidentified Theileria species (in 1 nymph), an undescribed Babesia sensu stricto (‘true’ Babesia, 2 adults, 2 nymphs), a Babesia sp. Coco (also a ‘true Babesia’, 1 adult, 1 nymph), as well as Babesia microti S837 (1 adult, 4 nymphs) were recovered. Babesia microti S837 is closely related to the human pathogen B. microti US-type. Additionally, a 132 bp sequence matching the cytochrome b locus of deer, Odocoileus virginanus, was obtained from 2 partially engorged H. longicornis. The diverse assemblage of piroplasms now associated with H. longicornis in the USA spans 3 clades in the piroplasm phylogeny and raises concerns of transmission amplification of veterinary pathogens as well as spillover of pathogens from wildlife to humans.

BosR and PlzA reciprocally regulate RpoS function to sustain Borrelia burgdorferi in ticks and mammals

The RNA polymerase alternative σ factor RpoS in Borrelia burgdorferi (Bb), the Lyme disease pathogen, is responsible for programmatic-positive and -negative gene regulation essential for the spirochete's dual-host enzootic cycle. RpoS is expressed during tick-to-mammal transmission and throughout mammalian infection. Although the mammalian-phase RpoS regulon is well described, its counterpart during the transmission blood meal is unknown. Here, we used Bb-specific transcript enrichment by tick-borne disease capture sequencing (TBDCapSeq) to compare the transcriptomes of WT and ΔrpoS Bb in engorged nymphs and following mammalian host-adaptation within dialysis membrane chambers. TBDCapSeq revealed dramatic changes in the contours of the RpoS regulon within ticks and mammals and further confirmed that RpoS-mediated repression is specific to the mammalian-phase of Bb's enzootic cycle. We also provide evidence that RpoS-dependent gene regulation, including repression of tick-phase genes, is required for persistence in mice. Comparative transcriptomics of engineered Bb strains revealed that the Borrelia oxidative stress response regulator (BosR), a noncanonical Fur family member, and the cyclic diguanosine monophosphate (c-di-GMP) effector PlzA reciprocally regulate the function of RNA polymerase complexed with RpoS. BosR is required for RpoS-mediated transcription activation and repression in addition to its well-defined role promoting transcription of rpoS by the RNA polymerase alternative σ factor RpoN. During transmission, ligand-bound PlzA antagonizes RpoS-mediated repression, presumably acting through BosR.

Genomic characterization of Francisella tularensis and other diverse Francisella species from complex samples

Francisella tularensis, the bacterium that causes the zoonosis tularemia, and its genetic near neighbor species, can be difficult or impossible to cultivate from complex samples. Thus, there is a lack of genomic information for these species that has, among other things, limited the development of robust detection assays for F. tularensis that are both specific and sensitive. The objective of this study was to develop and validate approaches to capture, enrich, sequence, and analyze Francisella DNA present in DNA extracts generated from complex samples. RNA capture probes were designed based upon the known pan genome of F. tularensis and other diverse species in the family Francisellaceae. Probes that targeted genomic regions also present in non-Francisellaceae species were excluded, and probes specific to particular Francisella species or phylogenetic clades were identified. The capture-enrichment system was then applied to diverse, complex DNA extracts containing low-level Francisella DNA, including human clinical tularemia samples, environmental samples (i.e., animal tissue and air filters), and whole ticks/tick cell lines, which was followed by sequencing of the enriched samples. Analysis of the resulting data facilitated rigorous and unambiguous confirmation of the detection of F. tularensis or other Francisella species in complex samples, identification of mixtures of different Francisella species in the same sample, analysis of gene content (e.g., known virulence and antimicrobial resistance loci), and high-resolution whole genome-based genotyping. The benefits of this capture-enrichment system include: even very low target DNA can be amplified; it is culture-independent, reducing exposure for research and/or clinical personnel and allowing genomic information to be obtained from samples that do not yield isolates; and the resulting comprehensive data not only provide robust means to confirm the presence of a target species in a sample, but also can provide data useful for source attribution, which is important from a genomic epidemiology perspective.

Capture Sequencing Enables Sensitive Detection of Tick-Borne Agents in Human Blood

Assay sensitivity can be a limiting factor in the use of PCR as a tool for the detection of tick-borne pathogens in blood. We evaluated the performance of Tick-borne disease Capture Sequencing Assay (TBDCapSeq), a capture sequencing assay targeting tick-borne agents, to test 158 whole blood specimens obtained from the Lyme Disease Biobank. These included samples from 98 individuals with signs and symptoms of acute Lyme disease, 25 healthy individuals residing in Lyme disease endemic areas, and 35 samples collected from patients admitted to the Massachusetts General Hospital or referred to the infectious disease clinic. Compared to PCR, TBDCapSeq had better sensitivity and could identify infections with a wider range of tick-borne agents. TBDCapSeq identified a higher rate of samples positive for Borrelia burgdorferi (8 vs. 1 by PCR) and Babesia microti (26 vs. 15 by PCR). TBDCapSeq also identified previously unknown infections with Borrelia miyamotoi, Ehrlichia, and Rickettsia species. Overall, TBDCapSeq identified a pathogen in 43 samples vs. 23 using PCR, with four co-infections detected versus zero by PCR. We conclude that capture sequencing enables superior detection of tick-borne agents relative to PCR.