# Complete Genome Sequences of Bordetella pertussis Clinical Isolate FR5810 and Reference Strain Tohama from Combined Oxford Nanopore and Illumina Sequencing

Valérie Bouchez, Sarah Louise Baines, Sophie Guillot, Sylvain Brisse

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

## To cite this version:

Valérie Bouchez, Sarah Louise Baines, Sophie Guillot, Sylvain Brisse. Complete Genome Sequences of Bordetella pertussis Clinical Isolate FR5810 and Reference Strain Tohama from Combined Oxford Nanopore and Illumina Sequencing. Microbiology Resource Announcements, 2018, 7 (19), pp.e01207-18. 10.1128/MRA.01207-18. pasteur-04838187

HAL Id: pasteur-04838187

https://pasteur.hal.science/pasteur-04838187v1

Submitted on 14 Dec 2024

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Complete Genome Sequences of *Bordetella pertussis* Clinical Isolate FR5810 and Reference Strain Tohama from Combined Oxford Nanopore and Illumina Sequencing

Valérie Bouchez,a,b Sarah Louise Baines,c Sophie Guillot,a,b Sylvain Brissea,b

**ABSTRACT** Here, we describe the complete genome sequences of two *Bordetella pertussis* strains, FR5810, a clinical isolate recovered from the respiratory tract of an infant, and Tohama, a key reference strain for the species. Sequences were obtained using a hybrid approach combining Oxford Nanopore Technologies MinION and Illumina NextSeq 500 sequence data.

Bordetella pertussis is a Gram-negative bacterial pathogen responsible for the vaccine-preventable disease whooping cough. Whole-genome sequencing (WGS) of clinical isolates of B. pertussis is useful to understand transmission at fine typing resolution scales. Circularization of the chromosomal sequence is needed to understand the dynamics of the numerous insertion sequences that may contribute to genetic adaptation of the pathogen populations toward escape of vaccine-induced immunity.

Here, we report the complete genome sequences obtained for clinical isolate FR5810 (1) and reference strain Tohama (2), using a hybrid assembly strategy combining long- and short-sequence reads. FR5810 was isolated in 2014 from the respiratory tract of a 44-day-old baby and was identified using API20E (bioMérieux) on the basis of morphological characteristics and biochemical properties as *B. pertussis*.

Isolates were first grown on Bordet Gengou Agar (Becton, Dickinson) supplemented with 15% defibrinated horse blood at 36°C for 3 days and then subcultivated for 24 h at 36°C. Bacteria were suspended in physiological water and pelleted for DNA extraction. For Illumina sequencing, bacteria were lysed in phosphate-buffered saline (PBS) 1×, with lysis buffer (Roche Diagnostics) and proteinase K (Roche Diagnostics), by incubating at 65°C for 10 min and then at 95°C for 10 min. Libraries were constructed using the Nextera XT DNA library preparation kit (Illumina), and sequenced on a NextSeq 500 platform (Illumina) using 2 × 150-bp chemistry. For Nanopore sequencing, genomic DNA was prepared by phenol-chloroform extraction using Phase Lock Gel tubes (Qiagen GmbH) and purified by ethanol precipitation. Libraries were prepared using a 1D ligation sequencing kit (SQK-LSK-108, Nanopore), without fragmentation, and sequenced on the MinION using a FLO-MIN-106 flow cell (Nanopore), following the protocol for 1D genomic DNA (gDNA) long reads without BluePippin (Nanopore).

NextSeq 500 sequencing generated 3,691,944 paired reads for FR5810 and 3,037,288 paired reads for Tohama. MinION sequencing generated 237,860 reads (5.2 Gb) for FR5810 and 49,416 reads (1.3 Gb) for Tohama. The average long-read read lengths were 21,785 bp for FR5810 and 26,212 bp for Tohama.

A hybrid Nanopore-Illumina de novo assembly was performed using Unicycler v0.4.4 (3) run under a normal assembly mode with default parameters, including a polishing

**Received** 30 August 2018 **Accepted** 19
October 2018 **Published** 15 November 2018

**Citation** Bouchez V, Baines SL, Guillot S, Brisse S. 2018. Complete genome sequences of *Bordetella pertussis* clinical isolate FR5810 and reference strain Tohama from combined Oxford Nanopore and Illumina sequencing. *Microbiol Resour Announc* 7:e01207-18. https://doi.org/10.1128/MRA.01207-18.

**Editor** Vincent Bruno, University of Maryland School of Medicine

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

Address correspondence to Sylvain Brisse, sylvain.brisse@pasteur.fr.

<footnote><sup>a</sup>Biodiversity and Epidemiology of Bacterial Pathogens Unit, Institut Pasteur, Paris, France</footnote>

<footnote><sup>b</sup>National Reference Center for Whooping Cough and Other Bordetella Infections, Institut Pasteur, Paris, France</footnote>

<footnote><sup>c</sup>Department of Microbiology and Immunology, The University of Melbourne at the Peter Doherty Institute for Infection and Immunity, Melbourne, Australia</footnote>
<figure><img src="image_3.png" /></figure>

