Applying the pathogen-targeted next-generation sequencing method to pathogen identification in cerebrospinal fluid.

BACKGROUND: The cerebrospinal fluid (CSF) culture is a widely used method for the diagnosis of meningitis, but its detection sensitivity is low. Several new methods have been developed for pathogen detection, including metagenomic next-generation sequencing (mNGS) and pathogen-targeted NGS (ptNGS). In this study, we aimed to evaluate the performance of ptNGS in pathogen detection in CSF.
METHODS: CSF specimens were acquired from 38 patients with meningitis who were diagnosed at Xuanwu Hospital, Capital Medical University between October 2020 and February 2021. DNA was extracted from the CSF samples, and pathogens were identified using both ptNGS and mNGS. SPSS 22.0 software was used to compare the pathogen detection performance of ptNGS and mNGS in CSF.
RESULTS: Among the 38 patients with meningitis, 14 had a non-infectious disease (NID) and 24 had an infectious disease (ID). Of the 38 samples, both ptNGS and mNGS detected 9 (23.7%) positive samples, and 12 (31.6%) negative samples. Thirteen (34.2%) samples were detected to be positive by ptNGS only, and 4 (10.5%) were detected to be positive by mNGS only. The positivity rate detected by ptNGS for the ID group was higher than that detected by mNGS (P=0.080), and the positivity rates detected by ptNGS and mNGS for the NID group were comparable. The positive predictive value (PPV) and negative predictive value (NPV) of diagnosing an ID by ptNGS were 77.3% and 56.3%, respectively. While, the PPV and NPV of diagnosing an ID by mNGS were 76.9% and 44.0%, respectively. ptNGS increased the sensitivity rate by approximately 70%. The sensitivity rate of ptNGS was higher than that of mNGS (70.8% vs. 41.7%), while the specificity rate of mNGS was higher than that of ptNGS (78.6% vs. 64.3%). Additionally, ptNGS required a shorter time for pathogen diagnosis (15 vs. 24 hrs) and had lower costs than mNGS.
CONCLUSIONS: ptNGS has a number of advantages over mNGS, including its sensitivity, timeliness, and economy, all factors that are important considerations in clinical use. 2021 Annals of Translational Medicine. All rights reserved.

Introduction

Meningitis is an inflammation of the fluid and meninges surrounding the brain and spinal cord. Meningitis, especially acute meningitis, represents a threat to people’s health. Clinically, the cause of acute meningitis cases is not identified in approximately 50% of patients (1-3). A failure to make a timely diagnosis in patients with meningitis contributes to poor patient outcomes. Accurate information about important etiological agents is necessary to ensure appropriate management. Viral infections, followed by bacterial infections, are the most common causes of meningitis. Conversely, fungal and parasitic infections rarely cause meningitis. Non-infectious causes, such as chemical reactions, drug allergies, cancer, and inflammatory diseases (e.g., sarcoidosis) can also cause meningitis.

The cerebrospinal fluid (CSF) culture has been widely used to diagnose meningitis (4). However, for cases with negative CSF cultures, the diagnosis of bacterial meningitis may be established by positive blood cultures (5). Due to the low sensitivity of CSF cultures, several alternative methods have been introduced for pathogen detection, including polymerase chain reaction (PCR), film array, immunoassay, microarray, and next-generation sequencing (NGS) (6-9). Among these methods, immunoassay shows best timeliness; PCR and film array show best accuracy. But the throughput of all these methods are low. NGS base methods show much higher throughput. It can also detects not only pathogens but also their drug-resistant genes in on run. Metagenomic NGS (mNGS) is an increasingly rapid and comparatively low-cost means of screening CSF in an unbiased manner for a broad range of human pathogens (10,11). mNGS is a single diagnostic test with promising potential for clinical application (10,11). Recently, the new method of pathogen-targeted NGS (ptNGS) was developed and has been used to identify pathogens in respiratory tract infection or mycobacterium infection cases (12,13). At present, it seems that ptNGS has the advantages of detection sensitivity not affected by human genome and background bacteria, detection sensitivity not affected by pathogen genome size, lower detection cost, reduced sample transportation requirements, and quantitative detection of pathogens. However, its pathogen detection performance in CSF remains unclear. In this study, ptNGS and mNGS were compared with conventional test results respectively, and then ptNGS and mNGS technologies based on the NGS platform were compared horizontally. This project is the first to apply ptNGS to the pathogenic diagnosis of CSF in patients with central nervous system infection.

