Tracing mumps transmission in Washington state
Louise H. Moncla*1, Allison Black*1,2, Misty Lang3, Nicholas R. Graff3, Ailyn C. Perez-Osorio3, Nicola F. Müller1, Dirk Haselow4, Chas DeBolt3,Trevor Bedford1,2
1Fred Hutchinson Cancer Research Center, Seattle, Washington, United States 1University of Washington, Seattle, Washington, United States 3Office of Communicable Disease Epidemiology, Washington State Department of Health, United States 4Arkansas Department of Health, Little Rock, Arkansas, United States
*these authors contributed equally to this work
Mumps has resurged among vaccinated individuals in recent years, sparking an outbreak in 2016 and 2017 that affected nearly 12,000. Burgeoning evidence demonstrates that mumps vaccine-induced immunity wanes over time, triggering post-vaccine era outbreaks to which young adults are particularly susceptible. During the 2016/2017 outbreak, Washington state experienced one of the highest incidence rates of mumps infection, despite high overall vaccination rates. Curiously, Washington cases were enriched among school-aged children and members of a particular close-knit ethnic community, suggesting that waning immunity alone may not explain transmission. To elucidate how the mumps outbreak in Washington originated and spread, we sequenced 110 near-complete genomes of mumps viruses collected throughout the outbreak, and trace mumps transmission into and within the state. Contrary to contact tracing data, we uncover that mumps was introduced into Washington at least 13 times, primarily from Arkansas, sparking overlapping transmission chains. We do not find evidence for antigenic evolution, or that vaccination status or age were the primary determinants of transmission. Instead, the outbreak was overwhelmingly sustained by transmission within a close-knit ethnic community. Community members were twice as likely to transmit mumps as non-community members, and transmission to non-community members lead to short, terminal transmission chains. These data suggest that transmission was inefficient outside of close-contact settings. Our data underscore the ability of genomic data to clarify epidemiologic factors driving transmission, and pinpoint contact networks as critical determinants of mumps transmission in Washington. Our results suggest that in contact patterns should be considered along with waning immunity when formulating outbreak response and 3rd dose vaccine recommendations.
This repo contains all of the code used to perform the analyses for this project tracing mumps transmission during an outbreak in Washington. Sequence data generated for this project can be found here, along with the protocols and pipelines for amplifying and sequencing mumps genomes from buccal swabs. We have deposited all of the raw sequencing reads (with human reads removed) onto the Short Read Archive under Project number PRJNA641715. Consensus genomes are currently being deposited into Genbank under accession numbers MT859507-MT859672, although they have not yet been approved.
Alignments and xml files used for beast analyses are available in the
structured-coalescent-analyses directories. We currently are not able to publish full metadata that links the sequence strain name to age, vaccination status, and Marshallese status. We have therefore removed these metadata fields from the uploaded alignments and xml files and replaced them with generic placeholder text, like
vaccination_status. If we receive permission from the Department of Health to release this metadata, we will update these files.
Similarly, we have a public metadata file with the date of sample collection, Ct value, specimen type, and vaccination status available here. If we are able to release further metadata, we will update this file with the appropriate fields.
Many of the plotting scripts in this repo require baltic. I have also included a version of baltic that will work with all of these jupyter notebooks within this repository.