Artículo original

Software for cerebrospinal fluid diagnostics and statistics

Software para el diagnóstico y estadísticas del líquido cefalorraquídeo

0000-0001-7210-5109Hansotto Reiber^{1  *} 

¹Georg August Universität. Göttingen, Germany.

ABSTRACT

Background:

In CSF analysis for diagnostics we have knowledge-based software for numerical and graphical data interpretation, but software programs for statistics are scarce. Free, stand-alone software programs that calculate all individual functions of CSF protein analysis and allow the statistical treatment of groups of diseases numerically and graphically are presented for relevant examples.

Methods:

Diagnosis of an intrathecal synthesis refers to the upper limit of the reference range, Qlim = Qmean +3SD, but statistical evaluation of its frequency is referred to Qmean+2SD. When quantifying intrathecal synthesis for statistics, either the absolute amount (IgG_loc) or the relative intrathecal fraction (IgG_IF) can be reported with reference to the mean reference curve, Qmean. The free software “CSF research Tool” for immunoglobulins allows diagnostic and statistic evaluations with Reibergrams and calculation of mean values and standard deviations from disease groups. The software “FLC-K statistics” for free light chains Kappa offers for diagnostics and statistics the numerical and graphical interpretation basis for statistical processing in exported Excel tables. A free “CSF-App” for Smartphones provides data calculation for diagnostics of single patients with examples of disease-related data patterns.

Results:

Patients with clinically isolated syndrome (CIS) who were later diagnosed as MS showed no immunological differences to patients initially diagnosed as MS (same mean quantity of intrathecal synthesis in CIS and MS detectable for IgG and FLC-K). The frequently claimed diagnostically higher sensitivity of the FLCK analysis compared to IgG, can be explained by the up to 3-fold higher mean intrathecal fraction of FLC-K, corresponding to a higher frequency in the detection of intrathecal synthesis with FLCK analysis.

Conclusions:

With a knowledge-based quantification in CSF analysis, supported by knowledge-based software programs, scientifically and diagnostically important results can be obtained.

Key words: Cerebrospinal fluid; Statistics; Software; CSF-App; Clinically Isolated Syndrome; Multiple Sclerosis; IgG; Free light chains Kappa

RESUMEN

Introducción:

Los programas de software gratuitos y autónomos que calculan todas las funciones individuales del análisis de proteínas del líquido cefalorraquídeo (LCR) y permiten el tratamiento estadístico de grupos de enfermedades de forma numérica y gráfica se presentan como ejemplos relevantes.

Métodos:

Cuando se cuantifica la síntesis intratecal para la estadística, se puede informar la cantidad absoluta (IgG_loc) o la fracción intratecal relativa (IgG_IF) con referencia a la curva de referencia media, Qmean. El software gratuito “CSF research Tool” para inmunoglobulinas permite realizar evaluaciones diagnósticas y estadísticas con Reibergrams y calcular los valores medios y las desviaciones estándar de los grupos de enfermedades. El software “FLC-K statistics” para Free light chains Kappa ofrece para el diagnóstico y la estadística la base de interpretación numérica y gráfica para el procesamiento estadístico en tablas exportadas de Excel. El programa “CSF-App” para teléfonos inteligentes es gratuito y ofrece el cálculo de datos para el diagnóstico de pacientes individuales con ejemplos de patrones de datos relacionados con enfermedades.

Resultados:

Los pacientes con síndrome clínico aislado (SCA) que posteriormente fueron diagnosticados como EM no mostraron diferencias inmunológicas con respecto a los pacientes inicialmente diagnosticados como EM (la misma cantidad media de síntesis intratecal en el síndrome clínico aislado y EM detectable para IgG y FLC-K). La sensibilidad más elevada que se afirma con frecuencia en el diagnóstico del análisis de FLC-K en comparación con la IgG, puede explicarse por la fracción intratecal media hasta tres veces mayor de FLC-K, que corresponde a una mayor frecuencia en la detección de la síntesis intratecal con el análisis de FLC-K.

