- Open Access
Periodicity of varicella-zoster virus in the presence of immune boosting and clinical reinfection with varicella
© Korostil et al.; licensee BioMed Central. 2015
Received: 2 December 2014
Accepted: 18 March 2015
Published: 11 April 2015
Clinical reinfection with varicella is normally ignored in mathematical transmission models as it is considered too rare to be important.
We apply basic bifurcation analysis to a simple mathematical model of varicella-zoster virus (VZV) transmission incorporating reinfection.
We demonstrate that under certain conditions this model can exhibit periodic behaviour as opposed to what is observed in VZV models that ignore the possibility of repeat varicella attacks. Periodicity can be induced by a combination of immune boosting and reinfection while the impact of zoster (shingles) recurrence on the onset of periodicity is negligible.
Our results suggest that mathematical models of VZV may benefit from inclusion of repeat varicella.
Varicella-zoster virus (VZV) is first manifested as varicella (also known as chickenpox), a rash that turns into open lesions. In the absence of vaccination, this usually happens when a person is young. While the symptoms clear within a week or so, VZV remains latent in the nervous system and may reactivate to cause a rash called herpes zoster (HZ) or shingles. Why this reactivation happens remains unclear, although it is generally assumed that onset is related to reduced levels of cellular immunity [1-3].
Clinical reinfection with varicella is thought to be rare and is typically not accounted for in epidemiological models of VZV transmission . However, observational studies suggest that this may be inaccurate, with a number of literature reports suggesting that reinfection occurs with moderate frequency. In particular, two US studies [5,6] based on active surveillance data indicate that a substantial proportion of varicella cases (between 4.5% and 25%) report a prior history of varicella infection. A smaller study demonstrated laboratory evidence of prior infection in 8 out of 115 (7%) patients with clinical varicella infection . There are also case reports of multiple varicella attacks in healthy children . Hence, it appears that prior infection does not always result in sustained protective immunity against subsequent varicella reinfection and there is a case for considering this process more closely in models.
Repeat HZ attacks in immunocompetent individuals have also typically been excluded in epidemiological modelling studies [4,9,10], under the assumption that repeat attacks primarily occur in immunocompromised individuals. However, second attacks were observed by Hope-Simpson  with 8 out of 192 cases recorded in his series being a repeat event. Repeat cases were also observed by Wilson in a rural UK practice with 6 out of 151 cases identified as repeat events . More recently, a US study found that 95 out of 1,669 individuals with a medically documented episode of HZ had two or more attacks (8 individuals had three or more attacks) . While immunocompromised patients had a higher incidence of recurrence, the rate of recurrence in immunocompetent patients was similar to the rate of first attacks reported in other studies.
The Hope-Simpson hypothesis that the immune response of someone who has been exposed to varicella can be strengthened (boosted) via subsequent exposure to VZV  has been supported by both clinical and modelling studies [3,13] although there are studies which failed to observe evidence of boosting . The degree of exposure that would be required to provide an immune boost remains unknown. The VZV modelling studies of which we are aware implement boosting by assuming that individuals with weakened immunity are boosted to full immunity at a fraction of the rate of at which they would be becoming infected with varicella (this rate is often referred to as a force of infection) [4,10,15-17]. In a recent work on the transmission of pertussis, the potential for the rate of boosting to exceed the force of infection was explored and shown to induce both cyclical behaviour and the potential for disease recurrence following immunisation . It was suggested that immune boost could occur even if the exposure dose is too small to be likely to cause clinical infection. This is supported by the observation that immune memory cells respond more rapidly and to lower doses of antigen when primed through previous exposure than in naive individuals [19,20].
In this paper, we consider a simple mathematical model of varicella-zoster virus transmission incorporating clinical reinfection with varicella, immune boosting requiring either more or less exposure to VZV than primary infection, and HZ recurrence. Our aim is to investigate whether the behaviour of this model qualitatively changes when its parameters are varied within ranges selected based on current knowledge of VZV. Specifically, we are interested in identifying conditions whereby there is switching from stable equilibrium behaviour (that modelling studies normally focus on) to periodic behaviour. While varicella is known to exhibit periodicity [3,21], the evidence for HZ periodicity is still unclear, with studies both supporting periodically-varying incidence [22,23] and failing to find evidence of it [24,25]. The majority of epidemiological modelling studies have not explicitly incorporated periodic VZV transmission (e.g. transmission rates modelled as sinusoid) or observed VZV-related periodicity, although extensive investigations of this periodicity have been conducted in the context of chaotic dynamics (see e.g. ).
Base value, refs.
Varicella transmission coefficient
Rate of becoming infectious with varicella
the base value corresponds to the average latent period of ≈ 14 days
Varicella recovery rate
the base value corresponds to the average infectious period of ≈ 7.3 days
Loss of full immunity rate
the base value corresponds to the average duration of “full” immunity of 10 years
Varicella force of infection reduction coefficient for those
who become reinfected following clearance of varicella
Varicella force of infection reduction coefficient for those
who become reinfected following HZ attack
the base value corresponds to the average life expectancy of ≈ 82 years
HZ reactivation rate
HZ recovery rate
the base value corresponds to the average infectious period of ≈ 7 days
Rate of HZ reactivation after loss of immunity following recovery from HZ
Rate of loss of immunity to VZV after recovery from HZ
the base value corresponds to immunity lasting 3 months on average
Relative VZV infectiousness of HZ
Hence, using the base values for ω,μ and γ provided in Table 1, we see that β≈325 year−1. The rate of loss of full immunity σ is unknown but attempts to estimate it have been made previously, hence we cite three relevant modelling studies. We set a base value for σ corresponding to the average duration of full immunity of 10 years and let σ take a range of possible values from 0.0125 (“full” immunity for 80 years) to 4.0 (“full” immunity for only 3 months). The boosting coefficient κ can be interpreted as a measure of how much exposure is required for immune boosting as compared with that required for primary infection. We allow it to vary from 0 to 20. The birth/death rate μ has been set to correspond to the current Australian life expectancy at birth . There are no readily available data to directly inform parameter p, which characterises the contribution of HZ to the varicella force of infection. Some VZV modelling studies have assumed that this contribution is so small that it can be ignored  while others attempted to estimate it. In particular, p=5.4×10−7 was estimated in  and , while p=3.11×10−6 was estimated in . A substantially larger value of 0.05 was used in . Hence, we allow p to vary between 0 and 0.1. The rate of HZ reactivation after loss of immunity following recovery from HZ δ is allowed to take values between 0 (no reactivation) and 0.1. Loss of immunity to VZV after recovery from HZ is characterised by parameter θ which ranges from 0 (no loss of immunity) to 12.0 (immunity lost after only 1 month). Finally, parameters ζ and α have been introduced to allow for reinfection with varicella. In the base case we assume that α=0.90 (i.e. the force of varicella infection on those who have recovered from HZ is reduced by only 10% as compared to that on someone susceptible to primary varicella infection) and allow α to take values from 0 (those who have had HZ cannot be reinfected with varicella) to 1. The base value for α is high not due to epidemiological considerations but because we determined that for substantially lower values of α system (1) does not deviate from stable equilibrium behaviour no matter what values the other parameters take. This is evident from the results presented in the next section. We set the base case value for ζ at 0.1, which means that the force of varicella infection on those who have recovered from varicella and have not had HZ is reduced by 90% as compared to that on individuals susceptible to primary infection.
In order to detect the periodic dynamics we systematically searched for Hopf points using the XPPAUT  bifurcation analysis software. A Hopf point is a point where a bifurcation from a branch of stable equilibrium to a branch of periodic oscillation occurs [33,34]. This point can be detected by varying a single model parameter while the other parameters are kept fixed. Having found the Hopf point, we can investigate where it would be if another parameter were given a different value. In this way we obtain a set of Hopf points which form a so-called Hopf curve.
We found that for certain parameters it was easy to obtain periodic behaviour for most or all of their values taken from within the ranges given in Table 1. Specifically, for all values of β under 400 we could normally observe oscillations, and this was also the case for β under 500 when η and/or σ were increased or μ and/or ζ were decreased as compared to their base values. This is illustrated in Figure 2a, where one can pick any value of β and for that value there is always a range of values of ζ to select to ensure periodic behaviour. For the rate of becoming infectious with varicella, ω, oscillations were possible for all values except when the HZ recovery rate ρ was at the bottom of its range (around 46-48) and δ was approaching 0.017. Similarly, the entire range for the varicella recovery rate γ values allowed for oscillatory solutions except when δ and θ were varied. The HZ recovery rate, ρ, behaved in the same way. Note that β,ω,γ and ρ are considered well informed.
Another group of parameters turned out to be unnecessary to support oscillatory behaviour. One of these parameters was ζ: it usually had to be under 0.02, i.e. the force of varicella infection on those who cleared varicella previously had to be reduced by at least 98%. This means that ζ could as well be zero, as shown in Figure 2a. Parameter δ (the rate of HZ reactivation following recovery from HZ) tended to be under 0.07 unless κ,μ,ζ and σ were manipulated. As shown in Figure 2b, δ can be as large as 0.1 if σ is somewhere between 0.2 and 0.55, but it can also be zero. Similarly, p was consistent with oscillations when small (of the order 10−3) or equal to zero in line with the values estimated in [4,16,31] but much less than p=0.05 used in .
At last, there were six parameters that allowed for oscillatory behaviour when the ranges of values they could vary in were narrowed to some extent. Specifically, the boosting parameter κ had to exceed 2.5. Figure 2c illustrates that in the grey area where the values of κ and θ corresponding to oscillations are, κ is clearly much larger than 1. Parameter μ was restricted to under 0.015 (the corresponding life expectancy of at least 66.7 years) unless η increased and then μ could even be 0.04. The loss of full immunity rate, σ, had to exceed 0.1 (duration of immunity under 10 years) but was usually limited by 0.8 (or 1.1 if κ was approaching 8). The values of θ under about 2 or a little less were not compatible with oscillations (see Figure 2c) and neither was the HZ reactivation rate η under 0.05, which is evident from Figure 2d. At last, α had to be over 0.7 in most cases - a notable exception is shown in Figure 2d where α can approach 0.3 assuming that η is around 0.5.
It is pertinent to note that periodic solutions typically corresponded to repeat varicella incidence that we would call substantial. For example, varying ζ while all other model parameters are fixed results in detection of a Hopf point at ζ=0.01314. This point corresponds to the yearly incidence of repeat varicella periodically reaching a maximum of 1.05% of the modelled population. The yearly HZ incidence is also periodic in this circumstance averaging 0.00215% and peaking at around ±0.0023%.
In this study, we developed a simple VZV transmission model with novel features to explore the influence of boosting and disease recurrence on the potential for periodicity in varicella and zoster disease incidence. We relied upon bifurcation analysis of this model to assess the relationship between key parameters and the existence of periodic solutions.
In comparison with other published VZV transmission models, our model has four additional features: repeated HZ attacks and rates of boosting exceeding the force of infection are allowed, while we also include the possibility of clinical reinfection with varicella from a waned immunity state and in individuals who have experienced zoster.
reinfection with varicella for individuals who recovered from varicella and are in the weakly immune state W (as mentioned in Results, ζ can be zero);
reactivation of HZ following recovery from zoster (δ=0 is possible);
contribution of HZ to varicella force of infection (p can be zero).
duration of full immunity to VZV is no longer than 10 years (σ>0.1);
duration of immunity to VZV following recovery from HZ should not exceed 6 months (θ>2)
immune boosting requires less exposure than primary varicella infection (κ≥2.5);
unless the average HZ reactivation rate η is only slightly under 0.2 (that is 4 times its base value), those who have had HZ and are in contact with individuals infected with varicella have little or no added protection against varicella as compared with individuals never infected with VZV (α>0.7).
It is not immediately clear how to assess plausibility of these requirements. The existing data can only give some indications regarding the average time needed to return from I to I or progress from I to IZ. Since the immunity in question is the full immunity corresponding to the average time needed to progress from I to W and the latter is included to implement the Hope-Simpson hypothesis (i.e. the very existence of this state is not an established fact), we should note that there is no data to suggest that σ>0.1 may be unrealistic. There is also not enough data to clarify the average duration of staying in the RZ compartment.
The exposure required for boosting VZV immunity is not well understood but the Hope-Simpson hypothesis is widely applied in modelling studies with κ≤1. Nonetheless, we are not aware of any strong evidence in support of this restriction which precludes κ≥2.5 (a relevant discussion on why κ>1 is likely to be valid for pertussis can be found in ).
Based on the Hope-Simpson hypothesis, our model makes individuals recovered from HZ immune but we do not know how long this immunity should last on average. Our analysis shows that we can obtain periodicity if this duration is under 6 months which we can not consider implausible until there is evidence to prove otherwise. One would also expect that being in the post-zoster state individuals may be able to strengthen their immunity when exposed to varicella and our results do not contradict this expectation. However, they indicate that the onset of oscillatory behaviour would occur if this protection, if gained at all, were insignificant: α>0.7 means that someone who have had HZ, then was immune for a while and then became susceptible to VZV is becoming reinfected with varicella at an average rate reduced by up to 30% at best or much less, as compared with a rate this individual would be becoming infected with varicella at if he or she was exposed to varicella for the first time. As those who have had HZ are most often elderly and their immunity is weakened by age, there is a reason to believe that they may become susceptible to varicella reinfection via contact with infectious varicella cases and the immune boosting mechanism may not work for these individuals very well. The absence of a strong evidence base for such cases does not necessarily preclude them being common. For example, routine reporting schemes are frequently recording numerous “unspecified VZV” diagnoses  which could be either HZ or varicella. Separating varicella and HZ diagnoses can be difficult. Laboratory testing is infrequent, and the ability to differentiate re-infection from reactivation is therefore based on the disease syndrome rather than confirmation through viral testing. As discussed by Volpi et al.  the clinical diagnosis of HZ is subject to error, with the Shingles Prevention Study reporting that 24% of clinical HZ diagnoses were not confirmed to have HZ by laboratory testing  and a UK study found 17% (48 out of 278) of GP clinical HZ diagnoses to be incorrect . Thus, while reinfection with varicella is suspected to be infrequent, confirmation of this fact is difficult given errors in reporting and diagnosis.
Hence, while other sources of periodicity, including seasonally varying transmissibility may appear more likely explanations of varicella periodicity, outcomes from our model appear to accord with data on recurrence and reinfection, suggesting a role for these processes in enhanced models of VZV transmission.
We must acknowledge substantial limitations in relation to the realism exhibited by our model. It is generally accepted that age-dependent transmission is important in the context of varicella and HZ transmission and disease incidence but we do not include this here. Instead we include only a basic demographic process involving a constant birth/death rate μ as our approach aimed to isolate the impacts of recurrence and reinfection in a prototype model via minimising the total number of model parameters (a similar approach was used in, for example, ). A more realistic model would incorporate age-dependent transmission rates and age-dependent risks of HZ reactivation, with potential to consider more subtle differences such as gender related risks of disease. Other models also implement the assumption that breakthrough varicella is milder and less contagious than primary infection, based on experience with vaccine trials and implementation. Consideration of the impact of vaccination is also ignored here, given that the focus is on qualitative behaviours in a simple model. Note that as a result of the mentioned limitations, varicella incidence in the entire population produced by our model is higher than zoster incidence (see Figure 3). This is typical in the absence of vaccination, while in the post-vaccination period one often observes zoster incidence increasing and exceeding that of varicella (see, for example, [17,38]). The increase in zoster could be further facilitated by the ageing of population as discussed in .
Periodicity of varicella and HZ as such do not seem to be of major consequences to public health. This is partly because periodicity of zoster has yet to be firmly substantiated. For example, an Italian study found an annual cycle in HZ notifications with significant peaks in May-June . On the other hand, HZ periodicity was not detected by a number of population-based studies and long-term analyses of data [24,25]. However, there has been a steady interest in gaining a better understanding of the reasons for VZV reactivation - should that be successful, practical implications are likely to follow. It is possible that seasonal variations in the incidence of immunologically related diseases may be related to circannual variations in immune response  or, in the case of HZ, to the effects of ultraviolet (UV) . We suggest that varicella and zoster periodicity may be potentially explained solely by VZV natural history and including clinical reinfection with varicella and the possibility that only a limited exposure is needed for immune boosting in modelling studies may result in more accurate predictions.
The Kirby Institute receives funding from the Australian Government Department of Health and Ageing. The views expressed in this publication do not necessarily represent the position of the Australian Government.
- Hope-Simpson RE. The nature of herpes zoster: a long-term study and a new hypothesis. Proc R Soc Med. 1965; 58:9–20.PubMed CentralPubMedGoogle Scholar
- Ogunjimi B, Van Damme P, Beutels P. Herpes Zoster Risk Reduction through Exposure to Chickenpox Patients: A Systematic Multidisciplinary Review. PLoS ONE. 2013; 8(6):66485.View ArticleGoogle Scholar
- Thomas SL, Hall AJ. What does epidemiology tell us about risk factors for herpes zoster?Lancet Infect Dis. 2004; 4(1):26–33.View ArticlePubMedGoogle Scholar
- Brisson M, Edmunds WJ, Gay NJ, Law B, De Serres G. Modelling the impact of immunization on the epidemiology of varicella zoster virus. Epidemiol Infect. 2000; 125(3):651–69.View ArticlePubMed CentralPubMedGoogle Scholar
- Hall S, Maupin T, Seward J, Jumaan AO, Peterson C, Goldman G, et al.Second varicella infections: are they more common than previously thought?Pediatrics. 2002; 109(6):1068–73.View ArticlePubMedGoogle Scholar
- Marin M, Watson TL, Chaves SS, Civen R, Watson BM, Zhang JX, et al.Varicella among adults: data from an active surveillance project, 1995–2005. J Infect Dis. 2008; 197(Suppl 2):94–100.View ArticleGoogle Scholar
- Gershon AA, Steinberg SP, Gelb L. Clinical reinfection with varicella-zoster virus. J Infect Dis. 1984; 149(2):137–42.View ArticlePubMedGoogle Scholar
- Junker AK, Angus E, Thomas EE. Recurrent varicella-zoster virus infections in apparently immunocompetent children. Pediatr Infect Dis J. 1991; 10(8):569–75.View ArticlePubMedGoogle Scholar
- Brisson M, Melkonyan G, Drolet M, De Serres G, Thibeault R, De Wals P. Modeling the impact of one- and two-dose varicella vaccination on the epidemiology of varicella and zoster. Vaccine. 2010; 28(19):3385–97.View ArticlePubMedGoogle Scholar
- van Hoek AJ, Melegaro A, Zagheni E, Edmunds WJ, Gay N. Modelling the impact of a combined varicella and zoster vaccination programme on the epidemiology of varicella zoster virus in England. Vaccine. 2011; 29(13):2411–20.View ArticlePubMedGoogle Scholar
- Wilson JB. Thirty one years of herpes zoster in a rural practice. Br Med J (Clinical Research ed.) 1986; 293(6558):1349–51.View ArticleGoogle Scholar
- Yawn B, Wollan P, Kurland M, St Sauver J, Saddier P. Herpes Zoster Recurrences More Frequent Than Previously Reported. Mayo Clinic Proc. 2011; 86(2):88–93.View ArticleGoogle Scholar
- Oxman MN, Levin MJ, Johnson GR, Schmader KE, Straus SE, Gelb LD, et al.Shingles Prevention Study Group, A vaccine to prevent herpes zoster and postherpetic neuralgia in older adults. N Engl J Med. 2005; 352(22):2271–84.View ArticlePubMedGoogle Scholar
- Donahue JG, Kieke BA, Gargiullo PM, Jumaan AO, Berger NR, McCauley JS, et al.Herpes Zoster and Exposure to the Varicella Zoster Virus in an Era of Varicella Vaccination. Am J Pub Health. 2010; 100(6):1116–22.View ArticleGoogle Scholar
- Brisson M, Gay NJ, Edmunds WJ, Andrews NJ. Exposure to varicella boosts immunity to herpes-zoster: implications for mass vaccination against chickenpox. Vaccine. 2002; 20(19–20):2500–7.View ArticlePubMedGoogle Scholar
- Gidding HF, Brisson M, MacIntyre CR, Burgess MA. Modelling the impact of vaccination on the epidemiology of varicella zoster virus in Australia. Austr N Z J Pub Health. 2005; 29(6):544–51.View ArticleGoogle Scholar
- Poletti P, Melegaro A, Ajelli M, del Fava E, Guzzetta G, Faustini L, et al.Perspectives on the Impact of Varicella Immunization on Herpes Zoster. A Model-Based Evaluation from Three European Countries. PLoS ONE. 2013; 8(4):60732.View ArticleGoogle Scholar
- Lavine JS, King AA, Bjørnstad ON. Natural immune boosting in pertussis dynamics and the potential for long-term vaccine failure. Proc Nat Acade Sci USA. 2011; 108(17):7259–64.View ArticleGoogle Scholar
- Janeway CJ, Travers P, Walport M. Immunobiology: The Immune System in Health and Disease, 5th edn. New York: Garland Science; 2001.Google Scholar
- Zinkernagel RM. Regulation of the immune response by Antigen. Science. 2001; 293(5528):251–3.View ArticlePubMedGoogle Scholar
- Heininger U, Seward JF. Varicella. The Lancet. 2006; 368(9544):1365–76.View ArticleGoogle Scholar
- Gallerani M, Manfredini R. Seasonal variation in herpes zoster infection. Br J Dermatol. 2000; 142(3):588–9.View ArticlePubMedGoogle Scholar
- Toyama N, Shiraki K, and Members of the Society of the Miyazaki Prefecture Dermatologists. Epidemiology of herpes zoster and its relationship to varicella in Japan: A 10-year survey of 48,388 herpes zoster cases in Miyazaki prefecture. J Med Virol. 2009; 81(12):2053–8.View ArticlePubMedGoogle Scholar
- Chant KG, Sullivan EA, Burgess MA, Ferson MJ, Forrest JM, Baird LM, et al.Varicella-zoster virus infection in australia. Aust N Z J Pub Health. 1998; 22(4):413–8.View ArticleGoogle Scholar
- Chapman RS, Cross KW, Fleming DM. The incidence of shingles and its implications for vaccination policy. Vaccine. 2003; 21(19–20):2541–7.View ArticlePubMedGoogle Scholar
- Olsen LF, Schaffer WM. Chaos versus noisy periodicity: alternative hypotheses for childhood epidemics. Science. 1990; 249(4968):499–504.View ArticlePubMedGoogle Scholar
- Nardone A, de Ory F, Carton M, Cohen D, van Damme P, Davidkin I, et al.The comparative sero-epidemiology of varicella zoster virus in 11 countries in the European region. Vaccine. 2007; 25(45):7866–72.View ArticlePubMedGoogle Scholar
- Diekmann O, Heesterbeek JAP. Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation. New York: Wiley; 2000.Google Scholar
- Australian Bureau of Statistics. 4102.0 – Life expectancy trends. 2011. www.abs.gov.au Accessed 29 March 2015.
- Guzzetta G, Poletti P, Del Fava E, Ajelli M, Scalia Tomba GP, et al.Hope-Simpson’s Progressive Immunity Hypothesis as a Possible Explanation for Herpes Zoster Incidence Data. Am J Epidemiol. 2013; 177(10):1134–42.View ArticlePubMedGoogle Scholar
- Karhunen M, Leino T, Salo H, Davidkin I, Kilpi T, Auranen K. Modelling the impact of varicella vaccination on varicella and zoster. Epidemiol Infect. 2009; 138(04):469.View ArticlePubMedGoogle Scholar
- Ermentrout B. Simulating, Analyzing, and Animating Dynamical Systems: a Guide to XPPAUT for Researchers and Students vol. 14. Philadelphia: Siam; 2002.View ArticleGoogle Scholar
- Kuznetsov Y. Elements of Applied Bifurcation Theory, 2nd edn. New York: Springer; 1998.Google Scholar
- Seydel R. Practical Bifurcation and Stability Analysis, 3rd edn. New York: Springer; 2009.Google Scholar
- Australian Government. Department of Health. National Notifiable Diseases Surveillance System - Number of notifications of Varicella zoster (Unspecified). www.health.gov.au. Accessed 29 March 2015.
- Volpi A, Gatti A, Pica F. Frequency of herpes zoster recurrence. Mayo Clinic Proc. 2011; 86(6):586–5867.View ArticleGoogle Scholar
- Scott FT, Leedham-Green ME, Barrett-Muir WY, Hawrami K, Gallagher WJ, Johnson R, et al.A study of shingles and the development of postherpetic neuralgia in East London. J Med Virol. 2003; 70(S1):24–30.View ArticleGoogle Scholar
- Wu P-Y, Wu H-DI, Chou T-C, Sung F-C. Varicella Vaccination Alters the Chronological Trends of Herpes Zoster and Varicella. PLoS ONE. 2013; 8(10):77709.View ArticleGoogle Scholar
- Levi FA, Canon C, Touitou Y, Reinberg A, Mathe G. Seasonal modulation of the circadian time structure of circulating T and natural killer lymphocyte subsets from healthy subjects. J Clin Investig. 1988; 81(2):407.View ArticlePubMed CentralPubMedGoogle Scholar
- Zak-Prelich M, Borkowski JL, Alexander F, Norval M. The role of solar ultraviolet irradiation in zoster. Epidemiol Infect. 2003; 129(03):593.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.