In the present study, we aimed to evaluate the performance of ptNGS in pathogen detection in CSF. CSF samples from 38 patients with meningitis, who had been diagnosed using the traditional diagnostic method of CSF, were examined using the mNGS and ptNGS methods. The performance of the ptNGS and mNGS was then compared. We present the following article in accordance with the MDAR reporting checklist (available at https://dx.doi.org/10.21037/atm-21-5488).

Methods

Patients and samples

CSF specimens were acquired from 38 patients who were diagnosed with meningitis using the traditional diagnostic method of CSF at Xuanwu Hospital, Capital Medical University between October 2020 and February 2021. Based on their final diagnoses, the patients were categorized into an infectious disease (ID) group and a non-infectious disease (NID) group. All procedures performed in this study involving human participants were in accordance with the Declaration of Helsinki (as revised in 2013). This study was approved by the Ethics Committee of Xuanwu Hospital, Capital Medical University [No. (2020)103]. All participants provided written informed consent. The CSF specimens were stored at −20 to −80 °C awaiting use.

DNA extraction

The CSF samples were centrifuged at 13,000 rpm for 10 min and were ground on a Grinding Mill (Tiss-24, Jingxin, Shanghai, China) at 60 HZ for 10 min. The ground samples were then used for DNA/RNA extraction and purification (ZymoBIOMICS DNA Miniprep Kit, R2002) in accordance with the manufacturer’s instructions. The extracted DNA was used to construct a library. Reduce host DNA contamination from clinical CSF samples: (I) before nucleic acid extraction: osmotic dissolution lyse human cells to retain the intact pathogen; the released human DNA is treated and degraded with DNase or chemical reagents (such as monoaziridin), and then the microbial nucleic acid is extracted (QIAamp DNA Microbiome Kit); (II) after nucleic acid extraction: after the nucleic acid of the specimen is extracted, the anti-methylated DNA specific binding protein antibody is used to selectively bind and remove the methylated DNA of CpG island (NEBNext Microbiome DNA Enrichment Kit).

ptNGS

A total of 544 meningitis associated pathogens were included in the identification model for meningitis diagnosis. First, the DNA sequences used in the identification model were selected as the targeted fragments. The primers were designed and synthesized at Pathogeno Biotech and Sangon Biotech (Shanghai, China), respectively. The targeted gene sequences were amplified and enriched by multiplex PCR (as described previously) (11). A sequencing library was constructed by adding the sequencing connectors to the purified PCR products and using DNA purification magnetic beads for purification. Targeted gene sequencing (300 cycles) was performed on a MiSeq system (Illumina, Inc., San Diego, CA, USA) using the MiSeq reagent kit v2. FastQ files were generated by the MiSeq Reporter software. Offline data generated by the MiSeq system were identified and counted by an adapter, and reads with a double-end length >60 bp were retained. Among the high-quality data, reads with <60 bp at either end, single-end primer recognition, or non-specific primer binding, were re-examined and deleted. Clean read pairs were obtained for identification and sequence alignment. Before aligning to pathogen reference sequences, clean read pairs were first aligned to human reference genome. Read pairs aligned to human genome were dropped to reduce host DNA contamination. Remained read pairs were aligned to pathogen reference sequences. Read paired counts for each pathogen were generated for further analysis.

ptNGS quality control: (I) nucleic acid purity control: The purity of the extracted nucleic acid is measured by Nanodrop. Generally, A260/A280 is greater than 1.8, and A260/A230 greater than 2.0 is considered a pure sample; (II) error-proof label control: every A unique error-proof label is added to each sample at the beginning of nucleic acid extraction. After the sequencing is completed, the correct and stable error-proof label signal in each sample is analyzed to monitor sample confusion or extraction/amplification failure; (II) library quality control: the normal library fragment size is about 350 bp, there is no obvious dimer or non-specific band, and the library concentration is greater than 1 ng/μL; (IV) sequencing quality control: after filtering and comparing and analyzing the total data volume of a single library, the effective data volume should not be less than 0.08 M reads.

mNGS

The mNGS method has been described previously (8). Briefly, DNA libraries were constructed using an end-repair method; adapters were added overnight, and amplified by PCR before the analysis using an Ion Torrent Proton Sequencer (Life Technologies, Carlsbad, CA, USA). Before sequencing, the quality of the DNA libraries was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) in conjunction with a quantitative PCR. The qualified DNA libraries were prepared by emulsion PCR in the OneTouch system and sequenced on the Ion Torrent Proton (Life Technologies, South San Francisco, CA, USA) sequencing platform. High-quality sequencing data were obtained by removing low-quality and short reads (length <35 bp), and Burrows-Wheeler alignment was then used to subtract human host sequences. After removing low-complexity reads, the remaining data were then classified by aligning them to 4 microbial genome databases (including viruses, bacteria, fungi and parasites) simultaneously. RefSeq contains 4,189 whole-genome sequences of viral taxa, 2,328 bacterial genomes or scaffolds, 199 fungi associated with human infection, and 135 parasites related to human diseases.