Conclusiones:

Con la cuantificación en el análisis del LCR se pueden obtener resultados importantes desde el punto de vista científico y diagnóstico.

Palabras-clave: líquido cefalorraquídeo; estadística; software; CSF-App; síndrome clínico aislado; esclerosis múltiple; IgG; cadenas ligeras libres kappa

Introduction

The importance of qualified CSF diagnostics has been described many times for the differential diagnosis of neurological diseases of the individual patient.^1,2,3,4,5 For this purpose, there is also good laboratory software that combines on-line and off-line data into an integrated report of the patient's findings.^3,6,7,8. The graphical representation of the immune response in quotient diagrams (Reibergram) as part of this report⁵ allows the rapid recognition of typical patterns of disease. The theoretical principles for numerical and graphical evaluation with nonlinear reference range were also reported in detail.^9,10,11,12

However, what has hardly been described since then is the statistical evaluation of groups of diseases on this theoretically sound basis for CSF analysis. However, free software programs have been available for some time¹³ which have recently been expanded. Without such a statistical treatment of patient cohorts, an important research tool would not be available. A knowledge-based quantification of the intrathecal synthesis can be helpful for the comparison of the immune response of different disease groups with controls, but also in clinical studies where the effectiveness of therapeutic measures is important. To improve diagnosis, it may be important to compare the diagnostic sensitivity by the relative intrathecal fractions of different parameters, such as IgG and free light chains.¹² As an application in psychiatry,¹⁴ it was helpful to exclude an inflammatory cause, e.g., in schizophrenia or bipolar disorders.¹⁴

The still frequently used diagnostic interpretation aid, the linear IgG index,¹⁵ does not allow to make a trustworthy statistic, this is only possible on the basis of a non-linear, hyperbolic reference range.^10,12 For the evaluation with nonlinear reference range, two free software packages for immunoglobulins and for the free light chains Kappa are described in this review.¹⁶ The application will be demonstrated by means of two current problems of multiple sclerosis research: The diagnosis of definite MS and clinically isolated syndrome (CIS)¹⁷ as well as the analytical sensitivity of free light chain kappa (FLC-K) with that of IgG to detect an intrathecal immune response.¹²

The numerical and graphical programs for statistics have recently been supplemented by a so-called CSF app, an interpretation software for smartphones.¹⁶

Methods

CSF statistics

The intrathecal synthesis in disease groups can be quantitatively characterized and compared by two different methods.

The number of patients with intrathecal synthesis per group

In Fig. 1 the limits of the reference range for IgG are shown in a linear quotient diagram. For the diagnosis the upper limit Qlim of the reference range with Qmean +3SD (s = standard deviation) is used. For statistical purposes, however, it is advisable to refer to an upper limit of the reference range with Qmean +2 SD. The reference range with this upper limit (dashed line in Fig.1) covers only 96% of the controls but results in a higher sensitivity. An example from the software-based calculation for a group of MS patients¹⁸ is shown in Table 1 to illustrate the difference in sensitivity and specificity between the two limits. The unit of measurement is either a dimensionless fraction in n/N or in percent n/N x100 [%].

Fig. 1 Hyperbolic reference range of blood-derived proteins in CSF. The linear diagram shows the upper diagnostic limit (Qlim) and the mean function Qmean of the hyperbolic reference range.¹⁰ Qlim is defined by QLim (3s) = Qmean + 3 SD. The statistically relevant upper border line (dashed line) is defined by Qlim(2s) = Qmean + 2SD. Differences in the diagnostic sensitivity of the different Qlim lines are shown in Table 1. 

Table 1 Reference range related frequency of intrathecal IgG, IgA and IgM synthesis in a group of MS patients 

Reference	IgG	IgA	IgM
Qlim (3s)	23	2	6
Qlim (2s)	28	5	8
N	29	29	29

