Original Research ARTICLE
Physical and Chemical Processes and the Morphofunctional Characteristics of Human Erythrocytes in Hyperglycaemia
- 1Department of Biotechnology, Bioengineering and Biochemistry, Faculty of Biotechnology and Biology, Federal State-Financed Academic Institution of Higher Education, National Research Ogarev Mordovia State University, Saransk, Russia
- 2Department of Biochemistry, Biotechnology and Bioengineering, Faculty of Biology, Federal State-Funded Educational Institution of Higher Professional Education, Samara State University, Samara, Russia
- 3Division of Cancer Research, Ninewells Hospital and Medical School, University of Dundee, Dundee, United Kingdom
- 4CMCBR, Abertay University, Dundee, United Kingdom
Background: This study examines the effect of graduated hyperglycaemia on the state and oxygen-binding ability of hemoglobin, the correlation of phospholipid fractions and their metabolites in the membrane, the activity of proteolytic enzymes and the morphofunctional state of erythrocytes.
Methods: Conformational changes in the molecule of hemoglobin were determined by Raman spectroscopy. The structure of the erythrocytes was analyzed using laser interference microscopy (LIM). To determine the activity of NADN-methemoglobinreductase, we used the P.G. Board method. The degree of glycosylation of the erythrocyte membranes was determined using a method previously described by Felkoren et al. Lipid extraction was performed using the Bligh and Dyer method. Detection of the phospholipids was performed using V. E. Vaskovsky method.
Results: Conditions of hyperglycaemia are characterized by a low affinity of hemoglobin to oxygen, which is manifested as a parallel decrease in the content of hemoglobin oxyform and the growth of deoxyform, methemoglobin and membrane-bound hemoglobin. The degree of glycosylation of membrane proteins and hemoglobin is high. For example, in the case of hyperglycaemia, erythrocytic membranes reduce the content of all phospholipid fractions with a simultaneous increase in lysoforms, free fatty acids and the diacylglycerol (DAG). Step wise hyperglycaemia in incubation medium and human erythrocytes results in an increased content of peptide components and general trypsin-like activity in the cytosol, with a simultaneous decreased activity of μ-calpain and caspase 3.
Conclusions: Metabolic disorders and damage of cell membranes during hyperglycaemia cause an increase in the population of echinocytes and spherocytes. The resulting disorders are accompanied with a high probability of intravascular haemolysis.
Various physiological, pathological and nutritional conditions such as physical activity, large amounts of sweet food, emotional stress, metabolic syndrome, and diabetes are accompanied by high level of glucose in blood plasma. The high content of glucose in plasma accelerates the probability of non-enzymatic glycosylation of proteins, which induce damage to the cell membrane due to nonspecific aggregation of protein molecules and changes in protein-protein and protein-lipid interactions (Vasilyeva, 2005).
Taken together, these changes initiate the rapid aging of cells and the human organism. Metabolic syndrome significantly accelerates the development of atherosclerotic vascular damage and provokes earlier disability and death. During metabolic syndrome, which is currently the most common pathology of metabolic disorders, glycosylation of erythrocytic membrane proteins induces the impairment of rheological parameters of blood, low deformability and mobility of erythrocytes, high aggregation of erythrocytes and thrombocytes, high blood viscosity, and arterial hypertension (Shilov et al., 2008).
In addition, glycosylation of erythrocytic membrane proteins and hemoglobin during hyperglycaemia increases adhesion to endothelial cells, resulting in membrane destabilization (change in the asymmetry of membrane phospholipids), changes in viscoelastic properties of cells and their morphology (Riquelme et al., 2005). Taken together, these changes can impair the oxygen-transport function of erythrocytes and reduce erythrocyte lifespan. Furthermore, the number of damaged circulating cells and, aging erythrocytes will increase (Lang et al., 2006; Mindukshev et al., 2010).
The biochemical mechanisms of impaired growth in human erythrocytes during the development of hyperglycaemia have not been sufficiently investigated. In particular, there are scarce data on the composition and status of the lipid phase of the membranes, the relationship of these processes with the activity of methemoglobin formation and the activity of apoptotic enzymes. Moreover, there is a lack of data in the literature on the effect of these processes on the morphofunctional state of erythrocytes and their oxygen-transport properties.
Therefore, we aimed to perform a comprehensive study of the effects of graduated hyperglycaemia on the composition of phospholipids, the activity of proteolytic enzymes, and, the consequent effect of ongoing processes on the morphofunctional state of erythrocytes and the oxygen-transport properties of hemoglobin.
Materials and Methods
This study was performed using pure erythrocyte fractions isolated from freshly obtained, anticoagulated donor blood from the regional station of blood transfusion. The average age of the donors was 28.4 ± 1.7 years. The research was approved by the Local Ethics Board at Mordovia State University in accordance with the principles of Good Clinical Practice (protocol number 12 of 17 September 2014). Informed consent statements were signed by all donors participating in the experiment.
Erythrocytes were washed 4 times with PBS until pure fractions were formed and centrifuged at 600 g for 10 min at temperature +4°C. Incubation of the erythrocytes was performed in a Ringer's solution at a ratio of 2:1 at 37°C for 30 min with gentle mixing every 5 min (to avoid haemolysis). For incubation, a Ringer's solution was used containing as follow (in mM/l): NaCl–136.75; KCl–to 2.68; NaHCO3–11.9 per; Na2HPO4–0.342; MgCl2–0.105; CaCl2–1.8, glucose–5 (normoglycaemia).
Hyperglycaemic conditions were characterized by glucose content, in which the graduated concentration was accompanied by an equimolar decrease in the content of sodium chloride: 10 mM/l glucose–NaCl: 131.75 mM/l; 15 mM/l glucose–NaCl: 126.75 mm/l; 20 mM/l NaCl; 121.75 mM/L. The incubation medium was separated by centrifugation in the same manner, and erythrocytes were analyzed using spectrophotometry (Zavodnik and Lapshina, 1996) to determine the correlation of the hemoglobin forms.
Conformation and Properties of Hemoglobin
Conformational changes in the hemoglobin molecule were determined by Raman scattering on In Via Renishaw (UK) microscope with a short-focus high-luminosity monochromator (focus distance–250 mm).
For excitation of the Raman spectra we used a laser with a wavelength of 532 nm, max radiation power of 100 mW, and 100 × objective. The data recorder was a –CCD detector (1,024 × 256 pixels with Peltier cooling to −70°C) with a grating of 1800 lines/mm. The digitized spectra were processed in WIRE 3.3 (part of the unit's software). To analyse the conformation of globin haematoporphyrin, we used specific bands of the RS spectrum, which enabled an estimate of the relative amount of oxyhemoglobin and the ability of hemoglobin to bind to and isolate ligands, including oxygen (Brazhe et al., 2009).
Laser Interference Microscopy
The structure of the erythrocytes was analyzed using laser interference microscopy (LIM) in vitro (Byazhe et al., 2006; Yusipovich et al., 2011; Revin et al., 2016) using MII-4 sysytem (Russia). The measurements were performed at room temperature, and a suspension of erythrocytes in the incubation medium (1:2) was placed on mirror glass. The smear was prepared and covered with cover glass. Images of 10 sites, with a monolayer arrangement of cells in an interference channel, were obtained, with reflected light in each sample. The images were processed using FIJI (Schindelin et al., 2012). The structure of the erythrocytes was assessed by registering the average value of the optical path difference (OPD) and phase image area using at least 100 cells from each sample. The phase volume of the erythrocyte was calculated using the following formula:
Where Fmean is the mean value of the optical path difference, proportional to the thickness of erythrocyte; S is the phase image area of the cells; ncell is the refractive index of erythrocyte, equal to 1.405; nm is the refractive index of the surrounding solution (1.333).
Determination of the Activity of NADN-Methemoglobinreductase
To determine the activity of NADN-methemoglobinreductase, we used the P.G. Board method (Board, 1981). The degree of glycosylation of the erythrocyte membranes was determined using a method previously described by Felkoren et al. (1991). Before defining the proteolytic enzymes, the suspension of erythrocytes was haemolysed by the addition of a buffer of 20 mM three-HCl containing 2 mM EDTA, pH 7.5 in a ratio of 1:9. The haemolysate was maintained at a temperature 2-4°C for 15 min and centrifuged at 16,000 × g for 40 min.
The supernatant was used to determine the activity of μ-calpain with the release of enzyme using ion exchange chromatography (column 3 × 15; DEAE-cellulose) and eluted by a gradient of 0.1–0.4 M NaCl. Next, the mean calpain activity of the fractions was determined using incubation medium previously described by Sorimachi et al. (1997) (imidasole buffer, 4% casein, 50 mM of CaCl2, 50 mM of cysteine) (Stroev et al., 1991; Elce John, 2000; Sorimachi et al., 2000). Then, the general proteolytic activity, in ascending order of the absorption level at a wavelength of 280 nm, was measured following incubation of the haemolysate and protein deposition by 5% trichloroacetic acid (TCA) (Bazarnova et al., 2008).
The release of fractions with calpain activity occurred during movement in a reverse gradient from 0.2 to 0.1 M of sodium chloride immediately after the release of hemoglobin. The activity of μ-calpain was calculated as the difference between activity with and without the inhibitor in the incubation medium (phenylmetylsulfonylfluoride, 2 mM PMSF, Sigma, USA). The content of peptides in the incubation medium and erythrocytes was determined using the Lowry method, with the Bio-Rad DC Protein Assay a set of protein assay Dc reagents (Bio-Rad). The content of peptides in the erythrocytes was determined after deproteinisation by 5% of TCA. The active concentration of caspase-3 in erythrocytes was recorded using an enzyme immunoassay (BD Biosciences, USA) with Stat Fax 3200 microplate reader (USA).
Analysis of Membrane Phospholipids
To analyse the state of the membrane phospholipids, we isolated membranes from the haemolysate using 5 mM NaH2PO4+ 0.5 mM PMSF (phenylmetylsulfonylfluoride) solution, cooled to 0°C, pH 8.0 in a ratio of 1:20. The mixture was incubated for 10 min at 4°C and then centrifuged at 20,000 × g for 40 min (0°C). The supernatant was removed and the residue was resuspended in lysis solution and centrifuged in the same manner. The sample was washed three times. Lipid extraction was performed using the Bligh and Dyer method (Bligh and Dyer, 1959).
To separate the phospholipid fractions we used one-dimensional chromatography combined with chloroform/methanol/glacial acetic acid/water in a ratio of 60/50/1/4 (Evans et al., 1990). Chromatographic separation was performed in a thin layer of silica gel deposited on a glass plate. Standard plates from HPTLC Silicagel 60 F254 (Merck. Germany) were used. To separate DAG and FFA we used the combination of heptane/diethyl ether/glacial acetic acid (60/40/2 by volume). Detection of the phospholipids was performed using V. E. Vaskovsky method (Vaskovsky et al., 1975).
The amount of phospholipid fractions and free fatty acids in the erythrocytic membranes was determined using a densitometer TCL Scanner 3 (Gamag, Switzerland) at an absorption wavelength of 360 nm using a deuterium lamp and winCATC software.
The statistical processing was performed using the Statistika 0.06 software package. To determine the significance levels, the Student's test was used. Repeats in the variation series of different indicators were from 8 to 20 values.
Hyperglycaemic conditions are characterized by a low affinity of hemoglobin to oxygen, which is manifested as a parallel decrease in the content of hemoglobin oxyform and the growth of desoxyform (Table 1).
Table 1. Correlation of hemoglobin formation (%), content of glycosylated hemoglobin (g/l) and the degree of glycosylation of membrane proteins (rel. units).
Upon a 30-min incubation of the cells, an increase in the glycosylation of hemoglobin occurred only when the glucose concentration was 20 mM/l, whereas a small but reliable growth in methemoglobin formation, increase in membrane-bound form of hemoglobin and degree of glycosylation of proteins in erythrocytic membranes has been previously observed when the glucose level exceeded a two-fold level (10 mm/l) (Table 1).
The data obtained using the spectrometric measurements were confirmed by Raman scattering (Table 2).
We found a reliable decrease in the relative content of oxygenated hemoglobin, with an average of 2%. Moreover, we obtained similar spectrophotometric percentages of hemoglobin forms (Tables 1, 2). The conformational changes observed indicate an increased ability of hemoglobin to bind to ligands and the reliable occurrence of low hemoglobin affinity to oxygen in the case of 20-mM/l hyperglycaemia, at approximately 11%. Furthermore, there is a low intensity of symmetric and asymmetric vibrations of pyrrole rings, which indicated a less conformational mobility of the haeme structures and an impairment of the effective binding of ligands (Table 2).
At 15- and 20-mM/l glucose hyperglycaemia, we detected a high level of lysophosphatidylcholine (LFH) in the erythrocyte membrane, and the content of phosphatidylcholine (PFS) was reliably low under in conditions of increasing amounts of glucose, i.e., up to 10 mM/l. In addition, we recorded a low level of sphingomyelin (SM) and fractions of phosphatidylinositol + phosphatidylserine (PI + PS). Strong hyperglycaemia (20 mm/l) is characterized by a low content of all fractions with a simultaneous increase of in LPC (Table 3).
Table 3. Phospholipid composition of erythrocytic membranes under normoglycaemia and graduated hyperglycaemia.
The percentage of increased lysophosphatides under strong 20-mm/l hyperglycaemia was 160%. The SM level was 2 to 3 times lower, with PC–3 times; PI+PS–2.35%; and PEA–4 times. The low content of phospholipids and high content of lysophosphatides indicated the activation of phospholipase A and C, which was accompanied by loosening of the membrane and higher ion permeability (Mills and Needham, 2005).
Analysis of the content of free fatty acids (FFA) and the product of the reaction catalyzed by phospholipase C and diacylglycerol (DAG) showed growth that was proportional to the level of hyperglycaemia, and the percentage of growth with 20 mM/l hyperglycaemia is 133 and 182%, respectively (Table 4).
These results clearly demonstrate the activity of phospholipase.
The normal correlation of phospholipid fractions pre-determines the effective regulation of active and passive transport of substances, the sensitivity of cells to the action of ligands, and the activity of membrane-bound enzymatic systems. The optimal operation of ion-carrying systems (Ca2+ and Na+, K+-ATP-PS) is facilitated by stable PEA content. During severe hyperglycaemia (20 mm/l), a sharp decrease in PEA content will be accompanied by a defect in ion-carrying processes (Table 3). Moreover, it will trigger a general change in the structure of the membrane phospholipid bilayer, as PEA is a structural phospholipid (Delaunay, 2002).
Importantly, the low content of PEA causes a disorder in endoglobular homeostasis and inhibition of the antioxidant activity of cells in (Afanasiev et al., 2007). The low content of phospholipids with polyunsaturated fatty acids (PC, PI, PEA) triggers defatting of membranes and a high correlation of cholesterol/PL, which is accompanied by changes in physico-chemical properties, namely, high microviscosity. These findings indicate that hyperglycaemia is followed by a disruption in membrane permeability.
In hyperglycaemia, human erythrocytes demonstrate a high content of peptide components, which was specifically high when the content of glucose in the incubation medium was 10 mm/l (Table 5). A further increase in the level (degree) of hyperglycaemia was accompanied by an increased release of protein compounds into the incubation medium due to higher cell membrane permeability (Table 5).
Table 5. Content of protein and peptide compounds in the incubation medium and peptides in erythrocytes (during normo- and hyperglycaemia).
Measurement of the proteolytic activity revealed interesting patterns. The release of μ-calpain fractions and identification of the average activity of the enzyme showed a decrease with a parallel increase in the degree of hyperglycaemia, which was statistically significant when the concentration of glucose reached 20 mM/l. A low activity of caspase 3 was detected, which was also strongest under conditions of severe hyperglycaemia (Table 6). Increased proteolytic activity can be detected only by the identification of the overall trypsin-like activity of cytosol with respect to the high content of peptides in the haemolysate, without the isolation of enzymes and addition of third proteolytic substrates (Table 6).
Table 7. Morphological parameters of erythrocytes in normo- and hyperglycaemia, according to LIM results.
During normoglycaemia, erythrocytes appeared similar to discocytes, with a normal distribution of hemoglobin inside the cells. During hyperglycaemia, the estimated parameters were high overall. In addition to the relative areas of the erythrocytes, all indicators were high in conditions of 20-mm/l hyperglycaemia (Table 7).
With a high content of glucose in the incubation medium, echinocytes accumulated in the population of cells, and the cell profile demonstrated visible outgrowths on the erythrocyte surface. At up to 15-mm/l glucose, spherocytes appeared, and the cells appeared more swollen (Figure 3). At 20-mm/l glucose, the cells were slightly compressed compared with 15-mm/l glucose; however, other changes became more robust (Table 7, Figure 4).
Disorder in the metabolic and functional state of human erythrocyte membranes at a high concentration of glucose in the incubation medium occurs primarily due to the general activation of biochemical and physicochemical processes in cells. The shift in the dissociation curve in the direction of deoxygenation causes increased methemoglobin formation as DOHb, which is less resistant to auto-oxidation compared tooxyform (Zavodnik and Lapshina, 1996; Ivanov, 2001). The gradual reduction of the conformational mobility of haeme structures (Table 2) is most likely caused by the high degree of the glycosylated protein component of hemoglobin, its permeability in the membrane, followed by a gradual increase in membrane-bound hemoglobin (Table 1).
The observed increase in the degree of glycosylation of membrane proteins induces changes in the state of membrane phospholipids. The change in protein-lipid interactions results in the activation of membrane phospholipase, causing an increase in LPC and DAG, and results in a high content of FFA in erythrocytic membranes. Membrane permeability under these conditions increases, and the peptide and protein compounds are released into the incubation medium. Some reduction in their release at 20-mm/l hyperglycaemia may be accounted for by the low activity of μ-calpain and caspase 3 due to less affinity to the glycosylated form of hemoglobin. However, their content in the erythrocytes and incubation medium was higher compared to normoglycaemia, most likely due to the continued high trypsin-like activity of the cell cytosol.
Metabolic disorders and injury to the cell membrane trigger morphological changes in erythrocytes. Previous research findings have shown a gradual accumulation of structurally damaged cells in the population. First, in 10-mM/l hyperglycaemia, the number of echinocytes increases (reversible morphological changes) and the phase volume and thickness of the cells increases. Next, stomacytes and spherocytes appear in the population. LIM microscopic results indicated strong cell swelling. In 20-mm/l hyperglycaemia, the number of spherocytes increased. Furthermore, the cells showed some shrinkage, but the ratio of surface/volume values reached a critical level. Thus, haemoglytic destruction is very likely.
The destruction of erythrocytes in the bloodstream in hyperglycaemia will result in a chronic inflammation response, which can damage the body's blood vessels.
Changes in the erythrocytic membrane not only result in morphological disorders of the whole erythrocyte, but they can also consequently result in conformational changes in hemoglobin as well as its oxygen-binding and oxygen-transport function. The low content of oxyhemoglobin, weak symmetric and asymmetric oscillations of pyrrole rings and increase in hemoglobin affinity to ligands, including oxygen, has been demonstrated.
Taken together, we demonstrate the significant risk of the prolonged availability of high glucose concentrations in plasma. Nonspecific glycosylation of membrane proteins and erythrocytesin hemoglobin results in a weak affinity of hemoglobin to oxygen and its loss by cells enroute to tissues in the human body. Moreover, the resulting damage to membranes and cell metabolism increase the probability of an accumulation of functionally defective aging erythrocytes in the circulating population. More rapid spherization of erythrocytes, in the absence of physiological mechanisms of apoptosis, introduces the possibility of necrotic cell death in the bloodstream, resulting in the gradual development of chronic states of hypoxia and inflammation.
VR: designed experiments, interpreted data, wrote manuscript. NK, NG, IG, IS, AT, AP, AS, and ER: designed and performed experiments wrote manuscript. KP: performed experiments and analyzed results. NZ and JCB: made substantial contribution to the analysis and interpretation of the data, revised and critically reviewed the manuscript.
The authors are grateful for the support of the Russian Science Foundation, Grant no. 15-15-10025.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Afanasiev, S. A., Rebrova, T. Y., and Kondratieva, D. S. (2007). Characteristics of phospholipid composition of erythrocyte membranes under conditions of postinfarction cardiosclerosis. Biomed. Chem. 53, 541–546. doi: 10.15690/vramn612
Brazhe, N. A., Abdali, S., Brazhe, A. R., Luneva, O. G., Bryzgalova, N. Y., Parchina, E. Y., et al. (2009). New insight into erythrocyte through in vivo surface-enhanced raman spectroscopy. Biophys. J. 97, 3206–3214. doi: 10.1016/j.bpj.2009.09.029
Byazhe, N. A., Brazhe, A. R., Pavlov, A. N., Erokhva, L. A., Yusipovich, A. I., and Maksimov, G. V. (2006). Unraveling cell processes: interference imaging interwoven with data analysis. J. Biol.Phys. 32, 191–208. doi: 10.1007/s10867-006-9012-1
Evans, W. S., Morro, D. D., O'Braytman, E., Sandelius, A. S., Beylez, E. M., Richardson, P. G., et al. (1990). The Biological Membrane. Methods. Biological Membranes. A Practical Approach. (Moscow: Mir), 424.
Mills, K., and Needham, D. (2005). Lysolipid incorporation in dipalmitoylphosphatidylcholine bilayer membranes enhances the ion permeability and drug release rates at the membrane phase transition. Biochim. Biophys. Acta 1716, 77–96. doi: 10.1016/j.bbamem.2005.08.007
Mindukshev, I. V., Dobrylko, I. A., Krivchenko, A. I., Krivoshlyk, V. V., Goncharov, N. V., Vivulanets, E. V., et al. (2010). Abnormalities of elastic and transporting properties of red blood cells under development of apoptosis. Biochemistry (Moscow) 4, C.22–C.31. doi: 10.1134/S1990747810010046
Revin, V. V., Gromova, N. V., Revina, E. S., Martynova, M. I., Seikina, A. I., Revina, N. V., et al. (2016). Role of membrane lipids in the regulation of erythrocytic oxygen-transport function in cardiovascular diseases. Biomed. Res. Int. 2016, 11. doi: 10.1155/2016/3429604
Riquelme, B., Foresto, P., D'Arrigo, M., Valverde, J., and Rasia, R. (2005). A dynamic and stationary rheological study of erythrocytes incubated in a glucose medium. J. Biochem. Biophys. Methods 62, 131–141. doi: 10.1016/j.jbbm.2004.10.004
Schindelin, J., Arganda-Carrera, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, S., et al. (2012). FIJI: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682. doi: 10.1038/nmeth.2019
Shilov, A. M., Avshalumov, A. S., Markovsky, V. B., Sinitsyna, E. N., and Poleshchuk, O. I. (2008). Changes of blood rheological properties in patients with metabolic syndrome. Russ. Med. J. 4, 200–204.
Stroev, E. A., Ryazanova, E. A., and Tavintsev, V. D. (1991). Method for Determining the Activity of Calpains in Biological Material. Certificate of authorship 1668950, Registered in the State Register of inventions of the USSR 04/08/. Application ◦4699987, priority of 05.06.1989.
Yusipovich, A. I., Zagubizhenko, M. V., Levin, G. G., Platonova, A., Parchina, E. Y., Grygorzcyk, R., et al. (2011). Laser interference microscopy of amphibian erythrocytes: impact of cell volume and refractive index. J. Micro. 244:223–229. doi: 10.1111/j.1365-2818.2011.03516.x
Keywords: human erythrocytes, hyperglycaemia, hemoglobin system, membrane lipids, proteolytic activity
Citation: Revin VV, Klenova NA, Gromova NV, Grunyushkin IP, Solomadin IN, Tychkov AY, Pestryakova AA, Sadykhova AV, Revina ES, Prosnikova KV, Bourdon J-C and Zhelev N (2017) Physical and Chemical Processes and the Morphofunctional Characteristics of Human Erythrocytes in Hyperglycaemia. Front. Physiol. 8:606. doi: 10.3389/fphys.2017.00606
Received: 31 May 2017; Accepted: 07 August 2017;
Published: 30 August 2017.
Edited by:Sergei G. Sokolovski, Aston University, United Kingdom
Reviewed by:Natalia Pavlovna Akentieva, Institute of Problems of Chemical Physics (RAS), Russia
Irma Borisovna Mosse, Institute of Genetics and Cytology at Belorussian Academy of Sciences, Belarus
Copyright © 2017 Revin, Klenova, Gromova, Grunyushkin, Solomadin, Tychkov, Pestryakova, Sadykhova, Revina, Prosnikova, Bourdon and Zhelev. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.