Summary

Identificazione delle malattie legate Patterns covarianza spaziale utilizzando i dati Neuroimmagine

Published: June 26, 2013
doi:

Summary

Tecniche multivariate tra cui l'analisi delle componenti principali (PCA) sono stati utilizzati per identificare i modelli delle firme di cambiamento regionale in immagini funzionali del cervello. Abbiamo sviluppato un algoritmo per identificare i biomarcatori di rete riproducibili per la diagnosi di patologie neurodegenerative, la valutazione della progressione della malattia, e la valutazione obiettiva degli effetti del trattamento in popolazioni di pazienti.

Abstract

Il modello in scala sottoprofilo (SSM) 1-4 è un algoritmo basato multivariata PCA che identifica le principali fonti di variazione paziente e gruppo di controllo dati immagine cerebrali e respingendo componenti minori (Figura 1). Applicato direttamente ai dati di covarianza voxel-by-voxel di immagini multimodali di stato stazionario, un intero set di immagini di gruppo può essere ridotto a un paio di significativi modelli di covarianza linearmente indipendenti e corrispondenti punteggi soggetto. Ogni modello, definito un gruppo sottoprofilo invariante (GIS), è un componente principale ortogonale che rappresenta una rete spazialmente distribuito di regioni cerebrali correlate funzionalmente. Grandi effetti globali scalari medi che possono oscurare più piccoli di rete specifici contributi vengono rimossi dalla conversione logaritmica intrinseca e significano centratura dei dati 2,5,6. Soggetti esprimono ciascuno di questi modelli in misura variabile rappresentata da un semplice punteggio scalare che può correlare con clin indipendentedescrittori di iCal o psicometrica 7,8. Utilizzando l'analisi di regressione logistica dei punteggi soggetti (ossia i valori di espressione del modello), coefficienti lineari possono essere derivate per combinare più componenti principali in singole malattie legate modelli di covarianza spaziale, cioè le reti composite con una migliore discriminazione dei pazienti provenienti da soggetti sani di controllo 5,6. Convalida incrociata all'interno del set di derivazione può essere eseguita utilizzando tecniche di ricampionamento bootstrap 9. Convalida in avanti è facilmente confermata dalla valutazione diretta punteggio dei modelli derivati ​​a potenziali dataset 10. Una volta convalidati, modelli correlati alla malattia possono essere utilizzati per segnare singoli pazienti rispetto ad un campione di riferimento fisso, spesso l'insieme di soggetti sani che è stato utilizzato (con il gruppo malattia) nel modello originale derivazione 11. Questi valori standardizzati possono a loro volta essere usati per aiutare nella diagnosi differenziale 12,13 e per valutare la malattiaprogressione e il trattamento degli effetti a livello di rete 7,14-16. Vi presentiamo un esempio di applicazione di questa metodologia di FDG PET dati dei pazienti con malattia di Parkinson e dei controlli normali utilizzando il nostro software in-house per ricavare una caratteristica covarianza modello biomarker di malattia.

Introduction

