SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

1 Southeastern European Medical Journal, 2021; 5(1) 
 

Original article 

Imported Infections Versus Herd Immunity Gaps; A Didactic 
Demonstration of Compartment Models Through the Example of a 
Minor Measles Outbreak in Hungary 1 

Katalin Böröcz *1, Ákos Markovics 2, Zsuzsanna Csizmadia 1, Joseph Najbauer 1, Timea Berki 1, Peter 
Németh 1 

1 Department of Immunology and Biotechnology, Clinical Centre, University of Pécs Medical School, Pécs, 
Hungary 

2 Department of General and Physical Chemistry, Faculty of Natural Sciences, University of Pécs, Pécs, 
Hungary 

*Corresponding author: Katalin Böröcz, borocz.katalin@pte.hu 

                                                      

Received: Feb 1, 2021; revised version accepted: Mar 15 2021; published: Apr 28, 2021 
  
KEYWORDS: MMR, vaccine, humoral antibody, epidemics, SEIR model 
 

Abstract 
Introduction: In Hungary, where MMR vaccine coverage is 99%, in 2017, a minor measles epidemic 
started from imported cases due to two major factors – latent susceptible cohorts among the 
domestic population and the vicinity of measles-endemic countries. Suspended immunization 
activities due to the COVID-19 surge are an ominous precursor to a measles resurgence.  This 
epidemiological demonstration is aimed at promoting a better public understanding of 
epidemiological data. 
Materials and Methods: Our previous MMR sero-epidemiological measurements (N of total measles 
cases = 3919, N of mumps cases = 2132, and N of rubella cases = 2132) were analyzed using open-
source epidemiological data (ANTSZ) of a small-scale measles epidemic outbreak (2017, Hungary). 
A simplified SEIR model was applied in the analysis. 
Results: In case of measles, due to a cluster-specific inadequacy of IgG levels, the cumulative 
seropositivity ratios (measles = 89.97%) failed to reach the herd immunity threshold (HIT Measles = 
92–95%). Despite the fact that 90% of overall vaccination coverage is just slightly below the HIT, 
unprotected individuals may pose an elevated epidemiological risk. According to the SEIR model, 
≥74% of susceptible individuals are expected to get infected. Estimations based on the input data of 
a local epidemic may suggest an even lower effective coverage rate (80%) in certain clusters of the 
population. 
Conclusion: Serological survey-based, historical and model-computed results are in agreement. A 
practical demonstration of epidemiological events of the past and present may promote a higher 
awareness of infectious diseases. Because of the high R0 value of measles, continuous large-scale 
monitoring of humoral immunity levels is important. 
 

(Böröcz K, Markovics Á, Csizmadia Z, Najbauer J, Berki T, Németh P. Imported Infections Versus Herd 
Immunity Gaps; A Didactic Demonstration of Compartment Models Through the Example of a Minor 
Measles Outbreak in Hungary. SEEMEDJ 2021; 5(1); 1-16) 

 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

2 Southeastern European Medical Journal, 2021; 5(1) 
 

Introduction 

Testing of acquired immunity and effectiveness 
of vaccination against infectious diseases has 
been increasingly important in the design of 
preventive public health strategies. Resurgence 
in measles cases in the United States and across 
Europe has occurred, including in individuals 
vaccinated with two doses of the vaccine (1). 
According to the Centers for Disease Control 
and Prevention (CDC), World Health 
Organization (WHO) and United Nations 
International Children’s Emergency Fund 
(UNICEF), measles has already been a global 
issue and now it has been aggravated by 
disrupted immunization protocols due to the 
COVID-19 pandemic (2–5). All six WHO regions 
have reported disrupted immunization activities, 
with major adverse effects on routine 
immunization and mass vaccination campaigns 
(4). According to CDC reports, in 2020, more than 
117 million children were at the risk of missing 
out on measles vaccines as a consequence of 
the COVID-19 surge (2). Measles immunity gaps 
resulting from suspended immunization 
activities are an ominous precursor to a measles 
resurgence (4). In Ukraine, one of Hungary’s 
neighboring countries that was already endemic 
for measles, vaccination has been interrupted in 
many regions (3). Regarding Europe, ECDC 
surveillance data have indicated an 
exceptionally high number of measles cases in 
2018, 2019 and 2020 in EU/EEA countries.  

Vaccination remains one of the safest and most 
effective interventions available in public health 
for the primary prevention of infectious diseases, 
resulting in both direct and indirect immunity in 
individuals vaccinated (herd immunity) (6–8). 
Even though a safe and effective two-dose 
measles/MMR vaccination schedule has been 
available in Europe since the 1960s, maintaining 
high vaccine coverage is still difficult, despite 
the fact that in Hungary, the MMR vaccine is 
mandatory and consequently the vaccine 
coverage is estimated to be at 98-99%. 
According to our previous publications (9,10) and 
in agreement with the results obtained by our 
colleagues (11), there are latent immunization 
gaps in certain age (or immunization) clusters of 

the Hungarian population, with predominance of 
the ~35-45-year-old adults. These are 
individuals who form a significant portion of the 
active labor force of the country, for instance 
health care workers (HCWs).  

Between January 2017 and May 2019, there were 
76 reported measles cases in Hungary (12), 54 of 
which were reported between 21 February and 
22 March 2017 (13). Because of the recent 
outbreaks worldwide, not only of measles, but 
also mumps and rubella (MMR) infections, and 
because of waning of immunity over time after 
vaccination (14–17), the importance of 
continuous MMR seroepidemiological screening 
is evident. 

Suboptimal vaccine effectiveness in certain 
clusters of the population has a negative impact 
on overall vaccination coverage. Small-scale 
outbreaks suggest that certain measles 
vaccines – applied during the early phases of the 
Hungarian vaccination history – failed to elicit 
the desired immunological response. The 
resulting immunization gap(s) raise the concern 
of potential further outbreaks (9,11). The 2020 
COVID-19 outbreak called attention to the 
importance of mathematical modelling of 
epidemics (18). Based on a reliable model, the 
timescale and economic impact of the disease 
can be predicted and preventive 
countermeasures can be taken (19). Through the 
example of the measles epidemic in Makó (2017, 
southeast Hungary), we demonstrated that, in 
possession of key epidemiological data (e.g. R0 
value, estimated vaccination coverage of a given 
population, number of infected and recovered 
individuals and duration of the epidemic), a 
simple open-source mathematical model can 
give a good approximation of the course of an 
infection and may provide better general 
compliance with protective measures. 

Materials and Methods  

Experimental work 

In this seroepidemiological survey, we 
combined the data from our previous findings 
with recent measurements, including anti-



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

3 Southeastern European Medical Journal, 2021; 5(1) 
 

