METHOD OF DIAGNOSING OR CATEGORIZING DISORDERS FROM BIOCHEMICAL PROFILES

Abstract: A method for diagnosing disorders in living organisms, in which fluid samples from normal and afflicted (abnormal) individuals are analyzed to generate patterns representative of molecular constituents of said samples. A data base of frequency distribution patterns of constituents of samples from organisms having known categories of disorders and controls is created, and the unknown sample analysis is compared for conformity to said frequency distribution patterns. The invention has particular applicability to diagnosing degenerative diseases such as Alzheimer's Disease, Parkinson's Disease, Huntington's Disease, schizophrenia, Progressive Supernuclear Palsy, amyotrophic lateral sclerosis and senile dimential, which are given as exemplary. The invention also may be advantageously employed to diagnose diseases such as tumors, carcinomas, cardiovascular abnormalities and other disorders, or for selection of the therapy based on categories of known vs. unsuccessful outcomes. Moreover, both treatment protocols and new pharmaceuticals may be evaluated.

Description: METHOD OF DIAGNOSING OR CATEGORIZING DISORDERS FROM BIOCHEMICAL PROFILES This invention relates to analytical and mathematical methods for diagnosing or categorizing disorders. The invention has particular utility for diagnosing or categorizing disorders in living animals from analysis profiles of biologically active materials such as neurotransmitters and other neurochemical substances in brain tissue, cerebrospinal fluid, plasma, serum, saliva, nasal mucosa, urine and the like, such as catecholamines, their precursors, cofactors and their metabolites. The invention is uniquely capable of differentiating a large number of compounds of biological, diagnostic and/or pharmaceutical significance and of using such differential for diagnosing disorders and will be described in connection with such utility although other uses are contemplated. There is an extensive body of literature relating abnormalities in neurotransmitters, precursors, and metabolites to degenerative, neuropsychiatric and behavioral disorders, hypertension and certain carcinomas. See, for example, Schildkraut et al in The Brain, Biochemistry and Behavior, Proceedings of the Sixth, Arnold 0. Beckman Conference in Clinical Chemistry, pages 47-68. Although the potential role of these compounds in a number of significant disorders has been established, their routine analysis has not yet achieved widespread clinical use. Two problems in the clinical utility of neurotransmitter measurements are related to the economic and technical limitations of current technology. First, there is felt to be a high degree of interlaboratory and intersample uncertainty in quantitative values. Second, it has been difficult to measure enough of the known metabolically related compounds of a particular neurotransmitter to fully describe its biochemical significance in an individual sample, or to detect, identify and measure unusual neurotransmitters - an important aspect of basic research in various disease states that is presently very expensive and specialized. While a number of interlaboratory technique intercomparisons for a variety of neurotransmitters have been carried out, there has been no comprehensive study within and among different techniques and laboratories for neurotransmitters in typical samples of interest. In the absence of such studies, given the complexity of the analytical problem and the historically wide variation whenever an analyte has been subjected to rigorous interlaboratory testing, the current values for normal and abnormal neurotransmitter levels must be take with unspecified and probably wide limits of confidence. Although the analysis of single neurotransmitters or metabolites from a complex biochemical pathway has been shown to correlate with a number of disorders utilizing statistical analysis over a large number of samples, the analytical level of a single neurotransmitter in an individual sample, with a few exceptions, has had relatively low clinical diagnostic utility. Essentially the state of the field of biochemical correlates of disorders is that while between large populations of normal and abnormal individuals a correlation generally can be determined for a particular biochemical, the scatter that results from both analytical and biochemical phenomena typically does not permit the level of a particular biochemical to be utilized diagnostically for a particular single individual. Nor may a single biochemical value be utilized for the rational prescription or development of a pharmaceutical for that individual. This is not particularly surprising in that both the levels and effects of a particular neurotransmitter are modified by a number of other neurotransmitters, in the same, or parallel metabolic pathways. If, for instance, 5-HT (serotonin) is to be used as a diagnostic tool for depression, suicidal tendencies, or schizophrenia, it would be necessary and perhaps provide a route to definitive diagnosis and pharmaceutical specification or development, to simultaneously determine the approximately 40 other compounds that derive from tryptophan and significantly effect the indolaminergic neuronal system's activity. In recent years, LCEC (Liquid Chromatography with Electrochemical Detection) has become a common tool for the determination of catecholamines biogenic amines and their metabolites in biological fluids. Because of sensitivity limitations (typically 20-50 pg) and the complexity of biological samples, both separation and concentration steps typically have been necessary. Heretofore, plasma catecholamine analysis typically required three steps. First, the sample is collected and the catecholamines separated and concentrated, for example, using the alumina extraction procedure of Anton and Sayre (See A.H. Anton and D.F. Sayre, J. Pharmacol, Exp. Ther., 138 (1962), p. 360-375). The analytes, norepinephrine, epinephrine and dopamine, along with the internal standard DHBH (dihydroxybenzylamine) , then are separated chromatographically, and finally detected electrochemically. Typical sample size requirements are 1.0 ml plasma or serum. In routine clinical use, there have been numerous problems with conventional techniques (alumina absorption, ion exchange and extraction), due to a large number of poorly understood variables, in the overall analysis system of sample acquisition, storage, preparation and sensor response. These problems have quite likely confused the relationships that may exist between levels and distribution of the catecholamines and various physiological and behavioral phenomena and disease states. In the analysis of complex biological materials such as blood, serum and cerebrospinal fluids which may contain numerous different constituents, the important (e.g. abnormal) metabolites such as neurotransmitters to be identified may be present in only parts per trillion. While a chromatographic column can achieve macro separation of the various constituents, it may not provide adequate spatial (in time) separation of the extremely small portion of metabolites of interest from the much larger percentage of the many other compounds coeluted from the column at the same time as the metabolites of interest. Many of these interfering coeluted materials are electrochemically active but electrochemically irreversible, while many metabolites such as neurotransmitters are both electrochemically active and electrochemically reversible. It has been found that the analytical problems of reliable measurements of neurochemicals and related compounds are complicated by the fact that interferences with conventional or prior technologies are disorder related. This was discussed in my prior publication, (Matson et al. Clinical Chemistry, Vol. 30, No. 9, 1984) (see U.S. Patent 4,511,659) for dopamine, dopac and seratonin measurements in directly analyzed brain extract and cerebrospinal fluid for normal, schizophrenics and Alzheimers. Recent work has indicated that even for the widely used and accepted technique of alumina extraction for plasma catecholamines that inferences may be disorder specific. Higher values for Norepinephrine (NE) and Epinephrine (EP) were observed following alumina extraction and analysis of a single energy conventional electrochemical detector than for a three cell redox detector on neonatal stress blood samples. Analysis of the neonate extracts on the sixteen channel chemical imaging system revealed several unexpected compounds that are potential interferences including dihydroxyphenylacetic acid (DOPAC) , 3 hydroxykynurenamine (3-OHKYA) and 3-hydroxy-anthranilic acid (3-OHAN). These compounds have not been detected in normal adult plasma alumina extracts. In my aforesaid U.S. Patent No. 4,511,659, there is provided an electrochemical detection system comprising a plurality of coulometrically efficient electrochemical cells, in series, for sequentially oxidizing and reducing selected substances in a sample solution under controlled conditions prior to measurement on a downstream testing electrode or electrodes. More specifically, in accordance with the invention provided in my aforesaid U.S. Patent No. 4,511,659, a sample solution (e.g. a body fluid) is passed through a suitable chromatographic column and the eluant is streamed in contact with a series of electrochemically isolated, in-line coulometric electrodes operated under conditions so as to establish a series of "gates" for the sequential oxidation and reduction of substances in the sample solution whereby to screen (remove) selected interfering and electrochemically irreversible substances contained in the sample solution, while passing selected electrochemically reversible products for detection and measurement on a downstream electrode. The gate electrode series is follows in-line by one or more, preferably an array of six or more coulometric measuring electrodes, each formed of porous electrode base material such as fritted graphite, fritted carbon or other conductive fritted material, for detecting and measuring the electrochemically reversible compounds of interest (e.g. neurotransmitters). As reported in my aforesaid U.S. Patent No. 4,511,659, there are several beneficial effects of this approach to electrochemical analysis. Long-term drift in response is effectively eliminated by acquiring essentially 100% of the signal. The capability of analyzing essentially 100% of a material allows the assay of compounds of unknown purity by relating them to the basic principles of electrochemical reaction embodied in Faraday's Law. Poisoning of the electrode, a dominant problem with electrochemical sensors, is effectively eliminated by the use of a much larger relative surface area for reaction. And, finally, and most important to the eventual development of array and gate cells, a coulometric electrode by virtue of its essentially 100% efficiency allows sequential oxidation and/or reduction of compounds at successive-in-line detectors. The improved sensitivity of the detection system as discussed in my aforesaid U.S. Patent No. 4,511,659, particularly where two or more active testing electrodes follow the screening electrodes has given the ability to do direct injections of serum filtrates and has also allowed the generation of reproducible patterns of compounds with catecholamine like electrochemical behavior of a large number of resolvable components. This provides the possibility of performing pattern recognition for the diagnosis or perhaps even predictive diagnosis, of various disorders or disease states. In my copending application Serial No. 797,615 and its parent U.S. Patent 4,863,873, I describe a system for resolving and detecting hundreds of compounds in a single sample at femtogram levels whereby to provide a small molecule inventory or metabolic pathway pattern of an individual. As taught in my aforesaid U.S. Patent 4,863,873, the small molecule inventory may be considered to reflect the underlying activity and distribution of the enzymatic pathways of an individual and hence reflect an operational measure of the genone determining those enzymes. The small molecule inventory of an individual may thus be used to determine the health state of the individual and/or to diagnose disease states. Correlation of the patterns from a plurality of individuals provides an understanding of the mechanisms of disorders or disease states or conditions and, in turn, provides a rational route to pharmacological development leading to treatment, cure or suppression of such disorders, disease states or conditions. The present invention is an improvement in the invention described in my aforesaid U.S. Patent 4,863,873. More particularly, in the practice of my invention as described in my U.S. Patent 4,863,873, I have observed that the biochemical analysis profiles of "normal" or healthy individuals may vary quite widely, while the biochemical profile analysis data of individuals having disorders is far less chaotic. More particularly, I have observed that the frequency distribution of certain biochemical compounds or ratios of compounds in individuals suffering from a disorder are far less chaotic than "normal" or healthy individuals. This leads to a general protocol for diagnosing, categorizing or differentiating individuals based on comparisons of biochemical analytical data of small molecule inventory against data bases of known or previously diagnosed cases. By way of example the process of the present invention may advantageously be employed in the differentiation of neurological degenerative dementing or affective disorders such as Alzheimer's disease, Huntington's disease, Parkinson's disease, schizophrenia. Progressive Supernuclear Palsy, amyotrophic lateral sclerosis and senile dementias from each other and neurologically normal controls. Moreover, by suitable selection of variables, the process of the present invention also is applicable to classification of tumors, carcinomas, cardiovascular abnormalities and other disorders. Similarly, the process of the present invention advantageously may be utilized to select therapy based on categories of known successful vs. unsuccessful outcomes. While not wishing to be bound by theory, the two fundamental hypotheses underlying the process of the present invention are: 1. The underlying genetic makeup or predisposition of an individual will reflect through the proteins, enzymes, and other factors it determines in patterns of small molecules. Individual components within these patterns will be affected by environmental effects such as diet, stress or chemical inset; however, the overall pattern of relationships will reflect the underlying operation of the genome or the interference of a particular disorder. .Among the small molecules are the transmitters, cofactors and metabolites that regulate neuronal and endocrine functions and the interactions of somatic and central nervous system processes. Thus, the compounds such as purines, tyrosine and tryptophan derived neurotransmitters, peptides, pterin and vitamin cof ctors are highly " relevant to the effect or etiology of neurological disorders, cardiovascular disfunction and certain tumors or carcinomas. 2. The relationships of these biochemical patterns from a disorder are less chaotic or more regular than those from healthy controls. All of the biochemical systems of small molecules are interconnected and interrelated in a complex web of feedback and response. These interactions are highly nonlinear and thus, depending on subtle differences in initial conditions, the response of individual components in a biochemical pattern will be highly variable. The overall system will thus behave in a mathematically chaotic fashion. In a disorder, elements within the biochemical pattern are over or underregulated, thus reducing the degrees of freedom or overall variability. Consequently, the presence of a disorder implies more regulated or less chaotic variability of compounds or relationships among compounds in patterns from disordered individuals. These two fundamental hypotheses provide an approach to diagnostic categorization of disorders using frequency distributions of compounds and relationships from large data bases. Description of the Drawings For a fuller understanding of the nature and objects of the present invention, reference should be had to the following detailed description taken in combination with the accompanying drawings wherein: Figs. la-Id are frequency distribution graphs of variables of Alzheimer's Disease and controls; Fig. 2 is a geographical representation of scoring algorithms; Fig. 3 is a plot showing initial scoring of Alzheimer's Disease v. control; Fig. 4 is a plot similar to Fig. 3, with five Alzheimer's Disease cases which scored as controls removed from the Alzheimer's Disease scoring Distribution; Fig. 5 is a plot showing analog distribution measurements in accordance with the present invention; Fig. 6A is a plot of measurements of nasal mucosa swab samples at low gain; Fig. 6B is a plot of nasal mucosa swab samples at high gain; and Fig. 7 is a plot of 8-amyloid in saline suspension of mucosa swab. Detailed Description of the Invention EXAMPLE I Methodology for Sample Analysis and Data Base Creation 280 CSF samples from the categories Alzheimer's Disease - AD (61 samples), Parkinson's Disease - PD (60 samples), schizophrenia - SC (60 samples), Huntington's Disease - HD (20 samples), Supernuclear Palsy - PSP (13 samples) and neurologically normal controls - C (68 samples), were electrochemically analyzed in accordance with the teachings of my aforesaid U.S. Patent 4,863,873. Samples from normal and diseased individuals were prepared and flowed through a chromatographic column, and detected in a 16 sensor electrochemical cell using an NCA Chemical Analyzer, Model No. CEAS available from ESA, Inc., Bedford, Massachusetts. Sensor potentials ranged from T^ -600mv to T^ς +900mv in lOO v increments. All samples were from 7th or 8th mL aliquots of rostral caudal gradients. Pools were created for each category utilizing small subaliquots of the samples, and pools of all samples were created for analytical quality control and evaluation of unknowns. Samples were run under a variant of a standard reverse phase gradient procedure in the repetitive sequence Control Standard, Pool, 7 Samples, Control Standard, Pool, ... as set forth in Table 1: TABLE 1. CHARACTERISTICS OF METHOD 1. Real time set up review; Times and events in the methods chromatographic functions. Review of live method: Potentials and currents of channels 1-16 range functions temperature and limits. -10A- Table 1 Mobile Phase: A: 0.05 M LiPO pH 2.86, 3 mg/L Dodecyl Sulfonic Acid B: 30% MeOH, 0.05 M LiPO 4 , pH 2.86, 3 mg/L Dodecyl Sulfonic Acid Column: Dual 5u ODS, 4.6 ram x 15 cm ESA NBS -11- Samples were analyzed for 38 known components (listed in Table 2) and for 18 well-defined unknown peaks that were isolated in all pools shown in Table 3. (Asterisks denote components used in evaluation of regression and cluster analysis statistical procedures for categorization of groups. ) Table 2 Oracle compatible record showing retention times, digitized characteristic responses across channels, and set of control standard. Abbreviations are described in Table 6. -11A- Table 2 -11B- Table 2 (Cont 'd. ) Table 3. Unknown Pool Peak Matching. Table shows dominant or maximum channel height and retention times of unknown peaks in the pool CSF with default concentrations set at 100. The peaks are referenced to XAN and Tyrosine, concentration 1000; HV.A, concentration 200; Tryptophan, concentration 700. The table shows the match of the seventh pool in the study to the nineteenth pool in the study and this represents typical drift over a two week period. Indicate variables used in regression analysis. -11C Table 3 •Standard: Table: ?CG 1JΛ Study: CΞF1SΛ # COM pound:; identified: 22 ϋ Kit our. CoiπpoiJnd≤ ;iot, Found: « Unknown Peak Clusters: 74ϋ -12- The analysis records were linked by a unique identifier to clinical data of clinical diagnosis, diagnostic criteria, age, pharmaceutical history, sex and race. Pools analyzed as samples against standards for known values were utilized to assess the precision of known compound values in the data base. Standards analyzed against sequential standards were used as a measure of instrumental performance and pools analyzed against sequential pools were utilized as a measure of the precision of unknown peaks. Validation of the Data Control standards analyzed.against sequential control standards yielded precision values ranging from ±l%-±4% CV with no outlying values. Pools analyzed as samples gave precision values ranging from ±2-±7% CV for compounds present at the 0.5 ng/mL level or greater and typically ±25-30% for compounds present at 2X the detection limit of 0.02-0.03 ng/mL (e.g. 5HT, EPI). Pools analyzed against sequential pools for unknowns gave values of ±3-±15% coefficient of variation. Typically, the coefficient of variation of the pools was 5-25 fold less than the coefficient of variation of analytes in a group of samples. Essentially, the contribution of assay variability to the results is minimal. The data base, upon completion, contained 280 samples by 57 analytes (17,000 records). Of these, 163 were null either because no peaks were detected at the sensitivity limits of the assay, or because a signal detected did not meet the qualitative criteria for purity. Regression Analysis Linear regression analysis and stepwise regression analysis were used in a preliminary evaluation of the data. Both raw and mean corrected data was evaluated. Regression comparison of the AD group (61) vs. controls (60) setting AD=1 and C=0 gave a categorical separation regression equation with an S (standard error -13- of estimate) value=0.39 and p=.0041 for 27 of the most significant known compound variables identified in stepwise regression (labeled with asterisks in Table 2). Inclusion of 7 of the most significant variables (labeled with asterisks in Table 3) from the pool analyzed unknown peak data base gave values of S=0.382 and p=0.0037. Assuming a clinical diagnostic error rate in the order of 10%, seven AD samples with regression calculated values (from -1.2 to 0.01) were removed from the calculation. The regression characteristics were then S=0.352 and p=0.0031. Regression of the AD group with AD=1 vs. all others (219)=0 for the same variable group yielded an equation with S=0.481 and p=0.0013. Observations; Although the AD group is separated from mother groups with a high degree of probability, there is too high a degree of overlap for a simple linear regression algorithm to accurately categorize an individual sample. Cluster Analysis Procedures Cluster analysis procedures using nearest neighbor and furthest neighbor approaches were applied to the data base. With both these approaches, the AD group tended to cluster, but controls were scattered relatively evenly, both outside and inside the AD region. Thus, the cluster analysis approach is not suitable as a categorization tool for this type of data. Observations; The behavior of the data under cluster analysis protocols, and the observations that the standard deviations of compound values and of precursor/product ratios across metabolic pathways within a disorder group are smaller than within control groups is consistent with the hypothesis that the biochemical response of controls or normal individuals is more chaotic than that of disordered individuals. -14- Frequency Distribution Probability Analysis The observations on the nature of the data distributions coupled with the technical ability to run large numbers of samples and variables offers an approach to categorization based on differences in the frequency distributions of variables in different disorder categories. This approach relies on basic probability considerations without any assumptions on the shape of a distribution curve or linearity of relationships. The simplest question that I have investigated for the preliminary data base is that given an unknown sample, what is the probability (p) that that sample belongs in one group and not another. For one variable, the question takes the form: where F(V n ) A or f(V n ) B = the frequency with which an unknown sample value (V n ) occurs in category A or category B. For multiple compounds, the expression expands: f( l) A -f( 2 )A---f( n)A P = f(V 1 ) A -f(V) A ...f(V n ) A + f(V 1 ) B -f(V 1 ) B ...f(V n ) B If all frequencies are the same, the P value is 0.5 or a 50/50 chance that the unknown sample is A and not B. A positive answer compresses the expression to a 1 and a negative answer to 0. Like cluster procedures and unlike regression, the use of the algorithm is independent of the number of variables used. Implementation of the Procedure: The implementation of the procedure is by the following steps: 1. Frequency distributions (shown in Figs. la-Id) were created by using a smoothing algorithm based on a 3 -15- point polynomial expansion function that treats each point in the sparse data distribution with equal weighting as the means of a distribution with a width at half height proportional to its value. The use of smoothing functions is a necessary assumption until the n of samples in a particular category reaches approximately 300-400. The procedure used was to divide all data in the categories by the maximum value among categories X 85, apply the polynomial expression algorithm, and normalize the data distributions for the number of samples in each category. The frequency distributions in each category are then organized into look up tables for each variable (Table 4) is then inserted into the look up table. Individual values are divided by the range value X 85 and the frequencies for each variable for category A and B are sequentially calculated in the algorithm after subtracting the effect of that sample from the frequency table. The effect of a sequential calculation across a group of variables listed and described in Table 6 is shown in Figure 2 for 3 AD and 3 C cases from a group of 61 AD and 44 controls. One of the major features of the algorithm is that no single variable predominates as a differentiator among a large group of samples. -16- TABLE 4. Distribution Table. The table shows frequency distribution intervals 1-85 for two of the 71 variables in a look up table initially used in scoring ad and control samples. The tables are obtained by dividing all data by the highest value for the combined data base times 85. The raw distribution is then smoothed by a polynominal expansion which treats each value with equal weighting as the mean of a distribution whose scatter increases with the value. The overall distributions of groups are normalized to the same area. DISTRIBUTION TABLE -17- 2. A sample record such as that shown in Table 5 TABLE 5. DATA CASE RECORD. The data case record is spooled from the Oracle data base. In the example shown, the query requested biochemical data, the diagnosis "control" C and the I.D. "AT0022". The distribute range is created from the record by dividing each value by the maximum range value times 85 to indicate the look up range for the freqency of occurrance of the compound in the AD and C distributions. Table 5 DATACASE RECORD C DISTRIBUT AT0022 RANGE MT3 am 80 BO 1 OHAN3 0.018080 3 QHKY3 0.090402 4 30 D 1.970763 44 HVA 58.03809 56 MHPG 8.714754 62 12 2.911604 22 P01 57.25159 9 P02 135.9827 49 P03 1.808040 1 -18- TABLE 6. DESCRIPTION OF SCORING VARIABLES. The table shows the arbitrary sequency in which 71 variables were applied to the initial scoring of AD vs. C (See Figure 2) . The order of scoring has no effect on the final outcome. Acetaminophen which was included in the assays was not used as a scoring variable. The table also contains the names of the known compounds assayed and the probable moieties in the unknown peaks inferred from the chromatographic/electrochemical data and in vitro studies. Pathway ratios one to twelve are calculated as molar ratios. -18A- Table 6 DESCRIPTION OF SCORING VARIABLES -18B- Table 6 (Cont.) Name or Possible Characteristic Xanthine Methoxyhydroxybenzene - Cysteine or Methionine Di Peptide Catechol - Catechol - Cysteine Condensation or Peptide Cysteine Condensation or Peptide TRP or TYR Peptide TRP or TYR Peptide Indole - Indole - Methoxybenzene - Indole - Hydroxybenzeπe - TYR or TRP Peptide TYR or TRP Peptide Indole - Hydroxyindole - Hydroxyindole - Oxidative of Backwave of HVA Oxidative of Backwave of TRP Oxidative of Backwave of TYR Oxidative of Backwave of XAN TRP/5HIAA + 5HTOL + 5HT + 5HTP TRP/ OHAN + KYN TRP/KYN + HTP + HTOL + 5HIAA + 5HT HTP/5HIAA 5HIAA/5HTOL KYN/OHKY HVA/5HIAA TYR/HPLA TYR/HVA + LD + E + NE + MHPG + DA + 3MT + NMN HVA/MHPG TYR/HPLA + HVA + LD + E + NE + MHPG + DA + 3MT + NMN XAN/HVA -19- Testing the Algorithm on AD vs. Controls For an initial test, conditions were set up such that each individual sample was evaluated as if the data base were set up without its contribution. The results of the initial scoring are shown in Fig. 3. The scoring of five of the 61 AD cases as controls (p=less than 0.01 that the sample is an AD and not a control) is not surprising given the probable diagnostic error rate in AD. The scoring of 4 of the controls as AD are of concern. One possible explanation is that the AD data base is in effect contaminated by five cases that clearly do not match the overall AD profile and are probably not AD. When these five samples are removed from the AD data base and all samples, including the 5 removed, are scored, the control and AD groups are uniquely separated as shown in Fig. 4. The five samples that were removed from the AD scoring data group distribute in an equivocal region from 0.1 to 0.9. In subsequent application of the procedure and algorithm to AD samples vs. all other samples (PD, SC, HD, PSP and C) in the data base asking the question is this sample in the AD distribution and not in the distribution of all others yielded similar results scoring AD samples with p values = 0.98 or greater. The distribution of scores of all others was scattered from 0.001 to .8 including the 5 AD samples which previously scored in this region vs. controls. Testing PD vs. Controls and PD vs. AD Eighteen unknown peaks were initially selected for assay in all samples because of differences among preliminary pools of AD, PC, C and SC. Of these, six occurred only in one sample of the group and 11 had distributions of only slight differentiation capability. One peak, designated P05 shown in Figure 5, occurred predominantly in the PD group. Otherwise, values above 100% of pool values were seen in only six of the AD group. The mean values of P05 and S.D. relative to a pool value -20- of 100 were PD 415 + 200, Control 20 ± 15, AD 30 ± 80, and SC 28 ± 17. This one variable utilized alone in the distributin p scoring algorithm separates PD from controls with PD values from 0.4 - 0.999 and C values from 0.51 - 0.001. With all variables, the PD values were all greater than 0.994 and C values less than 0.003. Scoring AD and not PD without P05 separated 58 PD samples with scores of less than 0.02 and two with scores of 0.14 and 0.18. 55 AD samples scored above 0.98 with six in the region between 0.25 and 0.76, including the 5 cases which initially scored as controls. Including P05 in the scoring, all PD samples scored below 0.001. 57 of the AD samples scored above 0.98 and 4 of the initial 5 that matched the control group scored near the PD samples from 0.11 to 0.32. Example II Pathological Changes in the Olfactory System in Alzheimer's Disease In Alzheimer's Disease a number of studies have shown that the structures of the limbic system, including entorhinal cortex, hippocampal formation, basal forebrain and the amygdala are the most severely affected areas of the brain. Since the amygdala, entorhinal cortex and uncal hippocampus are strongly related to olfactory input (Pearson et al, 1985) - suggested that the olfactory pathways may be the initial site of pathology in the disease. In fact some authors have suggested that inhaled molecules, such as aluminosilicates, could contribute to the etiology of the disease (Roberts, 1987). - 1. Pearson R.C.A., Esiri M.M., Hiorns R.W., Wilcock G.K. Powell T.P.S. Anatomical correlates of the distribution of the pathological changes in the neocortex in Alzheimer's disease. Proc. Natl. Acad. Sci. 82: 4531-4524, 1985 2. Roberts, E. Alzheimer's disease may begin in the nose and may be caused by aluminosilicates. Neurobiol. Aging. 7:561-567, 1986; Shipley M.T. Transport of Footnote Cont. Next Page -21- The olfactory nerve cells in the olfactory epithelium project through the cribiform plate to the olfactory bulb (Kosel et al, 1981 ) - Primary olfactory fibers synapse in olfactory glomeruli with descending dendrites of mitral and tufted cells ; which are the primary output neurons of the bulb. Axons of mitral and tufted cells enter the olfactory tract and provide input to the anterior olfactory nucleus, as well as to the central projections of the olfactory system. The anterior olfactory nucleus, located in the center of the bulb, gives rise to a recurrent collateral to the bulb and to a crossed projection to the anterior commissure. The olfactory tract passes through the anterior perforated substance and projects to prepiriform cortex, corticomedial nuclei of the amygdala, entorhinal and perirhinal cortices, inferior surface of the frontal lobe, insula, temporal pole and basal forebrain (Haberty et al, 1977). 4 Several lines of evidence have shown that the olfactory system is affected in ALzheimer's Disease. Several groups have shown that Alzheimer's Disease patients show deficits in olfactory recognition (Doty et al, 1987). 5 A study which examined patients with very Footnote 2 Cont. molecules from nose to brain; transneuronal anterograde and retrograde labeling in the rat olfactory system by wheat germ agglutinin-horseradish peroxidase applied to the nasal epithelium. Brain Res. Bull. 15: 129-142, 1985; Peri D.P., Good P.F. Uptake of aluminum into central nervous system along nasal-olfactory pathways. Lancet 1: 1028, 1987. 3. Kosel K.C., Van Hoesen G.W., West J.R. Olfactory bulb projections to the parahippocampal area of the rat. J. Comp. Neurol. 198: 467-482, 1981. 4. Haberty L.B., Price J.L. The axonal projection patterns of the mitral and tufted cells in the olfactory bulb in the rat. Brain Res. 129: 152-157, 1977. 5. Doty R.L., Reys P.F., Gregor T. Presence of both odor identification and detection deficits in Alzheimer's Footnote Cont. Next Page -22- mild Alzheimer's Disease showed impairment on a task of identification of odors, however olfactory thresholds were normal (Koss et al, 1988). - Pathologic studies have recently documented changes in the primary olfactory receptor neurons in Alzheimer's Disease. Talamo et al, 1976^ reported nasal epithelium from Alzheimer's Disease patients showed increased immunoreactivity for phosphorylated neurofilaments, as well as Tau and Alz-50 positive neurites. These neurons however did not contain neurofibrillary tangles. Neurite formation in the olfactory epithelium did correlate with numbers of neurofibrillary tangles and senile plaques in the hippocampus of the Alzheimer's Disease brains. The possibility of using biopsy of nasal epithelium as a clinical marker of Alzheimer's Disease during life was suggested. Pathological studies of the olfactory bulb in Alzheimer's Disease have shown consistent changes. Esiri and Wilcock, (1984) 8 found neurofibrillary tangles in the olfactory bulb in tufted cells, outer granule cells, and Footnote 5 Cont. Disease. Brain Res. Bull. 18: 597-600, 1987; Green J.E., Songsanand P., Corkin S., Growdon J.H. Olfactory capacities in Alzheimer's Disease. Neurology 39, Suppl 1: 138, 1989; Rezek D.L. Olfactory deficits as a neurologic sign on dementia of the Alzheimer type. Arch. Neurol. 44: 1030-1032, 1987. 6. Koss E., Weiffenbach J.M. , Haxby J.V. , ^riedland R.P. Olfactory detection and identification performance are dissociated in early Alzheimer's Disease. Neurology 38: 1228-1232, 1988. 7. Talamo B.R., Rudel R.A. , Kosik K.S., Lee V.Y.M., Neff S., Adelman L. , Kauer J.S. Pathological changes in olfactory neurons inn patients with Alzheimer's Disease. Aging 7: 11-14, 1976. 8. Erisi M.M. , Wilcock G.K. The olfactory bulbs in Alzheimer's Disease. J. Neurol. Neurosurg. Psychiatr. 47: 56-60, 1984. -23- in the anterior olfactory nucleus. Hyman and collegues (1991)9 h as examined 10 control and 10 Alzheimer's Disease olfactory bulbs. Large numbers of neurofibrillary tangles were consistently found in the Alzheimer's Disease anterior olfactory nucleus, with only small numbers in the mitral and tufted cells. This is consistent with other primary sensory systems in Alzheimer's Disease in which higher order association areas show more severe pathology than primary sensory areas. The olfactory epithelial neurons and olfactory bulb have been shown to contain very high levels of carnosin (β-alanyl-L-histidine) , and several studies have shown that carnosine is released in the olfactory bulb in response to peripheral inputs ( acrides et al, 1983 ) . -® The external tufted cells in the olfactory bulb contain dopamine, substance P or both. The deeper tufted cells may use excitatory amino acids. Periglomerular cells contain GABA, enkephalin or dopamine. The olfactory bulb receives a very strong centrifugal projection including enkephalin, vasoactive intestinal polypeptide, LHRH, somatostatin and substance P fibers. A cholinergic projection originates mainly in the ventral nucleus of the diagonal band. Serotonergic projections originate from the dorsal raphe nuclei, while there is a noradrenergic projection from the locus ceruleus. The olfactory bulbs therefore receive a rich and variegated neurochemical input, including projections from the cholinergic, 9. Hyman B.T. , Arriagada P.V., Van Hoesen G.W. Pathological changes in the olfactory system in aging and Alzheimer's disease. Int. Study Group for the Pharmacology of Memory Disorders, in press. (1991) 10. Macrides F. , Davis B.L. The olfactory bulb. In: Emson P.C. (ed.) Chemical Neuroanatomy, Raven Press, New York, pp. 391-426, 1983. -24- noradrenergic and serotonergic nuclei which are known to be affected in Alzheimer's Disease (Macrides et al, 1983 ) . -- Several factors relating to the suitability of using nasal secretions for neurochemical analysis have been examined. These include sample acquisition, potential effects of bacterial or viral infection, the specific compounds which can be identified, and means for normalizing the data to compensate for sampling variability and sample size. Following the electrochemical analysis procedure above described for CFS all the compounds listed in Table 2 have been identified using nasal secretion samples although the concentration profiles are significantly different than for CSF. (See typical patterns in Figures 6A and 6B.) Notably, dihydroxyphenylacetic acid (DO AC) is present in greater concentrations than homovanillic acid (HVA) and serotonin (5HT) is relatively equivalent to 5-hydroxindoleacetic acid (5HIAA). Also notably different is the complexity of the region of small di-and tri- peptides with retention times of 30-45 min on channels 10-15. Essentially, the nasal mucosa appear to contain representative compounds from all metabolic pathways that occur in the neuronal projections to the olfactory bulb. Two characteristics make nasal secretions attractive for neurochemical analysis. The secretions are close to the site or origin of the compounds being measured and they are in a strong reducing environment. Oxidation potential measurements of saline suspensions of nasal secretions with platinum vs. Ag/Ag/Cl yield values from -0.230 to -0.300 mV compared to lumbar CSF values of -0.050 to -0.170. The reducing character of the tissue is consistent with observation of a large number of easily 11.Macrides F. , Davis B.L. The olfactory bulb. In: Emson P.C. (ed.) Chemical Neuroanatomy, Raven Press, New York, pp. 391-426, 1983 -25- oxidized peaks on the lower channels of the chromatographic pattern. The reducing character of the sample makes it highly attractive for clinical use because of the stability inferred for the target compounds. Indeed preliminary studies indicate no change in 5-hydroxyindoleacetic acid or homovanillic acid over 24 hours in perchloric acid extracts at room temperature. From the number of peaks in the 35-50 min region, the nasal mucosa also appear to be a stable matrix for small peptides, including fragments of 3-amyloid peptide. This hypothesis was tested by incubating β-amyloid in saline suspensions of nasal mucosa and observing the development over 6 hours of eleσtroactive peaks, five of which corresponded to original peaks in the suspension as shown in Figure 7. Sample Acquisition Preparation and Normalization Samples typically weigh approximately 100 mg. A preliminary evaluation of swabs taken from high in the nostril, approximately 2 cm, vs. a wipe of the nares at approximately 0.5 cm, indicate that all species of interest are present lower in the nose, but at significantly reduced concentrations overall. In an initial study from an n of 4 there were no significant changes in the ratio of DOPAC/HVA or in tryptophan (TRP) to kynurenine (KYN) , but the ratio of 5HIAA/5HT increased, consistent with reduced concentrations of 5HT in swabs obtained lower in the nose. Similarly, on an n of 4 , a relatively aggressive scrubbing the swab yields a pattern with higher concentrations than a more gentle scrub, but without significant changes in overall pattern. Swabs were extracted by cutting the tip and placing it in 300 1 of 0.1 M HCIO4 , they ere vortexed for 1 min. centrifuged and then reextracted in another 300μ 1 of 0.1 M HCIO . Further sequential extractions did not result in any improvement in recoveries. Following centrifugation the pellets are combined, and can be extracted with 40:60 -26- acetonitrilehexane for analysis of ubiquinones and larger peptides. Several brands of cotton swab were evaluated for blank effects. Of these, plastic handled Johnson and Johnson (TM) swabs proved clean enough to be used as received. Paper composite and wooden handles of various types showed a number of small peaks that have the potential to interfere with the target compounds of interest. The selection of appropriate divisors or normalizers of the data is a major consideration since sample size is not controllable. This issue was investigated in a preliminary study of 6 individuals from whose left and right nostril swabs were obtained. From the initial data analysis, xanthine (XAN), uric acid (UA), tyrosine (TYR), DOPAC, HVA, TRP, KYN, 5HIAA and 5HT have been measured. The 6 left and right individual samples as duplicate pairs coefficient of variation were from ±60 - ±80%. When all data was divided by xanthine, the variations were ±15 - ±30%. When precursor product ratios were used, the variations were further lowered (DOPAC/HVA mean 1.4 ± 7%; TRP/KYN mean 2.6 ± 14%; and 5HT/5HIAA mean l.l± 13%) Possible Effects of Bacterial and Viral Infections Two available cultures of Pseudomonas and Staphlococcus were evaluated to look for possible interferences from bacterial infection. Aliquots of approximately 10 mg wet weight of cultured bacteria isolated by centrifugation and washing were sonified in 200μL of 0.1 M HCIO4, centrifuged and the supernatant analyzed. Responses for XAN, KYN, TYR, TRP, guanosine and hypoxanthine were observed of approximately 50-100% of the response of a typical nasal swab. No dopaminergic, serotonergic, or noradrenergic metabolites were observed. Consequently, it is thought that bacterial contamination would have to be quite massive, constituting over 10% of a -27- typical sample to have a 10% or greater effect on the measurements. In the initial study three of the six subjects had colds and there were no significant differences between them and unaffected individuals. There was a slight, but not statistically significant increase in KYN/XAN and a decrease in TRP/KYN ratios. It thus appears that nasal secretions may advantageously be employed as samples for neurochemical analysis using electrochemical detection techniques. Similar results were found using platelets. Moreover, the invention advantageously may be employed for diagnosing a disease condition at an early stage, i.e. before observable physical manifestations. For example, in the case of Alzheimer's Disease, the exact etiology is unknown. However, there is strong evidence that genetic factors play a role (St. George-Hyslop et al, 1990). ^~ The means by which a genetic defect contributes to the pathologic features of the illness is unclear. A major feature of the pathology is the accumulation of the 3-amyloid protein in senile plaques, blood vessels, skin By correlating changes in neurochemical markers with changes in accumulation of the β -amyloide protein as the Alzheimer's Disease progresses, it is possible to provide an early diagnosis for Alzheimer's Disease. Also, identification of neurochemical markers for Alzheimer's Disease may provide a basis for prevention and/or treatment, i.e. by identifying precursors, progression of the disease may be slowed, stopped or even reversed. The invention has been described for use in diagnosing Alzheimer's Disease. It will be understood, however, that the invention advantageously may be used to diagnose and 12. St. George-Hyslop P.H. , Haines , J.L., Farrer L.A., et al. Genetic linkage studies suggest that Alzheimer's disease is not a single homogeneous disorder. Nature 347: 194-197, 1990. -28- characterize other neurological, degenerative or defective disorders such as Huntington's Disease, Parkinson's Disease, schizophrenia, progressive supernuclear palsy, amyotrophic lateral sclerosis and senile dementias. The invention also advantageously may be used to classify and diagnose tumors, carcinomas, cardiovascular abnormalities and other disorders, or for selection of therapy based on categories of known successful vs. unsuccessful outcomes. Moreover, both treatment protocols and new pharmaceuticals may be evaluated. Still other changes and advantages will be obvious to one skilled in the art. ■ 13. Joachim C.L. , Morris J.H., Selkoe D.J. Clinically diagnosed ALzheimer's disease: autopsy results in 150 cases. Ann. Neurol. 24: 50-56, 1988.

Claims: -29- CLAIMS 1. In a method for diagnosing disorders in which b iological samples from normal and afflicted (abnormal) individuals are analyzed to generate patterns representative of molecular constituents of said samples, the improvement which comprises creating a data base of frequency distribution patterns of constituents of samples from organisms having known categories of disorders and controls, and comparing the sample analysis for conformity to said frequency distribution patterns. 2. A method according to claim 1, wherein said samples comprise body fluid. 3. A method according to claim 2, wherein said body fluid is selected from cerebrospinal fluid, plasma, platelets, nasal mucosa, serum, saliva, and urine. 4. A method according to claim 1 or 2 or 3, wherein said samples comprise electrochemically active compounds, and including the steps of passing said fluid samples sequentially through a liquid chromatographic column for achieving time-space separation of the materials eluting from the column and an electrochemical detection apparatus whereby to generate electrochemical patterns of said electrochemically active compounds. 5. A method according to claim 1 or 2 or 3 or 4, wherein said disorder is selected from Alzheimer's Disease, Parkinson's Disease, Huntington's Disease, schizophrenia, Progressive Supernuclear Palsy, amyotrophic lateral sclerosis, senile dementis, tumors, carcinomas, and cardiovascular abnormalities. 6. A method according to claim 4, and including the step of separating said electrochemically active compounds by electrochemical characteristics in said electrochemical detection apparatus. 7. A method for diagnosing disorders in which biological samples from normal and afflicted (abnormal) individuals are analyzed to generate patterns -30- representative of selective parameters, the improvement which comprises examining said patterns for chaotic or non-linear values. 8. A method according to claim 7, wherein the parameters measured comprise biochemical patterns of small molecules. 9. A method according to claim 8, wherein said parameters are representative of at least one of the following selected from the group consisting of neurotransmitters, cofactors, precursors, metabolites, and associated compounds. 10. A method according to claim 7, wherein said disorder comprises Alzheimer's Disease, and the parameters measured comprise tyrosine and t'ryptophan peptide degradation fragments from beta amyloid. 11. A method according to claim 7, wherein said disorder comprises Parkinson's Disease, and the parameter measured comprises the amylyte designated P05.

Inventors: MATSON WAYNE R, MATSON, WAYNE, R.
Applicants: ESA INC, ESA, INC.
Application year: 1992
Main classification: G01N30
IPC codes: G01N27, G01N30, G01N33