Genetic characterization of mumps viruses associated with the resurgence of mumps in the United States: 2015-2017

Despite high coverage with measles, mumps, and rubella vaccine in the United States, outbreaks of mumps occur in close contact settings such as schools, colleges, and camps. Starting in late 2015, outbreaks were reported from several universities, and by the end of 2017, greater than 13,800 cases had been reported nation-wide. In 2013, the CDC and the Association of Public Health Laboratories contracted four Vaccine Preventable Diseases Reference Centers (VPD-RCs) to perform real-time reverse transcription PCR (RT-qPCR) to detect mumps RNA in clinical samples and to determine the genotype. Twelve genotypes of mumps virus are currently recognized by the World Health Organization, and the standard protocol for genotyping requires sequencing the entire gene coding for the small hydrophobic (SH) protein. Phylogenetic analysis of the 1862 mumps samples genotyped from 2015 through 2017 showed that the overall diversity of genotypes detected was low. Only 0.8 % of the sequences were identified as genotypes C, H, J, or K, and 0.5 % were identified as vaccine strains in genotypes A or N, while most sequences (98.7 %) were genotype G. The majority of the genotype G sequences could be included into one of two large groups with identical SH sequences. Within genotype G, a small number of phylogenetically significant outlier sequences were associated with epidemiologically distinct chains of transmission. These results demonstrate that molecular and epidemiologic data can be used to track transmission pathways of mumps virus; however, the limited diversity of the SH sequences may be insufficient for resolving transmission in all outbreaks.

Variability analysis and inter-genotype comparison of human respiratory syncytial virus small hydrophobic gene

Background

Small hydrophobic (SH) gene is one of the mostly diverse genomic regions of human respiratory syncytial virus (HRSV). Its coding region constitutes less than 50% of the complete gene length, enabling SH gene to be highly variable and the SH protein highly conserved. In standard HRSV molecular epidemiology studies, solely sequences of the second hypervariable region of the glycoprotein gene (HVR2) are analyzed. To what extent do the strains identical in HVR2 differ elsewhere in genomes is rarely investigated. Our goal was to investigate whether diversity and inter-genotype differences observed for HVR2 are also present in the SH gene.

Methods

We sequenced 198 clinical samples collected within a limited area and time frame. In this HRSV collection, rapid and significant changes in HVR2 occurred.

Results

Over 20% of strains from this pool (containing HRSV genotypes NA1, ON1, GA5, BA9 and BA10) would be incorrectly assumed to be identical to another strain if only the HVR2 region was analysed. The majority of differences found in SH gene were located in the 5′ untranslated region (UTR). Seven indels were detected, one was genotype GA5 specific. An in-frame deletion of 9 nucleotides (coding for amino acids 49–51) was observed in one of group A strains. Fifteen different SH protein sequences were detected; 68% of strains possessed the consensus sequence and most of others differed from the consensus in only one amino acid (only 4 strains differed in 2 amino acids). The majority of differing amino acids in group A viruses had the same identity as the corresponding amino acids in group B strains. When analysis was restricted to strains with identical HVR2 nucleotide sequences and differing SH protein sequences, 75% of differences observed in the SH ectodomain were located within region coding for amino acids 49–51.

Conclusions

Basing HRSV molecular epidemiology studies solely on HVR2 largely underestimates the complexity of circulating virus populations. In strain identification, broadening of the genomic target sequence to SH gene would provide a more comprehensive insight into viral pool versatility and its evolutionary processes.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-1020-9) contains supplementary material, which is available to authorized users.

Mumps virus F gene and HN gene sequencing as a molecular tool to study mumps virus transmission

Various mumps outbreaks have occurred in the Netherlands since 2004, particularly among persons who had received 2 doses of measles, mumps, and rubella (MMR) vaccination. Genomic typing of pathogens can be used to track outbreaks, but the established genotyping of mumps virus based on the small hydrophobic (SH) gene sequences did not provide sufficient resolution. Therefore, we expanded the sequencing to include fusion (F) gene and haemagglutinin-neuraminidase (HN) gene sequences in addition to the SH gene sequences from 109 mumps virus genotype G strains obtained between 2004 and mid 2015 in the Netherlands. When the molecular information from these 3 genes was combined, we were able to identify separate mumps virus clusters and track mumps virus transmission. The analyses suggested that multiple mumps virus introductions occurred in the Netherlands between 2004 and 2015 resulting in several mumps outbreaks throughout this period, whereas during some local outbreaks the molecular data pointed towards endemic circulation. Combined analysis of epidemiological data and sequence data collected in 2015 showed good support for the phylogenetic clustering.