# Complete Genome Sequences of Monkeypox Virus from a French Clinical Sample and the Corresponding Isolated Strain, Obtained Using Nanopore Sequencing

Charlotte Balière, Véronique Hourdel, Aurelia Kwasiborski, Quentin Grassin, Maxence Feher, Damien Hoinard, Jessica Vanhomwegen, Fabien Taieb, Paul-Henri Consigny, Jean-Claude Manuguerra, et al.

► To cite this version:

Charlotte Balière, Véronique Hourdel, Aurelia Kwasiborski, Quentin Grassin, Maxence Feher, et al. Complete Genome Sequences of Monkeypox Virus from a French Clinical Sample and the Corresponding Isolated Strain, Obtained Using Nanopore Sequencing. Microbiology Resource Announcements, 2023, 12 (4), pp.e0000923. 10.1128/mra.00009-23. hal-04127243

HAL Id: hal-04127243

https://hal.science/hal-04127243

Submitted on 30 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Complete Genome Sequences of Monkeypox Virus from a French Clinical Sample and the Corresponding Isolated Strain, Obtained Using Nanopore Sequencing

Charlotte Balière,ª Véronique Hourdel,ª Aurelia Kwasiborski,ª Quentin Grassin,ª Maxence Feher,ª Damien Hoinard,ª Jessica Vanhomwegen,ª Fabien Taieb,ᵇ Paul-Henri Consigny,ᵇ Jean-Claude Manuguerra,ª India Leclercq,ª Christophe Batéjat,ª Valérie Caroª

**ABSTRACT** We report the whole-genome sequences of a monkeypox virus from the skin lesion of a French patient and the corresponding isolated viral strain. Both viral genomic sequences were successfully obtained by applying shotgun metagenomics using the Oxford Nanopore Technologies sequencing approach.

Mopox, an emerging disease of zoonotic origin, is caused by *Monkeypox virus* (MPXV), belonging to the genus *Orthopoxvirus* of the family *Poxviridae*. Traditionally endemic to Western and Central African countries, since May 2022, MPXV has caused the largest multicountry outbreak described so far (1).

As of 3 January 2023, a total of 83,974 confirmed monkeypox cases and 75 deaths had been reported to the World Health Organization, with 4,114 cases in France (2). MPXV harbors a large double-stranded DNA (dsDNA) genome of around 200 kb (3), which has already been fully sequenced using the MinION device (4). Here, as part of a global MPX surveillance effort, we directly sequenced a French clinical sample and the corresponding isolated strain.

A swab was collected from a pustular lesion on the leg of a 36-year-old French man who presented for an infectious disease consultation at the Medical Center, Institut Pasteur, in June 2022. Vero E6 cells were inoculated with the clinical specimen and incubated for 3 days, and the culture supernatant was harvested and tested for the presence of MPXV by PCR. The clinical sample and the corresponding isolate were found positive with respective threshold cycle values of 20 and 15, using a generic MPXV real-time PCR assay (G2R_G), as previously described (5). The lab investigations described in this article were carried out in accordance with the General Data Protection Regulation (regulation [EU] 2016/679 and directive 95/46/EC) and the French data protection law (law 78-17 on 6 January 1978 and decree 2019-536 on 29 May 2019), which do not require a review by an ethics committee for the secondary use of samples collected for healthcare purposes. In such cases, secondary use for research is authorized, as the persons have been informed of such secondary use (article L.1211-2 of the French public health code).

Genomic DNA (gDNA) was extracted from 150 µL each of the clinical sample (CIBU-220035) and isolate (CMIP-P0) using the NucleoSpin Dx virus kit (Macherey-Nagel). Using the Qubit dsDNA high-sensitivity assay kit (Invitrogen), the concentration of gDNA was estimated at 0.261 ng/µL and 6.18 ng/µL for the clinical sample and the isolated strain, respectively.

For Nanopore sequencing, libraries were prepared using the rapid PCR barcoding kit (SQK-RPB004; Oxford Nanopore Technologies). The library for CMIP-P0 was prepared from 5 ng of gDNA, but due to the low gDNA concentration of the CIBU-220035 sample, its library was prepared from the maximal volume input (i.e., 3 µL). Each library was

**Editor** Simon Roux, DOE Joint Genome Institute

**Copyright** © 2023 Balière et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Valérie Caro, valerie.caro@pasteur.fr.

The authors declare no conflict of interest.

**Received** 20 January 2023

**Accepted** 13 February 2023

**Published** 27 March 2023

<footnote>ªInstitut Pasteur, Université Paris Cité, Environment and Infectious Risks Unit, Laboratory for Urgent Response to Biological Threats (CIBU), Paris, France</footnote>