Malattie neurodegenerative sono stati ampiamente studiati con tecniche che localizzano e quantificare le anomalie del metabolismo del cervello così come i metodi non-inferenziali che studiano le interazioni regionali 17. Data-driven strategie di analisi multivariata come l'analisi delle componenti principali (PCA) 1,2,4,18 e analisi delle componenti indipendenti (ICA) 19,20, così come le tecniche di supervisione come quadrati parziali (PLS almeno) 21 e tendenze ordinali canoniche variates analisi (Ort / CVA) 22 può rivelare modelli caratteristici o "reti" di attività correlate. Le basi di procedure multivariate, in particolare il modello in scala sottoprofilo (SSM) 1,2,4-6,18 sono stati precedentemente descritti in JoVE 3. Questo approccio PCA-based è stato originariamente sviluppato per esaminare anormali rapporti covarianza funzionali tra le regioni del cervello in stato stazionario di singole immagini di volume di flusso sanguigno cerebrale e il metabolismo acquired nello stato di riposo di modalità come la PET e la SPECT che presentano caratteristiche di alto rapporto segnale-rumore. Modelli SSM malattie specifiche sono biomarcatori che riflettono differenze complessive nella topografia regionale nei pazienti rispetto ai soggetti normali 7,16 e può riflettere un singolo processo di rete o l'assimilazione di diverse complesse funzioni anormali 23. Reti del cervello modello di covarianza metabolici sono associati a valori di espressione (soggetto punteggi) in grado di distinguere tra i gruppi normali e di malattia e di fornire misure di network-based che correlano con le valutazioni cliniche di gravità della malattia. In genere, soggetti punteggi per tali modelli aumentano con la progressione della malattia e possono anche essere espressi prima della comparsa dei sintomi 14,24. Anzi, correlati alla malattia biomarker rete sono state caratterizzate per patologie neurodegenerative come il morbo di Parkinson 10 (PD), la malattia di Huntington 25 (HD) e malattia di Alzheimer 8 </sup> (AD). È importante sottolineare che le malattie legate topografie metabolici sono stati individuati per atipiche disturbi motori parkinsoniani, quali atrofia multisistemica (MSA) e la paralisi sopranucleare progressiva (PSP). Questi modelli sono stati utilizzati in concerto per la diagnosi differenziale di individui con sindromi clinicamente simili "look-alike" 12,13,26.

Al contrario, i metodi univariati voxel-based fMRI tipici valutare la significatività delle differenze tra pazienti e controlli in ammassi cerebrali isolate. Più recentemente, sono stati sviluppati metodi per misurare connettività funzionale tra le regioni del cervello variamente definiti 27-29. Questa definizione di connettività funzionale è limitato a soggetto e regione interazioni e devia specifici dal concetto originale SSM / PCA che si riferisce alla interconnettività sezione trasversale di intrinseche spazialmente distribuiti cerebrali aree di rete 1,2,23,30. A loro vantaggio, piattaforme MRI unre facilmente installato, ampiamente disponibile, non invasivo e in genere richiedono tempi di scansione più breve di modalità di imaging radiotracciante tradizionali come il PET o SPECT con un conseguente aumento di potenziali metodologie descritte nella letteratura recente. Tuttavia, i segnali fMRI tempo-dipendenti risultanti forniscono misure indirette di attività neurale locale 31,32. Gli algoritmi analitici generalmente complessi impiegati sono stati limitati dalle grandi dimensioni di set di dati, il rumore fisiologico insito nei segnali fMRI, così come l'elevata variabilità dell'attività del cervello che esiste tra temi e regioni 19,23. Anche se le informazioni interessanti per quanto riguarda l'organizzazione del cervello può essere dedotta dalle proprietà delle "reti" fMRI, non sono stati sufficientemente stabile per essere utilizzato come biomarker di malattia affidabili. Inoltre, le topografie risultanti rete non sono necessariamente equivalenti a quelli individuati utilizzando metodologie di imaging funzionale affermati come SSM / PCA. Per The la maggior parte, rigorosa convalida incrociata delle topografie risultanti fMRI è stata carente, con alcuni esempi di applicazione di successo avanti di modelli derivati ​​nei dati di scansione prospettici da casi singoli.

Un vantaggio di covarianza analisi PCA consiste nella sua capacità di identificare le principali fonti di variazione dati nei primi componenti principali ma è inefficace se gli autovettori rappresentano importanti fattori di disturbo casuale anziché risposta effettiva rete intrinseca. Selezionando solo le prime autovettori e limitando a quelli che mostrano differenze significative nei punteggi di controllo dei pazienti contro il normale, abbiamo notevolmente ridurre l'influenza di elementi di disturbo. Tuttavia, per l'approccio di base qui descritto, tali misure possono non essere sufficienti a generare stimatori robusti in un tipico insieme di dati fMRI ad eccezione delle modalità descritte di seguito.