Statistical analysis

SPSS 22.0 software was used for the data analysis. A Pearson chi-squared (χ2) test or a Fisher’s exact test was used to compare the frequency of the categorical data. A t-test (normal distribution) or a non-parametric test (Mann-Whitney test, non-normal distribution) was used to compare the quantitative data between the two groups. For all presented data, a P value <0.05 was considered statistically significant.

Results

Patient characteristics

The clinicopathologic information of 38 patients with meningitis is listed in Table S1. Notably, 14 patients had a NID, and 24 had an ID. Specifically, there was 18 cases of viral meningoencephalitis (VM), 3 cases of purulent meningitis (PM), 2 cases of tuberculous meningitis (TBM), and 1 case of fungal meningitis (FM) patients. Six (42.9%) and 12 (50.0%) male patients had a NID and an ID, respectively. The median age of patients with a NID and an ID was 39 (range, 11–77) and 46 (range, 15–68) years old, respectively.

Consistency between ptNGS and mNGS detection

Both ptNGS and mNGS were used to detect the pathogens in each patient. The pathogen results obtained by ptNGS and mNGS are presented in Table S1. As shows, of the 38 cases, both ptNGS and mNGS identified 9 (2 NID and 7 ID; 23.7%) positive samples and 12 (8 NID and 4 ID; 31.6%) negative samples. ptNGS identified 13 (34.2%) positive samples only (3 NID and 10 ID), and mNGS identified 4 (10.5%) positive samples only (1 NID and 3 ID). The consistency between ptNGS and mNGS among all the enrolled patients was estimated, and the Kappa value was 0.148 (P=0.490; see ). In the ID group, ptNGS and mNGS both identified positive samples in 7 patients (2 results were completely matched and 5 were partially matched) (see Table S2 and ).

Figure 1

Overview of sample classification and comparison. (A) Based on diagnoses obtained using the traditional diagnostic method of CSF, the samples were categorized into ID or NID groups. All samples were examined by ptNGS and mNGS for the concordance analysis, and ID and NID patients were used to evaluate their diagnostic performance. (B) For the double-positive subset in the ID group (n=7), 2 patients had completely matched results and 5 had partially matched results. CSF, cerebrospinal fluid; ID, infectious disease; NID, non-infectious disease; ptNGS, pathogen-targeted next-generation sequencing; mNGS, metagenomic next-generation sequencing.

Table 1

The consistency between ptNGS and mNGS

	ptNGS (+)	ptNGS (−)	Total	Kappa	P value
mNGS (+)	9	4	13	0.148	0.490
mNGS (−)	13	12	25
Total	22	16	38

The symbol “+” stands for “positive”, and the symbol “−” stands for “negative”. ptNGS, pathogen-targeted next-generation sequencing; mNGS, metagenomic next-generation sequencing.

Comparison of the diagnostic performance of ptNGS and mNGS

In the ID group (n=24), 17 (70.8%) patients were diagnosed with confirmed pathogens by ptNGS, while only 10 (41.7%) were diagnosed by mNGS (see ).

Figure 2

Positivity rate comparison between ptNGS and mNGS. (A) The differences between ptNGS and mNGS in infectious disease (n=24). (B) The differences between ptNGS and mNGS in non-infectious disease (n=14). ptNGS, pathogen-targeted next-generation sequencing; mNGS, metagenomic next-generation sequencing; ID, infectious disease; NID, non-infectious disease.

In the ID group, the positivity rate detected by ptNGS for was higher than that detected by mNGS (P=0.080). In the NID group (n=14), 5 (35.7%) patients were diagnosed with confirmed pathogens by ptNGS and 3 (21.4%) by mNGS. The positivity rates for the NID group between ptNGS and mNGS were comparable (see ).

The performance of ptNGS and mNGS in pathogen detection was compared (see ). The positive predictive value (PPV) and the negative predictive value (NPV) of diagnosing ID by ptNGS were 77.3% and 56.3%, respectively. Conversely, the PPV and the NPV of diagnosing ID by mNGS were 76.9% and 44.0%, respectively. ptNGS increased the sensitivity rate by approximately 70%. The sensitivity rate of ptNGS was higher than that of mNGS (70.8% vs. 41.7%), while the specificity rate of mNGS was higher than that of ptNGS (78.6% vs. 64.3%).

