Genome AnnouncGenome AnnouncgagaGAGenome Announcements2169-8287American Society for Microbiology1752 N St., N.W., Washington, DC243568393868863genomeA01075-1310.1128/genomeA.01075-13ProkaryotesGenome Sequences of 28 Bordetella pertussis U.S. Outbreak Strains Dating from 2010 to 2012HarvillEric T.aGoodfieldLaura L.aIvanovYuryaMeyerJessica A.aNewthChristopherbCassidayPamelacTondellaMaria LuciacLiaoPattybZimmermanJerrybMeertKathleenbWesselDavidbBergerJohnbDeanJ. Michaelb,dHolubkovRichardb,dBurrJerib,dLiuTeresab,dBrinkacLaureneKimMariaeLosadaLilianaeThe Pennsylvania State University, University Park, Pennsylvania, USANICHD Collaborative Pediatric Critical Care Research Network, Salt Lake City, Utah, USACenters for Disease Control and Prevention, Atlanta, Georgia, USAUniversity of Utah Data Coordinating Center, Salt Lake City, Utah, USAJ. Craig Venter Institute, Rockville, Maryland, USAAddress correspondence to Liliana Losada, llosada@jcvi.org.19122013Nov-Dec201316e01075-131311201326112013Copyright © 2013 Harvill et al.2013Harvill et al.This is an open-access article distributed under the terms of the Creative Commons Attribution 3.0 Unported license.

Despite the availability of highly effective vaccines, Bordetella pertussis incidence has been rapidly rising in highly vaccinated populations. Recent outbreaks have received media attention, feeding concerns about the emergence of dangerous new strains with increased virulence or that escape vaccine-induced immunity. To accelerate the study of this reemerging pathogen, we sequenced the genomes of 28 B. pertussis strains isolated during outbreaks from 2010 through 2012, making both strains and sequence data available to the scientific community.

cover-dateNovember/December 2013access-typefree
GENOME ANNOUNCEMENT

The Centers for Disease Control and Prevention classifies whooping cough as a reemerging disease, documenting increasing numbers of cases nearly every year, from hundreds per year in the 1970s to >41,000 in 2012 (1). The recent high-profile epidemics in California (in 2010) and Washington (in 2012) contribute to growing concerns and feed speculation about the ongoing evolution of Bordetella pertussis. To address these concerns, we collected 28 strains of B. pertussis from these and other outbreaks and made them available at the Biodefense and Emerging Infections Research Resources Repository for further study by others.

