Laser en Autismo
Laser en Autismo
Laser en Autismo
Clinical
Medicine
Research
Advances in Experimental Medicine
and Biology
Clinical and Experimental Biomedicine
Volume 1116
Subseries Editor
Mieczyslaw Pokorski
More information about this series at http://www.springer.com/series/16003
Mieczyslaw Pokorski
Editor
Clinical Medicine
Research
Editor
Mieczyslaw Pokorski
Opole Medical School
Opole, Poland
This Springer imprint is published by the registered company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Contents
v
vi Contents
Bioprogressive Paradigm
in Physiotherapeutic and Antiaging
Strategies: A Review
1
2 M. Pokorski et al.
avoid the risk of falling or to begin a rehabilitation frameworks often occur in adulthood, and
protocol (Israely et al. 2018; Zur et al. 2018; behaviors are generally modifiable through a
Aalen et al. 2014; Arbeev et al. 2009). series of rehabilitative interventions (Atkins
A bioprogressive multisensory rehabilitation et al. 2014; Barzilai et al. 2012; Bergman et al.
approach has proven the most suitable to face a 2007). Thus, nutrition, frequency of physical
physiological reduction in stability. The activity, and other crucial factors need to be
approach, supported by modern technology, examined to develop the appropriated lifestyle,
enables the creation of a changing training envi- which would also provide succor to metabolic
ronment, with the aim to increase proprioceptive activity, cerebral functionality, and postural con-
perception and balance reaction (Bellomo et al. trol (Si and Liu 2014; Bosma-den Boer et al.
2009). Improvements can be monitored in a mod- 2012; Misirli et al. 2012; Agnoli et al. 2011; Hu
ern gait analysis laboratory to describe variations 2002).
in balance reactions, which enables the therapist Studying the inter-system correlations and the
to periodically adapt the training program to rules that underlie the perceptive-motor pattern
changes in the patient condition to reach increas- holds a central role in the understanding of
ingly important therapeutic goals (Basile et al. relationships between movement, function, envi-
2016). Parkinson’s patients treated with a ronment context, and psycho-behavioral
bioprogressive rehabilitation paradigm, consisting characteristics of subjects. The human being
of a combination of dynamic antigravity postural forms a unique biosystem whose management
system (SPAD), auditory cue system (metro- requires a holistic approach, personalized for spe-
nome), high-intensity vibratory system (VISS), cific patient’s characteristics (Saggini et al. 2012).
and active stretching based on virtual reality reha- Posture is identified as the morphological, senso-
bilitation system (VRRS) for several months, rial, and motor expression of the evolved identity.
experience a distinct improvement in balance sta- The somatosensory system, visual and auditory
bility and gait pattern, with a reduction in falls structures, and vestibular receptors are essential
and overall enhancement of quality of life for maintaining the brain-regulated body motion
(Bellomo et al. 2014). Modern assistance and balance, as well as a feel of one’s comfort
methods enable the elderly to access the training (Zimmermann et al. 2013). A potential alteration
facilities at any time and to use them in a safe of sensory afferentation can be compensated
environment while continuing to provide through the redundancy characteristics of a sys-
instructions for rehabilitation and monitoring of tem, which also respond to appropriate
progress (Akushevich et al. 2013). Some innova- myofascial adaptations (Reghem et al. 2014).
tive examples for the preventive and rehabilitative The aging process underlies a progressive deteri-
intervention, which illustrate this increasingly oration of the delicate mechanisms of neurologi-
diversified scope and ever-changing future per- cal osteomuscular-fascial optimization,
spective of the intervention program, have been manifesting as an altered “afferent stimulation”
described by Young (1983). that determines “postural dysfunction” (Barassi
Aging concerns physical, cognitive, and social et al. 2018a, b, c, d).
spheres of individuals. The cardiovascular sys- Already at gestational age week 5, a trilaminar
tem, muscle-skeletal mass, osteoarthritic struc- embryo may be recognized, with the interposition
ture, immuno-humoral metabolism, and central of mesoderm between ectoderm and endoderm.
nervous system disorders predominate in the The epithelial to mesenchymal cell transition
aged, altering the organism’s functionality emerges, during which the structure assumes its
(Israely et al. 2018; Carmeli et al. 2016). Diabetes function and form. In this mechanism, muscular
and cardiovascular diseases, closely linked to chains and fascial continuity appear to be the
inflammation and oxidative stress, belong to at evolutionary expressions of the spinal cord
least partly modifiable risk factors, with a proper (Barassi et al. 2018b). Paraxial mesodermal
diet and correct lifestyle. Such pathologic elements of the skeletal muscular system are
Bioprogressive Paradigm in Physiotherapeutic and Antiaging Strategies: A Review 3
formed as the splanchnic mesodermal derivations, 2018). There appears a plausibility that such
in particular smooth muscles, derived from the structures play a central role in the body’s self-
occipital-to-sacral region somites and head-region management of edema to recreate homeostatic
somitomeres (Saggini et al. 1996). harmony (Johnson 1977). Thus, contrarily to the
It is known that fascial connective tissue is old physiological tenet of functional neutrality of
able to pick up nociceptive and proprioceptive connective tissue, purportedly serving the crea-
information thanks to the presence of mechano- tion of a parenchymal scaffold for cell arrange-
and nociceptors and its spatial organization (Van ment, connective cells appear to host a multitude
der Wal 2009). Studies also demonstrate that of active functions, inter alia, inter-cell communi-
connective tissue creates a network between the cation, proprioceptive signal transmission, regu-
musculoskeletal, neural, endocrine, circulatory, lation of extracellular space fluid, and others
respiratory, and immune systems, creating a (Bordoni and Zanier 2015). Nowadays, these
milieu for inter-tissue cooperation (Vecchiet cells appear to underlie the effectiveness of the
et al. 1999). In particular, fascia is composed of bioprogressive rehabilitation paradigms that are
cells that have the ability to contract and commu- set to affect their function to restore a proper body
nicate with one another, ensuring the supervision neuro-myo-sensory integration, as schematically
of body activity, such as the transmission of ten- depicted in Fig. 1.
sion created by muscles or management of inter-
stitial fluids (Bordoni and Zanier 2015; Schleip
2003). One of the latest studies performed on 2 Stretch Exercise
fascia tissue demonstrates, using a confocal laser
endomicroscopy, the presence of “reticular One of the most common exercise techniques is
structures” in visceral submucosa and dermal muscular-fascial stretch. Stretching technique
periarterial stroma, which are hypothesized to be modifies the Golgi apparatus activity which, via
an extension of the interstitial space (Benias et al. the spinal cord, decreases the firing rate of alpha
BODY DYNAMICS
MESODERM
(Connective Tissue)
Microgravity Environment
Tissue Remodelling Rehabilitation Rehabilitation
Action
GOLGI OSTEO-MYO-FASCIAL
APPARATUS SYSTEM
STRETCHING MAV
Mechano-Acoustic Vibration Terrestrial Acquatic
BODY DYNAMICS
OPTIMIZATION
Fig. 1 Schematic diagram of the bioprogressive philosophy of physiotherapeutic interventions, pointing to the central
role of connective tissue in maintaining the neuro-myo-sensory integration
4 M. Pokorski et al.
motor neurons, and, in turn, muscle tone, for held in that position, while rehabilitant contracts
instance, as occurs in Yoga posture (Barassi the muscle isometrically to its maximum for a
et al. 2018c; Schleip 2003). There exist various certain amount of time, followed by a brief
types of stretching technique. The unresolved phase of passive stretching. CRAC is an exten-
issue with physical stretch is what it should be sion of the former technique in which after the
like in the setting of health-care practitioners, maximum isometric contraction of the agonistic
particularly those engaged in sports training or muscle, rehabilitant contracts the corresponding
rehabilitative framework in the elderly (Caplan antagonistic muscle for a certain amount of time.
et al. 2009). There are three types of stretch exer- Ultimately, it seems that a mix of stretch modes,
cise: static stretching, dynamic stretching, and in particular static stretch followed by dynamic
proprioceptive neuromuscular facilitation (PNF) stretch, optimizes the benefits gained from physi-
(Fig. 1). A combination of these approaches cal rehabilitation, simultaneously reducing the
guarantees autogenic inhibition, mutual inhibi- number of exercise-induced muscle impairments
tion, stress relaxation, and a reduction of pain or injuries (Rowlands et al. 2003). To this end,
consequently to the gate control theory, all of meditation techniques consisting of a variable
which leads to improvements in the range of motion ingredient and resembling the CRAC
motion of joints (Nelson et al. 2005; Wicke sequence of stretches, such as is Tai Chi and the
et al. 2014; Hindle et al. 2012). Static stretching like, advance in popularity of late (Hindle et al.
consists of holding muscles in an elongated posi- 2012). An additional benefit of meditation is the
tion for an extended period. Contrarily, dynamic improvement in mood, emotional stability, and
stretch consists of moving joints through their cognition (Pokorski and Suchorzynska 2018;
range of motion. There are pros and cons of Lam et al. 2014).
both types of stretch. Static isometric exercise,
although being physically motionless, is a strong
respiratory stimulant that provides succor to arte- 3 Aging and Rehabilitation
rial blood and, consequently, tissue oxygenation.
It also increases heart rate and blood pressure Old age should not be seen as a disease or pathol-
(Pokorski et al. 1990). These effects, thought of ogy but as a stage of life that can be fully explored
as typical and beneficial from the standpoint of to pursue personal and social goals. Human func-
physical rehabilitation or training, may jeopardize tionality has an outstanding plasticity and is capa-
the frail elderly due to abruptly increased meta- ble of compensating for the lost structural
bolic rate. Static stretch may also lead to muscle elements caused by the age-related dampening
weakness, which hampers subsequent physical of cellular metabolism. That may be exemplified,
performance. Dynamic stretch, on the other for instance, by the adaptive hyperventilatory
hand, gives proneness to injuries, which may responses to strenuous stimuli, such as exercise,
actually hamper the range of motion that already in old age, despite a severe structural deterioration
declines with age but is important for daily life of the lungs and airways, limiting the alveolar gas
activities (Konrad and Tilp 2014). Dynamic exchange area and thus oxygen diffusion
stretch, however, increases the nerve conduction (Pokorski and Marczak 2003). Likewise, cogni-
velocity, muscle compliance, enzymatic turnover, tive brain functions are maintained due to a
and the central control running down from the greater neuronal synchronization and a flexible
brain to muscles (Wicke et al. 2014). The PNF use of cognitive reserve or recruitment of alterna-
combines the stretching and contraction of tive neuronal pathways, normally constituting a
reciprocally related muscle groups in a cyclical dormant reserve. The antiaging strategy is a pro-
pattern. “Contract-relax” (CR) and “contract- cess that not only involves combating the disease
relax-antagonist-contract” (CRAC) methods are and its effects but also a large psychological com-
the representative techniques of PNF. In the CR ponent that is closely linked to behavior, person-
technique, the target muscle is lengthened and ality traits, and other factors (Angel and Hogan
Bioprogressive Paradigm in Physiotherapeutic and Antiaging Strategies: A Review 5
1992). Exercise is an antiaging strategy linked to regulation of genes that are responsible for
physical and cognitive rehabilitation, develop- expression of sarcomeric proteins and for recruit-
ment of positive stress, and lifestyle changes ment of satellite cells. Integrated treatment with
(Rebelo-Marques et al. 2018). The mind-to-body focused mechano-acoustic vibrations has a bene-
connection is indispensable for a motivated ficial effect on bone mineral density (BMD) and
engagement in antiaging behaviors, such as exer- T-score and on muscle strength and quality of life
cise, and in physical rehabilitation in case of a of osteoporotic subjects (Kwak et al. 2017;
disease-induced handicap (Wrycza and Baudisch Saggini et al. 2017).
2012; Depp and Jeste 2009). Exercise has a pro- The body’s agility to react to external forces
tective effect on cognitive functions and helps plays a central role in rehabilitation. Acting on
combat the deleterious effects of stress and this mechanism enables the most individual adap-
aging. The practice of proper physical exercise tation of a therapeutic protocol and thus the use of
can improve physiological and functional the entirety of the rehabilitation power and also
responses, notably cardiovascular fitness and the exposition of compensative reactions and
respiratory function, and thus physical health rebalancing of human body network (Fig. 1). To
and quality of life of the elderly (Bellomo et al. this end, microgravity environment is particularly
2017). Exercise also plays a key role in the assess- effective in reversing functional detriments and
ment and prevention of frailty in the elderly limits concerning daily activities in patients with
(Singh 2002). In adulthood, aerobic skills and muscle and balance disorders (Kaneda et al.
general physical activity are associated with 2008). The exemplary is the aquatic microgravity
lower mortality and morbidity from cardiovascu- environment, produced by immersion in water,
lar diseases. Benefits of physical activity also which improves body stability due to a decrease
include reductions in a risk for stroke, diabetes in neural system activation caused by a central
mellitus, cancer, and osteoporosis (Ruiz et al. redistribution of the blood volume (Barassi et al.
2008; Lee and Skerrett 2001; Bijnen et al. 1999; 2018a). A decrease in sympathetic activity,
Hakim et al. 1998; Paffenbarger et al. 1993). accompanied by increased parasympathetic activ-
Aging is inherently involved with a loss of ity, promoting vasodilation and blood circulation,
muscle mass and function leading to sarcopenia increases blood flow, accelerates cellular metabo-
and, consequently, to postural dysfunction lism and removal of unwanted materials, and
(Carmeli 2017; Lin and Bhattacharyya 2012). decreases pain sensitivity (Forestier and Françon
Sarcopenia is underlain by a reduction in protein 2008). The autonomic nervous system
synthesis and an increase in muscle protein deg- remodeling positively affects the perception pain
radation (Delli Pizzi et al. 2017). Further, the and fatigue associated with musculoskeletal
ability of muscle regeneration is severely disorders due mostly to decreasing muscle ten-
compromised by aging, which can lead to disabil- sion (Barassi et al. 2018d; Bellomo et al. 2012;
ity, particularly in patients with concomitant Yasui et al. 2010). The rebalancing of gait pattern
diseases or organ failures. Skeletal muscle dys- involves the development of relatively stable
function can be ameliorated by well-suited exer- changes in spinal and supraspinal networks that
cise that improves the function of muscle are modeled by proprioceptive and somatosen-
mitochondria (Barbieri et al. 2015). Beneficial sory sensory information specific for the task.
effects of focal vibrations on the well-being of The development of specialized training
ailing elderly patients are beyond doubt when protocols for assisted exercise in the microgravity
analyzing the results of studies on proprioception, environment hinges on the realization of the
muscle tone, and quality of gait (Toosizadeh et al. importance of locomotor learning and how the
2018; Saggini and Bellomo 2015). Vibrations of use of specific devices alters the demand stem-
300 Hz effectively increase muscle strength, act- ming from motor tasks (Wolpert et al. 2001). The
ing on neuromuscular receptors, Golgi tendon literature unambiguously shows that a specific
organs, Pacini corpuscles, as well as through the reeducation protocol in relieving body weight
6 M. Pokorski et al.
can improve the static and dynamic balance, lead- of biology and function into the realm of rehabil-
ing to a better overall motor performance itation therapy (Fig. 1). It enables the optimiza-
(Moreno et al. 2013; Threlkeld et al. 2003). In tion of redundant biological circuits and
our daily practice, we implement the single- arrangements through personalized therapeutic
photon avalanche detector (SPAD), a photon interventions rather than application of rigidly
detecting system for optical communications. fixed therapy schemes (Trieschmann 1987). Fur-
SPAD is a device able to compose a microgravity ther research should scrutinize the physiological
environment during gait motion, which optimizes underpinnings of the bioprogressive approach to
the body tridimensional alignment in the body better understand its function and potential
weight support modality. It has a dual action: benefits to be gained in physiotherapy.
mechanical, which allows for a cortical-
subcortical neuromotor retraining aimed at the Ethical Approval This review article has been written in
reacquisition of a body scheme in three planes compliance with the ethical standards of the institutional
and/or national research committees and with the 1964
of the space, and functional, which minimizes the
Helsinki declaration and its later amendments. The article
energy consumption and increases the proprio- does not contain any studies with human participants or
ceptive signal inflow to cortical areas during gait animals performed by any of the authors.
(Delli Pizzi et al. 2017; Saggini et al. 2004). Thus,
SPAD improves postural stability through an
increase in exteroceptive and proprioceptive sen-
References
sitivity. The technique is also useful in therapy of
sarcopenic patients for muscle strengthening and Aalen OO, Valberg M, Grotmol T, Tretli S (2014) Under-
pain resolution, along with a postural and gait standing variation in disease risk: the elusive concept
rehabilitation. This innovative therapeutic of frailty. Int J Epidemiol 44(4):1408–1421
approach combines a motor and sensory feedback Agnoli C, Krogh V, Grioni S, Sieri S, Palli D, Masala G,
Sacerdote C, Vineis P, Tumino R, Frasca G, Pala V,
task with creation of feed-forward stimulation in Berrino F, Chiodini P, Mattiello A, Panico S (2011) A
the microgravity environment (Di Pancrazio et al. priori–defined dietary patterns are associated with
2013). reduced risk of stroke in a large Italian cohort. J Nutr
In elderly patients, clinical relevance of 141(8):1552–1558
Akushevich I, Kravchenko J, Ukraintseva S, Arbeev K,
exercise-based rehabilitation and nutritional Kulminski A, Culminskaya IV, Yashin AI (2013) Mor-
interventions, such as increased ingestion of bidity risks among older adults with pre-existing
essential nutrients, results in the long-term age-related diseases. Exp Gerontol 48(12):1395–1401
improvements in strength and muscle mass and Angel JL, Hogan DP (1992) The demography of minority
aging populations. J Fam Hist 17(1):95–115
in functional capacity and decreases a risk of Arbeev KG, Akushevich I, Kulminski AM, Arbeeva LS,
chronic metabolic disorders (Snijders et al. Akushevich L, Ukraintseva SV, Yashin AI (2009)
2009). This integrated, personalized therapeutic Genetic model for longitudinal studies of aging, health,
and precautionary approach to health and disease, and longevity and its potential application to incom-
plete data. J Theor Biol 258(1):103–111
with characteristics of adaptability and biological Atkins JL, Whincup PH, Morris RW, Lennon LT,
progressiveness, improves quality of life in adult- Papacosta O, Wannamethee SG (2014) Sarcopenic
hood. This is a holistic strategy, contained within obesity and risk of cardiovascular disease and mortal-
the biopsychosocial model, developed by Engels ity: a population-based cohort study of older men. J
Am Geriatr Soc 62(2):253–260
(1977), which emphasizes, but is not limited to, Barassi G, Bellomo RG, Porreca A, Giannuzzo G, Irace G,
biological determinants for maintaining health, Trivisano L, Saggini R (2018a) Rehabilitation of
and does not overlook the overall importance of neuromotor disabilities in aquatic microgravity envi-
psychological or socioeconomic factors. All of ronment. Adv Exp Med Biol. https://doi.org/10.1007/
5584_2018_164
them may integrate to act as physical traumas to Barassi G, Bellomo RG, Di Giulio C, Giannuzzo G,
the body system and compromise the homeostatic Irace G, Barbato C, Saggini R (2018b) Effects of
harmony. The neuro-myo-sensory integration of manual somatic stimulation on the autonomic nervous
the bioprogressive philosophy combines aspects system and posture. Adv Exp Med Biol 1070:97–109
Bioprogressive Paradigm in Physiotherapeutic and Antiaging Strategies: A Review 7
Barassi G, Bellomo RG, Di Iulio A, Lococo A, Porreca A, Bosma-den Boer MM, van Wetten ML, Pruimboom L
Di Felice PA, Saggini R (2018c) Preoperative rehabili- (2012) Chronic inflammatory diseases are stimulated
tation in lung cancer patients: yoga approach. Adv Exp by current lifestyle: how diet, stress levels and medica-
Med Biol 1096:19–29 tion prevent our body from recovering. Nutr Metab 9
Barassi G, Bellomo RG, Porreca A, Di Felice PA, (1):32
Prosperi L, Saggini R (2018d) Somato-visceral effects Caplan N, Rogers R, Parr MK, Hayes PR (2009) The
in the treatment of dysmenorrhea: neuromuscular man- effect of proprioceptive neuromuscular facilitation
ual therapy and standard pharmacological treatment. J and static stretch training on running mechanics. J
Altern Complement Med 24(3):291–299 Strength Cond Res 23(4):1175–1180
Barbieri E, Agostini D, Polidori E, Potenza L, Guescini M, Carmeli E (2017) Frailty and primary sarcopenia: a review.
Lucertini F, Annibalini G, Stocchi L, De Santi M, Adv Exp Med Biol 1020:53–68
Stocchi V (2015) The pleiotropic effect of physical Carmeli E, Bachar A, Rom O, Aizenbud D (2016) Oxida-
exercise on mitochondrial dynamics in aging skeletal tive stress and nitric oxide in sedentary older adults
muscle. Oxidative Med Cell Longev 2015:917085 with intellectual and developmental disabilities. Adv
Barzilai N, Huffman DM, Muzumdar RH, Bartke A Exp Med Biol 884:21–27
(2012) The critical role of metabolic pathways in Christensen K, Doblhammer G, Rau R, Vaupel JW (2009)
aging. Diabetes 61(6):1315–1322 Ageing populations: the challenges ahead. Lancet 374
Basile C, Magnani S, Cifani AE, Galasso P, Zucchi S, (9696):1196–1208
Bray A, Saggini R (2016) From gait to walk analysis: Delli Pizzi S, Bellomo RG, Carmignano SM, Ancona E,
the walk analysis lab. In: V Congresso Gruppo Franciotti R, Supplizi M, Barassi G, Onofrj M,
Nazionale Di Bioingegneria (GNB 2016), Naples, Bonanni L, Saggini R (2017) Rehabilitation program
Italy, pp 246–247 based on sensorimotor recovery improves the static
Bellomo RG, Iodice P, Savoia V, Saggini A, Vermiglio G, and dynamic balance and modifies the basal ganglia
Saggini R (2009) Balance and posture in the elderly: an neurochemistry. Medicine 96(50):e8732–e8732
analysis of a sensorimotor rehabilitation protocol. Int J Depp CA, Jeste DV (2009) Definitions and predictors of
Immunopathol Pharmacol 22(3 Suppl):37–44 successful aging: a comprehensive review of larger
Bellomo RG, Barassi G, Iodice P, Di Pancrazio L, quantitative studies. Focus 7(1):137–150
Megna M, Saggini R (2012) Visual sensory disability: Di Pancrazio L, Bellomo RG, Franciotti R, Iodice P,
rehabilitative treatment in an aquatic environment. Int J Galati V, D’Andreagiovanni A, Bifolchetti S,
Immunopathol Pharmacol 25(1 Suppl):17–22 Thomas A, Onofrj M, Bonanni L, Saggini R (2013)
Bellomo RG, Di Pancrazio L, Khodor H, Saggini R, Combined rehabilitation program for postural instabil-
Barassi G, Carmignano SM, Di Stefano A (2014) ity in progressive supranuclear palsy. NeuroReh-
Global bioprogressive rehabilitation program and pos- abilitation 32(4):855–860
tural instability in Parkinson’s disease. Eur Sci J Engel GL (1977) The need for a new medical model: a
2:310–321 challenge for biomedicine. Science 196
Bellomo RG, Tripodi D, D’Ercole S, Pecoraro I, Bosna C, (4286):129–136
Barassi G, Melle G, Porreca A, Saggini R (2017) Forestier R, Françon A (2008) Crenobalneotherapy for
Quality of life: effects of physical activity in an anthro- limb osteoarthritis: systematic literature review and
pometric, cognitive and psychosocial background, and methodological analysis. Joint Bone Spine 75
variation of odontoiatric parameters. Sci Philos 5 (2):138–148
(2):3–18 Hakim AA, Petrovitch H, Burchfiel CM, Ross GW,
Benias PC, Wells RG, Sackey-Aboagye B, Klavan H, Rodriguez BL, White LR, Yano K, Curb JD, Abbott
Reidy J, Buonocore D, Miranda M, Kornacki S, RD (1998) Effects of walking on mortality among
Wayne M, Carr-Locke DL, Theise ND (2018) Struc- nonsmoking retired men. N Engl J Med 338(2):94–99
ture and distribution of an unrecognized interstitium in Hindle K, Whitcomb T, Briggs W, Hong J (2012) Propri-
human tissues. Sci Rep 8(1):4947 oceptive neuromuscular facilitation (PNF): its
Bergman A, Atzmon G, Ye K, MacCarthy T, Barzilai N mechanisms and effects on range of motion and mus-
(2007) Buffering mechanisms in aging: a systems cular function. J Hum Kinet 31:105–113
approach toward uncovering the genetic component Hu FB (2002) Dietary pattern analysis: a new direction
of aging. PLoS Comput Biol 3(8):e170 in nutritional epidemiology. Curr Opin Lipidol
Bijnen FC, Feskens EJ, Caspersen CJ, Nagelkerke N, 13(1):3–9
Mosterd WL, Kromhout D (1999) Baseline and previ- Israely S, Leisman G, Carmeli E (2018) Improvement in
ous physical activity in relation to mortality in elderly hand trajectory of reaching movements by error-
men: the Zutphen Elderly Study. Am J Epidemiol 150 augmentation. Adv Exp Med Biol 1070:71–84
(12):1289–1296 Johnson WC (1977) Case of suspected connective tissue
Bordoni B, Zanier E (2015) Understanding fibroblasts in disease presenting as chronic edema or swelling of
order to comprehend the osteopathic treatment of the face, neck, and upper trunk. Cutis 20(5):670–674
fascia. Evid Based Complement Alternat Med Kaneda K, Sato D, Wakabayashi H, Hanai A, Nomura T
2015:860934 (2008) A comparison of the effects of different water
8 M. Pokorski et al.
exercise programs on balance ability in elderly people. Ruiz JR, Sui X, Lobelo F, Morrow JR, Jackson AW,
J Aging Phys Act 16(4):381–392 Sjöström M, Blair SN (2008) Association between
Konrad A, Tilp M (2014) Increased range of motion after muscular strength and mortality in men: prospective
static stretching is not due to changes in muscle and cohort study. BMJ 337:a439
tendon structures. Clin Biomech 29(6):636–642 Saggini R, Bellomo RG (2015) Integration to focal vibra-
Kwak K, Kim H, Kim D (2017) A study on the improve- tion in neurorehabilitation. Eur J Phys Rehabil Med 51
ment of walking characteristics of the elderly with (4):508
vibration stimuli applied to the tibialis anterior tendon. Saggini R, Giamberardino MA, Gatteschi L, Vecchiet L
Biomed Res Int 2017:5342485 (1996) Myofascial pain syndrome of the peroneus
Lam LC, Chan WM, Kwok TC, Chiu HF (2014) Effec- longus: biomechanical approach. Clin J Pain 12
tiveness of Tai Chi in maintenance of cognitive and (1):30–37
functional abilities in mild cognitive impairment: a Saggini R, Cancelli F, Di Bonaventura V, Bellomo RG
randomised controlled trial. Hong Kong Med J 20 (2004) Efficacy of two micro-gravitational protocols to
(3 Suppl 3):20–23 treat chronic low back pain associated with discal
Lee IM, Skerrett PJ (2001) Physical activity and all-cause lesions: a randomized controlled trial. Eur J Phys
mortality: what is the dose-response relation? Med Sci Rehabil Med 40(4):311
Sports Exerc 33(6 Suppl):S459–S471 Saggini R, Iodice P, Saggini A, Scuderi N, Bellomo RG
Lin HW, Bhattacharyya N (2012) Balance disorders in the (2012) Biomechanical modification and sense motor
elderly: epidemiology and functional impact. Laryngo- control of body posture after plastic surgery. Int J
scope 122(8):1858–1861 Immunopathol Pharmacol 25:1 Suppl):1–1 Suppl):8
Misirli G, Benetou V, Lagiou P, Bamia C, Saggini R, Ancona E, Carmignano SM, Supplizi M,
Trichopoulos D, Trichopoulou A (2012) Relation of Barassi G, Bellomo RG (2017) Effect of combined
the traditional Mediterranean diet to cerebrovascular treatment with focused mechano-acoustic vibration
disease in a Mediterranean population. Am J and pharmacological therapy on bone mineral density
Epidemiol 176(12):1185–1192 and muscle strength in post-menopausal women. Clin
Moreno JC, Barroso F, Farina D, Gizzi L, Santos C, Cases Miner Bone Metab 14(3):305
Molinari M, Pons JL (2013) Effects of robotic guid- Schleip R (2003) Fascial plasticity–a new neuro-
ance on the coordination of locomotion. J Neuroeng biological explanation part 2. J Bodyw Mov Ther
Rehabil 10(1):79 7(2):104–116
Nelson AG, Kokkonen J, Arnall DA (2005) Acute muscle Si H, Liu D (2014) Dietary antiaging phytochemicals and
stretching inhibits muscle strength endurance perfor- mechanisms associated with prolonged survival. J Nutr
mance. J Strength Cond Res 19(2):338 Biochem 25(6):581–591
Paffenbarger RS Jr, Hyde RT, Wing AL, Lee IM, Jung Singh MA (2002) Exercise comes of age: rationale and
DL, Kampert JB (1993) The association of changes in recommendations for a geriatric exercise prescription.
physical-activity level and other lifestyle J Gerontol A Biol Sci Med Sci 57(5):M262–M282
characteristics with mortality among men. N Engl J Snijders T, Verdijk LB, van Loon LJ (2009) The impact of
Med 328(8):538–545 sarcopenia and exercise training on skeletal muscle
Pokorski M, Marczak M (2003) Ventilatory response to satellite cells. Ageing Res Rev 8(4):328–338
hypoxia in elderly women. Ann Hum Biol 30:53–64 Threlkeld AJ, Cooper LD, Monger BP, Craven AN,
Pokorski M, Suchorzynska A (2018) Psychobehavioral Haupt HG (2003) Temporospatial and kinematic gait
effects of meditation. Adv Exp Med Biol 1023:85–91 alterations during treadmill walking with body weight
Pokorski M, Masuda A, Paulev PE, Sakakibara Y, Ahn B, suspension. Gait Posture 17(3):235–245
Takaishi S, Nishibayashi Y, Honda Y (1990) Ventila- Toosizadeh N, Ehsani H, Miramonte M, Mohler J (2018)
tory and cardiovascular responses to hypoxic and Proprioceptive impairments in high fall risk older
hyperoxic static handgrip exercise in man. Respir adults: the effect of mechanical calf vibration on pos-
Physiol 81:189 202 tural balance. Biomed Eng Online 17(1):51
Rebelo-Marques A, De Sousa Lages A, Andrade R, Fontes Trieschmann RB (1987) Aging with a disability. Demos
Ribeiro C, Mota Pinto A, Carrilho F, Espregueira- Publications, New York, 148 pp
Mendes J (2018) Aging hallmarks: the benefits of Van der Wal J (2009) The architecture of the connective
physical exercise. Front Endocrinol 9:258 tissue in the musculoskeletal system—an often
Reghem E, Chèze L, Coppens Y, Pouydebat E (2014) The overlooked functional parameter as to proprioception
influence of body posture on the kinematics of prehen- in the locomotor apparatus. Int J Ther Massage Body-
sion in humans and gorillas (Gorilla gorilla). Exp Brain work 2(4):9
Res 232(3):1047–1056 Vecchiet L, Vecchiet J, Giamberardino MA (1999)
Rowlands AV, Marginson VF, Lee J (2003) Chronic flex- Referred muscle pain: clinical and pathophysiologic
ibility gains: effect of isometric contraction duration aspects. Curr Rev Pain 3(6):489–498
during proprioceptive neuromuscular facilitation Wicke J, Gainey K, Figueroa M (2014) A comparison
stretching techniques. Res Q Exerc Sport 74(1):47–51 of self-administered proprioceptive neuromuscular
Bioprogressive Paradigm in Physiotherapeutic and Antiaging Strategies: A Review 9
facilitation to static stretching on range of motion and hemodynamics during thermotherapy on the neck.
flexibility. J Strength Cond Res 28(1):168–172 Auton Neurosci 156(1–2):96–103
Wolpert DM, Ghahramani Z, Flanagan JR (2001) Young EA (1983) Nutrition, aging, and the aged. Med
Perspectives and problems in motor learning. Trends Clin North Am 67(2):295–313
Cogn Sci 5(11):487–494 Zimmermann M, Toni I, de Lange FP (2013) Body posture
Wrycza TF, Baudisch A (2012) How life expectancy modulates action perception. J Neurosci 33
varies with perturbations in age-specific mortality. (14):5930–5938
Demogr Res 27:365–376 Zur O, Berner Y, Ohel Y, Carmeli E (2018) Two-year
Yasui H, Takamoto K, Hori E, Urakawa S, Nagashima Y, follow-up of fall prediction among older adults in an
Yada Y, Ono T, Nishijo H (2010) Significant correla- independent-living community. Adv Exp Med Biol
tion between autonomic nervous activity and cerebral 1040:63–71
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 11–17
https://doi.org/10.1007/5584_2018_243
# Springer International Publishing AG, part of Springer Nature 2018
Published online: 21 August 2018
11
12 M. Bujar-Misztal and A. Chciałowski
pain when breathing, chest mobility impairments, groups: PNF (aged 62.3 8.4; n ¼ 22) and
cough reflex weakness, retention of secretions in standard rehabilitation (aged 62.0 7.4;
the bronchial tree, and sometimes lower lung n ¼ 20). The surgery was performed on a beating
lobes collapse, all of which may manifest in heart, without the use of extracorporeal circula-
lung dysfunction and arterial blood gas disorders tion. The principles of the rehabilitation methods
(El-Kader 2011; Wynne and Botti 2004; are presented in Table 1.
Westerdahl et al. 2001). Pulmonary rehabilita- During the 2-month period after surgery,
tion, whose task is to restore the patient’s opti- patients in both groups continued cardiac rehabil-
mum breathing efficiency, prevents itation according to the rules established during
complications mostly arising from physical inac- the hospital stay. Spirometry and body
tivity or immobilization (Pack et al. 2013; Leon plethysmography examinations were performed
et al. 2005), particularly after surgery. The using BodyBox 5500 (Medisoft, Limassol,
existing methods of physiotherapy take into Cyprus), according to the guidelines the
account the scope of hospital and posthospital European Respiratory Society (ERS) and Ameri-
treatment and the possibility of home-based reha- can Thoracic Society (ATS) (Curtis et al. 1994;
bilitation (Perk et al. 2012; Dylewicz et al. 2004). Efthimiou et al. 1991; Berrizbeitia et al. 1989).
Proprioceptive neuromuscular facilitation Table 2 shows the schedule of pulmonary func-
(PNF), known as Kabat-Kaiser rehabilitation, is tion examinations conducted in the study.
currently used in daily physiotherapy practice. In Data of respiratory variables are as the percent-
the method “P” stands for proprioception, i.e., age means SD of predicted values. The Shapiro-
sensory receptor-mediated collection of informa- Wilk test was used to assess normality of data
tion about the executed motion and spatial body distribution. A paired t-test or Wilcoxon’s test
positon, “N” for activation of neuromuscular was used to assess differences in each variable
structures, and “F” for facilitation. PNF between the time points of measurement in either
represents a holistic approach taking into account group patients of post-surgery rehabilitation. A
the therapeutic influence on the whole person, p-value <0.05 defined the statistically significant
including mental and social spheres, rather than differences. Calculations were carried out with R
just disease symptoms. The method is 2.13.0 software for Windows.
underpinned by a belief that bodily parts are inti-
mately interconnected and the treatment efforts
ought to be addressed to the whole (Adler et al. 3 Results
2008). Therefore, the aim of this study was to
determine the influence on lung function of 3.1 Spirometry
PNF vs. standard physiotherapy in patients after
CABG surgery. We found no signs of obstruction in any of the
patients subjected to either PNF or standard type
of cardiological rehabilitation after CABG sur-
2 Methods gery. The FEV1%FVC index was around 80%
and 90% in the PNF and standard rehabilitation
The study was approved by the Ethics Committee groups, respectively, and changed irrelevantly
for Scientific Research of the Academy of Physi- from the baseline level after 2-month-long post-
cal Education in Warsaw and by the Ethics Com- surgery rehabilitation (Table 3). In general, all the
mittee of the Military Institute of Medicine in variables were modestly reduced after 2-month-
Warsaw, Poland, and the patients gave informed long rehabilitation with either method; the reduc-
consent for study procedures. Forty-two patients, tion was more distinct in the PNF group in terms
aged 48–85, diagnosed with stable coronary of statistical significance, albeit the variables
artery disease, scheduled for elective CABG, remained in the normal level with the use of
were randomly divided into two physiotherapy both rehabilitative methods.
Influence of Proprioceptive Neuromuscular Facilitation on Lung Function in. . . 13
Table 3 Changes in spirometry variables from the base- standard method in patients with coronary artery bypass
line level after 2-month-long post-surgery rehabilitation graft (CABG)
using proprioceptive neuromuscular facilitation (PNF) or
PNF method Standard method
Variable Pre-surgery baseline 2 months post-surgery Pre-surgery baseline 2 months post-surgery
FVC 99.11 3.99 95.91 3.97 *** 104.98 3.77 96.98 3.65 **
FEV1 96.73 4.01 89.45 3.80 *** 99.54 4.05 90.73 3.66
FEV1%FVC 79.51 8.25 80.81 8.81 90.55 9.33 89.51 7.45
PEF 100.16 4.96 97.26 5.27 *** 104.96 4.02 98.69 4.84
FEF25 89.67 5.68 71.87 6.71* 71.75 8.00 58.33 7.19
FEF50 89.23 7.53 72.74 5.49 *** 90.34 8.73 76.92 6.22 *
FEF75 66.11 5.67 68.56 5.83* 98.45 6.90 93.59 5.74
Data are presented as % predicted value; *p < 0.05; **p < 0.01; ***p < 0.001
Table 4 Plethysmography variables in patients treated with proprioceptive neuromuscular facilitation (PNF) and
standard rehabilitation method at baseline pre-surgery and on the third day post-surgery
PNF method Standard method
Variable Pre-surgery baseline Day 3 post-surgery Pre-surgery baseline Day 3 post-surgery
VC 88.38 20.27 49.35 13.53 100.54 20.44 60.75 15.06 *
TGV 112.36 20.69 96.09 20.29 *** 126.82 31.36 101.72 28.34 **
RV 111.39 23.35 90.91 34.44 120.74 40.91 99.47 39.50 *
TLC 93.89 16.79 63.01 16.12 ** 99.89 20.26 69.21 18.76
RV%TLC 113.21 19.09 133.75 26.51 114.04 26.68 133.80 24.49
Data are presented as % predicted value; *p < 0.05; **p < 0.01; ***p < 0.001
Table 5 Plethysmography variables in patients treated with proprioceptive neuromuscular facilitation (PNF) and
standard rehabilitation method on the third and the seventh day post-surgery
PNF method Standard method
Variable Day 3 post-surgery Day 7 post-surgery Day 3 post-surgery Day 3 post-surgery
VC 49.35 13.53 64.31 20.94 * 60.75 15.06 71.68 16.54 ***
TGV 96.09 20.29 101.00 27.70 *** 101.72 28.34 115.56 17.17 **
RV 90.91 34.44 107.12 34.75 *** 99.47 39.50 120.00 35.33 **
TLC 63.01 16.12 77.94 20.77 *** 69.21 18.76 84.57 17.74 ***
RV%TLC 133.75 26.51 130.32 23.41 133.80 24.49 134.20 21.28 **
Data are presented as % predicted value; *p < 0.05; **p < 0.01; ***p < 0.001
Table 6 Plethysmography variables in patients treated with proprioceptive neuromuscular facilitation (PNF) and
standard rehabilitation method on the seventh day and then 2 months post-surgery
PNF method Standard method
Variable Day 3 post-surgery 2 months post-surgery Day 3 post-surgery 2 months post-surgery
VC 64.31 20.94 88.51 21.66 *** 71.68 16.54 97.15 16.61
TGV 101.00 27.70 113.23 29.07 *** 115.56 17.17 118.47 23.49 ***
RV 107.12 34.75 109.02 39.53 *** 120.00 35.33 110.75 29.06
TLC 77.94 20.77 92.33 20.77 *** 84.57 17.74 97.17 16.59 *
RV%TLC 130.32 23.41 109.30 30.12 *** 134.20 21.28 108.17 19.50
Data are presented as % predicted value; *p < 0.05; ***p < 0.001
The results of body plethysmography are post-surgery in both groups. The examination
presented in detail in Tables 4, 5, 6, and 7. conducted 2 months post-surgery showed a sub-
There were a number of significant differences stantial increase in the variables, with an advan-
in specific variables, beyond those in the RV% tage for the PNF treatment.
TLC index above outlined. However, RV%TLC
is a reflection of lung air trapping after CABG
surgery, and thus it is of key importance for lung
4 Discussion
health.
It is worth noticing that the plethysmography
Rehabilitation is an important part of the treat-
variables were generally lower in patients treated
ment for patients after cardiovascular and thoracic
with PNF than those with standard rehabilitation.
surgeries, as it helps improve physical and respi-
The CABG caused a significant decrease in the
ratory capacity. Breathing exercises lead to a
value of individual variables, particularly on Day
reduction or complete elimination of potential
3 post-surgery, with a gradual rebound on Day 7
pulmonary complications resulting from impaired
Influence of Proprioceptive Neuromuscular Facilitation on Lung Function in. . . 15
Table 7 Plethysmography variables in patients treated with proprioceptive neuromuscular facilitation (PNF) and
standard rehabilitation method at baseline pre-surgery and 2 months post-surgery
PNF method Standard method
Variable Pre-surgery baseline 2 months post-surgery Pre-surgery baseline 2 months post-surgery
VC 88.38 20.27 88.51 21.66 *** 100.54 20.44 97.15 16.61
TGV 112.36 20.69 113.23 29.07 *** 126.82 31.36 118.47 23.49 **
RV 111.39 23.35 109.02 39.53 120.74 40.91 110.75 29.06
TLC 93.89 16.79 92.33 20.77 ** 99.89 20.26 97.17 16.59
RV%TLC 113.21 19.09 109.30 30.12 114.04 26.68 108.17 19.50
Data are presented as % predicted value; **p < 0.01; ***p < 0.001
respiratory mechanics, increasing pain, cough, changes in the lungs. In a study of Shenkman
phlegm retention, or diaphragm impairment et al. (1997), impaired lung function, expressed
(Efthimiou et al. 1991). as the 25% reductions in FVC, FEV1, and PEF,
The present study enabled to relate functional was observed 3–4 months after surgery, while
breathing variables in patients after coronary reductions in FEF50 and FEF75 were sustained
artery bypass graft surgery to the pre-surgery for a longer time. A decline in lung function is
baseline and to compare the effects on lung func- due likely to a modified breathing pattern after
tion of PNF and standard rehabilitation methods chest opening. Damage to ribs limits the chest
in the extended post-surgery period of up to wall motion, and respiratory muscles weakness
2 months. We, generally, noticed reductions in usually persist for several months. In addition,
spirometry indices 2 months after revasculariza- atelectatic changes in the basal segments of
tion, pointing to impairment of airflow in large lower lung lobes and pleural damage contribute
and small bronchi, irrespective of rehabilitation. to respiratory dysfunction (Montes et al. 2004).
The reduction was somehow less in the The diaphragm dysfunction caused by a transient
PNF-rehabilitated patients. The FEV1/FVC phrenic nerve damage as a result of heart hypo-
index did not change significantly at the post- thermia also is a serious adverse effect that occurs
surgery time marks studied, which was due likely in about 22–31% of patients, even as late as up to
to the proportional reduction of both components. 1 year after surgery (Westerdahl et al. 2003, 2005;
Body plethysmography also revealed a decline in Curtis et al. 1994; Efthimiou et al. 1991). A
individual variables on the third and the seventh weaker cardiac mechanical function with low
day post-surgery. The decline was associated with ejection fraction may predispose to the develop-
increased value of RV%TLC that became ment of foci of pulmonary edema and an abnor-
normalized 2 months later but only in patients mal ventilation/perfusion ratio, followed by
undergoing PNF rehabilitation. hypoxia (Staton et al. 2005). Vaidya et al.
The PNF method apparently has never before (1996) have shown that spirometry variables
been used for rehabilitation after CABG surgery, decrease after mitral valve surgery for up to
so there is no comparative reference to the present 3 months, despite improved hemodynamic heart
findings. However, the general impairment of function.
lung function after CABG operations is in line Staton et al. (2005) have shown significant
with the results of other studies. A decline, some- reductions in FEV1 and FVC, with unchanged
times reaching 40–50%, in FVC and FEV1 has FEV1/FVC ratio, compared with the pre-surgery
been observed in the first 3 days post-surgery, baseline. However, the adoption of FEV1/FVC as
gradually normalizing after a few months the exponent of a possible restriction has been due
(Wynne and Botti 2004; Dubach et al. 1995; largely to the lack of TLC measurement.
Berrizbeitia et al. 1989). Rachwalik et al. (2006) Rachwalik et al. (2006, 2007), after a triple exam-
have observed a slight increase in FEV1/FVC, ination of patients, on the fifth day and 2.5 and
which points to the trend toward restrictive 6 months after cardiopulmonary bypass surgery,
16 M. Bujar-Misztal and A. Chciałowski
Nichols M, Townsend N, Scarborough P, Rayner M Shenkman Z, Shir Y, Weiss Y, Bleiberg B, Gross D (1997)
(2014) Cardiovascular disease in Europe 2014: epide- The effects of cardiac surgery on early and late pulmo-
miological update. Eur Heart J 35(42):2950–2959 nary functions. Acta Anaesthesiol Scand
Pack QR, Goel K, Lahr BD, Greason KL, Squires RW, 41:1193–1199
Lopez-Jimenez F, Zhang Z, Thomas RJ (2013) Partici- Staton GW, Williams WH, Mahoney EM, Hu J, Chu H,
pation in cardiac rehabilitation and survival after coro- Duke PG, Puskas JD (2005) Pulmonary outcomes of
nary artery bypass graft surgery: a community based off-pump vs on-pump coronary artery bypass surgery
study. Circulation 128(6):590–597 in a randomized trial. Chest 127(3):892–901
Perk J, DeBacker G, Gohlke H, Graham I, Reiner Z, Vaidya R, Husain T, Ghosh PK (1996) Spirometric
Verschuren WM (2012) European guidelines on car- changes after open mitral surgery. J Cardiovasc Surg
diovascular disease prevention in clinical practice. The 37(3):295–300
fifth joint task force of the European Society of Cardi- Westerdahl E, Lindmark B, Almgren SO, Tenling A
ology and Other Societies on cardiovascular disease (2001) Chest physiotherapy after coronary artery
prevention in clinical practice. Eur Heart J bypass graft surgery – a comparison of three different
33:1635–1701 deep breathing techniques. J Rehab Med 33:79–84
Pinheiro AC, Magalhaes Novais MC, Neto MG, Herbst Westerdahl E, Lindmark B, Bryngelsson I, Tenling A
Rodrigues MV, de Souza Rodrigues E Jr, Aras R Jr, (2003) Pulmonary function 4 months after coronary
Carvalho VO (2011) Estimation of lung vital capacity artery bypass graft surgery. Respir Med 97:317–322
before and after coronary artery bypass grafting sur- Westerdahl E, Lindmark B, Eriksson T, Friberg Ö,
gery: a comparison on incentive spirometry and Hedenstierna G, Tenling A (2005) Deep-breathing
ventilometry. J Cardiothorac Surg 6:70 exercises reduce atelectasis and improve pulmonary
Rachwalik M, Lysenko L, Kustrzycki W, Pelczar M, function after coronary artery bypass surgery. Chest
Wachnik J (2006) Comparison of spirometry in 128:3482–3488
patients after on-pump and off-pump coronary artery Wynne R, Botti M (2004) Postoperative pulmonary dys-
bypass grafting. Anaesthesiol Intensive Ther 1:6–10 function in adults after cardiac surgery with cardiopul-
Rachwalik M, Pelczar M, Gozdzik A (2007) Postopera- monary bypass: clinical significance and implications
tive, early and late pulmonary functions – off– for practice. Am J Crit Care 13(5):384–393
pump vs. on–pump coronary artery bypass surgery.
Pol J Thoracic Cardiovasc Surgery 4(3):252–260
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 19–25
https://doi.org/10.1007/5584_2018_282
# Springer Nature Switzerland AG 2018
Published online: 29 September 2018
19
20 M. Wojciechowska et al.
1 Introduction 2 Methods
Table 2 Remote ischemic preconditioning (RICP) and incidence of contrast-induced nephropathy (CIN)
RICP (n ¼ 62) Control (n ¼ 61) p
Serum creatinine 24-B, mg/dL 0.01 0.09 0.04 0.10 0.087
Serum creatinine 72-B, mg/dL 0.03 0.09 0.04 0.08 0.623
Serum cystatin C 24-B, mg/dL 0.00 0.13 0.01 0.11 0.827
CIN (SC) 3/61 (4.9%) 7/58 (12.1%) 0.197
CIN (CC) 1/58 (1.7%) 2/57 (3.5%) 0.618
CIN (CC) contrast-induced nephropathy assessed by serum cystatin C criteria, CIN (SC) contrast-induced nephropathy
assessed by serum creatinine criteria, Serum creatinine 24-B difference between serum creatinine concentration
measured after 24 h and the baseline (B) level, Serum creatinine 72-B difference between serum creatinine concentration
measured after 72 h and the baseline level, Serum cystatin C 24-B difference between serum cystatin C concentration
measured after 24 h and the baseline level
Remote Ischemic Preconditioning in Renal Protection During Elective. . . 23
basis of SC. There was no coincidence between been demonstrated to prevented contrast-induced
the diagnosis of contrast-induced nephropathy nephropathy in patients with ST-elevation
based on SC and CC (McNemar test: 0.02, myocardial infarction (Yamanaka et al. 2015) or
κ ¼ 0.28). Serum creatinine measured twice in coexisting chronic kidney disease (Igarashi
(24 and 72 h after PCI) appeared to more accu- et al. 2013; Er et al. 2012), with a high Mehran
rately identify all patients with contrast-induced risk score (Mehran et al. 2004). After a thorough
nephropathy than CC measured once (from analysis of all the studies, along with the present
10 cases of contrast-induced nephropathy based results, we conclude that RIPC may be of less
on SC, only 3 met criteria to recognize contrast- importance in case of elective PCI and low-risk
induced nephropathy based on CC). patients, but it may be particularly useful in unsta-
ble patients with a high to very high risk of
contrast-induced nephropathy.
3.4 One-Year Follow-Up The mechanisms underlying contrast-induced
nephropathy are multifactorial, and it is possible
One-year follow-up was possible in case of that they vary depending on comorbidities and
95 patients: 9 (90%) with contrast-induced indications for PCI (elective vs. urgent). In
nephropathy (CIN+) and 86 (79%) without it lower-risk patients, the plausible mechanism of
(CIN–). The RIPC had no appreciable influence contrast-induced nephropathy may be a direct
on the occurrence of major adverse cardiac and toxicity of contrast on tubular cells, caused or
cerebrovascular events during the 1-year follow- aggravated by patient’s dehydration (Sendeski
up, nor did it significantly change the MACCE 2011). In case of higher-risk patients, the
index (Table 3). pathomechanism may be more complicated as
hypoxic injury may be a predominant factor.
RIPC may prevent hypoxic injury and thus may
4 Discussion have a renoprotective effect in higher-risk
patients, but not in the general population.
In this study, general risk for developing contrast- Urea-to-creatinine (U-to-SC) ratio is used to
induced nephropathy was 8.4% (10/119 patients), determine the patient hydration. An abnormal
and it was compatible with the predicted proba- value that could indicate hypovolemia and a
bility based on the Mehran risk score (Mehran reduced renal perfusion is more than 40 (Higgins
et al. 2004). The study demonstrates that RIPC 2016). In the study, baseline mean U-to-SC ratio
applied prior to elective PCI failed to provide an in the patients of both groups was 45.1 1.3, and
added benefit to hydration and N-acetylcysteine it was over 40 in 63.4% of patients. That indicates
in decreasing the incidence of contrast-induced that more than 60% of patients could be
nephropathy. There was a trend toward lower hypovolemic on admission. We demonstrate a
numerical incidence of contrast-induced nephrop- positive association between higher baseline U-
athy, assessed on the basis of both SC and CC, to-SC ratio and contrast-induced nephropathy
and also lower increase in SC and CC levels after (55.1 14.8 in CIN+ vs. 43.9 10.0 in CIN–;
PCI in the RIPC group; the differences were p < 0.002), which means that more adequate fluid
however insignificant. As the risk for developing therapy should be implemented in selected
contrast-induced nephropathy was generally low patients for preventing nephropathy. In this
in our patients (Mehran risk score about 3 points; study, patients did not interrupt treatment with
Table 1), these results are in line with similar potentially nephrotoxic drugs, such as ACEI and
studies in the literature, where RIPC failed to AT1 blockers, so we indirectly demonstrate that
improve kidney function after contrast adminis- such drugs did not increase the risk of
tration in patients with low-to-moderate risk of nephropathy.
contrast-induced nephropathy (Gholoobi et al. CC is considered an early marker of kidney
2015; Menting et al. 2015). However, RIPC has injury. SC typically goes up between 48–72 h
24 M. Wojciechowska et al.
Table 3 Effects on the occurrence of major adverse cardiac and cerebrovascular events of remote ischemic
preconditioning (RIPC) during 1-year follow-up
RICP (n ¼ 62) Control (n ¼ 61)
n ¼ 47 (75.8%) n ¼ 48 (78.7%) p
TIA/stroke n (%) 1 (2.1%) 1 (2.1%) 1.000
ACS (UA/MI) n (%) 4 (8.5%) 5 (10.4%) 1.000
MACCE n (%) 4 (8.5%) 6 (12.5%) 0.740
TIA Transient ischemic attack, ACS (UA/MI) acute coronary syndrome (unstable angina/myocardial infarct), MACCE
composite endpoint defined as acute coronary syndrome, transient ischemic attack, or stroke
after the use of contrast, whereas CC reaches the after percutaneous coronary intervention does not
highest value 12–24 h after PCI (Herget- detect all cases of contrast-induced nephropathy.
Rosenthal et al. 2004; Dharnidharka et al. 2002; Remote ischemic preconditioning before percuta-
Mussap et al. 2002). As patients are usually neous coronary intervention does not influence
discharged from the hospital 24 h after PCI, CC 1-year follow-up outcome.
might be a better option to identify contrast-
induced nephropathy. The present results indi- Conflicts of Interest The author declares no conflicts of
cate, however, that SC measured twice, 24 and interest in relation to this article.
72 h after PCI, may be more accurate to detect
nephropathy as we found ten CIN+ cases based Ethical Approval All procedures in this study were
performed in accordance with the ethical standards of the
on SC and three on CC.
institutional research committee and with the 1964
In previous studies, RIPC performed before Helsinki Declaration for Human Research and its later
PCI has been shown to improve clinical outcomes amendments.
during the following 4 (Yamanaka et al. 2015) or
6 weeks (Er et al. 2012). In this study, RIPC, Informed Consent Written informed consent was
performed before elective PCI, failed to influence obtained from all individual participants included in the
study.
the incidence of contrast-induced nephropathy or
the MACCE index during in 1-year follow-up.
Nonetheless, we cannot rule out the appearance
of early post-PCI advantages of RIPC, which References
would wane with the time of follow-up. The
Bartorelli AL, Marenzi G (2008) Contrast-induced
main limitation of this study was a small number nephropathy. J Interv Cardiol 21:74–85
of patients investigated. It cannot be ruled out that Dharnidharka VR, Kwon C, Stevens GA (2002) Serum
in a larger population of patients the differences cystatin C is superior to serum creatinine as a marker of
kidney function: a meta-analysis. J Kidney Dis
between RIPC-treated and control groups would
40:221–226
become significant. Er F, Nia AM, Dopp H, Hellmich M, Dahlem KM,
Caglayan E, Kubacki T, Benzing T, Erdmann E,
Burst V, Gassanov N (2012) Ischemic preconditioning
for prevention of contrast medium-induced nephropa-
5 Conclusions thy: randomized pilot RenPro Trial (Renal Protection
Trial). Circulation 126:296–303
Gholoobi A, Sajjadi SM, Shabestari MM, Eshraghi A,
We believe we have demonstrated that in patients
Shamloo AS (2015) The impact of remote ischemic
referred for elective percutaneous coronary inter- pre-conditioning on contrast-induced nephropathy in
vention, remote ischemic preconditioning does patients undergoing coronary angiography and angio-
not provide statistically significant benefits in plasty: a double-blind randomized clinical trial. Elec-
tron Physician 7:1557–1565
decreasing the incidence of contrast-induced
Goldenberg I, Matetzky S (2005) Nephropathy induced by
nephropathy. Cystatin C is considered an early contrast media: pathogenesis, risk factors and preven-
marker of kidney injury, but measuring it 24 h tive strategies. CMAJ 172:1461–1471
Remote Ischemic Preconditioning in Renal Protection During Elective. . . 25
Hausenloy DJ, Yellon DM (2008) Remote ischaemic Mussap M, Dalla Vostra M, Fioretto P, Saller A,
preconditioning: underlying mechanisms and clinical Varagnolo M, Nosadini R, Plebani M (2002) Cystatin
application. Cardiovasc Res 79:377–386 C is a more sensitive marker than creatinine for the
Herget-Rosenthal S, Marggraf G, Husing J, Göring F, estimation of GFR in type 2 diabetic patients. Kidney
Pietruck F, Janssen O, Philipp T, Kribben A (2004) Int 61:1453–1461
Early detection of acute renal failure by serum cystatin Raingruber B, Kirkland-Walsh H, Chahon N, Kellermann
C. Kidney Int 66:1115–1122 M (2011) Using the Mehran risk scoring tool to predict
Heyman SN, Rosen S, Brezis M (1994) Radiocontrastne- risk for contrast medium–induced nephropathy in
phropathy: a paradigm for the synergism between toxic patients undergoing percutaneous angiography. Crit
and hypoxic insults in the kidney. Exp Nephrol Care Nurse 31:e17–e22
2:153–157 Rihal CS, Textor SC, Grill DE, Berger PB, Ting HH, Best
Higgins C (2016) Urea and creatinine concentration, the PJ, Singh M, Bell MR, Barsness GW, Mathew V,
urea: creatinine ratio. www.acutecaretesting.org. Garratt KN, Holmes DR Jr (2002) Incidence and prog-
Accessed 24 Aug 2018 nostic importance of acute renal failure after percuta-
Igarashi G, Iino K, Watanabe H, Ito H (2013) Remote neous coronary intervention. Circulation
ischemic pre-conditioning alleviates contrast-induced 105:2259–2264
acute kidney injury in patients with moderate chronic Sendeski MM (2011) Pathophysiology of renal tissue
kidney disease. Circ J 77:3037–3044 damage by iodinated contrast media. Clin Exp
McCullough PA, Adam A, Becker CR, Davidson C, Pharmacol Physiol 38:292–299
Lameire N, Stacul F, Tumlin J, CIN Consensus Work- Wong GT, Irwin MG (2007) Contrast-induced nephropa-
ing Panel (2006) Risk prediction of contrast-induced thy. Br J Anaesth 99:474–483
nephropathy. Am J Cardiol 98:27K–36K Yamanaka T, Kawai Y, Miyoshi T, Mima T, Takagaki K,
Mehran R, Aymong ED, Nikolsky E, Lasic Z, Iakovou I, Tsukuda S, Kazatani Y, Nakamura K, Ito H (2015)
Fahy M, Mintz GS, Lansky AJ, Moses JW, Stone GW, Remote ischemic preconditioning reduces contrast-
Leon MB, Dangas G (2004) A simple risk score for induced acute kidney injury in patients with
prediction of contrast-induced nephropathy after per- ST-elevation myocardial infarction: a randomized con-
cutaneous coronary intervention: development and ini- trolled trial. Int J Cardiol 178:136–141
tial validation. J Am Coll Cardiol 44:1393–1399 Zimmerman RF, Ezeanuna PU, Kane JC, Cleland CD,
Menting TP, Sterenborg TB, de Waal Y, Donders R, Kempananjappa TJ, Lucas FL, Kramer RS (2011)
Wever KE, Lemson MS, van der Vliet JA, Wetzels Ischemic preconditioning at a remote site prevents
JF, SchultzeKool LJ, Warlé MC (2015) Ischemic acute kidney injury in patients following cardiac sur-
preconditioning to reduce contrast-induced nephropa- gery. Kidney Int 80:861–867
thy: a randomized controlled trial. Eur J Vasc
Endovasc Surg 50:527–532
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 27–36
https://doi.org/10.1007/5584_2018_238
# Springer International Publishing AG, part of Springer Nature 2018
Published online: 29 June 2018
27
28 S. V. Tabatabaei et al.
other malignancies; (2) lack of data about prog- Scores 5 out of the 9 possible were judged to
nostic parameters in the title/abstract; (3) insuffi- indicate a high risk of bias.
cient data to enable a statistical comparison of The pooled risk ratios (RR) and 95% confi-
survival between ECE+ and ECE– patients; dence intervals (CI) of all-cause mortality
(4) NSCLC stage IV; (5) palliative therapy of a between ECE+ and ECE patients were calcu-
subgroup of patients without sufficient data to lated using DerSimonian–Laird random-effect
exclude them in a secondary analysis; models (DerSimonian and Laird 1986). Hetero-
(6) in vitro or animal studies; (7) neo-adjuvant geneity across studies was assessed by the I2
chemoradiotherapy; and (8) nodal status N3. metric and χ2 statistics (Higgins and Thompson
2002).
Poor Vascular
Mean/Median differentiation invasion Squamous cell carcinoma Pneumonectomy R0–Resection
age (year) Male (%) T4 (%) N2 (%) (%) (%) (%) Smoking (%) Stage IIIA (%) (%) (%)
Study ECE ECE ECE ECE ECE ECE ECE ECE ECE
Author country NOS + ECE + ECE + ECE + ECE + ECE ECE+ ECE + ECE + ECE + ECE ECE+ ECE + ECE
Duran Spain 7 NR NR NR 33 NR NR NR NR 37 NR NR
Cantolla et al.
(1994)
NOS Newcastle–Ottawa Scale score, ECE+ with extracapsular tumor extension, ECE without extracapsular tumor extension, T4 locally advanced lung cancer, N2 cancer is in lymph nodes,
Stage IIIA cancer spread to lymph nodes on the same side of the chest as the tumor, R0-Resection a microscopically margin-negative resection, NR not reported
Prognostic Impact of Extracapsular Lymph Node Invasion on Survival in. . . 31
Fig. 1 Summary of survival risk ratios (RR) in original studies and the pooled estimate; RR risk ratio
32 S. V. Tabatabaei et al.
reserved only for pN2 and pN3 status. Since the This result, in conjunction with our finding,
findings of the present review suggest a decreased suggests that ECE plays an important role in the
survival in the ECE+ group, we hypothesize that prognosis of solid tumors as well as NSCLC, and
PORT could benefit the subgroup of pN1 patients consequently it should be incorporated into the
with extra-nodal lymph node invasion. It would TNM staging criteria.
be interesting whether this approach increases It is important to know that the penetration of
survival compared to no pN1-PORT in NSCLC the lymph node capsule can be investigated after
(Borghetti et al. 2016; Zaric et al. 2013), similar surgical lymph node dissection (Gu et al. 2017).
to pN1 breast cancer where PORT has been Since ECE is very focal, the surrounding fatty
shown to benefit ECE+ patients (Belkacemi tissue has to be completely included in surgical
et al. 2015). Lynch et al. (2014) have shown that resection of lymph nodes, in order to make a valid
in patients with squamous cell carcinoma of the statement about tumor extension (Veronese et al.
tonsil, ECE status is predictive of contralateral 2016). In lung cancer patients, endobronchial
nodal recurrence. Consequently, prophylactic ultrasound-guided transbronchial needle aspira-
contralateral radiotherapy has been suggested to tion (EBUS-TBNA) has recently been widely
improve overall survival in these patients. The used for mediastinal lymph nodes staging and as
same could also be true for NSCLC patients a minimal invasive staging tool has replaced
with lymph node status N1 (Borghetti et al. video-mediastinoscopy. EBUS-TBNA provides
2016). More data is required to discuss this issue both cytologic and histologic materials, but it
in NSCLC, particularly with lymph node status cannot provide reliable information regarding
N2 and N3. ECE status.
Furthermore, Fig. 3 represents a graphical In the present systematic review and meta-
reconstruction of summarized all-cause mortality analysis, we examined three prospective and two
data from studies included in three systematic retrospective observational studies involving
reviews published in 2016 and 2007 (Luchini 828 patients affected by NSCLC with metastatic
et al. 2016; Veronese et al. 2016; Wind et al. lymph nodes. Of these, 298 patients had ECE+ of
2007). As seen, ECE+ is associated with a signif- lymph node metastasis, while 530 showed only
icantly higher all-cause mortality in solid tumors. intra-nodal metastasis. While the results of this
Prognostic Impact of Extracapsular Lymph Node Invasion on Survival in. . . 33
Fig. 3 Pooled estimate of all-cause mortality in different types of cancer. The data is summarized using a random-effect
model; RR risk ratio
value of the number and laterality of metastatic ingui- with non–small–cell lung carcinoma (NSCLC) who
nal lymph nodes in vulvar cancer: revisiting the FIGO have undergone curative treatment. Arch
staging system. Eur J Surg Oncol 39(7):780–785 Bronconeumol 30(4):196–201 (Article in Spanish)
Belkacemi Y, Kaidar-Person O, Poortmans P, Ozsahin M, Eckel SP, Cockburn M, Shu YH, Deng H, Lurmann FW,
Valli MC, Russell N, Kunkler I, Hermans J, Kuten A, Liu L, Gilliland FD (2016) Air pollution affects lung
van Tienhoven G, Westenberg H, Breast Working cancer survival. Thorax 71(10):891–898
Party of the EORTC Radiation Oncology Group Ettinger DS, Wood DE, Akerley W, Bazhenova LA,
(ROG) (2015) Patterns of practice of regional nodal Borghaei H, Camidge DR, Cheney RT, Chirieac LR,
irradiation in breast cancer: results of the European D’Amico TA, Demmy TL, Dilling TJ, Govindan R,
Organization for Research and Treatment of Cancer Grannis FW Jr, Horn L, Jahan TM, Komaki R, Kris
(EORTC) NOdal radiotherapy (NORA) survey. Ann MG, Krug LM, Lackner RP, Lanuti M, Lilenbaum R,
Oncol 26(3):529–535 Lin J, Loo BW Jr, Martins R, Otterson GA, Patel JD,
Bollen EC, Theunissen PH, van Duin CJ, Drenth BM, van Pisters KM, Reckamp K, Riely GJ, Rohren E, Schild S,
Noord JA, Blijham GH (1994) Clinical significance of Shapiro TA, Swanson SJ, Tauer K, Yang SC,
intranodal and extranodal growth in lymph node Gregory K, Hughes M (2014) Non-small cell lung
metastases of non-small cell lung cancer. Scand J cancer, version 1.2015. J Natl Compr Cancer Netw
Thorac Cardiovasc Surg 28(3–4):97–102 12(12):1738–1761
Borghetti P, Barbera F, Bonù ML, Trevisan F, Ciccarelli S, Fons G, Hyde SE, Buist MR, Schilthuis MS, Grant P,
Vitali P, Maddalo M, Triggiani L, Pasinetti N, Burger MP, van der Velden J (2009) Prognostic value
Pedretti S, Bonetti B, Pariscenti G, Tironi A, of bilateral positive nodes in squamous cell cancer of
Caprioli A, Buglione M, Magrini SM (2016) Resected the vulva. Int J Gynecol Cancer 19(7):1276–1280
pN1 non-small cell lung cancer: recurrence patterns Fujii T, Yajima R, Tatsuki H, Kuwano H (2016) Prediction
and nodal risk factors may suggest selection criteria of extracapsular invasion at metastatic sentinel nodes
for post-operative radiotherapy. Radiol Med 121 and non-sentinel lymph nodal metastases by FDG-PET
(9):696–703 in cases with breast cancer. Anticancer Res 36
Brabender J, Bollschweiler E, Hölscher AH, Strobel K, (4):1785–1789
Gutschow C, Prenzel K, Grimminger P, Drebber U, Gatzinsky P, Berglin E, Dernevik L, Larsson I, William-
Schröder W, Metzger R, Vallböhmer D (2012) The Olsson G (1985) Resectional operations and long–term
prognostic impact of extracapsular lymph node results in carcinoma of the esophagus. J Thorac
involvement in rectal cancer patients: implications for Cardiovasc Surg 89(1):71–76
staging and adjuvant treatment strategies. Oncol Lett 3 Goldstraw P, Crowley J, Chansky K, Giroux DJ, Groome
(4):825–830 PA, Rami-Porta R, Postmus PE, Rusch V, Sobin L,
Bradley JD, Paulus R, Graham MV, Ettinger DS, International Association for the Study of Lung Cancer
Johnstone DW, Pilepich MV, Machtay M, Komaki R, International Staging Committee, Participating
Atkins J, Curran WJ, Groupet RTO (2005) Phase II Institutions (2007) The IASLC lung cancer staging
trial of postoperative adjuvant paclitaxel/carboplatin project: proposals for the revision of the TNM stage
and thoracic radiotherapy in resected stage II and IIIA groupings in the forthcoming (seventh) edition of the
non-small-cell lung cancer: promising long-term TNM classification of malignant tumours. J Thorac
results of the radiation therapy oncology group-- Oncol 2(8):706–714
RTOG 9705. J Clin Oncol 23(15):3480–3487 Gomez DR, Komaki R (2012) Postoperative radiation
Chen WC, Lai CH, Fang CC, Yang YH, Chen PC, Lee CP, therapy for non-small cell lung cancer and thymic
Chen MF (2016) Identification of high–risk subgroups malignancies. Cancers (Basel) 4(1):307–322
of patients with oral cavity cancer in need of postoper- Gu Y, Shi H, Su C, Chen X, Zhang S, Li W, Wu F, Gao G,
ative adjuvant radiotherapy or chemo–radiotherapy. Wang H, Chu H, Zhou C, Zhou F, Ren S (2017) The
Medicine (Baltimore) 95(22):e3770 role of endobronchial ultrasound elastography in the
Dai C, Ren Y, Xie D, Zheng H, She Y, Fei K, Jiang G, diagnosis of mediastinal and hilar lymph nodes.
Chen C (2016) Does lymph node metastasis have a Oncotarget 8(51):89194–89202
negative prognostic impact in patients with NSCLC Higgins, Thompson SG (2002) Quantifying heterogeneity
and M1a disease? J Thorac Oncol 11(10):1745–1754 in a meta-analysis. Stat Med 21(11):1539–1558
DerSimonian R, Laird N (1986) Meta-analysis in clinical Hui Z, Dai H, Liang J, Lv J, Zhou Z, Feng Q, Xiao Z, Chen
trials. Control Clin Trials 7(3):177–188 D, Zhang H, Yin W, Wang L (2015) Selection of
Di Giorgio A, Botti C, Sammartino P, Mingazzini P, proper candidates with resected pathological stage
Flammia M, Stipa V (1991) Extracapsular lymph IIIA-N2 non-small cell lung cancer for postoperative
node metastases in the staging and prognosis of gastric radiotherapy. Thorac Cancer 6(3):346–353
cancer. Int Surg 76(4):218–221 Igai H, Matsuura N, Tarumi S, Chang SS, Misaki N, Go T,
Duran-Cantolla J, González Macías JJ, Agüero Balbín R, Ishikawa S, Yokomise H (2011) Clinicopathological
Carbajo Carbajo M, Ortega Morales FJ, Hernández study of p-T1aN0M0 non-small-cell lung cancer, as
Alonso M, Ondiviela Gracia R (1994) The prognosis defined in the seventh edition of the TNM
of mediastinal lymph node infiltration (N2) in patients
Prognostic Impact of Extracapsular Lymph Node Invasion on Survival in. . . 35
classification of malignant tumors. Eur J Cardiothorac Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gøtzsche
Surg 39(6):963–967 PC, Ioannidis JP, Clarke M, Devereaux PJ, Kleijnen J,
Jing X et al (2016) Delineation of clinical target volume Moher D (2009) The PRISMA statement for reporting
for postoperative radiotherapy in stage IIIA–pN2 non– systematic reviews and meta-analyses of studies that
small–cell lung cancer. Onco Targets Ther 9:823–831 evaluate healthcare interventions: explanation and
Jonnalagadda S, Smith C, Mhango G, Wisnivesky JP elaboration. BMJ 339:b2700
(2011) The number of lymph node metastases as a Liu W, Shao Y, Guan B, Hao J, Cheng X, Ji K, Wang K
prognostic factor in patients with N1 non-small cell (2015) Extracapsular extension is a powerful prognos-
lung cancer. Chest 140(2):433–440 tic factor in stage IIA–IIIA non–small cell lung cancer
Kelder W, Ebrahimi A, Forest VI, Gao K, Murali R, Clark patients with completely resection. Int J Clin Exp
JR (2012) Cutaneous head and neck squamous cell Pathol 8(9):11268–11277
carcinoma with regional metastases: the prognostic Luchini C, Nottegar A, Solmi M, Sergi G, Manzato E,
importance of soft tissue metastases and extranodal Capelli P, Scarpa A, Veronese N (2016) Prognostic
spread. Ann Surg Oncol 19(1):274–279 implications of extranodal extension in node-positive
Kerr KM, Nicolson MC (2013) Prognostic factors in squamous cell carcinoma of the vulva: a systematic
resected lung carcinomas. EJC Suppl 11(2):137–149 review and meta-analysis. Surg Oncol 25(1):60–65
Khosrotehrani K, van der Ploeg AP, Siskind V, Hughes Lupattelli M, Maranzano E, Bellavita R, Tarducci R,
MC, Wright A, Thomas J, Barbour A, Allan C, Latini R, Castagnoli P, Bufalari A, Corgna E,
Bayley G, Eggermont A, Verhoef C, Smithers BM, Pinaglia D, Rossetti R, Ribacchi R, Latini P (2001)
Green AC (2014) Nomograms to predict recurrence Adjuvant radiochemotherapy in high–risk rectal cancer
and survival in stage IIIB and IIIC melanoma after results of a prospective non–randomized study. Tumori
therapeutic lymphadenectomy. Eur J Cancer 50 87(4):239–247
(7):1301–1309 Lynch J, Lal P, Schick U, Nutting CM, Newbold K,
Komori K, Kanemitsu Y, Kimura K, Sano T, Ito S, Abe T, Harrington K, Bhide S (2014) Multiple cervical
Senda Y, Shimizu Y (2013) Detailed stratification of lymph node involvement and extra–capsular extension
TNM stage III rectal cancer based on the presence/ predict for contralateral nodal recurrence after ipsilat-
absence of extracapsular invasion of the metastatic eral radiotherapy for squamous cell carcinoma of the
lymph nodes. Dis Colon Rectum 56(6):726–732 tonsil. Oral Oncol 50(9):901–906
Komuta K, Okudaira S, Haraguchi M, Furui J, Kanematsu Maurizi G, D’Andrilli A, Venuta F, Rendina EA (2016)
T (2001) Identification of extracapsular invasion of the Reconstruction of the bronchus and pulmonary artery.
metastatic lymph nodes as a useful prognostic sign in J Thorac Dis 8(Suppl 2):S168–S180
patients with resectable colorectal cancer. Dis Colon Metindir J, Bilir Dilek G (2008) Evaluation of prognostic
Rectum 44(12):1838–1844 significance in extracapsular spread of pelvic lymph
Kuo SW, Chen JS, Huang PM, Hsu HH, Lai HS, Lee JM node metastasis in patients with cervical cancer. Eur J
(2014) Prognostic significance of histologic differenti- Gynaecol Oncol 29(5):476–478
ation, carcinoembryonic antigen value, and Metzger R, Drebber U, Baldus SE, Mönig SP, Hölscher
lymphovascular invasion in stage I non-small cell AH, Bollschweiler E (2009) Extracapsular lymph node
lung cancer. J Thorac Cardiovasc Surg 148 involvement differs between squamous cell and adeno-
(4):1200–1207.e3 carcinoma of the esophagus. Ann Surg Oncol 16
Lagarde SM, ten Kate FJ, de Boer DJ, Busch OR, (2):447–453
Obertop H, van Lanschot JJ (2006) Extracapsular Misthos P, Sepsas E, Kokotsakis J, Skottis I, Lioulias A
lymph node involvement in node–positive patients (2008) The significance of one-station N2 disease in
with adenocarcinoma of the distal esophagus or gas- the prognosis of patients with nonsmall-cell lung can-
troesophageal junction. Am J Surg Pathol 30 cer. Ann Thorac Surg 86(5):1626–1630
(2):171–176 Moretti L, Yu DS, Chen H, Carbone DP, Johnson DH,
Le Péchoux C (2011) Role of postoperative radiotherapy Keedy VL, Putnam JB Jr, Sandler AB, Shyr Y, Lu B
in resected non-small cell lung cancer: a reassessment (2009) Prognostic factors for resected non-small cell
based on new data. Oncologist 16(5):672–681 lung cancer with pN2 status: implications for use of
Lee YC, Wu CT, Kuo SW, Tseng YT, Chang YL (2007) postoperative radiotherapy. Oncologist 14
Significance of extranodal extension of regional lymph (11):1106–1115
nodes in surgically resected non–small cell lung can- Naito Y, Goto K, Nagai K, Ishii G, Nishimura M,
cer. Chest 131(4):993–999 Yoshida J, Hishida T, Nishiwaki Y (2010) Vascular
Lerut T, Coosemans W, Decker G, De Leyn P, Ectors N, invasion is a strong prognostic factor after complete
Fieuws S, Moons J, Nafteux P, Van Raemdonck D, resection of node-negative non-small cell lung cancer.
Leuven Collaborative Workgroup for Esophageal Car- Chest 138(6):1411–1417
cinoma (2003) Extracapsular lymph node involvement Nakamura K, Ozaki N, Yamada T, Hata T, Sugimoto S,
is a negative prognostic factor in T3 adenocarcinoma Hikino H, Kanazawa A, Tokuka A, Nagaoka S (2005)
of the distal esophagus and gastroesophageal junction. Evaluation of prognostic significance in extracapsular
J Thorac Cardiovasc Surg 126(4):1121–1128
36 S. V. Tabatabaei et al.
spread of lymph node metastasis in patients with gas- Torre LA, Siegel RL, Jemal A (2016) Lung cancer statis-
tric cancer. Surgery 137(5):511–517 tics. Adv Exp Med Biol 893:1–19
Nakano S, Baba M, Shimada M, Shirao K, Noguchi Y, Ueno H, Mochizuki H, Tamakuma S (1998) Prognostic
Kusano C, Natsugoe S, Yoshinaka H, Fukumoto T, significance of extranodal microscopic foci discontinu-
Aikou T (1999) How the lymph node metastases ous with primary lesion in rectal cancer. Dis Colon
toward cervico–upper mediastinal region affect the Rectum 41(1):55–61
outcome of patients with carcinoma of the thoracic van der Steen S, de Nieuwenhof HP, Massuger L, Bulten J,
esophagus. Jpn J Clin Oncol 29(5):248–251 de Hullu JA (2010) New FIGO staging system of
Paraf F, Fléjou JF, Pignon JP, Fékété F, Potet F (1995) vulvar cancer indeed provides a better reflection of
Surgical pathology of adenocarcinoma arising in prognosis. Gynecol Oncol 119(3):520–525
Barrett’s esophagus. Analysis of 67 cases. Am J Surg Van Schil PE, Balduyck B, De Waele M, Hendriks JM,
Pathol 19(2):183–191 Hertoghs M, Lauwers P (2013) Surgical treatment of
Peters JL, Sutton AJ, Jones DR, Abrams KR, Rushton L early-stage non-small-cell lung cancer. EJC Suppl 11
(2008) Contour–enhanced meta–analysis funnel plots (2):110–122
help distinguish publication bias from other causes of Vansteenkiste JF, De Leyn PR, Deneffe GJ, Stalpaert G,
asymmetry. J Clin Epidemiol 61(10):991–996 Nackaerts KL, Lerut TE, Demedts MG (1997) Survival
Pinto AP, Schlecht NF, Pintos J, Kaiano J, Franco EL, and prognostic factors in resected N2 non–small cell
Crum CP, Villa LL (2004) Prognostic significance of lung cancer: a study of 140 cases. Leuven Lung Cancer
lymph node variables and human papillomavirus DNA Group Ann Thorac Surg 63(5):1441–1450
in invasive vulvar carcinoma. Gynecol Oncol 92 Veronese N, Nottegar A, Pea A, Solmi M, Stubbs B,
(3):856–865 Capelli P, Sergi G, Manzato E, Fassan M, Wood LD,
Puppa G, Maisonneuve P, Sonzogni A, Masullo M, Scarpa A, Luchini C (2016) Prognostic impact and
Capelli P, Chilosi M, Menestrina F, Viale G, Pelosi G implications of extracapsular lymph node involvement
(2007) Pathological assessment of pericolonic tumor in colorectal cancer: a systematic review with meta-
deposits in advanced colonic carcinoma: relevance to analysis. Ann Oncol 27(1):42–48
prognosis and tumor staging. Mod Pathol 20 Wang S, Wu Y, Ou W (2001) Skip metastasis to the
(8):843–855 mediastinal lymph nodes in non-small cell lung cancer.
Riquet M (1996) Capsular rupture of N2 resected lymph Zhonghua Zhong Liu Za Zhi 23(3):259–261
nodes and surgery of bronchial cancer. Rev Mal Respir Wind J, Lagarde SM, Ten Kate FJ, Ubbink DT, Bemelman
13(6):535–537. (Article in French) WA, van Lanschot JJ (2007) A systematic review on
Riquet M, Manac’h D, Saab M, Le Pimpec-Barthes F, the significance of extracapsular lymph node involve-
Dujon A, Debesse B (1995) Factors determining sur- ment in gastrointestinal malignancies. Eur J Surg
vival in resected N2 lung cancer. Eur J Cardiothorac Oncol 33(4):401–408
Surg 9(6):300–304 Wisnivesky JP, Arciniega J, Mhango G, Mandeli J, Halm
Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson EA (2011) Lymph node ratio as a prognostic factor in
GD, Rennie D, Moher D, Becker BJ, Sipe TA, Thacker elderly patients with pathological N1 non-small cell
SB (2000) Meta-analysis of observational studies in lung cancer. Thorax 66(4):287–293
epidemiology: a proposal for reporting. Meta-analysis Zacho A, Fischermann K, Soresen BL (1963) Prognostic
of observational studies in epidemiology (MOOSE) role of breach of lymph node capsule in nodal
group. JAMA 283(15):2008–2012 metastates from gastric carcinoma. Acta Chir Scand
Suzuki H, Fujii T, Asao T, Tsutsumi S, Wada S, Araki K, 125:365–369
Kubo N, Watanabe A, Tsukagoshi M, Kuwano H Zaric B, Stojsic V, Tepavac A, Sarcev T, Zarogoulidis P,
(2014) Extracapsular lymph node involvement is Darwiche K, Tsakiridis K, Karapantzos I, Kesisis G,
associated with colorectal liver metastases and impact Kougioumtzi I, Katsikogiannis N, Machairiotis N,
outcome after hepatectomy for colorectal metastases. Stylianaki A, Foroulis CN, Zarogoulidis K, Perin B
World J Surg 38(8):2079–2088 (2013) Adjuvant chemotherapy and radiotherapy in the
Tachikawa D, Inada S, Kotoh T, Futami K, Arima S, treatment of non–small cell lung cancer (NSCLC). J
Iwashita A (1999) An evaluation of malignancy and Thorac Dis 5(Suppl 4):S371–S377
prognostic factors based on mode of lymph node Zhang B, Zhao L, Yuan Z, Pang Q, Wang P (2016) The
metastasis in esophageal carcinoma. Surg Today 29 influence of the metastasis pattern of mediastinal
(11):1131–1135 lymph nodes on the postoperative radiotherapy’s effi-
Taghavi N, Yazdi I (2015) Prognostic factors of survival cacy for the IIIA-pN2 non-small-cell lung cancer: a
rate in oral squamous cell carcinoma: clinical, histo- retrospective analysis of 220 patients. Onco Targets
logic, genetic and molecular concepts. Arch Iran Med Ther 9:6161–6169
18(5):314–319
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 37–50
https://doi.org/10.1007/5584_2018_264
# Springer Nature Switzerland AG 2018
Published online: 23 September 2018
37
38 W. Krajewski et al.
1 Introduction
2 Methods
The majority of urinary bladder cancers are
non-muscle invasive cancer at the time of diagno- The study was approved by the Ethics Committee
sis. This type of cancer in which tumor resides in of Wroclaw Medical University in Wroclaw,
the mucosa or submucosal layer of the bladder Poland, and was performed in accordance with
contrasts with muscle invasive cancers that infil- the recommendation of the Declaration of
trate or pass through the muscular layer (Babjuk Helsinki for Human Research. Informed consent
et al. 2013). The standard of care in non-muscle was obtained from all participants enrolled into
invasive bladder cancer consists of transurethral the study. Further, the investigating urologist
tumor resection (TURB) with or without addi- explained the purpose of the study to each patient,
tional intravesical therapy. TURB surgery is a the protection of participant’s confidentiality, and
procedure that allows the visualization of the the participant’s freedom to drop out of the study
urethra and the bladder walls, and the removal at any time.
of any suspicious foci. It is commonly performed
under spinal anesthesia with the use of a rigid,
hollow, steel tube (resectoscope) equipped with 2.1 Patients and General Procedures
an optic system and cutting tools, which are
introduced through the urethra. We conducted a prospective observational analy-
Treatment results of a surgical procedure are sis of bladder cancer patients who underwent
evaluated by both clinical outcomes and patient’s TURB surgery. The patients were recruited dur-
reported outcomes; the latter described by the ing 12 months, from October 2016 to October
patient himself and devoid of any interpretation 2017 in the Department of Urology and Urologic
by a third party. The vast majority of available Oncology of Wroclaw Medical University;
studies concerning the patient’s reported Department of General, Oncological and Func-
outcomes in bladder cancer relate to the muscle tional Oncology of Warsaw Medical University,
invasive cancer and its extensive mutilating Warsaw; and the Department of Urology and
effects (Feuerstein and Goenka 2015). Despite Oncologic Urology of Lower Silesian Specialty
that TURB is one of the most commonly Hospital in Wrocław, Poland. Three hundred
performed surgical procedures in urology and sixty-seven male patients scheduled for TURB
the belief of the operation safety, reports were initially considered for the study. Patients
concerning the self-reported outcomes in those with a clinical suspicion of a bladder tumor in
who survive the malignancy are lacking. It is whom resection with cutting current was
noteworthy that patients after TURB surgery performed were included into the study. Those
have a significantly higher readmission rate com- in whom only removal of a tumor with cold
pared with other urological procedures due to a forceps or cold biopsy with fulguration was
relatively high occurrence of postoperative lower performed were excluded. Other exclusion
urinary tract symptoms, bleeding, and bladder criteria were the following: age under 18 years
perforation risk (Rambachan et al. 2014). These old, inability to understand the information about
factors raise the possibility that TURB might have the study or to cooperate with psychological
a considerable impact on the patient’s reported evaluations, and the lack of contact after
Influence of Transurethral Resection of Bladder Cancer on Sexual Function. . . 39
where the higher the value, the better the erectile logistic regression for binary dependent variables.
function and intercourse satisfaction. The IIEF-5 Finally, non-regression taxonomy was used for
score above 21 is considered a normal erectile selected variables to detect the patient types
function; the lower the score, the worse erectile which were compared to psychological and sex-
function. Erectile dysfunction is classified into ual performance outcomes. The computation was
four categories: severe (1–7), moderate (8–11), performed in the R platform (Team 2017).
moderate (12–16), and mild (17–21) (Rhoden
et al. 2002). The tool has been validated in vari-
ous settings (Tang et al. 2014; Deveci et al. 2008). 2.5 Taxonomy
Taxonomic Analysis Dendrograms of taxonomic related to changes in the level IIEF-5 and depres-
distances grouped patients of a given type into sion (one-way ANOVA, Table 5). A rough inter-
aggregates, forming hierarchical super-groups pretation of the clinical patterns of patients is as
(Figs. 1 and 2). The groups were defined follow- follows: a population with a relatively average
ing the increment of the joining heights between BMI 30 kg/m2, aged >45, with nonsmoking
objects. After pruning the dendrogram along the status generates the most negative IIEF-5 changes
observed cut-off value (red line), six groups of after surgery, compared to the remaining types of
patients have been formed. These collections patients (type 1 in the right-hand part of Table 5).
were used for further statistical comparisons. Further, all the mean differences between the
The assigned taxonomic types of patients were patient groups are significant.
42 W. Krajewski et al.
Table 2 Questionnaire results at baseline (pre-TURB) and after transurethral resection of bladder (post-TURB) tumor
Psychometric tool/level Pre-TURB Post-TURB
HADS – anxiety (points) 7.9 3.7 (1–19) 8.8 4.4 (1–21)
HADS – anxiety qualitative borderline (n; %) 80 (31.7) 58 (23.0)
HADS – anxiety qualitative definitive (n; %) 56 (22.2) 83 (32.9)
HADS – depression (points) 5.7 4.1 (0–18) 7.3 4.9 (0–20)
HADS – depression qualitative borderline (n; %) 36 (14.3) 45 (17.9)
HADS – depression qualitative definitive (n; %) 35 (13.9) 61 (24.2)
SSQ sexual satisfaction (points) 29.1 5.3 (13–40) 27.3 5.5 (11–40)
IIEF-5 (points) 16.2 5.0 (2–25) 13.4 5.4 (5–25)
IIEF-5 without ED (n; %) 28 (14.3) 13 (6.7)
IIEF-5 mild ED (n; %) 72 (36.8) 44 (22.4)
IIEF-5 mild to moderate ED (n; %) 61 (31.1) 65 (33.2)
IIEF-5 moderate ED (n; %) 24 (12.2) 44 (22.4)
IIEF-5 severe ED (n; %) 11 (5.6) 30 (15.3)
Data are mean point scores SD (range) or number (%); HADS Hospital Anxiety and Depression Scale, SSQ Sexual
Satisfaction Questionnaire, IIEF-5 simplified International Index of Erectile Function, ED erectile dysfunction. All of the
pre-TURB/post-TURB differences were significant at p < 0.05
Intercept 4.9 0.3 <0.0001 Anxiety Intercept 9.2 0.3 <0.0001 Intercept 0.41 0.23 0.0688
Depression 0.5 0.1 <0.0001 TURB 0.3 0.2 0.0426 SSQ change 0.19 0.06 0.0021
Intercept 11.9 0.9 <0.0001 Intercept 8.1 0.4 <0.0001 Depression Intercept 0.65 0.30 0.0307
IIEF-5 0.3 0.1 <0.0001 Tumor size 0.4 0.1 0.0070 SSQ change 0.03 0.01 0.0002
Depression Intercept 9.1 1.0 <0.0001 Intercept 4.5 0.4 <0.0001 Intercept 1.13 0.19 <0.0001
IIEF-5 0.2 0.1 <0.0001 Depression 0.6 0.04 <0.0001 Cath. Time 0.01 0.00 <0.0001
Intercept 1.4 1.4 0.3340 Intercept 12.3 0.8 <0.0001 Intercept 1.22 0.18 <0.0001
Age 0.1 0.0 <0.0001 Sec. IIEF-5 0.3 0.1 <0.0001 SSQ change 0.16 0.05 0.0016
IIEF-5 Intercept 28.9 2.0 <0.0001 Depression Intercept 1.2 1.7 0.4990 Intercept 1.27 0.14 <0.0001
Age 0.2 0.0 <0.0001 Age 0.1 0.0 <0.0001 Anxiety 0.41 0.05 <0.0001
change
43
(continued)
Table 3 (continued)
44
Fig. 1 Dendrogram of taxonomic analysis of depressive symptoms. Red line presents the cut-off level for statistically
different patients groups
and depressive symptoms were positively results, was clearly influenced by surgery, cathe-
associated with each other and depressive terization, and the occurrence of complication.
symptoms were influenced by the patient’s age That may be related to psychological stress stem-
and presurgery urine abnormalities. That lends ming from the hospitalization and confirmation of
support to previous studies that have a neoplasm, bothersome and demeaning catheter-
demonstrated that urinary infection and lower ization, and physical complaints caused by
urinary tract symptoms (LUTS) might influence LUTS, bleeding, or a post-anesthesia syndrome,
quality of life (Jeong et al. 2015). which often occurs after surgery (Ghoneim and
We demonstrated that TURB procedure O’Hara 2016). However, in some cases, particu-
caused an increase in the level of anxiety and larly in TURB naïve patients, we noticed a sub-
depressive symptoms. The second post-surgery stantial decrease in post-surgery anxiety; the
of HADS survey, which showed worsened reason could be a “falsely” elevated anxiety
46 W. Krajewski et al.
Fig. 2 Dendrogram of taxonomic analysis of erectile function based on the International Index of Erectile Function
(IIEF-5). Red line presents the cut-off level for statistically different patients groups
Acknowledgments This research did not receive any material discursive intra–psychic approach. Cult
specific grants from public, commercial, or nonprofit Health Sex 15:881–895
agencies. Guercio C, Mehta A (2018) Predictors of patient and
partner satisfaction following radical prostatectomy.
Sex Med Rev 6(2):295–301
Conflicts of Interest The authors declare no conflicts of
Hatzimouratidis K (2007) Epidemiology of male sexual
interest in relation to this article.
dysfunction. Am J Mens Health 1:103–125
Hinz A, Brahler E (2011) Normative values for the hospi-
tal anxiety and depression scale (HADS) in the general
German population. J Psychosom Res 71:74–78
References Hinz A, Krauss O, Hauss JP, Höckel M, Kortmann RD,
Stolzenburg JU, Schwarz R (2010) Anxiety and
Asgari MA, Safarinejad MR, Shakhssalim N, depression in cancer patients compared with the gen-
Soleimani M, Shahabi A, Amini E (2013) Sexual func- eral population. Eur J Cancer Care (Engl) 19:522–529
tion after non–nerve–sparing radical cystopros- Ibrahim AI, El-Malik EM, Ismail G, Rashid M, Al Zahrani
tatectomy: a comparison between ileal conduit AB (2002) Risk factors associated with sexual dys-
urinary diversion and orthotopic ileal neobladder sub- function after transurethral resection of the prostate.
stitution. Int Braz J Urol 39:474–483 Ann Saudi Med 22:8–12
Babjuk M, Burger M, Zigeuner R, Shariat SF, van Rhijn Jaidane M, Arfa NB, Hmida W, Hidoussi A, Slama A,
BW, Compérat E, Sylvester RJ, Kaasinen E, Böhle A, Sorba NB, Mosbah F (2010) Effect of transurethral
Palou Redorta J, Rouprêt M, European Association of resection of the prostate on erectile function: a prospec-
Urology (2013) EAU guidelines on non–muscle–inva- tive comparative study. Int J Impot Res 22:146–151
sive urothelial carcinoma of the bladder: update 2013. Jeong SM, Suh B, Jang SH, Jin HS, Kim N, Kwon H,
Eur Urol 64:639–653 Cho B, Park JH (2015) Depression and its severity are
Bowen DK, Butcher MJ, Botchway A, McVary KT (2015) strongly associated with both storage and voiding
Counseling on sexual side effects from TURP. Can J lower urinary tract symptoms independently of pros-
Urol 22:8063–8068 tate volume. J Korean Med Sci 30:1646–1651
Clavien PA, Dindo D (2007) Surgeon’s intuition: is it Kowalkowski MA, Chandrashekar A, Amiel GE, Lerner
enough to assess patients’ surgical risk? World J Surg SP, Wittmann DA, Latini DM, Goltz HH (2014)
31:1909–1911 Examining sexual dysfunction in non–muscle–inva-
Corona G, Lee DM, Forti G, O’Connor DB, Maggi M, sive bladder cancer: results of cross–sectional mixed–
O’Neill TW, Pendleton N, Bartfai G, Boonen S, methods research. Sex Med 2(3):141–151
Casanueva FF, Finn JD, Giwercman A, Han TS, Kulpa M, Kosowicz M, Stypula-Ciuba BJ, Kazalska D
Huhtaniemi IT, Kula K, Lean ME, Punab M, Silman (2014) Anxiety and depression, cognitive coping
AJ, Vanderschueren D, Wu FC, EMAS Study Group strategies, and health locus of control in patients with
(2010) Age–related changes in general and sexual digestive system cancer. Prz Gastroenterol 9:329–335
health in middle–aged and older men: results from Kulpa M, Zietalewicz U, Kosowicz M, Stypula-Ciuba B,
the European Male Ageing Study (EMAS). J Sex Ziolkowska P (2016) Anxiety and depression and cog-
Med 7:1362–1380 nitive coping strategies and health locus of control in
Deveci S, O’Brien K, Ahmed A, Parker M, Guhring P, patients with ovary and uterus cancer during anticancer
Mulhall JP (2008) Can the international index of erec- therapy. Contemp Oncol (Pozn) 20:171–175
tile function distinguish between organic and psycho- Laumann EO, Waite LJ (2008) Sexual dysfunction among
genic erectile function? BJU Int 102:354–356 older adults: prevalence and risk factors from a nation-
Feuerstein MA, Goenka A (2015) Quality of life outcomes ally representative U.S. probability sample of men and
for bladder cancer patients undergoing bladder preser- women 57–85 years of age. J Sex Med 5:2300–2311
vation with radiotherapy. Curr Urol Rep 16:75 Lindau ST, Gavrilova N (2010) Sex, health, and years of
Ficarra V, Righetti R, D’Amico A, Pilloni S, Balzarro M, sexually active life gained due to good health: evidence
Schiavone D, Malossini G, Mobilio G (2000) General from two US population based cross sectional surveys
state of health and psychological well–being in patients of ageing. BMJ 340:c810
after surgery for urological malignant neoplasms. Urol Lindau ST, Schumm LP, Laumann EO, Levinson W,
Int 65:130–134 O’Muircheartaigh CA, Waite LJ (2007) A study of
Gerharz EW, Mansson A, Mansson W (2005) Quality of life sexuality and health among older adults in the United
in patients with bladder cancer. Urol Oncol 23:201–207 States. N Engl J Med 357:762–774
Ghoneim MM, O’Hara MW (2016) Depression and post- Maechler M (2017) Package ‘cluster’. A language and
operative complications: an overview. BMC Surg 16:5 environment for statistical computing. In: Computing
Gilbert E, Ussher JM, Perz J, Wong WK, Hobbs K, Mason RFfS, version 206, Vienna, Austria
C (2013) Men’s experiences of sexuality after cancer: a
50 W. Krajewski et al.
Marczewski E, Steinhaus H (1958) On a certain distance Sighinolfi MC, Micali S, De Stefani S, Mofferdin A,
of sets and the corresponding distance of functions. Ferrari N, Giacometti M, Bianchi G (2007) Bacille
Colloq Math 6:319–327 Calmette–Guérin intravesical instillation and erectile
Mazzola CR, Deveci S, Heck M, Mulhall JP (2012) function: is there a concern? Andrologia 39:51–54
Androgen deprivation therapy before radical prostatec- Sobin LH, Gospodarowicz MK (2010) TNM classification
tomy is associated with poorer postoperative erectile of malignant tumors. In: Wittekind C (ed) International
function outcomes. BJU Int 110:112–116 union against cancer, 7th edn. Wiley–Blackwell,
Mohamed NE, Diefenbach MA, Goltz HH, Lee CT, New York
Latini D, Kowalkowski M, Philips C, Hassan W, Hall Tang Z, Li D, Zhang X, Yi L, Zhu X, Zeng X, Tang Y
SJ (2012) Muscle invasive bladder cancer: from diag- (2014) Comparison of the simplified International
nosis to survivorship. Adv Urol 2012:142135 Index of Erectile Function (IIEF–5) in patients of erec-
Nomejko A, Dolińska-Zygmunt G (2014) The sexual sat- tile dysfunction with different pathophysiologies.
isfaction questionnaire – psychometric poperties. Pol J BMC Urol 14:52
Appl Psychol 12:105–112 Team RC (2017) A language and environment for statisti-
Nomejko A, Dolińska-Zygmunt G, Zdrojewicz Z (2012) cal computing. Vienna, Austria
Quality of life and satisfaction of sexual life – original Thomas HN, Hess R, Thurston RC (2015) Correlates of
results. Seksuologia Polska 10:54–60 sexual activity and satisfaction in midlife and older
Palapattu GS, Haisfield-Wolfe ME, Walker JM, women. Ann Fam Med 13:336–342
BrintzenhofeSzoc K, Trock B, Zabora J, Schoenberg van der Aa MN, Bekker MD, van der Kwast TH, Essink-
M (2004) Assessment of perioperative psychological Bot ML, Steyerberg EW, Zwarthoff EC, Sen FE,
distress in patients undergoing radical cystectomy for Elzevier HW (2009) Sexual function of patients under
bladder cancer. J Urol 172(5 Pt 1):1814–1187 surveillance for bladder cancer. BJU Int 104(1):35–40
Placer J, Salvador C, Planas J, Trilla E, Lorente D, Watrowski R, Rohde A (2014) Validation of the polish
Celma A, López MÁ, Morote J (2015) Effects of version of the hospital anxiety and depression scale in
holmium laser enucleation of the prostate on sexual three populations of gynecologic patients. Arch Med
function. J Endourol 29:332–339 Sci 10:517–524
Rambachan A, Matulewicz RS, Pilecki M, Kim JY, Wittmann D, Northouse L, Foley S, Gilbert S, Wood DP
Kundu SD (2014) Predictors of readmission following Jr, Balon R, Montie JE (2009) The psychosocial
outpatient urological surgery. J Urol 192:183–188 aspects of sexual recovery after prostate cancer treat-
Rhoden EL, Teloken C, Sogari PR, Vargas Souto CA ment. Int J Impot Res 21(22):99–106
(2002) The use of the simplified International Index Wylie K, Kenney G (2010) Sexual dysfunction and the
of Erectile Function (IIEF–5) as a diagnostic tool to ageing male. Maturitas 65:23–27
study the prevalence of erectile dysfunction. Int J Yoshimura K, Utsunomiya N, Ichioka K, Matsui Y,
Impot Res 14:245–250 Terai A, Arai Y (2005) Impact of superficial bladder
Schmidt S, Frances A, Lorente Garin JA, Juanpere N, cancer and transurethral resection on general health–
Lloreta Trull J, Bonfill X, Martinez-Zapata MJ, related quality of life: an SF–36 survey. Urology
Morales Suarez-Varela M, de la Cruz J, Emparanza 65:290–294
JI, Sánchez MJ, Zamora J, Pijoan JI, Alonso J, Ferrer Zak–Łykus A, Nawrat M (2013) Sexual satisfaction: mar-
M (2015) Quality of life in patients with non–muscle– ital and partnership. Family Forum 3:171–186 (Article
invasive bladder cancer: one–year results of a in Polish)
multicentre prospective cohort study. Urol Oncol 33 Zigmond AS, Snaith RP (1983) The hospital anxiety and
(1):19.e7–19.e15 depression scale. Acta Psychiatr Scand 67:361–370
Selvin E, Burnett AL, Platz EA (2007) Prevalence and risk
factors for erectile dysfunction in the US. Am J Med
120:151–157
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 51–62
https://doi.org/10.1007/5584_2018_265
# Springer Nature Switzerland AG 2018
Published online: 29 September 2018
51
52 P. Naumczyk et al.
(information or vocabulary subtest of the The structural brain aging theories seek the
Wechsler Adult Intelligence Scale-Revised origins of pathological aging in noxious processes
(WAIS-R)), and executive functioning (Color in different areas of the brain. Namely, the frontal
Traits Test). The findings unraveled several hypothesis focuses on the damage in the anterior
interdependencies. The higher the intelligence, cortical and subcortical regions (West 1996),
the thicker was the grey matter in nine regions whereas the neurodegenerative hypothesis
of both hemispheres, but also some paradoxi- assumes that aging follows roughly the
cal reversed associations were found in four Alzheimer’s disease pattern of brain damage,
areas; all of them were localized along differ- but at a slower rate (Swerdlow 2007), which
ent sections of the cingulate gyrus in both points to the possible brain pathology in the tem-
hemispheres. An inverse association was poral lobes. The expansion of in vivo neuroimag-
found between crystallized intelligence and ing techniques, such as magnetic resonance
the thickness of the pars opecularis of the imaging (MRI), sheds some light on this theoreti-
right hemisphere. The better the executive cal debate. In the studies on normal healthy aging,
functioning, the thicker was the grey matter brain volume shrinkage has been repeatedly
of a given region. The better the motor perfor- reported (Fjell et al. 2013; Fjell et al. 2009a; Raz
mance, the thicker was the grey matter of the et al. 2010; Raz et al. 2005), with the annual
rostral middle frontal area of the left hemi- percentage change of cortical volume up to
sphere and the lingual gyrus of both 0.8% (Fjell et al. 2014). Yet, current data give
hemispheres. In conclusion, the associations more support to the frontal hypothesis. Although
unraveled demonstrate that the neural cortical atrophy is widespread and present both in
mechanisms underlying healthy aging are frontal and temporal areas (Tamnes et al. 2013;
complex and heterogenic across different cog- Enzinger et al. 2005), it seems that cortical thin-
nitive domains and neuroanatomical regions. ning is more prevalent in the anterior regions
No brain aging theory seems to provide a (Fjell et al. 2009b; Resnick et al. 2003). Nonethe-
suitable interpretative framework for all the less, the frontal lobes dysfunction does not
results. A novel, more integrative approach to explain the diversity of functional changes
the brain aging should be considered. accompanying normal aging. Healthy old people
increasingly seek clinical advice on cognitive
Keywords problems. It has been shown that a range of cog-
Aging process · Cortical thickness · Grey nitive functions deteriorate with age (Josefsson
matter · Healthy aging · Inhibition control · et al. 2012; Finkel et al. 2007). The deterioration
Intelligence · Motor speed primarily concerns working memory (Mattay
et al. 2006; Balota et al. 2000), inhibitory control
(Park and Reuter-Lorenz 2009), the pace of infor-
mation processing (Salthouse 2000), and a long-
1 Introduction
term memory – specifically in the episodic mem-
ory retrieval domain (Reuter-Lorenz and Park
There are four major hypotheses describing the
2010; Grön et al. 2003).
neuronal underpinnings of aging: the frontal
In the research focused on these four areas,
hypothesis, the neurodegenerative hypothesis, the
elderly participants perform worse than the
acceleration in right hemisphere’s aging hypothe-
young ones. Unfortunately, those results are not
sis (the right-hemi theory), and the hemispheric
easily explained in the framework of the frontal
asymmetry reduction in older adults model
hypothesis. Although patients with the anterior
(HAROLD). The first two are grounded in struc-
brain damage experience deficits in the
tural studies, whereas the latter two are rooted in
information processing, inhibition control, and
neuropsychological data on healthy aging in which
working memory, their dysfunction usually far
the subject functions well in everyday life.
Cognitive Predictors of Cortical Thickness in Healthy Aging 53
exceeds what we attribute to the typical pattern of thickness with motor speed, inhibition control,
cognitive aging. To account for the variability in and fluid and crystallized intelligence. Whereas
elderly individuals’ outcomes, the hypothesis of crystallized intelligence is, roughly speaking,
accelerated aging in the right hemisphere’s, a unaffected by normal aging, the other three
right-hemi model has been proposed (Dolcos factors decline with age. Our goal was to compare
et al. 2002). This model suggests an asymmetry the predictive value of all these factors in the
in the age-related vulnerability between the estimation of the grey matter thickness and its
hemispheres. The right-hemi theory was rooted variations across different brain regions in healthy
in the neuropsychological findings of a corre- elderly. In the study, we remained attentive to the
spondence between cognitive performance of patterns supporting or undermining the
healthy aging individuals and patients suffering hypotheses of the brain aging outlined above.
from lateralized brain damage (Taconnat et al.
2007; Uekermann et al. 2006; Royall et al.
2005; Joy et al. 2001). Yet, the studies of struc- 2 Methods
tural brain do not support a notion of asymmetry
in cortical atrophy associated with aging. In the 2.1 Participants
light of such findings, the HAROLD model was
put forward (Cabeza 2002). It postulates that The study protocol was approved by the Ethics
aging is reflected mostly in the functional archi- Committees of the Universities involved in the
tecture of the brain in the form of a reduction of project (permissions NKBBN/296/2017 and
prior specialization of cerebral regions. The 7/2014), and all participants gave written
HAROLD model is, to an extent, consistent with informed consent for the study procedures. The
previous structural findings (see Tamnes et al. participants consisted of 46 pensioners
2013) suggesting that higher-order cortical (37 women and 9 men) over 60 years of age.
systems are especially prone to atrophy in elderly. Individuals with any psychoneurological
Nonetheless, the current understanding of the incidents, brain pathology diagnosed in magnetic
mechanisms underlying age-related brain resonance scan images, dementia symptoms
changes and healthy cognitive aging remains based on the Mini Mental State Examination
inadequate. There are many factors, unaccounted (MMSE), metal bodily implants, or overt neuro-
for, which have a protective or harmful effect on degenerative condition excluded from the study.
the brain’s structure and function during aging. Detailed demographic data of participants are
On the other hand, there are no data providing a depicted in Table 1.
direct contrast between the hypotheses mentioned
above, which would require a difficult task of
collecting and matching both anatomical and neu- 2.2 MRI Acquisition and Processing
ropsychological data. Rather, the two structural or
the two functional models are compared with The participants were examined in two equivalent
each other (Fjell et al. 2009a; Dolcos et al. 3.0 T MRI scanners (Philips Achieva TX;
2002). Further, some previous studies that try to Andover, MA and GE 750 Discovery; Chicago,
link cortical thickness and cognitive functioning IL). The protocol included standard anatomical
in the elderly have failed to refer to the sequences, as brain pathology reference, and a
hypotheses described (Lee et al. 2016; Persson high-resolution T1-weighted sequence, later
et al. 2016; Gautam et al. 2015; Fjell et al. 2006). used for the estimation of cortical thickness. The
The present study seeks to contribute to the T1 sequence parameters for the Philips were:
debate on the theories of aging. We focused on repetition time (TR) ¼ 8.10 ms, echo time
defining neuropsychological correlates of grey (TE) ¼ 3.7 ms, voxel size 111 mm, flip
matter thickness in the elderly. Specifically, we angle 8 ; and for the GE were: TR ¼ 8.19 ms,
examined the relationships between cortical TE ¼ 3.18 ms, voxel size 1x1x1 mm, flip
54 P. Naumczyk et al.
angle ¼ 12 . All anatomical scans were processed thickness of each region was measured by aver-
with FreeSurfer v.6.0 software (Athinoula aging all individual metrics of vertices for a given
A. Martinos Center for Biomedical Imaging at region. The evaluation resulted in 68 cortical
Massachusetts General Hospital, Boston, MA) thickness values for every subject (66 regions
following the recommended processing stream and 2 mean cortical thickness scores).
(http://surfer.nmc.mgh.harvard.edu/), previously
described in detail works (Dale et al. 1999; Fischl
et al. 2004a, b; Fischl et al. 2002; Fischl et al.
2.3 Neuropsychological Assessment
1999). In short, the processing includes transla-
tion into Talairach space and removal of
All of the participants underwent a thorough neu-
non-brain tissue from the images, followed by a
ropsychological assessment that took place about
conversion of the voxel-based images into the
10 days after the MRI session. The neuropsycho-
surface-based 3D representations of triangles
logical protocol included fluid intelligence assess-
(vertices), basing on white and gray matter
ment with the Raven Progressive Matrices
boundaries. The procedure allows for a precise
(RPM-S) standard version, crystallized intelli-
subject-to-subject and subject-to-template
gence assessment with the Wechsler Adult Intel-
registrations upon individual gyral and sulcal
ligence Scale-Revised (WAIS-R) or vocabulary
topology. Anatomical landmarks derived from
subtest of the WAIS-R, and the executive func-
the Desikan-Killiany brain atlas (Desikan et al.
tioning assessment with the Color Traits Test
2006) were fitted into the data resulting in
(CTT).
66 regions of interest, 33 per hemisphere, for
As normative fluid intelligence scores, both
each participant. The atlas is presented in Fig. 1.
percentile and sten scores, were non-normally
The closest distance from the grey matter bound-
distributed (Shapiro-Wilk (46) ¼ 0.787,
ary to the pial boundary at each vertex represents
p < 0.001 and Shapiro-Wilk (46) ¼ 0.901,
cortical thickness at a given vertex. The cortical
Cognitive Predictors of Cortical Thickness in Healthy Aging 55
Fig. 1 The Desikan-Killiany brain atlas used for automated cortical parcellations. (Reproduced from Klein and
Tourville (2012) under the terms of the Creative Commons Attribution License)
p ¼ 0.001, respectively), the normally distributed 1 completion, divided by the time needed for
raw scores of the RPM-S were used for further part 1 completion. The interference score is 0, if
evaluation. For crystallized intelligence, a calcu- there is no difference in the performance between
lated normative score of a given subtest was cho- the test parts. The higher the interference score,
sen to account for differences in the subtests used. the worse the executive functioning, as more time
As for the CTT, three metrics were derived from is needed to account for the increased
the test: the time taken to complete part 1 – inhibition load.
marker of visual attention and motor speed, the
time taken to complete part 2 – marker of execu-
tive control and visual attention, and the interfer-
2.4 Statistical Evaluation
ence index – marker of the flexibility of voluntary
attention, especially over the inhibition processes.
The evaluation focused on the relationship
The last marker was computed as a difference
between the neuropsychological metrics
between the time needed for part 2 and part
obtained, i.e., fluid IQ raw score, crystallized IQ
56 P. Naumczyk et al.
significant linkages in both hemispheres, 4 signif- Age-related effects were notable only in the
icant corrected linkages). The factor that right hemisphere, mostly in the temporal lobe.
accounted for most of them (13 uncorrected and This result, particularly the mean thickness score
2 Bonferroni corrected) was fluid intelligence. prediction, is in line with the right-hemi model
Higher fluid intelligence was associated with described in the introduction. Despite the tempo-
thicker grey matter in nine regions of both ral location, it does not support the neurodegener-
hemispheres, but also some paradoxical reversed ative hypothesis, as the Alzheimer’s pathologies
associations were found in four areas; all of them affect primarily parahippocampal and entorhinal
were localized along different sections of the gyri (Braak and Braak 1991). Still, this is a
cingulate gyrus in both hemispheres. The same finding contradictory to previous studies, where
reversed direction was significant for crystallized no profound asymmetry of age-related effects has
intelligence and the thickness of the pars been found (Lee et al. 2016; Zhou et al. 2013;
opecularis of the right hemisphere. Fjell et al. 2009b; Raz et al. 2005). There are a
There was an inverse relationship between the few explanations possible. The first refers to the
interference index and cortical thickness; the gender inequality present in this study. Our group
higher the interference index, the worse the exec- had 37 women and only 9 men, and there is
utive performance. This relationship indicates research suggesting that age-related brain atrophy
that the better participants’ executive functioning, is less pronounced or even statistically insignifi-
i.e., the smaller the interference index, the thicker cant in females, when compared to peer males
the grey matter of a given region. Thus, executive (Thambisetty et al. 2010; Hutton et al. 2009).
functioning explained underlying differences in Another distinguishing feature is the demo-
variance in five regions of both hemispheres graphic structure. Our participants were all
(one region after Bonferroni correction). Also, elderly, with age ranging from 62 to 86 years,
the better was the participants’ motor perfor- which was a rather low age disparity compared to
mance, the thicker the grey matter of the rostral other studies in which the spread was as wide as
middle frontal area of the left hemisphere and the 18–93 years (Fjell et al. 2009a), 20–79 years (Lee
lingual gyrus of both hemispheres (the left et al. 2016), 22–60 years (Hutton et al. 2009), or
hemisphere’s association was significant after 5–59 years (Zhou et al. 2013). It might be that in
Bonferroni correction). The results, with the our group the difference in brain atrophy was not
corresponding measures of the effect size, are so much apparent within a 24-year span between
provided in Table 2. the eldest and youngest participants. In the main,
our group also was older than some of the others,
which could affect regional variability in the cor-
4 Discussion tical damage observed. An interesting possibility,
however, cannot be ruled out that the method we
This study defined the predictive value of age and employed in this study enabled a direct compari-
cognitive variables for cortical thickness in differ- son of the predictive value of age with variables
ent anatomical regions. Specifically, fluid and associated with cognitive functioning; the latter is
crystallized intelligence, motor speed, and inhibi- evincibly a better predictor of brain atrophy in
tion control were chosen as factors with declining multiple regions.
versus roughly invariable developmental The relationship between fluid and crystallized
trajectories. After correcting for multiple intelligence, and grey matter thickness, on the
comparisons, significant “predictor-cortical thick- other side, we unraveled was twofold. For some
ness” pairs were found in five regions. It is worth associative regions throughout the frontal, parie-
noting that the uncorrected results were in line tal, and temporal lobes of both hemispheres, the
with the ones that passed the rigorous threshold; linkage was positive, meaning that higher intelli-
therefore, those trends will be briefly discussed gence accompanied thicker grey matter, whereas
as well. for others (localized in different sections of
58 P. Naumczyk et al.
cingulate cortex for the fluid intelligence and in 2013), the negative ones raise more interest. It is
the pars opecularis of the right hemisphere for intriguing that such paradoxical associations
crystallized intelligence), it was negative, mean- between the global and/or crystallized intelli-
ing that better intelligence related to thinner grey gence and grey matter thickness, although previ-
matter. As the positive linkages are in line with ously reported, have been left out of discussions
previous studies (Lee et al. 2016; Menary et al. by the authors (Fjell et al. 2006). Although some
Cognitive Predictors of Cortical Thickness in Healthy Aging 59
of our results should be taken with caution, as whereas crystallized intelligence is based on
they were not corrected in the multiple compari- declarative knowledge and the ability to utilize
son analysis, those previous reports of the para- facts (Reber and Reber 2001). Therefore, the for-
doxical associations lend support for interpreting mer has a strong hereditary component (Gray and
them as non-incidental findings. Thompson 2004), while the latter is acquired
There are a variety of possible mechanisms through learning and then shaped by experience
underlying this reversed (paradoxical) relation- (Cattell 1963). Our results suggest that a stronger
ship (higher intelligence – thinner cortical grey prognostic is congenital traits.
matter) found in this and previous studies. One The present results also identified five regions
possible explanation comes from research on the in which the core predictor of grey matter thick-
developing brain. During childhood and adoles- ness was the performance of inhibitory control in
cence, cortical thickness metrics follow a reversed the elderly. Those were localized mostly in asso-
U-shape trajectory. Starting at some age (different ciative regions of the frontal, parietal, and occipi-
between regions and individuals), the thickness of tal lobes. Worse inhibitory control was related to
grey matter continuously decreases greater cerebral atrophy in all the linkages found.
(Khundrakpam et al. 2016). That is interpreted This is a novel finding, as in previous studies on
as a healthy process of synaptic pruning, which neuroanatomy-cognition correspondence in the
plays a key role in regional functional specializa- elderly, the executive domain has either been
tion of the brain (Haartsen et al. 2016). Menary left out of investigation (Raz and Rodrigue
et al. (2013) has provided support to this hypoth- 2006; Persson et al. 2016) or few significant
esis, demonstrating a reversed relationship results have been found (Fjell et al. 2006). Also,
between lexicon scores and grey matter thickness an additional finding in the present study was the
of frontal areas, meaning that the youth with predictive value of the motor speed and visual
extensive vocabulary had thinner cortices. Thus, attention outcome for thickness of three areas of
it may be that thinner grey matter of the cingulate the cerebral cortex. This factor also was not of
cortex and pars opecularis we found, rather than primary interest in previous studies. Lee et al.
atrophy, represents a greater higher degree of (2016) have identified a principal component
specialization of those regions with increased they attributed to motor performance, which is
participants’ intelligence. associated with cortical thickness of multiple
The mechanism above outlined lends support regions, such as inferior temporal, olfactory,
for the HAROLD model. This theory puts empha- parahippocampal, and inferior frontal orbital
sis on the dedifferentiation of previously areas. In contrast, our findings were few, which
specialized brain regions in triggering the cogni- may be a result of using univariate scores of only
tive performance changes observed in the elderly. one test as an indicator in the study.
As fluid and crystallized intelligence are known Nonetheless, the predictive power of the
preservative cognitive factors, their predictive elderly’s executive and motor performance is
value for a greater specialization of certain brain hardly interpretable under any of the theoretical
regions would suggest that these variables play a contexts provided. The results do not support the
key role in restraining the asymmetry reduction frontal hypothesis (as only two linkages in ante-
with aging. From this perspective, it is interesting rior regions were found), nor do they confirm a
that it was fluid, not crystallized, intelligence that higher effect in the temporal areas (as the neuro-
accounted for most of the associations unraveled degenerative model postulates). Also, no obvious
in this study. There are noteworthy differences lateralization was present or absent in the data;
between the two constructs that may contribute thus, it cannot be interpreted as enhanced deterio-
to this result. Fluid intelligence, often synony- ration of the right hemisphere, nor as a reduction
mous with abstract reasoning, is based on the in asymmetry between the hemispheres
ability to solve novel problems creatively, (as HAROLD model assumes). It seems rather
60 P. Naumczyk et al.
Selective increase of cortical thickness in high- Mattay VS, Fera F, Tessitore A, Hariri AR, Berman KF,
performing elderly – structural indices of optimal cog- Das S, Meyer-Lindenberg A, Goldberg TE, Callicott
nitive aging. NeuroImage 29(3):984–994 JH, Weinberger DR (2006) Neurophysiological
Fjell AM, Walhovd KB, Fennema-Notestine C, McEvoy correlates of age-related changes in working memory
LK, Hagler DJ, Holland D, Brewer JB, Dale AM capacity. Neurosci Lett 392(1–2):32–37
(2009a) One-year brain atrophy evident in healthy Menary K, Collins PF, Porter JN, Muetzel R, Olson EA,
aging. J Neurosci 29(48):15223–15231 Kumar V, Steinbach M, Lim KO, Luciana M (2013)
Fjell AM, Westlye LT, Amlien I, Espeseth T, Reinvang I, Associations between cortical thickness and general
Raz N, Agartz I, Salat DH, Greve DN, Fischl B, Dale intelligence in children, adolescents and young adults.
AM, Walhovd KB (2009b) High consistency of Intelligence 41(5):597–606
regional cortical thinning in aging across multiple Park DC, Reuter-Lorenz P (2009) The adaptive brain:
samples. Cereb Cortex 19(9):2001–2012 aging and neurocognitive scaffolding. Ann Rev
Fjell AM, McEvoy L, Holland D, Dale AM, Walhovd KB Psychol 60(1):173–196
(2013) Brain changes in older adults at very low risk Persson N, Ghisletta P, Dahle CL, Bender AR, Yang Y,
for Alzheimer’s disease. J Neurosci 33(19):8237–8242 Yuan P, Daugherty AM, Raz N (2016) Regional brain
Fjell AM, Westlye LT, Grydeland H, Amlien I, shrinkage and change in cognitive performance over
Espeseth T, Reinvang I, Raz N, Dale AM, Walhovd two years: the bidirectional influences of the brain and
KB, Initiative ADN (2014) Accelerating cortical thin- cognitive reserve factors. NeuroImage 126:15–26
ning: unique to dementia or universal in aging? Cereb Raz N, Rodrigue KM (2006) Differential aging of the
Cortex 24(4):919–934 brain: patterns, cognitive correlates and modifiers.
Gautam P, Anstey KJ, Wen W, Sachdev PS, Cherbuin N Neuroscie Biobehav Rev 30(6):730–748
(2015) Cortical gyrification and its relationships with Raz N, Lindenberger U, Rodrigue KM, Kennedy KM,
cortical volume, cortical thickness, and cognitive per- Head D, Williamson A, Dahle C, Gerstorf D, Acker
formance in healthy mid-life adults. Behav Brain Res JD (2005) Regional brain changes in aging healthy
287:331–339 adults: general trends, individual differences and
Gray JR, Thompson M (2004) Neurobiology of intelli- modifiers. Cereb Cortex 15(11):1676–1689
gence: science and ethics. Nat Rev Neurosci 5 Raz N, Ghisletta P, Rodrigue KM, Kennedy KM,
(6):471–482 Lindenberger U (2010) Trajectories of brain aging in
Grön G, Bittner D, Schmitz B, Wunderlich AP, middle-aged and older adults: regional and individual
Tomczak R, Riepe MW (2003) Variability in memory differences. NeuroImage 51(2):501–511
performance in aged healthy individuals: an fMRI Reber AS, Reber ES (2001) Dictionary of psychology.
study. Neurobiol Aging 24(3):453–462 Penguin books, London
Haartsen R, Jones EJH, Johnson MH (2016) Human brain Resnick SM, Pham DL, Kraut MA, Zonderman AB,
development over the early years. Curr Opin Behav Sci Davatzikos C (2003) Longitudinal magnetic resonance
10:149–154 imaging studies of older adults: a shrinking brain. J
Hutton C, Draganski B, Ashburner J, Weiskopf N (2009) Neurosci 23(8):3295–3301
A comparison between voxel-based cortical thickness Reuter-Lorenz PA, Park DC (2010) Human neuroscience
and voxel-based morphometry in normal aging. and the aging mind: a new look at old problems. J
NeuroImage 48(2):371–380 Gerontol B Psychol Sci Soc Sci 65(4):405–415
Josefsson M, de Luna X, Pudas S, Nilsson LG, Nyberg L Royall DR, Palmer R, Chiodo LK, Polk MJ (2005) Normal
(2012) Genetic and lifestyle predictors of 15-year lon- rates of cognitive change in successful aging: the free-
gitudinal change in episodic memory. J Am Geriatr dom house study. J Int Neuropsychol Soc 11
Soc 60(12):2308–2312 (7):899–909
Joy S, Fein D, Kaplan E, Freedman M (2001) Quantifying Salthouse TA (2000) Aging and measures of processing
qualitative features of block design performance speed. Biol Psychol 54(1–3):35–54
among healthy older adults. Arch Clin Neuropsychol Swerdlow RH (2007) Is aging part of Alzheimer’s disease,
16(2):157–170 or is Alzheimer’s disease part of aging? Neurobiol
Khundrakpam BS, Lewis JD, Zhao L, Chouinard-Decorte- Aging 28(10):1465–1480
F, Evans AC (2016) Brain connectivity in normally Taconnat L, Clarys D, Vanneste S, Bouazzaoui B,
developing children and adolescents. NeuroImage Isingrini M (2007) Aging and strategic retrieval in a
134:192–203 cued-recall test: the role of executive functions and
Klein A, Tourville J (2012) 101 labeled brain images and a fluid intelligence. Brain Cogn 64(1):1–6
consistent human cortical labeling protocol. Front Tamnes CK, Walhovd KB, Dale AM, Østby Y,
Neurosci 6:171 Grydeland H, Richardson G, Westlye LT, Roddey
Lee S, Habeck C, Razlighi Q, Salthouse T, Stern Y (2016) JC, Hagler DJ Jr, Due-Tønnessen P, Holland D, Fjell
Selective association between cortical thickness and AM, Alzheimer’s Disease Neuroimaging Initiative
reference abilities in normal aging. NeuroImage (2013) Brain development and aging: overlapping
142:293–300 and unique patterns of change. NeuroImage 68:63–74
62 P. Naumczyk et al.
Thambisetty M, Wan J, Carass A, An Y, Prince JL, West RL (1996) An application of prefrontal cortex func-
Resnick SM (2010) Longitudinal changes in cortical tion theory to cognitive aging. Psychol Bull 120
thickness associated with normal aging. NeuroImage (2):272–292
52(4):1215–1223 Zhou D, Lebel C, Evans A, Beaulieu C (2013) Cortical
Uekermann J, Channon S, Irene Daum I (2006) Humor thickness asymmetry from childhood to older adult-
processing, mentalizing, and executive function in nor- hood. NeuroImage 83:66–74
mal aging. J Int Neuropsychol Soc 12(2):184–191
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 63–73
https://doi.org/10.1007/5584_2018_274
# Springer Nature Switzerland AG 2018
Published online: 5 October 2018
63
64 E. Witkowska-Sędek et al.
is a further difficulty in the exact evaluation of the children. In the 41 strong group of GH-deficient
relationships between the GH/IGF-1 axis and the children, 18 of them were treated with GH for
RANKL/RANK/OPG pathway. 12 months, 14 were treated for 6 months, and in
This study seeks to define serum OPG and 9 children only baseline data were available. All
total sRANKL concentrations in short children, GH-deficient children received recombinant
with different GH secretory status, and to evalu- human GH subcutaneously once daily at bedtime.
ate the associations between the GH/IGF-1 axis Mean GH doses are reported in Table 2. Height,
and the OPG/RANKL pathway in GH-deficient weight, and body mass index (BMI) were
children during GH replacement therapy. expressed as standard deviation scores (SDS)
according to standards of the Institute for Mother
and Child in Warsaw, Poland. BMI was calcu-
2 Methods lated as weight in kilograms divided by height in
square meters (kg/m2). Height was normalized for
2.1 Subjects and General Protocol chronological age, while weight and BMI were
normalized for height-age. Height-age was
This study was conducted at the Department of defined as the age that corresponds to the patient’s
Pediatrics and Endocrinology of the Medical Uni- height when plotted at the 50th percentile on a
versity of Warsaw, Poland, and included 41 growth chart. Baseline height velocity (HV) was
GH-deficient and 20 GH-sufficient short children calculated in all children using data from the
and adolescents recruited prospectively. All the period of 6–18 months before the initiation of
children underwent GH testing in order to evalu- GH treatment or taking baseline anthropometric
ate the GH secretory status. GH deficiency was and biochemical measurements in GH-sufficient
diagnosed according to the following criteria: children. Bone age was evaluated using the
height below the 3rd percentile for sex and age standards of Greulich and Pyle (1969). Pubertal
according to Polish growth charts, height velocity status was assessed according to staging system
less than –1 SD below mean for sex- and of Tanner (1962).
age-matched Polish population, bone age delay,
exclusion of causes of short stature other than GH
deficiency, and the GH level below 10 ng/mL in a 2.2 Biochemical Analysis
test of spontaneous nocturnal growth hormone
secretion (five measurements of GH Serum OPG, total sRANKL, IGF-1, PTH,
concentrations made every 30 min during the 25-hydroxyvitamin D (25(OH)D), calcium and
first 2 h after falling asleep) and in two pharma- phosphorus content, total alkaline phosphatase
cological tests with different stimuli such as clo- (ALP) activity, and 24-h urine calcium and phos-
nidine, insulin, glucagon, or arginine. Children phorus excretion were measured in all children at
with GH release equalling or exceeding 10 ng/ baseline and at 6 and 12 months of GH replace-
mL were considered as GH-sufficient. The maxi- ment therapy in the GH-deficient group. Serum
mum release of GH (GHmax) in a particular content of OPG and total sRANKL were
patient was defined as the highest GH concentra- measured with ELISA (DRG Instruments
tion in any measurement in three tests in GmbH, Marburg, Germany) and those of GH,
GH-deficient children and in any measurement IGF-1, and PTH with an immunoassay using the
in the test or tests performed in GH-sufficient IMMULITE 2000 Xpi analyzer (Siemens,
children. Erlangen, Germany). The serum content of 25
The protocol included baseline anthropometric (OH)D was measured by immunoassay using
and biochemical evaluation of both GH-deficient the ARCHITECT analyzer (Abbott Diagnostics,
and GH-sufficient children and anthropometric Abbott Park, IL). The serum ALP activity, cal-
and biochemical evaluation at 6 and 12 months cium and phosphorus content, and urine calcium
of GH replacement therapy in GH-deficient and phosphorus excretion in 24-h urine samples
66 E. Witkowska-Sędek et al.
were measured by dry chemistry system using the the Spearman correlation analysis for nonpara-
VITROS 5600 chemistry analyzer (Ortho Clinical metric data. The independent predictors of OPG
Diagnostics, Raritan, NJ). IGF-1 concentrations elevation during the first 6 months of GH treat-
were normalized for sex and bone age and were ment were analyzed using backward stepwise
expressed as SDS according to the normative data regression beginning with the three variable
provided by the manufacturer (Siemens regression models. A p-value <0.05 defined sta-
Healthcare Diagnostics Inc., Deerfield, IL). tistically significant differences. The analysis was
performed using a commercial Statistica v13.1
package (StatSoft, Tulsa, OK).
2.3 Statistical Evaluation
Table 1 Anthropometric and biochemical parameters in GH-deficient and GH-sufficient children at baseline
Parameter GH-deficient children GH-sufficient children p
Number of patients 41 20 –
Prepubertal/pubertal status) 31/10 13/7 –
Age (years) 10.3 3.5 10.8 3.4 ns
Bone age (years) 8.2 3.5 9.7 3.7 ns
Height SDS 2.59 0.58 2.60 0.57 ns
Weight SDS for height-age 0.09 0.72 0.36 0.57 ns
BMI SDS for height-age 0.08 1.00 0.45 0.75 ns
HV (cm/year) 5.3 1.4 5.2 1.4 ns
IGF-1 SDS 0.37 1.04 0.20 1.41 ns
OPG (pmol/L) 5.6 (3.6–10.1) 9.4 (4.8–11.5) ns
Total sRANKL (pmol/L) 226.3 (164.5–317.3) 243.8 (179.9–337.2) ns
Total ALP (U/L) 202.0 (171.0–231.0) 193.5 (160.0–236.5) ns
PTH (pg/mL) 20.7 (12.2–24.1) 17.6 (11.1–22.9) ns
25(OH)D (ng/mL) 25.8 (22.2–34.0) 27.6 (24.3–34.6) ns
Serum Ca (mg/dL) 10.0 (9.8–10.3) 10.1 (10.0–10.2) ns
Serum P (mg/dL) 5.0 0.46 5.1 0.53 ns
Urine Ca excretion (mg/kg/24 h) 2.9 (1.5–4.0) 2.5 (1.0–4.2) ns
Urine P excretion (mg/kg/24 h) 16.8 (14.0–20.1) 19.5 (16.1–23.1) ns
Data are presented as means SD or medians with interquartile ranges (IR) as appropriate; ns not significant, GH growth
hormone, SDS standard deviation score, BMI body mass index, HV height velocity, IGF-1 insulin-like growth factor-1,
OPG osteoprotegerin, sRANKL receptor activator of nuclear factor kappa B ligand, ALP alkaline phosphatase, PTH
parathormone, 25(OH)D 25-hydroxyvitamin D, Ca calcium, P phosphorus
Osteoprotegerin, Receptor Activator of Nuclear Factor Kappa B Ligand, and. . . 67
(p < 0.001), which led to an increase in height significantly higher than that at baseline after
SDS observed after the first 6 months of treatment 12 months of GH treatment (p < 0.05), but its
(p < 0.001). Weight SDS and BMI SDS for 24-h urine excretion did not change significantly
height-age did not change significantly during during therapy. Nor did serum phosphorus con-
the first 12 months of GH treatment. The mean tent change during GH treatment. The anthropo-
IGF-1 SDS increased significantly after 6 months metric and biochemical parameters at baseline
of GH treatment compared to the baseline value and at 6 and 12 months of GH replacement ther-
(p < 0.001) but then declined back toward the apy are presented in Table 2.
baseline value after another 6 months of treat-
ment, remaining insignificantly higher than that
at baseline. The OPG content and total ALP
3.2 Associations Between OPG or
activity increased significantly after the first
Total sRANKL
6 months of GH treatment (p < 0.01 and
and Anthropometric
p < 0.05, respectively) and then reached a plateau,
and Biochemical Parameters
remaining significantly higher than that at base-
line after 12 months of treatment (p < 0.001 and
Baseline associations between OPG or total
p < 0.05, respectively). Total sRANKL, PTH, and
sRANKL and the anthropometric and biochemi-
25(OH)D did not change significantly during GH
cal parameters investigated, taken as independent
treatment. Serum calcium content was
variables, in GH-deficient and GH-sufficient
Table 2 Anthropometric and biochemical parameters at baseline and during GH treatment in GH-deficient children
Parameter Baseline 6 months 12 months pa pb pc
Number of patients 41 32 18 – – –
Age (years) 10.3 3.5 11.3 3.4 11.7 3.6 – – –
Bone age (years) 8.2 3.5 – 9.9 3.5 – <0.001 –
Height SDS 2.59 0.58 2.23 0.59 2.01 0.64 <0.001 <0.001 <0.01
Weight SDS for height-age 0.09 0.72 0.17 0.72 0.06 0.72 ns ns ns
BMI SDS for height-age 0.08 1.00 0.18 0.98 0.03 1.01 ns ns ns
HV (cm/year) 5.3 1.4 – 9.8 1.03 – <0.001 –
GH dose (mg/kg/week) – 0.19 0.01 0.19 0.01 – – –
IGF-1 SDS 0.37 1.04 2.02 2.51 0.94 1.06 <0.001 ns <0.05
OPG (pmol/L) 5.55 7.57 9.70 <0.01 <0.001 ns
(3.56–10.05) (4.31–12.16) (6.36–12.42)
Total sRANKL (pmol/L) 226.3 260.9 316.4 ns ns ns
(164.5–317.3) (155.4–354.0) (129.3–511.4)
Total ALP (U/L) 202.0 244.0 264.0 <0.05 <0.05 ns
(171.0–231.0) (212.0–291.0) (225.0–282.0)
PTH (pg/mL) 20.7 (12.2–24.1) 19.2 (11.2–30.3) 16.6 (12.6–22.5) ns ns ns
25(OH)D (ng/mL) 25.8 (22.2–34.0) 27.3 (23.1–31.2) 28.0 (23.1–31.8) ns ns ns
Serum Ca (mg/dL) 10.0 (9.8–10.3) 10.2 (9.9–10.3) 10.2 (10.0–10.5) ns <0.05 ns
Serum P (mg/dL) 5.00 0.46 5.30 0.36 5.40 0.80 ns ns ns
Urine Ca excretion (mg/kg/ 2.9 (1.5–4.0) 2.3 (2.0–3.2) 2.0 (1.2–4.4) ns ns ns
24 h)
Urine P excretion (mg/kg/ 16.8 (14.0–20.1) 20.1 (14.5–26.1) 17.8 (14.3–22.3) ns ns ns
24 h)
Data are presented as means SD or medians with interquartile ranges (IR) as appropriate; pa p-value baseline vs. 6 months,
pb p-value baseline vs. 12 months, pc p-value 6 vs. 12 months, ns not significant, GH growth hormone, SDS standard
deviation score, BMI body mass index, HV height velocity, IGF-1 insulin-like growth factor-1, OPG osteoprotegerin,
sRANKL receptor activator of nuclear factor kappa B ligand, ALP alkaline phosphatase, PTH parathormone, 25(OH)D
25-hydroxyvitamin D, Ca calcium, P phosphorus
68 E. Witkowska-Sędek et al.
Table 3 Spearman’s correlation coefficients (r) between baseline serum OPG (pmol/L) or total sRANKL (pmol/L)
concentrations and baseline anthropometric and biochemical parameters in GH-deficient and GH-sufficient children
GH-deficient children GH-sufficient children
Parameter OPG Total sRANKL OPG Total sRANKL
Age (years) ns ns ns 0.45 (<0.05)
Bone age (years) ns ns ns 0.46 (<0.05)
Height SDS ns ns ns ns
Weight SDS for height-age 0.33 (<0.05) 0.46 (<0.01) ns ns
BMI SDS for height-age ns 0.41 (<0.01) ns 0.47 (<0.05)
HV (cm/year) ns ns ns ns
GHmax (ng/mL) ns ns ns ns
IGF-1 SDS ns ns 0.55 (< 0.05) ns
Total ALP (U/L) ns ns ns ns
PTH (pg/mL) ns ns ns ns
25(OH)D (ng/mL) ns 0.34 (<0.05) ns ns
Serum Ca (mg/dL) ns ns 0.46 (<0.05) ns
Serum P (mg/dL) ns ns ns ns
Urine Ca excretion (mg/kg/24 h) ns ns ns ns
Urine P excretion (mg/kg/24 h) ns ns ns ns
GHmax maximum growth hormone release, SDS standard deviation score, BMI body mass index, HV height velocity,
IGF-1 insulin-like growth factor-1, OPG osteoprotegerin, sRANKL receptor activator of nuclear factor kappa B ligand,
ALP total alkaline phosphatase, PTH parathormone, 25(OH)D 25-hydroxyvitamin D, Ca calcium, P phosphorus, ns not
significant
children, are presented in Table 3. There were no with baseline height SDS (r ¼ 0.49; p < 0.01),
significant associations between either maximum weight SDS and BMI SDS for height-age
GH release, baseline height SDS or HV and either (r ¼ 0.54; p < 0.01 and r ¼ 0.58; p < 0.001,
baseline OPG or total sRANKL levels in GH- respectively), and with baseline IGF-1 SDS
deficient and GH-sufficient children. In (r ¼ 0.46; p < 0.01). Total sRANKL content did
GH-deficient children, the strongest associations not change significantly during GH treatment.
were found between the total sRANKL content Likewise, we did not find any associations of
and the nutritional parameters such as weight significance between OPG concentrations after
SDS for height-age and BMI SDS for height-age the first 6 months of GH treatment and any of
(r ¼ 0.46; p < 0.01 and r ¼ 0.41; p < 0.01, the anthropometric or biochemical parameters
respectively). In GH-sufficient children, OPG measured after the first 6 months of GH treatment.
content was adversely associated with IGF-1 The increment in IGF-1 during the first 6 months
SDS (r ¼ 0.55, p < 0.05), while total sRANKL of therapy (ΔIGF-16-month-baseline) and GH doses
was adversely associated with age (r ¼ 0.45; administered in that period were not associated
p < 0.05), bone age (r ¼ 0.46; p < 0.05), and with OPG concentrations after 6 months of ther-
BMI SDS for height-age (r ¼ 0.47; p < 0.05). apy either. ΔOPG6-month-baseline was associated
There were no statistically significant with height SDS (r ¼ 0.46; p < 0.01), weight
associations between OPG or total sRANKL SDS, and BMI SDS for height-age (r ¼ 0.61;
levels at 6 months of GH treatment and baseline p < 0.001 and r ¼ 0.50; p < 0.01, respectively)
anthropometric and biochemical parametres. Fur- and total ALP activity (r ¼ 0.39; p < 0.05) and
ther analysis revealed that the increment in OPG was adversely associated with 25(OH)D
concentration in the first 6 months of treatment concentrations (r ¼ 0.45; p < 0.05) evaluated
(ΔOPG6-month-baseline) was adversely associated after the first 6 months of GH treatment. Total
with GHmax (r ¼ 0.41; p < 0.05) and positively sRANKL concentration after the first 6 months of
Osteoprotegerin, Receptor Activator of Nuclear Factor Kappa B Ligand, and. . . 69
GH treatment was associated only with total ALP months of therapy. Backward stepwise regression
activity at 6 months of therapy (r ¼ 0.42; analysis was not performed at 12 months of GH
p < 0.05). treatment because there were too few patients
To evaluate the relationship between the who could be included in such analysis.
ΔOPG6-month-baseline, considered as a dependent
variable, and the independent study variables,
measured at baseline (models 1, 2, and 3) and at 4 Discussion
6 months of GH treatment (models 4, 5, and 6),
backward stepwise regression analysis was used. Several authors have investigated the relationships
Regression models included the following inde- between the GH/IGF-1 axis and the RANKL/
pendent variables: model 1, baseline height SDS, RANK/OPG pathway in experimental and clinical
weight SDS, and BMI SDS for height-age; model studies (Xia et al. 2015; Guerra–Menéndez et al.
2, GH max, baseline bone age, and IGF-1 SDS for 2013; Flint et al. 2009; Wang et al. 2006; Lanzi
bone age; model 3, baseline total ALP, PTH, and et al. 2003; Ueland et al. 2003), but such studies in
25(OH)D; model 4, height SDS, weight SDS, and GH-deficient patients, especially in children and
BMI SDS for height-age, all at 6 months of GH adolescents, are scarce. IGF-1 influences OPG and
treatment; model 5, IGF-1 SDS for bone age, RANKL concentrations in a dose- and time-
ΔIGF-16-month-baseline, and GH dose administered dependent manner in vitro and in vivo, which
in the first 6 months of treatment; model 6, total confirms its role as a coupling agent in direct
ALP, PTH, and 25(OH)D, all at 6 months of GH activation of bone formation and indirect, through
treatment. The variables, baseline BMI SDS for OPG/RANKL, activation of bone resorption
height-age (β ¼0.42; p < 0.05), baseline ALP (Rubin et al. 2002). It is hypothesized that GH
activity (β ¼0.36; p < 0.05), weight SDS for induces OPG production and is able to modulate
height-age at 6 months of GH treatment (β bone remodeling by directly influencing
¼1.86; p < 0.01), and total ALP activity at osteoblast-osteoclast cross talk (Mrak et al.
6 months of GH treatment (β ¼ 0.48; p < 0.01), 2007). The results of studies evaluating OPG and
were identified as independent predictors of RANKL concentrations in children and
ΔOPG6-month-baseline in models 1, 3, 4, and adolescents are divergent. Further multicenter
6, respectively. Neither GH max, nor IGF-1 studies, including large cohorts of children and
SDS or increase in IGF-1 SDS after the initiation adolescents, are needed to establish valuable refer-
of GH treatment, nor GH dose, nor bone age was ence data for pediatric population and to foster the
found to be significant independent predictors of clinical use of the RANKL/RANK/OPG pathway.
OPG elevation after the initiation of therapy. Buzi et al. (2004) have found that the OPG level in
After the first 12 months of GH treatment, we healthy children is gender independent but changes
failed to find any significant associations between with age and is significantly higher in the first
the OPG content or ΔOPG12-month-baseline and any 4 years of life than at 4–14 years of age, while
of the anthropometric and biochemical RANKL content slightly declines with age. Those
parameters. GH doses administered in the first authors have also reported that serum OPG content
12 months of replacement therapy were adversely in children affected by Turner’s syndrome
associated with OPG at 12 months (r ¼ 0.60; (TS) and early or precocious puberty (PP) is com-
p < 0.05). We also noticed that there was a ten- parable to those in healthy peers, while that in
dency for greater elevation in OPG in the first children with rheumatoid arthritis (RA) is signifi-
12 months of therapy in children treated with cantly higher than in controls. Simultaneously,
lower GH doses, but that result was statistically they have also found that serum RANKL does
insignificant (r ¼ 0.49; p ¼ 0.06). Total not differ significantly between TS-, PP-, and
sRANKL content at 12 months of GH treatment RA-afflicted children and age-matched healthy
was not associated with any of the anthropometric population. On the other hand, Wasilewska et al.
and biochemical parameters evaluated at 12 (2009) have suggested that in healthy children and
70 E. Witkowska-Sędek et al.
and nutritional status independently predicted the disorders in GH-deficient adults. They demon-
elevation of OPG after the initiation of GH treat- strate that plasma OPG content is higher in
ment. We also found that GH doses were GH-deficient adults than in healthy controls and
adversely associated with the increase in OPG in is associated with age, BMI, waist and hip circum-
the first year of replacement therapy. The plausi- ference, waist-hip ratio, and the content of fasting
ble association between GH administration and insulin, total cholesterol, triglycerides, high sensi-
OPG content ought to be confirmed in a larger tivity C-reactive protein, and interleukin-6. The
group of patients. Meazza et al. (2014) have also authors suggest that OPG could be a potential
investigated the effects of 12-month GH treat- therapeutic target in cardiovascular disorders.
ment on OPG content and on various metabolic Those results suggest that the RANKL/RANK/
parameters, including adipokine profile, in OPG pathway is associated not only with bone
GH-deficient children, and report, in opposition turnover but also with metabolic status and could
to the present study, a lack of significant changes be used as a marker of impaired metabolic status
after the initiation of GH treatment. Those authors leading to increased cardiovascular risk. Although
have found that height SDS, both at baseline and those observations have been made in
after 12 months of GH treatment, was adversely GH-deficient adults, it is possible that the same
associated with the OPG content in GH-deficient associations could also apply for GH-deficient
children. They also found an adverse association children and adolescents. Erol et al. (2016) have
between the baseline OPG and IGF-1 content found that OPG content is lower in obese than
before treatment. The discrepancy between the nonobese children, but they do not find any
present findings and those above outlined points associations between OPG and insulin resistance.
to the need of further exploration of OPG content In conclusion, we found that serum OPG or
changes due to GH treatment. sRANKL was independent from GH secretory
The results of studies in adult populations par- status and IGF-1 content in short children. After
tially correspond to those reported in children. the initiation of GH replacement treatment, serum
Ueland et al. (2001) have reported that serum OPG content increased significantly, while total
OPG content does not depend on GH secretory sRANKL remained unchanged. The increase in
status and does not differ between GH-deficient OPG during GH treatment was positively
and healthy adults. Lanzi et al. (2003), who associated with total ALP activity and nutritional
investigated the effects of 6-month-long GH treat- status in GH-deficient children. The exploration of
ment in both childhood and adult-onset GH defi- mutual relationships between the RANKL/RANK/
ciency, have reported a significant increase in OPG pathway and both bone turnover and
plasma OPG content after the first 6 months of glucolipid metabolism could lead to new avenues
GH treatment, which is accompanied by signifi- of clinically applicable research concerning GH
cant increases in total serum calcium and urinary deficiency and GH replacement treatment.
calcium excretion with unchanged PTH and 25
(OH)D levels. In the same study, the authors also Conflicts of Interest The authors declare no conflicts of
evaluate the associations between the OPG level interest in relation to this article.
during GH treatment and two bone turnover
markers: osteocalcin, considered a bone formation Ethical Approval All procedures performed in this study
were in accordance with the ethical standards of the insti-
marker, and urine deoxypyridinoline excretion,
tutional research committee and with the 1964 Helsinki
considered a bone resorption marker. They have Declaration and its later amendments. The study protocol
reported an inverse association of OPG increase was approved by the Bioethics Committee of the Medical
during GH treatment with the increase in University of Warsaw (permit: KB/93/2016).
osteocalcin and with deoxypyridinoline excretion.
Xia et al. (2015) have revealed that OPG, despite Informed Consent Informed written consent was
obtained from all individual participants included in the
its associations with bone metabolism, might also
study or from their parents and legal guardians.
play a role in the development of cardiovascular
72 E. Witkowska-Sędek et al.
Rubin J, Ackert–Bicknell CL, Zhu L, Fan X, Murphy TC, members of the IGF–1 family in cortical bone from
Nanes MS, Marcus R, Holloway L, Beamer WG, women with adult onset growth hormone deficiency–
Rosen CJ (2002) IGF–I regulates osteoprotegerin relationship with bone turn–over. Bone 33:638–645
(OPG), and receptor activator of nuclear factor–Kb Vega D, Maalouf NM, Sakhaee K (2007) The role of
ligand in vitro and OPG in vivo. J Clin Endocrinol receptor activator of nuclear factor–κB (RANK)/
Metab 87(9):4273–4279 RANK Ligand/Osteoprotegerin: clinical implications.
Simonet WS, Lacey DL, Dunstan CR, Kelley M, Chang J Clin Endocrinol Metab 92:4514–4521
MS, Lüthy R, Nguyen HQ, Wooden S, Bennett L, Wada T, Nakashima T, Hiroshi N, Penninger JM (2006)
Boone T, Shimamoto G, DeRose M, Elliott R, RANKL–RANK signaling in osteoclastogenesis and
Colombero A, Tan HL, Trail G, Sullivan J, Davy E, bone disease. Trends Mol Med 12:17–25
Bucay N, Renshaw–Gegg L, Hughes TM, Hill D, Wang Y, Nishida S, Elalieh HZ, Long RK, Halloran BP,
Pattison W, Campbell P, Sander S, Van G, Tarpley J, Bikle D (2006) Role of IGF–1 signaling in regulating
Derby P, Lee R, Boyle WJ (1997) Osteoprotegerin: a osteoclastogenesis. J Bone Miner Res 21:1350–1358
novel secreted protein involved in the regulation of Wasilewska A, Rybi–Szuminska AA, Zoch–Zwierz W
bone density. Cell 89(2):309–319 (2009) Serum ostoprotegerin (OPG) and receptor acti-
Siomou E, Challa A, Printza N, Giapros V, vator of nuclear factor kappa B (RANKL) in healthy
Petropopulou F, Mitsioni A, Papachristou F, Stefanidis children and adolescents. J Pediatr Endocrinol Metab
CJ (2011) Serum osteoprotegerin, RANKL and fibro- 22(12):1099–1104
blast growth factor–23 in children with chronic kidney Witkowska–Sędek E, Kucharska A, Rumińska M, Pyrżak
disease. Pediatr Nephrol 26:1105–1114 B (2016) Relationship between 25(OH)D and IGF–I in
Szulc P, Hofbauer LC, Delmas PD (2001) OPG serum children and adolescents with growth hormone defi-
levels in men: correlation with age, estrogen and tes- ciency. Adv Exp Med Biol 912:43–49
tosterone status. J Clin Endocrinol Metab Witkowska–Sędek E, Stelmaszczyk–Emmel A,
86:3162–3165 Kucharska A, Demkow U, Pyrżak B (2017) Associa-
Takayanagi H (2007) Osteoimmunology: shared tion between vitamin D and carboxy–terminal cross–
mechanisms and crosstalk between the immune and linked telopeptide of type I collagen in children during
bone systems. Nat Rev Immunol 7:292–304 growth hormone replacement therapy. Adv Exp Med
Tanner JM (1962) Growth at adolescence. Blackwell, Biol 1047:53–60
Oxford Witkowska–Sędek E, Stelmaszczyk–Emmel A,
Theoleyre S, Wittrant Y, Tat SK, Fortun Y, Redini F, Majcher A, Demkow U, Pyrżak B (2018) The relation-
Heymann D (2004) The molecular triad OPG/RANK/ ship between alkaline phosphatase and bone alkaline
RANKL: involvement in the orchestration of patho- phosphatase activity and the growth hormone/insulin–
physiological bone remodelling. Cytokine Growth like growth factor–1 axis and vitamin D status in
Factor Rev 15:457–475 children with growth hormone deficiency. Acta
Ueland T, Bollerslev J (2005) Bone metabolism and Biochim Pol 65(2):269–275
growth hormone deficiency. In: Jorgensen JOL, Xia J, Ren W, Zheng X, Liu C, Li J, Chen T, Li X,
Christiansen JS (eds) Frontiers in hormone research. Wang L, Hu Y (2015) Correlation of increased plasma
Karger, Basel osteoprotegerin and cardiovascular risk factors in
Ueland T, Bollerslev J, Godang K, Muller FS, Aukrust P patients with adult growth hormone deficiency. Int J
(2001) Increased serum osteoprotegerin in disorders Clin Exp Med 8(3):3184–3192
characterized by persistence autoimmune activation Yasuda H, Shima N, Nakagawa N, Yamaguchi K,
or glucocorticoid excess – possible role in bone Kinosaki M, Mochizuki S, Tomoyasu A, Yano K,
homeostasis. Eur J Endocrinol 145:685–690 Goto M, Murakami A, Tsuda E, Morinaga T,
Ueland T, Bollerslev J, Flyvbjerg A, Hansen TB, Vahl N, Higashio K, Udagawa N, Takahashi N, Suda T
Mosekilde L (2002) Effects of 12 months of GH treat- (1998) Osteoclast differentiation factor is a ligand for
ment on cortical and trabecular bone content of IGFs osteoprotegerin/osteoclastogenesis–inhibitory factor
and OPG in Adults with acquired GH deficiency: a and is identical to TRANCE/RANKL. Proc Natl
double–blind, randomized, placebo–controlled study. Acad Sci U S A 95(7):3597–3602
J Clin Endocrinol Metab 87:2760–2763 Yates L, Górecki DC (2006) The nuclear factor–kappa B
Ueland T, Odgren PR, Yndestad A, Godang K, (NF–κB): from a versatile transcription factor to a
Schreiner T, Marks SC Jr, Bollerslev J (2003) Growth ubiquitous therapeutic target. Acta Biochim Pol 53
hormone substitution increases gene expression of (4):651–662
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 75–79
https://doi.org/10.1007/5584_2018_266
# Springer Nature Switzerland AG 2018
Published online: 23 September 2018
Inhibition of Cross-Reactive
Carbohydrate Determinants in Allergy
Diagnostics
75
76 M. Grzywnowicz et al.
makes CCD different from other carbohydrate positive levels of anti-CCD antibodies. To this
residues is their antigenicity and ability to stimu- end, we compared the asIgE content in relation
late the production of allergen specific immuno- to 30 allergens using multiparametric
globulin E (asIgE). The term “cross-reactive” Polycheck® tests with and without a CCD
determinants was first used more than 30 years inhibitor.
ago (Aalberse et al. 1981). Subsequent research
has demonstrated that anti-CCD antibodies are of
no clinical relevance, since their presence is not
2 Methods
associated with any symptoms (Mertens et al.
2010; Muller et al. 2009; Mari et al. 2008; Ebo
The study was approved by the institutional
et al. 2004; Mari 2002). While basophil activation
Ethics Committee and was performed according
tests have demonstrated bioactivity of anti-CCD
to the guidelines of the Helsinki Declaration for
asIgE, no connection has been substantiated to
Human Research. There were 24 children of
clinical symptoms in patients (Mertens et al.
either genders (12 boys and 12 girls) aged 3–10
2010; Foetisch et al. 2003).
(average 5 years) enrolled into the study. They
Despite the seeming clinical irrelevance, anti-
were hospitalized with chronic upper and lower
CCD antibodies pose a serious problem in case of
respiratory tract symptoms and skin symptoms in
the in vitro asIgE diagnostics with respect to the
the Regional Branch of Warsaw Institute of
allergens in which they appear, making it harder to
Tuberculosis and Lung Diseases in the mountain-
correctly diagnose the allergy source. That is
ous city of Rabka-Zdrój from July to August
caused by the structural heterogeneity of CDD
2017. The goal of hospitalization was to conduct
and their resemblance in various allergen sources.
the diagnostics procedures and to determine ther-
The most important type of CCD, having a
apy. Seventeen children were diagnosed with
sequence of N-glycan with alpha(1,3)-fructose
asthma, associated with allergic rhinitis in 8 of
attached to the polypeptide chain through nitrogen
them, and 7 children were diagnosed with moder-
atom in the asparagine radical, is present in both
ately severe atopic dermatitis. Anti-CCD
plant allergens, including inhalable allergens and
antibodies were found in the sera of all children
various food allergens, hymenoptera venom, and
in a range of 0.43–25 kU/L (median 1.5 kU/L).
latex (Altmann 2016). Therefore, when an exam-
The Polycheck® Atopic 30-I (Biocheck
ined person produces anti-CCD antibodies, they
GmbH, Münster, Germany), an allergy test rou-
can bind with allergens which contain CDDs dur-
tinely used in allergy screening, was the diagnos-
ing the test, giving a false positive or falsely ele-
tic method applied. This test is a multiparametric,
vated result in relation to the actual asIgE content.
quantitative test with a lower limit of detection of
Taking into account that the prevalence of anti-
0.15 kU/L. As recommended by the producer, it
CCD antibodies in people subjected to the in vitro
was assumed that positive results start from the
asIgE testing can reach 25%, they pose a serious
level of 0.35 kU/L. The allergens included in the
diagnostic problem in the allergy testing (Muller
assay were the following: cow’s milk, alpha-
et al. 2009; Ebo et al. 2004; Mari et al. 2008). The
lactalbumin, beta-lactoglobulin, casein, egg
prevalence of anti-CCD antibodies is even higher
white, egg yolk, codfish, peanut, cacao, soybean,
among young people, amounting to 35%.
apple, carrot, tomato, flour mix, chicken meat,
One way to circumvent or eliminate the prob-
citrus mix, rice, grass mix, rye, dog epithelia, cat
lem posed by anti-CCD antibodies in the in vitro
epithelia, Cladosporium herbarum, Alternaria
asIgE diagnostics is to perform the component-
alternata, Aspergillus fumigatus,
resolved diagnostics, based on the antigens free
Dermatophagoides pteronyssinus, D. farinae,
from CCD or utilizing a CCD inhibitor. The
hazel pollen, birch pollen, and mugwort pollen.
premise behind the present study was to evaluate
The test also includes a CCD marker, pineapple
the influence of a CCD inhibitor on the asIgE
bromelain.
content in individuals previously diagnosed with
Inhibition of Cross-Reactive Carbohydrate Determinants in Allergy Diagnostics 77
A Polycheck® CCD inhibitor was used Table 1 Anti-CCD antibody content assayed without and
according to the producer’s recommendations, with the CCD inhibitor (kU/L)
added at a ratio of 1:50 to each tested serum. CCD inhibitor Difference
This synthetic inhibitor is composed of polyva- Sample Without With %
lent CCD structures of bromelain, with up to four 1 25.00 6.30 75
amino acids left after the proteolysis and purifica- 2 0.47 0.22 53
tion processes, which guarantees the lack of pep- 3 0.78 0.41 47
tide epitopes. The glycopeptide obtained in such a 4 1.00 0.56 44
way is attached to a protein carrier, which is 5 2.40 1.40 42
6 1.90 1.20 37
highly purified human serum albumin.
7 8.90 6.10 31
Data were presented in the tabular form as two
8 0.89 0.62 30
nominal variables (sample number and the 9 0.80 0.59 26
absence or presence of CDD inhibitor) and one 10 0.82 0.63 23
measurement variable (content of anti-CCD anti- 11 0.76 0.60 21
body). Differences in the content of anti-CCD 12 2.70 2.30 15
antibody assayed without and with the addition 13 0.56 0.49 13
of a CCD inhibitor, expressed as the percentage 14 2.80 2.50 11
of the former, were statistically analyzed using 15 0.87 0.80 8
the Wilcoxon signed-rank test. A p-value <0.05 16 3.40 3.30 3
defined a statistically significant difference. Sta- 17 1.50 1.50 0
tistical analysis was performed using a commer- 18 2.00 2.00 0
19 2.20 2.20 0
cial GraphPad Prism v7 package (San Diego,
20 1.50 1.60 +7
CA).
21 1.00 1.10 +10
22 0.43 0.48 +12
23 2.40 2.90 +21
3 Results 24 6.30 7.90 +25
antibodies to plant allergens, except for the the limit of detection, completely mitigating the
antibodies related to tree and grass pollens false positive. On average, the level of the entire
which remained at a level exceeding 100 kU/L. anti-CCD antibody content was lowered by 72%,
despite a 23% decrease in the CCD marker. That
clearly illustrates an advantage of the use of a
CCD inhibitor, instead of testing different levels
4 Discussion
of a CCD marker. Unfortunately, in the patients
investigated in this study, despite positive values
A CCD inhibitor is used in the asIgE detection
for CCD, only one sample displayed numerous
test to limit false positive reactions, by competi-
positive results for asIgE. It is worth emphasizing
tively binding anti-CCD antibodies in the sample
that the use of multiparametric tests makes it
assayed before they become bound with allergens
much easier for a physician to identify the pres-
containing a CCD residue. In general, CCD inhib-
ence of anti-CCD IgE antibodies than when
itor, as opposed to allergen-extract molecules, is
interpreting individual results (Altmann 2016;
composed of a mix of molecules that are unasso-
Holzweber et al. 2013).
ciated with the solid phase of the test, which is
The CCD inhibitor is a tool for in vitro asIgE
why the CCD inhibitor-anti-CCD IgE complexes
testing, which helps reduce the effect of anti-CCD
can be washed out along with the other unbound
antibodies on the test results, and thereby
elements during the procedure, and thus do not
increases the specificity of the method. The appli-
influence the test result. The CCD inhibitor used
cation of this reagent can contribute to the better
in this study has been tested before in a group of
and faster diagnostics in cases of multiple
more than 1,300 CCD-positive patients, effec-
allergies tested in vitro. The present findings
tively reducing the number of false positives
lead to a conclusion that the use of CCD inhibi-
(Holzweber et al. 2013).
tion is a viable option in a clinical setting.
The present study shows that the use of a CCD
Multiparametric tests in the diagnostics of
inhibitor lowers the content of anti-CCD
IgE-dependent allergies facilitate the identifica-
antibodies directly analyzed in the test, which
tion of an influence of anti-CCD antibodies on
complies with the previous observations
the results for other allergens. However, due to
(Altmann 2016). The absence of an inhibitory
additional costs stemming from the need to retest
effect in some samples we analyzed remains
with CCD inhibitor, the use of a CCD inhibitor
unexplained. That could stem from the lack of
should be considered on a case-by-case basis.
duplicate measurements. Also, a small quantity
of the material tested or unique patient’s antibody
Conflicts of Interest The authors declare no conflicts of
characteristics may constitute some limitation. A interest in relation to this article.
low content of asIgE found in the present sample
of patients (median 1.5 kU/L, with no inhibitor)
can be yet another factor making the statistical
References
evaluation difficult, based on the CCD marker
alone. Be it as it may, other authors have also Aalberse RC, Koshte V, Clemens JG (1981) Immunoglob-
noticed the lack of a uniformly inhibitory effect ulin E antibodies that crossreact with vegetable foods,
on the anti-CCD antibody content of a CCD pollen, and Hymenoptera venom. J Allergy Clin
Immunol 68:356–364
inhibitor, although there was a marked reduction
Altmann F (2016) Coping with cross-reactive carbohy-
in false positives (Holzweber et al. 2013). drate determinants in allergy diagnosis. Allergo J Int
As exemplified by a sample showing, in 25:98–105
in vitro testing, different plant-based allergies, Ebo DG, Hagendorens MM, Bridts CH, De Clerck LS,
Stevens WJ (2004) Sensitization to cross-reactive car-
the action of the CCD inhibitor was clearly visible
bohydrate determinants and the ubiquitous protein
in case of reduced specific asIgE (Table 2). In profilin: mimickers of allergy. Clin Exp Allergy
case of the tomato allergen, asIgE dropped below 34:137–144
Inhibition of Cross-Reactive Carbohydrate Determinants in Allergy Diagnostics 79
81
82 I. Stanisławska et al.
regulation of cellular redox state (Kalinina et al. After weaning at 4–5 weeks, animals obtained
2015; Hultberg and Hultberg 2006). GSH also feed consisting of 16% protein for another 7 days.
plays a protective role against oxidative stress The composition and nutritional value of the feed
and damage in the hepatocytes, scavenging a are presented in Table 1. Then, the animals were
variety of free radicals (Wang et al. 2015). GSH divided into the glutathione- and 0.9% NaCl-
is notably present in the hepatocyte cytosol treated (control) groups, consisting of 15 animals
(approx. 90%), where it is synthesized in a each. Over the next 7 days, the mice received
two-step reaction, and in the mitochondrial matrix reduced glutathione (Sigma-Aldrich; St. Louis,
(10%) to which it passes from the cytosol MS) in a dose of 100 μg/g dissolved in 12 μL/g
(Koehler et al. 2006). body mass of 0.9% NaCl or saline alone. All
Lysosomes are the membrane-bound cellular injections were made i.p., twice daily, between
entities that contain hydrolytic enzymes whose 8.15–8.45 a.m. and 4.15–4.45 p.m. 6 h after the
main role, although not explicitly limited to, is last injection, under brief diethyl ether anesthesia,
to degrade and dispose of unwanted by the cell the mice were euthanized by cervical dislocation
material and its debris (Settembre et al. 2013). and were decapitated. Fragments of liver tissue
Oxidative stress contributes to the destabilization were immediately isolated, perfused with 0.9%
of lysosomal membranes causing a release of NaCl solution at 5 C, and placed in a 0.25 M
hydrolases into the cytosol, with a subsequent saccharose solution containing 2 mM EDTA, at a
induction of cellular apoptosis or necrosis (Bidere concentration of 1 g tissue per 6 ml solution. The
et al. 2003; Jóźwik et al. 2003). fragments were homogenized at 200 rpm in a
Recently, GSH supplementation is a heady Potter-Elvehjem tissue grinder with a Teflon
trend in the general public. The peptide is consid- pestle (Omni International; Kennesaw, GA),
ered a food additive with the presumed enhancing according to Beaufay’s (1972) method.
role in the cellular clearance and energy pro- Homogenates were subjected to sequentially
cesses, which also makes it of use in the sports increasing centrifugations: 8 min at 5,500 rpm in
field to increase physical efficacy. However, the a K26D centrifuge, 20 min at 14,000 rpm (Sorvall
influence of GSH on the activity of specific lyso- RC-5C), and 60 min at 10,000 rpm (Beckman
somal processes is largely unknown. The present L7–66; GMI, Ramsey, MN), according to the
study seeks to examine the effect of exogenous method of Marzella and Glaumann (1980).
GSH, applied in vivo, on the in vitro activity of The lysosomal and microsomal sediments
hydrolytic enzymes in the lysosomal, micro- were suspended in 100 mM phosphate-buffered
somal, and cytosolic fractions of the mouse liver. saline (PBS), pH 6.0, with the addition of 0.1%
Triton X100. They were then frozen and thawed was determined using a modified method of
twice and centrifuged for 10 min at 700 rpm in a Kirschke et al. (1972). All substrates were pur-
K-26D centrifuge. After another double freezing chased from SERVA Feinbiochemica GmbH &
and thawing, the cytosolic fractions were then Co (Heidelberg, Germany).
assayed. Data were presented as means SD.
Differences between group means of the glutathi-
one and control groups were compared using
analysis of variance ANOVA. A p-value <0.05
2.2 Cellular Enzymes
defined statistically significant differences. A
commercial statistical package of SAS/STAT
Activity of cellular enzymes was determined in
software was used (SAS Institute Inc., Cary, NC).
each of the subcellular fraction obtained during
the preparation procedures described in the sub-
section above, i.e., lysosomal, microsomal, and
cytosolic fractions, using the appropriate 3 Results
substrates and methods depicted in Table 2. A
Nicolet Evolution 300 BB spectrophotometer In the glutathione group, activity of β-GlcUr sig-
setup (Thermo Electron Corporation, Rugby, nificantly decreased, on average, in the lysosomal
Warwickshire, UK) was used for the enzyme fraction of hepatocytes (to 64.4% of that in the
activity determination. control group) and increased in microsomal
Activities of β-GlcUr, β-Gal, β-Glu, Hex, AcP, (to 174.9% of control) and cytosolic (to 119.9%
EL, and LL were determined using 4-nitrophenyl of control) fractions. The activity of β-Gal and
substrates at 420 nm wavelength according to β-Glu increased in all the fractions studied, while
Barrett and Heath’s (1977) micromethod. Cath the activity of Hex markedly increased in lyso-
D and Cath L were determined using 6% somal and cytosolic fractions only (to 130.5%
azocasein as a substrate at 366 nm wavelength and 136.2% of control, respectively) (Table 3).
according to Langner et al.’s (1973) method. GSH-fortified animal feed distinctly enhanced the
LeuAP and AlaAP were determined using Fast activity of all remaining hydrolytic enzymes in
Blue BB salt (4-benzoylamino-2,5-diethoxyben- the subcellular fractions of mouse liver, with the
zenediazonium chloride) substrates at 540 nm exception of esterase, leucine and alanine
wavelength according to McDonald and Barrett’s aminopeptidases, and cathepsins in the cytosol
(1986) method. The enzyme activity was and lipase in the microsomal fraction.
expressed in nmol/mg protein/h. The protein
Table 3 Activity of glycosidases in lysosomal (L), microsomal (M), and cytosolic (C) fractions of liver tissue in the
in vivo glutathione (GSH)-treated and control untreated mice
Enzyme Fraction Control GSH % of Control
L 303.1 30.0 195.1 17.0** 64.4
β-GlcUr M 94.0 9.2 164.4 16.0** 174.9
C 62.4 5.0 74.8 6.1* 119.9
L 347.8 35.0 450.5 40.0** 129.5
β-Gal M 59.5 5.8 112.9 11.0** 189.7
C 325.5 27.0 394.6 32.0** 121.2
L 117.1 12.0 215.3 19.0** 183.9
β-Glu M 23.9 2.2 34.2 3.3** 143.1
C 134.0 11.0 157.2 12.0* 117.3
L 842.1 84.0 1098.7 99.0** 130.5
Hex M 287.1 29.0 245.5 24.0 85.5
C 225.0 19.0 306.5 25.0** 136.2
L 1979.4 198.0 2571.4 233.0*** 129.9
AcP M 795.7 79.0 1226.9 122.0*** 154.2
C 1520.0 127.0 2444.5 203.0*** 160.8
L 1434.1 143.0 1932.0 175.0*** 134.7
E M 864.9 86.0 1017.7 101.0* 117.78
C 320.5 27.0 302.2 25.0 94.3
L 574.8 58.0 728.8 66.0*** 126.8
L M 592.6 59.0 576.4 57.0 97.3
C 134.2 11.0 157.6 12.0 117.4
L 95.6 9.7 143.7 12.0*** 150.3
AlaAP M 208.2 21.0 328.7 32.0*** 157.9
C 134.1 11.0 132.3 10.0 98.7
L 146.8 15.0 236.3 21.0*** 161.0
LeuAP M 232.5 23.0 359.8 35.0*** 154.8
C 75.1 6.1 80.7 6.6 107.5
L 22.7 2.4 33.7 2.9*** 148.5
Cath. D & L M 3.0 0.1 3.9 0.2*** 130.0
C 0.89 0.2 0.55 0.1*** 61.8
β-GlcUr β-glucuronidase, β-Gal β-galactosidase, β-Glu β-glucosidase, Hex N-acetyl-hexosaminidase, AcP acid phos-
phatase, E esterase, L lipase, AlaAP alanine aminopeptidase, LeuAP leucine aminopeptidase, Cath A & L cathepsins D
and L. Enzyme activity was expressed in nmol/mg protein/h; *p 0.05; **p 0.01, and ***p 0.001 for the intergroup
difference
Table 4 Synopsis of in vivo GSH supplementation on the in vitro activity of all hydrolases investigated in the cellular
fractions of mouse hepatocytes
Activity of hydrolases
Hepatocyte fraction Increase (n) No change (n) Decrease (n)
Lysosomal 9 0 1
Microsomal 8 2 0
Cytosolic 5 4 1
Total; n (%) 22 (73.4) 6 (20.0) 2 (6.6)
might thus play a protective role in liver dysfunc- Hum S, Koski KG, Hoffer LJ (1992) Varied protein intake
tional states by developing adaptive responses alters glutathione metabolism in rats. J Nutr 122
(10):2010–2018
aimed at rebalancing of a distorted “redox state” Jóźwik A, Śliwa-Jóźwik A, Fronczyk W, Kołątaj A (2003)
through the enhancement of cellular antioxidant The level of in brain young and old mice. J Physiol
capacity. The findings of this study imply a poten- Pharmacol 2:64–65
tial benefit of GSH supplementation for decreas- Kalinina EV, Chernov NN, Novichkova MD (2015) Role
of glutathione, glutathione transferase, and
ing oxidative damage, particularly in tissues of glutaredoxin in regulation of redox-dependent pro-
high metabolic demand, at the cellular level. cesses. Biochem Mosc 79(13):1562–1583
There is an urgent need to explore the safety of Kirschke H, Langner J, Wiederanders B, Ansorge S,
GSH-fortified diet, which requires alternative Bohley P (1972) Intracellular protein degradation.
IV. Isolation and characterization of peptidases from
study designs. rat liver lysosomes. Acta Biol Med Ger 28(2):305–322
Article in German
Conflicts of Interest The authors declare no conflict of Koehler CM, Beverly KN, Leverich EP (2006) Redox
interest in relation to this article. pathways of the mitochondrion. Antioxid Redox Sig-
nal 8(5–6):813–822
Langner J, Wakil A, Zimmermann M, Ansorge S,
Ethical Approval All procedures performed in this study Bohley P, Kirschke H, Wiederanders B (1973)
were in accordance with the ethical standards of the insti- Aktivitätsbestiung proteolytischer Enzym mit
tution at which the studies were conducted and with the Azokasein als substrat. Acta Biol Med Ger 31:1–18
1964 Helsinki Declaration and its later amendments. All (Article in German)
applicable international, national, and/or institutional Marzella L, Glaumann H (1980) Increased degradation in
guidelines for the care and use of animals were followed. rat liver induced by vinblastine. I. Biochemical charac-
The experiments were carried out with the approval of the terization. Lab Investig 42:8–17
Bioethics Committee of the Świętokrzyska Chamber of McCarty MF, DiNicolantonio JJ (2015) An increased need
Physicians in the city of Kielce in Poland (permit for dietary cysteine in support of glutathione synthesis
no. 49/2016). may underlie the increased risk for mortality associated
with low protein intake in the elderly. Age (Dordr) 37
(5):96
McDonald JK, Barrett AJ (1986) Exopeptidases. In:
References McDonald JK (ed) Mammalian proteases: a glossary
and bibliography. Academic, London, pp 111–144
Anthonsen HW, Baptista A, Drabløs F, Martel P, Petersen Morini P, Casalino E, Sblano C, Landriscina C (1991) The
SB, Sebastião M, Vaz L (1995) Lipases and esterases: response of rat liver lipid peroxidation, antioxidant
a review of their sequences, structure and evolution. enzyme activities and glutathione concentration of the
Biotechnol Annu Rev 1:315–371 thyroid hormone. Int J Biochem 23:1025–1030
Barrett AJ, Heath MF (1977) Lysosomal enzymes. In: Mosharov E, Cranford MR, Banerjee R (2000) The quan-
Dingle JT (ed) Lysosomes. A laboratory handbook. titatively important relationship between homocysteine
North-Holland Publishing Company, metabolism and glutathione synthesis by the transsul-
Amsterdam/New York/Oxford, pp 19–145 furation pathway and its regulation by redox changes.
Beaufay H (1972) Methods for isolation of lysosomes. In: Biochemistry 39:13005–13011
Dingle JT (ed) Lysosomes. A laboratory handbook. Ollinger K, Brunk U (1995) Cellular injury induced by
North-Holland Publishing Company, Amsterdam, pp oxidative stress is mediated through lysosomal dam-
1–30 age. Free Radic Biol Med 19:565–574
Bidere N, Lorenzo HK, Carmona S, Laforge M, Harper F, Radtke KK, Coles LD, Mishra U, Orchard PJ, Holmay M,
Dumont C, Senik A (2003) Cathepsin D triggers Bax Cloyd JC (2012) Interaction of N-acetylcysteine and
activation, resulting in selective apoptosis-inducing cysteine in human plasma. J Pharm Sci 101:4653–4659
factor (AIF) relocation in T lymphocytes entering the Seo AY, Xu J, Servais S, Hofer T, Marzetti E,
early commitment phase to apoptosis. J Biol Chem Wohlgemuth SE, Knutson MD, Chung HY,
278:31401–31411 Leeuwenburgh C (2008) Mitochondrial iron accumu-
Fernandez-Checa JC, Kaplowitz N (2005) Hepatic mito- lation with age and functional consequences. Aging
chondrial glutathione: transport and role in disease and Cell 7(5):06–716
toxicity. Toxicol Appl Pharmacol 204:263–273 Settembre C, Fraldi A, Medina DL, Ballabio A (2013)
Hultberg M, Hultberg B (2006) The effect of different Signals from the lysosome: a control centre for cellular
antioxidants on glutathione turnover in human cell clearance and energy metabolism. Nat Rev Mol Cell
lines and their interaction with hydrogen peroxide. Biol 14(5):283–296
Chem Biol Interact 163:192–198
Effects of Glutathione on Hydrolytic Enzyme Activity in the Mouse Hepatocytes 87
Śliwa-Jóźwik A, Jóźwik A, Fronczyk W, Guszkiewicz A, level and severity of portal hypertension in cirrhotic
Kołątaj A (2004) Effect of reduced glutathione (GSH) patients. Int J Clin Exp Med 8(7):11007–11013
on activity of lysosomal system in subcellular fractions Zhu J, Wu Y, Tang Q, Leng Y, Cai W (2014) The effects
of mouse kidney. Anim Sci Paper Rep 22(2):237–245 of choline on hepatic lipid metabolism, mitochondrial
Valencia E, Marin A, Hardy G (2001) Glutathione – function and antioxidative status in human hepatic
nutritional and pharmacologic viewpoints: part C3A cells exposed to excessive energy substrates.
IV. Nutrition 17:783–784 Nutrients 6:2552–2571
Wang SL, Zhu XY, Zhang DW, Zhang ZJ, Gao HJ, Yang
CQ (2015) Relevance of plasma malondialdehyde
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 89–109
https://doi.org/10.1007/5584_2018_261
# Springer Nature Switzerland AG 2018
Published online: 5 October 2018
89
90 M. A. Rieger et al.
iron (III) chloride and lime) combined with low 2.2.1 Inclusion Criteria
concentrations of other biological hazards (e.g., and Sociodemographic Data
total bacteria: 105–106 CFU/m3 and total molds: Students The following inclusion criteria were
<103 CFU/m3), and no relevant exposures to applied for students: healthy male non-smokers,
ammonia (Rieger and Diefenbach 2007; without any clinically manifested allergy and not
Steinberg 2001). Based on those results and fur- living or working in a nearby farm, sewage treat-
ther exploratory measurements, we selected four ment plant, or waste industry. The students were
sewage treatment plants of a regional association aged 24.7 (range 22–29) years and had a body
for sewage treatment in North Rhine-Westphalia mass index (BMI) of 23.7 (range 20.7–26.6) kg/m2
(Germany) for the present study. Data were (Table 1). After inclusion in the study, a screening
assessed in 2002–2003. test was performed at least 1 month prior to expo-
To get insight into the adaptation to endotoxin sure. The test consisted of a skin prick test with
exposure, we compared chronically exposed common aero-allergens. In addition, specific
workers on the sewage treatment plants to just immunoglobulin E (sIgE) in response to a variety
acutely exposed voluntary subjects (students) of ubiquitous aero-allergens was assessed in the
regarding the possible changes in various serum (atopy screen sx1 Phadiatop®, CAP-FEIA;
outcomes during one working shift. The workers ThermoFisher Scientific, Waltham, MA). Allergen
and students were matched in pairs, i.e., the sIgE 0.35 kU/L was considered positive. Based
workers were escorted by students during their on the screening results and medical history, 4 out
ordinary work on the plant. The students of the 10 students were classified as having atopy.
documented the workers’ activities in a log. To On the screening day, students also were trained by
investigate a possible dose-response relationship, a laboratory’s specialist to gather NALF.
students were asked to participate in the study
several times with intervals of several months On the day of examination, sociodemographic
between the repeated exposures. Four students data, general symptoms, and the judgment of the
were purposefully exposed to high concentrations stay on the sewage treatment plant were assessed
of endotoxin on two following days to assess a by means of a standardized questionnaire. None of
possible cumulative effect on inflammatory the students indicated any concurrent possible
markers. Lastly, specimens were obtained from exposure to biological hazards due to occupation
five students not only during their stay on sewage or leisure time activities. Concerning the occupa-
treatment plants but also after 3 days free of any tional history, only two students indicated former
specific exposure in offices and lecture rooms at occupational exposure to airborne biological
the university. Thus, variability of parameters hazards (agricultural work years before the study)
within 1 day could be compared with changes and none had any other previous inhalation expo-
evoked by specific exposure. sure. Concerning the judgment of stress due to the
working environment of the sewage treatment
plant, most of the inquired aspects were rated as
2.2 Study Sample little-to-moderate (heat, dust, chemicals, humidity,
cold, and time pressure) and four aspects as
As described above, a major goal of the study was moderate-to-high (noise, humidity, infectious
a comparison of chronically exposed workers in agents, and odor), which was assessed on a
sewage treatment plants (n ¼ 8) to just acutely 3-point Likert scale: 1 ¼ little stress and
exposed voluntary subjects, represented by 3 ¼ high stress. The self-reported students’ health
students (n ¼ 10). In addition, a possible dose- status was rated to be very good (mean 1.4), based
response relationship of the immunological on the following scale: 1 ¼ very good, 2 ¼ good,
reactions and a cumulative effect of exposure to and 6 ¼ insufficient. The evaluation of general
endotoxin on two following days were complaints and medical history showed no specific
investigated.
92 M. A. Rieger et al.
findings. Spirometry revealed no abnormalities in chemicals, humidity, cold, and time pressure and
the morning of the examination day. as moderate-to-high due to infectious agents and
odor. The majority of workers (7/8) indicated a
Workers The workers were recruited among the high work satisfaction. The evaluation of general
staff of selected sewage treatment plants working complaints and medical history showed no spe-
at the chamber filter press or sand trap. Only cific findings. The subjective ability to work was
workers fit for work on the day of examination rated slightly better (mean 2.3, range 1–4) than
were asked to participate in the study. All workers the own health status (mean 2.8, range 2–4)
were men, 3/8 workers were current and 2/8 were (1 ¼ very good, 2 ¼ good, and 6 ¼ insufficient).
former smokers. The mean age was 40.9 (range In the morning of the examination day, spirome-
24–56) years. The BMI was 28.4 (range try resulted in regular findings in all but one
23.6–37.6) kg/m2, with one worker being obese worker; the latter presenting with expiratory
according to the definitions of WHO. The mean flow limitation. Based on sIgE (in all but one
working duration on the sewage treatment plant worker) and the medical history, one worker
was 13.4 (range 4.0–33.8) years and the workers was described as being atopic.
indicated to work >6 h/day in the chamber filter As the recruitment of sewage treatment plants
press hall (Table 2). Beyond the work on sewage turned out difficult, the time interval between the
treatment plants, only two workers reported a measurements in single plants was between 1 and
previous occupational exposure to airborne 9 months. That made it possible that four students
biological hazards (agricultural work or sewage (s4, s7, s8, and s9) could be examined on two
sludge treatment of about 7 years’ duration) and different plants and one student (s5) on three
one worker reported a former inhalation exposure plants. According to the usual size of the sewage
from mining of 6 years’ duration. In a treatment plant, only two or three workers per
standardized questionnaire, the workers judged plant could be recruited for participation in the
the stress as little-to-moderate by such aspects of study.
the work environment as heat, noise, dust,
Table 2 Sociodemographic data of workers
Daily work in Work on Subjective evaluation of Work satisfaction Subjective ability to
SHSb (1 ¼ very chamber filter the sewage work environment on (single item; 1 ¼ little, work (school grades,
Age BMI Atopya good press hall plant sewage plantc (range 9–27 2 ¼ medium, 1 ¼ very good
Workers (years) (kg/2) (yes/no) 6 ¼ insufficient) (hours) (years) points) 3 ¼ high) 6 ¼ insufficient)
w1 38 28.1 No 4 >6 8.3 15 3 3
w2 43 27.8 No 4 >6 6.0 17 3 2
w3 38 26.2 No 2 >6 8.3 12 3 2
w4 24 23.6 No 2 >6 4.0 17 3 2
w5 49 29.1 No 3 >6 19.8 15 3 4
w6 44 26.3 Yes 3 >6 17.7 15 2 2
w7 56 37.6 No 2 <3 33.8 10 3 1
w8 35 28.7 No 2 Highly variable 9.5 16 3 2
All Mean Mean Yes: 2.8 > 6 h: 6/8 Mean 13.4 Mean 14.6 Mean 3 Mean 2.3
workers 40.9 28.4 1/8
a
Assessed by antibodies against inhalative allergens (IgE, sx1-test) and medical history
b
SHS Subjective health status
Adaptation to Occupational Exposure to Moderate Endotoxin Concentrations: A. . .
c
Working environment was judged by the workers regarding nine aspects: heat, noise, dust, chemicals, infectious agents, odor, humidity, cold, and time pressure (1 ¼ little stress,
2 ¼ moderate stress, 3 ¼ high stress)
93
94 M. A. Rieger et al.
Table 3 Endotoxin exposure on sewage treatment plants during the days of examination assigned to study participants
EU/m3
Plant no. (date) Students Workers Two measurements/daya; average
1 s1 w1 90/46; 68.0
(23/4/2002) s2 w2 17/16; 16.5
s3 w3 12/7; 9.5
2 s4 w4 32/41; 36.5
(22/5/2002) s5 w4 32/41; 36.5
s6 w5, (w6)b 99/91; 95.0
s7 – 121/47; 84.0
3c s7 w7 10/7; 8.5
(12/2/2003) s8 w7 10/7; 8.5
s9 w8 4/5; 4.5
s5 w8 4/5; 4.5
(13/2/2003) s10 w7 58/13; 35.5
s5 w8 12/33; 22.5
s9 – 51/8; 29.5
4c,d s4 318/234; 276.0
(3/11/2003) s5 798/885; 841.5
s8 – (No personal sampling); 465.0e
s9 394/162; 278
(4/11/2003) s4 536/643; 589.5
s5 231/1039; 635.0
s8 – (No personal sampling); 659.7e
s9 866/643; 754.5
a
Exposure was assessed by personal sampling during routine work on sewage treatment plants (especially in chamber
filter press hall and sand trap) before and after lunch break, sampling equipment usually carried by the worker(s), duration
of measurements: 112–265 min (mean 166 min, median 178 min)
b
W6 no full exposure due to part-time work on the day of examination
c
Plant 3 and plant 4: measurements on two following days
s student, w worker, EU endotoxin unit
d
Plant 4: students only, no routine work but high exposure near the sand trap
e
Student s8 no personal sampling but exposure assigned as mean value of exposure of students s4, s5, and s9 (each
personal sampling)
following soluble markers were measured: inter- PeliKine™-Tool sets (CLB, Amsterdam,
leukin (IL)-1β, IL-5, IL-6, IL-8, soluble (s)CD14, Netherlands) in a standard range of 0.4–300 pg/
tumor necrosis factor (TNF)-α, eosinophilic cat- mL for IL-1β and a standard range of
ionic protein (ECP), NO, total protein content, 1.4–1000 pg/mL for TNF-α. IL-8 was measured
and albumin. The present evaluation of cellular with the OptEIA™ ELISAs (BDBiosciences
and soluble inflammatory markers in NALF Pharmingen, Heidelberg, Germany) in a standard
focuses on total cell count, proportion of range of 3–200 pg/mL and most of the samples
neutrophils, and the concentrations of total pro- must be diluted to reach the standard range. IL-6
tein, IL-8, IL-6, TNF-α, IL-1β, and sCD14. and sCD14 were determined with the DuoSet™
Therefore the methods of these selected soluble ELISA (R&D Systems, Wiesbaden, Germany) in
markers are described in more detail. a standard range of 3–600 pg/mL for IL-6 and of
Interleukins, sCD14, and TNF-α were 250–16,000 pg/mL for sCD14. Total protein con-
measured using purchased monoclonal “sand- tent was determined according to the method of
wich” enzyme immunoassay (EIA) kits. IL-1β Bradford (1976) with bovine serum albumin as a
and TNF-α were measured with the respective standard solution (range 10–100 μg/mL). For all
96 M. A. Rieger et al.
soluble markers, concentrations below the limit of For this major goal of the study, 12 pairs of
quantification (LOQ) were set at 2/3 of LOQ. workers and students were combined and
investigated on three sewage treatment plants
2.4.3 Spirometry (see Table 4)
Vital capacity (VC), forced vital capacity (FVC),
and forced expiratory volume in 1 s (FEV1) as 2. Variability of inflammatory markers without
well as peak expiratory flow (PEF), peak inspira- specific exposure, possible dose-response rela-
tory flow (PIF), and maximal instantaneous tionship, and cumulative effects of the endo-
forced expiratory flow when x% of the FVC toxin exposure:
remains to be expired (MEFx%; 25%, 50%, and
75%) were measured pre- and post-shift, as well Five students (s4, s5, s7, s8, and s9) were
as in the next morning with a portable spirometer investigated without any specific exposure with
(Jaeger SpiroPro®; Viasys Würzburg, Germany) a time interval of at least 3 weeks after the last
using the software V1.0. Measurements were exposure to endotoxin (s4, s5, s8, and s9:
performed according to the guidelines of the 3 weeks, s7: 9 months). This examination was
American Thoracic Society (ATS 1995), but repeated twice each with a 1-week interval.
repeated measurements were not possible at all Dose-related changes of inflammatory markers
time points due to organizational constraints on concentrations were investigated based on
the sewage treatment plants. Reference values 20 values as some students participated repeat-
were chosen according to Quanjer et al. (1993). edly in the study. The students’ values were
grouped according to three categories of endo-
2.4.4 Questionnaires toxin concentration: <50 EU/m3 – eight students;
The occupational and medical history (including 50–200 EU/m3 – three students; >200 EU/m3 –
smoking), the current occupation, potential other four students. Possible cumulative effects could
sources of exposure (e.g., leisure time activities, be investigated in four students (s4, s5, s8, and s9)
or place of living), and the subjective stress who were exposed to high endotoxin
caused by the working conditions on the sewage concentrations on plant 4 on two following days
treatment plant were assessed by means of a (see Table 5).
standardized questionnaire and by the subjective
evaluation of health and workability. In addition,
the study participants were asked to indicate the 2.6 Statistical Analysis
extent of symptoms in line with a mucous mem-
brane irritation syndrome (e.g., symptoms from Data were expressed as median with min–max
the eyes, nose, and throat), general effects of range (tables) or stem and leaf plots (box plots)
exposure to endotoxin conforming to an organic illustrating the interquartile range (75th percentile
dust toxic syndrome (e.g., limb pain, fever, and – 25th percentile) and minimum and maximum
feeling like getting a flu), or other symptoms (e.g., (figures). As the concentration of inflammatory
nausea and diarrhea) at time point A, B, and C, markers in NALF showed no normal distribution,
using a Likert scale (0 ¼ not at all, 4 ¼ very logarithmic (total cells, total protein, IL-8, IL-1β,
strong). and sCD14) or dichotomous values (neutrophils
categorized below/equal or above 40%, IL-6 and
TNF-α categorized as below/equal or above the
2.5 Sample Size for Single Research respective limit of quantification) were used for
Issues statistical analysis comparing workers and
students (research issue 1). In this context, after
1. Comparison of chronically and only acute univariate analysis (logarithmic values: regres-
exposed subjects: sion analysis and analysis of variance, with a
UNIANOVA procedure; dichotomous values:
Table 4 Cellular and soluble inflammatory markers in nasal lavage fluid (NALF) before exposure (time point A) in chronically exposed (workers) and so far not exposed
subjects (students). These data were assessed on plant 1–3 (only Day 1)
Students (n ¼ 10; n ¼ 12 values) Workers (n ¼ 8; n ¼ 8 values) Students vs. workersa
Number of values >LOQ (%; Number of values >LOQ (%; Univariate
NALF 100% ¼ 12) Median (range) 100% ¼ 8) Median (range) analysis GEE
Total cell count 12 (100%) 6.0 (4.0–12.0) 8 (100%) 10.0 (4.0–24.0) p ¼ 0.090 p ¼ 0.079
(104)
Neutrophils (%) 12 (100%) 78.8 (27.5–93.3) 8 (100%) 35.8 (0.0–94.0) ns/nt ns/nt
Total protein 12 (100%) 58.5 (18.1–160.3) 8 (100%) 61.3 (10.0–370.3) ns/nt ns/nt
IL-8 (pg/mL) 12 (100%) 147.1 (38.3–802.3) 8 (100%) 231.9 (59.1–1162.1) p ¼ 0.024 p ¼ 0.128
IL-6 (pg/mL) 1 (8%) 2.0 (2.0–43.2) 2 (25%) 2.0 (2.0–18.4) ns/nt ns/nt
TNF-α (pg/mL) 8 (67%) 1.6 (0.9–90.9) 5 (63%) 1.6 (0.9–20.2) ns/nt ns/nt
IL-1β (pg/mL) 10 (83%) 1.5 (0.3–6.3) 8 (100%) 2.1 (0.8–36.6) p ¼ 0.019 p ¼ 0.085
sCD14 (pg/mL) 12 (100%) 3461.0 (1657.8–12083.8) 8 (100%) 3907.6 (1266.3–12876.2) ns/nt ns/nt
a
This statistical analysis is based on the values of all three time points; at time point A, a marked difference was found for the total cell count only
ns no statistically significant difference between students and workers (i.e., p > 0.05)
Adaptation to Occupational Exposure to Moderate Endotoxin Concentrations: A. . .
nt no tendency regarding significant differences between students and workers (i.e., p > 0.100)
LOQ limit of quantification, GEE general estimating equation
97
98
Table 5 Cellular and soluble inflammatory markers in nasal lavage fluid (NALF) in students at time point A (morning) and B (afternoon) according to endotoxin exposure during
the day
Time point A
(morning) Time point B (afternoon after endotoxin exposure)b
(n ¼ 10 students/ No specific exposurea < 50 EU/m3 (n ¼ 8 50–200 EU/m3 (n ¼ 3 > 200 EU/m3 (n ¼ 4/n ¼ 8 values); median (range)
n ¼ 31 values); (n ¼ 5 students/n ¼ 14 students/n ¼ 9 values)c; students/n ¼ 3 values)c;
NALF median (range) valuesd); median (range) median (range) median (range) Day 1c Day 2c
Total cellsa 8.0 (4.0–16.0) 4.0 (4.0–12.0) 8.0 (4.0–12.0) 8.0 (8.0–8.0) 4.0 (4.0–8.0) 4.0 (4.0–4.0)
104
Neutrophils 67.7 (0.0–93.3) 38.7 (0.0–84.7) 56.1 (9.0–94.0) 27.0 (16.0–81.3) 12.2 (1.0–55.0) 60.3 (7.0–66.0)
(%)
Total 73.6 (6.7–359.1) 74.7 (21.6–372.7) 37.5 (21.3–150.9) 42.9 (21.5–43.3) 49.6 (44.2–164.7) 47.9 (33.6–133.5)
protein (μg/
mL)
IL-8 246.7 373.6 (96.2–2982.0) 137.1 (38.8–572.2) 105.5 (89.0–414.2) 379.9 (169.2–636.9) 558.2 (51.7–1324.2)
(pg/mL) (38.3–1590.9)
IL-6 2.0 (2.0–43.2) 2.0 (2.0–33.6) 2.0 (2.0–26.2) 2.0 (2.0–2.0) 2.0 (2.0–2.0) 2.0 (2.0–5.7)
(pg mL)
TNF-α 0.9 (0.9–90.9) 1.7 (0.9–3.2) 1.5 (0.9–9.5) 6.6 (0.9–7.6) 0.9 (0.9–7.7) 0.9 (0.9–1.6)
(pg/mL)
IL-1β 2.8 (0.3–17.7) 1.2 (0.3–15.7) 0.9 (0.3–4.2) 0.6 (0.5–1.7) 1.8 (1.4–5.7) 1.4 (0.9–3.4)
(pg/mL)
sCD14 3548.8 3337.4 (1563.8–18374.0) 2968.0 (943.2–7415.5) 2495.9 (2046.5–3473.3) 3108.2 (1937.3–6620.7) 2426.8 (1943.2–7741.1)
(pg/mL) (275.3–13574.1
The number of values is based on repeated participation of students in the study (see Table 3)
a
No specific exposure: offices and lecture rooms of the university
b
Endotoxin exposure derived from the students’ stay on a sewage treatment plant (plant 1–3: escorting workers during their routine working day)
c
Concentrations per category: <50 EU/m3: median 17 EU/m3 (4–36 EU/m3), 50–200 EU/m3: median 84 EU/m3 (70–95 EU/m3); >200 EU/m3: Day 1 median 356 EU/m3
(162–885 EU/m3; mean: 465 EU/m3), Day 2 median 643 EU/m3 (231–1039 EU/m3; mean 659.7 EU/m3)
d
Student 8: did not participate on the third day without specific exposure
M. A. Rieger et al.
Adaptation to Occupational Exposure to Moderate Endotoxin Concentrations: A. . . 99
logistic regression) regarding the possible effects work shift. In addition, workers w3 and w4
either due to group (chronic or only acute expo- showed a tendency of exposure-induced changes
sure), time points (A, B, and C), or the in lung function, whose significance remains
interactions of group and time points, generalized unclear due to methodological reasons i.e. no
linear models (GLM) with generalized estimating possibility of repeated measurements due to orga-
equations (GEE) were calculated, as some nizational constraints on the sewage treatment
students participated several times in the study. plants. All students had normal spirometry results
In addition, the possible effect of group and plant at all time points and showed no signs of
on the values gathered at time point A was exposure-induced changes in lung function.
analyzed by univariate analysis (logarithmic
values) or logistic regression (dichotomous
values). Due to a limited number of subjects, we 3.2 Research Issue 1: Effects
additionally reported results ignoring the cluster- of Chronic Endotoxin Exposure
ing caused by repeated measurements of the same
subjects. A two-sided significance level of 0.05 The emphasis of the present study was on the
was chosen for all tests. According to small num- possible effects of the endotoxin exposure on
bers, no statistical analysis was possible for cellular and soluble inflammatory markers in
research issue 2. Data were analyzed using a NALF. For the comparison of chronically
commercial SPSS software v22 and v23 (IBM exposed workers with acutely exposed students,
Co, Armonk, NY). 12 pairs of workers and students based on data of
8 workers and 10 students could be examined on
3 sewage treatment plants. There were significant
3 Results differences between both groups in IL-8 (univari-
ate analysis: p ¼ 0.024, GEE: p > 0.128) and
3.1 General Findings IL-1β (univariate analysis: p ¼ 0.019, GEE:
p ¼ 0.085), each with a lower concentration in
None of the healthy volunteers cancelled the students than workers. Similarly, total cell count
exposure on the sewage treatment plant. Only tended to be lower in students than workers (uni-
after the exposure to a high endotoxin concentra- variate analysis: p ¼ 0.090, GEE: p ¼ 0.079).
tion on plant 4, some students indicated TNF-α was the only inflammatory marker with a
symptoms in line with a mucous membrane irrita- significant concentration difference at the three
tion in the standardized questionnaire (increase time points (univariate analysis: p ¼ 0.079,
from “no symptoms” to “little” or “clear” irrita- GEE: p ¼ 0.016) with fewer values below LOQ
tion especially of the nose) which were no longer at time point A, compared to time point C. For
prevalent on the next morning. One student TNF-α, no significant effect or interaction with
reported stomach ache, diarrhea, and nausea on the variable “group” could be detected. For the
the morning of the second day of exposure on other inflammatory markers, no effect or interac-
sewage plant 4, but the symptoms eased during tion regarding “group” or “time point” could be
the day. The sewage workers did not report any seen in the statistical analysis.
changes of symptoms during the working day in Concerning the possible effect of
the standardized questionnaire. chronic vs. acute exposure (variable “group”) or
the effect of plant on the concentration resulting
3.1.1 Lung Function at time point A, a tendency toward a significant
All but one worker (w8) had normal values in the difference was found for the variable “group”
spirometric examination before the exposure and only in the total cell count in NALF (p ¼ 0.064)
showed no clinically relevant changes after the (lower concentration in students) and the variable
work shift. Worker w8 had a peripheral flow “plant” turned out to be significantly associated
limitation which seemed to increase during a with a higher concentration of this marker on
100 M. A. Rieger et al.
plant 3 than on plant 1 (p ¼ 0.034). Likewise, endotoxin concentration than that without spe-
logistic regression showed a significant effect of cific exposure (median: 12.2% (exposure >200
the variable “plant” on TNF-α (p ¼ 0.027), with EU/m3) vs. 53.3% (no specific exposure)). Gen-
more results below the limit of detection on plant erally, concentrations of inflammatory markers
1 (five samples) than on plant 2 or 3 (only one without specific exposure included most other
sample each). Concerning both markers, there values gathered after endotoxin exposure in the
was no significant interaction between the samples investigated in the present study (see
variables “group” and “plant”. For the other Table 5).
inflammatory markers gathered at time point A, Concerning the pattern of inflammatory
no such effects of the two variables could be seen marker concentrations, a higher endotoxin con-
in this univariate analysis. The concentrations of centration seemed associated with an increased
inflammatory markers in students and workers at concentration of total protein, IL-8, IL-1β, and
time point A are given in Table 4. sCD14, and with a decrease in the percentage of
neutrophils and in the concentration of IL-6 in
descriptive analysis. Differently, comparing the
inflammatory marker concentrations in the after-
3.3 Research Issue 2: Different
noon of the first and second days with high endo-
Patterns of Inflammatory
toxin exposure on plant 4 (mean endotoxin
Markers in NALF After Exposure
concentration 465 EU/m3 (Day 1) and 659.7
to a Moderate Endotoxin
EU/m3 (Day 2) (for details regarding the exposure
Concentration as Compared
see Table 3 and legend to Table 5)), marked
to Nonspecific Exposure. Dose-
differences between the values gathered at time
Response Relationship of Acute
point B on the first and second days seemed
Exposure, Including
detectable regarding the total cells (lower on
a Cumulative Exposure
Day 2), neutrophils (higher on Day 2), IL-8
(higher), IL-6 (higher), and sCD14 (lower) (see
The inflammatory markers in NALF at time point
last columns in Table 5 and Fig. 1). Compared to
B for all students participating in the study were
the initial findings on the morning of Day 1 (time
analyzed according to the respective endotoxin
point A1), the concentrations found on the morn-
exposure during the day (plant 3: values of the
ing of Day 2 (C1/A2) and Day 3 (C2) during the
first day of exposure only) and compared to time
cross-shift examination showed almost complete
point B values without specific exposure and to
recovery regarding the neutrophils (C1 and C2)
all values gathered before the exposure (i.e., time
but slight changes regarding IL-6 and especially
point A) (Table 5). The values gathered from five
IL-1β (C1 and C2), and an increase in variability
students at time point B, on 3 days without any
especially at time point C1 regarding the total
specific exposure (see Table 5 – column “no
cells, total protein, TNF-α, and sCD14 (Fig. 1).
specific exposure”), and the respective time
Due to a small sample size, combined with a
point A- and C-values (data not shown) showed
great variability in the concentration of all
a great variability and there was no marked dif-
markers investigated, no statistical analysis was
ference between A-, B-, and C-values. Due to this
possible to compare the values at time point B,
great variability, comparison of NALF markers
according to different kinds of exposure or to
gathered on sewage plant 4 with those gathered in
compare between all of the time points (A1 to
the same four students during 3 days without
C2), or between the three time points of each day
specific exposure revealed no significant differ-
(A1 to C1 and A2 to C2), with or without specific
ence between both conditions (see Table 5). Only
exposure. The study’s main finding was the great
the percentage of neutrophils at time point B1
variability of inflammatory marker concentrations
tended to be lower after exposure to a high
in naïve subjects without specific exposure.
12 100
NALF total cells (10e4)
4 0
A1 B1 C1/A2 B2 C2 A1 B1 C1/A2 B2 C2
15 14
10
10
8
6
4
5
2
0
A1 B1 C1/A2 B2 C2 A1 B1 C1/A2 B2 C2
Time of collection (Day 1 & Day 2) Time of collection (Day 1 & Day 2)
350
NALF total protein (mg/mL)
1000
300
NALF IL-8 (pg/mL)
250
200
100
150
100
50
0 10
A1 B1 C1/A2 B2 C2 A1 B1 C1/A2 B2 C2
6 12000
10000
NALF sCD14 (pg/mL)
NALF IL-1b (pg/mL)
5 8000
6000
4
4000
3
2
2000
1
0
A1 B1 C1/A2 B2 C2 A1 B1 C1/A2 B2 C2
Time of collection (Day 1 & Day 2) Time of collection (Day 1 & Day 2)
Fig. 1 Cellular and soluble inflammatory markers in nasal lavage fluid (NALF) in four students at time point A (before
exposure), time point B (afternoon of Day 1 and Day 2), and time point C (the morning after exposure during Day 1 and
Day 2); mean endotoxin exposure of 465 EU/m3 (Day 1) and 659.7 EU/m3 (Day 2)
102 M. A. Rieger et al.
Concerning condition “no specific exposure”, the latter question has not yet been addressed in
we decided not to perform any exposure assess- the occupational settings, where endotoxin is the
ment, as all office and lecture rooms, where the leading biological hazard, as opposed to studies
students stayed during the 3 days of the investi- examining workers and naïve controls in occupa-
gation, were in a good condition and had a good tional settings with more mixed biological
indoor climate with no apparent risk of elevated exposures, e.g., waste collection (Wouters et al.
concentrations of endotoxin, as reported to occur 2002) or pig farming (Sahlander et al. 2012; von
due to dampness in schools (Jacobs et al. 2014; Essen and Romberger 2003). As the workers on
Jacobs et al. 2013). the sewage treatment plants were only men, we
recruited no female students. Thus, we could not
Inflammatory Markers in NALF The noninva- examine any gender effect despite this might have
sive collection of NALF turned out well applica- been interesting according to the findings of a
ble for this field study, and it has been applied in Taiwanese study addressing the recovery in lung
many studies addressing the reaction to endotoxin function after cessation of exposure to organic
exposure in both field and laboratory settings dust (Lai et al. 2015).
(Raulf et al. 2017; Sigsgaard et al. 2015; Wouters There were some differences in the
et al. 2002; Douwes et al. 2000). According to the sociodemographic data of workers and students.
concept of “united airway disease” in the occupa- The workers were older, some were smokers, and
tional setting, inflammatory markers assessed in some showed signs of peripheral flow limitation
NALF may not only point to the inflammatory before (one worker) or after (three workers) the
processes in the upper but also lower airways work shift. This is in line with cross-sectional
(Quirce et al. 2010). As the students were trained studies in wastewater industries (Steinberg
in collecting their NALF at the study onset, the 2001), waste collectors (Raulf et al. 2017;
recovery rate and quality of specimens were gen- Schantora et al. 2015), compost workers
erally high. The time points for NALF collection (Hoffmeyer et al. 2014), or agricultural workers
(B about 8 h and C about 24 h after time point A) (Spierenburg et al. 2015) presenting with respira-
were chosen according to other studies, where tory symptoms during their employment. How-
NALF was gathered at 1 h (Danuser et al. ever, with respect to the ATS recommendations,
2000), 3 h (Sigsgaard et al. 2000), 6/7 h methodological limitations of our spirometry
(Doreswamy et al. 2011; Wang et al. 1997), or assessment have to be quoted. In general, there
23/24 h or later after endotoxin exposure was no possibility of repeated measurements due
(Danuser et al. 2000). In most of those studies, to organizational constraints on the sewage treat-
the exposure-related effect in inflammatory ment plants. On the other hand, students were
markers has occurred 3 or 6 h after exposure, only non-smokers and had normal spirometry
with a recovery approximately 24 h after a single findings.
exposure. Thus, time points B and C seemed to be Atopy was determined in 1/8 workers and 4/10
chosen properly in the present study. students by objective parameters. Students rated
their general health better than workers
Study Sample The constitution of pairs of (1.4 vs. 2.8; 1 ¼ very good, 6 ¼ insufficient)
workers (i.e., chronically exposed subjects) and and indicated a higher stress than workers due to
acutely exposed healthy volunteers (students), some aspects of the working environment on the
escorting the workers, was designed to compare: sewage treatment plants, e.g., noise 2.7
(a) the baseline concentration of inflammatory (students) vs. 1.8 (workers) or humidity 2.4
markers in NALF in the morning and (b) the (students) vs. 1.4 (workers) (1 ¼ little, 3 ¼ high
patterns of inflammatory markers in NALF after stress), whereas both groups indicated a similarly
exposure to endotoxin in the afternoon and on the high stress due to odor and infectious hazards.
next morning in both groups. To our knowledge,
104 M. A. Rieger et al.
The fact that the workers, despite partly was not reflected in statistical analysis (n ¼ 12
longstanding work on sewage treatment plants, values from students, n ¼ 8 values from workers).
indicated a rather good health status and even The difference between workers and only
better ability to work in conjunction with a high acutely exposed students is similar to the findings
work satisfaction might point to a healthy worker comparing current and former compost workers
survival effect. Yet, one out of eight workers to white collar workers. There, current workers
showed pathologic findings in spirometry before show a higher concentration of inflammatory
exposure, and he and two additional workers markers in induced sputum than former compost
showed a tendency of exposure-induced changes workers or nonexposed subjects (Raulf et al.
in lung function at time point B. According to the 2015). Similarly, before exposure, workers from
questionnaire, some students experienced slight biofuel plants have a higher IL-8 concentration in
symptoms of mucous membrane irritation, but no NALF compared to controls from conventional
severe complaints. To contrast the findings fuel plants (Zheng et al. 2014). Pig farmers com-
derived from the high exposure, some students pared to naïve controls show signs of local and
were also asked to participate in the assessment systemic inflammation associated with altered
during 3 days without specific exposure. Yet, no innate immunity, which has been ascribed to
data of chronically exposed subjects were avail- adaptation to chronic endotoxin exposure
able in either condition. (Sahlander et al. 2012; von Essen and Romberger
2003). So far, this effect has not been described
for workers in wastewater industries.
Other than in studies comparing farm workers
4.2 Research Issue 1: Effects
to naïve subjects in swine confinement (Sahlander
of Chronic Endotoxin Exposure
et al. 2012; von Essen and Romberger 2003), in
the present study there was no definable higher
We aimed to answer the question of whether there
concentration of inflammatory markers in
was a difference between chronically exposed
students than in workers at time point B or C, as
(workers of sewage treatment plants) and acute
a result of exposure to endotoxin on plant 1, 2, or
exposed subjects (students) regarding the markers
3. This might be due to a rather low endotoxin
in NALF before and after exposure to endotoxin
concentration associated with routine work on the
during a routine work on a plant. The results of
days of examination and the fact that the sample
univariate analysis and GEE pointed to a certain
size was too limited to detect differences between
difference between students and workers, which
workers and students.
can be interpreted as the effect of chronic expo-
sure: workers showed higher concentrations of
some inflammatory markers in NALF (i.e., total
4.3 Research Issue 2: Different
cell count, IL-8, and IL-1β) across all three time
Patterns of Inflammatory
points but without significant clinical symptoms
Markers in NALF After Exposure
or findings due to the exposure. Focusing only on
to a Moderate Endotoxin
the values gathered at time point A, this differ-
Concentration as Compared
ence occurred in the total cell count. To summa-
to Nonspecific Exposure. Dose-
rize, the effect of chronic exposure may result in
Response Relationship of Acute
higher concentrations of inflammatory markers in
Exposure, Including
NALF. Yet, these findings are based on only a
a Cumulative Exposure
small sample size (n ¼ 10 students, n ¼ 8
workers). In addition, regarding the comparison
We found a great intra- and interindividual
of concentrations at time point A, the uneven
variability concerning the concentrations of
distribution due to repeated participation of
inflammatory markers in NALF gathered in five
some students has to be accounted for, as this
healthy students in the morning and afternoon of
Adaptation to Occupational Exposure to Moderate Endotoxin Concentrations: A. . . 105
3 days without any specific exposure. Due to a comparing time point B values of the students
small sample size we could not perform any sta- exposed to different endotoxin concentrations,
tistical analysis, but descriptive analysis suggests higher endotoxin concentrations seemed
that there were no circadian differences in the associated with increased concentration of total
inflammatory markers investigated. protein, IL-8, IL-1β, and sCD14, and with
In the literature, we could not find any study decreased percentage of neutrophils and IL-6
focusing on the variability of inflammatory concentration. A cumulative exposure to the
markers in NALF (or other specimens) in naïve endotoxin concentration > 200 EU/m3 (Day 1:
subjects without a specific exposure. Yet, hints median 356 EU/m3 (162–885 EU/m3; mean:
for a great variability can be driven from studies 465 EU/m3), Day 2: median 643 EU/m3
where the results of (sometimes repeated) (231–1039 EU/m3; mean 660 EU/m3)), for
measurements before exposure are reported. A about 8 h each, seemed associated with a differ-
broad concentration range of inflammatory ence in the concentration of time point B values
markers in NALF has been reported, e.g., in for- collected on the first and second day. The differ-
mer waste collectors, especially in never smokers ence, particularly, were present in total cell count
(Raulf et al. 2017), in former compost or white (lower on Day 2), percentage of neutrophils
collar workers in induced sputum (Raulf et al. (higher on Day 2), and concentrations of IL-8
2015), or in NALF and exhaled breath condensate (higher), IL-6 (higher), and sCD14 (lower). Yet
in workers from conventional fuel plants (Zheng these effects should be considered cautiously
et al. 2014). A repeated collection of NALF with- against the background of a great intra- and inter-
out exposure (0 min, 30 min, and 2, 6, and 24 h) individual variability of inflammatory markers in
in six healthy volunteers or pre-exposure in NALF specimens collected at time point B in five
16 healthy volunteers (Raulf-Heimsoth et al. students on 3 days without any specific exposure,
2010) has shown a broad range of inflammatory as well as compared to a broad range of values
markers comparable to the present findings. collected at time point A in the students.
Unfortunately, other authors have not reported Beside a small sample size combined with
the range of repeated measurements of inflamma- great variability of concentrations, a reason for
tory markers in healthy subjects, but just the the lack of a distinguished dose-dependent reac-
means or medians, e.g., in 20 healthy volunteers tion in NALF inflammatory markers in the pres-
with repeated measurements (Kitz et al. 2008) or ent study could be that the endotoxin exposure
in 9 healthy volunteers with rather “similar” during routine work (plant 1–3) or at specific
concentrations (Michel et al. 1997). Therefore, places known to result in a high endotoxin expo-
no comparison is possible between the present sure (plant 4) was too small to induce a detectable
results and those reported previously. Interest- response. So far, there have been several studies
ingly, the median concentration of inflammatory where the dose-dependent changes of inflamma-
markers at different time points after saline instil- tory markers are investigated, mainly with the
lation reported by Danuser et al. (2000) do not goal to describe the effect of the lowest possible
show such high variability as found in the present level of endotoxin exposure in healthy volunteers
study. Thus, our findings seem rather new as a or chronically exposed workers (Danuser et al.
great variability of inflammatory markers has so 2000; Michel et al. 1997).
far been described explicitly only in response to Soluble inflammatory markers in NALF are
organic dusts (Sigsgaard et al. 2005). the main outcomes in an experimental study of
In the present study, the concentration of Danuser et al. (2000), where seven subjects were
inflammatory markers without a specific endo- exposed to 0, 10, and 40 μg lipopolysaccharide
toxin exposure in the students entailed the major- (LPS) by nasal instillation for 10 s. The instilla-
ity of values collected after the exposure (see tion of 10 μg LPS increased IL-6 over 2-fold and
Table 5). When possible dose-related changes in that of 40 μg caused an increase in nasal
the marker concentration were investigated by symptoms and a 13-fold increase in IL-6
106 M. A. Rieger et al.
compared to saline treatment, with a peak at 6 h. are akin to the present findings concerning
As there were no repeated measurements in single neutrophils and sCD14, at variance with IL-8.
individuals, the variability of the effects is To wrap it up, there still seems to be a lack of
unknown. Yet, comparing the present findings replicated inhalation studies in healthy volunteers
to those of that study suggests that the exposure describing a consistent dose-response relation-
on plant 4 was high enough to produce a distinct ship or cumulative effects due to endotoxin expo-
reaction in the NALF inflammatory markers. sure. That holds true for both work-related (Thorn
Michel et al. (1997) published the results of a 2001) and environmental exposures (Walser et al.
study on the dose-response relationship to inhaled 2015). Concerning the present study, it remains
endotoxin, focusing on clinical symptoms, blood unclear whether exposure to a far higher endo-
and induced sputum indices, and lung function in toxin concentration than found on plant 4 would
nine healthy volunteers after a weekly bronchial have induced changes in inflammatory markers
challenge with increasing doses of 0.5, 5, and distinct from the greatly variable concentrations
50 μg LPS. They noticed significant increases in found without any specific exposure in the
the number of neutrophils and monocytes in the students. Inversely, it could be that the apparent
sputum after inhalation of 5 μg LPS, which was changes described in the experimental studies
significantly greater after 50 μg LPS, compared to cited would not come to sight if those studies
saline. No repeated measurements were reported. accounted for the variability of findings without
Thus, variability of the effects in the nine specific exposure.
volunteers is unknown. The authors have set a
threshold value for no response in sputum inflam-
matory cells at 0.5 μg inhaled LPS, which 5 Conclusions
corresponds to the exposition for 10 h a day and
lung ventilation of 1 m3 per hour of air containing We replicated previous findings of a great
50 ng LPS. This concentration corresponds to variability of the concentration of inflammatory
about 500 EU/m3, which was exceeded by the markers in NALF in healthy subjects without
mean concentration the majority of students specific exposure at different time points, i.e.,
were exposed to for about 8 h on plant 4 in the within 1 day and between several days. The con-
present study. Yet, in the cited study, significant centration range was even broader than the appar-
effects in sputum inflammatory cells were ent reaction after exposure to a rather high
observed only at a concentration of 5 μg LPS endotoxin concentration (mean: 465 EU/m3
and higher, corresponding to a rather high expo- (Day 1)/660 EU/m3 (Day 2)). Therefore, in our
sure of about 5000 EU/m3. Thus, exposure even study we could not replicate findings on a dose-
on sewage treatment plant 4 possibly was too low response relationship or cumulative effects
to evoke distinct effects. obtained by other authors.
Concerning the effect of repeated exposure Based on our findings, it seems necessary to
(Day 1 and Day 2 on plant 4) at time point B on carefully assess individual ranges of concentra-
the total cell count (lower on Day 2), neutrophils tion for single inflammatory markers without any
(higher on Day 2), IL-8 (higher), IL-6 (higher), specific exposure before investigating a reaction
and sCD14 (lower on Day 2), our present results to a specific exposure such as endotoxin or other
can be compared with a study of Doreswamy hazardous substances. In addition, a great
et al. (2011). Those authors have reported a sig- variability of findings in single individuals
nificant increase in neutrophils and a decrease in underscores the importance of a sample size,
sCD14 and IL-8-concentration in NALF after a which should be large enough for the use of
6-hour experimental exposure (instillation of non-parametric tests (i.e., 20 subjects) or corre-
2 μg LPS in the nostril) on the third day com- spond to the sample size indicated in the test-
pared to pre-exposure or 6 h after just a single specific table accounting for the predictable
exposure in 15 healthy volunteers. Those results great variability of findings (Wilcoxon 1945).
Adaptation to Occupational Exposure to Moderate Endotoxin Concentrations: A. . . 107
Concerning the adaptation to chronic endo- Bradford MM (1976) A rapid and sensitive method for the
toxin exposure, our results suggest a higher con- quantitation of microgram quantities of protein
utilizing the principle of protein-dye binding. Anal
centration of inflammatory markers in chronically Biochem 72:248–254
versus acutely exposed subjects. Yet, as in our Castellan RM, Olenchock SA, Kinsley KB, Hankinson JL
study the exposure during routine work was not (1987) Inhaled endotoxin and decreased spirometric
high enough to induce specific effects, further values. An exposure-response relation for cotton dust.
N Engl J Med 317:605–610
research is needed to compare this pattern of Danuser B, Rebsamen H, Weber C, Krueger H (2000)
concentration changes of inflammatory markers. Lipopolysaccharide-induced nasal cytokine response:
a dose-response evaluation. Eur Arch Otorhino-
Acknowledgments The study was financed by a grant of laryngol 257:527–532
the Lieselotte und Dr. Karl Otto Winkler Stiftung für DECOS (Dutch Expert Committee on Occupational
Arbeitsmedizin and own research funds of the institutions Standards) (1998) Endotoxins: health based
involved (Department of Applied Physiology, Occupa- recommended exposure limit. In: A report of the
tional Medicine and Infectiology, University of Health Council of the Netherlands, publication no
Wuppertal, Germany, and Department of Allergology/ 1998/03WGD. Health Council of the Netherlands,
Immunology, Institute for Prevention and Occupational Rijswijk, pp 1–82
Medicine, Institute of the German Social Accident Insur- Doreswamy V, Alexis NE, Zhou H, Peden DB (2011)
ance, Ruhr University Bochum, Germany). The work of Nasal PMN response to repeated challenge with endo-
the Institute of Occupational and Social Medicine and toxin in healthy volunteers. Inhal Toxicol 23:142–147
Health Services Research (IASV) is supported by an unre- Douwes J, Wouters I, Dubbeld H, van Zwieten L,
stricted grant of the employers’ association of the metal Steerenberg P, Doekes G, Heederik D (2000) Upper
and electric industry Baden-Württemberg airway inflammation assessed by nasal lavage in com-
(Südwestmetall). post workers: a relation with bio-aerosol exposure. Am
The authors thank the Emscher-Genossenschaft/ J Ind Med 37:459–468
Lippeverband (EGLV) for giving access to their sewage Health Council of the Netherlands (2010) Endotoxins.
treatment plants, all volunteers for study participation, and Health-based recommended occupational exposure
the involved staff of IPA, IASV, and University limit. Health Council of the Netherlands, The Hague
Wuppertal, especially Monika Lahr, Heiko Diefenbach, Hoffmeyer F, van Kampen V, Taeger D, Deckert A,
and Nicole Blomberg for technical assistance in collection Rosenkranz N, Kaßen M, Schantora AL, Brüning T,
of specimens and on the sewage treatment plants Raulf M, Bünger J (2014) Prevalence of and relation-
and Gerda Borowitzki, Susanne Freundt, and Sabine ship between rhinoconjunctivitis and lower airway
Bernard for technical assistance in analysis of the NALF diseases in compost workers with current or former
samples. exposure to organic dust. Ann Agric Environ Med
21:705–711
Jacobs JH, Krop EJ, de Wind S, Spithoven J, Heederik DJ
Conflicts of Interest The authors declare no conflicts of (2013) Endotoxin levels in homes and classrooms of
interest in relation to this article. Dutch school children and respiratory health. Eur
Respir J 42:314–322
Jacobs JH, Krop EJ, Borras-Santos A, Zock J, Taubel M,
Hyvarinnen A, Pekkanen J, Doekes G, Heederik DJ
References (2014) Endotoxin levels in settled airborne dust in
European schools: the HITEA school study. Indoor
ATS (1995) American Thoracic Society Standardization Air 24:148–157
of spirometry, 1994 update. Am Rev Respir Dis 152 Kitz R, Rose MA, Placzek K, Schulze J, Zielen S,
(3):1107–1136 Schubert R (2008) LPS inhalation challenge: a new
Bakirci N, Kalaca S, Francis H, Fletcher AM, Pickering C, tool to characterize the inflammatory response in
Anthony C, Tumerdem N, Cali S, Oldham L, Niven R humans. Med Microbiol Immunol 197:13–19
(2007) Natural history and risk factors of early Lai PS, Hang J, Valeri L, Zhang F, Zheng B, Mehta AJ,
respiratory responses to exposure to cotton dust in Shi J, Su L, Brown D, Eisen EA, Christiani DC (2015)
newly exposed workers. J Occup Environ Med Endotoxin and gender modify lung function recovery
49:853–861 after occupational organic dust exposure: a 30-year
BGIA-Arbeitsmappe (1989) Verfahren zur Bestimmung study. Occup Environ Med 72:546–552
der Endotoxinkonzentration in der Luft am Liebers V, Brüning T, Raulf-Heimsoth M (2006) Occupa-
Arbeitsplatz (Kennzahl 9450). In: Berufsgenos- tional endotoxin-exposure and possible health effects
senschaftliches Institut für Arbeitssicherheit (ed) on humans. Am J Ind Med 49:474–491
BIA-Arbeitsmappe Messung von Gefahrstoffen, 28. Liebers V, Raulf-Heimsoth M, Brüning T (2008) Health
Lfg. IV/02. Erich Schmidt Verlag, Bielefeld (Article effects due to endotoxin inhalation (review). Arch
in German) Toxicol 82:203–210
108 M. A. Rieger et al.
Michel O, Nagy A, Schroeven M, Duchateau J, Nève J, recommended threshold limit values for biological
Fondu P, Sergysels R (1997) Dose-response relation- hazards in Germany. J Agric Saf Health 11:185–191
ship to inhaled endotoxin in normal subjects. Am J Rylander R (1995) Endotoxins in the environment. In:
Respir Crit Care Med 156:1157–1164 Lewin J, Alving CR, Munford RS, Redl H (eds) Bac-
Quanjer PH, Tammeling GJ, Cotes JE, Pedersen OF, terial endotoxins, lipopolysaccharids from genes to
Peslin R, Yernault JC (1993) Lung volumes and forced therapy. Proceedings of the third conference of the
ventilatory flows. Eur Respir J 6(Suppl 16):5–40 international endotoxin society. Wiley-Liss Inc,
Quirce S, Lemière C, de Blay F, del Pozo V, van Gerth New York, pp 79–90
WR, Maestrelli P, Pauli G, Pignatti P, Raulf-Heimsoth- Rylander R (2002) Endotoxin in the environment – expo-
M, Sastre J, Storaas T, Moscato G (2010) Noninvasive sure and effects. J Endotoxin Res 8:241–252
methods for assessment of airway inflammation in Rylander R, Bake B, Fischer JJ, Helander IM (1989)
occupational settings. Allergy 65:445–458 Pulmonary function and symptoms after inhalation of
Radon K, Garz S, Riess A, Koops F, Monso E, Weber C, endotoxin. Am Rev Respir Dis 140:981–986
Danuser B, Iversen M, Opravil U, Donham K, Sahlander K, Larsson K, Palmberg L (2012) Daily expo-
Hartung J, Pedersen S, Nowak D (2003) Respiratory sure to dust alters innate immunity. PLoS One 7:
diseases in European farmers-II. Part of the European e31646
farmers’ project. Atemwegs- und Lungenerkrankungen Schantora AL, Casjens S, Deckert A, van Kampen V,
in der Europäischen Landwirtschaft – Teil 2: Neumann HD, Brüning T, Raulf M, Bünger J,
Ergebnisse der europäischen Landwirtschaftsstudie. Hoffmeyer F (2015) Prevalence of work-related
Pneumologie (Stuttgart, Germany) 57:510–517 rhino-conjunctivitis and respiratory symptoms among
Raulf M, Hoffmeyer F, van Kampen V, Deckert A, domestic waste collectors. Adv Exp Med Biol
Brüning T, Bünger J (2015) Cellular and soluble 834:53–61
inflammatory markers in induced sputum of Sigsgaard T, Bonefeld-Jørgensen EC, Kjaergaard SK,
composting plant workers. Adv Exp Med Biol Mamas S, Pedersen OF (2000) Cytokine release from
858:19–29 the nasal mucosa and whole blood after experimental
Raulf M, van Kampen V, Neumann HD, Liebers V, exposures to organic dusts. Eur Respir J 16:140–145
Deckert A, Brüning T, Bünger J, Hoffmeyer F (2017) Sigsgaard T, Bonefeld-Jørgensen EC, Hoffmann HJ,
Airway and blood inflammatory markers in waste Bønløkke J, Krüger T (2005) Microbial cell wall
collectors. Adv Exp Med Biol 1021:1–14 agents as an occupational hazard. Toxicol Appl
Raulf-Heimsoth M, Wirtz C, Papenfuss F, Baur X (2000) Pharmacol 207:310–319
Nasal lavage mediator profile and cellular composition Sigsgaard T, Thorne PS, Schlünssen V, Bønløkke J,
of nasal brushing material during latex challenge tests. Riddervold IS, Hoppe KA, Andersen NT, Mackenzie
Clin Exp Allergy 30:110–121 NM (2015) The change in nasal inflammatory markers
Raulf-Heimsoth M, Hoffmeyer F, van Thriel C, after intranasal challenges with particulate chitin and
Blaszkewicz M, Bünger J, Brüning T (2010) Assess- lipopolysaccharide: a randomized, double-blind, pla-
ment of low dose effects of acute sulphur dioxide cebo-controlled, crossover study with a positive con-
exposure on the airways using non-invasive methods. trol. Int Forum Allergy Rhinol 5:716–723
Arch Toxicol 84:121–127 Spierenburg EA, Smit LA, Heederik D, Robbe P, Hylkema
Raulf-Heimsoth M, Pesch B, Kendzia B, Spickenheuer A, MN, Wouters I (2015) Healthy worker survivor analy-
Bramer R, Marczynski B, Merget R, Brüning T (2011) sis in an occupational cohort study of Dutch agricul-
Irritative effects of vapours and aerosols of bitumen on tural workers. Int Arch Occup Environ Health
the airways assessed by non-invasive methods. Arch 88:1165–1173
Toxicol 85(Suppl 1):S41–S52 Steinberg R (2001) Belastungen und Beanspruchungen
Reynolds SJ, Milton DK, Heederik D, Thorne PS, von Beschäftigten in Arbeitsbereichen der Abwasser-
Donham KJ, Croteau EA, Kelly KM, Douwes J, wirtschaft unter besonderer Berücksichtigung der
Lewis D, Whitmer M, Connaughton I, Koch S, Gefährdung durch biologische Arbeitsstoffe. Gesell-
Malmberg P, Larsson B, Deddens J, Saraf A, Larsson schaft für Abwassertechnik e.V. (GFA), Hennef (Arti-
L (2005) Interlaboratory evaluation of endotoxin cle in German)
analyses in agricultural dusts-comparison of LAL Thorn J (2001) The inflammatory response in humans after
assay and mass spectrometry. J Environ Monit inhalation of bacterial endotoxin: a review. Inflamm
7:1371–1377 Res 50:254–261
Rieger MA, Diefenbach H (2007) Biological stress and Visser MJ, Spaan S, Arts HJ, Smit LA, Heederik DJ
strain in the field of waste water management – how (2006) Influence of different cleaning practices on
important is endotoxin? Gefahrstoffe – Reinhaltung endotoxin exposure at sewage treatment plants. Ann
der Luft 67:363–366 Occup Hyg 50(7):731–736
Rieger MA, Lohmeyer M, Nübling M, Neuhaus S, von Essen S, Romberger D (2003) The respiratory inflam-
Diefenbach H, Hofmann F (2005) A description of matory response to the swine confinement building
the standardized measurement procedures and environment: the adaptation to respiratory exposures
Adaptation to Occupational Exposure to Moderate Endotoxin Concentrations: A. . . 109
in the chronically exposed worker. J Agric Saf Health Wilcoxon F (1945) Individual comparisons by ranking
9:185–196 methods. Biom Bull 1:80–83
Walser SM, Gerstner DG, Brenner B, Bünger J, Wouters IM, Hilhorst SK, Kleppe P, Doekes G, Douwes J,
Eikmann T, Janssen B, Kolb S, Kolk A, Nowak D, Peretz C, Heederik D (2002) Upper airway inflamma-
Raulf M, Sagunski H, Sedlmaier N, Suchenwirth R, tion and respiratory symptoms in domestic waste
Wiesmüller G, Wollin K, Tesseraux I, Herr CE (2015) collectors. Occup Environ Med 59(2):106–112
Evaluation of exposure-response relationships for Zheng Y, Schlünssen V, Bønløkke J, Madsen AM,
health effects of microbial bioaerosols – a systematic Skov S, Sigsgaard T (2014) Change in airway inflam-
review. Int J Hyg Environ Health 218(7):577–589 matory markers in Danish energy plant workers during
Wang Z, Larsson K, Palmberg L, Malmberg P, Larsson P, a working week. Ann Agric Environ Med 21
Larsson L (1997) Inhalation of swine dust induces (3):534–540
cytokine release in the upper and lower airways. Eur
Respir J 10:381–387
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 111–130
https://doi.org/10.1007/5584_2018_234
# Springer International Publishing AG, part of Springer Nature 2018
Published online: 29 June 2018
111
112 G. Leisman et al.
of the above, we investigated behavioral and cog- their parents seeking treatment for the symptoms
nitive changes in ASD as a consequence of the of ASD and from among the individuals who
delivery of red LLLT. responded to recruitment flyers and print ads.
Participants did not receive any financial compen-
sation for participation.
2 Methods
2.2 Inclusion and Exclusion Criteria
2.1 Participants
Each eligible participant satisfied each of the fol-
The study received an approval from the Helsinki
lowing inclusion criteria and none of the exclu-
Committee of the Institute for Neurology and
sion criteria. Males or females aged 5–17 years
Neurosurgery in Havana, Cuba, and was
meeting the criteria of the Diagnostic and Statis-
registered with the FDA (Identifier:
tical Manual of Mental Disorders, 4th Edition, for
NCT03379662). Informed written consent was
ASD, are diagnosed by a qualified medical pro-
obtained from the parent or guardian of each
fessional. Diagnosis was confirmed by the
participant after a full explanation of the
Autism Diagnostic Interview (ADI-R)
procedures to be undertaken. The informed con-
(Constantino et al. 2003). Each participant
sent forms, research protocol, and approvals are
demonstrated “irritable” behaviors such as
available for inspection in the Office of Research
tantrums, aggression, self-injurious behavior, or
Integrity at the Institute of Neurology and Neuro-
a combination of thereof. The participant’s Aber-
surgery in Havana, Cuba.
rant Behavior Checklist (ABC) irritability sub-
The participants consisted of 40 individuals
scale (Rojahn et al. 2003) score was >18; the
distributed at baseline in the fashion described in
Clinical Global Impressions Severity Scale
Tables 1, 2, and 3. Participants in both groups
(CGI-S) (Berk et al. 2008) score was 4
spanned 5–16 years, with the mean participant
(4 ¼ moderately ill). The participants’ therapeu-
age of approx. 8 years. A t-test for independent
tic/intervention plan had been consistent/stable
samples revealed no statistically significant dif-
for 3 months. They abstained from undertaking
ference in age between test and placebo group
new treatments during the study time.
participants (μa-μb ¼ 0.13; t ¼ 0.13;
Exclusion criteria were the following: history
df ¼ 38; p(two-tailed) ¼ 0.90 ( p > 0.05)).
of a primary or concurrent diagnosis of another
All study participants were recruited from
disorder, including neurological, use of a psycho-
among the investigator’s normal pool of new
tropic drug, or any participation in a research
and existing patients who voluntarily came with
study within 30 days prior to the commencement
of the current study.
Table 1 Age of study participants by the procedure group
Age (years) Testa (n ¼ 21) Placebo (n ¼ 19)
Mean 8.2 8.4 2.3 Baseline Concomitant
Standard deviation 3.0 3.2 Medications and Low-Level Light
Range 5–16 5–16 Therapy (LLLT)
a
Test group was the active treatment group
Table 4 presents the OTC and prescription
Table 2 Participant gender breakdown for the test and medications reported by participants at the base-
placebo group participants line evaluation used in the past by both test and
Gender Testa (n ¼ 21) Placebo (n ¼ 19) placebo group participants. No participant took
Male 16 14 any OTC medication during the study time.
Female 5 5 Table 5 lists the non-drug therapies routinely
a
Test group was the active treatment group used by both test and placebo group participants
114 G. Leisman et al.
Table 4 Prior medication used to treat the symptoms of per week, 3–4 days apart at the investigator’s test
autistic disorder by the procedure group site. A pulsed laser of 635 nm with a power output
Medicationa Testb (n ¼ 21) Placebo (n ¼ 19) of 15 mW and a red 635 nm LED were used as
Carbamazepine 10 7 treatment and placebo, respectively.
Risperdal 9 8 Participants were required to maintain their
Valproate 2 3 regular medication schedule and treatment
Methylphenidate 3 1 regimens, as reported at the baseline evaluation,
Haloperidol 1 3 to treat symptoms related to autistic disorder
Chlorpromazine 1 2 throughout the study time. All of them complied
Risperidone 1 1 with this requirement.
Conductosa 1 1
Thioridazine 0 2
Levomepromazine 0 1
a 2.5 Outcome Measures
Some participants previously used several medications;
therefore, medication use by participant group adds up to
greater than the participant sample size per group; btest Pre-post treatment outcome measures included
group was the active treatment group the Aberrant Behavior Checklist (ABC). The
global score and the five subscale scores
to treat the symptoms of ASD at the time of entry consisted of (a) irritability and agitation,
to the study. Therapy used to treat symptoms of (b) lethargy and social withdrawal,
ASD was minimal and equal between test and (c) stereotypic behavior, (d) hyperactivity and
placebo group participants. noncompliance, and (e) inappropriate speech.
Medications routinely used by both test and The ABC is an informant rating instrument
placebo group participants at the time of study empirically derived by principal component anal-
enrollment for indications other than to treat the ysis. The global score for the ABC has not been
symptoms of autistic disorder are listed in psychometrically derived and is not statistically
Table 6. valid (Farmer and Aman 2012).
The second outcome measure consisted of the
Clinical Global Impressions (CGI) scale includ-
2.4 Procedure ing a severity-of-illness scale (CGI-S) and a
global improvement/change scale (CGI-C).
All the 40 participants completed the course Leucht and Engel (2006) have noted that both
according to the protocol. Twenty-one of them the CGI and the Brief Psychiatric Rating Scale
were randomized to the test (active treatment) (BPRS) are often employed in drug trials. Those
procedure group, and 19 were randomized to the authors have found that by calculating the effect
placebo procedure group. size and its 95% confidence interval for both
Participants received eight 5-min laser light continuous (standardized mean differences) and
applications to the base of the skull and temporal dichotomous (odds ratio) outcomes, there was no
areas with the Erchonia® EAL Laser (active or significant differences between tests, which
sham) across a 4-week period: two applications indicates a good inter-test reliability.
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 115
Table 5 Therapies used to treat the symptoms of autistic disorder by the procedure group
Therapy Testa (n ¼ 21) Placebo (n ¼ 19)
Logopedic therapy for language 2 2
Hydrotherapy for hyperactivity 1 0
a
Test group was the active treatment group
The ABC and CGI-S and CGI-C were 3.2 Primary Efficacy Outcome
administered prior to treatment and then 2 (mid- Analysis
point) and 4 weeks (endpoint) of treatment and
finally 8 weeks after treatment (post-procedure The evaluation time point at which study success
evaluation). was analyzed (study endpoint) was predetermined
as 4 weeks following baseline (pretreatment)
evaluation. The study was predetermined to be
considered successful if, using the intent-to-treat
3 Results
(ITT) last observation carried forward (LOCF)
analysis, the primary endpoint was statistically
3.1 Pre-procedure Measures
significant at the 0.05 level. The primary efficacy
outcome measure was predefined as the mean
change from baseline to endpoint (end of the
Aberrant Behavior Checklist (ABC) Mean and
4-week treatment period) in the ABC irritability
standard deviation baseline (pre-procedure) of
and agitation subscale score, with a minimum
the ABC global and subscale scores for test and
mean difference of 8.5 points between the test
placebo group participants are provided in
and placebo groups.
Table 7.
The primary outcome measure was evaluated
for the following two populations:
T-tests for independent samples revealed no sta-
tistically significant differences in the baseline
1. Intent-to-Treat (ITT). ITT analysis included all
ABC global score or in any of the five ABC
participants with valid measurements at base-
subscale scores between test and placebo group
line, randomized to a procedure group.
participants ( p > 0.05) as shown in Table 8.
Dropouts and terminated participants were
handled by carrying forward the last observa-
Clinical Global Impressions Severity (CGI-S)
tion for all time points following the dropout
Score Table 9 shows the CGI-S baseline score
(LOCF). If a participant was not a dropout but
for the test and placebo group participants. The
had no data in a relative day range, the last
majority of participants in both groups had a
observation prior to the time point being
baseline CGI-S score of 6, corresponding to
analyzed was employed.
“severely ill”.
2. Per-Protocol Population. Per-protocol analy-
sis of results intended to corroborate
116 G. Leisman et al.
Table 8 T-test results for differences in ABC global and subscale scores between the two procedure group
(test vs. placebo)
μa-μb t df p
Global score 0.95 +0.44 38 0.66
Irritability and agitation 1.64 0.68 38 0.50
Lethargy and social withdrawal 1.40 +0.91 38 0.37
Stereotypic behavior 4.09 1.64 38 0.11
Hyperactivity and noncompliance 0.77 +0.69 38 0.49
Inappropriate speech 2.60 +0.33 38 0.74
conclusions drawn from the ITT analysis. This adjusted means for the absolute change in ABC
analysis excluded the participants with major irritability subscale score from baseline to end-
protocol deviations and incompletes point, adjusting for the covariate of baseline ABC
(dropouts, non-compliant participants, irritability subscale scores.
disqualified participants, etc.). All participants The adjusted mean difference in the baseline to
in this study completed all study visits and endpoint change in the ABC between the test and
procedures and had all study measurements placebo groups was 15.2, almost double the
recorded through to the study endpoint evalu- pre-established study success criteria of 8.5
ation. Therefore LOCF was not applied. points between the procedure groups, in favor of
the test group (F ¼ 99.3; p < 0.0001). Table 11
shows the observed mean and standard deviation
at baseline and endpoint for the ABC score and
change between the evaluations by the procedure
3.3 General Statistical Evaluation
group.
The ABC score decreased 14.8 points for
For the ABC irritability subscale score, changes
participants in the test group, while it increased
were evaluated by analysis of covariance
by 0.32 points for participants in the placebo
(ANCOVA), with baseline measures as a covari-
group, resulting in an observed difference
ate and a procedure group (active or placebo) as
between the two groups in the ABC irritability
the main effect. Table 10 shows the observed and
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 117
Table 10 Change in ABC irritability subscale score from baseline to endpoint, adjusting for baseline ABC irritability
subscale score
ABC irritability subscale score Testa (n ¼ 21) Placebo (n ¼ 19)
Observed mean 14.81 0.32
Adjusted mean 14.83 0.34
a
Test group was the active treatment group
Table 11 ABC irritability subscale score from baseline to endpoint by the procedure group
Testa (n ¼ 21) Placebo (n ¼ 19)
ABC irritability subscale score Mean SD Mean SD
Baseline 30.5 6.7 29.6 6.8
Endpoint 15.7 9.9 29.9 6.6
Change 14.8 6.4 0.3 1.4
a
Test group was the active treatment group
subscale score from baseline to endpoint evalua- on the participant’s primary caregiver’s rating and
tion of 15.1. (b) CGI-I rating of 1 (very much improved) or
2 (much improved) at study endpoint as deter-
T-Test for Two Correlated Samples A correlated mined by the clinician’s evaluation.
sample t-test was used to compare the difference Participants attaining a 25% reduction
in the mean change in ABC score from study from baseline to endpoint in the ABC irritability
baseline to endpoint. The evaluation for partici- subscale score. Eighteen (86%) of the 21 test
pant groups separately found the mean change to group participants attained a minimum of 25%
be significant for the test group (t ¼ +10.6; reduction in the ABC score from baseline to end-
df ¼ 20; p < 0.001), but not for the placebo point, while none of the 19 placebo participants
group (t ¼ 1.0; df ¼ 18, p ¼ 0.33). did. Table 12 shows the mean and standard devi-
ation of the percent (%) change in ABC scores
T-Test for Independent Samples A t-test for two from baseline to endpoint for the test and placebo
independent samples, used to compare the mean groups.
ABC score change from baseline to endpoint
between the test and placebo groups, found that ABC Irritability Subscale A t-test for two inde-
the mean difference of 15.1 was significant in pendent samples was performed to compare the
favor of the participants in the test group mean difference in percent change in ABC scores
(t ¼ 10.1; df ¼ 38; p < 0.0001). from baseline to endpoint between the test group
and placebo groups. The mean difference of
-52.5% was significant in favor of the test group
participants (t ¼- 9.79; df ¼ 38; p < 0.0001).
3.4 Secondary Efficacy Outcome
Analysis
CGI-I Ratings Seventeen (81%) of the 21 test
group participants received a CGI-I rating of
3.4.1 Positive Responder Rate (PRR)
1 or 2 at the study endpoint evaluation, while
The evaluation of the difference in PRR between
none of the 19 placebo participants did.
the procedure groups was performed as a
predetermined secondary efficacy outcome mea-
Positive Responder Rate Successes Table 13
sure to provide support for the primary efficacy
shows the number and percentage of test and
analysis. The PRR was defined as satisfaction
placebo group participants meeting the dual
with both of the following: (a) 25% reduction
criteria for PRR. There was a difference of 81%
from baseline to endpoint in the ABC score based
118 G. Leisman et al.
Table 12 Percentage change in ABC irritability subscale score by the procedure group
% change in ABC irritability score Testa (n ¼ 21) Placebo (n ¼ 19)
Mean SD 51.0 22.4 1.5 6.6
a
Test group was the active treatment group
Table 14 Comparison of the proportion of positive responder rates between the procedure groups
2 2 table Positive responder rate (PRR) Negative responder rate (NRR)
Testa group 17 4
Placebo group 0 19
p < 0.00001 17 23
a
Test group was the active treatment group
in the proportion of participants who met the PRR score as a covariate and the group as a main effect
between procedure groups: 81% of test group (Tables 16 and 17). The F values were statistically
participants met the PRR criteria compared with greater for the mean change for each of the ABC
none in the placebo group. Fisher’s exact test for scores between groups ( p < 0.0001). With
two independent proportions was conducted to removal of individual differences in baseline
compare the proportion of PRR between the scores, the two adjusted means significantly dif-
groups (Table 14). The difference in this pro- fered ( p < 0.0001), and the means were signifi-
portion was significant ( p < 0.00001). The cant for the test but not for placebo group
greater PRR from baseline to endpoint for the participants (Tables 18 and 19).
test group relative to the placebo group also was A t-test analysis for two independent samples
significant. comparing differences in the mean change in the
ABC scores from baseline to endpoint between
3.4.2 Aberrant Behavior Checklist (ABC) the two groups is provided in Table 20; the dif-
Global and Subscale Scores ference was significantly greater for the test than
The evaluation of the mean change from baseline placebo group participants.
to endpoint for each of the ABC global score and The evaluation of the mean change in each of
the remaining four ABC subscale scores (leth- the ABC global score and the five subscale sores
argy/social withdrawal, stereotypic behavior, (irritability/agitation, lethargy/social withdrawal,
hyperactivity/noncompliance, and inappropriate stereotypic behavior, hyperactivity/noncompli-
speech) between the two procedure groups was ance, and inappropriate speech) across the study
performed as a predetermined secondary efficacy duration between the procedure groups was
outcome measure to provide support for the pri- performed as a predetermined supportive measure
mary efficacy analysis outcome. Table 15 shows providing support for the primary efficacy analy-
the mean and standard deviation of baseline and sis. The ABC scores were recorded at the follow-
endpoint ABC global and subscale scores and the ing four evaluation points: baseline
changes for the test and placebo participants (pre-procedure), week 2 (interim procedure
(Aman et al. 1985a, 1985b). The ABC global administration), week 4 (study endpoint), and
and subscale scores were evaluated with week 8 (i.e., 4 weeks post-procedure). Table 21
ANCOVA, with the baseline measure for each shows these data for the test and placebo groups.
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 119
Table 15 ABC global and subscale scores from baseline to endpoint by the procedure group
Testa (n ¼ 21) Placebo (n ¼ 19)
Mean SD Mean SD
Global score Baseline 107.3 20.3 104.7 28.7
Endpoint 63.8 30.5 105.4 28.4
Change 43.5 19.1 0.7 2.6
Lethargy and social withdrawal Baseline 23.1 9.3 24.7 5.1
Endpoint 13.8 8.8 24.7 5.1
Change 9.3 5.8 0.1 0.2
Stereotypic behavior Baseline 13.7 4.1 12.3 5.6
Endpoint 8.2 5.1 12.3 5.6
Change 5.5 4.0 0.0 0.0
Hyperactivity and noncompliance Baseline 32.8 7.8 36.9 7.9
Endpoint 21.1 9.5 37.3 7.4
Change 11.7 7.5 0.4 1.1
Inappropriate speech Baseline 7.2 3.1 6.4 4.0
Endpoint 4.9 2.4 6.4 3.9
Change 2.3 2.3 0.0 0.3
a
Test group was the active treatment group
Table 16 Change in ABC scores from baseline to endpoint, adjusted for baseline score
ABC global and subscales Mean scores Testa (n ¼ 21) Placebo (n ¼ 19)
Global score Observed 43.53 0.74
Adjusted 43.60 0.82
Lethargy and social withdrawal Observed 9.29 0.05
Adjusted 9.44 0.22
Stereotypic behavior Observed 5.47 0
Adjusted 5.43 0.05
Hyperactivity and noncompliance Observed 11.67 0.37
Adjusted 11.95 0.68
Inappropriate speech Observed 2.29 0
Adjusted 2.21 0.08
a
Test group was the active treatment group
Table 17 Significance of the mean baseline to endpoint change in the ABC global and subscale scores between the test
and placebo groups
Observed difference Adjusted difference F p
Global score 44.3 44.4 99.0 <0.0001
Lethargy and social withdrawal 9.7 9.3 58.1 <0.001
Stereotypic behavior 5.5 5.4 33.8 0.000001
<0.00005
Hyperactivity and noncompliance 12.0 12.6 50.3 <0.0001
Inappropriate speech 2.3 2.1 19.3 0.000091
<0.0001
Figures 1, 2, 3, 4, 5, and 6 show the mean indicate that each ABC score in the test (active
changes in each of the ABC global and five sub- treatment group) decreased significantly from the
scale scores across the study duration by the pro- baseline level across the three evaluation points;
cedure group. The figures as well as Table 21 the decrease was progressively augmented over
120 G. Leisman et al.
Table 18 t-test analysis for the mean change in the ABC global and subscale scores for the test group participants
t df p
Global score +10.45 20 <0.0001
Lethargy and social withdrawal +7.34 20 <0.0001
Stereotypic behavior +6.34 20 <0.0001
Hyperactivity and noncompliance +7.14 20 <0.0005
Inappropriate speech +4.59 20 <0.0001
Table 19 t-test analysis for the mean change in the ABC placebo groups was performed as a
global and subscale scores for the placebo group predetermined secondary efficacy outcome mea-
participants
sure providing support for the primary efficacy
t df p analysis outcome. Table 24 shows the number of
Global score 1.25 18 0.23 participants by the CGI-S rating at baseline and
Lethargy and social withdrawal 1.00 18 0.33 endpoint for both groups of participants. All
Stereotypic behavior –0.0 18 –1.00
21 test (active treatment) participants showed a
Hyperactivity and noncompliance 1.44 18 0.17
one-category or greater improvement in CGI-S
Inappropriate speech 0.0 18 1.00
ratings from baseline to endpoint. In contradis-
tinction, the majority of placebo participants
Table 20 t-test analysis for the mean change in the ABC (17) showed no change in CGI-S rating, one
scores between the test and placebo groups demonstrated a worsening of one category, and
t df p another demonstrated an improvement of
Global score 10.01 38 <0.0001 one-category rating from baseline to endpoint. A
Lethargy and social 7.01 38 <0.0001 2 3 Fisher’s exact test for matched categorical
withdrawal data was conducted to compare the proportion of
Stereotypic behavior 6.02 38 <0.0001 participants whose CGI-S rating improved,
Hyperactivity and 6.92 38 <0.0005 showed no change, or worsened from baseline
noncompliance
to endpoint in the test and placebo groups
Inappropriate speech 4.32 38 <0.0001
(Table 25).
A difference in the proportion of participants
time, including a 4-week follow-up period during whose CGI-S rating changed one or more
which no further treatment occurred. Conversely, categories from baseline to endpoint between the
there was no such change in the placebo group. test and placebo groups was significant
One-way ANOVA and Tukey’s HSD post hoc ( p < 0.00001), with a greater proportion of test
analyses were performed to evaluate the mean group participants who had improved CGI-S
change in each ABC global and subscale scores ratings from baseline to endpoint relative to the
across the four study evaluation points in the test placebo group (Table 26). In both test and pla-
and placebo groups. Table 22 shows that the score cebo groups, all participants (100%) were rated at
improvement in each of the ABC subscales was baseline in the top three severity-of-condition
significant across almost all the comparative ratings, i.e., marked, severe, and most extreme.
evaluations in the test group. Conversely, score All placebo group participants (100%) retained
changes in the placebo participants were not sta- these top CGI-S ratings (categories 5–7 in
tistically significant (Table 23). Table 26) by week 8, while only 3 of the 21 of
active test treatment (14%) received the same top
3.4.3 Clinical Global Impressions ratings. This continuously progressive and sub-
Severity-of-Illness (CGI-S) Ratings stantial improvement for the test over placebo
The evaluation of the change in CGI-S ratings group participants was illustrated in Fig. 7 which
from baseline to endpoint between the test and indicates the percentage of participants who
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 121
Table 21 ABC global and subscale scores across the study duration by the test and placebo groups
Testa (n ¼ 21) Placebo (n ¼ 19)
Mean SD Mean SD
Global score Baseline 107.3 20.3 104.7 28.7
Week 2 89.3 24.1 110.0 19.4
Week 4 63.8 30.5 105.4 28.4
Week 8 48.8 25.1 105.3 31.6
Irritability Baseline 30.5 6.7 29.6 6.8
Week 2 23.4 7.4 29.9 7.5
Week 4 15.7 9.9 29.9 6.5
Week 8 11.8 10.0 29.9 6.6
Lethargy and social withdrawal Baseline 23.1 9.3 24.7 5.1
Week 2 18.2 9.3 24.7 5.1
Week 4 13.8 8.8 24.7 5.1
Week 8 10.7 9.8 24.7 5.1
Stereotypic behavior Baseline 13.7 4.1 12.3 5.6
Week 2 11.6 4.4 12.3 5.6
Week 4 8.2 5.1 12.3 5.6
Week 8 6.9 5.1 12.6 5.7
Hyperactivity and noncompliance Baseline 32.8 7.8 36.9 7.9
Week 2 30.3 8.2 36.7 7.8
Week 4 21.1 9.5 37.3 7.4
Week 8 16.0 10.5 37.3 7.4
Inappropriate speech Baseline 7.2 3.1 6.4 4.0
Week 2 5.8 2.6 6.4 4.0
Week 4 4.9 2.4 6.4 3.9
Week 8 3.5 2.2 6.5 3.9
a
Test group was the active treatment group
received the 5–7 categories of CGI-S ranking relative to baseline, with the majority (13) receiv-
across the study duration by the procedure group. ing a rating of “much improved”. Conversely, no
placebo group participant was rated as
3.4.4 Clinical Global Impressions demonstrating symptom improvement at end-
Improvement/Change (CGI-C) point relative to baseline, with the majority
Ratings (17) rated as “no change” and the remaining
The evaluation of CGI-C ratings at study end- 2 as “minimally worse”. A proportion of
point between the procedure groups was participants in both test and placebo groups
performed as a predetermined secondary efficacy whose CGI-C rating at study endpoint, relative
outcome measure to provide support for the pri- to baseline, showed “improvement” (CGS-C
mary efficacy analysis outcome. The number of rating of 1, 2, and 3), “no change” (CGI-C rating
participants by the CGI-C rating at endpoint for of 4), or “worsening” (CGI-C rating of 5, 6, and
the test and placebo groups is shown in Table 27. 7) in symptoms, evaluated with 2 x 3 Fisher’s
Twenty out of the 21 test group participants exact test, is shown in Table 28. A difference in
showed some degree of improvement in the this proportion between the two groups was sig-
autism-related symptoms at study endpoint nificant ( p < 0.00001). There was a greater
122 G. Leisman et al.
Points
group participants 80
60
40
Baseline Week 2 Week 4 Week 8
Test Placebo
Test Placebo
10
Baseline Week 2 Week 4 Week 8
Test Placebo
6
Baseline Week 2 Week 4 Week 8
Test Placebo
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 123
Points
noncompliance subscale 30
score in the test (active
25
treatment) and placebo
group participants 20
15
Baseline Week 2 Week 4 Week 8
Test Placebo
Test Placebo
Table 22 One-way ANOVA and Tukey’s HSD post hoc follow-up results for the change in the Aberrant Behavior
Checklist (ABC) score across the study duration in the test (active treatment) group participants
ABC scale/subscale F p Tukey’s HSD
Global score 101.6 <0.0001 Changes between each possible evaluation
Combination was significant at p < 0.01
Irritability 114.5 <0.0001 Changes between each possible evaluation
Combination was significant at p < 0.01
Lethargy and social 49.6 <0.0001 Changes between each possible evaluation
withdrawal Combination was significant at p < 0.01, except week 4–8 which was
significant at p < 0.05
Stereotypic behavior 38.2 <0.0001 Baseline to week 2 significant at p < 0.05
All others significant at p < 0.01, except week 4–8 which was not
significant
Hyperactivity and 59.4 <0.0001 All significant at p < 0.01, except baseline to week 2 which was not
noncompliance significant
Inappropriate speech 26.9 <0.0001 All significant at p < 0.01, except week 2–4 which was not significant
124 G. Leisman et al.
Table 23 One-way ANOVA results for the change in Aberrant Behavior Checklist (ABC) score across the study
duration in the placebo group participants
ABC scale/subscale F p
Global score 0.4 0.789
Irritability 0.4 0.739
Lethargy and social withdrawal 1.0 0.399
Stereotypic behavior 1.0 0.399
Hyperactivity and noncompliance 2.4 0.081
Inappropriate speech 0.4 0.739
Table 24 Clinical Global Impressions severity-of-illness (CGI-S) ratings at baseline and endpoint by the procedure
group
Testa (n ¼ 21) Placebo (n ¼ 19)
CGI-S category Baseline Endpoint Baseline Endpoint
7: Most extreme of all autism patients 3 – – –
6: Severe 14 2 15 16
5: Marked 4 1 4 3
4: Moderate – 13 – –
3: Mild – 4 – –
2: Borderline – 1 – –
1: Normal – – – –
a
Test group was the active treatment group
Table 25 Fisher’s exact test for matched categorical data to compare the proportion of participants whose Clinical
Global Impressions severity-of-illness (CGI-S) ratings improved
2 3 table Improved Unchanged Worsened n
Testa group 21 0 0 21
Placebo group 1 17 1 19
Total 22 17 1 40
a
Test group was the active treatment group; p < 0.00001
Table 26 Clinical Global Impressions severity-of-illness (CGI-S) ratings across the study duration by the procedure
groups
Testa (n ¼ 21) Placebo (n ¼ 19)
CGI-S Baseline Week 2 Week 4 Week 8 Baseline Week 2 Week 4 Week 8
7: Most extreme 3 2 – 1 – – – –
6: Severe 14 1 2 2 15 15 16 15
5: Marked 4 15 1 – 4 4 3 4
4: Moderate – 3 13 9 – – – –
3: Mild – – 4 6 – – – –
2: Borderline – – 1 3 – – – –
1: Normal – – – – – – – –
a
Test group was the active treatment group; p < 0.00001
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 125
%
100
80
60
40
20
0
Test Placebo
Fig. 7 Percentage of participants receiving a 5–7 rating in the severity-of-illness subscale of the Clinical Global
Impressions Scale (CGI-S) across the study duration
Table 27 Clinical Global Impressions improvement/change subscale (CGI-C) ratings at the study endpoint by the
procedure group
CGI-C category Testa (n ¼ 21) Placebo (n ¼ 19)
1: Very much improved 4 –
2: Much improved 13 –
3: Minimally improved 3 –
4: No change 1 17
5: Minimally worse – 2
6: Much worse – –
7: Very much worse – –
a
Test group was the active treatment group; p < 0.00001
Table 28 Fisher’s exact test comparing the proportion of participants with Clinical Global Impressions improvement/
change subscale (CGI-C) ratings of ASD symptoms at the study endpoint
2 3 table Improved Unchanged Worsened n
Test groupa 20 1 0 21
Placebo group 0 17 2 19
Total 20 18 2 40
a
Test group was the active treatment group; p < 0.00001
Table 29 Clinical Global Impressions improvement/change subscale (CGI-C) ratings across the study duration by the
procedure group
Testa (n ¼ 21) Placebo (n ¼ 19)
CGI-C rating Week 2 Week 4 Week 8 Week 2 Week 4 Week 8
1: Very much improved – 5 10 – – –
2: Much improved 9 12 8 – – –
3: Minimally improved 11 3 3 – – –
4: No change 1 1 – 18 17 17
5: Minimally worse – – – 1 2 2
6: Much worse – – – – – –
7: Very much worse – – – – – –
a
Test group was the active treatment group; p < 0.00001
indicated a continuous progressive improvement LLLT is more effective than placebo in effecting
in symptoms. As early as 2 weeks into the 4-week a positive change in the ASD symptoms.
procedure, 20 out of the 21 participants received All participants of the active treatment showed
CGI-C ratings of improvement, relative to base- a one-category or greater improvement in CGI-S
line, with 11 showing “minimal improvement” from baseline to endpoint, while 17 of the placebo
and 9 “much improvement”. By week 8, the pre- participants showed no CGI-S rating change. The
sentation of autistic symptoms and behaviors was difference in the proportion of participants whose
reported as “very much improved” for almost half CGI-S rating changed one or more categories
(48%) of active treatment participants. Con- from baseline to endpoint between the two groups
versely, there was essentially no change in the was significant ( p < 0.00001). This difference is
CGI-C ratings across the study duration for the attributable to the efficacy of LLLT compared
placebo participants, indicating no observable with placebo. All 21 active treatment participants
symptomatic improvement. showed a one-category or greater improvement in
CGI-S from baseline to endpoint.
Twenty out of the 21 active treatment
participants showed some degree of improvement
4 Discussion
in autism-related symptoms at endpoint relative to
baseline. The majority (13) received ratings of
Baseline ABC irritability subscale scores
“much improved”. No placebo group participant
indicated that LLLT was significantly more effec-
demonstrated improvement in symptoms at end-
tive than placebo in treating symptoms of ASD in
point relative to baseline. The majority
children and adolescents independent of baseline
(17) demonstrated “no change”, and the
ABC scores, which accounted for 6% of the vari-
remaining 2 placebo participants rated “mini-
ance in the ABC score change from baseline to
mally worse”. We found that the ABC global
endpoint. The difference in the proportion of PPR
and five subscale scores decreased progressively
between the LLLT and placebo groups was sta-
and significantly from baseline across each of the
tistically significant ( p < 0.00001).
three successive evaluation points; the decrease
Considering the baseline score as a covariate,
progressed over time, including a 4-week follow-
F values were statistically significant
up during which no further LLLT occurred. Con-
( p < 0.0001) for the mean change from baseline
versely, the placebo group demonstrated no sig-
to endpoint between the LLLT and placebo
nificant change across the study duration,
groups and within participant groups for both
demonstrating the effectiveness of LLLT in
ABC global and subscale scores. Removing indi-
reducing ASD-associated symptoms.
vidual differences in the baseline score, the
For both test and placebo group, all
adjusted means significantly differed
participants were rated at baseline in the top
( p < 0.0001), supporting the conclusion that
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 127
three severity conditions. By week 8, all placebo 2013). Abnormal anatomical connectivity and
participants retained a 5–7 CGI-S score, while performance of hub regions have been
only 3 of the 21 active treatment (14%) received hypothesized to be associated with cognitive and
this rating, which speaks for a progressive behavioral dysfunction in numerous neurological
improvement in the symptom presentation across and psychiatric disorders (Bullmore and Sporns
the study time. As early as 2 weeks into the 2012; Seeley et al. 2009). For example, analyses
4-week active treatment, 20 out of the 21 test of structural and functional connectivity in
group participants received improved CGI-C schizophrenia have shown reduced frontal hub
ratings, with 11 participants showing “minimal” connectivity (Fornito et al. 2012).
and 9 demonstrating “much improvement”. By We have recently noted that altered functional
week 8, symptom presentation of ASD was connectivity, i.e., synchronous brain activity,
reported as being “very much improved” for might be associated with the deficits characteris-
almost half (48%) of the test group participants, tically found in ASD (Machado et al. 2015). Of
whereas no significant change in the CGI-C rating specific importance is the integrity of functional
was noted for the placebo group participants. connectivity in the default mode network (DMN),
There is a well-established literature on a network active during inactive resting states,
photobiomodulation, supporting improvement in and in cognitive functions linked to the
dysfunctional neuronal activity with the use of ASD-related social dysfunction. Assaf et al.
low intensity red and near-infrared (NIR) light (2010) have found a decreased functional connec-
(Hiwaki and Miyaguchi 2018; Naeser et al. tivity between the precuneus/medial/prefrontal/
2014). Ten-year experience in animal studies anterior/cingulate cortices and default mode net-
has indicated that transcranial photobiomo- work core areas, with the degree of functional
dulation has a positive effect in animal models connectivity in these regions inversely correlated
of traumatic brain injury (Wu et al. 2012a; Oron with the ASD communication and social deficits.
et al. 2007), Alzheimer’s (Wu et al. 2012b), These results support the hypothesis that default
depression (Purushothuman et al. 2015), and mode networks’ under-connectivity contributes
stroke (Hiwaki and Miyaguchi 2018; Lapchak to the core deficits seen in ASD.
and deTaboada 2010; Oron et al. 2006), while Of particular interest is the association of
human studies have included traumatic brain default mode networks with ASD (Buckner
injury (Naeser et al. 2011), depression (Schiffer et al. 2008). Numerous studies have found vari-
et al. 2009), and stroke (Zivin et al. 2009). Fur- ous “hubs” with the effectiveness of the entire
ther, low-level light energy has been found safe network highly dependent on the hub status. The
for humans in the stroke studies, providing a high most common ASD-related hubs are those
benefit-to-risk ratio, with no reported side effects reported by Raichle et al. (2001) and Greicius
of LLLT. et al. (2003) being associated with the ventrome-
Almost all aspects of cognitive function dial and dorsomedial prefrontal, posterior cingu-
require integration of widely distributed neural late, precuneus, lateral parietal, and entorhinal
activity. Network analysis of human brain con- cortices (Raichle 2015). Assaf et al. (2010) have
nectivity has reliably classified sets of neocortical examined three default mode sub-networks
areas essential for supporting optimized neuronal obtained from the resting fMRI scans of
signaling and communication. Optimal “brain 16 ASD individuals and 16 matched controls,
hubs” exist in the networks of effective connec- using independent component analysis. The
tivity, supporting the notion of the hub’s function ASD individuals demonstrate a reduced func-
in a wide range of cognitive tasks and a dynamic tional connectivity between the medial-prefron-
coupling in and across effective functional tal/anterior-cingulate/precuneus cortices and
networks. A high level of brain hub centrality other default mode sub-networks areas. The
renders these networks vulnerable to dysfunction degree of functional connectivity in these regions
and disconnection (van den Heuvel and Sporns
128 G. Leisman et al.
is inversely associated with the severity of social employed in this study. The authors would like to thank
and communication deficits. Ms. Elvira Walls for her assistance with the statistical
analysis and research methodology.
LLLT promotes cell and neuronal repair
(Dawood and Salman 2013) and brain network
Conflicts of Interest None of the authors have any fidu-
rearrangement (Erlicher et al. 2002) in many neu- ciary interests in the research reported herein or with the
rologic disorders identified with lesions in the source of funding.
hubs of default mode networks (Buckner et al.
2008). LLLT facilitates a fast-track wound-
healing (Dawood and Salman 2013) as
mitochondria respond to light in the red and References
near-infrared spectrum (Quirk and Whelan
Aman MG, Singh NN, Stewart AW, Field CJ (1985a) The
2011). On the other hand, Erlicher et al. (2002) aberrant behavior checklist: a behavior rating scale for
have demonstrated that weak light directs the the assessment of treatment effects. Am J Ment Defic
leading edge of growth cones of a nerve. There- 89:485–491
Aman MG, Singh NN, Stewart AW, Field CJ (1985b)
fore, when a light beam is positioned in front of a
Psychometric characteristics of the aberrant behavior
nerve’s leading edge, the neuron will move in the checklist. Am J Ment Defic 89:492–502
direction of the light and grow in length (Black Assaf M, Jagannathan K, Calhoun VD, Miller L, Stevens
et al. 2013; Quirk and Whelan 2011). Nerve cells MC, Sahl R, O'Boyle JG, Schultz RT, Pearlson GD
(2010) Abnormal functional connectivity of default
appear to thrive and grow in the presence of
mode sub-networks in autism spectrum disorder
low-energy light, and we think that the effect patients. NeuroImage 53:247–256
seen here is associated with the rearrangement Barrett LF, Satpute AB (2013) Large-scale brain networks
of connectivity. in affective and social neuroscience: towards an inte-
grative functional architecture of the brain. Curr Opin
Reports are now emerging that LLLT and
Neurobiol 23:361–372
photobiomodulation significantly upregulate Behrmann M, Plaut DC (2013) Distributed circuits, not
brain-derived neurotrophic factor (BDNF), a fac- circumscribed centers, mediate visual recognition.
tor highly associated with dendritic sprouting, Trends Cogn Sci 17:210–219
Berk M, Ng F, Dodd S, Callaly T, Campbell S,
neuroplasticity, and brain reconnectivity
Bernardo M, Trauer T (2008) The validity of the CGI
(Hamblin 2018). In particular, photobiomo- severity and improvement scales as measures of clini-
dulation influences the ubiquitous transcription cal effectiveness suitable for routine clinical use. J Eval
factor CREB (cAMP-response element-binding Clin Pract 14:979–983
Black B, Mondal A, Kim Y, Mohanty SK (2013) Neuronal
protein) to increase BDNF and consequently den-
Beacon. Opt Soc Amer Optics Lett 38:2174–2176
drite growth (length, spine density, and branching Brovelli A, Ding M, Ledberg A, Chen Y, Nakamura R,
in the hippocampal neurons). Functional mag- Bressler SL (2004) Beta oscillations in a large-scale
netic resonance studies demonstrate the activity sensorimotor cortical network: directional influences
revealed by granger causality. Proc Nat Acad Sci
modulation in intrinsic brain networks, including
USA 101:9849–9854
default mode networks, likely to be dysfunctional Buckner RL, Andrews-Hanna JR, Schacter DL (2008) The
in a variety of conditions. All of that supports the brain’s default network: anatomy, function, and rele-
notion that LLLT can modify functional and vance to disease. Ann N Y Acad Sci 1124:1–38
Bullmore E, Sporns O (2012) The economy of brain net-
effective connectivity in the neocortex. Research
work organization. Nat Rev Neurosci 13:336–349
on photobiomodulation reveals the beneficial Constantino JN, Davis SA, Todd RD, Schindler MK,
effects of LLLT for a rapidly expanding list of Gross MM, Brophy SL, Metzger LM, Shoushtari CS,
conditions, making this method increasingly Splinter R, Reich W (2003) Validation of a brief quan-
titative measure of autistic traits: comparison of the
accepted by the mainstream medicine, even
social responsiveness scale with the autism diagnostic
though its mechanisms of action remains by far interview-revised. J Autism Dev Disord 33
an area of limited knowledge. (4):427–433
Cotler HB, Hamblin MR, Carroll J (2015) The use of low
level laser therapy (LLLT) for musculoskeletal pain.
Acknowledgments The project was funded by Erchonia
MOJ Orthop Rheumatol 2(5):pii:00068
Corp., Melbourne, FL, the producers of the equipment
Effects of Low-Level Laser Therapy in Autism Spectrum Disorder 129
Dawood MS, Salman SD (2013) Low level diode laser Leucht S, Engel RR (2006) The relative sensitivity of the
accelerates wound healing. Lasers Med Sci clinical global impressions scale and the brief psychi-
28:941–945 atric rating scale in antipsychotic drug trials.
deFreitas LF, Hamblin MR (2016) Proposed mechanisms Neuropsychopharmacology 31:406–412
of photobiomodulation or low–level light therapy. Machado C, Estévez M, Leisman G, Melillo R,
IEEE J Select Topics Quant Electronics 22:7000417 Rodríguez R, DeFina P, Hernández A, Pérez–
deTaboada L, Ilic S, Leichliter-Martha S, Oron U, Oron A, Nellar J, Naranjo R, Chinchilla M, Garófalo N,
Streeter J (2006) Transcranial application of Vargas J, Beltrán C (2015) qEEG spectral and coher-
low-energy laser irradiation improves neurological ence assessment of autistic children in three different
deficits in rats following acute stroke. Lasers Surg experimental conditions. J Autism Dev Disord
Med 38:70–73 45:406–424
Ding S, Meng Han Y, Xue Y (2017) A review on feature McPartland J, Volkmar FR (2012) Autism and related
binding theory and its functions observed in perceptual disorders. Handb Clin Neurol 106:407–418
process. Cognit Comput 9:194–206 Melillo R, Leisman G (2009) Autistic spectrum disorders
Erlicher A, Betz T, Stuhtmann B, Koch D, Milner V, as functional disconnection syndrome. Rev Neurosci
Raizen J (2002) Guiding neuronal growth with light. 20:111–131
Proc Nat Acad Sci U S A 99:16024–16028 Melillo R, Leisman G (2010) Neurobehavioral disorders
Farivar S, Malekshahabi T, Shiari R (2014) Biological of childhood: an evolutionary perspective. Springer,
effects of low level laser therapy. J Lasers Med Sci New York
5:58–62 Naeser MA, Saltmarche A, Krengel MH, Hamblin MR,
Farmer C, Aman MG (2012) Aberrant behavior checklist. Knight JA (2011) Improved cognitive function after
In: Volkmar F (ed) Encyclopedia of autism Spectrum transcranial, light emitting diode treatments in chronic,
disorders. Springer, New York traumatic brain injury: two case reports. Photomed
Fornito A, Zalesky A, Pantelis C, Bullmore ET (2012) Laser Surg 29:351–358
Schizophrenia, neuroimaging and connectomics. Naeser MA, Zafonte R, Krengel MH, Martin PI, Frazier J,
NeuroImage 62:2296–2314 Hamblin MR, Knight JA, Meehan WP 3rd, Baker EH
Friederici AD, Gierhan SM (2013) The language network. (2014) Significant improvements in cognitive perfor-
Curr Opin Neurobiol 23:250–254 mance post–transcranial, red/near–infrared light–emit-
Fuster JM (1997) Network memory. Trends Neurosci ting diode treatments in chronic, mild traumatic brain
20:451–459 injury: open–protocol study. J Neurotrauma
Gadow KD, Drabick DA (2012) Symptoms of autism and 31:1008–1017
schizophrenia spectrum disorders in clinically referred Oron A, Oron U, Chen J, Eilam A, Zhang C, Sadeh M,
youth with oppositional defiant disorder. Res Dev Lampl Y, Streeter J, DeTaboada L, Chopp M (2006)
Disabil 33:1157–1168 Low level laser therapy applied transcranially to rats
Goldman-Rakic PS (1988) Topography of cognition: par- after induction of stroke significantly reduces long-
allel distributed networks in primate association cortex. term neurological deficits. Stroke 37:2620–2624
Ann Rev Neurosci 11:137–156 Oron A, Oron U, Streeter J, deTaboada L,
Greicius MD, Krasnowm B, Reissm AL, Menon V (2003) Alexandrovich A, Shohami E (2007) Low–level laser
Functional connectivity in the resting brain: a network therapy applied transcranially to mice following trau-
analysis of the default mode hypothesis. Proceedings matic brain injury significantly reduces long–term neu-
of the Nat Acad Sci U S A 100:253–258 rological deficits. J Neurotrauma 24:651–656
Griffiths TD, Kumar S (2017) Driving working memory. Pessoa L (2012) Beyond brain regions: network perspec-
Neuron 94:5–6 tive of cognition–emotion interactions. Behav Brain
Grover F Jr, Weston J, Weston M (2017) Acute effects of Sci 35:158–159
near infrared light therapy on brain state in healthy Power JD, Petersen SE (2013) Control–related systems in
subjects as quantified by qEEG measures. Photomed the human brain. Curr Opin Neurobiol 23:223–228
Laser Surg 35:136–141 Purushothuman S, Johnstone DM, Nandasena C, van
Hamblin MR (2018) Photobiomodulation for traumatic Eersel J, Ittner LM, Mitrofanis J, Stone J (2015) Near
brain injury and stroke. J Neurosci Res 96:731–743 infrared light mitigates cerebellar pathology in trans-
Hiwaki O, Miyaguchi H (2018) Noninvasive measurement genic mouse models of dementia. Neurosci Lett
of dynamic brain signals using light penetrating the 591:155–159
brain. PLoS One 13:e0192095 Quirk BJ, Whelan HT (2011) Near–infrared irradiation
Lapchak PA, deTaboada L (2010) Transcranial near infra- photobiomodulation: the need for basic science.
red laser treatment (NILT) increases cortical adeno- Photomed Laser Surg 29:143–144
sine-50 -triphosphate (ATP) content following embolic Raichle ME (2015) The brain’s default mode network.
strokes in rabbits. Brain Res 1306:100–105 Ann Rev Neurosci 38:413–427
Leisman G, Machado C, Melillo R, Mualem R (2012) Raichle ME, MacLeod AM, Snyder AZ, Powers WJ,
Intentionality and “free–will” from a neurodeve- Gusnard DA, Shulman GL (2001) Default mode of
lopmental perspective. Front Integr Neurosci 6:36 brain function. Proc Nat Acad Sci U S A 98:676–682
130 G. Leisman et al.
Rojahn J, Aman MG, Matson JL, Mayville E (2003) The Vissers ME, Cohen MX, Geurts HM (2012) Brain connec-
aberrant behavior checklist and the behavior problems tivity and high functioning autism: a promising path of
inventory: convergent and divergent validity. Res research that needs refined models, methodological
Develop Disabil 24:391–404 convergence, and stronger behavioral links. Neurosci
Salo E, Salmela V, Salmi J, Numminen J, Alho K (2017) Biobehav Rev 36:604–625
Brain activity associated with selective attention, Wang L, Hu L, Grygorczyk R, Shen X, Schwarz W (2015)
divided attention and distraction. Brain Res Modulation of extracellular ATP content of mast cells
1664:25–36 and DRG neurons by irradiation: studies on underlying
Scherman M, Mishina OS, Lombardi P, Giacobino E, mechanism of low–level–laser therapy. Med Inf
Laurat J (2012) Enhancing electromagnetically– 2015:630361
induced transparency in a multilevel broadened Wang X, Tian F, Soni SS, Gonzalez–Lima F, Liu H (2016)
medium. Opt Express 20:4346–4351 Interplay between up–regulation of cytochrome–c–
Schiffer F, Johnston AL, Ravichandran C, Polcari A, oxidase and hemoglobin oxygenation induced by
Teicher MH, Webb RH, Hamblin MR (2009) Psycho- near–infrared laser. Scientif Rep 6:30540
logical benefits 2 and 4 weeks after a single treatment Weis S, Rivière R, Deléglise S, Gavartin E, Arcizet O,
with near infrared light to the forehead: a pilot study of Schliesser A, Kippenberg TJ (2010) Optomechanically
10 patients with major depression and anxiety. Behav induced transparency. Science 330:1520–1523
Brain Funct 5:46 Wu Q, Xuan W, Ando T, Xu T, Huang L, Huang YY,
Seeley WW, Crawford RK, Zhou J, Miller BL, Greicius Dai T, Dhital S, Sharma SK, Whalen MJ, Hamblin MR
MD (2009) Neurodegenerative diseases target large– (2012a) Low–level laser therapy for closed–head trau-
scale human brain networks. Neuron 62:42–52 matic brain injury in mice: effect of different
Shanks S, Leisman G (2018) Perspective on broad acting wavelengths. Lasers Surg Med 44:218–226
clinical physiological effects of photobiomodulation. Wu X, Alberico SL, Moges H, deTaboada L, Tedford CE,
Adv Exp Med Biol. https://doi.org/10.1007/5584_ Anders JJ (2012b) Pulsed light irradiation improves
2018_188 behavioral outcome in a rat model of chronic mild
South M, Newton T, Chamberlain PD (2012) Delayed stress. Lasers Surg Med 44:227–232
reversal learning and association with repetitive behav- Zappella M (2012) Reversible autism and intellectual dis-
ior in autism spectrum disorders. Autism Res 5 ability in children. Am J Med Genet C Semin Med
(6):398–406 Genet 160C(2):111–117
Testa–Silva G, Loebel A, Giugliano M, de Kock CP, Zivin JA, Albers GW, Bornstein N, Chippendale T,
Mansvelder HD, Meredith RM (2012) Hypercon- Dahlof B, Devlin T, Fisher M, Hacke W, Holt W,
nectivity and slow synapses during early development Ilic S, Kasner S, Lew R, Nash M, Perez J, Rymer M,
of medial prefrontal cortex in a mouse model for men- Schellinger P, Schneider D, Schwab S, Veltkamp R,
tal retardation and autism. Cereb Cortex 22:1333–1342 Walker M, Streeter J, NeuroThera Effectiveness and
van den Heuvel MP, Sporns O (2013) Network hubs in the Safety Trial-2 Investigators (2009) Effectiveness and
human brain. Trends Cog Sci 17:683–696 safety of transcranial laser therapy for acute ischemic
stroke. Stroke 40:1359–1364
Adv Exp Med Biol - Clinical and Experimental Biomedicine (2018) 3: 131–138
https://doi.org/10.1007/5584_2018_239
# Springer International Publishing AG, part of Springer Nature 2018
Published online: 4 July 2018
Epidemiology of Granulomatosis
with Polyangiitis in Poland, 2011–2015
131
132 K. Kanecki et al.
GPA can occur at all ages, and it has been respectively. The incidence was stable throughout
reported from infancy to old age. the study period. There was an increase in the
Results from the population-based surveys annual prevalence from 28.8 per million in 1990
support the existence of differential international to 64.8 per million in 2005 (Watts et al. 2009).
variations in the occurrence of GPA. It has Another British study has reported 462 GPA
become a well-accepted notion that the geo- cases diagnosed in 1997–2013 (Pearce et al.
graphic distribution of GPA decreases with the 2017). In that study, the overall estimate of GPA
decreasing latitude, both north and south of the incidence was 11.8 (95% CI, 10.7–12.9) per mil-
equator. In the Nordic countries, the incidence or lion person-years. The incidence was 0.88 (95%
prevalence of GPA is roughly two- to sevenfold CI, 0.40–1.96) in children (aged <16 years) and
higher than that in the southern Europe. A similar 14.0 (95% CI, 12.8–15.4) in adults. The preva-
discovery was made in the southern hemisphere, lence in 2013 was 134.9 (95% CI, 121.3–149.6)
with a threefold higher incidence of GPA, con- per million population.
trolled for age, sex, and ethnicity, in the southern- In a study from northern Germany in
most part of New Zealand compared with the 1998–2005, GPA incidence ranged between
northernmost part (O’Donnell et al. 2007). Ethnic 6 and 12 new cases per million per year (median
origin is also thought to play a part in the devel- 8.5) (Herlyn et al. 2008). In another German
opment of GPA. In a multiethnic population of study in 1994–2006, the prevalence of GPA rose
France, GPA is twice as prevalent among popula- from 58 to 98 per million population (Herlyn
tion of the European versus non-European et al. 2014). In a Swedish study conducted until
descent; and non-European cases of GPA are 30 June 2008, using hospital records and serology
rather rare (Mahr et al. 2004). Relatives of database, the annual incidence for GPA was
patients with GPA are at increased risk of being estimated at 9.8 (95% CI, 7.4–12.2) per million
diagnosed with other autoimmune and inflamma- (Mohammad et al. 2009). In a retrospective Nor-
tory diseases, indicating a shared susceptibility to wegian cohort study, using hospital discharge
those diseases (Knight et al. 2010a). records, the annual incidence increased from 5.2
HLA-DPB1*0401 is a GPA susceptibility allele, (95% CI, 2.7–9.0) during 1984–1988 to 12.0
and the population frequency of the allele may (95% CI, 8.0–17.3) during 1994–1998 per mil-
help explain variations in the GPA incidence lion. The point prevalence increased from 30.4
described in the literature (Watts et al. 2015). (95% CI, 16.6–51.0) to 95.1 (95% CI,
69.1–129.0) per million population (Koldingsnes
and Nossent 2000). A Finish study has shown an
1.1 Incidence and Prevalence annual increase in the incidence from 1.9 (95%
of Granulomatosis CI, 1.4–2.6) during 1981–1985 to 9.3 (95% CI,
with Polyangiitis (GPA) 8.1–10.6) during 1996–2000 per million (Takala
et al. 2008). In the region of Lugo of northwest
GPA is a rare disease. Therefore, collection of a Spain, a study based on the data from 1988–2001,
large population of patients to determine the inci- coming from a single reference hospital and
dence and prevalence takes time, which makes regarding all patients aged 15 or older, has
studies rather scarce. In a study conducted in reported the annual GPA incidence of 2.95
1990–2005, using of the General Practice (95% CI, 1.44–6.05) per million population
Research Database in the UK, 295 patients with (Gonzalez-Gay et al. 2003). In the Costa del Sol
newly diagnosed GPA have been identified. region of southern Spain in 2010, the annual
The overall annual incidence of GPA was 8.4 incidence and prevalence were 2.1 and 15.8 per
(95% confidence interval (95% CI), 7.5–9.4) million, respectively (Romero-Gómez et al.
per million. In that study, the annual rate was 2015). In northern Italy, the mean annual GPA
8.1 (95% CI, 6.8–9.6) per million and 8.8 (95% incidence increased from 1.7 during 1995–1999
CI, 7.4–10.3) per million in women and men, to 3.4 per million per year during 2005–2009.
Epidemiology of Granulomatosis with Polyangiitis in Poland, 2011–2015 133
The point prevalence increased from 17.8 (95% was determined on the basis of the International
CI, 7.7–35.1) in 1999 to 34.3 (95% CI, Classification of Diseases (M31.3) for the dis-
20.3–54.2) per million in 2009 (Catanoso et al. charge diagnosis using the Integrated South
2014). In a French urban multiethnic population, Australian Activity Collection (Hissaria et al.
GPA prevalence was estimated at 23.7 (95% CI, 2008). A different study also shows that discharge
16–31) per million adults (Mahr et al. 2004). In records might be a useful element in the epidemi-
Lithuania, the annual GPA incidence is reported ological evaluation of eosinophilic
at 2.1 (95% CI, 1.1–4.1) (Dadoniene et al. 2005), granulomatosis with polyangiitis, a necrotizing
whereas in Turkey at 4.8 per million population vasculitis (Kanecki et al. 2017).
(Pamuk et al. 2016). Between 2000 and 2004,
GPA prevalence in the Australian Capital Terri-
tory and the surrounding rural region was 2 Methods
148.0 (95% CI, 125.1–173.9) per million, and
the disease-specific incidence was 8.4 per million This study was approved by the institutional
per year (Ormerod and Cook 2008). A 5-year inci- board for clinical studies. It is a retrospective,
dence of GPA in South Australia was 56.0 (95% population-based study based on a review of hos-
CI, 44.1–68.4) per million over the period pital discharge files reporting the diagnosis of
2001–2005 (Hissaria et al. 2008). A 5-year preva- GPA in Poland in 2011–2015. The incidence of
lence of GPA varies in New Zealand between GPA was estimated on the basis of data obtained
29–75 cases per million population (de Zoysa from the Polish hospital morbidity study carried
2013). Using the National Inpatient Sample, the out by the National Institute of Public Health.
largest all-payer inpatient database in the USA, Data on all inpatients, discharged and those who
the annual hospitalization rate for patients with passed away, excluding psychiatric and military
the principal diagnosis of GPA increased by 24%, hospitals, were sent to the National Institute of
from 5.1 to 6.3 per 1 million US population, in the Public Health, on a monthly basis. The
period 1993–2011 (p < 0.0001 for trend) (Wallace anonymized data included the information on
et al. 2017). The annual incidence of GPA in hospitalization with ICD10-coded diagnoses,
Japan is 2.1 per million adult population (95% date of admittance and discharge, birth date, gen-
CI, 0.6–3.7 per million) (Fujimoto et al. 2011). der, and the place of residence. We evaluated the
Another study performed in a southern hemisphere patients with the first-time hospitalization for
region shows a 5-year GPA prevalence of GPA. The incidence rate was calculated using
152 (95% CI, 117–186) per million population the number of GPA patients and the
(Gibson et al. 2006). corresponding census data. In addition, demo-
In Poland, the average annual GPA incidence graphic data for the general Polish population
has been estimated at 4.9 per million in the gen- were obtained from the Central Statistical Office
eral population, 5.8 per million in the adult popu- in Poland (CSO 2016).
lation, and 1 per million in the population below
18 years of age (Kanecki et al. 2014). In view of
the paucity of epidemiological data, the present 2.1 Statistical Analysis
study seeks to define a current burden of GPA in
Poland. The updated information would also help Categorical data were expressed as percentages.
compare recent trends in GPA morbidity in Assuming the Poisson distribution of the cases
Poland with those in other geographical regions. analyzed, 95% confidence intervals (CI) were cal-
There is a growing interest in the development culated. Nonparametric tests were used when the
of surveillance ways for GPA, using the existing normality assumptions were unmet. A p-value of
resources such as hospital morbidity databases. In <0.05 defined statistically significant differences.
a New Zealand study, a 5-year incidence of GPA The incidence rate was calculated using the num-
in South Australia over the period 2001–2005 ber of disease cases divided by the population
134 K. Kanecki et al.
figure. Linear regression was used for trend anal- The age distribution of patients at their first
ysis. All analyses were performed using Statistica hospitalization with GP diagnosis is presented in
software v2013 and WINPEPI software Fig. 3. The median and mean age of all patients
(Abramson 2011). were 55 and 52 years, respectively. The average
annual incidence of GPA in Poland was
estimated at 7.7 per million population (95%
CI, 4.1–11.4), and the point prevalence was
3 Results
38.4 per million at the end of 2015. GPA inci-
dence decreased significantly over time in
There were 6995 cases of hospitalization with
Poland between 2011 and 2015 (12.9 vs. 6.7
GPA diagnosis in 2011–2015, out of which 1491
per million, respectively; p < 0.001). Of note,
cases represented patients with the first-time GPA
there was a significant seasonal decreasing trend
diagnosis. The majority consisted of adult patients,
in the first-time GPA incidence from January to
1455 (97.5%; 18 years of age). The whole
December, presented in Fig. 4 (p < 0.01). There
patient group comprised 742 (49.8%) males and
were no significant differences in age according
749 (50.2%) females. The number of GPA patients
to gender at the time of GPA diagnosis. In all
by gender per year is presented in Fig. 1. The
6,995 hospitalized cases of GPA, there were
number of GPA patients in rural and urban regions
108 deaths during the study period (54 females
per year is presented in Fig. 2. There were signifi-
and 54 males; mean age 60 years, range
cant differences in the incidence rate between peo-
19–86 years). Fifty-one out of the 108 deaths
ple in urban (995 patients) and rural areas
were observed in patients hospitalized with
(496 patients); 67% vs. 33%, p < 0.001.
newly diagnosed GPA (27 females and
250 Male
polyangiitis (GPA) by
gender in 2011–2015 200
150
100
50
0
2011 2012 2013 2014 2015
Calendar Year
Fig. 3 Age distribution of patients diagnosed with granulomatosis with polyangiitis (GPA) in 2011–2015
24 males). During the study period, A 1-year Mohammad et al. 2009; Watts et al. 2009; Herlyn
survival rate for GPA was 94%. et al. 2008; Takala et al. 2008; Dadoniene et al.
2005; Gonzalez-Gay et al. 2003; Koldingsnes and
Nossent 2000). Moreover, the average annual
4 Discussion incidence in 2011–2015 was higher than that
reported in a previous epidemiological study
On the basis of the inpatient discharge files, we from Poland (Kanecki et al. 2014). The GPA
estimated that the average annual incidence of incidence was also higher than that previously
GPA in Poland was 7.7 per million population. reported in a recent British study (Pearce et al.
This finding is comparable with data from other 2017). The GPA prevalence in northern Germany
European countries (Pearce et al. 2017; Romero- rose from 58 to 98 per million population in
Gómez et al. 2015; Catanoso et al. 2014; 1994–2006 (Herlyn et al. 2014). However in
136 K. Kanecki et al.
another study, it amounted to 6–12 new cases per distributed across age groups. Population-based
million per year in the period of 1998–2005 studies indicate that GPA may occur in a wide
(Herlyn et al. 2008). A high GPA incidence in range of age (Lane et al. 2003; Koldingsnes and
2011 in Poland and a significant drop in succes- Nossent 2000). A wide range of age at the first-
sive years may be related to environmental, time GPA hospitalization was also noticed in the
infectious, or other unknown factors. present study, the median and mean age of 55 and
Additionally, we noticed a significantly 52 years, respectively. The wide age ranges may
decreasing trend in GPA incidence month-to- imply that stochastic events are at work in the
month, with the highest rate of newly diagnosed disease pathogenesis (Gibelin et al. 2011). In the
GPA in January. It is a long-held belief that the UK, median age of patients with newly diagnosed
etiology of systemic vasculitides could be infec- GPA was 59 years (interquartile range
tious. For vasculitides with a frequent pulmonary 47–70 years), with a mean of 52 years (Watts
involvement, this hypothesis is all the more rele- et al. 2009). In other studies, mean age at the
vant because lung diseases are notably initiated diagnosis was 64 years (Herlyn et al. 2014) and
by airborne infections. The most compelling 55–65 years (Herlyn et al. 2008). In Finland,
observation to support the infectious cause mean age at the diagnosis was 53 years (Takala
comes from the circumstantial evidence showing et al. 2008), and in Norway it was 50 (range
that the onset of these diseases is preceded by 10–84 years). The highest incidence rate occurred
infectious events and from the studies showing in Norwegian men aged 65–74 (Koldingsnes and
that GPA first becomes evident during the fall and Nossent 2000).
winter months, the seasons with frequent In the present study, male to female ratio was
respiratory infection epidemics. However, there 0.99. Likewise, 51–52% of GPA patients in the
is a controversy regarding the concept of season- British studies were male (Wallace et al. 2016;
ality concerning the onset of disease, with other Watts et al. 2009). Yet another British study has
studies suggesting either no GPA seasonality also reported that GPA incidence is lower in
(Mahr et al. 2006a) or its increased occurrence females (adjusted incidence rate ratio ¼ 0.68;
in the summer (Mahr et al. 2006b). 95% CI, 0.56–0.81), with the greatest incidence
We found that there were significantly more in those aged 55–69 years (adjusted incidence
patients who lived in urban rather than rural rate ratio ¼ 6.8, 95% CI, 4.9–9.6; reference
areas. Likewise, a greater GPA prevalence was group 16–39 years of age) (Pearce et al. 2017).
observed in urban areas in northern Germany In Norway, 62% of GPA patients are male
(Herlyn et al. 2014). An Australian study has (Koldingsnes and Nossent 2000), whereas only
also reported a trend toward a higher incidence 49% of male patients have been reported in
of GPA in rural areas (Ormerod and Cook 2008). Finland (Takala et al. 2008).
In a multicultural population of New Zealand, During the present study period, 108 deaths in
GPA predominates in Europeans (de Zoysa all GPA patients were observed, with a 1-year
2013). A long-standing exposure to silica may survival rate for GPA was 94%. All the deaths
be associated with GPA onset (Hogan et al. were observed in adults. Among the dead, there
2007). Other positive relationships were were 51 patients with newly diagnosed GPA. In a
suggested for GPA and farming, particularly Swedish study, 5-year survival rate for GPA was
due to a contact with livestock and organic 83% (Mohammad et al. 2009). In the UK, mortal-
solvents (Lane et al. 2003). The role of occupa- ity was reported at 13.6% for the 1-year rate
tional exposures as determinants of GPA risk is (Pearce et al. 2017). In New Zealand, 1- and
at variance with a recent study that demonstrates 10-year probability of survival in GPA patients
no such relationship between (Knight et al. were 91% and 62%, respectively (Khan et al.
2010b). 2012). A meta-analysis of GPA cases indicates a
GPA occurs primarily during adulthood, but similar mortality risk, with a standardized mortal-
its age-specific incidence rates are not uniformly ity ratio of 2.63 (95% CI 2.02–3.43) (Tan et al.
Epidemiology of Granulomatosis with Polyangiitis in Poland, 2011–2015 137
2017). A population-based British study suggests prevalence at 36 per million at the end of 2015.
that the survival rate of GPA patients has The incidence rate of GPA in Poland is similar to
improved considerably over the past two decades, the data reported by other European countries. A
affirming the benefits of recent trends in the man- higher incidence of GPA was reported in urban
agement of GPA and its complications (Wallace than rural regions of Poland. During the study
et al. 2016). Using the National Inpatient Sample, period, a 1-year survival rate for GPA was 94%.
the largest all-payer inpatient database in the Trends in GPA incidence observed in this study
USA, it has been reported that inhospital deaths may be related to the environmental, infectious,
in patients with the principal diagnosis of GPA or other unknown factors.
declined by 73%, from 9.1% to 2.5% (p < 0.0001
for trend) from 1993 to 2011 (Wallace et al. Conflicts of Interest Authors declare no conflict of inter-
2017). Patients with GPA have a ninefold est in relation to this article.
increased risk of death in the first year of the
disease, attributed to infection, active vasculitis,
and renal failure. The risk is at its lowest between References
the second and the eighth year of the disease,
although higher than that in the control popula- Abramson JH (2011) WINPEPI updated: computer
tion. There is an increased mortality from the programs for epidemiologists, and their teaching
ninth year onward which remains unexplained potential. Epidemiol Perspect Innov 8(1):1
Catanoso M, Macchioni P, Boiardi L, Manenti L,
(Luqmani et al. 2011). Tumiati B, Cavazza A, Luberto F, Pipitone N,
This study has several limitations. The data- Salvarani C (2014) Epidemiology of granulomatosis
base which was used provides general informa- with polyangiitis (Wegener’s granulomatosis) in north-
tion regarding hospitalization of GPA patients, ern Italy: a 15-year population-based study. Semin
Arthritis Rheum 44:202–207
thereby allowing the authors to present an epide- CSO (2016) Central Statistical Office of Poland. http://stat.
miological report. However, it did not allow gov.pl. Accessed on 13 Oct 2016
examining risk factors or other correlates of the Dadoniene J, Kirdaite G, Mackiewicz Z, Rimkevicius A,
disease in individual patients. Another limitation Haugeberg G (2005) Incidence of primary systemic
vasculitides in Vilnius: a university hospital population
is that the evaluation relied on discharge files based study. Ann Rheum Dis 64:335–336
from inpatient hospitalizations, which excludes de Zoysa JR (2013) Vasculitis in New Zealand. Clin Exp
patients who were seen at the outpatient basis. Nephrol 17:712–714
This inaccuracy may result in the underestimation Fujimoto S, Watts RA, Kobayashi S, Suzuki K, Jayne
DRW, Scott DGI, Hashimoto H, Nunoi H (2011)
of the incidence rate. On the other hand, this type Comparison of the epidemiology of anti-neutrophil
of vasculitis is a systemic necrotizing disease with cytoplasmic antibody-associated vasculitis between
multi-organ involvement and usually requires Japan and the U.K. Rheumatol Oxf Engl
treatment or advanced diagnostic procedures to 50:1916–1920
Gibelin A, Maldini C, Mahr A (2011) Epidemiology and
be performed at short-notice in the inpatient etiology of wegener granulomatosis, microscopic
setting. In addition, the first-time GPA diagnosis polyangiitis, churg-strauss syndrome and Goodpasture
in the inpatient discharge file may not necessarily syndrome: vasculitides with frequent lung involve-
be the very time of the disease occurrence. A ment. Semin Respir Crit Care Med 32:264–273
Gibson A, Stamp LK, Chapman PT, O’Donnell JL (2006)
longer duration of the observation may improve The epidemiology of Wegener’s granulomatosis and
the estimation. microscopic polyangiitis in a Southern Hemisphere
region. Rheumatol Oxf Engl 45:624–628
Gonzalez-Gay MA, Garcia-Porrua C, Guerrero J,
Rodriguez-Ledo P, Llorca J (2003) The epidemiology
5 Conclusions of the primary systemic vasculitides in Northwest
Spain: implications of the Chapel Hill consensus con-
The average annual incidence of GPA in Poland ference definitions. Arthritis Rheum 49:388–393
was estimated at 7.7 per million population (95% Herlyn K, Hellmich B, Gross WL, Reinhold-Keller E
(2008) Stable incidence of systemic vasculitides in
confidence interval, 4.1–11.4) and the point
138 K. Kanecki et al.
Schleswig-Holstein, Germany. Dtsch Arztebl Int Mahr A, Artigues N, Coste J, Aouba A, Pagnoux C,
105:355–361 Guillevin L, French Vasculitis (2006b) Study group sea-
Herlyn K, Buckert F, Gross WL, Reinhold-Keller E (2014) sonal variations in onset of Wegener’s granulomatosis:
Doubled prevalence rates of ANCA-associated increased in summer? J Rheumatol 33:1615–1622
vasculitides and giant cell arteritis between 1994 and Mohammad AJ, Jacobsson LTH, Westman KWA,
2006 in northern Germany. Rheumatology Sturfelt G, Segelmark M (2009) Incidence and survival
53:882–889 rates in Wegener’s granulomatosis, microscopic
Hissaria P, Cai FZJ, Ahern M, Smith M, Gillis D, Roberts- polyangiitis, Churg-Strauss syndrome and polyarteritis
Thomson P (2008) Wegener’s granulomatosis: epide- nodosa. Rheumatol Oxf 48:1560–1565
miological and clinical features in a south Australian O’Donnell JL, Stevanovic VR, Frampton C, Stamp LK,
study. Intern Med J 38:776–780 Chapman PT (2007) Wegener’s granulomatosis in
Hogan SL, Cooper GS, Savitz DA, Nylander-French LA, New Zealand: evidence for a latitude-dependent inci-
Parks CG, Chin H, Jennette CE, Lionaki S, Jennette dence gradient. Intern Med J 37:242–246
JC, Falk RJ (2007) Association of silica exposure with Ormerod AS, Cook MC (2008) Epidemiology of primary
anti-neutrophil cytoplasmic autoantibody small-vessel systemic vasculitis in the Australian capital territory
vasculitis: a population–based, case-control study. Clin and south-eastern New South Wales. Intern Med J
J Am Soc Nephrol 2:290–299 38:816–823
Kanecki K, Życińska K, Moskalewicz B, Tyszko P (2014) Pamuk ÖN, Dönmez S, Calay{r GB, Pamuk GE (2016)
Granulomatosis with polyangiitis in Poland – epidemi- The epidemiology of antineutrophil cytoplasmic
ological study. Reumatologia 52:99–104 antibody-associated vasculitis in northwestern Turkey.
Kanecki K, Nitsch-Osuch A, Gorynski P, Tarka P, Tyszko Clin Rheumatol 35:2063–2071
P (2017) Hospital morbidity database for epidemiolog- Pearce FA, Grainge MJ, Lanyon PC, Watts RA, Hubbard
ical studies on Churg-Strauss syndrome. Adv Exp Med RB (2017) The incidence, prevalence and mortality of
Biol 980:19–25 granulomatosis with polyangiitis in the UK Clinical
Khan AR, Chapman PT, Stamp LK, Wells JE, O’Donnell Practice Research Datalink. Rheumatol Oxf
JL (2012) Wegener’s granulomatosis: treatment and 56:589–596
survival characteristics in a high-prevalence southern Romero-Gómez C, Aguilar-García JA, García-de-Lucas
hemisphere region. Intern Med J 42:23–26 MD, Cotos-Canca R, Olalla-Sierra J, García-Alegría
Knight A, Sandin S, Askling J (2010a) Increased risk of JJ, Hernández-Rodríguez J (2015) Epidemiological
autoimmune disease in families with Wegener’s study of primary systemic vasculitides among adults
granulomatosis. J Rheumatol 37:2553–2558 in southern Spain and review of the main epidemiolog-
Knight A, Sandin S, Askling J (2010b) Occupational risk ical studies. Clin Exp Rheumatol 33(2 Suppl
factors for Wegener’s granulomatosis: a case-control 89):11–18
study. Ann Rheum Dis 69:737–740 Takala JH, Kautiainen H, Malmberg H, Leirisalo-Repo M
Koldingsnes W, Nossent H (2000) Epidemiology of (2008) Incidence of Wegener’s granulomatosis in
Wegener’s granulomatosis in northern Norway. Arthri- Finland 1981–2000. Clin Exp Rheumatol 26(3 Suppl
tis Rheum 43:2481–2487 49):81–85
Lane SE, Watts RA, Bentham G, Innes NJ, Scott DGI Tan JA, Dehghan N, Chen W, Xie H, Esdaile JM, Avina-
(2003) Are environmental factors important in primary Zubieta JA (2017) Mortality in ANCA-associated vas-
systemic vasculitis? A case-control study. Arthritis culitis: a meta-analysis of observational studies. Ann
Rheum 48:814–823 Rheum Dis 76:1566–1574
Langford C (2012) Clinical features and diagnosis of small- Wallace ZS, Lu N, Unizony S, Stone JH, Choi HK (2016)
vessel vasculitis. Cleve Clin J Med 79(Suppl 3):3–7 Improved survival in granulomatosis with polyangiitis:
Luqmani R, Suppiah R, Edwards CJ, Phillip R, Maskell J, a general population-based study. Semin Arthritis
Culliford D, Jayne D, Morishita K, Arden N (2011) Rheum 45:483–489
Mortality in Wegener’s granulomatosis: a bimodal pat- Wallace ZS, Lu N, Miloslavsky E, Unizony S, Stone JH,
tern. Rheumatol Oxf Engl 50:697–702 Choi HK (2017) Nationwide trends in hospitalizations
Mahr A, Guillevin L, Poissonnet M, Aymé S (2004) and in-hospital mortality in granulomatosis with
Prevalences of polyarteritis nodosa, microscopic polyangiitis (Wegener’s). Arthritis Care Res
polyangiitis, Wegener’s granulomatosis, and Churg- 69:915–921
Strauss syndrome in a French urban multiethnic popu- Watts RA, Al-Taiar A, Scott DGI, Macgregor AJ (2009)
lation in 2000: a capture-recapture estimate. Arthritis Prevalence and incidence of Wegener’s granulomatosis
Rheum 51:92–99 in the UK general practice research database. Arthritis
Mahr AD, Neogi T, Merkel PA (2006a) Epidemiology of Rheum 61:1412–1416
Wegener’s granulomatosis: lessons from descriptive Watts RA, MacGregor AJ, Mackie SL (2015) HLA allele
studies and analyses of genetic and environmental variation as a potential explanation for the geographi-
risk determinants. Clin Exp Rheumatol 24(2 Suppl cal distribution of granulomatosis with polyangiitis.
41):82–91 Rheumatol Oxf Engl 54:359–362