measles, -mumps and -rubella antibody level 
(IgG) determination. Measurements were 
performed on the automated Siemens BEP 2000 
Advance® platform (Siemens AG, Germany), 
using our self-developed ELISA assays 
validated by well-established commercial kits, 
as previously described (9,10). Indirect 
immunofluorescent microscopy was used a 
reference (Euroimmun, Germany).  

In case of large-scale seroepidemiological 
measurements, a serum bank consisting of 
anonymous patient sera was used (N of total 
measles cases = 3919, N of mumps cases = 2132, 
and N of rubella cases = 2132) (Ethical License 
number 2015/5726). Nationally representative 
samples included randomly selected clinical 
residual samples, with the exclusion criteria of 
neonates, children under the vaccination age 
and severely immunocompromised patients. 
Samples were collected from the Department of 
Laboratory Medicine (University of Pécs, Clinical 
Centre). Serum samples were from all listed age 
groups participating in this study and they were 
categorized based on past changes introduced 
in measles and MMR immunization schedules. 
Age group determination was based on the 
landmarks in the history of measles and MMR 
vaccination schedules in Hungary (Figure 1). 
Human sera were stored in the accredited 
laboratory of the Department of Immunology 
and Biotechnology (University of Pécs, Medical 
School, Pécs, Hungary) according to quality 
assurance criteria (ISO 17025). 

Population-level result evaluation and 
seropositivity ratio assessment was performed 
in relation to the concept of herd immunity 
threshold (HIT) values (HIT Measles = 92–95%, 
HIT Mumps = 85–90%, HIT Rubella = 83–86). The 
study relies on the full virus antigen repertoire-
based indirect ELISA method. Therefore, it must 
be considered a good surrogate, rather than an 
absolute correlate marker for immunity – as far 
as Plotkin’s nomenclature is considered 
normative (20–22). We examined vaccination 
group-related infection- and vaccine-induced 
antibody titres using the following software: 
SPSS, Origin Pro, Excel. 

SEIR model example and input data  

A small-scale measles outbreak in Hungary in 
2017 raised questions about the vaccination 
coverage rate in the country. Experimental 
results supported the theory of ineffective 
vaccines, as previously mentioned (9). In spite of 
its limitations, it seemed reasonable to set up a 
SEIR model calculation in order to see whether 
a few percent decrease in effective vaccination 
could result in a local epidemic. To demonstrate 
the disease spread in a well-immunized 
population where latent immunity gaps may be 
present, input data were based on the data of 
the 2017 measles outbreak in Makó, southeast 
Hungary. The following parameters were used 
to perform the calculations: 

Population (N): The epidemic was linked to the 
small-town hospital. During that year, 65 
physicians were responsible for medical 
attendance of the estimated 30,000 inhabitants 
of Makó and the surrounding villages. In our 
model, a population of N = 400-1,000 people 
was assumed, including patients, health-care 
workers and family members. 

Number of infected individuals (I): A total of 29 
cases were reported. 

Incubation time and contagious period: The 
incubation time for measles ranges from 10 to 12 
days on average, an infected person can be 
contagious even 1-2 days before the first 
characteristic symptoms are visible, up to 4 days 
after the rashes appear. In our model, the 
incubation time (Tinc) was assumed to be 10 
days, whereas the contagious period (Tcont) 

parameters were determined by equations (5) 
and (6). 

Reproduction rate ranging from 12 to 18 can be 
found in the literature and both values were 
tested. The higher value is applicable to 
communities where no social distancing is 
present and the ratio of vaccinated or 
immunized inhabitants is low. In Central Europe, 
the use of the lower value seems more rational, 
although this specific epidemic was kept mainly 
in a hospital, where circumstances promote the 
spread of the infection. In this case, the start of 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

4 Southeastern European Medical Journal, 2021; 5(1) 
 

the outbreak was defined as the possible first 
day of the first patient’s infection, while the 
model was set to stop after the recovery of the 
last infected person. When it comes to large-
scale epidemics, a different approach is used. If 
no new cases are found after a certain period, 
the outbreak is over. This time period is usually 
determined by the incubation time, with a 
calculation method suggested by the WHO. 

Based on the ELISA antibody measurements, it 
can be assumed that only ~90% of the Hungarian 
population has effective immunization, which is 
under the theoretical 92-94% of HIT. In the 
model, 90% of seroprevalence was assumed, 
but lower values were also tested subsequently. 
No additional vaccinated (V) compartment was 
created and immunized individuals were treated 
as recovered. Vital dynamics was disregarded 
due to the short period of the epidemic. Death 
rate was not taken into consideration either, as 
no fatalities were observed during the 
Hungarian outbreak. Calculations were 

performed using Microsoft Excel Visual Basic 
Application (VBA), but the graphs were plotted in 
Origin. VBA is a built-in feature of the Microsoft 
Office Suite with several limitations, but its 
prevalence and the user-friendly computer 
language makes it suitable for educational 
purposes. 

Results 