Here, we report the genome sequences of these 28 clinical isolates derived from whooping cough hospital cases that occurred between 2010 and 2012. Genomic DNA was prepared using a phenol-chloroform extraction method and ethanol precipitation (2). A combination of 3- or 5-kb mate pair (~30× coverage) and 100-bp Illumina paired-end reads (~50× coverage) were used for genome sequence determination. After quality trimming, all reads were used to generate assemblies with Celera assembler 6.1 (3) or Velvet assembler (4). After improvement, all genomes had between 4 and 186 scaffolds containing 25 to 285 contigs. Underlying consensus sequences and gaps were improved using custom scripts developed at the J. Craig Venter Institute (JCVI). The overall G+C content in all cases was ~67%, with genome sizes ranging from 3.83 Mb to 4.15 Mb (average, 4.0 Mb). Up to 218 copies of insertion element IS481 and 95 tripartite transcarboxylate transporter (TTT) element were found in any one B. pertussis genome, amounting to >215 kb on average. All of the B. pertussis isolates, irrespective of geographic location, belong to the same multilocus sequence type and are nearly identical to each other at the protein level (average 99% identity) in conserved regions of the chromosome. However, the B. pertussis genomes were subject to massive genome rearrangements and different gene losses that account for the majority of the diversity between strains. All the genomes were annotated using JCVI’s annotation pipeline (http://www.jcvi.org) and were predicted to have between 3,750 and 4,193 genes. As expected, mobile elements were overrepresented in these genomes compared to in published Bordetella bronchiseptica genomes due to the expansion of the repetitive elements. Over 1,000 core B. bronchiseptica genes were absent in all of the B. pertussis genomes. Many of these genes are involved in capsule biosynthesis, alternate respiration, nutrient acquisition, type VI secretion, and antibiotic resistance. All B. pertussis isolates encoded pertussis toxin, filamentous hemagglutinin, type III secretion system, adenylate cyclase, and other virulence factors.

The findings from this study suggest that currently circulating B. pertussis isolates in the United States are derived from a single genetic background. A full analysis of the virulence genes and evolution of B. pertussis is under way and will be published in a subsequent report.

Nucleotide sequence accession numbers.

The B. pertussis whole-genome shotgun projects have been deposited for each isolate at DDBJ/EMBL/GenBank as described in Table 1. The version described in this paper is the first version.

Characteristics of the 28 B. pertussis strains

B. pertussis strain nameNCBI accession no.No. of contigsN50 (bp)Total length (bp)G+C %
2250905AXSU0000000028520,9083,940,18067.76
2356847AXST0000000017931,1983,907,01967.86
2371640AXSS0000000019333,4623,904,60567.88
STO1-SEAT-0006AXSR0000000049139,3624,028,42167.77
STO1-SEAT-0007AXSQ0000000026343,1554,053,85167.78
STO1-CHLA-0011AXSP0000000025419,2744,083,34967.69
H897AXSO0000000017336,6023,938,67567.81
H918AXSN0000000060150,3974,073,71267.71
H921AXSM0000000048164,7854,064,08667.73
H939AXSL0000000043177,1964,018,81067.74
H973AXSK0000000062112,3674,014,25067.81
STO1-SEAT-0004AXSJ0000000012152,1703,839,59667.9
I002AXSI0000000063102,5224,047,56767.74
I036AXSH0000000046189,6164,057,28567.74
I176AXSG0000000033266,5264,045,24667.78
STO1-CHLA-0006AXSF0000000065134,3874,044,28367.79
CHLA-15AXSD0000000013846,1583,844,20667.89
CHLA-13AXSE000000009574,0094,029,82267.79
CHLA-20AXSC0000000013146,1403,828,49967.9
CHLA-26AXSB0000000041264,5094,045,40867.8
STO1-CHOC-0016AXSA0000000054122,8554,033,84967.79
STO1-CHOC-0017AXRZ0000000010845,1974,050,20867.8
STO1-CHOC-0018AXRY0000000012743,3914,051,57767.81
STO1-CHOC-0019AXRX0000000012052,2554,150,26267.72
STO1-CHOC-0021AXRW0000000056115,7234,027,48067.8
STO1-CHOC-0008AXRV0000000012450,5073,834,60267.9
STO1-CHOM-0012AXRU0000000037222,9724,065,44367.71
STO1-CNMC-0004AXSV0000000014946,1593,855,04267.88

Citation Harvill ET, Goodfield LL, Ivanov Y, Meyer JA, Newth C, Cassiday P, Tondella ML, Liao P, Zimmerman J, Meert K, Wessel D, Berger J, Dean JM, Holubkov R, Burr J, Liu T, Brinkac L, Kim M, Losada L. 2013. Genome sequences of 28 Bordetella pertussis U.S. outbreak strains dating from 2010 to 2012. Genome Announc. 1(6):e01075-13. doi:10.1128/genomeA.01075-13.

ACKNOWLEDGMENTS

We acknowledge the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) Collaborative Pediatric Critical Care Research Network (CPCCRN) for assisting in collecting B. pertussis isolates.

This project has been funded in part with federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under contract no. HHSN272200900007C, and by cooperative agreements from the Eunice Kennedy Shriver National Institute of Child Health and Human Development, Department of Health and Human Services (no. U10HD050096, U10HD049981, U10HD049983, U10HD050012, U10HD063108, U10HD063106, U10HD063114, U10HD049945, U10HD050009, and U01HD049934), and the National Vaccine Program Office at the United States Department of Health and Human Services. All B. pertussis isolates were collected in agreement with each institution’s IRB requirements for patient safety and confidentiality.

REFERENCES CDC 2013 Provisional pertussis surveillance report. MMWR Morb. Mortal. Wkly. Rep. 61:719732 ChomczynskiPSacchiN 1987 Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal. Biochem. 162:1561592440339 MillerJRDelcherALKorenSVenterEWalenzBPBrownleyAJohnsonJLiKMobarryCSuttonG 2008 Aggressive assembly of pyrosequencing reads with mates. Bioinformatics 24:2818282418952627 ZerbinoDRBirneyE 2008 Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 18:82182918349386