Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 22, Number 3—March 2016
Research

Whole-Genome Sequencing to Determine Origin of Multinational Outbreak of Sarocladium kiliense Bloodstream Infections

Figures
Tables
Article Metrics
36
citations of this article
EID Journal Metrics on Scopus
Author affiliations: Centers for Disease Control and Prevention, Atlanta, Georgia, USA (K.A. Etienne, R.M. Smith, S. Vallabhaneni, A.P. Litvintseva); Georgia Institute of Technology, Atlanta (K.A. Etienne, F.O. Vannberg); Translational Genomics Research Institute, Flagstaff, Arizona, USA (C.C. Roe, J. Gillece, D. Lemmer, E.M. Driebe, D.M. Engelthaler); Instituto Nacional de Salud, Bogota, Colombia (C. Duarte, P. Escandón, E. Castañeda); Universidad del Rosario, Bogota (B.L. Gómez); Corporación para Investigaciones Biológicas, Medellín, Colombia (B.L. Gómez, C. de Bedout, L.F. López); Instituto de Salud Publica de Chile, Santiago, Chile (V. Salas, L.M. Hederra, J. Fernández, P. Pidal, J.C. Hormazabel); Ministry of Health, Santiago (F. Otaíza-O’Ryan)

Cite This Article

Abstract

We used whole-genome sequence typing (WGST) to investigate an outbreak of Sarocladium kiliense bloodstream infections (BSI) associated with receipt of contaminated antinausea medication among oncology patients in Colombia and Chile during 2013–2014. Twenty-five outbreak isolates (18 from patients and 7 from medication vials) and 11 control isolates unrelated to this outbreak were subjected to WGST to elucidate a source of infection. All outbreak isolates were nearly indistinguishable (<5 single-nucleotide polymorphisms), and >21,000 single-nucleotide polymorphisms were identified from unrelated control isolates, suggesting a point source for this outbreak. S. kiliense has been previously implicated in healthcare-related infections; however, the lack of available typing methods has precluded the ability to substantiate point sources. WGST for outbreak investigation caused by eukaryotic pathogens without reference genomes or existing genotyping methods enables accurate source identification to guide implementation of appropriate control and prevention measures.

Despite modern advances in technology to control fungal contamination in clinical settings, fungi are continuously implicated in clusters or outbreaks of infections, particularly among immunosuppressed patients (1). The sources of fungal nosocomial outbreaks often are difficult to assess because of the widespread prevalence of fungi in the environment. Specifically, differentiating between fungal infections originating from a single contaminated point source and those independently acquired from the environment frequently is difficult. Molecular typing methods for discriminating strains have been an essential tool to identify potential source(s) of fungal infections in outbreaks. Small-scale DNA-based typing methods, such as variable-number tandem-repeat typing and multilocus sequence typing (MLST), use genomic similarity to assist in determining epidemiologic relatedness of fungal isolates in an outbreak investigation (2). However, the robustness and accuracy of such genotyping tools depend largely on the discriminatory power of the genotyping method and availability of reference data, which often are inconsistent or incomplete for many fungi.

The advent of whole-genome sequence typing (WGST) has made fungal genotyping feasible for outbreak investigations (37), especially for fungi for which conventional genotyping methods do not exist. With WGST, genetic relationships among isolates are determined by the phylogenetic analysis of single-nucleotide polymorphism (SNP) differences among analyzed genomes a population: typically, the fewer the number of SNPs observed between the strains, the more closely the strains are related, and the more likely they are to have a point source, provided supporting epidemiologic evidence exists. Point source outbreaks are typically clonal, and the resulting isolate genomes display few to no SNP differences. Environment-linked outbreaks might have 1 or multiple source populations that also are identifiable and distinct from background and control strains (5).

In January 2014, the Chilean Ministry of Health contacted the Mycotic Diseases Branch (in the Division of Foodborne, Waterborne, and Environmental Disease, National Center for Emerging and Zoonotic Infectious Diseases), US Centers for Disease Control and Prevention (CDC), about a cluster of 67 cases of Sarocladium kiliense (formerly Acremonium kiliense) bloodstream infections (BSI). This cluster was identified at 8 different hospitals in Santiago, Chile, by the National Infection Control Program. The infections occurred during June 2013–January 2014; 39 infections were in children and 2 in adults, all of whom were undergoing chemotherapy. The Chilean Ministry of Health initiated an epidemiologic investigation with technical assistance from CDC and the Pan American Health Organization. An environmental source was considered unlikely because of the spread of these infections among multiple locales in Chile; however, the possibility of environmental contamination could not be excluded.