Così, a causa della relazione diretta stabile di g regionalelucose metabolismo e l'attività sinaptica 33, questa metodologia è stata applicata principalmente per l'analisi dei dati dello stato FDG PET riposo. Tuttavia, dato che il flusso ematico cerebrale (CBF) è strettamente accoppiato ad attività metabolica in stato di riposo 10,11,34, SPECT 35,36 e più recentemente arteriosa spin labeling (ASL) metodi di imaging di perfusione MRI 37.38, sono stati utilizzati per valutare l'attività metabolica anomala nei singoli casi. Detto questo, la derivazione di affidabili modelli di covarianza spaziale con fMRI stato di riposo (rsfMRI) è come in precedenza non rilevato semplice 31,32. Anche così, SSM / PCA analisi preliminare dei dati rsfMRI da PD pazienti e soggetti di controllo ha rivelato alcune omologie topografiche tra i modelli di relazione con la malattia identificati utilizzando le due modalità, PET e ampiezza delle oscillazioni a bassa frequenza (ALFF) di BOLD fMRI 39,40 . Infine, si segnala, inoltre, che questo approccio è stato applicato con successo in voxel base morfometria (VBM) dati MRI strutturali 41,42, rivelando particolari modelli di covarianza spaziale associati alla perdita di volume legata all'età e in ulteriori confronti di VBM e modelli ASL nelle stesse materie 43. Il rapporto tra SSM / PCA topografie di covarianza spaziale e le reti cerebrali analoghe identificate utilizzando diversi approcci analitici e piattaforme di imaging è un argomento di indagine in corso.

Protocol

1. Raccolta dati e pre-elaborazione Il metodo di SSM / PCA può essere applicata a immagini di volume singole ottenuti da varie fonti e modalità. In particolare, per la on-site di imaging PET del metabolismo, preparare un adeguato tracciante radionuclide ad esempio [18 F]-fluorodeossiglucosio (FDG) e somministrare a ciascun paziente. I pazienti sono di solito sottoposti a scansione a riposo con gli occhi aperti, a seguito di un digiuno di almeno 12 ore, fuori farmaci. Scansione di ogni o…

Representative Results

Una semplice applicazione di analisi multivariata SSM / PCA per derivare un modello biomarcatore neuroimaging per PD di seguito è illustrato. PET FDG immagini di dieci clinicamente diagnosticati pazienti PD (6M/4F, 59y 7y ± sd) di durata variabile malato (9 anni ± 5A sd) e dieci età e sesso abbinati controlli normali (6M/4F, 58y 7y ± sd) sono stati analizzati utilizzando nostra routine ssmpca. I venti corrispondenti spazialmente immagini pre-normalizzati sono stati inizialmente selezionati nelle categorie …

Discussion

Il modello SSM / PCA originariamente presentato da Moeller et al. 4 è evoluta 1-3 in una tecnica semplice e robusto per l'analisi dei dati di neuroimaging. Tuttavia, ci sono stati ambiguità nell'applicazione di questa metodologia che abbiamo tentato di chiarire qui e in pubblicazioni precedenti 5-7,10. Alcuni di questi problemi sono stati affrontati nel testo, ma sono reemphasized qui per la loro importanza. Come specificato nell'introduzione, SSM / PCA è soprattu…

Disclosures

The authors have nothing to disclose.

Acknowledgements

Questo lavoro è stato sostenuto da Grant No. P50NS071675 (Morris K. Udall Centro di Eccellenza nella ricerca sulle malattie di Parkinson presso il Feinstein Institute for Medical Research) di DE presso l'Istituto nazionale dei disordini neurologici e colpo. Il contenuto è di esclusiva responsabilità degli autori e non rappresentano necessariamente il punto di vista ufficiale del National Institute of Neurological Disorder and Stroke o il National Institutes of Health. Lo sponsor non ha giocato un ruolo nel disegno dello studio, la raccolta, l'analisi e l'interpretazione dei dati, la scrittura del report o nella decisione di presentare la carta per la pubblicazione.

Materials

