Volume 26, Number 6—June 2020
Fatal Encephalitis Caused by Cristoli Virus, an Emerging Orthobunyavirus, France
We report the discovery of a new orthobunyavirus, Cristoli virus, by means of shotgun metagenomics. The virus was identified in an immunodepressed patient with fatal encephalitis. Full-length genome sequencing revealed high-level expression of a virulence factor, possibly explaining the severity of the infection. The patient’s recent history suggests circulation in France.
The Orthobunyavirus genus of the Peribunyaviridae family contains numerous viruses, usually transmitted by mosquitoes (1). New members are regularly discovered through mosquito screening campaigns, but most of them are not pathogenic (2–4). Global warming implies changes in the distribution of their vectors, possibly exposing humans to the onset of new diseases.
Broad-spectrum diagnostic tests, such as metagenomics methods, are useful to discover new viruses, especially in patients with encephalitis of unknown etiology (2,5). We have developed an original method based on shotgun metagenomics (MetaMIC) for the diagnosis of bacterial, viral, fungal, and parasitic infections from any human fluid or tissue (6). We used it to discover a previously unknown member of the Orthobunyavirus genus of the Peribunyaviridae family, Cristoli virus, in a patient with fatal encephalitis. The study was approved by the Créteil Institutional Review Board (Créteil, France).
A 58-year-old woman living in the Paris, France, area was hospitalized in September 2018 for isolated fever resistant to amoxicillin/clavulanic acid. She had a history of complete remission of non-Hodgkin lymphoma, decompensated cirrhosis related to autoimmune hepatitis treated with sirolimus, possibly congenital hypogammaglobulinemia associated with B-lymphopenia, and episodes of lower respiratory tract and digestive infections in the years preceding her admission. She last traveled to Italy and on a Mediterranean cruise in summer 2017. Although her fever declined on ofloxacin, she was admitted again in October 2018 for deterioration of her general health, anorexia, and psychomotor retardation. Her neurologic symptoms worsened during the following 6 months. Multiple electroencephalograms showed nonspecific slowing of background activity. Magnetic resonance imaging revealed limbic and cerebellar abnormalities compatible with inflammatory or infectious encephalitis (Figure). Results of testing for multiple cerebrospinal fluid samples showed only an increase in α interferon; no microbial agent was identified. A cerebral biopsy was performed in March 2019 for shotgun metagenomics testing. The patient’s condition continued to deteriorate, and she died in the intensive care unit on March 27, 2019.
We ground the cerebral biopsy material in a sterile isotonic solution. We performed preextraction using bead beating combined with chemical cell disruption, then a combined DNA/RNA extraction using QIAsymphony (QIAGEN, https://www.qiagen.com). We tested a negative environmental control (isotonic sterile solution) in parallel. We prepared DNA libraries with Nextera XT DNA, and RNA libraries with RNA Human RiboZero TruSeq Stranded Total RNA Library Prep Kit (Illumina, https://www.illumina.com). We pair-end sequenced libraries, 2 × 150 bp, with High Output Kit version 2 on a NextSeq500 device (Illumina). We analyzed the sequences with MetaMIC software (6).
Identification of nonhuman sequences uses different databases, including cleaned National Center for Biotechnology Information nucleotide and nonredundant databases (GenBank release 229, December 2018) that contain all known microorganisms. For new species, we conducted verification using an algorithm combining de novo assembly with Meta-SPAdes version 3.12.0 (http://cab.spbu.ru/software/meta-spades) and contig alignment with the ViPR database (https://www.viprbrc.org). After manual curation of full-length viral sequences, we performed alignment and phylogenetic analyses with the closest known viral species using Muscle version 3.8.31 (http://www.drive5.com/muscle) and IQ-Tree version 188.8.131.52 (http://www.iqtree.org). We used the general time reversible plus gamma 4 plus invariate sites model of nucleotide substitution and 10,000 full maximum-likelihood bootstrap replicates.
Shotgun metagenomics revealed the presence of sequences from a new member of the Orthobunyavirus genus of the Peribunyaviridae family; the sequences were not present in the negative environmental control. This new virus was named Cristoli virus after the city of Créteil where it was described. We reconstructed the full-length sequence of the viral genome de novo in silico. The large (L), medium (M), and small (S) segments were individualized and could be aligned with similar sequences from known orthobunyaviruses (1) (Appendix 1 Figure 1, panels A, B). Cristoli virus segregated within serogroup Turlock; the closest related virus was Umbre virus. These viruses cluster on a different branch from all other known orthobunyaviruses with a 100% bootstrap (Appendix 1 Figure 1, panel C). Cristoli virus differed from Umbre virus by 11.4% of L, 12.3% of M, and 7.0% of S segments (differences with other orthobunyaviruses) (Appendix 2).
We identified a 6,783-bp open reading frame (ORF), coding for a protein of 2,261 aa in the L segment, which, by analogy with other orthobunyaviruses, corresponded to the RNA-dependent RNA polymerase (Appendix 1 Figure 1). The conserved H….PD….DxK orthobunyavirus catalytic domain (7) was present between positions 37 and 101, identical to that of Umbre virus. Umbre and Cristoli viruses differed by 44 aa in the RNA polymerase (Table).
In segment M, we identified a 4,336-bp ORF coding for a protein of 1,445 aa, corresponding to a polyprotein processed into the 2 envelope glycoproteins n and c (Gn and Gc) and nonstructural protein m (NSm) in orthobunyaviruses (Appendix 1 Figure 1). Cristoli and Umbre viruses differed by 78 aa in the M segment, including 2 deletions (P758del and R759del) (Table). R759del corresponds to a trypsin cleavage site used for bunyavirus characterization (8). The conserved RxxR cleavage motif at the C-terminus of the Gn domain was conserved in Cristoli virus at positions 298–301. All conserved cysteine positions (9) were also conserved.
We identified 2 overlapping ORFs in the S segment, one corresponding to the nucleocapsid (N) protein (714 bp, 238 aa), the other to nonstructural protein s (NSs, 240 bp, 80 aa). We found few differences with Umbre virus for these proteins (Table).
The amounts of viral RNAs we semiquantified by shotgun metagenomics were 5.4 log viral genomes/mg of biopsy for segment L, log viral genomes/mg of biopsy 5.0 for segment M, and 6.2 log viral genomes/mg of biopsy for segment S. For confirmation, we designed specific reverse transcription PCRs for Cristoli virus segments L, M, and S. Brain extracts were positive for all three, with amplicon sequences identical to those in the full-length genome (Appendix 1 Figure 2). We deposited the annotated genome sequence into Genbank (accession nos. MN488996 [S segment], MN488997 [M segment], and MN488998 [L segment]).
The etiologic diagnosis of encephalitis is difficult because >25% of cases remain of unknown cause (5). Diagnostic methods without prior knowledge of the microorganisms sought, such as metagenomics, detect infectious agents when other techniques have failed (10). We used an original shotgun metagenomics approach and in-house software MetaMIC to discover a new virus belonging to the Orthobunyavirus genus of the Peribunyaviridae family in an immunodepressed patient with slowly progressive fatal encephalitis. Cristoli virus is close to Umbre virus, a member of the Turlock serogroup not previously associated with human disease. Because the patient had not traveled during the year preceding the onset of her symptoms, our findings suggest that Cristoli virus is endemic to France. Like other members of its family, it could have been transmitted through a mosquito bite. Changes in the distribution of mosquito species related to climate change may make such transmissions more frequent in the future. In this context, the metagenomics method we developed is particularly useful to detect unusual or emerging infectious agents in patients with unassigned infectious syndromes.
Full-length sequence analysis of Cristoli virus revealed a classical Peribunyaviridae genome organization with conserved motifs. The genome segment coding for the NSs protein was found in an amount 10 times higher than that of the other segments. Because NSs has been described as a virulence factor (11,12), its high-level expression in a context of immune depression could explain the unusual severity of the disease. If efficacious antiviral treatments become available (13), a rapid etiologic diagnosis of orthobunyavirus infections will become necessary. The shotgun metagenomics method we developed fulfills the criteria for such routine diagnosis.
Dr. Rodriguez is an assistant professor in the Department of Microbiology at Henri Mondor Hospital, University of Paris-Est-Créteil, Créteil, France, and the head of the institution’s genomics platform. He has implemented a clinical metagenomics technique in diagnostic routine and is working on the exploration of complex infectious diseases using this tool.
The authors declare that the planning, conduct, and reporting of the case was in line with the Declaration of Helsinki and French laws for biomedical research. It only involved the use of archived, residual samples that were sent for routine diagnosis and sequence characterization with nonopposition for leftover sample to be used in other assays.
C.R., G.G., V.D. and J.-M.P. are co-inventors of a patent protecting the shotgun metagenomics method and MetaMIC software.
Author contributions: C.R and J.-M.P conceived the study. C.R., G.G., M.N. and V.D. performed sequence analyses, developed the algorithms, and ran computational analyses. S.D., S.S, S.B., D.B., and V.C. provided and interpreted the clinical, radiologic, and virologic data. L.P. developed the confirmation RT-PCR methods. C.R., P.-L.W., and J.-M.P. analyzed the findings and wrote the manuscript. All authors discussed the results and contributed to the final manuscript.
- Elliott RM. Orthobunyaviruses: recent genetic and structural insights. Nat Rev Microbiol. 2014;12:673–85. DOIPubMedGoogle Scholar
- Li D. A highly pathogenic new bunyavirus emerged in China. Emerg Microbes Infect. 2013;2:
e1. DOIPubMedGoogle Scholar
- Huang YS, Higgs S, Vanlandingham DL. Emergence and re-emergence of mosquito-borne arboviruses. Curr Opin Virol. 2019;34:104–9. DOIPubMedGoogle Scholar
- Huang B, Allcock R, Warrilow D. Newly characterized arboviruses of northern Australia. Virol Rep. 2016;6:11–7. DOIGoogle Scholar
- Granerod J, Ambrose HE, Davies NW, Clewley JP, Walsh AL, Morgan D, et al.; UK Health Protection Agency (HPA) Aetiology of Encephalitis Study Group. Causes of encephalitis and differences in their clinical presentations in England: a multicentre, population-based prospective study. Lancet Infect Dis. 2010;10:835–44. DOIPubMedGoogle Scholar
- Rodriguez C, Jary A, Hua C, Woerther PL, Bosc R, Desroches M, et al.; Multidisciplinary Necrotizing Fasciitis Study Group. Pathogen identification by shotgun metagenomics of patients with necrotizing soft-tissue infections. Br J Dermatol. 2019;
bjd.18611. DOIPubMedGoogle Scholar
- Reguera J, Weber F, Cusack S. Bunyaviridae RNA polymerases (L-protein) have an N-terminal, influenza-like endonuclease domain, essential for viral cap-dependent transcription. PLoS Pathog. 2010;6:
e1001101. DOIPubMedGoogle Scholar
- Fazakerley JK, Gonzalez-Scarano F, Strickler J, Dietzschold B, Karush F, Nathanson N. Organization of the middle RNA segment of snowshoe hare Bunyavirus. Virology. 1988;167:422–32. DOIPubMedGoogle Scholar
- Briese T, Rambaut A, Lipkin WI. Analysis of the medium (M) segment sequence of Guaroa virus and its comparison to other orthobunyaviruses. J Gen Virol. 2004;85:3071–7. DOIPubMedGoogle Scholar
- Wilson MR, Sample HA, Zorn KC, Arevalo S, Yu G, Neuhaus J, et al. Clinical metagenomic sequencing for diagnosis of meningitis and encephalitis. N Engl J Med. 2019;380:2327–40. DOIPubMedGoogle Scholar
- Bridgen A, Weber F, Fazakerley JK, Elliott RM. Bunyamwera bunyavirus nonstructural protein NSs is a nonessential gene product that contributes to viral pathogenesis. Proc Natl Acad Sci U S A. 2001;98:664–9. DOIPubMedGoogle Scholar
- Weber F, Bridgen A, Fazakerley JK, Streitenfeld H, Kessler N, Randall RE, et al. Bunyamwera bunyavirus nonstructural protein NSs counteracts the induction of alpha/beta interferon. J Virol. 2002;76:7949–55. DOIPubMedGoogle Scholar
- Ter Horst S, Conceição-Neto N, Neyts J, Rocha-Pereira J. Structural and functional similarities in bunyaviruses: Perspectives for pan-bunya antivirals. Rev Med Virol. 2019;29:
e2039. DOIPubMedGoogle Scholar
TableCite This Article
Original Publication Date: May 04, 2020
1These authors contributed equally to this article.
Table of Contents – Volume 26, Number 6—June 2020
|EID Search Options|
|Advanced Article Search – Search articles by author and/or keyword.|
|Articles by Country Search – Search articles by the topic country.|
|Article Type Search – Search articles by article type and issue.|
Please use the form below to submit correspondence to the authors or contact them at the following address:
Christophe Rodriguez, Department of Microbiology, Hôpital Henri Mondor, 51 avenue du Maréchal de Lattre de Tassigny, 94010 Créteil CEDEX, France