Ricostruzioni su larga scala e indipendente, imparziale Clustering Sulla base morfologici metrica per classificare i neuroni in selettivi Popolazioni

Riepilogo

Questo protocollo descrive ricostruzioni su larga scala di popolazioni neuronali selettivi, etichettati dopo l'infezione retrograda con esprimendo un virus della rabbia modificato marcatori fluorescenti, e le analisi di cluster indipendente, imparziale che consentono la caratterizzazione completa delle metriche morfologiche tra sottoclassi neuronali distinti.

Abstract

Questo protocollo delinea le ricostruzioni su larga scala di neuroni in combinazione con l'uso di cluster indipendente e imparziale analisi per creare un'indagine completa delle caratteristiche morfologiche osservate in una popolazione neuronale selettiva. La combinazione di queste tecniche costituisce un nuovo approccio per la raccolta e l'analisi dei dati neuroanatomici. Insieme, queste tecniche consentono di larga scala, e quindi più completo, il campionamento delle popolazioni neuronali selettivi e stabilire metodi quantitativi per la descrizione imparziali classi neuronali morfologicamente unici all'interno di una popolazione.

Il protocollo descrive l'uso di virus della rabbia modificato per etichettare selettivamente i neuroni. G-cancellato atti virus della rabbia come tracciante retrogrado seguenti iniezione stereotassica in una struttura di destinazione cerebrale di interesse e serve come veicolo per l'erogazione e l'espressione di EGFP in neuroni. Un gran numero di neuroni sono infettate con questotecnica e di esprimere la GFP per tutta la loro dendriti, la produzione di "Golgi-like" riempimenti completi dei singoli neuroni. Di conseguenza, il metodo di analisi retrograda virus-mediata migliora le tecniche di tracciamento retrogrado tradizionali a base di coloranti con la produzione di riempimenti intracellulari completi.

I singoli neuroni ben isolate che abbracciano tutte le regioni della zona del cervello in fase di studio sono selezionati per la ricostruzione al fine di ottenere un campione rappresentativo di neuroni. Il protocollo delinea le procedure per ricostruire corpi cellulari e completare i modelli arborizzazione dendritiche di neuroni marcati che abbracciano molteplici sezioni di tessuto. dati morfologici, comprese le posizioni di ogni neurone all'interno della struttura del cervello, vengono estratti per ulteriori analisi. funzioni di programmazione standard sono stati utilizzati per eseguire gruppo indipendente analisi e valutazioni cluster basato su metriche morfologiche. Per verificare l'utilità di queste analisi, valutazione statistica di un perfo cluster analysisconfermate su 160 neuroni ricostruite nel nucleo reticolare del talamo del talamo (TRN) del macaco è stato fatto. Sia la cluster analysis originale e le valutazioni statistiche svolte qui indicano che i neuroni TRN sono separati in tre sottopopolazioni, ognuno con caratteristiche morfologiche uniche.

Introduzione

Neuroanatomy è uno dei fondamenti delle neuroscienze ¹ e recente interesse "connectomics" ha rinnovato entusiasmo per la comprensione della diversità morfologica delle popolazioni neuronali e le connessioni tra i neuroni specifici ^2. Metodi per l'etichettatura e neuroni ricostruendo hanno notevolmente migliorato con le recenti innovazioni, tra cui genetica e virus-mediata circuito di tracciamento approcci ^{3, 4,} consentendo studi morfologici più complete delle popolazioni neuronali ^5. Oltre a miglioramenti nella etichettatura singoli neuroni, sono anche emerse tecniche di analisi quantitativa dei dati che consentono di classificazione indipendente e imparziale dei neuroni in sottopopolazioni distinte in base a dati morfologici ^{5, 6.} Queste tecniche imparziali sono un miglioramento su più traditional metodi di classificazione qualitativi che sono stati lo standard nel settore per più di un secolo. L'obiettivo di questo studio è quello di delineare, passo dopo passo, la combinazione di etichettatura virus-mediata di neuroni all'interno di una popolazione selettivo, ricostruzioni su larga scala di un campione globale di questi neuroni, e l'analisi quantitativa dei dati basata su cluster indipendente con valutazione statistica. Grazie alla combinazione di questi metodi, delineiamo un nuovo approccio verso la raccolta e l'analisi dei dati neuroanatomici da facilitare la campionatura completa e imparziale la classificazione dei tipi di neuroni morfologicamente unici all'interno di una popolazione neuronale selettiva.