step with Pilon. This produced one circular chromosome of 4,102,412 bp for Tohama (67.7% G+C content). For FR5810, a long-read only assembly was performed with the first 30 fastq files (out of the 59 files produced during the run), corresponding to 124,000 reads, using Unicycer v0.4.4 (3). This produced one circular chromosome of 4,108,173 bp (67.7% G+C content). This sequence was then polished by mapping the Illumina short reads against the chromosome using Snippy v3.2 (https://github.com/tseemann/snippy), correcting variants, and then remapping in an iterative process until no variants were identified. Both chromosomes were annotated using Prokka v1.9 (4).

The obtained sequence of strain Tohama was compared to its initial sequence (NCBI accession number NC_002929) using progressiveMauve snapshot-2015-02-25 (5), which revealed a very good synteny. The clinical isolate FR5810 was characterized as harboring the following vaccine antigen alleles: ptxP3 and ptxA1 for the promoter and subunit 1 of pertussis toxin, respectively, alleles fim2-1 and fim3-1 for fimbria-encoding genes, and a PRN2-type pertactin. In addition, comparison to the Tohama genome sequence using progressiveMauve snapshot-2015-02-25 (5) revealed genome insertions, deletions, and multiple genomic rearrangements, as previously observed (6, 7).

This work demonstrates, for the first time to our knowledge, the feasibility of achieving high-quality *Bordetella pertussis* genome sequences from a hybrid approach combining Nanopore and Illumina sequencing.

**Data availability.** The complete genome sequences were deposited in GenBank under the accession numbers CP031787 and CP031788, corresponding to the sample accession numbers SAMN09862822 and SAMN09862835, respectively. Illumina raw reads files were previously published (1) under BioSample accession numbers ERS1869884 and ERS1869859, and Oxford Nanopore Technologies MinION raw reads files were deposited in SRA under the accession numbers SRX4744160 and SRX4744159. The versions described in this paper are the first versions.

## ACKNOWLEDGMENTS

This study was supported financially by the INCEPTION program (French Government Investissement d'Avenir grant ANR-16-CONV-0005). S.L.B. is supported by an Australian Government Research Training Program Scholarship and a Victorian Fellowship award by the Victorian State Government, Australia.

Illumina NextSeq 500 sequencing was performed at the Mutualized Platform for Microbiology (P2M) of Institut Pasteur (Paris, France). MinION sequencing was performed in the Biodiversity and Epidemiology of Bacterial Pathogens Unit of Institut Pasteur.

## REFERENCES

1. Bouchez V, Guglielmini J, Dazas M, Landier A, Toubiana J, Guillot S, Criscuolo A, Brisse S. 2018. Genomic sequencing of *Bordetella pertussis* for epidemiology and global surveillance of whooping cough. Emerg Infect Dis 24:988-994. https://doi.org/10.3201/eid2406.171464.

2. Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, Holden MTG, Churcher CM, Bentley SD, Mungall KL, Cerdeño-Tárraga AM, Temple L, James K, Harris B, Quail MA, Achtman M, Atkin R, Baker S, Basham D, Bason N, Cherevach I, Chillingworth T, Collins M, Cronin A, Davis P, Doggett J, Feltwell T, Goble A, Hamlin N, Hauser H, Holroyd S, Jagels K, Leather S, Moule S, Norberczak H, O'Neil S, Ormond D, Price C, Rabbinowitsch E, Rutter S, Sanders M, Saunders D, Seeger K, Sharp S, Simmonds M, Skelton J, Squares R, Squares S, Stevens K, Unwin L, Whitehead S, Barrell BG, Maskell DJ. 2003. Comparative analysis of the genome sequences of *Bordetella pertussis*, *Bordetella parapertussis* and *Bordetella bronchiseptica*. Nat Genet 35:32-40. https://doi.org/10.1038/ng1227.

3. Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycer: resolving bacterial

genome assemblies from short and long sequencing reads. PLoS Comput Biol 13:e1005595. https://doi.org/10.1371/journal.pcbi.1005595.

4. Seemann T. 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068-2069. https://doi.org/10.1093/bioinformatics/btu153.

5. Darling AE, Mau B, Perna NT. 2010. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 5:e11147. https://doi.org/10.1371/journal.pone.0011147.

6. Weigand MR, Peng Y, Loparev V, Batra D, Bowden KE, Burroughs M, Cassiday PK, Davis JK, Johnson T, Juieng P, Knipe K, Mathis MH, Pruitt AM, Rowe L, Sheth M, Tondella ML, Williams MM. 2017. The history of *Bordetella pertussis* genome evolution includes structural rearrangement. J Bacteriol 199:e00806-16. https://doi.org/10.1128/JB.00806-16.

7. Bowden KE, Weigand MR, Peng Y, Cassiday PK, Sammons S, Knipe K, Rowe LA, Loparev V, Sheth M, Weening K, Tondella ML, Williams MM. 2016. Genome structural diversity among 31 *Bordetella pertussis* isolates from two recent U.S. whooping cough statewide epidemics. mSphere 1:e00036-16. https://doi.org/10.1128/mSphere.00036-16.