Next Article in Journal
Neonatal SARS-CoV-2 Infection: Practical Tips
Next Article in Special Issue
Molecular Characterization of 4/91 Infectious Bronchitis Virus Leading to Studies of Pathogenesis and Host Responses in Laying Hens
Previous Article in Journal
Existence of Similar Leptospira Serovars among Dog Keepers and Their Respective Dogs in Mwanza, Tanzania, the Need for a One Health Approach to Control Measures
Previous Article in Special Issue
Molecular Epidemiological Survey of Canine Parvovirus Circulating in China from 2014 to 2019
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Brief Report

Direct Diagnosis of Echovirus 12 Meningitis Using Metagenomic Next Generation Sequencing

1
IHU Méditerranée Infection, 13005 Marseille, France
2
Institut de Recherche pour le Développement (IRD), Microbes Evolution Phylogeny and Infections (MEPHI), Aix-Marseille-Université, IHU Méditerranée Infection, 13005 Marseille, France
3
Service de Pédiatrie, Hôpital de la Timone, Assistance Publique à Marseille, 13005 Marseille, France
*
Author to whom correspondence should be addressed.
Submission received: 16 March 2021 / Revised: 2 May 2021 / Accepted: 12 May 2021 / Published: 17 May 2021
(This article belongs to the Collection Feature Papers in Viral Pathogens)

Abstract

:
The current point-of-care diagnosis of enterovirus meningitis does not identify the viral genotype, which is prognostic. In this case report, more than 81% of an Echovirus 12 genome were detected and identified by metagenomic next-generation sequencing, directly from the cerebrospinal fluid collected in a 6-month-old child with meningeal syndrome and meningitis: introducing Echovirus 12 as an etiological agent of acute meningitis in the pediatric population.

1. Introduction

The highly diverse viral genus Enterovirus encompasses more than 300 genotypes [1]. Some Enterovirus members are responsible for central nervous system (CNS) infections, which clinical and epidemiological characteristics and prognosis, vary according to the precisely identified enterovirus, and some genotypes have been associated with a particular clinical severity and mortality [2]. The Enterovirus genotype is not routinely determined by molecular diagnosis assays detecting enterovirus RNA at the core and point-of-care (POC) laboratories [3,4,5,6]. Echovirus strains belonging to Enterovirus B species preferentially infect infants and young children [1,4] and are frequently involved in aseptic meningitis and encephalitis. Indeed, Enterovirus genotyping is most commonly based on partial VP1 gene Sanger sequencing using a generic standard protocol which is not commonly applied during the time of care [7].
In order to challenge Enterovirus genotyping at the POC laboratory, we herein developed a unique protocol to diagnose and genotype Enterovirus CNS infection directly from the cerebrospinal fluid (CSF) using metagenomic Next-Generation Sequencing (mNGS). This diagnosis approach is here illustrated by the diagnosis of Echovirus 12 meningitis in a child, a rarely reported situation in such setting [8,9].

2. Case Presentation and Methods