Come esempio di questi metodi, descriviamo la nostra analisi di una vasta popolazione di neuroni all'interno di un singolo settore del nucleo reticolare del talamo (TRN) del macaco. Questi dati provengono da uno studio preliminare ^7. Metodi per l'etichettatura selettivamente i neuroni TRN proiettano al Latera dorsalil genicolato nucleo del talamo (dLGN) con iniezione chirurgica del virus della rabbia modificato codifica EGFP ^{4, 8} (vedi tabella di specifici materiali / attrezzature, fila 2) sono delineati. Questo virus della rabbia modificato manca il gene che codifica una proteina essenziale cappotto, eliminando movimento trans-sinaptica del virus. Una volta che il virus entra terminali degli assoni nel sito di iniezione, si comporta come un tracciante tradizionale retrograda con l'importante vantaggio di guidare espressione EGFP durante l'intero arborization dendritica dei neuroni infettati ^{5, 9, 10.} Di conseguenza, questo virus della rabbia G-cancellato può essere utilizzato per infettare selettivamente etichettare qualsiasi popolazione neuronale dopo l'iniezione e trasporto retrogrado.

Per eseguire un'analisi completa di una specifica popolazione neuronale, è importante campionareun'ampia distribuzione dei neuroni all'interno della popolazione. Poiché la tecnica di etichettatura virus-mediata produce completa intracellulare, "Golgi-like" riempie di molti neuroni con assoni nel sito di iniezione virus, è possibile ricostruire un campione molto ampio di neuroni nella misura massima di una struttura cerebrale. Inoltre, poiché il virus della rabbia modificato è così efficace a infettare ed etichettatura gran numero di neuroni, è possibile ricostruire centinaia di neuroni per animale. Procedure di campionamento 160 neuroni tutto il settore visivo del TRN ¹¹ per generare un campione completo di dLGN sporgenti neuroni TRN sono delineati. Il processo di ricostruzione neuroni individuali utilizzando un sistema di ricostruzione neurone compreso un microscopio, fotocamera e software di ricostruzione è descritto. Anche descritto sono metodi per determinare le posizioni di singoli neuroni all'interno di una struttura del cervello (in questo caso all'interno del TRN) e per verificare il virus iniezione site dal volume e all'interno di una struttura (in questo caso all'interno del dLGN) utilizzando ricostruzioni contorno volumetriche. I passaggi per esportare i dati morfologici ed eseguire gruppo indipendente analisi basate su metriche morfologiche misurati per ogni neurone sono delineati. Ci sono limitazioni per metodi di clustering e ci sono anche una varietà di diversi algoritmi di clustering disponibili. Di conseguenza, sono descritte le opzioni ei vantaggi di alcuni degli algoritmi più comunemente usato. L'analisi dei cluster non prevede la verifica statistica della unicità dei cluster. Pertanto, passaggi aggiuntivi sono delineati verificare raggruppamento ottimale nonché le relazioni tra i dati morfologici all'interno e tra cluster. Metodi statistici per la valutazione cluster per il set di dati TRN per confermare che i neuroni TRN sono raggruppati in tre gruppi unici basato su 10 parametri morfologici indipendenti sono descritte.

Così, delineando passi per selettivamente l'etichettatura, Ricostruendo, e l'analisi dei dati morfologici da una specifica popolazione neuronale, si descrivono i metodi per quantificare le differenze morfologiche tra i neuroni all'interno di una popolazione. risultati precedenti di tipi neuronali distinti all'interno del settore visivo del TRN scimmia macaco sono confermate con separati metodi di valutazione statistica. Insieme, ci auguriamo che queste tecniche saranno ampiamente applicabile a set di dati neuroanatomiche e aiutare a stabilire la classificazione quantitativa della diversità delle popolazioni neuronali attraverso il cervello.

Protocollo

Nota: Il tessuto esaminato in questo studio è stato preparato come parte di uno studio separato ^5. Pertanto, tutti i metodi sperimentali che prevedono l'uso di animali sono stati descritti in dettaglio nella sezione Metodi sperimentali di Briggs et al. (2016). Tutte le procedure che coinvolgono gli animali condotti come parte dello studio preliminare sono stati approvati dai comitati cura degli animali e uso istituzionale. sono descritti brevemente di seguito nelle sezioni 1 La procedura per l'iniezione del virus nel dLGN e l'elaborazione istologica del tessuto cerebrale - 2.

1. Stereotassica Iniezione

1. Eseguire tutte le procedure chirurgiche in un ambiente sterile usando tecniche asettiche. Steam-sterilizzare tutti gli strumenti chirurgici di metallo, garze e materiale drappeggio in autoclave. Sterilizzare tutti materiali che potrebbero essere danneggiati da vapore con la sterilizzazione chimica (per esempio ossido di etilene).
2. Indurre e mantenere l'anestesia in base al animale- e prequisiti specifici ROTOCOLLO.
   1. Per l'iniezione di virus nelle scimmie, indurre l'anestesia con ketamina (10 mg / kg) e mantenere gli animali in anestesia chirurgica complete con isoflurano (1-3% inalato in ossigeno) monitoraggio scaduto CO _2, ECG, frequenza respiratoria, e la temperatura costantemente e periodicamente verifica per il tono della mascella per garantire la corretta profondità dell'anestesia.
3. Posizionare animale in un telaio stereotassico per stabilizzare la posizione della testa e per consentire l'uso di coordinate stereotassica per localizzare la struttura di iniezione di interesse (ad es dLGN). Se la misurazione risposte visive, posto collirio (1% occhio atropina scende se la dilatazione della pupilla è desiderata, altrimenti collirio soluzione salina) e poi le lenti a contatto in entrambi gli occhi per prevenire la secchezza. Se non misurare le risposte visive, posizionare pomata oftalmica negli occhi.
4. Effettuare una incisione sulla linea mediana del cuoio capelluto e ritrarre la pelle ei muscoli.
5. Secondo stereotassica coordinate per la struttura di interesse (dLGN) in uno hemisphere, fare una piccola craniotomia.
6. Abbassare gradualmente un elettrodo di registrazione (ad esempio tungsteno o platino / iridio elettrodo (vedi Tabella di specifici materiali / attrezzature, riga 3) fissata ad una micromanipolatore che è posizionato manualmente) attraverso la craniotomia nel cervello. Bassa impedenza (<1 MW) e leggermente più spesso (~ 250 micron di diametro) elettrodi possono essere utilizzati per penetrare la dura, se desiderato.
7. Registrare la posizione del manipolatore in superficie corticale e nel punto in cui vengono misurate risposte visive. risposte visive sono udibili, tempo-locked risposte neuronali alla luce balenò negli occhi con un oftalmoscopio o piccolo LED / torcia elettrica.
8. Determinare la posizione e lo spessore del dLGN notando posizioni manipolatore all'inizio, metà e fine risposte visive e sottraendo la posizione della superficie corticale da questi valori. Registrare la profondità per ciascun sito di iniezione progettato all'interno del dLGN.
   NOTA: Operazi simileres possono essere usate per localizzare strutture non visivi utilizzando appropriati stimoli sensoriali. Inoltre, i sistemi di registrazione / iniezione accoppiato può essere utilizzato anche per indirizzare le strutture cerebrali piccoli.
9. Una volta che sono noti località iniezione ottimali, rimuovere l'elettrodo di registrazione e inserire un pipetta di iniezione (pipetta di vetro o siringa) nella stessa posizione stereotassica ed inferiore alla profondità appropriata per ogni iniezione, iniziando con il sito di iniezione più profonda e procedendo alla superficiale, in attesa di almeno 1 minuto dopo l'iniezione prima di cambiare la posizione degli elettrodi. pipette di vetro con diametro esterno di ~ 50 micron possono ridurre il riflusso del virus rispetto alle punte siringa (~ 200 micron di diametro).
10. L'utilizzo di un sistema di iniezione (vedi Tabella di materiali specifici / Equipmen t, fila 4), picospritzer o pompa a siringa, iniettare piccoli volumi (1 - 5 ml, fare riferimento alle istruzioni del produttore per le procedure di iniezione) del virus della rabbia modificato su più (3 - 10) Depthsad una velocità di ~ 100 nL / min all'interno della struttura di destinazione (dLGN). Nota: penetrazioni di iniezione separati possono essere realizzati in strutture bersaglio più grandi (ad esempio scimmia dLGN). Regolare il volume di iniezione totale in base alle dimensioni struttura obiettivo.
11. Pausa almeno 5 minuti prima di ritrarre la pipetta di iniezione. Quando la pipetta di iniezione viene rimosso, riempire la craniotomia con materiale protettivo (pezzo di osso colla in luogo, applicare la cera osso o Gelfoam), poi suturare il muscolo e pelle di nuovo insieme.
12. Somministrare analgesici (ad esempio, ketoprofene) e gli antibiotici prima della fine di un intervento chirurgico e monitorare animale continuamente fino a quando animale è ambulatoriale.
13. Per almeno 3 giorni e fino a 10 giorni dopo l'intervento, il monitoraggio degli animali quotidianamente per garantire suture siano pulite, asciutte, e intatta e animale non mostra segni di dolore o fastidio. Somministrare analgesici e antibiotici ogni giorno, come richiesto. Inizia edilizia sociale di animali post-chirurgica dopo animale ha pienamente recuperato da un intervento chirurgico.

2. Tissue raccolta, sezionamento e colorazione

1. Consentire 7 - 14 giorni per il trasporto retrogrado di virus, quindi euthanize l'animale per overdose di soluzione di eutanasia (per esempio Euthasol) e profumato l'animale transcardially con 0,1 M tampone fosfato salino, 4% paraformaldeide, quindi paraformaldeide 4% con il 10% di saccarosio.
2. Rimuovere il cervello (vedi ¹² per i passaggi analoghi in un modello di roditore) e il luogo in 20 - 30% di saccarosio con paraformaldeide al 4% in 0.1 M tampone fosfato e conservare in frigorifero per 1 - 2 giorni.
3. Quando il cervello ha toccato il fondo del contenitore, rimuoverlo e tagliare il tessuto in blocchi contenenti la regione del cervello con neuroni retrogrado etichettati (es corteccia visiva) e la struttura porta di iniezione (es dLGN). Tagliare gli stessi blocchi di tessuto dall'emisfero non per iniezione - queste serviranno come sezioni di controllo. Per ottenere risultati ottimali, iniziare le procedure istologiche immediately sul tessuto appena sezionato.
4. Sezione di tessuto coronale con uno spessore di 50 micron per sezione con un microtomo di congelamento (vedi tabella di specifici materiali / attrezzature, riga 5).
5. Macchiare tutte le sezioni per l'attività del citocromo ossidasi (vedi Tabella di materiali specifici / Attrezzature, righe 5 - 8) Per visualizzare strati corticali e strutture sottocorticali ^13. Nota: Le sezioni possono essere memorizzati durante la notte in un frigorifero in risciacquo finale dopo la macchia citocromo ossidasi.
6. Etichettare tutte le sezioni con un anticorpo primario contro la GFP (1: 1000 diluizione, ~ 12 ore di incubazione, si veda la Tabella di specifici materiali / attrezzature, fila 9) seguito da un anticorpo secondario che corrisponde al host primario e contrassegnati con biotina ⁵ (1: 500 diluizione 2 - 4 h di incubazione, vedi Tabella di specifici materiali / attrezzature, riga 10). Si raccomanda di incubare le sezioni in una notte anticorpo secondario.
7. Sezioni Label con un complesso avidina / biotina poi reagiscono con DAB e perossido di macchiare in modo permanente tutti i neuroni marcati ^5.
8. Montare le sezioni su vetrini Subbed e lasciare asciugare per una notte.
9. Sezioni sgrassare utilizzando una serie di alcol e xilene risciacqui quindi coprire slittamento ^14.

3. Ricostruzione neuronale

NOTA: Tutte le ricostruzioni di esperimenti originali sono stati realizzati utilizzando un sistema di ricostruzione dei neuroni costituito da un microscopio (vedi tabella di specifici materiali / attrezzature, riga 14), fotocamera collegata (vedi Tabella di specifici materiali / attrezzature, riga 13), e il software di ricostruzione pacchetto (vedi tabella di materiali specifici / Attrezzature, righe 11 - 12). Software-assistita ricostruzione neuronale permette la visualizzazione dei vetrini di tessuto sovrapposti con disegni basati su computer di processi neuronali. È importante sottolineare che il software digitalizza ricostruzione morfologicaDati zione in tre dimensioni, che permette l'estrazione di informazioni morfologiche specifiche posizioni. Il programma di estrazione di dati associati (vedi Tabella di specifici materiali / attrezzature, riga 12) consente l'estrazione di un ricco set di dati morfologici di ogni ricostruzione salvato.

1. Mettere una diapositiva nel supporto vetrino da microscopio e concentrarsi su una singola sezione di interesse utilizzando un obiettivo a basso ingrandimento (ad es 2X o 10X). Assicurarsi che l'immagine della telecamera è visibile nella vista software posizionando correttamente l'otturatore della fotocamera.
2. Scegliere neuroni marcati all'interno della struttura cerebrale di interesse (ad es TRN) che sono ragionevolmente isolato per ricostruire inequivocabilmente tanto del arborization dendritica possibile. Preferibilmente scegliere neuroni se il corpo cellulare è completamente all'interno di una singola sezione al fine di catturare la più grande estensione di ogni corpo cellulare e stimare con precisione la sua area e rotondità. Utilizzare 10X obiettivo microscopio per identificare cellule Body nella sezione casa.
3. Una volta che un ben colorato, ben isolato-neurone viene identificato, aggiungere la sezione "casa" che contiene il corpo cellulare alla sezione Manager facendo clic sul pulsante "Nuova sezione". Inserire il numero di sezioni da includere (raccomandare 2 - 3 sezioni per cominciare). Assegnare un spessore della sezione di 50 micron quando richiesto (vedere il punto 2.3).
4. Tracciare i contorni delle strutture cerebrali pertinenti nella sezione casa contenente il corpo cellulare. Utilizzare il 2X obiettivo / ingrandimento per il contorno tracciato.
   1. In primo luogo, selezionare "Contorno" dal menu a discesa sulla barra degli strumenti superiore e quindi selezionare il tipo di contorno dal menu a discesa (ad esempio "LGN" o "TRN").
   2. Tracciare contorni della struttura di destinazione (es TRN) così come le strutture adiacenti, se desiderato, facendo clic sui punti lungo il contorno con il mouse come strumento tiraggio.
   3. Selezionare "Chiudi Contour" facendo clic destro del mouse e selezionando questa opzionedal menu per completare e racchiudere ogni contorno tracciato.
5. Posizionare un marcatore al centro della struttura di destinazione facendo clic sul simbolo marcatore desiderato sulla barra degli strumenti a destra e poi cliccando sul punto desiderato nella ricostruzione per posizionare il marcatore in quella posizione. Utilizzare lo stesso tipo di marcatore per contrassegnare il centro della struttura bersaglio in tutte le ricostruzioni.
6. Una volta che i contorni sono completi, tracciare il contorno del corpo cellulare utilizzando 40 - 60X obiettivi / ingrandimento. Per fare ciò, prima selezionare "Neuron" dal menu a discesa sulla barra degli strumenti superiore e quindi selezionare la struttura neurone per tracciare, in questo caso "Cell Corpo". Successivo tracciare il corpo cellulare come in 3.4.
7. Posizionare un marcatore stile diverso al centro del corpo cellulare, seguendo la procedura descritta in 3.5.
8. Tracciare tutti i dendriti che iniziano al corpo cellulare.
   1. In primo luogo, selezionare "Dendrite" dal menu a discesa. Poi tracciare ogni dendrite a partire dalla CdA cellularey e utilizzando il mouse come strumento tiraggio. Assicurarsi di regolare la profondità Z tutto il tracciato per catturare con precisione l'angolo e la direzione del dendrite.
   2. Posizionare un nodo in ogni nodo lungo il dendrite facendo clic destro del mouse e selezionando "biforcano nodo" o "triforcazione nodo" dal menu a discesa. Utilizzare 40 - 60X obiettivo / ingrandimento per tutte le ricostruzioni dendriti.
   3. Alla fine di ogni dendrite, fare clic con il mouse e selezionare "Fine" dal menu a discesa.
9. Identificare terminazioni dendritiche che possono continuare nella sezione adiacente e portare questi a fuoco al appropriata z approfondita l'immagine al microscopio a. Includere i principali punti di riferimento nelle vicinanze, come i vasi sanguigni o facilmente riconoscibili i modelli dendriti o fasci l'immagine al microscopio a.
   1. Ridurre l'ingrandimento di 20X o 10X sul microscopio (scegliere l'ingrandimento che permette la visualizzazione più grande punto di riferimento), e prendereuna foto del immagine al microscopio utilizzando una macchina fotografica digitale (ad esempio telefono cellulare, tablet) tenuto in mano allo schermo del computer.
10. Spostare alla sezione adiacente della diapositiva e allineare i contorni della sezione precedentemente tracciato con i confini del dLGN e TRN nella nuova sezione.
11. Ritorno il campo di vista per la zona del corpo cellulare (basato su contorni e grandi punti di riferimento) e portare l'ingrandimento indietro fino a 10 - 20X.
12. Utilizzare la foto delle terminazioni dendrite dalla sezione precedentemente tracciata per aiutare ad allineare le terminazioni con l'inizio dei dendriti nella nuova sezione.
    1. Per ruotare il tracciato e spostarlo per allineare dendriti, utilizzare lo strumento freccia per selezionare la ricostruzione, tasto destro del mouse e selezionare "Sposta" dal menu a discesa.
    2. In alternativa, utilizzare lo strumento "Match" per abbinare terminazioni dendritiche a inizi. Selezionare lo strumento Partita dalla barra degli strumenti superiore e, quando richiesto, fare clic sul finalenella ricostruzione e il corrispondente punto di prosecuzione nell'immagine. Ripetere per 3 o più finali.
13. Una volta che il tracciato della sezione precedente è allineato con i dendriti nella nuova sezione, assicurarsi che la nuova sezione corrispondente viene selezionata nel gestore sezione facendo clic sulla sezione corrente. Oppure, aggiungere una nuova sezione alla sezione Manager, come sopra descritto in 3.3, regolando la profondità Z e la posizione di ogni sezione nuova rispetto alla sezione d'origine e impostando ogni spessore di sezione 50 mm (vedi punto 2.3).
14. Aumentare l'ingrandimento di 40 - 60X e tracciare le continuazioni dendrite cliccando col tasto destro del mouse sulla fine del dendrite dalla ricostruzione e selezionando "Aggiungi a porre fine" dal menu a discesa, quindi tracciare il dendrite utilizzando il mouse come disegnare strumento. Seguire richiesta di specificare se continuazione è nella nuova sezione.
15. Seguire dendriti attraverso almeno 3 sezioni adiacenti (uno su ciascun lato della sezione di casa) a seguito di passaggi 3,8-3,15. Aggiungere alla ricostruzione fino almeno 3 sezioni totali sono tracciati per il neurone o finché i dendriti non possono più essere seguite o trovato. Tracciare i contorni delle sezioni limitrofe seguito i passaggi precedenti, se lo si desidera.

4. indipendente Clustering

Nota: Gruppo indipendente di analisi consentono analisi imparziali di grandi insiemi di dati, multi-dimensionali che altrimenti potrebbero essere difficili da visualizzare e, soprattutto, fornire una valutazione quantitativa della diversità morfologica. Una piattaforma di programmazione matriciale è molto utile per l'analisi di set di dati multi-dimensionali e consente sofisticate manipolazioni di dati e analisi statistiche. Le funzioni elencate nei punti 4 - 6 sono definiti nella piattaforma di programmazione elencati nella tabella di specifici materiali / attrezzature, riga 15.

1. Estrarre i dati morfologici da ogni singolo ricostruzione neuronale utilizzando un extprogramma raction associata al sistema di ricostruzione dei neuroni (vedi tabella di materiali specifici / Attrezzature, righe 11 - 12).
   1. In primo luogo, scegliere i dati morfologici desiderati per ogni ricostruzione tra cui: le informazioni di contorno, le informazioni marcatore, dati morfologici del corpo cellulare e arborizzazione dendritica, ecc.
   2. Dal menu a discesa Modifica nella barra degli strumenti superiore, selezionare "Seleziona tutti gli oggetti". Quindi selezionare "Analisi struttura ramificata" dal menu a discesa Analysis nella barra degli strumenti in alto e cliccare sulla scheda e selezionare le opzioni di analisi desiderate in ogni scheda.
   3. Estrarre tutti i dati desiderati facendo clic sul pulsante "OK" nella finestra di analisi e salvare in un formato foglio di calcolo facendo clic destro sulle finestre di uscita e selezionando "Esporta in Excel".
2. Compilare un foglio di calcolo Master (vedi Tabella di specifici materiali / attrezzature, riga 16), in cui ogni variabile morfologica / meTric è rappresentata da una singola osservazione numerica per neurone. In questo foglio di master, tenere un registro di identificatore riga per ogni neurone.
3. Verificare che tutte le variabili incluse nel cluster analysis (es metriche morfologiche) sono indipendenti l'uno dall'altro. Ad esempio, il numero di nodi sarà correlata con il numero di rami in un arborization dendritica e assonale. Di conseguenza, solo una di queste due variabili dovrebbero essere inclusi in un cluster analysis.
4. Assicurarsi che ogni misura contiene una osservazione per tutte le variabili, vale a dire ogni neurone (di misura), deve avere un punto di dati (osservazione) corrispondente ad ogni metrica morfologica (variabile). Se le osservazioni di una particolare variabile (per esempio, lunghezza assone) sono carenti per un sottoinsieme di neuroni, questa variabile (lunghezza assone) non può essere incluso nell'analisi cluster.
5. Organizzare i dati sul foglio di calcolo master in una singola matrice in cui ogni riga rappresenta un neurone e each colonna contiene dati numerici per ciascuna metrica morfologica (Figura 1). Ad esempio, un insieme di dati di cui 100 neuroni per i quali vi sono osservazioni per 5 variabili indipendenti sarebbe rappresentato da una matrice 100 x 5 (riga per colonna). Non è necessario includere qualsiasi identificazione di ogni neurone all'interno della matrice - l'indice di riga servirà come identificatore univoco di ogni neurone. Ci dovrebbe essere anche senza lacune o s "nan" nella matrice. Salvare la matrice come una singola variabile (ad esempio, dati = [100 x 5 matrice]; vedi file di codice supplementare per il codice di esempio).
6. Scegliere un algoritmo per calcolare le distanze tra i punti che rappresentano singoli neuroni in uno spazio n-dimensionale, dove n è definito dal numero di variabili. La funzione 'pdist' permette flessibilità nel calcolo distanze tra neuroni. Il calcolo della distanza di default è la distanza euclidea ed è stato utilizzato in precedenza per le analisi simili di morphologica neuronalel dati ^{5, 6.}
7. Con il tra-neurone allontana uscita della funzione 'pdist', utilizzare la funzione di 'linkage' per formare gruppi. Anche in questo caso, sono disponibili diverse opzioni di clustering. Il metodo di Ward e la distanza baricentro hanno dato risultati simili nelle analisi di set di dati morfologici ^{5, 6.}
8. Se lo si desidera, utilizzare la funzione 'Kmeans' come un approccio di clustering alternativa al 'pdist' e funzioni 'linkage'. '' Kmeans comuni utilizza squadrati distanza euclidea per calcolare distanze tra neuroni e quindi la distanza baricentro per assegnare cluster.
9. Per visualizzare un albero di cluster gerarchica, utilizzare la funzione 'dendrogramma' sull'uscita dell'operazione 'linkage'. Questa funzione ha un valore predefinito che limita la visualizzazione a 30 misurazioni, ma può essere sovrascritto specificando il numero di misure ( ad esempio i neuroni) visualizzati. Il dendrogramma illustra le distanze di collegamento sul l'asse y per ciascuno dei neuroni, identificato dal loro indice di riga sul asse x.
   NOTA: Le uscite di 'linkage' e 'dendrogramma' non definiscono il numero ottimale di cluster. L'uscita di '' Kmeans non assegna misure ai cluster, ma non verifica se questi gruppi sono ottimali. confronti statistici separate e la verifica di raggruppamento ottimale possono essere eseguiti (vedi sotto).

5. Verifica dei Clustering

Nota: Come indicato in precedenza, l'analisi cluster stesso non fornisce direttamente una valutazione statistica dei se i cluster illustrati nel dendrogramma cluster siano unici e rappresentativo del campione. I metodi per la verifica cluster dal dendrogramma sono stati proposti ^15, tuttavia questi non forniscono verifica statistica di raggruppamento ottimale. Ci sono multimetodi PLE per verificare il clustering ottimale.

1. Metodo 1: Valutare il clustering:
   1. Utilizzare la funzione '' evalclusters specificando il metodo di calcolo cluster, come ad esempio 'Kmeans' o 'legame'. Una gaussiana il risultato del modello miscela può anche essere valutata (vedi punto 5.2; vedere file di codice supplementare per il codice di esempio).
   2. Specificare il criterio di valutazione da una serie di opzioni (ad esempio 'CalinskiHarabasz' - per la descrizione delle opzioni di criterio, consultare il menu Guida in linea, alla ricerca di "evalclusters").
   3. Specificare un intervallo di numeri di cluster ottimali per testare (es [1: 6] per verificare se 1, 2, 3, 4, 5, o 6 cluster sono ottimali).
      Nota: L'uscita di 'evalclusters' è un numero ottimale di cluster dato l'algoritmo di clustering e il criterio specificato.
2. Metodo 2: gaussiana modello mistura di clustering:
   NOTA: miscela gaussiana modello di clustering Algorithms (MGM) supporre che le osservazioni per ogni variabile provengono da una miscela di distribuzioni gaussiana che definiscono ogni cluster putativo, ognuno con la propria media e covarianza. Il GMM quindi utilizza un algoritmo em assegnare probabilità a posteriori ad ogni misurazione, indicando la probabilità di appartenere ad uno specifico cluster ^16.
   1. Implementare una GMM utilizzando la funzione 'fitgmdist' inserendo il numero putativo di cluster osservabili nei dati. In alternativa, per evitare l'assegnazione a priori il numero di cluster, utilizzare analisi delle componenti principali (PCA) con una serie di diversi numeri di cluster ottimale utilizzando la procedura descritta qui (vedi file di codice supplementare per il codice di esempio).
   2. Utilizzare la funzione 'PCA' sulla matrice dei dati originali e salvare i punteggi dei componenti principali come uscita.
   3. In un ciclo a partire da 1 e passando attraverso un numero di cluster ottimali putativi, usare thfunzione 'fitgmdist' e sui punteggi dei componenti principali e generare MGM per ogni numero di cluster ottimali putativi.
   4. Valutare ogni GMM esaminando la probabilità negativo registro, Akaike criterio di informazione, e di Bayes criterio di informazione. Il GMM con i criteri più bassi è ottimale.
   5. Se lo si desidera, verificare se i valori medi per ogni cluster sono significativamente differenti quando il numero di cluster è ottimale (mezzi di ogni cluster sono memorizzati nel "mu" uscita di ciascun GMM).

6. Analisi statistica dei dati in cluster

1. Una volta che il numero ottimale di cluster è determinato utilizzando il clustering gerarchico e valutata, tornare i dati originali, i neuroni separati in gruppi, e determinare quali parametri morfologici contribuiscono al raggruppamento di neuroni in classi uniche.
2. Ordinare i dati metrici morfologici originali, che i neuroni sono raggruppati in base alla assegnazione cluster.
3. Per esaminare le relazioni tra i dati morfologici tra neuroni (per tutti i neuroni o per i neuroni separati in cluster), eseguire la regressione lineare si adatta a 2 e 3 vie comparazioni delle osservazioni metriche morfologiche. Questi attacchi possono essere stimati utilizzando la funzione di 'forma' e valutati utilizzando la funzione 'fitlm', in cui la bontà di uscite adatte comprendono R ² e p valori per attacchi di regressione.
4. Per esaminare le relazioni statistiche tra i neuroni in ogni cluster, utilizzare non parametrico due campioni (Wilcoxon rank test-sum) o multipla-campionari di confronti (ANOVA), a seconda del numero di cluster. È importante sottolineare che l'uso di più campioni ANOVA assicura che i valori p sono corretti per confronti multipli.

Risultati

Abbiamo dimostrato in precedenza che le ricostruzioni su larga scala di neuroni all'interno di una popolazione selettivo sta seguendo iniezione fattibile del virus della rabbia modificato nel dLGN ^5. Recentemente, lo stesso tessuto è stato utilizzato per ricostruire 160 neuroni nel settore visivo del TRN (Bragg et al, in revisione;. Figure 2A-B) seguendo le fasi metodologiche dettagliate qui sopra descritti. Nello studio TRN, tre gruppi ...

Discussione

studi neuroanatomici sono rimasti un pilastro di neuroscienze e recente interesse connectomics e le relazioni struttura-funzione ha rinnovato entusiasmo per la caratterizzazione morfologica dettagliata delle popolazioni di neuroni selettivi. Tradizionalmente, gli studi neuroanatomici hanno fatto affidamento su classificazioni qualitative di neuroni in classi morfologicamente distinte di neuroni definiti da neuroanatomisti esperti. Con i progressi nelle tecniche per la ricostruzione neuroni e l'estrazione dei dati mo...

Materiali

Name	Company	Catalog Number	Comments
SADΔG-EGFP	E.M. Callaway Laboratory, Salk Institute		Prepared by Dr. F. Osakada. G-deleted rabies virus available through the Salk Institute Viral Core
Recording electrode: platinum/iridium or tungsten	FHC	UEPSGGSE1N2M	Visit website (www.fh-co.com) for alternative order specifications
Nanoject II	Drummod Scientific	3-000-204, 110V	Alternatives: picospritzer, Hamilton syringe
Freezing microtome	Thermo Scientific
DAB	Sigma Aldrich	D5905-50TAB	3,3'-Diaminobenzidine tetrahydrochloride, tablet, 10 mg substrate per tablet. Caution: carcinogen - must be bleached before discarding
Cytochrome C	Sigma Aldrich	C2037-100MG
Catalase	Sigma-Aldich	C9322-5G
Rabbit anti-GFP	Life Technologies/Thermo Fisher	#A-11122	Primary antibody
Biotinylated goat anti-rabbit	Vector Laboratories	#BA-1000	Secondary antibody
Neurolucida System	MicroBrightField		Software for neuron tracing and analysis. http://www.mbfbioscience.com/neurolucida
Neurolucida Explorer	MicroBrightField		Data export software
Microfire Camera	Optronics		2-Megapixel true color microscope camera. http://www.simicroscopes.com/pdfs/microfire.pdf
Nikon E800 Microscope	Nikon Instruments Inc.		Biological research microscope. http://www.microscopyu.com/museum/eclipseE800.html
Matlab	The MathWorks Inc.		Matrix-based computational mathematics software. http://www.mathworks.com
Microsoft Office Excel	Microsoft		Spreadsheet program