Name of Equipment Company Catalog Number Comments
Image Acquisition
PET Scanner GE Medical Systems GE Advance Any PET, PET/CT and PET/MRI Scanners from GE, Siemens and Philips
PC Workstations Lenovo Any http://www.lenovo.com/us/en/
Radiopharmaceuticals
[18F]-fluorodeoxyglucose Feinstein Institute for Medical Research Routine Production Also distributed by Cardinal Health http://www.cardinal.com/
Software
ScanVP Feinstein Institute for Medical Research Version 5.9.1, Version 6.2, To be released www.feinsteinneuroscience.org
SPM The UCL Institute of Neurology spm99-spm8 http://www.fil.ion.ucl.ac.uk/spm
Windows Microsoft Any
Matlab Mathworks Matlab Version 7.0, 7.3 http://www.mathworks.com/
JMP SAS Version 5 http://www.jmp.com/

References

  1. Moeller, J. R., Strother, S. C. A regional covariance approach to the analysis of functional patterns in positron emission tomographic data. J. Cereb. Blood Flow Metab. 11 (2), 121-135 (1991).
  2. Alexander, G. E., Moeller, J. R. Application of the scaled subprofile model to functional imaging in neuropsychiatric disorders: A principal component approach to modeling brain function in disease. Hum. Brain Mapp. 2, 1-16 (1994).
  3. Habeck, C. G. Basics of multivariate analysis in neuroimaging data. J. Vis. Exp. (41), e1988 (2010).
  4. Moeller, J. R., Strother, S. C., Sidtis, J. J., Rottenberg, D. A. Scaled subprofile model: a statistical approach to the analysis of functional patterns in positron emission tomographic data. J. Cereb. Blood Flow Metab. 7 (5), 649-658 (1987).
  5. Spetsieris, P. G., Eidelberg, D. Scaled subprofile modeling of resting state imaging data in Parkinson’s disease: methodological issues. Neuroimage. 54 (4), 2899-2914 (2011).
  6. Dhawan, V., Tang, C. C., Ma, Y., Spetsieris, P., Eidelberg, D. Abnormal network topographies and changes in global activity: Absence of a causal relationship. Neuroimage. 63 (4), 1827-1832 (2012).
  7. Eidelberg, D. Metabolic brain networks in neurodegenerative disorders: a functional imaging approach. Trends Neurosci. 32 (10), 548-557 (2009).
  8. Habeck, C., Foster, N. L., et al. Multivariate and univariate neuroimaging biomarkers of Alzheimer’s disease. Neuroimage. 40 (4), 1503-1515 (2008).
  9. Efron, B., Tibshirani, R. . An introduction to the bootstrap. , (1994).
  10. Ma, Y., Tang, C., Spetsieris, P., Dhawan, V., Eidelberg, D. Abnormal metabolic network activity in Parkinson’s disease: test-retest reproducibility. J. Cereb. Blood Flow & Metab. 27 (3), 597-605 (2007).
  11. Ma, Y., Eidelberg, D. Functional imaging of cerebral blood flow and glucose metabolism in Parkinson’s disease and Huntington’s disease. Mol. Imaging Biol. 9 (4), 223-233 (2007).
  12. Tang, C. C., Poston, K. L., et al. Differential diagnosis of parkinsonism: a metabolic imaging study using pattern analysis. Lancet Neurol. 9 (2), 149-158 (2010).
  13. Spetsieris, P. G., Ma, Y., Dhawan, V., Eidelberg, D. Differential diagnosis of parkinsonian syndromes using PCA-based functional imaging features. Neuroimage. 45 (4), 1241-1252 (2009).
  14. Tang, C. C., Poston, K. L., Dhawan, V., Eidelberg, D. Abnormalities in metabolic network activity precede the onset of motor symptoms in Parkinson’s disease. J. Neurosci. 30 (3), 1049-1056 (2010).
  15. Mure, H., Hirano, S., et al. Parkinson’s disease tremor-related metabolic network: characterization, progression, and treatment effects. Neuroimage. 54 (2), 1244-1253 (2011).
  16. Niethammer, M., Eidelberg, D. Metabolic brain networks in translational neurology: concepts and applications. Ann. Neurol. , (2012).
  17. Petersson, K. M., Nichols, T. E., Poline, J. B. Statistical limitations in functional neuroimaging. I. Non-inferential methods and statistical models. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 354 (1387), 1239-1260 (1999).
  18. Habeck, C., Stern, Y. Multivariate data analysis for neuroimaging data: overview and application to Alzheimer’s disease. Cell Biochem. Biophys. 58 (2), 53-67 (2010).
  19. McKeown, M. J., Hansen, L. K., Sejnowsk, T. J. Independent component analysis of functional MRI: what is signal and what is noise. Current Opinion in Neurobiology. 13 (5), 620-629 (2003).
  20. Stone, J. V. Independent component analysis: an introduction. Trends Cogn. Sci. 6 (2), 59-64 (2002).
  21. McIntosh, A. R., Bookstein, F. L., Haxby, J. V., Grady, C. L. Spatial pattern analysis of functional brain images using partial least squares. Neuroimage. 3 (3 Pt. 1), 143-157 (1996).
  22. Habeck, C., Krakauer, J. W., et al. A new approach to spatial covariance modeling of functional brain imaging data: ordinal trend analysis. Neural Comput. 17 (7), 1602-1645 (2005).
  23. Habeck, C., Moeller, J. R. Intrinsic functional-connectivity networks for diagnosis: just beautiful pictures. Brain Connect. 1 (2), 99-103 (2011).
  24. Huang, C., Tang, C., et al. Changes in network activity with the progression of Parkinson’s disease. Brain. 130, 1834-1846 (2007).
  25. Feigin, A., Tang, C., et al. Thalamic metabolism and symptom onset in preclinical Huntington’s disease. Brain. 130, 2858-2867 (2007).
  26. Poston, K. L., Tang, C. C., et al. Network correlates of disease severity in multiple system atrophy. Neurology. 78 (16), 1237-1244 (2012).
  27. Biswal, B., Yetkin, F. Z., Haughton, V. M., Hyde, J. S. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magnetic Resonance in Medicine: Official Journal of the Society of Magnetic Resonance in Medicine / Society of Magnetic Resonance in Medicine. 34 (4), 537-541 (1995).
  28. Greicius, M. D., Krasnow, B., Reiss, A. L., Menon, V. Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proceedings of the National Academy of Sciences of the United States of America. 100 (1), 253-258 (2003).
  29. Friston, K. J., Frith, C. D., Liddle, P. F., Frackowiak, R. S. Functional connectivity: the principal component analysis of large (PET) data sets. J. Cereb. Blood Flow Metab. 13, 5-14 (1993).
  30. Strother, S. C., Kanno, I., Rottenberg, D. A. Commentary and opinion: I. Principal component analysis, variance partitioning, and “functional connectivity”. Journal of Cerebral Blood Flow and Metabolism: Official Journal of the International Society of Cerebral Blood Flow and Metabolism. 15 (3), 353-360 (1995).
  31. Ekstrom, A. How and when the fMRI BOLD signal relates to underlying neural activity: the danger in dissociation. Brain Research Reviews. 62 (2), 233-244 (2010).
  32. Logothetis, N. K. What we can do and what we cannot do with fMRI. Nature. 453 (7197), 869-878 (2008).
  33. Lin, T. P., Carbon, M., et al. Metabolic correlates of subthalamic nucleus activity in Parkinson’s disease. Brain. 131 (Pt. 5), 1373-1380 (2008).
  34. Hirano, S., Asanuma, K., et al. Dissociation of metabolic and neurovascular responses to levodopa in the treatment of Parkinson’s disease. J. Neurosci. 28 (16), 4201-4209 (2008).
  35. Feigin, A., Antonini, A., et al. Tc-99m ethylene cysteinate dimer SPECT in the differential diagnosis of parkinsonism. Mov. Disord. 17 (6), 1265-1270 (2002).
  36. Eckert, T., Van Laere, K., et al. Quantification of Parkinson’s disease-related network expression with ECD SPECT. Eur. J. Nucl. Med. Mol. Imaging. 34 (4), 496-501 (2007).
  37. Ma, Y., Huang, C., et al. Parkinson’s disease spatial covariance pattern: noninvasive quantification with perfusion MRI. J. Cereb. Blood Flow Metab. 30 (3), 505-509 (2010).
  38. Melzer, T. R., Watts, R., et al. Arterial spin labelling reveals an abnormal cerebral perfusion pattern in Parkinson’s disease. Brain. 134 (Pt. 3), 845-855 (2011).
  39. Skidmore, F., Spetsieris, P., et al. Diagnosis of Parkinson’s disease using resting state fMRI. , LB22 (2011).
  40. Peng, S., Wu, T., et al. A comparison study of Parkinson’s disease-related patterns between FDG PET and fMRI at rest state. Neuroimage. 61, 5610 (2012).
  41. Brickman, A. M., Habeck, C., Zarahn, E., Flynn, J., Stern, Y. Structural MRI covariance patterns associated with normal aging and neuropsychological functioning. Neurobiol. Aging. 28 (2), 284-295 (2007).
  42. Bergfield, K. L., Hanson, K. D., et al. Age-related networks of regional covariance in MRI gray matter: reproducible multivariate patterns in healthy aging. Neuroimage. 49 (2), 1750-1759 (2010).
  43. Steffener, J., Brickman, A. M., Habeck, C. G., Salthouse, T. A. Cerebral blood flow and gray matter volume covariance patterns of cognition in aging. Human Brain Mapping. , (2012).
  44. Cangelosi, R., Goriely, A. Component retention in principal component analysis with application to cDNA microarray data. Biol. Direct. 2, 2 (2007).
  45. Akaike, H. A new look at the statistical model identification. IEEE Transactions on Automatic Control. 19 (6), 716-723 (1974).
  46. Eidelberg, D., Moeller, J. R., et al. Assessment of disease severity in parkinsonism with fluorine-18-fluorodeoxyglucose and PET. J. Nucl. Med. 36 (3), 378-383 (1995).
  47. Ma, Y., Tang, C., Moeller, J. R., Eidelberg, D. Abnormal regional brain function in Parkinson’s disease: truth or fiction. Neuroimage. 45 (2), 260-266 (2009).
  48. Habeck, C., Steffener, J., Rakitin, B., Stern, Y. Can the default-mode network be described with one spatial-covariance network. Brain Res. 1468, 38-51 (2012).
  49. Joliffe, I. T. Principal Components Analysis. Springer Series in Statistics. , (2002).
  50. Limpert, E., Stahel, W. A., Abbt, M. Log-normal distributions across the sciences: keys and clues. BioScience. 51 (5), 341-352 (2001).
  51. Huang, C., Mattis, P., et al. Metabolic abnormalities associated with mild cognitive impairment in Parkinson disease. Neurology. 70 (16 Pt. 2), 1470-1477 (2008).
  52. Mattis, P. J., Tang, C. C., Ma, Y., Dhawan, V., Eidelberg, D. Network correlates of the cognitive response to levodopa in Parkinson disease. Neurology. 77 (9), 858-865 (2011).
  53. Feigin, A., Kaplitt, M. G., et al. Modulation of metabolic brain networks after subthalamic gene therapy for Parkinson’s disease. Proc. Natl. Acad. Sci. U.S.A. 104 (49), 19559-19564 (2007).
  54. Mure, H., Tang, C. C., et al. Improved sequence learning with subthalamic nucleus deep brain stimulation: evidence for treatment-specific network modulation. The Journal of Neuroscience: The Official Journal of the Society for Neuroscience. 32 (8), 2804-2813 (2012).
check_url/50319?article_type=t

Play Video

Cite This Article
Spetsieris, P., Ma, Y., Peng, S., Ko, J. H., Dhawan, V., Tang, C. C., Eidelberg, D. Identification of Disease-related Spatial Covariance Patterns using Neuroimaging Data. J. Vis. Exp. (76), e50319, doi:10.3791/50319 (2013).

View Video