A detailed review of medication administration records revealed that all patients received 4 intravenous medications: ondansetron, heparin, saline, and potassium. Heparin, saline, and potassium were products used in many different patient populations within the hospital. However, ondansetron, an antinausea medication, was used as standard protocol among oncology patients, and all patients with S. kiliense BSIs received ondansetron from a single source, pharmaceutical company A, in Colombia. Three lots of ondansetron, manufactured by pharmaceutical company A, were investigated in Chile Drug National Agency (ANAMED, ISP). In accordance with Mycotic Diseases Branch recommendations (S. Vallabhaneni, pers. comm.) the Chilean Ministry of Health laboratory cultured 10% of unopened ondansetron vials of these 3 lots. Vials from 2 of the 3 available lots yielded S. kiliense on February 15, 2014. All isolates were identified by traditional methods and confirmed by DNA sequencing by ISP. After this finding, all ondansetron products made by this manufacturer were recalled in Chile; the Pan American Health Organization issued an international health alert on February 17, 2014.

Concurrently, the Corporación para Investigaciones Biológicas (Medellin, Colombia) and the Instituto Nacional de Salud (Bogota, Colombia) contacted CDC about an isolate identified as S. kiliense by conventional DNA sequencing methods and 16 isolates originally identified as a Fusarium spp. by phenotypic methods in Colombia dating to November 2013. Because of the findings in Chile, these isolates were reevaluated and confirmed as S. kiliense by conventional DNA sequencing methods (8). Further investigation by officials in Colombia showed that at least 14 of the 16 patients also received ondansetron manufactured by pharmaceutical company A in Colombia; culturing and conventional DNA sequence identification methods also confirmed that ondansetron was contaminated with S. kiliense.

A subset of isolates from patients and medication vials in this investigation were sent to CDC’s Mycotic Diseases Branch for further identification and molecular typing. To determine whether the contaminated lots of ondansetron harbored the same fungal strains as infected patients, we used WGST.

Methods

Isolates

On the basis of the availability of epidemiologic and patient information, we subjected a subset of patient isolates from each country to molecular analysis: 7 isolates from Chile and 11 isolates from Colombia, 1 isolate per patient. Additionally, 7 isolates from contaminated ondansetron vials were collected from both countries (Table). Eleven unrelated S. kiliense control isolates from the CDC culture collection, American Type Culture Collection, Centraalbureau voor Schimmelcultures, and Universitat Rovira I Virgili also were included for analysis. No background isolates of S. kiliense from the affected countries were available for analysis. Genomic DNA was extracted from cells grown on Sabouraud dextrose agar by using the DNeasy Blood and Tissue kit (QIAGEN, Hilden, Germany) as referenced in Litvintseva et al. (3).

Library Preparation and Illumina Sequencing

The 36 DNA samples were prepared for Illumina sequencing (Illumina, San Diego, CA, USA) by using the KAPA Biosystems Library Preparation with Standard PCR kit (KAPA Biosystems, Wilmington, MA, USA) protocol with a modified 8-bp index. Approximately 1 μg of double-stranded DNA was sheared by using a Sonicman sonicator (Brooks Automation, Spokane, WA, USA) to an average insert size of 650 bp, and DNA libraries were prepared for Illumina paired-end sequencing as described by the manufacturer. All 36 libraries were sequenced to a read length of 100 bp on the Illumina HiSeq 2500 system. Whole-genome sequence read files were deposited in the National Center for Biotechnology Information Sequence Read Archive under BioProject PRJNA291140.

Genome Assembly

We assembled raw sequencing reads from isolate CDC-B10657 using ABySS as a reference genome for S. kiliense (9). The assembly and its corresponding read data were subject to self-alignment to determine coverage. Contigs were filtered for length and coverage; contigs <200 bp were removed from samples assembled from HiSeq data. Contigs with coverage <20% of the average coverage of the 20 largest contigs per assembly were manually removed. We assessed the final assembly for erroneous sites using Pilon (10).

SNP Variant Detection and Phylogenetic Analysis