Changes and historical data regarding 
epidemics in the Hungarian measles/MMR 
vaccination schedule (23–25) have been plotted 
on a timeline in order to evaluate 
seroepidemiological data accordingly. Figure 1 
shows changes in measles and MMR vaccination 
schedules in Hungary since the introduction of 
the vaccine (1969). High age-specific attack rates 
characterizing major epidemics (1980-81 and 
1988-89) along with 93%-99% of vaccine 
coverage evidence insufficiencies of the early 
vaccination program. 

 
Figure 1. Measles and MMR vaccination schedules in Hungary 
(a) Vaccination against measles was introduced in Hungary in 1969. (b) From 1969 to 1974, a single dose of measles 
vaccine was administered in mass campaigns to persons aged 9-27 months. (c) After vaccination was 
implemented, the incidence rate decreased until 1973-74, when large epidemics occurred primarily in 
unvaccinated 6-9-year-olds. (d) The recommended age for vaccination was 10 months until 1978, when it was 
changed to 14 months. (e) After the 1980-81 epidemic, persons born between 1973 and 1977, who received vaccine 
when the recommended age was 10 months, were revaccinated. (f) The 1988-89 epidemic mainly affected persons 
aged 17-21, who had been targeted to receive vaccine during mass campaigns in the first years of the vaccination 
program in Hungary. After 1989, children were re-vaccinated at the age of 11 with a monovalent measles vaccine 
in a scheduled manner. Also, in 1989, the rubella vaccine was introduced. (g) In 1990, measles-rubella bivalent 
vaccines were introduced. (h) The administration of the first vaccine at the age of 14 months lasted from 1978 to 
1991. Also, in 1991, the measles-mumps-rubella trivalent vaccine was introduced. (i) In 1992, the administration of 
the first MMR vaccine was shifted to 15 months of age. (j) In 1996, the MERCK MMR II vaccine (Enders’ Edmonston 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

5 Southeastern European Medical Journal, 2021; 5(1) 
 

strain, live attenuated) was introduced. (k) In 1999, measles-mumps-rubella revaccination replaced the monovalent 
measles vaccine. Also, in 1999, the GSK PLUSERIX vaccine (Measles Schwarz Strain) was introduced. (l) In 2003, the 
GSK PRIORIX vaccine was introduced. (m) Between 2004 and 2005, the MERCK MMR II vaccine was used. (n) 
Between 2006 and 2010, the GSK PRIORIX vaccine was in use. (o) Starting from 2011, we have been using a Sanofi-
MSD product, MMRvaxPro (Measles virus Enders’ Edmonston strain, live, attenuated), for vaccination and 
revaccination of children. GSK PRIORIX is still on the market, commonly used for vaccination in adulthood. (p) 
Between January 2017 and December 2019, there were 76 reported measles cases in Hungary (according to ECDC 
Surveillance reports). (Source of information: MMWR Weekly October 06, 1989 / 38(39); 665-668, International 
Notes Measles – Hungary, http://www.vacsatc.hu, https://www.ecdc.europa.eu) 

 

Figure 2 shows the age or vaccination group-
specific seropositivity and seronegativity ratios 
for measles, mumps and rubella. The lowest 
seropositivity ratios in terms of anti-measles 
antibody titres (IgG) were observed in the groups 
‘Vaccinated between 1969-1977’ (87.56%) and 
‘Vaccinated between 1978-1987’ (78.48%). These 
results are further confirmed by the 

abovementioned vaccine insufficiencies of the 
relative periods, described in Figure 1. Regarding 
the mumps and rubella seroepidemiological 
survey, in terms of humoral antibody levels, all 
vaccination groups satisfied the requirements 
necessary for the achievement of herd 
immunity.

 
 

(a) Measles 

 

 
(b) Mumps 

 

 
 

(c) Rubella 

 

 

Figure 2. Measles, mumps and rubella 
seropositivity ratios according to vaccination 
groups 
Age / vaccination groups: (I) Individuals born before 
1969. (II) Individuals vaccinated between 1969 and 1977. 
(III) Individuals vaccinated between 1978 and 1987. (IV) 
Individuals vaccinated between 1988 and 1990. (V) 
Individuals vaccinated between 1991 and 1995. (VI) 
Individuals vaccinated between 1996 and 1998. (VII) 
Individuals vaccinated between 1999 and 2002. (VIII) 
Individuals vaccinated in 2003. (IX) Individuals 
vaccinated between 2004 and 2005. (X) Individuals 
vaccinated between 2006 and 2010 (XI) Individuals 
vaccinated after 2011. The lowest seropositivity ratio 
(78.48%) was observed in the anti-measles antibody 
titres (IgG) in the group ‘Vaccinated between 1978 and 
1987’. 

In case of measles, mumps and rubella 
cumulative results, the seropositivity ratios were 
89.97%, 91.60% and 92.58%, respectively, as 
shown in Figure 3. Due to previously detailed 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

6 Southeastern European Medical Journal, 2021; 5(1) 
 

cluster-specific inadequacy of humoral antibody 
levels, the cumulative anti-measles 

seropositivity ratios also failed to reach the herd 
immunity threshold (HIT Measles = 92–95%).

 
Measles Mumps Rubella 

   
 Overall seropositivity ratio  
 Overall seronegativity ratio  

𝑆𝑒𝑟𝑜𝑝𝑜𝑠𝑖𝑡𝑣𝑖𝑡𝑦 𝑟𝑎𝑡𝑖𝑜 = 𝑇𝑜𝑡𝑎𝑙 𝑛𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑠𝑎𝑚𝑝𝑙𝑒𝑠 −
𝛴 (𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑒 + 𝑒𝑞𝑢𝑖𝑣𝑜𝑐𝑎𝑙 𝑠𝑎𝑚𝑝𝑙𝑒𝑠)

𝑇𝑜𝑡𝑎𝑙 𝑛𝑢𝑚𝑏𝑒𝑟 𝑜𝑓 𝑠𝑎𝑚𝑝𝑙𝑒𝑠
∗ 100 

HIT Measles = 92–95% HIT Mumps = 85–90% HIT Rubella = 83–86% 

Figure 3. Overall seropositivity and seronegativity ratios 
N measles = 3,919; N mumps, rubella = 2,132. In case of measles, mumps and rubella cumulative results, the 
seropositivity ratios were 89.97%, 91.60% and 92.58%, respectively. The overall ratio of seropositive samples was the 
lowest in the ‘measles’ group, where it remained under the threshold value. Seropositivity ratios were calculated as 
follows: 

Using the seronegativity ratio of 89.97% (≈ 90%) 
obtained by the cumulative data representation 
of anti-measles (IgG) antibody levels, the model 
of possible outcomes of a measles outbreak in a 
hospital as a function of the vaccination 
coverage rate was investigated. The results of 
the VBA-based SEIR model of the 2017 epidemic 

in Hungary are summarized in Table 1. Three 
parameters – population of the sample, ratio of 
immunized individuals and reproduction rate of 
the virus – were set to different values. The 
effect of these adjustments was investigated 
and changes in the number of measles cases 
and timescale of the epidemic were observed. 

Table 1. SEIR model results for the 2017 measles epidemic in Makó, Hungary 

10%

90%

8%

92%

7%

93%

Population of the 
sample (N) 

Ratio of immunised 
(%) 

Total number of 
measles cases 

Duration of epidemic 

𝑹𝟎 = 𝟏𝟖 
1000 90 73 6 months 
400 90 29 4 months 
400 80 78 3 months 
150 80 29 2.5 months 

𝑹𝟎 = 𝟏𝟐 
1000 90 2 6 days 
400 90 2 6 days 
400 80 70 4 months 
150 80 26 3 months 
Empirical values 
? 90 29 2 months 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

7 Southeastern European Medical Journal, 2021; 5(1) 
 

At R_0=18 and N = 1000, assuming 90% effective 
vaccination, 100 susceptible individuals can be 
found in the population. The model estimates a 
total number of infected persons at 74 and the 
duration of the epidemic at half a year, which is 
more than double of the real values. By setting 
the population at N = 400, 30 infected individuals 
and 4 months were given by the model. This way 
the number of infected persons corresponds to 
the actual clinical data, but the duration is still 
longer compared to empirical findings.  

Timescale of the epidemic can be compressed 
by increasing the proportion of susceptible 
people. If the vaccination coverage rate is 
changed from 90% to 80%, the duration of the 
epidemic is reduced to 3 months, but the total 
number of infected individuals becomes higher. 
Based on this anomaly, it can be presumed that 
the total number of involved population might 
be even lower than 400. Unfortunately, the 
results of the contact tracing procedure were 
not available for a better approximation.  

An acceptable correspondence between the 
model calculations and the clinical data was 
observed by assuming N = 150 and 80% of 
vaccination coverage as input parameters. 

The results – 30 infections in a two-month period 
– are close to the official values. For a better 
comparison, modelling with R_0=12 was also 
performed. The less contagious the virus, the 
fewer cases are found. Using this lower 
reproduction rate, only isolated cases can occur 
at 90% of vaccination coverage (which is a value 
that resembles the HIT). By decreasing the 
vaccination rate, the number of cases increases 
and the timescale is shortened, similarly to 
previous test examples. 

 

Discussion 

MMR vaccination in Hungary  

In Hungary, MMR vaccine is mandatory. A single-
dose, live-virus combined measles-mumps-
rubella (MMR) vaccine is used to vaccinate 
infants of ≥15 months of age. A reminder vaccine 
is given to sixth year primary school students (~11 
years of age). PRIORIX (GSK), PRIORIX-TETRA 
(GSK), ProQuad (MERCK) and the M-M-
RVAXPRO (MSD Pharma) vaccines are currently 
used in Hungary for vaccination of children (at 15 
months and 11 years of age) and for adults (62). 
The vaccines contain live attenuated viruses (26). 
Regarding insufficient cumulative anti-measles 
seropositivity levels, we would like to emphasize 
that potential gaps in the population-level 
humoral immunity (IgG) are attributable to early 
vaccination periods and are not a general 
phenomenon relative to the current 
immunization practices. The susceptibility of 
certain cohorts is likely attributable to the 
thermal instability of the historical Leningrad-16 
vaccine, inefficient seroconversion owing to 
vaccination at a premature age (e.g. 9 months of 
age) and the questionable efficiency of the 
inoculum itself (9, 11, 25, 29, 30, 31). The 2017 
measles outbreak in Makó was provoked by 
imported cases. Some of our bordering 
countries are still endemic for measles (27–30). 
Supplementary Figure 1 shows the European 
measles cases in the time period relative to the 
epidemics in Makó and Szeged. COVID-19 is 
increasing the risk of measles outbreaks. 
According to CDC Global Measles Outbreak 
reports of January 2021, 41 countries may 
postpone their measles campaigns for 2020 or 
2021 due to the COVID-19 pandemic. This 
increases the risk of bigger outbreaks around 
the world (31).



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

8 Southeastern European Medical Journal, 2021; 5(1) 
 

 

Supplementary Figure 1. European measles cases in the time period relative to the epidemics in 
Makó and Szeged (ecdc.europa.eu) 
Between December 2016 and November 2017, numerous measles cases occurred in Europe, most of which were 
reported by Romania, one of Hungary’s neighbouring countries. Source of data: https://www.ecdc.europa.eu/ 

 

2017 measles epidemic in Hungary 

In 2017, according to the data of the national 
authorities, a total of 76 persons were infected 
with measles (corrected to 73 laboratory 
confirmed cases by ECDC Surveillance reports). 
The outbreak in the hospital of the small town of 
Makó involved 29 individuals and lasted from 
January 2017 to March 2017 (32,33).  In order to 
demonstrate the spread of virus in a well-
immunized population, where despite good 
vaccination coverage, latent immunization gaps 
(unprotected, seronegative cohorts) are present, 

we used an open-source epidemiological report 
of the Hungarian National Public Health and 
Medical Officer Service (ANTSZ) (17 March 2017): 
‘At the peak of the Hungarian measles 
epidemics during the spring of 2017, 52 cases 
with measles-specific symptoms were reported. 
Of these, 15 laboratory confirmed cases 
(National Reference Laboratory for Measles and 
Rubella, National Public Health Institute, 
Budapest, Hungary) were registered by 16 
March. Of these patients, 12 were health care 
workers (HCWs) and two were hospitalized 
patients. One of them was a foreigner, while the 

86
245
86
7
3

137
5
1

13
296

647
389

36
3

19
3.891

0
2
4
0

18
1
43
29

2.257
1
6

157
30
226

8.638 5755

0 2000 4000 6000 8000 10000 12000 14000 16000

Austria

Belgium

Bulgaria

Croatia

Cyprus

Czech Republic

Denmark

Estonia

Finland

France

Germany

Greece

Hungary

Iceland

Ireland

Italy

Latvia

Lithuania

Luxembourg

Malta

Netherlands

Norway

Poland

Portugal

Romania

Slovakia

Slovenia

Spain

Sweden

United Kingdom

Total

Measles cases by notification rate (ECDC)
(December 2016 - November 2017; cumulated results)

Laboratory confirmed case Cases



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

9 Southeastern European Medical Journal, 2021; 5(1) 
 

other one was a patient living in the vicinity of a 
HCW. The epidemic affected two health care 
institutions, the Hospital of Makó and the clinics 
of the University of Szeged. The first measles 
case was imported in mid-February 2017 to the 
Hospital of Makó. The epidemic affected the 
hospital staff and their contacts. By 17 March, a 
measles infection was confirmed in case of a 
patient who was presumed to be the original 
importer of the virus, in case of 11 HCWs and in 
case of one of the HCW’s contacts. At the time 
of this report, additional 11 cases (of which seven 
HCWs and three patients’ contacts) were still 
under investigation. At the clinics of the 
University of Szeged, two persons – a HCW and 
a patient – fell ill with measles. Another 11 
persons (six patients and five HCWs) were also 
suspected at the time of the report. Following 
the appearance of the abovementioned 
measles cases, in Csongrád County, a total of 
391 people were vaccinated against measles, 
mumps and rubella (MMR). As the first cases of 
this period had been revealed, the National Chief 
Physician ordered strict monitoring and 
reporting of suspected measles virus infections. 
Thus, another 15 suspected cases were 
registered in several other counties. At the time 
of the report, laboratory testing was still ongoing 
(12)’. 

The second group of imported cases was 
detected at the end of July 2017 in Nyíregyháza, 
Szabolcs-Szatmár-Bereg County, Hungary (11). 
Six unvaccinated Romanian children were 
hospitalised with clinical symptoms of measles. 
These cases were later laboratory confirmed 
(National Reference Laboratory for Measles and 
Rubella, National Public Health Institute, 
Budapest, Hungary). The subsequent disease 
spread among two additional HCWs (also 
laboratory confirmed) supports the 
susceptibility of certain clusters in the Hungarian 
population (11). 

Epidemiological model- a didactic representation 

In this section, we explain the spreading 
mechanism of infectious diseases for those who 
are not familiar with the computational 
background of modelling. To understand the 
basics of epidemic models, a simplified 
mathematical interpretation can be used. The 
spread of a disease can be described by S-
shaped sigmoid mathematical functions, similar 
to the well-known pH titration curve, or 
haemoglobin saturation curves. As infectious 
diseases spread from human to human, the 
number of susceptible persons is decreasing 
over time and it influences the propagation of 
the pathogenic agent. In the beginning of the 
outbreak, the damping effect of recovered 
patients is minimal; the curve is very close to 
exponential and the number of new cases 
increases rapidly. At a certain time, a kind of 
equilibrium follows, daily recoveries can balance 
new infections and the curve reaches its 
inflection point. Afterwards, in the saturation 
phase, the epidemic slows down and at the end, 
no new cases are found and the vast majority of 
the population has recovered (Figure 4). The 
curves represented in Figure 4 are a graphic 
interpretation of a commonly used method for 
epidemic modelling – the compartment model. 
In this model, the population is divided into 
compartments – well-defined categories based 
on their epidemiological properties. In a 
compartment, all individuals behave exactly the 
same, e.g. they are all infected, all vaccinated, all 
exposed, etc. The simplest among these 
compartment models is the SIR model, where 
the letters of the acronym stand for susceptible, 
infectious and recovered. 

.

 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

10 Southeastern European Medical Journal, 2021; 5(1) 
 

0 30 60 90 120 150 180

0

200

400

600

800

1000

P
o
p

u
la

ti
o

n

Time (days)

Susceptible 
Recovered

Infected

94%

6%

8%

Figure 4. SIR model curves of a hypothetical 
epidemic 
As the disease spreads, the number of susceptible 
individuals decreases. First they get infected (I), but 
later on they will progress to the recovered 
compartment (R). Approximately 6% of the population 
managed to avoid contact with infected individuals. 
The peak of infections could be observed almost three 
months after the first case was recorded, affecting 8% 
of the population at the same time. 

The progression of an individual in this model is 
easy to follow, each member of the population 
progresses from susceptible to infectious to 
recovered. 

𝑆
𝛽
⇒ 𝐼

𝛾
⇒ 𝑅 (1) 

Transition between compartments is described 

probability of transmitting the disease between 
a susceptible and an infectious person. In other 

individuals to whom an infectious person can 
pass the disease per day (18,39,40) For example, 
if the infection rate is 0.2, it will take five days on 
average to infect someone. If we assume that 
the patient is contagious for 10 days, two new 
infections are expected in this case.  

The overall efficacy of the epidemic can be 
described by the number of these secondary 
infections originated from the primary infection, 
our first patient. This important parameter is the 

basic reproduction number (R0). Each virus has 
its own average R0 value – 12-18 for measles 
and 3.3-5.7 for COVID-19, according to the 
literature. 

transition into the recovered compartment. For 
instance, if this rate is 0.1, the contagious period 
will last for 10 days.  

From a mathematical perspective, the 
transitions can be described by the following 
differential equations, where S, I and R are the 
number of individuals in the corresponding 
compartments, while N is the whole population. 

dS

dt
= −

βIS

N
  (2)   

dI

dt
=

βIS

N
− γI(3) 

dR

dt
= γI(4) 

Mathematical methods (such as the Runge-
Kutta method) are available for solving similar 
equations, but there is a simpler option. Using 
the built-in features of Microsoft Excel (or any 
equivalent spreadsheet application), it is 
possible to make calculations using an iterative 
method. Instead of solving the equations, the 
computer performs calculations that follow the 
daily changes in different compartments. 

R0 have to be defined. Based 
on the definition of the transition rates, it can be 
seen that the recovery rate can be determined 
by the number of contagious days (T_cont). 

𝛾 =
1

𝑇𝑐𝑜𝑛𝑡
 (5) 

 

The basic reproduction number can be given as 
follows: 

𝑅0 =
𝛽

𝛾
 (6) 

Let us assume that in a certain population 
measles can be transmitted from a single person 
to 12 others (R_0=12) and they stay contagious for 
6 days (T_cont=6). In this case: 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

11 Southeastern European Medical Journal, 2021; 5(1) 
 

𝛾 =
1

𝑇𝑐𝑜𝑛𝑡
=

1

6
 (7) 

𝛽 = 𝑅0𝛾 = 2 (8) 

Incubation time plays an important role in the 
spread of a disease. In a more sophisticated 
model (SEIR model), this can also be taken into 
consideration. A new compartment for the 
exposed part of the population can be 
generated. The susceptible person first gets 
exposed and will progress to the infectious state 
only after a certain time. 

𝑆
𝛽
⇒𝐸

𝛼
⇒𝐼

𝛾
⇒ 𝑅 (9) 

The parameter ‘α’ is a new transition rate, 
which can be determined by the incubation time 
(T_inc) , similarly to γ: 

𝛼 =
1

𝑇𝑖𝑛𝑐
 

New compartments can be added to the model 
anytime, such as the compartment M for 
individuals with maternal immunity or the 
compartment E for exposed individuals, who are 
already infected, but not infectious. Based on the 
characteristics of certain infectious diseases, 
further models have been developed, such as 
the SIS, MSIR, SEIR, SEIS, MSEIR and MSEIRS 
models. The second ‘S’ in the acronym indicates 
that after the infection, no permanent immunity 
can be reached and the individuals step to the S 
compartment again. In other models, the ratio of 
hospitalization, the ratio of mild and severe 
cases and epidemiological interventions can be 
included, with a more complex mathematical 
background. 

In the examples described above some 
important parameters are simply disregarded, 
although it is possible to perform a more 
detailed computation. Vital dynamics, the 
natural dynamics of birth and death, can be 
included by adding two further parameters.  

It is necessary to mention that compartment 
models have their well-known limitations and 
shortcomings. For instance, all individuals in the 
population are assumed to have an equal 
probability of coming in contact with others, 

although society is inhomogeneous from the 
perspective of social distancing. Another 
drawback is that the traditional compartment 
model cannot handle uncertainty in model 
parameters. Working with a smaller set of data 
increases this uncertainty, making predictions 
unreliable. To overcome this problem, it is usual 
to calculate the SIR model over a few possible 
values for each parameter. A more complex 
solution is to use distribution functions instead of 
single numbers and if real-time data is available 
(e.g. we are in the middle of a pandemic), a 
clinical dataset can be utilized to adjust these 
parameters (36–38). 

Regarding the SEIR model resembling the 2017 
measles outbreak in Makó (Figure 4), we would 
like to note that both the simplified 
mathematical method and the input data were 
unreliable. With more sophisticated models, 
many different parameters can be taken into 
consideration (37,39). Despite that fact, the 
calculated values correspond in order of 
magnitude to the available data on the epidemic 
and support the experimental results describing 
the vaccination gap.  

Model curves using a lower percentage of the 
population-level anti-measles protection rate 
are more fitting. This finding may indicate an 
even lower percentage of effectively vaccinated 
population than it was found previously (~90%).  

It is concluded that the importance of 
seroepidemiological surveys is confirmed by the 
recent outbreaks of measles, mumps and 
rubella infections in several countries 
(14,16,17,40–45). Considering the HIT values, 
suboptimal anti-measles seropositivity ratios 
were detected in certain clusters of the early 
vaccination era (78.48% of sufficient anti-
measles IgG antibody titres among individuals 
vaccinated between 1978 and 1987). This finding, 
which is in accordance with a recent study 
published by our colleagues (11) and historical 
literature data (46), suggests the existence of 
age-specific immunization gaps in the 
Hungarian population. For mumps and rubella, 
our preliminary data shows satisfactory 
immunity levels. Nowadays, in our country, the 
MMR vaccination coverage is ideal due to the 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

12 Southeastern European Medical Journal, 2021; 5(1) 
 

mandatory administration of safe and modern 
trivalent vaccines. Nevertheless, dubious 
immunization practices in some of our 
neighboring countries, aggravated by the 
detrimental effect of the COVID-19 pandemic 
and subsequent suspension of measles 
vaccination campaigns, may facilitate the 
occurrence of minor importation-related MeV 
outbreaks in susceptible cohorts. Using the 
example of the 2017 measles outbreak in Makó, 
it has been demonstrated that in possession of 
key epidemiological parameters (e.g. R0 value, 
estimated vaccination coverage of a given 
population, number of infected and recovered 
individuals, duration, etc.), a simple SEIR model 
can give a good approximation regarding the 
course of an infection.  

We believe that awareness may significantly 
reduce the extent of an epidemic (38,47). In the 
light of current disquieting epidemiological 
circumstances, we suggest the introduction of 
open-access mathematical and epidemiological 
models into modern natural science education 
of students. Today, online epidemic models are 
easily available for the public (35,36). Practical 
introduction to these plain calculation models 
could help students understand the rationale 
behind epidemiological data. We believe that a 
practical demonstration of epidemiological 
events can promote a better understanding of 
countermeasures and also allow for an easier 
adaptation to the current epidemiological 
regulations. 

Limitations of experimental work 

The diagnostic ability of our assay was 
calculated based on the results obtained by 
well-established kits capable of humoral 
antibody detection, rather than on neutralizing 
antibody titres that could serve as an absolute 
correlate of protection (48–50). It is important to 
emphasize that immunity to measles is a 

complex orchestration between the cellular and 
humoral immunity. For this reason, only 
antibody-based definitions of vaccine success 
and failure may be misleading, or at least 
simplistic and incomplete (51). 

Limitations of mathematical modelling 

Input data plays a key role in modelling of 
epidemics. Even when the number of cases is 
high – like in the 2020 COVID outbreak – the 
confidence of fitting is poor at the beginning of 
new cases vs. time graph. The first cases are 
usually unexpected, quarantine and social 
distancing protocols are not applied yet and if 
the disease has a low prevalence in the 
population, the accuracy of the diagnosis might 
be low. Besides that, atypical symptoms can be 
misleading for physicians. Furthermore, 
statistical values, such as basic reproduction 
number, incubation and recovery time, depend 
on other factors, such as social distancing and 
the health care system. 

Acknowledgement. None. 

Disclosure 
Funding. This work was realized with the 
financial support of the University of Pécs, 
Medical School (Pécs, Hungary) and the 
following grants: EFOP-3.6.1-16-2016-00004, 
2019 Higher Education Excellence Program 
(FIKP II PoP) and HUHR/1901/3.1.1/0032 
CABCOS3. Project no. TKP2020-IKA-08 has been 
implemented with the support provided from 
the National Research, Development and 
Innovation Fund of Hungary, financed under the 
2020-4.1.1-TKP2020 funding scheme.” 
Transparency declaration: We declare that we 
have no commercial or potential competing 
interests or any financial and personal 
relationships with other people or organizations 
that could inappropriately influence our 
network. 
 

References 

1.  Haralambieva IH, Kennedy RB, 
Ovsyannikova IG, Whitaker JA, Poland GA. 
Variability in Humoral Immunity to Measles 
Vaccine: New Developments. Trends Mol Med 

2015;21(12):789-801. doi: 
10.1016/j.molmed.2015.10.005.  

2.  COVID-19’s Impact on Measles 
Vaccination Coverage. Available from: 
https://www.cdc.gov/globalhealth/measles/n



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

13 Southeastern European Medical Journal, 2021; 5(1) 
 

ews/covid-impact-on-measles-
vaccination.html 

3.  Suspension of vaccination due to COVID-
19 increases the risk of infectious diseases 
outbreaks - UNICEF and WHO. Available from: 
https://www.unicef.org/ukraine/en/press-
releases/suspension-vaccination-due-covid-
19-increases-risk-infectious-diseases-
outbreaks (press release) 

4.  Durrheim DN, Andrus JK, Tabassum S, 
Bashour H, Githanga D, Pfaff G. A dangerous 
measles future looms beyond the COVID-19 
pandemic. Nat Med 2021; 27:360–361. 
https://doi.org/10.1038/s41591-021-01237-5 

5.  WHO Western Pacific Region. WHO 
Western Pacific Region: Countries and areas. 
Measles-Rubella Bull. 2021;15(1):1–10. 
https://iris.wpro.who.int/bitstream/handle/106
65.1/14713/Measles-Rubella-Bulletin-2021-Vol-
15-No-01.pdf 

6.  Torre G La, Saulle R, Unim B, Meggiolaro 
A, Barbato A, Mannocci A, Spadea A. The 
effectiveness of measles-mumps-rubella (MMR) 
vaccination in the prevention of pediatric 
hospitalizations for targeted and untargeted 
infections: A retrospective cohort study. Hum 
Vaccin Immunother. 2017;13(8):1879- 1883. doi: 
10.1080/21645515.2017.1330733  

7.  Rosian-Schikuta I, Fröschl B, Habl C, 
Stürzlinger H. The Measels-Mumps-Rubella 
Vaccination from a health political and 
economical point of view. GMS Health Technol 
Assess. 2007 Nov 28;3:Doc12.  

8.  Plans-Rubió P. Evaluation of the 
establishment of herd immunity in the 
population by means of serological surveys and 
vaccination coverage. Hum Vaccin Immunother. 
2012; 8(2):184–8.  

9.  Böröcz K, Csizmadia Z, Markovics, Farkas 
N, Najbauer J, Berki T, Nemeth P. Application of 
a fast and cost-effective “three-in-one” MMR 
ELISA as a tool for surveying anti-MMR humoral 
immunity: The Hungarian experience. Epidemiol 
Infect. 2020; 148: e17.  doi: 
10.1017/S0950268819002280 

10.  Böröcz K, Csizmadia Z, Markovics Á, 
Mészáros V, Farkas K, Telek V, Varga V, Ouma 
Maloba  G, Bodó K, Najbauer  J, Berki T, Németh  
P. Development of a robust and standardized 
immunoserological assay for detection of anti-
measles IgG antibodies in human sera. J 
Immunol Methods. 2019 Jan; 464:1-8. doi: 
10.1016/j.jim.2018.07.009.  

11.  Lengyel G, Marossy A, Ánosi N, Farkas SL, 
Kele B, Nemes-Nikodém É, et al. Screening of 
more than 2000 Hungarian healthcare workers’ 
anti-measles antibody level: results and 
possible population-level consequences. 
Epidemiol Infect. 2019 Sep 11;147:e7.  

12.  ANTSZ - Tájékoztatás a kanyarós 
megbetegedésekről. Available from: 
https://www.antsz.hu/felso_menu/temaink/ja
rvany/jarvany_archivum/kanyaro/tajekoztatas
kanyarosmegbetegedesekrol.html 

13.  Epidemiological update: Measles - 
monitoring European outbreaks, 7 July 2017. 
Available from: 
https://www.ecdc.europa.eu/en/news-
events/epidemiological-update-measles-
monitoring-european-outbreaks-7-july-2017 

14.  Lambert N, Strebel P, Orenstein W, 
Icenogle J, Poland GA. Rubella. Lancet. 2015; 
385(9984):2297–307.  

15.  Moss WJ. Measles. Lancet. 2017; 
390(10111):2490–502.  

16.  Hviid A, Rubin S, Mühlemann K. Mumps. 
Lancet. 2008;371(9616):932–44.  

17.  Lewnard JA, Grad YH. Vaccine waning 
and mumps re-emergence in the United States. 
Sci Transl Med. 2018; 10(433):eaao5945.  

18.  Peng L, Yang W, Zhang D, Zhuge C, Hong 
L. Epidemic analysis of COVID-19 in China by 
dynamical modeling. 2020. arXiv:2002.06563 [q-
bio.PE]  

19.  Kretzschmar M. Epidemic Model - an 
overview | ScienceDirect Topics. International 
Encyclopedia of Public Health (Second Edition). 
2017. p. 571–84. 
https://www.sciencedirect.com/topics/medici
ne-and-dentistry/epidemic-model 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

14 Southeastern European Medical Journal, 2021; 5(1) 
 

20.  Böröcz K, Csizmadia Z, Markovics Á, 
Mészáros V, Farkas K, Telek V, Varga V, Maloba  
GO, Bodó K, Najbauer J, Berki  T, Németh P. 
Development of a robust and standardized 
immunoserological assay for detection of anti-
measles IgG antibodies in human sera. J 
Immunol Methods. 2019 Jan;464:1-8. doi: 
10.1016/j.jim.2018.07.009. 

21.  Plotkin SA. Correlates of protection 
induced by vaccination. Clin Vaccine Immunol. 
2010; 17(7):1055–65. Available from: 
http://www.ncbi.nlm.nih.gov/pubmed/204631
05 

22.  Plotkin SA. Complex correlates of 
protection after vaccination. Clin Infect Dis. 
2013;56(10):1458–65.  

23.  Plotkin SA, Gilbert PB. Nomenclature for 
immune correlates of protection after 
vaccination. Clin Infect Dis [Internet]. 2012 Jun 
[cited 2019 Aug 30];54(11):1615–7. Available from: 
http://www.ncbi.nlm.nih.gov/pubmed/224372
37 

24.  Böröcz K, Csizmadia Z, Markovics Á, 
Farkas N, Najbauer J, Berki T, et al. Application of 
a fast and cost-effective “three-in-one” MMR 
ELISA as a tool for surveying anti-MMR humoral 
immunity: The Hungarian experience. Epidemiol 
Infect. 2020;  

25.  Agócs MM, Markowitz LE, Straub I, 
Dömök I. The 1988-1989 measles epidemic in 
Hungary: assessment of vaccine failure. Int J 
Epidemiol [Internet]. 1992 Oct [cited 2019 May 
5];21(5):1007–13. Available from: 
http://www.ncbi.nlm.nih.gov/pubmed/146883
7 

26.  MMWR Weekly. MMWR Publications | 
MMWR [Internet]. [cited 2019 May 5]. Available 
from: 
https://www.cdc.gov/mmwr/publications/ind
ex.html 

27.  Hungary V. Vacsatc - Oltásbiztonság | 
Nyitólap [Internet]. [cited 2020 Mar 31]. Available 
from: http://www.oltasbiztonsag.hu/ 

28.  Vacsatc - Oltásbiztonság | Szakmai 
oldalak | Virális oltóanyagok [Internet]. [cited 

2019 Dec 21]. Available from: 
http://www.vacsatc.hu/?Vir%E1lis-
olt%F3anyagok&pid=37 

29.  International Notes Measles -- Hungary 
[Internet]. Available from: 
https://www.cdc.gov/mmwr/preview/mmwrh
tml/00001472.htm 

30.  Böröcz K, Csizmadia Z, Markovics Á, 
Mészáros V, Farkas K, Telek V, et al. 
Development of a robust and standardized 
immunoserological assay for detection of anti-
measles IgG antibodies in human sera. J 
Immunol Methods. 2019;464.  

31.  Lengyel G, Ánosi N, Marossy A, Mátyus M, 
Bosnyákovits T, Orosz L. Az európai és a 
magyarországi kanyaróhelyzet összefoglalása 
és tanulságai. Orv Hetil [Internet]. 2019 May 12 
[cited 2019 Jul 31];160(20):767–73. Available from: 
https://www.akademiai.com/doi/10.1556/650.
2019.31397 

32.  Habersaat KB, Pistol A, Stanescu A, 
Hewitt C, Grbic M, Butu C, et al. Measles 
outbreak in Romania: understanding factors 
related to suboptimal vaccination uptake. Eur J 
Public Health [Internet]. 2020 Oct 1 [cited 2021 
Mar 5];30(5):986–92. Available from: 
https://academic.oup.com/eurpub/article/30
/5/986/5847956 

33.  Dascalu S. Measles Epidemics in 
Romania: Lessons for Public Health and Future 
Policy [Internet]. Vol. 7, Frontiers in Public Health. 
Frontiers Media S.A.; 2019 [cited 2021 Mar 5]. p. 
98. Available from: 
/pmc/articles/PMC6496956/ 

34.  The Lancet. Measles, war, and health-
care reforms in Ukraine. Vol. 392, The Lancet. 
Lancet Publishing Group; 2018. p. 711.  

35.  Wadman M. Measles epidemic in Ukraine 
drove troubling European year [Internet]. Vol. 
363, Science. American Association for the 
Advancement of Science; 2019 [cited 2021 Mar 
5]. p. 677–8. Available from: 
http://stm.sciencemag.org/content/scitransm
ed/10/433/eaao5945.full 



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

15 Southeastern European Medical Journal, 2021; 5(1) 
 

36.  Global Measles Outbreaks [Internet]. 
[cited 2021 Jan 31]. Available from: 
https://www.cdc.gov/globalhealth/measles/
data/global-measles-outbreaks.html 

37.  ECDC. Epidemiological update: Measles - 
monitoring European outbreaks, 7 July 2017 
[Internet]. [cited 2020 Mar 31]. Available from: 
https://www.ecdc.europa.eu/en/news-
events/epidemiological-update-measles-
monitoring-european-outbreaks-7-july-2017 

38.  ANTSZ. ANTSZ - Public information on 
measles epidemics [Internet]. [cited 2020 Mar 31]. 
Available from: 
https://www.antsz.hu/felso_menu/temaink/ja
rvany/jarvany_archivum/kanyaro/tajekoztatas
kanyarosmegbetegedesekrol.html 

39.  Stella R. Quah WC. International 
Encyclopedia of Public Health - Second Edition-
Google Books [Internet]. Elsevier. 2017 [cited 
2020 Apr 4]. Available from: 
https://books.google.hu/books?hl=hu&lr=&id=
WAnpCgAAQBAJ&oi=fnd&pg=PP1&dq=Measure
ment+and+Modeling:+Infectious+Disease+Modeli
ng+Mirjam+Kretzschmar,+in+International+Encycl
opedia+of+Public+Health+(Second+Edition),+2017
&ots=k662f-7bGN&sig=CQ-q-
rJ9YMRIG9BDwTm 

40.  Epidemic Calculator [Internet]. [cited 
2020 Apr 4]. Available from: 
https://gabgoh.github.io/COVID/index.html 

41.  Buckman SR, Glick R, Lansing KJ, 
Petrosky-Nadeau N, Seitelman LM, Francisco S. 
Replicating and Projecting the Path of COVID-19 
with a Model-Implied Reproduction Number. 
2020;2020–44.  

42.  Röst G, Wu J. SEIR EPIDEMIOLOGICAL 
MODEL WITH VARYING INFECTIVITY AND 
INFINITE DELAY [Internet]. Vol. 5, AND 
ENGINEERING. 2008 [cited 2020 Apr 4]. Available 
from: http://www.mbejournal.org/ 

43.  Li M, Wang M, Xue S, Ma J. The influence 
of awareness on epidemic spreading on random 
networks. J Theor Biol. 2020 Feb 7;486:110090. 
doi: 10.1016/j.jtbi.2019.110090.  

44.  Lipsitch M, Cohen T, Cooper B, Robins 
JM, Ma S, James L, et al. Transmission dynamics 
and control of severe acute respiratory 
syndrome. Science (80- ). 2003 Jun 
20;300(5627):1966–70.  

45.  Ujiie M. Rubella Resurgence in Japan 
2018–2019. J Travel Med [Internet]. 2019 Sep 2 
[cited 2021 Feb 1];26(6). Available from: 
https://academic.oup.com/jtm/article/doi/10.
1093/jtm/taz047/5523888 

46.  WHO | New measles surveillance data 
from WHO. WHO [Internet]. 2019 [cited 2019 Aug 
22]; Available from: 
https://www.who.int/immunization/newsroom
/new-measles-data-august-2019/en/ 

47.  Wadman M. Measles epidemic in Ukraine 
drove troubling European year. Science (80- ) 
[Internet]. 2019 Feb 15 [cited 2019 May 
16];363(6428):677–8. Available from: 
http://www.sciencemag.org/lookup/doi/10.11
26/science.363.6428.677 

48.  The Lancet. Measles, war, and health-
care reforms in Ukraine. Lancet [Internet]. 2018 
Sep 1 [cited 2019 May 16];392(10149):711. 
Available from: 
http://www.ncbi.nlm.nih.gov/pubmed/3019181
2 

49.  Ukraine: Measles Outbreak - Jan 2019 | 
ReliefWeb [Internet]. [cited 2019 Aug 22]. 
Available from: 
https://reliefweb.int/disaster/ep-2019-
000017-ukr 

50.  Balbi AM, Van Sant AA, Bean EW, Jacoby 
JL. Mumps: Resurgence of a once-dormant 
disease. J Am Acad Physician Assist [Internet]. 
2018 May 1 [cited 2021 Feb 1];31(5):19–22. 
Available from: 
https://pubmed.ncbi.nlm.nih.gov/29642091/ 

51.  Takahashi S, Sato K, Kusaka Y, Hagihara 
A. Public preventive awareness and preventive 
behaviors during a major influenza epidemic in 
Fukui, Japan. J Infect Public Health. 2017 Sep 
1;10(5):637–43.  

52.  Liu Y, Liu Z, Deng X, Hu Y, Wang Z, Lu P, 
et al. Waning immunity of one-dose measles-



SEEMEDJ 2021, VOL 5, NO. 1 Imported Infections Versus Herd Immunity Gaps 

16 Southeastern European Medical Journal, 2021; 5(1) 
 

mumps-rubella vaccine to mumps in children 
from kindergarten to early school age: a 
prospective study. Expert Rev Vaccines 
[Internet]. 2018 May 4 [cited 2019 May 
16];17(5):445–52. Available from: 
http://www.ncbi.nlm.nih.gov/pubmed/294783
47 

53.  Bankamp B, Hickman C, Icenogle JP, Rota 
PA. Successes and challenges for preventing 
measles, mumps and rubella by vaccination. 
Curr Opin Virol [Internet]. 2019 Feb [cited 2019 
May 15];34:110–6. Available from: 
http://www.ncbi.nlm.nih.gov/pubmed/308524
25 

54.  Kontio M, Jokinen S, Paunio M, Peltola H, 
Davidkin I. Waning antibody levels and avidity: 
Implications for MMR vaccine-induced 
protection. J Infect Dis. 2012;206(10):1542–8.  

55.  Ovsyannikova IG, Dhiman N, Jacobson 
RM, Vierkant RA, Poland GA. Frequency of 
measles virus-specific CD4+ and CD8+ T cells in 
subjects seronegative or highly seropositive for 
measles vaccine. Clin Diagn Lab Immunol 
[Internet]. 2003 May [cited 2019 Aug 2];10(3):411–
6. Available from: 
http://www.ncbi.nlm.nih.gov/pubmed/127386
40

 

1 Author contribution. Acquisition of data: Böröcz K, 
Markovics Á, Csizmadia Z  
Administrative, technical or logistic support: Böröcz K, 
Csizmadia Z 
Analysis and interpretation of data: Böröcz K, 
Markovics Á, Najbauer J, Németh P 
Conception and design: Böröcz K, Markovics Á, 
Csizmadia Z, Berki T, Németh P 
Critical revision of the article for important intellectual 
content: Böröcz K, Najbauer J, Berki T, Németh P 
Drafting of the article: Böröcz K, Markovics Á, 
Csizmadia Z, Najbauer J 
Final approval of the article: Najbauer J, Németh P. 
Guarantor of the study: Berki T, 
Obtaining funding: Berki T, Németh P 
Provision of study materials or patients: Böröcz K,  
Statistical expertise (statistical analysis of data): 
Markovics Á, 

http://www.ncbi.nlm.nih.gov/pubmed/12738640
http://www.ncbi.nlm.nih.gov/pubmed/12738640