<footnote>ᵇInstitut Pasteur, Medical Center (CMIP), Paris, France</footnote>
<table><caption>TABLE 1 Nanopore sequencing data for CIBU-220035 and CMIP-P0</caption><thead><tr><th rowspan="2">Feature</th><th colspan="2">Data for strain:</th></tr><tr><th>CIBU-220035</th><th>CMIP-P0</th></tr></thead><tbody><tr><td>Run data (Gb)</td><td>94.5</td><td>96.05</td></tr><tr><td>No. of total reads (million)</td><td>3.66</td><td>3.92</td></tr><tr><td>No. of mapped reads (million)</td><td>13,597</td><td>511,657</td></tr><tr><td>Mean length of reads (bp)</td><td>2,158</td><td>1,814</td></tr><tr><td>Mean quality score</td><td>11.99</td><td>11.68</td></tr><tr><td>Avg depth (×)</td><td>145</td><td>4,456</td></tr><tr><td>Avg coverage 20× (%)</td><td>99.6</td><td>99.9</td></tr><tr><td>Consensus sequence length (bp)</td><td>197,205</td><td>197,205</td></tr><tr><td>GC content (%)</td><td>33</td><td>33</td></tr></tbody></table>

loaded independently onto an R9.4 flow cell (FLO-MIN106) and sequenced on an Mk1C device for 24 h, yielding reads with a Phred quality score of 8. Default parameters were used for all software. MinKNOW v19.12.2 software was used to collect the raw data and the Guppy v3.4 tool to perform base calling. Subsequently, the “What’s in My Pot?” (WIMP) workflow (FASTQ [human + viral] v2021.11.26) was launched for real-time species classification, allowing the assignation of 0.46% and 18.5% of reads to monkeypox virus strain Zaire-96-I-16 (NCBI taxonomy ID 619591) for CIBU-220035 and CMIP-P0, respectively.

Using Minimap2 v2.17-r941 (6), Medaka v1.0.3, and bcftools v1.10.2, a consensus sequence of 197,205 bp was generated for CIBU-220035 and CMIP-P0 by mapping the reads to the reference genome found under GenBank accession number MT903345, with a high coverage depth along the total length (Table 1), and manually curating them using Tablet v1.21 (7). Both viral genomic sequences belonged to lineage B.1 and clade IIb, according to Nextclade (https://clades.nextstrain.org) (8). Comparison of the genomic sequences obtained for the clinical specimen and its respective isolate showed 99% nucleotide identity, with 5 nucleotide substitutions, resulting in 1 amino acid change (a G-to-A change at position 193366 [G193366A] encoded by the ankyrin-like gene) compared to the reference genome MT903345. Despite the limitations of the genome quality obtained using MinION, this result would suggest that the culture on Vero E6 cells induced some molecular changes, immediately after the primary viral isolation.

We report the successful and rapid sequencing of MPXV from a pustular lesion specimen and its corresponding strain using a direct sequencing approach with the Nanopore platform, despite the low-viral load sample.

**Data availability.** The consensus genome sequences were deposited in the GISAID EpiPox database (accession numbers EPI_ISL_16260351 and EPI_ISL_16260402) and at GenBank (accession numbers OQ249660 and OQ249661). The raw reads were deposited in the NCBI Sequence Read Archive under accession number PRJNA918303.

## ACKNOWLEDGMENT

This project received funding from the French Institut Carnot “Microbes-Santé” in collaboration with the “France Futur Elevage” Carnot of INRAE Institute, in the framework of the FIELD project.

## REFERENCES

1. World Health Organization. 2022–2023. Mpox (monkeypox) outbreak 2022. https://www.who.int/emergencies/situations/monkeypox-oubreak-2022. Accessed 7 November 2022.

2. Our World in Data, World Health Organization. 2022–2023. Mpox (monkeypox) data explorer: visualizing the data produced by the World Health Organization on the 2022 mpox outbreak. https://ourworldindata.org/explorers/monkeypox. Retrieved 7 November 2022.

3. Isidro J, Borges V, Pinto M, Sobral D, Santos JD, Nunes A, Mixão V, Ferreira R, Santos D, Duarte S, Vieira L, Borrego MJ, Núcio S, de Carvalho IL, Pelerito A, Cordeiro R, Gomes JP. 2022. Phylogenomic characterization

and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus. Nat Med 28:1569–1572. https://doi.org/10.1038/s41591-022-01907-y.

4. Vandenbogaert M, Kwasiborski A, Gonofio E, Descorps-Declère S, Selekon B, Nkili Meyong AA, Ouilibona RS, Gessain A, Manuguerra J-C, Caro V, Nakoune E, Berthet N. 2022. Nanopore sequencing of a monkeypox virus strain isolated from a pustular lesion in the Central African Republic. Sci Rep 12:10768. https://doi.org/10.1038/s41598-022-15073-1.

5. Li Y, Zhao H, Wilkins K, Hughes C, Damon IK. 2010. Real-time PCR assays for the specific detection of monkeypox virus West African and Congo Basin
strain DNA. J Virol Methods 169:223-227. https://doi.org/10.1016/j.jviromet.2010.07.012.

6. Li H. 2018. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34:3094-3100. https://doi.org/10.1093/bioinformatics/bty191.

7. Milne I, Stephen G, Bayer M, Cock PJA, Pritchard L, Cardle L, Shaw PD,

Marshall D. 2013. Using Tablet for visual exploration of second-generation sequencing data. Brief Bioinform 14:193-202. https://doi.org/10.1093/bib/bbs012.

8. Aksamentov I, Roemer C, Hodcroft EB, Neher RA. 2021. Nextclade: clade assignment, mutation calling and quality control for viral genomes. J Open Source Softw 6:3773. https://doi.org/10.21105/joss.03773.