Table 2

The pathogen detection performance of ptNGS and mNGS

Method	Sensitivity (%)	Specificity (%)	PPV (%)	NPV (%)
ptNGS	70.8	64.3	77.3	56.3
mNGS	41.7	78.6	76.9	44.0

ptNGS, pathogen-targeted next-generation sequencing; mNGS, metagenomic next-generation sequencing; PPV, positive predictive value; NPV, negative predictive value.

It took approximately 15 and 24 hrs for pathogen diagnosis by ptNGS and mNGS, respectively. mNGS has the ability to detect unknown pathogens; however, ptNGS only detects major pathogens with a clinical incidence of more than 95%. ptNGS can quantify pathogens at the copy number level, while mNGS can only perform relative quantifications. ptNGS can detect 50 copies of pathogens per milliliter of CSF sample, while the detection limit of mNGS depends on the context of the samples. provides a summary of the comparison of the technical advantages of ptNGS and mNGS.

Table 3

Comparison of the technical advantages of ptNGS and mNGS

Items	ptNGS	mNGS
Detection target	Targeted pathogens	Metagenome
Detection period	15 hrs	24 hrs
Detection of unknown pathogens	No	Yes
Limit of detection	50 copies per mL	Context dependent
Quantitative detection	Copy number level quantification	Relative quantification
Economic cost	Low	High

ptNGS, pathogen-targeted next-generation sequencing; mNGS, metagenomic next-generation sequencing.

Discussion

Meningitis may be infectious or non-infectious. The main types of infectious meningitis include VM, PM, TBM, and FM. Prognosis and treatment vary for different types of meningitis. Bacterial meningitis is a serious ID that can be fatal to both children and adults. The incidence and mortality rates of bacterial meningitis vary depending on the type of pathogen (14). To effectively treat bacterial meningitis, the microorganisms and their antibiotic susceptibility patterns need to be rapidly identified.

Currently, the CSF culture is the gold standard for diagnosing bacterial meningitis. However, the low bacterial growth rates cause high false negative results (15-17). Thus, new test methods need to be developed. mNGS, which is a rapid and high-throughput method for pathogen detection, has been applied to CSF samples in several studies (6-8). Miao et al. systematically compared mNGS and the CSF culture, and found that mNGS was advantageous in several aspects (10). mNGS has a higher sensitivity for pathogen identification and is less affected by previous exposure to antibiotics; thus, it has become a promising technology for the detection of IDs. However, it is difficult for mNGS to define specific microbial profiles that are diagnostic or predictive of disease development. Additionally, mNGS typically requires at least 20 million reads for each sample library and has high costs. Therefore, a new method is needed to solve these limitations. ptNGS requires only 0.08 million reads for a single sample library, which greatly improves the detection throughput and reduces the sequencing costs. Based on the superior performance of mNGS in pathogen identification, we evaluated the performance of ptNGS by comparing it to that of the mNGS.

Li et al. reported that compared with conventional tests, ptNGS (targeted amplicon sequencing assay) has better sensitivity and specificity for the diagnosis of alveolar lavage fluid bacterial and viral pathogens in severe community-acquired pneumonia patients (17). Chao et al. revealed that using pathogen target sequencing method to detect pathogenic bacteria in the sputum of patients with Acute lower respiratory infections is more accurate and sensitive than traditional Sputum culture detection (14). In this study, we use CSF samples to compare the differences between ptNGS and mNGS in the detection of meningitis patients. The analysis showed that ptNGS had a higher sensitivity rate than that of the mNGS (70.8% vs. 41.7%). However, the specificity of ptNGS was low (56.3%). This may be due to the false negative results of traditional diagnostic methods. Surprisingly, the positivity rates of ptNGS and mNGS were not higher than traditional diagnostic methods (68.4%/55.3% vs. 63.2%). This may be due to the relatively small size of the study cohort. Additionally, compared to mNGS, ptNGS took a shorter time to identify pathogens (15 vs. 24 hrs) and had lower economic costs. Although this is a small sample size study and more studies need to evaluate the performance of ptNGS in clinical use, it is proven that ptNGS has a promising potential method for guiding clinical detection.

Conclusions

ptNGS has a number of advantages over mNGS, including its sensitivity, timeliness, and economy, all of which are usually considered in clinical use.

The article’s supplementary files as