A 6-month-old girl born from twin pregnancy was admitted at the emergency department with fever (37.7 °C), cough and meningeal syndrome. She had an history of Respiratory Syncytial Virus bronchiolitis that led to hospitalization one month prior to meningitis. At the admission, CSF analysis after lumbar puncture showed a leukocyte count of 1 cell/mm3, protein at 0.14 g/L and glucose at 3.41 mmol/L. Microscopic analysis after Gram staining was negative. Investigation of the CSF at the POC laboratory using the Biofire FilmArray Meningitis/Encephalitis panel (bioMérieux, Marcy-l, Etoile, France) [6] was positive for Enterovirus.
In parallel, total RNA was manually extracted from 200 µL of CSF, following an in-house developed protocol, using QIAamp Viral RNA Mini Kit solutions (Qiagen, Hilden, Germany) for lysing, and washing steps, then total RNA was then purified using RNA specific magnetic Dynabeads (Life technology, Oslo, Norway). Briefly, 200 µL of CSF were incubated with 50 µL proteinase K (Qiagen) for 5 min at room temperature, then 300 µL AVL lysis buffer (Qiagen) were added and incubated for 15 min at room temperature. A 150-µL volume absolute ethanol (99%) were added to the lysis mix, 50-µL Dynabeads (40 mg/µL) were added and incubated for 15 min at room temperature. The Dynadeads/RNA complex washed twice with 850 µL AV1 solution (Qiagen), then two times with 450 µL AV2 solution (Qiagen) in the presence of 70% ethanol. After the second wash, Dynabeads were dried for 15 min at room temperature and eluted in a 60 µL-volume, then incubated for 3 min at 70 °C, followed by magnetic separation. Finally, 3 µL of RNaseOUT™ Recombinant Ribonuclease Inhibitor (Invitrogen, Illkirch, France were added to the purified RNA was and stored at −70 °C.
A 40-µL volume of total RNA was treated with ezDNase (Invitrogen, Illkirch, France) and concentrated with (Zymo Research, Irvine, CA, USA) kit, then eluted in 20 µL sterile water (Figure 1). The complementary DNA (cDNA) synthesis was performed using TaqMan kit according to the manufacturer protocol (Applied Biosystem, Foster City, CA, USA) in 50µL containing 19.25 µL eluted RNA, then 20 µL of the cDNA were used as a matrix for double strand synthesis, using 3 units of DNA Polymerase I, Large (Klenow) Fragment (BioLabs) in a 30 µL-volume. Double stranded DNA purified with Agencourt® AMPure beads (Invitrogen) and eluted in 17 µL of 1x-sterile Tris-EDTA solution. Finally, 1 ng of cDNA was used for metagenomics Next-Generation Sequencing (mNGS) library preparation (Supplementary Materials 3), using Illumina Nextera XT paired-end protocol (Illumina, San Diego, CA, USA), as previously described [10,11] and sequenced on iSeq 100 instrument in a single 17.5-h run providing 2 × 150-base pair (bp) long reads. The NGS generated sequences were assembled by Spades on-line software available on Galaxy/Europe bioinformatics (https://usegalaxy.eu, accessed on 17 June 2020) and mapped with CLC Genomics Workbench software version 7.5.0 (Qiagen).

3. Results and Discussion

BLAST analysis of the contigs generated by mNGS after assembling reads with Spades (https://usegalaxy.eu, accessed on 17 June 2020), yielded as best match the Echovirus 12 strain Travis 2–85 gene (GenBank accession number AF295499.1). This strain was originally isolated from a 6-year-old healthy American male, caused cytopathic effect in tissue culture, was not neutralized by poliomyelitis antiserum, and failed inducing disease in infant mice [12]. As the complete genome of this strain was not available in the GenBank database, the Echovirus 12 complete genome, prototype Travis (X77708.1) was used as reference sequence for mapping of total reads by CLC Genomics Workbench software. The iSeq sequencing generated 114.818 reads, and 76.284 (66.4%) reads could be mapped on the Echovirus 12 genome, generating 7 contigs (GenBank accession number; PRJEB39568) covering 6,127 bp, hence 81.7% of this genome (Figure 1). Phylogenetic analysis based on VP1 gene and 3D polymerase encoding genes identified that these sequences belong to an Echovirus 12 as supported by bootstrap values of 98% and 97%, respectively (Figure 2, Supplementary Materials).
Using this approach of whole genome sequencing, Echovirus 12 was obtained in one-shot protocol directly from CSF sample. To confirm the NGS result, two Echovirus 12-specific primers were designated and used to target a 291-bp Echovirus 12 genome fragment. The same strain of Echovirus 12 was identified by sequencing of the amplified fragment at 98.44% sequence identity (not published data).
Echovirus 12 has been detected in patients with diarrhea and aseptic meningitis [9,13]. In this study, we described for the first time a clinical case of Echovirus 12 meningitis diagnosed by near whole genome sequencing of an Echovirus 12 directly from CSF by mNGS. This strategy proves high sensitivity in enterovirus detection, which warrants its introduction for routine diagnostic of enterovirus meningitis in addition to viral genomic surveillance and may even be considered for POC laboratories if the very fast Oxford Nanopore Technology is used. Current routine diagnostic targets a short sequence covering around 7% of the Enterovirus genome, so that is not indicative for Enterovirus genotyping [11]. The availability of the genome instead of a gene fragment necessarily provides improved information regarding typing of the viral strain and identifying mutations and recombinations, and correlating these genotypic features with epidemiological and/or clinical ones. Another important benefit of mNGS upon qPCR or Sanger sequencing is its versatility. Indeed, this is not a targeted approach but instead it is an opened approach that can potentially detect sequences from any virus or microorganism provided these are in sufficient amount; this is thus of particular interest in cases when no infectious agent could be diagnosed. In addition, NGS cost per clinical sample is currently in the same order of magnitude than that of Sanger sequencing. Indeed, we have estimated that current cost per nucleotide is EUR 0.3 for Sanger sequencing and EUR 0.03 for NGS sequencing with Illumina technology; such cost cannot be directly extrapolated to other laboratories due to highly variable cost components, among which the commercial policy of suppliers or the infrastructure of the laboratory. Overall, the present case exemplifies the powerfulness of mNGS in the setting of the routine diagnosis of meningitis.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/pathogens10050610/s1, Supplementary Material 1: Phylogenetic tree based on 3D polymerase encoding gene analysis. Supplementary Material 2: DNase treatment and cDNA synthesis. Supplementary Material 3: NGS Library preparation.

Author Contributions

M.M.: contributed to experimental design, realization of the work, data analysis, interpretation, and writing. C.Z., performed routine analyses and collected clinical information and writing. A.M., clinical diagnostic, and CSF sampling. P.C., Bioinformatic data analysis, interpretation, and writing. M.D., contributed to critically reviewing the manuscript, data interpretation, coordinated and directed the work. All authors have read and agreed to the published version of the manuscript.

Funding

Madjid Morsli is PhD student supported by “Fondation Meditérranée Infection”. This work (Lab material and reactifs) was supported by the French Government under the Investissements d’Avenir (Investments in the Future) programme managed by the Agence Nationale de la Recherche (ANR,fr: National Agency for Research) [reference: Méditerranée Infection 10-IAHU-03]. This work was also supported by Région Le Sud (Provence Alpes-Côte d’Azur) and European funding [FEDER PA 0000320 PRIMMI].

Institutional Review Board Statement

All data were generated as part of routine work at the Assistance Publique-Hôpitaux de Marseille (Marseille University hospitals), and this study is the result of routine clinical management, validated by the IHU Méditerranée Infection Ethics Committee under the following numbers (N: 2021–004). No specific clinical sampling has been done in this study.

Informed Consent Statement

Patient tutors’ consent was waived as the study is reporting a diagnostic activity which did not require any specific clinical sample, which result did not modify the medical management of the patient; who cannot be identified through the case report.

Data Availability Statement

The identified Echovirus 12 genome available on NCBI GenBank accession number; PRJEB39568.

Acknowledgments

This study was supported by Fondation Méditerranée Infection, IHU Méditerranée Infection, Marseille, France. MM benefits from a PhD grant by Fondation Méditerranée Infection.

Conflicts of Interest

The authors have no conflict of interest to declare.

References

  1. Genus: Enterovirus—Picornaviridae—Picornavirales—International Committee on Taxonomy of Viruses (ICTV). Available online: https://talk.ictvonline.org/ictv-reports/ictv_online_report/positive-sense-rna-viruses/picornavirales/w/picornaviridae/681/genus-enterovirus (accessed on 31 December 2020).
  2. Knoester, M.; Helfferich, J.; Poelman, R.; Van Leer-Buter, C.; Brouwer, O.F.; Niesters, H.G. Twenty-nine Cases of Enterovirus-D68-associated Acute Flaccid Myelitis in Europe 2016: A Case Series and Epidemiologic Overview. Pediatr. Infect. Dis. J. 2019, 38, 16–21. [Google Scholar] [CrossRef] [Green Version]
  3. Khumalo, J.; Nicol, M.; Hardie, D.; Muloiwa, R.; Mteshana, P.; Bamford, C. Diagnostic accuracy of two multiplex real-Time poly-merase chain reaction assays for the diagnosis of meningitis in children in a resource-limited setting. PLoS ONE 2017, 12, e0173948. [Google Scholar] [CrossRef] [PubMed]
  4. Dumaidi, K.; Al-Jawabreh, A.; Samarah, F.; Zraiqi, A.; Yaseen, D. Genetic diversity of the enteroviruses detected from cerebrospinal fluid (CSF) samples of patients with suspected aseptic meningitis in northern West Bank, Palestine in 2017. PLoS ONE 2018, 13, e0202243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Archimbaud, C.; Mirand, A.; Chambon, M.; Regagnon, C.; Bailly, J.-L. Improved diagnosis on a daily basis of enterovirus meningitis using a one-step real-time RT-PCR assay. J. Med Virol. 2004, 74, 604–611. [Google Scholar] [CrossRef] [PubMed]
  6. Vincent, J.-J.; Zandotti, C.; Baron, S.; Kandil, C.; Levy, P.-Y.; Drancourt, M.; Raoult, D.; Ninove, L. Point-of-care multiplexed diagnosis of meningitis using the FilmArray® ME panel technology. Eur. J. Clin. Microbiol. Infect. Dis. 2020, 39, 1573–1580. [Google Scholar] [CrossRef] [PubMed]
  7. Oberste, M.S.; Maher, K.; Williams, A.J.; Dybdahl-Sissoko, N.; Brown, B.A.; Gookin, M.S.; Peñaranda, S.; Mishrik, N.; Uddin, M.; Pallansch, M.A. Species-specific RT-PCR amplification of human enteroviruses: A tool for rapid species identification of uncharacterized enteroviruses. J. Gen. Virol. 2006, 87, 119–128. [Google Scholar] [CrossRef]
  8. Schiff, G.M.; Stefanovic’, G.M.; Young, E.C.; Sander, D.S.; Pennekamp, J.K.; Ward, R.L. Studies of Echovirus-12 in Volunteers: Determination of Minimal Infectious Dose and the Effect of Previous Infection on Infectious Dose. J. Infect. Dis. 1984, 150, 858–866. [Google Scholar] [CrossRef]
  9. Rosenwirth, B.; Eggers, H.J. Structure and Replication of Echovirus Type 12. 1. Analysis of the Polypeptides and RNA of Echovirus 12 Particles. JBIC J. Biol. Inorg. Chem. 1978, 92, 53–60. [Google Scholar] [CrossRef] [PubMed]
  10. Diop, A.; Khelaifia, S.; Armstrong, N.; Labas, N.; Fournier, P.-E.; Raoult, D.; Million, M. Microbial culturomics unravels the halophilic microbiota repertoire of table salt: Description of Gracilibacillus massiliensis sp. nov. Microb. Ecol. Heal. Dis. 2016, 27. [Google Scholar] [CrossRef]
  11. Morsli, M.; Vincent, J.-J.; Milliere, L.; Colson, P.; Drancourt, M. Direct next-generation sequencing diagnosis of echovirus 9 meningitis, France. Eur. J. Clin. Microbiol. Infect. Dis. 2021. [Google Scholar] [CrossRef]
  12. Hammon, W.M.; Ludwig, E.H.; Pavia, R.A.; McCloskey, L.W.; Sather, G.E. Problems raised by certain ECHO viruses in the attempted laboratory detection of poliomyelitis virus infection. Ann. N. Y. Acad. Sci. 1957, 67, 304–310. [Google Scholar] [CrossRef]
  13. Liu, H.; Zhang, J.; Zhao, Y.; Zhang, H.; Lin, K.; Sun, H.; Huang, X.; Yang, Z.; Ma, S. Molecular characterization of echovirus 12 strains isolated from healthy children in China. Sci. Rep. 2018, 8, 11716. [Google Scholar] [CrossRef]
Figure 1. Echovirus 12 whole genome sequencing direct diagnosis workflow. The total duration of the handling process is less than 24 h including al steps. (A) CSF sampling by lumbar puncture, followed by specific point-of-care (POC) enterovirus (EVs) diagnosis in FilmArray® CSF direct test (BioFire Diagnostics, Salt Lake City, UT, USA). The viral RNA was extracted and purified following an in-house developed protocol, using QIAamp Viral RNA Mini Kit (Qiagen) and RNA was purified by magnetic DynaBeads (Thermo Fisher, Waltham, MA, USA). The EVs POC diagnostic was confirmed by real-time PCR using LightCycler Multiplex RNA Virus Master kit (Roche Diagnostics®, Mannheim, Germany). (B) The extracted RNA was treated by Turbo DNase (Thermo Fisher) and purified with (Zymo Research) kit. RT-cDNA synthesis reaction was performed in 50 µL-volume using kit TaqMan (Thermo Fisher), followed by double strand synthesis using DNA Polymerase I, Large (Klenow) Fragment (BioLabs). (C) The double stranded cDNA sequenced in 150-cycle iSeq Illumina instrument following Nextera NGS library preparation iSeq protocol. Finally, NGS data analysis was performed using CLC Genomics Workbench software version 7.5.0 (Qiagen), and more than 81.68% of Echovirus 12 genome was obtained directly by next-generation sequencing.
Figure 1. Echovirus 12 whole genome sequencing direct diagnosis workflow. The total duration of the handling process is less than 24 h including al steps. (A) CSF sampling by lumbar puncture, followed by specific point-of-care (POC) enterovirus (EVs) diagnosis in FilmArray® CSF direct test (BioFire Diagnostics, Salt Lake City, UT, USA). The viral RNA was extracted and purified following an in-house developed protocol, using QIAamp Viral RNA Mini Kit (Qiagen) and RNA was purified by magnetic DynaBeads (Thermo Fisher, Waltham, MA, USA). The EVs POC diagnostic was confirmed by real-time PCR using LightCycler Multiplex RNA Virus Master kit (Roche Diagnostics®, Mannheim, Germany). (B) The extracted RNA was treated by Turbo DNase (Thermo Fisher) and purified with (Zymo Research) kit. RT-cDNA synthesis reaction was performed in 50 µL-volume using kit TaqMan (Thermo Fisher), followed by double strand synthesis using DNA Polymerase I, Large (Klenow) Fragment (BioLabs). (C) The double stranded cDNA sequenced in 150-cycle iSeq Illumina instrument following Nextera NGS library preparation iSeq protocol. Finally, NGS data analysis was performed using CLC Genomics Workbench software version 7.5.0 (Qiagen), and more than 81.68% of Echovirus 12 genome was obtained directly by next-generation sequencing.
Pathogens 10 00610 g001
Figure 2. Phylogenetic analysis based on VP1 gene. The Echovirus 12 near full-length genome obtained here by mNGS (GenBank accession no LR862408; indicated by a white bold font and a black background, and 9 hit blast Echovirus 12 recovered from the NBCI GenBank nucleotide sequence database (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/nucleotide/, accessed on 29 April 2021), indicated by a bold font, were incorporated in the phylogeny reconstruction in addition to Echovirus 1 sequence recovered from GenBank database. The sequence obtained in the present study is most similar to Echovirus 12 Travis stains (X79047, X77709, AF295499) and are clustered with this sequence, confirming the BLAST result. The evolutionary history was inferred in the MEGA 7 software version 7.0.2. This analysis involved 11 nucleotide sequences. There were a total of 548 positions in the final dataset. The tree was performed by applying the neighbor-joining method and the Kimura 2-parameter method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches. The tree was drawn to scale, with branch lengths in the same units as the evolutionary distances used to infer the phylogenetic tree; the scale bars represent the corresponding number of nucleotide substitutions per site. Bootstrap values ≥ 75% are indicated at the nodes.
Figure 2. Phylogenetic analysis based on VP1 gene. The Echovirus 12 near full-length genome obtained here by mNGS (GenBank accession no LR862408; indicated by a white bold font and a black background, and 9 hit blast Echovirus 12 recovered from the NBCI GenBank nucleotide sequence database (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/nucleotide/, accessed on 29 April 2021), indicated by a bold font, were incorporated in the phylogeny reconstruction in addition to Echovirus 1 sequence recovered from GenBank database. The sequence obtained in the present study is most similar to Echovirus 12 Travis stains (X79047, X77709, AF295499) and are clustered with this sequence, confirming the BLAST result. The evolutionary history was inferred in the MEGA 7 software version 7.0.2. This analysis involved 11 nucleotide sequences. There were a total of 548 positions in the final dataset. The tree was performed by applying the neighbor-joining method and the Kimura 2-parameter method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches. The tree was drawn to scale, with branch lengths in the same units as the evolutionary distances used to infer the phylogenetic tree; the scale bars represent the corresponding number of nucleotide substitutions per site. Bootstrap values ≥ 75% are indicated at the nodes.
Pathogens 10 00610 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Morsli, M.; Zandotti, C.; Morand, A.; Colson, P.; Drancourt, M. Direct Diagnosis of Echovirus 12 Meningitis Using Metagenomic Next Generation Sequencing. Pathogens 2021, 10, 610. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10050610

AMA Style

Morsli M, Zandotti C, Morand A, Colson P, Drancourt M. Direct Diagnosis of Echovirus 12 Meningitis Using Metagenomic Next Generation Sequencing. Pathogens. 2021; 10(5):610. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10050610

Chicago/Turabian Style

Morsli, Madjid, Christine Zandotti, Aurelie Morand, Philippe Colson, and Michel Drancourt. 2021. "Direct Diagnosis of Echovirus 12 Meningitis Using Metagenomic Next Generation Sequencing" Pathogens 10, no. 5: 610. https://0-doi-org.brum.beds.ac.uk/10.3390/pathogens10050610

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop