Genetic heterogeneity of L-Zagreb mumps virus vaccine strain
Virology Journal volume 5, Article number: 79 (2008)
The most often used mumps vaccine strains Jeryl Lynn (JL), RIT4385, Urabe-AM9, L-Zagreb and L-3 differ in immunogenicity and reactogenicity. Previous analyses showed that JL, Urabe-AM9 and L-3 are genetically heterogeneous.
We identified the heterogeneity of L-Zagreb throughout the entire genome. Two major variants were defined: variant A being identical to the consensus sequence of viral seeds and vaccine(s) and variant B which differs from variant A in three nucleotide positions. The difference between viral variants in L-Zagreb strain is insufficient for distinct viral strains to be defined. We demonstrated that proportion of variants in L-Zagreb viral population depends on cell substrate used for viral replication in vitro and in vivo.
L-Zagreb strain should be considered as a single strain composed of at least two variant viral genomes.
Mumps virus (MuV) genome consists of a 15 384 nt long non-segmented single-stranded negative sense RNA. The genomic RNA contains seven genes which encode nine open reading frames: NP (nucleoprotein), P (phosphoprotein, V protein, I protein), M (matrix protein), F (fusion protein), SH (small hydrophobic protein), HN (haemagglutinin-neuraminidase) and L (large protein) [1, 2].
Mumps virus causes an acute systemic infection involving glandular, lymphoid and nervous tissues. Prior to the introduction of live attenuated virus vaccines, mumps virus was a leading cause of the virus-induced CNS disease .
Live attenuated mumps vaccines have been used worldwide since late 1960s [3, 4]. Nowadays, the most often used vaccine strains are Jeryl Lynn (JL), RIT 4385, Urabe-AM9, L-Zagreb and Leningrad-3 (L-3) [5, 6].
Although at the time of their development the knowledge of molecular content of mumps vaccines was not the issue, recently it has become obvious that the molecular consistency of vaccine production is not a trivial matter. Sauder et al  showed that the change in genetic heterogeneity at the specific genome sites can have a profound effect on neurovirulent phenotype of Urabe-AM9 strain. The RNA viral population consists of virus particles that differ from the consensus sequence in one or more nucleotides (quasisipecies), the feature that arises because of the high mutation rate of RNA-dependent RNA polymerase (RdRp) (10-3 to 10-5 errors per nucleotide site and replication cycle) [8, 9]. Given that all mumps vaccines are quasiespecies populations, an adequate description of the vaccine virus genome should include not only the consensus sequence, but also the quantitative assessment of the existing viral variants.
Previous analyses confirmed that mumps vaccine strains JL, Urabe-AM9 and L-3 are genetically heterogeneous. JL is composed of a mixture of two distinct viral strains (JL5 and JL2) [10, 11] while Urabe-AM9 represents a quasispecies mix [12, 13].
L-3 vaccine strain was prepared from five mumps virus isolates combined into a single strain in 1953 . It was characterized as heterogenic on the basis of plaque morphology  and a sequence autoradiogram with several ambiguities in P and F genes  but precise vaccine composition of L-3 was never published.
L-Zagreb vaccine strain was developed by further subcultivation of L-3 mumps vaccine strain in primary culture of chicken embryo fibroblast (CEF) . Genetic stability at the level of the consensus sequence of the L-Zagreb vaccine strain in the course of the production process was demonstrated .
Here, we analyzed the detailed genetic composition of L-Zagreb vaccine strain. Due to mixture of mumps virus isolates in L-3 production and its heterogeneity we wonder about the composition of L-Zagreb strain.
By two independent cloning experiments we showed that L-Zagreb vaccine strain contains only one viral strain. However, numerous nucleotide positions showed to be heterogenic and indicating a quasispecies nature of this strain. We successfully isolated two types of viral clones: identical to consensus sequence, named as variant A, and with the nucleotide sequences different from the consensus sequence (quasispecies). The most abundant quasispecies, named variant B, was detected in all analyzed L-Zagreb samples.
Finally, we demonstrated that the heterogenic composition of L-Zagreb strain strongly depends on the number of passages and the type of the cell culture that the virus is replicating on.
Results and discussion
Heterogenic nucleotide positions in the L-Zagreb vaccine strain genome
The strategy for defining heterogenic positions in the L-Zagreb vaccine strain genome involved cloning of eleven overlapping PCR fragments into pUC19 plasmid vector and sequencing of resulting plasmid clones. For each fragment, two independent cloning experiments were performed in order to avoid misinterpretation of artificial heterogeneity arisen from the error of Pfu DNA polymerase used for fragment amplification . Twenty and ten clones were analyzed in the first and the second experiment, respectively. Cloned genome fragments were compared to the consensus sequence of the L-Zagreb strain [GenBank: AY685920] in order to select clones with changed nucleotides.
As a result, 88 and 49 nucleotides different from the consensus sequence were identified throughout a complete genome of L-Zagreb strain, in the first and the second experiment respectively (Fig 1). The distribution of the heterogenic positions seem to be at random except for the region between approx. 2000 and 3000 nt which corresponds to almost a complete coding region of P gene (which spans region 1979–3152 nt) where no heterogeneity was found (Fig 1).
By comparing changed nucleotides of the first and the second cloning experiment, alterations of six nucleotides (1059, 1073, 1996, 5261, 11345 and 13054) were identified in both cloning experiments in one or more clones. Alterations G1059A was identified in 1/20 and 1/10 sequenced plasmid clones, G1073T in 5/20 and 1/10; A1996C in 1/20 and 1/10, C5261T in 7/20 and 4/10, G11345T in 5/20 and 3/10 plasmid clones and alteration C13054A was identified in 2/20 and 1/10 plasmid clones (Fig 1).
Other 125 changed nucleotides were found within a single clone in one of the experiments (Fig 1). They could easily be considered as possible heterogeneities in L-Zagreb strain, although it should not be ruled out that some of them originated from a low error rate of the enzymes involved in fragment amplification.
Heterogeneity of isolated viral clones
Based on the above results it may be predicted that there are different viral variants constituting L-Zagreb strain. However, solely by defining heterogenic positions in cloned fragments it was not possible to identify the variants. Therefore, viral clones of the L-Zagreb strain were plaque isolated in the Vero cell culture and additionally subcultivated once in the same cell culture. Twelve viral clones were identified by sequencing in regions which included six heterogenic positions identified in cloning experiments: 1059, 1073, 1996, 5261, 11345 and 13054. The analysis of the picked viral clones indicated two viral variants, herein named variant A and variant B. Variant A was shown to be identical to the consensus sequence in all 6 positions: 1059G, 1073G, 1996A, 5261C, 11345G, and 13054C. Variant B differed from the consensus sequence in positions 1073T, 5261T and 11345T.
Seven out of 12 analyzed viral clones represented variant A and five out of 12 represented variant B.
In the process of determination of the viral variants we used only 6 positions identified as heterogenic by the cloning experiments. As that could limit the criteria for defining the variants, a complete genome of one viral clone belonging to variant B was sequenced and compared to the consensus sequence. Again, the differences were only in positions 1073, 5261 and 11345, what confirmed our criteria for defining variant B as adequate. Also it confirmed the presence of Ts at positions 1073, 5261 and 11345 as the genetic marker of variant B versus variant A.
Three positions (1073, 5261 and 11345) which differentiate variants A and B caused diversity in amino acid composition of the N (310 aa, A→S), F (239 aa, T→I) and L (970 aa, V→L) protein. Although these amino acids are not part of any known functional domains it is difficult to predict biological impact of these differences since the molecular structures of mumps virus proteins are not well characterized.
Heterogenic positions 1059, 1996 and 13054 were not detected in any of the twelve analyzed viral clones. Since they were identified in both cloning experiments, they can not be considered as misinterpretation due to the erroneous nucleotide incorporation in the amplification process. They should rather be considered as differential positions of viral variants which coexist with variant A and B, but at low amount. Although heterogenic positions represented by low amount seem at random and not relevant for the virus, the in vivo impact of these minor variants should not be minimized. Previous study  indicated the existence of multiple genetic markers and a need for evaluation of a total viral population instead of putting the relevance of a consensus sequence in foreground.
Since L-Zagreb vaccine strain originated from L-3 vaccine strain , we compared the nucleotide sequences of both variants to the consensus sequence of L-3 [GenBank: AY508995] (Table 1). Comparison of the entire genome of L-3 strain with the variant A and variant B genomes showed difference in five and eight nucleotides, respectively (Table 1). Since the heterogeneity of L-3 strain is not completely resolved it is not possible to define when the heterogenic positions in L-Zagreb strain occurred: were they newly created by mutational process in the course of adaptation in primary CEF culture or was it the selection of preexisting minor quasispecies of L-3 strain which occurred when Japanese quail embryo fibroblast culture was replaced by primary CEF.
High level of difference between JL5 and JL2 (414 nt, 87 aa), demonstrated JL vaccine as a mixture of two distinct viral strain . In contrast to JL, the two most abundant viral variants in L-Zagreb strain differ in three nucleotide positions (three aa) what is an insufficient variation for two distinct viral strains to be defined. L-Zagreb strain should be considered as a single strain composed of a number of quasispecies viral genomes. Thus, L-Zagreb strain is similar to Urabe-AM9 vaccine strain which consists of at least two viral variants with minor genetic changes [19, 12, 13].
Screening of heterogeneity in the L-Zagreb seeds and vaccine lots
Here, the same sample was used for both cloning experiments and plaque isolation. Due to the high genetic plasticity of RNA viruses one could easily assume that the identified heterogenic positions reflect only the composition of that sample and is not the intrinsic feature of L-Zagreb vaccine strain.
Therefore the heterogeneity of L-Zagreb strain was analyzed in vaccine seeds (master and working) and two final vaccine batches. The T in heterogenic position 5261 in the F gene is located within the SspI restriction site while the C in the same position eliminates this restriction site. The existence or the absence of restriction site facilitated the use of PCR-RFLP assay as an adequate method for distinguishing the two variants, A and B. A 321 bp uncleaved fragment indicated variant A while a 219 bp cleaved fragment indicated variant B.
Heterogenic positions 1073 and 11345 which are also genetic markers of variant B were unsuitable for RFLP assay. However, it was proved above that these three positions represent the marker of variant B.
Both variants were detected in all four viral samples, but at different proportions (Table 2, Fig 2A). The master seed consisted of 93.0 ± 1.8% variant A and 7.0 ± 1.8% variant B. In the working seed the proportion of variant B increased to 9.7 ± 1.1% and the proportion of variant A decreased to 90.3 ± 1.1%. The two final vaccine batches contained even lower proportions of variant A (80.4 ± 2.6% and 79.4 ± 2.9%, respectively) and higher proportions of variant B (19.6 ± 2.6% and 20.6 ± 2.9%, respectively) in comparison to the viral seeds.
Altogether, these data indicate that different L-Zagreb vaccine strain samples are heterogenic in the same nucleotide positions. Although the stability of the L-Zagreb consensus sequence of the master seed and a final vaccine batch was confirmed , it is clear that the quasispecies content is changing by the passage number.
Influence of different cell culture on heterogeneity of L-Zagreb
Different production stages of L-Zagreb vaccine analyzed above originated through increasing number of passages in primary chicken embryo fibroblasts (CEF). Analysis of the proportion of variants A and B in the vaccine samples clearly shows that CEF somewhat favor replication of variant B over variant A (Fig 2A). The influence of selected cell line on replication efficiency of viral clones and thus a genetic heterogeneity of the entire viral sample was also reported previously .
A limited number of serial passages of L-Zagreb vaccine strain in Vero and SH-SY5Y cell line further confirmed the influence of the cell culture selection on the genetic composition of L-Zagreb strain (Fig 2B and 2C). Cell supernatants from each of five passages were analyzed by PCR-RFLP assay in position 5261. The original sample (p0) contained 77.7 ± 1.5% of variant A and 22.3 ± 1.5% of variant B. The first passage decreased the proportion of variant B in both Vero and SH-SY5Y cells to 12.2 ± 1.5% and 10.0 ± 3.6%, respectively (Fig 2B and 2C, respectively). Surprisingly, the proportion of variant B was diminished to an undetectable level already in the following passage and remained undetectable for the next three passages in both cell types (Fig 2B and 2C). This clearly shows that both Vero and SH-SY5Y cells, in contrast to CEF, promote the replication of variant A leading to the loss of variant B, the second most abundant genomic variant in L-Zagreb strain.
However, when propagated as a genetically homogenous viral clone, variant B was able to replicate in both cell types up to the same extent as variant A (data not shown) indicating that variant A either uses up the cell capacity to produce viral particles faster then variant B or it actively suppresses the replication of variant B by an unknown mechanism.
Our data confirm heterogeneous nature of the L-Zagreb mumps vaccine strain and permanent existence of minor variant B whose replication is favored by the production cell culture. Variant B differed from consensus sequence in only tree nucleotides what give us reason to conclude that vaccine strain L-Zagreb is composed of only one mumps viral strain. Also we detected low represented variant changed in nt positions 1059, 1996 and/or 13054.
Sauder et al.  showed that changes in the neurovirulent phenotype were merely associated with the changes of the level of genetic heterogeneity. In addition to cell substrate, serial passaging of mumps virus could reduce or increase neurovirulent phenotype of the virus.
Due to its heterogenic profile, examined L-Zagreb vaccine lot could serve as comparator in investigations of genetic profile of L-Zagreb postvaccinal mumps virus  or L-Zagreb horizontally transmitted mumps virus .
Materials and methods
Viral material and cell culture
L-Zagreb master and working seeds were stored at -80°C until used for RNA extraction. Freeze-dried L-Zagreb vaccine lots 1 and 2 were reconstituted in 250 μl of water prior to RNA extraction or in 500 μl of sterile water prior to plaque assay. All viral materials were produced at the Institute of Immunology Inc, Zagreb, Croatia.
Vero cell culture (African green monkey kidney cells) was obtained from the American Type Culture Collection (USA) and cultured in minimal essential medium (MEM-H) (AppliChem, Germany) supplemented with 10% fetal calf serum (FCS) (Moregate, Australia) and neomycin 50 μl/ml (Gibco-BRL, USA). Human neuroblastoma cell culture, SH-SY5Y was obtained from European Collection of Cell Cultures, (UK) and cultured in D-MEM medium supplemented with 10% FCS and neomycin 50 μl/ml.
Viral RNA was extracted from 250 μl of viral seeds or reconstituted vaccine as reported by Chomczynski and Mackey .
cDNA synthesis and amplification of the complete genome segmented in eleven overlapping fragments, were performed as described in Ivancic et al . Briefly, RNA was reverse transcribed with random hexamers and MuLV (Applied Biosystems, USA) for 1 h at 42°C. Amplification was performed with the whole reverse transcription mix containing 2.4 U Pfu DNA polymerase (Promega, USA) in a total volume of 100 μl.
Successful amplification of fragments of the expected size was confirmed by electrophoresis in 1% agarose gels. The fragments were purified using a Qiaquick kit (Qiagen, Germany) and cloned into pUC19 vector (New England Biolabs, USA) using protocol for blunt-end cloning described in pMOSBlue Blunt Ended Cloning Kit (GE Healthcare, UK) with modification.
Briefly, plasmid DNA was linearised with SmaI and dephosphorylated using alkaline phosphatase (Roche, Germany). DNA fragments were converted into blunt, phosphorylated products in a one step reaction using PK enzyme mix from pMOSBlue Blunt Ended Cloning Kit (GE Healthcare, UK). Followed by brief heat incubation the product was ligated overnight at 16°C into dephosphorylated blunt ended pUC19.
Escherichia coli strain DH5amcrAB (Life Technologies Ltd, USA) was transformed with the ligation mix. Transformed bacteria were picked out using blue-white selection. Plasmid DNA was isolated from overnight cultures by alkaline lysis .
Segments of mumps virus genome were sequenced either upon cloning into pUC19 vector or directly as PCR products of isolated viral clones by using M13 FSP (5'gctggcgaaagggggatgtg3') and M13 RSP (5'cactttatgcttccggctcg3') primers or mumps virus specific primers . Sequencing was performed on a 3130 Genetic Analyzer (Applied Biosystems, USA) using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, USA) according to the protocol recommended by the manufacturer. Obtained sequences were analyzed by CloneManager Suite software (Scientific&Educational Software, USA).
Isolation of viral clones
The 8 × 105 Vero cells were grown in six-well plates for 24 h. Virus was diluted in MEM-H with 2% FCS and neomycin. Aspirated cell monolayer was infected with 0.5 ml of viral suspension. After 1 h at 37°C viral suspension was aspirated, cells were washed twice with PBS and 3 ml of overlay I (1 v/v 2 × MEM-H with 10% FBS without phenol red and 1 v/v 1.4% Noble agar (Sigma, USA)). Plates were incubated at 37°C in a humidified atmosphere of 5% CO2. After five days 1 ml of overlay II (0.02% Neutral red (Sigma, USA) plus overlay I) was added. Plates with solidified overlay II were incubated at 37°C for the next 24 h in a humidified atmosphere of 5% CO2. Viral clones were cut out and transferred on Vero cells for one additional passage in order to prepare viral suspension for the viral identification by DNA sequencing.
Virus passages in cell cultures
Vero and SH-SY5Y were grown in six-well plates for 24 h. Cells were infected with monovalent L-Zagreb vaccine at m.o.i. 0.05 and incubated at 37°C for 1 h. Cells were then washed twice with PBS. Infected Vero and SH-SY5Y cells were incubated in MEM-H or D-MEM, respectively, with 2% FCS and neomycin at 37°C in a humidified atmosphere of 5% CO2. After 48 h the medium was collected and used for further passage and viral identification by DNA sequencing. Each consecutive passage was done with 1/10 of collected culture medium. Five consecutive passages in both cell lines were performed.
PCR products were purified by using QIAquick kit (Qiagen, Germany) and cleaved in a reaction mixtures consisting of 0.05 μg of PCR product, 10 U of SspI (GE Healthcare, UK) and 1× cleavage buffer in a total volume of 25 μl. The reaction was carried out overnight at 37°C. Cleaved PCR product was diluted 1:10 with water, and 1 μl was mixed with 0.5 μl of GS-500 LIZ size standard (Applied Biosystems, USA) and 9.5 μl of HiDi formamide (Applied Biosystems, USA). The mixture was denaturated at 95°C for 2 min followed by cooling on ice. Electrophoresis of cleaved PCR products was performed on 3130 Genetic Analyzer (Applied Biosystems, USA) using POP7 polymer (Applied Biosystems, USA). Analyses of PCR products were done by GeneMapper (Applied Biosystems, USA) software using peak area data .
Carbone KM, Wolinsky JS: Mumps Virus. Fields-Virology. Edited by: Fields BN, Howley PM, Griffin DE, Lamb RA, Martin MA, Roizman B, Straus SE, Knipe DM. 2001, 1381-1400. 4
Elango N, Varsanyi TM, Kovamees J, Norrby E: Molecular cloning and characterization of six genes, determination oF-gene order and intergenic sequences and leader sequence of mumps virus. J Gen Virol. 1988, 69: 2893-2900. 10.1099/0022-1317-69-11-2893.
Hilleman MR: The development of live attenuated mumps virus vaccine in historic perspective and its role in the evolution of combined measles-mumps-rubella. Vaccinia, vaccination and vaccinology: Jenner, Pasteur and their successors. Edited by: Plotkin S, Fantini B. 1996, 283-292.
Smorodintsev AA: New live vaccines against virus diseases. AJPH. 1960, 50: 40-45.
WHO: Mumps virus vaccines. WER. 2001, 45: 346-356.
WHO: Mumps virus vaccines. WER. 2007, 82: 51-60.
Sauder CJ, Vandenburgh KM, Iskow RC, Malik T, Carbone KM, Rubin SA: Changes in mumps virus neurovirulence phenotype associated with quasispecies heterogeneity. Virology. 2006, 350: 48-57. 10.1016/j.virol.2006.01.035.
Domingo E, Holland JJ: RNA virus mutations and fitness for survival. Annu Rev Microbiol. 1997, 51: 151-178. 10.1146/annurev.micro.51.1.151.
Holland JJ, de la Torre HC, Clarke D, Duarte E: Quantitation of Relative Fitness and Great Adaptability of Clonal Populations of RNA Viruses. J Virol. 1991, 65: 2960-2967.
Afzal MA, Pickford AR, Forsey T, Heath AB, Minor PD: The Jeryl Lynn vaccine strain of mumps virus is a mixture of two distinct isolates. J Gen Virol. 1993, 74: 917-920. 10.1099/0022-1317-74-5-917.
Amexis G, Rubin S, Chizhikov V, Pelloquin F, Carbone K, Chumakov K: Sequence diversity of Jeryl Lynn strain of mumps virus: quantitative mutant analysis for vaccine quality control. Virology. 2002, 300: 171-179. 10.1006/viro.2002.1499.
Amexis G, Fineschi N, Chumakov K: Correlation of genetic variability with safety of mumps vaccine Urabe AM9 strain. Virology. 2001, 287: 234-241. 10.1006/viro.2001.1009.
Wright KE, Dimock K, Brown EG: Biological characteristic of genetic variants of Urabe AM9 mumps vaccine virus. Virus Res. 2000, 67: 49-57. 10.1016/S0168-1702(00)00129-5.
Boriskin YS, Kaptsova TI, Booth JC: Mumps virus variants in heterogeneous mumps vaccine. Lancet. 1993, 341: 318-319. 10.1016/0140-6736(93)92683-K.
Boriskin YS, Yamada A, Kaptsova TI, Skvortsova OI, Sinitsyna OA, Takeuchi K, Tanabayashi K, Sugiura A: Genetic evidence for variant selection in the course of dilute passaging of mumps vaccine virus. Res Virol. 1992, 143: 279-283. 10.1016/S0923-2516(06)80116-0.
Beck M, Weisz-Maleèek R, Meško-Prejac M, Radman V, Juzbašić M, Rajninger-Miholić M, Prislin-Musklic M, Dobrovsak-Sourek V, Smerdel S, Stainer DW: Mumps vaccine L-Zagreb, prepared in chick fibroblasts. I. Production and field trials. J Biol Stand. 1989, 17: 85-90. 10.1016/0092-1157(89)90031-0.
Ivancic J, Gulija TK, Forcic D, Baricevic M, Jug R, Mesko-Prejac M, Mazuran R: Genetic characterization of L-Zagreb mumps vaccine strain. Virus Res. 2005, 109: 95-105. 10.1016/j.virusres.2004.11.013.
Malet I, Belnard M, Agut H, Cahour A: From RNA to quasispecies: a DNA polymerase with proofreading activity is highly recommended for accurate assessment of viral diversity. J Virol Methods. 2003, 109: 161-170. 10.1016/S0166-0934(03)00067-3.
Brown EG, Dimock K, Wright KE: The Urabe AM9 mumps vaccine is a mixture of viruses differing at amino acid 335 of the hemaglutinin-neuraminidase gene with one form associated with disease. J Infect Dis. 1996, 174: 619-622.
Santos-Lopez G, Cruz C, Pazos N, Vallejo V, Reyes-Leyva J, Tapia-Ramirez J: Two clones obtained from Urabe AM9 mumps virus vaccine differ in their replicative efficiency in neuroblastoma cells. Microbes Infect. 2006, 8: 332-339. 10.1016/j.micinf.2005.06.031.
Tesović G, Lesnikar V: Aseptic meningitis after vaccination with L-Zagreb mumps strain – virologically confirmed cases. Vaccine. 2006, 24: 6371-6373. 10.1016/j.vaccine.2006.06.025.
Tesović G, Poljak M, Lunar MM, Kocjan BJ, Seme K, Vukić BT, Sternak SL, Cajić V, Vince A: Horizontal transmission of the Leningrad-Zagreb mumps vaccine strain: A report of three cases. Vaccine. 2008, 26: 1922-1925.
Chomczynski P, Mackey K: Single-step method of total RNA isolated by acid guanidine phenol extraction. Cell Biology: A Laboratory Handbook. Edited by: Celis JE. 1998, 221-224. 2
Sambrook J, Russell DW: Molecular cloning: a laboratory manual. 2001, Cold pring Harbor Laboratory Press, New York, third
Ivancic-Jelecki J, Baricevic M, Santak M, Forcic D: Restriction enzyme cleavage of fluorescently labeled DNA fragments-analysis of the method and its usage in examination of digestion completeness. Anal Biochem. 2006, 349: 277-284. 10.1016/j.ab.2005.11.041.
This work was supported by the Ministry of science, education and sports of the Republic of Croatia, grant 021-0212432-3123 (to M.S.). We thank J. Ivancic-Jelecki for critical reading of the manuscript.
The authors declare that they have no competing interests.
TKG participated in the conception of the study, performed the majority of the experiments and wrote the manuscript. DF helped in the conception of the study, its design and coordination and helped to draft the manuscript. MS participated in the molecular cloning, the nucleic acid sequencing and sequence alignment and helped to draft the manuscript. SML and AR prepared all mumps virus samples. RM participated in study design and helped to draft the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Kosutic-Gulija, T., Forcic, D., Šantak, M. et al. Genetic heterogeneity of L-Zagreb mumps virus vaccine strain. Virol J 5, 79 (2008). https://doi.org/10.1186/1743-422X-5-79