Note: The data originate from an earlier study.¹⁸ The cohort of patients with definite MS had an age of 28-29 years at first diagnostic puncture. According to Fig 1 either the diagnostic reference limit QLim = Qmean + 3 s or the statistics reference with Qlim (2s) = Qmean + 2s was used. s = Standard deviation, SD. Calculated with the software, CSF research Tool (http://www.albaum.it).

The mean intrathecal synthesized amount of protein.

For the statistically usable quantification of the mean quantity, the mean value of the reference range Qmean must be referred to. Qmean has the formula in (Eq. 1).

(Eq. 1)

The values for a, b, c depend on the molecule size with examples in Table 2.^10,11,12

Table 2 Parameters of the mean hyperbolic reference function for IgG , IgA , IgM¹⁰ and free light chain kappa, FLC-K.¹²  

Parameters	a/b	b^2 (x 10^6)	c (x 10^3)
QIgG (mean)	0,65	8	1,4
QIgA (mean)	0,47	27	2,1
QIgM (mean)	0,33	306	5,7
QKappa (mean)	1.95	33	4,85

Note: The parameters chararacterize with s = a/b the slope and the ordinate interval (c) of the asymptote of the hyperbolic functions with the parameters a/b , b², c.

The interpretation of the intrathecal quantity can refer either to the absolute intrathecal quantity or to the relative fraction of the total concentration in CSF.

Absolute amount synthesized

The Qmean related local (intrathecal) synthesis as absolute concentration, e.g. IgG_loc (mean) or K_loc (mean) for the Free Light Chain Kappa, is calculated with equation (Eq 2) for each individual patient and then the statistical mean is calculated for the whole collective.

(Eq. 2)

Relative amount synthesized

For the diagnostic sensitivity of a parameter, the relative proportion of its total cerebrospinal fluid concentration (K_CSF), the intrathecal fraction, e.g. IgG_IF or K_IF, calculated with Eq 3, may be of interest. Here, too, the statistical averaging in the collective for IgG_IF(mean) or K_IF(mean) is found from the calculated individual values.

Other

(Eq.3)

Software

The software “CSF research Tool/Reibergrams”^13,16 for immunoglobulins calculates the IgGloc (mean) and IgGIF (mean) values for the data transferred from Excel tables and offers their statistical evaluation to a selected group with mean values and coefficients of variation directly (Table 3). This is linked to a number of other statistical evaluations, such as the frequency of a barrier disturbance, etc. With this software, the data are also displayed graphically in rubbing programs (Fig. 2).

Fig. 2 Quantitation of protein data in Quotient diagrams (Reibergrams) for IgG (left) and freelight chains Kappa (right). Filled circles represent patients diagnosed with definite MS. Open circles represent patients initially diagnosed as clinically isolated syndrome (CIS), later diagnosed as definite MS. The data for IgG and FLC-K are from the same kohorts of CIS or MS patients. Corresponding statistics data are shown in Table 3. The graphs are developed with the software “CSF research tool” and “FLC-K statistics and graphic program” (http://www.albaum.it). 

Table 3 Statistics data of cohorts of MS and CIS patients. Definite MS is diagnosed according to the Mc Donald criteria 

Note: CIS patients have been diagnosed later in the Course of disease as a definite MS. Data were calculated with the free software “CSF research tool” and “FLC-K statistics and graphic program” from http://www.albaum.it

The software “FLC-K statistics and graphic program”¹⁶ for Free Light Chains Kappa, which is intended as an extension, calculates the Kloc (mean) and KIF (mean) values for the data taken from Excel tables. After export to Excel tables, the mean values can be calculated. The graphical results are shown in Fig. 2 as a friction program for FLCK.

The “CSF- App”¹⁶ was developed for Android-based smartphones. It offers calculation of albumin, IgG, IgA, IgM and FLC-K quotients for a complete CSF data report with graphical results in Reibergrams.⁵ The Antibody Index calculation for the MRZ -Antibody reaction¹⁹ are provided as well. As a tutorial this program offers examples for the possible data interpretations and an archive of immunoglobulin patterns from 15 different neurological diseases with comments for differential diagnosis.

Patients

The patients data and diagnosis originate from the laboratory of Zeman D. and Kušnierová P. as described in their earlier publications.¹² Multiple Sclerosis diagnosis was made according to the current modified Mc Donald criteria.²⁰

Results

Multiple sclerosis and clinically isolated syndrome

Patients with clinically based suspicion of multiple sclerosis who do not meet the McDonald criteria²⁰ are diagnosed as clinically isolated syndrome (CIS).¹⁷ For various authors, the difference between CIS and multiple sclerosis (MS) is associated with the idea of immunological conversion.²¹ This can be investigated by quantifying the intrathecal immune response in both groups. In addition to the detection of intrathecal IgG, the more sensitive FLC-K analysis in CSF is increasingly used.^12,17 However, quantification can also provide the specific information why the quantitative detection of FLCK is more sensitive than that of IgG. The CIS group investigated here represents patients who later became definite MS after initial CIS diagnosis.

By comparing the mean quantitative intrathecal synthesis amount (IgGloc mean and K loc mean in mg/l) between MS and CIS (Table 3) it becomes clear that no significant quantitative differences in the sense of an immunological conversion from CIS to MS can be detected. It rather gives the impression that the CIS group has slightly higher values on average.

With the analysis of the mean intrathecal synthesized fraction (IgGIF(mean) and KIF (mean) in %) we gain a further insight. The mean intrathecal fraction of FLC-K is almost three times larger than for IgGIF (Table 3 and Fig 2). This means that the relative fraction synthesized in the brain has a higher percentage of the total amount in cerebrospinal fluid for FLC-K than for IgG. The use of FLC-K instead of IgG makes a significant difference in diagnostic sensitivity. However, an important prerequisite is to refer to the nonlinear reference range.¹²

However, we have another analytical criterion, which is the frequency of synthesis in the two groups with reference to the statistical upper limit, Qlim (2s) (Table 3 and Fig 2). For IgG, only 25/45 patients (MS) or 37/50 (CIS) are above the statistically relevant limit ( QIgG > Qlim(2)), whereas for the free light chains 42/45 (MS) or 49/50 (CIS) are above the limit (QK > Qlim (2)). This confirms the theoretically superior sensitivity of FLC-K to IgG for the diagnosis of an inflammatory process. These ratios are documented numerically in Table 3 and graphically in the rubbing programs of Fig. 2.

Conclusion: The quantification of intrathecal protein synthesis provides important information for diagnosis and pathophysiology.

Discussion

Numerical and graphical CSF data presentation

The graphical representation of the immunoglobulin response as part of an integrating CSF data report was introduced to identify at a glance possible disease-typical immunoglobulin patterns.^1,5 The “CSF-App”¹⁶ offers a collection of such disease-typical patterns with differential diagnostic comments.

The reference range for serum proteins in CSF in the quotient diagrams (Reibergram) is non-linear,^10,11,12 i.e. it is characterized by empirically and theoretically based hyperbolic functions (Fig 2). Unless these diagrams are automatically produced with the laboratory software,^6,7,8 it is difficult for many users to create such a diagram with a reference range. This has led to the fact that many authors still prefer the evaluation with linear index, although false positive and false negative interpretations are unavoidable.¹ In order to remedy this, the stand-alone solutions of the described software programs were developed.¹⁶⁾ Data for individual patients as well as for groups of patients can be created in a publication-ready form. The software programs also offer data processing to obtain all diagnostically relevant values (e.g., intrathecal synthesis quantity) numerically. This also raises the question of the statistical use of these data of disease groups. Since not every parameter that can be obtained is suitable for statistical treatment (see Methods), all statistically relevant data are additionally calculated and can be exported as Excel tables for extended statistics.

Reference range for statistics

The limit used in the rubbing programs, Qlim (Fig 1 and 2), provides the highest specificity to diagnose an inflammatory process (smallest number of false positive interpretations). At this limit 99% of the negative controls are integrated in the reference range, but a correspondingly lower sensitivity for inflammatory processes is accepted. In statistics, a variance of only mean ± 1 SD is usually used for group comparisons. A comparison of mean ± 2 SD increases the accuracy for smaller groups or groups with a high analytical imprecision (error propagation in quotients). The reference range for this statistical reference curve selected here (dashed line in Fig 1) still covers 96% of the controls (i.e. only max 2% false positive controls).

Intrathecal fraction

The intrathecal fraction, IF, was introduced to compare the intrathecal synthesis in the quotient diagram of different proteins with different total amounts in serum.^10,11 This value was related to Qlim (Fig. 1). These calculations in the Reibergrams represent the theoretical minimum of the intrathecal synthesis of the different immunoglobulin classes. For the statistical comparison of groups, the calculation of mean synthesis quantities must of course be related to Qmean, the group mean value (Fig. 1). From the IF data (Table 3) it is clear that the smaller the normal serum concentration and accordingly the normal CSF concentration is (100-fold difference between IgG and FLC-K), the larger the relative concentration of the intrathecal fraction (Eq.3). This also means that the detection sensitivity for intrathecal synthesis increases (e.g. by almost 3 times the difference between IgG and FLCK in Table 3).

Extended statistics

The individual data for Igloc(mean) or Kloc(mean) calculated in the software can be exported into an Excel table and are thus available for any desired form of statistics with appropriate significance analysis. The “CSF research tool”¹⁶ also provides the data for testing the null hypothesis.

It should be noted that only Igloc(mean) and Kloc(mean) values are suitable as absolute concentrations for group comparisons. The IgIF (mean) and KIF (mean) values as relative values are only suitable for parameter comparisons.

Immune response in MS and CIS

The Gold Standard for the detection of an intrathecal IgG synthesis is the qualitative analysis of oligoclonal IgG bands. But this method does not allow a quantitative discrimination between different disease groups like CIS and MS. The quantitative analysis in Quotient diagrams allows to characterize at least the extent of intrathecal synthesis (intrathecal fraction in %) for the individual patient (like dominance of immunoglobulin glasses or difference in IgG, IgA, IgM patterns ^1,5). To compare reliably differences or similarities between groups we need statistics. For the chosen example of CIS and MS we get results from which we can conclude that the integration of the intrathecal fractions of FLC-K in the Mc Donald criteria could lead to an earlier more sensitive diagnosis of MS albeit the most specific parameter for detection of an intrathecal immune response, the MRZ antibody synthesis¹⁹ is increasingly investigated for MS diagnosis.²²

References

1.  Reiber H, Peter JB. Cerebrospinal fluid analysis - disease-related data patterns and evaluation programs. J Neurol Sci. 2001;184:101-22.

2.  Dorta-Contreras AJ. Reibergramas: elemento esencial en el análisis inmunológico del líquido cefalorraquídeo. Rev Neurol. 1999;28:996-8.

3.  Regeniter A, Kuhle J, Mehling M, Möller H, Wurster U, Freidank H, et al. A modern approach to CSF analysis: Pathophysioly, clinical application, proof of concept and laboratory reporting. Clin Neurol Neurosurg. 2009;111(4):313-8.

4.  Akaishi T, Narikawa K, Suzuki Y, Mitsuzawa S, Tsukita K, Kuroda H, et al. Importance of the quotient of albumin, quotient of immunoglobulin G and Reibergram in inflammatory neurological disorders with disease-specific patterns of blood-brain barrier permeability. Neurol Clin Neurosci. 2015;3(3):94-100. DOI:10.1111/ncn3.158

5.  Reiber H. Cerebrospinal fluid data compilation and knowledge-based interpretation of bacterial, viral, parasitic, oncological, chronic inflammatory and demyelinating diseases: Diagnostic patterns not to be missed in Neurology and Psychiatry. Arq Neuropsiquiatr. 2016;74(4):337-50.

6.  Albaum W. CSF-Lab. Labororganisation für das Liquorlabor; 2020 [access: 21/01/2020]. Available at: Available at: http://www.albaum.it 6.

7.  Wormek A. Desktop application as knowledge-based system for cerebrospinal fluid (CSF) analysis; 2020 [access: 21/01/2020]. Available at: Available at: http://www.wormek.org 7.

8.  Protis system. Automated CSF assessment; 2020 [access: 21/01/2020]. Available at: Available at: http://www.siemens-healthineers.com/diagnostics-it 8.

9.  Reiber H. Knowledge-base for interpretation of Cerebrospinal fluid data patterns - Essentials in Neurology and Psychiatry. Arq Neuropsiquiatr. 2016;74(6):501-512.

10.  Reiber H. Flow rate of cerebrospinal fluid (CSF) - a concept common to normal blood-CSF barrier function and to dysfunction in neurological diseases. J Neurol Sci. 1994;122:189-203.

11.  Reiber H. Proteins in cerebrospinal fluid and blood: Barriers, CSF flow rate and source-related dynamics. Rest Neurol Neurosci. 2003;21:79-96.

12.  Reiber H, Zeman D, Kušnierová P, Mundwiler E, Bernasconi L. Diagnostic relevance of Free light chains in Cerebrospinal fluid - The hyperbolic reference range for reliable data interpretation in quotient diagrams. Clin Chim Acta. 2019;497:153-62.

13.  Reiber H, Albaum W. Statistical evaluation of intrathecal protein synthesis in CSF/Serum quotient diagrams. Acta Neuropsychiatrica 2008; 20 (Supplement 1): 48-49

14.  Bechter K, Reiber H, Herzog S, Fuchs D, Tumani H, Maxeiner HG. Cerebrospinal fluid analysis in affective and schizophrenic spectrum disorders. Identification of subgroups with immune responses and blood-CSF barrier dysfunction. Journal Psychiatric Research. 2010;44:321-30.

15.  Link H, Tibbling G. Principles of albumin and IgG disorders. Evaluation of IgG synthesis within the central nervous system in multiple sclerosis. Scan J Clin Lab Invest. 1977;37:397-401.

16.  Albaum W. CSF -App; CSF Research Tool/Reibergrams; FLC-K Statistics and Graphic program; 2020 [access: 12/04/2020]. Available at: Available at: http://www.albaum.it 16.

17.  Pesslauer S, Milosavljevic D, Huebl W, Aboulenein-Djamshidian F, Krugluger W, Deisenhammer F, et al. Validation of kappa free light chains as a diagnostic biomarker in multiple sclerosis and clinically isolated syndrome: A multicenter study. Mult Scler. 2016; 22:502-10.

18.  Reiber H, Teut M, Pohl D, Rostasy K.M, Hanefeld F. Paediatric and adult Multiple Sclerosis: Age Related Differences and time course of the neuroimmunological response in Cerebrospinal fluid. Multiple Sclerosis. 2009;15:1466-80.

19.  Reiber H. Polyspecific antibodies without persisting antigen in multiple sclerosis, neurolupus and Guillain-Barré syndrome: immune network connectivity in chronic diseases. Arq Neuropsiquiatr. 2017;75(8):580-8. DOI: 10.1590/0004-282X20170081

20.  Polman CH, Reingold SC, Banwell B, et al. Diagnostic criteria for multiple sclerosis: 2010 revisions to the McDonald criteria. Ann Neurol. 2011;69(2):292-302.

21.  Villar LM, Espiño M, Costa-Frossard L, Muriel A, Jiménez J, Álvarez-Cermeño JC. High levels of cerebrospinal fluid free kappa chains predict conversion to Multiple sclerosis. Clin Chim Acta. 2012;413:813-6.

22.  Jarius S, Eichhorn P, Franciotta D, Petereit HF, Akman-Demir G, Wick M, Wildemann B. The MRZ reaction as a highly specific marker of multiple sclerosis: re-evaluation and structured review of the literature. J Neurol. 2017;264(3):453-66. DOI: 10.1007/s00415-016-8360-4

Received: April 15, 2020; Accepted: April 30, 2020

^*Correo electrónico: [e-mail]

The author declares no conflict of interest.