We applied a reference-based analysis method to determine genetic relatedness among all isolates. Read data of all samples were aligned to a reference by using Novoalign 3.00.03 (Novocraft Technologies, Selangor, Malaysia), and SNPs were identified by using the Genome Analysis Toolkit version 2.4 (11) from a custom pipeline, NASP (Northern Arizona SNP Pipeline, http://tgennorth.github.io/NASP/). The SNP calls were filtered in the final step of the pipeline and were included in the final matrix if they were not identified in repetitive regions, found in <90% of the base calls, and had a minimum read depth coverage of 10×. We excluded reads that mapped to multiple locations within the genome, as well as insertions and deletions. Only SNP loci present in all 36 isolates were included. This final matrix was created from NASP output and converted to FASTA format by using an in-house script NASP. We conducted genomewide SNP-based phylogenetics analyses using the simple parsimony algorithm in PAUP version 4.0b10 (12).

Results

The total length of the de novo assembled reference contigs we showed the genome size of S. kiliense was ≈38 MB. The average sequencing depth for the reference strain was 354×, which resulted in 98% coverage across the genome. The assembly contained 1,092 contigs ranging in length from 500 bp to 1.7 million bp (average contig length 35,528 bp; N50 591,374 bp) (13).

Figure

Thumbnail of Whole-genome single-nucleotide polymorphism (SNP) typing of Sarocladium kiliense strains, Chile and Colombia, 2013–2014. All patient (clinical) and drug (vial) isolates from these 2 countries differed by &lt;5 SNPs, and &gt;21,000 SNPs were identified for the control isolates (≈117,000 total SNPs, ≈73,000 parsimoniously informative SNPs). Scale bar indicates nucleotide substitutions per site.

Figure. Whole-genome single-nucleotide polymorphism (SNP) typing of Sarocladium kiliense strains, Chile and Colombia, 2013–2014. All patient (clinical) and drug (vial) isolates from these 2 countries differed by <5 SNPs, and >21,000...

Reference-based phylogenetic analysis identified ≈117,000 shared SNPs, of which ≈60% were parsimoniously informative (Figure). The cladogram showed 1 distinct clade that comprised all outbreak isolates and included isolates from the patients and the contaminated drugs. No more than 5 SNPs were detected between any patient and drug isolates from Chile and Colombia, demonstrating that these isolates had nearly indistinguishable genomes. Conversely, all control isolates from different sources clustered separately from the outbreak clade.

The genomewide SNP analysis showed greater diversity among the control isolates. Although the closest control isolate, B10764, differed from the outbreak clade by 501 SNPs, the remaining majority of the control isolates diverged from the outbreak clade by ≈21,000 SNPs. Two control isolates, B5504 and B5505, isolated from a 1996 cluster of S. kiliense infections in a US hospital, were genetically indistinguishable from each other.

Discussion

S. kiliense is primarily a saprobic soil organism; it can cause opportunistic infections in humans that typically occur after inoculation of the fungus (14). Here we report the results of molecular epidemiologic investigation of an outbreak of S. kiliense BSI that affected >50 patients in 2 South American countries. The epidemiologic investigation suggested contaminated antinausea medication as a possible source of this infection. With no existing genotyping methods for this uncommon pathogen, we used WGST to understand the genetic relationships among the isolates and identify a potential source of this outbreak.

The use of WGST to investigate fungal outbreaks has become integral to epidemiologic investigations (35,15). During nosocomial outbreak investigations, ascertaining potential source(s) of infection based solely on the descriptive epidemiologic findings often is difficult. Molecular genotyping frequently is needed to test hypotheses generated by epidemiologic investigation: the presence of a single strain or dominant clone usually suggests a point source, whereas the presence of multiple strains is usually consistent with an environmental exposure or exposure to mixed populations from a single source. For example, in a recent outbreak of Curvularia spp. (formerly Bipolaris spp.) among cardiac surgery patients, initial epidemiologic investigation suggested a point source; however, a molecular epidemiology analysis demonstrated multiple strains, consistent with an environmental source (16,17). Similarly, in an outbreak of Fusarium spp. associated with the use of contact lenses, initial epidemiologic investigation linked the infections with a particular lens cleaning solution, suggesting a point source; however, molecular analysis indicated multiple sources (18). In our current investigation, the WGST analysis demonstrated that the patient isolates from Chile and Colombia were nearly genetically indistinguishable (<5 SNPs) from those recovered from the medication vials, indicating the likely presence of a single-source infection. Conversely, the control isolates clustered differently from the outbreak clade by ≈21,000 SNPs, and except for 2 strains from the same cluster, thousands of SNPs separated any 2 control strains (Figure).

Although BSIs with this organism are rare, case reports/clusters of S. kiliense fungal infections have been reported in the literature (1922). In a cluster of S. kiliense infections of endophthalmitis and catheter-related BSIs, an environmental source was strongly suggested but could not be confirmed because of the lack of available typing methods (21). Isolates from this cluster of infections served as controls in our study: specifically, isolates B5504 and B5505 from 2 patients who underwent cataract extraction with intraocular lens implantation differed from each other by 1 SNP, indicating that the 2 genomes were nearly indistinguishable and suggesting a common source of infection.

In this investigation, the use of WGST identified a likely source of S. kiliense BSI in oncology patients in both countries by linking these infections to the receipt of contaminated medication. Results from this outbreak are consistent with those from other outbreaks of fungal pathogens in which WGST was used and a common source was hypothesized. For example, Litvintseva et al. investigated Exerohilium rostratum infections associated with the injection of contaminated methylprednisolone acetate by using whole-genome SNP typing; they found that all outbreak isolates were genotypically indistinguishable: no more than 2 SNPs separated the strains in the outbreak clade (3).

The genetic diversity among the S. kiliense control strains in our study was congruent to the level of diversity among control isolates found in other studies. We identified comparable levels of genetic similarities to those found among isolates of Saprochaete clavata from patients from a multicenter outbreak in France (15). Conversely, >28,000 and up to 1.2 million SNPs were observed among the control isolates of Apophysomyces trapeziformis infections associated with a tornado in Joplin, Missouri, USA (5). Although we were able to include a variety of controls from different sources, a major limitation in this study is the lack of background isolates of S. kiliense from the affected countries to assess genomic similarity among the unrelated isolates from South America.

Traditional molecular strain typing methods, such as MLST, are limited to analysis of specific genomic regions, typically protein-coding regions. However, the conservation in these regions might be insufficient to discriminate strains of certain fungi or provide the resolution needed to identify a likely source in a fungal outbreak. In the absence of established methods for fungi, the rapid development of traditional typing methods is often necessary. In this investigation, before deploying WGST, we evaluated MLST using 7 protein-coding regions to determine the relatedness among all isolates (17). Across 3,651 nt amplified, only a 1-nt polymorphism was identified in the β-tubulin gene of one of the control isolates, which did not provide sufficient resolution between outbreak and control isolates, indicating that the use of MLST was not informative in this outbreak investigation. Although conventional typing methods are rapid, the conclusions drawn from them might be confounded by either their discriminatory power or their character state conflict (i.e., homoplasy). Typically, genotyping methods with high resolution (e.g., variable-number tandem-repeat typing) have higher degrees of homoplasy, whereas methods with lower resolution (e.g., MLST) have less homoplasy; both might be at risk of misidentifying a common source. Conversely, WGST enables accurate deduction of genetic relationships among strains in fungal outbreaks for which current typing methods are ineffective or nonexistent. Furthermore, with the improvement of next-generation sequencing technologies, investigating fungal outbreaks in real time will soon be possible.

Contamination of medical products, particularly with rare fungi, poses growing concern and a public health threat, especially in vulnerable populations. Fungal clusters/outbreaks of common environmental species (rather than classical clinical pathogens) have been associated with medical products. For example, the 2002 outbreak of fungal meningitis was caused by injected steroids contaminated with Exophiala dermatitidis, and the 2012–2015 outbreak of fungal meningitis and other infections resulted from methylprednisolone contaminated by Exserohilum rostratum (3,23). Similarly, in 2009, intestinal zygomycosis resulted from ingestion of allopurinol tablets contaminated with Rhizopus microsporus (24). In 2012, an endophthalmitis outbreak was associated with use of an ophthalmologic dye contaminated with Fusarium sp. in the Fusarium incarnatum-equiseti species complex and triamcinolone contaminated with Curvularia hawaiiensis (Bipolaris hawaiiensis) (16,25,26). Increased vigilance and the use of advanced technologies are needed to rapidly identify the likely source(s) of infection to efficiently guide epidemiologic investigations and initiate appropriate control measures.

In summary, our study highlights the utility of advanced molecular methods to investigate outbreaks involving rare fungi. Next-generation sequencing and bioinformatics analyses will remain critical molecular epidemiology tools in such epidemiologic investigations.

Ms. Etienne is a public health microbiologist at CDC’s Mycotic Diseases Branch, Division of Foodborne, Waterborne, and Environmental Disease, National Center for Emerging and Zoonotic Infectious Diseases, in Atlanta and a Ph.D. student in the Bioinformatics program at Georgia Institute of Technology. Her research interests include integrating advanced laboratory methods, epidemiology and public health strategies to the study of fungal infections.

Top

Acknowledgment

We thank Shawn Lockhart, Joyce Peterson, and Cau Pham for assistance in the microbiological aspects and initial molecular analysis for this outbreak investigation.

Top

References

  1. Repetto  EC, Giacomazzi  CG, Castelli  F. Hospital-related outbreaks due to rare fungal pathogens: a review of the literature from 1990 to June 2011. Eur J Clin Microbiol Infect Dis. 2012;31:2897904. DOIPubMedGoogle Scholar
  2. Guinea  J, Garcia de Viedma  D, Pelaez  T, Escribano  P, Munoz  P, Meis  JF, Molecular epidemiology of Aspergillus fumigatus: an in-depth genotypic analysis of isolates involved in an outbreak of invasive aspergillosis. J Clin Microbiol. 2011;49:3498503. DOIPubMedGoogle Scholar
  3. Litvintseva  AP, Hurst  S, Gade  L, Frace  MA, Hilsabeck  R, Schupp  JM, Whole-genome analysis of Exserohilum rostratum from an outbreak of fungal meningitis and other infections. J Clin Microbiol. 2014;52:321622. DOIPubMedGoogle Scholar
  4. Gillece  JD, Schupp  JM, Balajee  SA, Harris  J, Pearson  T, Yan  Y, Whole genome sequence analysis of Cryptococcus gattii from the Pacific Northwest reveals unexpected diversity. PLoS ONE. 2011;6:e28550. DOIPubMedGoogle Scholar
  5. Etienne  KA, Gillece  J, Hilsabeck  R, Schupp  JM, Colman  R, Lockhart  SR, Whole genome sequence typing to investigate the Apophysomyces outbreak following a tornado in Joplin, Missouri, 2011. PLoS ONE. 2012;7:e49989. DOIPubMedGoogle Scholar
  6. Engelthaler  DM, Hicks  ND, Gillece  JD, Roe  CC, Schupp  JM, Driebe  EM, Cryptococcus gattii in North American Pacific Northwest: whole-population genome analysis provides insights into species evolution and dispersal. MBio. 2014;5:e0146414. DOIPubMedGoogle Scholar
  7. Litvintseva  AP, Brandt  ME, Mody  RK, Lockhart  SR. Investigating fungal outbreaks in the 21st century. PLoS Pathog. 2015;11:e1004804. DOIPubMedGoogle Scholar
  8. Petti  CA, Bosshard  PP, Brandt  ME, Clarridge  JE, Feldblyum  TV, Foxall  P, Interpretive criteria for identification of bacteria and fungi by DNA target sequencing: approved guideline. MM18-A. Wayne (PA): Clinical and Laboratory Standards Institute; 2008.
  9. Simpson  JT, Wong  K, Jackman  SD, Schein  JE, Jones  SJ, Birol  I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:111723. DOIPubMedGoogle Scholar
  10. Walker  BJ, Abeel  T, Shea  T, Priest  M, Abouelliel  A, Sakthikumar  S, Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE. 2014;9:e112963. DOIPubMedGoogle Scholar
  11. McKenna  A, Hanna  M, Banks  E, Sivachenko  A, Cibulskis  K, Kernytsky  A, The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297303 and. DOIPubMedGoogle Scholar
  12. Wilgenbusch  JC, Swofford  D. Inferring evolutionary trees with PAUP*. Curr Protoc Bioinformatics. 2003;Chapter 6:Unit 6.4
  13. Miller  JR, Koren  S, Sutton  G. Assembly algorithms for next-generation sequencing data. Genomics. 2010;95:31527. DOIPubMedGoogle Scholar
  14. Summerbell  RC, Gueidan  C, Schroers  HJ, de Hoog  GS, Starink  M, Rosete  YA, Acremonium phylogenetic overview and revision of Gliomastix, Sarocladium, and Trichothecium. Stud Mycol. 2011;68:13962. DOIPubMedGoogle Scholar
  15. Vaux  S, Criscuolo  A, Desnos-Ollivier  M, Diancourt  L, Tarnaud  C, Vandenbogaert  M, Multicenter outbreak of infections by Saprochaete clavata, an unrecognized opportunistic fungal pathogen. MBio. 2014;5:e0230914. DOIPubMedGoogle Scholar
  16. Manamgoda  DS, Rossman  AY, Castlebury  LA, Crous  PW, Madrid  H, Chukeatirote  E, The genus Bipolaris. Stud Mycol. 2014;79:22188. DOIPubMedGoogle Scholar
  17. Pham  CD, Purfield  AE, Fader  R, Pascoe  N, Lockhart  SR. Development of a multilocus sequence typing system for medically relevant Bipolaris species. J Clin Microbiol. 2015;53:323946. DOIPubMedGoogle Scholar
  18. Chang  DC, Grant  GB, O’Donnell  K, Wannemuehler  KA, Noble-Wang  J, Rao  CY, Multistate outbreak of Fusarium keratitis associated with use of a contact lens solution. JAMA. 2006;296:95363.PubMedGoogle Scholar
  19. Pastorino  AC, Menezes  UP, Marques  HH, Vallada  MG, Cappellozi  VL, Carnide  EM, Acremonium kiliense infection in a child with chronic granulomatous disease. Braz J Infect Dis. 2005;9:52934. DOIPubMedGoogle Scholar
  20. Ioakimidou  A, Vyzantiadis  TA, Sakellari  I, Arabatzis  M, Smias  C, Douka  V, An unusual cluster of Acremonium kiliense fungaemias in a haematopoietic cell transplantation unit. Diagn Microbiol Infect Dis. 2013;75:3136. DOIPubMedGoogle Scholar
  21. Fridkin  SK, Kremer  FB, Bland  LA, Padhye  A, McNeil  MM, Jarvis  WR. Acremonium kiliense endophthalmitis that occurred after cataract extraction in an ambulatory surgical center and was traced to an environmental reservoir. Clin Infect Dis. 1996;22:2227. DOIPubMedGoogle Scholar
  22. Das  S, Saha  R, Dar  SA, Ramachandran  VG. Acremonium species: a review of the etiological agents of emerging hyalohyphomycosis. Mycopathologia. 2010;170:36175. DOIPubMedGoogle Scholar
  23. Centers for Disease Control and Prevention. Exophiala infection from contaminated injectable steroids prepared by a compounding pharmacy—United States, July–November 2002. MMWR Morb Mortal Wkly Rep. 2002;51:110912 .PubMedGoogle Scholar
  24. Cheng  VC, Chan  JF, Ngan  AH, To  KK, Leung  SY, Tsoi  HW, Outbreak of intestinal infection due to Rhizopus microsporus. J Clin Microbiol. 2009;47:283443. DOIPubMedGoogle Scholar
  25. Mikosz  CA, Smith  RM, Kim  M, Tyson  C, Lee  EH, Adams  E, Fungal endophthalmitis associated with compounded products. Emerg Infect Dis. 2014;20:24856. DOIPubMedGoogle Scholar
  26. Centers for Disease Control and Prevention. Multistate outbreak of postprocedural fungal endophthalmitis associated with a single compounding pharmacy—United States, March–April 2012. MMWR Morb Mortal Wkly Rep. 2012;61:3101 .PubMedGoogle Scholar

Top

Figure
Table

Top

Cite This Article

DOI: 10.3201/eid2203.151193

Table of Contents – Volume 22, Number 3—March 2016

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.

Top

Comments

Please use the form below to submit correspondence to the authors or contact them at the following address:

Kizee A. Etienne, Centers for Disease Control and Prevention, 1600 Clifton Rd NE, Mailstop G11, Atlanta, GA 30329-4027, USA

Send To

10000 character(s) remaining.

Top

Page created: February 17, 2016
Page updated: February 17, 2016
Page reviewed: February 17, 2016
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external