NOTITIE: Alle software die hier wordt gebruikt, wordt weergegeven in de materiaaltabel. De gegevens die in deze studie voor demonstratiedoeleinden werden gebruikt, waren afkomstig van het Human Connectome Project (HCP: http://www. humanconnectome.org)15. Alle experimentele procedures werden goedgekeurd door de Institutional Review Board (IRB) van de Washington University. Beeldvormingsgegevens in de HCP-dataset werden verkregen met behulp van een gemodificeerde 3T Siemens Skyra-scanner met een 32-kanaals kopspoel. Andere parameters voor beeldacquisitie worden gedetailleerd beschreven in een eerder artikel16. Er werden minimale voorbewerkte gegevens gedownload voor de demonstratie, die de volgende voorbewerkingsstappen had voltooid: gradiëntvervormingscorrectie, bewegingscorrectie, voorbewerking van veldkaarten, correctie van ruimtelijke vervorming, ruimtelijke normalisatie naar de ruimte van het Montreal Neurological Institute (MNI), intensiteitsnormalisatie en verwijdering van biasveld. Rusttoestand fMRI-gegevens van de projecten van onderzoekers kunnen ook worden gebruikt.
1. Voorbewerking van de gegevens
- Controleer de gegevenskwaliteit en sluit deelnemers uit met ontbrekende hertestgegevens en overmatige hoofdbewegingen (3 mm in translatie en 3° in rotatie).
OPMERKING: Vijf deelnemers werden verwijderd en 38 jongvolwassenen (22-35 jaar oud) werden opgenomen in de hoofdanalyses.
- Open de graph theoretical network analysis (GRETNA) toolbox17 in MATLAB om verdere voorbewerkingsstappen uit te voeren. Klik op de batch van FC Matrix Construction. Selecteer het pad van de functionele gegevensset om de NIFTI-documenten te laden en voer de volgende stappen uit, zoals weergegeven in de pijplijnoptie in afbeelding 1:
- Verwijder de eerste 10 afbeeldingen door te dubbelklikken op Tijdstip Te verwijderen nummer in Eerste afbeeldingen verwijderen en 10 in te voeren.
- Ruimtelijk glad (volledige breedte op de helft van het maximum [FWHM] = [4 4 4] door te dubbelklikken op FWHM (mm) in Ruimtelijk glad en [4 4 4]).
- Regressie uit covariabelen. Kies Witte stofsignalen, CSF-signalen en Hoofdbeweging als WAAR. Selecteer het juiste masker op basis van de werkelijke voxelgrootte, bijvoorbeeld masker met 2 mm hier, en kies Friston-24-parameters voor Head Motion.
- Tijdelijk filteren. Voer de waarde van TR in op basis van de herhalingstijd van de MRI-scan (bijv. 720 ms hier) en verwijder hoogfrequente en laagfrequente ruis door te dubbelklikken op Band (Hz) en [0.01 0.1] in te voeren.
OPMERKING: Resultaten met en zonder regressie van signalen van het hele brein worden hieronder weergegeven. Bij het gebruik van niet-voorbewerkte gegevens worden ook gevestigde pijplijnen aanbevolen, zoals fMRI-prep18 en Data Processing Assistant for Resting-State fMRI (DPARSF)19 .

Figuur 1: Schatting van het voorproces van Rs-fMRI en functionele netwerkconnectiviteit. De instellingen van het voorproces (het verwijderen van de eerste 10 beelden, ruimtelijk gladstrijken met FWHM van 4 mm, lineair temporeel detrending, regressie van witte stofsignalen, cerebrospinale vloeistof (CSF) signalen, en hoofdbeweging met 24 parameters, filteren van de band van 0,01-0,1 HZ) en de statische correlatie met fisher' Z getransformeerd. Afkortingen: Rs-fMRI = functionele magnetische resonantiebeeldvorming in rusttoestand; FWHM = volledige breedte op half-maximum; CSF = hersenvocht. Klik hier om een grotere versie van deze figuur te bekijken.
2. Netwerkopbouw en -analyses
OPMERKING: De algemene workflow voor de aanleg en analyse van het navigatienetwerk wordt samengevat in drie hoofdstappen (figuur 2).

Figuur 2: Algemene workflow voor de opbouw en analyse van het navigatienetwerk. (A) Kies navigatie als de term waarnaar moet worden gezocht in de Neurosynth-database. (B) Er kan een lijst met activeringscoördinaten worden gegenereerd. (C) Voer een meta-analyse uit met behulp van functies van de Neurosynth om verschillende hersenkaarten te krijgen. (D,E) Door de meta-analytische kaart en een integraal-brein-parcellatieatlas (AICHA) op te nemen, kunnen knooppunten (ROI) worden gegenereerd. (F) De aanleg van een navigatienetwerk met behulp van de resulterende navigatieknooppunten en hun functionele connectiviteit (connectiviteitsraming en netwerkanalyse). Afkortingen: ROI = regio van belang; AICHA = atlas van intrinsieke connectiviteit van homotope gebieden. Klik hier om een grotere versie van deze figuur te bekijken.
- Definitie van netwerkknooppunten
- Download de nieuwste Neurosynth-database (neurosynth.org)20 door het commando in Python te typen:
>Neurosynth importeren als NS
>ns.dataset.download (path='./', unpack = True)
OPMERKING: Het datasetarchief ('current_data.tar.gz') bevat twee bestanden: 'database.txt' en 'features.txt'. Deze bevatten alle activeringscoördinaten van neuroimaging-artikelen en meta-analysetags die respectievelijk met een hoge frequentie in dat artikel voorkomen.
- Genereer een nieuw exemplaar van de gegevensset op basis van de database.txt en voeg functies toe aan deze gegevens door de opdracht te typen:
> van neurosynth.base.dataset gegevensset importeren
> dataset = Dataset('data/database.txt')
> dataset.add_features('gegevens/features.txt')
- Voer een meta-analyse uit met de term van belang (d.w.z. 'navigatie') door het commando te typen:
> ids = dataset.get_ids_by_features ('navigatie', drempelwaarde=0,01)
> ma = meta. MetaAnalyse (dataset, id's)
> ma.save_results('.', 'navigatie')
OPMERKING: De meta-analyse resulteert in verschillende hersenkaarten in NIFTI-formaat. Er werd een FDR-drempel (False Discovery Rate) van 0,01 toegepast om het percentage fout-positieven te controleren. Bij deze stap is gearchiveerde kennis nodig om ervoor te zorgen dat de algemeen gerapporteerde regio's worden opgenomen in de meta-analytische kaart. Soortgelijke stappen kunnen worden toegepast om meta-analyses uit te voeren voor andere cognitieve functies zoals taal en geheugen.
- Definieer interessante clusters door de meta-analytische kaart en een parcellatieatlas voor het hele brein op te nemen door het commando van FSL te typen:
>fslmaths navigation_0.01.nii.gz -bin navi_bin.nii.gz
>fslmaths navi_bin.nii.gz -mul AICHA/AAL.nii.gz navi_label_aicha/aal.nii.gz
>fslmaths navi_label_aicha/aal.nii.gz -thr n -uthr n label _n.nii.gz
>cluster -i label _n.nii.gz -t 0,2 -o cluster_n.nii.gz
>fslmaths cluster_n.nii.gz -thr m -uthr m cluster_n_m.nii.gz
>fslmaths cluster_n_m.nii.gz -bin -mul x node_x.nii.gz
>fslmaths node_1.nii.gz -add ... -add node_x.nii.gz navi_AICHA/AAL_mask.nii.gz
LET OP: Hier zijn twee atlassen gebruikt: AAL en AICHA. De AAL is de atlas die in de oorspronkelijke studie werd gebruikt voor de knooppuntdefinitie12. Deze atlas is gemaakt op basis van de anatomische profielen. De atlas van intrinsieke connectiviteit van homotope gebieden (AICHA)21 heeft een veel fijnere resolutie en een hogere functionele homogeniteit. Aan de hand van elk van de atlassen hebben we de interessegebieden gedefinieerd.
- Typ scripts in Python om de grootte van elke regio op de kaart te controleren:
>For I in np.arange(n)+1:
>____region_list.toevoegen(i)
>____size1_list.append(np.sum(img_dat==i))
>____size2_list.append(np.sum(aicha_img_dat==i))
>____pct_list.append(np.sum(img_dat==i)/np.sum(aicha_img_dat==i))
OPMERKING: Het gehele getal n in het script geeft het totale aantal regio's binnen de AICHA- en AAL-parcellatie aan (respectievelijk 384 en 128). Om de effecten van valse clusters te voorkomen, wordt gesuggereerd dat clusters met relatief kleine afmetingen (bijv. 100 voxels) kunnen worden verwijderd. De AICHA-atlas die hier wordt gebruikt, wordt gegenereerd met behulp van functionele connectiviteitsgegevens, waarbij elke regio homogeniteit van functionele temporele activiteit in zichzelf vertoont.
- Schatting van de netwerkconnectiviteit
OPMERKING: De GRETNA-toolbox wordt gebruikt voor het schatten van connectiviteit en netwerkanalyse.- Klik op de batch van FC Matrix Construction. Laad de voorbewerkte rs-fMRI-gegevens door het pad van de functionele gegevensset te selecteren. Klik op de optie statische correlatie . Upload het knooppunt dat in de vorige stap is verkregen als een atlas om de statische correlatie van rs-fMRI-signalen van elk paar regio's te berekenen en deze over te brengen naar Fisher's z-scores voor het verbeteren van de normaliteit.
NOTITIE: De gedetailleerde werking wordt weergegeven in afbeelding 1. De navigatienetwerkmatrices van N × N (N staat voor het aantal knooppunten) voor elke deelnemer zouden in .txt formaat worden verkregen.
- Krijg een positief en gewogen netwerk met de volgende stappen, zoals weergegeven in figuur 3.
- Klik op de batch Netwerkanalyse. Voeg de netwerkmatrixen toe aan het venster Brain Connectivity Matrix en kies een uitvoermap voor voorbereiding.
- Selecteer voor de pijplijnoptie van Netwerkconfiguratie positief in de tekenmatrix, waardoor negatieve verbindingen in de functieverbindingsmatrix op 0 worden gezet en dubbelzinnige verbindingen22 worden geëlimineerd. Kies het netwerktype als gewogen om het ongerichte gewogen netwerk te krijgen.
OPMERKING: Naast de gewogen netwerken zou men de netwerken ook kunnen binariseren om binaire netwerken te creëren voor latere analyses (met verschillende benaderingen), maar de gewogen netwerken wordt vaak beschouwd als een hogere betrouwbaarheid 23,24.
- Netwerk analyse
- Voeg kleine wereld, globale efficiëntie, clustercoëfficiënt, kortste padlengte, mate centraliteit en lokale efficiëntie toe aan de metrische analysepijplijn van het GRETNA-netwerk, zoals weergegeven in afbeelding 3.
OPMERKING: Kleine wereld en wereldwijde efficiëntie zijn twee wereldwijde netwerkstatistieken. In het bijzonder kan het netwerk met kleinschaligheid de efficiëntie van informatieoverdracht maximaliseren tegen relatief lage bedradingskosten. Globale efficiëntie weerspiegelt de transmissie-efficiëntie parallelle informatie in het transportnetwerk. Voor metrische gegevens van knooppuntnetwerken meet de mate van centraliteit het aantal verbindingen dat met een knooppunt is verbonden. De kortste weglengte, zoals de naam al aangeeft, is een basis voor het meten van integratie. De clustercoëfficiënt geeft aan in welke mate de buren van de knooppunten met elkaar verbonden zijn. Lokale efficiëntie is de efficiëntie van de communicatie met het knooppunt en zijn buren (de gedetailleerde formule en het gebruik worden weergegeven in deze papers) 17,25. Brain connectivity toolbox (BCT)25 en andere toolboxen kunnen ook worden gebruikt voor de berekening van de netwerkstatistieken.
- Selecteer Network Sparsity in de drempelmethode om de verstorende effecten van onechte verbindingen uit te sluiten en voer een reeks drempelreeksen in (d.w.z. 0,05, 0,1, 0,15, 0,2, 0,25, 0,3, 0,35, 0,4, 0,45, 0,5 wordt hier gebruikt) om de juiste drempelwaarde verder te bepalen op basis van de statistische resultaten.
OPMERKING: De verhouding tussen randen en het maximale aantal randen in een netwerk met het steekproefaantal knooppunten staat bekend als de schaarsheidsdrempel. Een sparsiteitsdrempel garandeert dat verschillende individuen hetzelfde aantal randen hebben. We hebben ervoor gekozen om verschillende drempels voor validatie te onderzoeken, wat nuttige gegevens zou kunnen opleveren voor het kiezen van een optimale drempel in toekomstige studies.
- Stel het willekeurige netwerknummer in op 1.000 om willekeurige netwerken te genereren met behulp van een Markov-bedradingsalgoritme26. Klik op Uitvoeren om de pijplijn in GRETNA uit te voeren nadat alle stappen zijn ingesteld.
OPMERKING: Net als bij echte hersennetwerken, behouden de willekeurige netwerken hetzelfde aantal knooppunten, randen en graadverdeling. Om te bepalen of ze topologisch aanzienlijk niet-willekeurig zijn geconstrueerd, worden ze vergeleken met hersennetwerken. Na het uitvoeren van de pijplijn zou een groep scores voor de netwerkstatistieken voor elk van de drempels worden verkregen voor verdere statistische analyses.
- Bepaal in vier stappen het optimale aantal modules in het netwerk.
- Bereken het gemiddelde navigatienetwerk. Klik op het bad van Metrische vergelijking en kies Verbinding. Laad de hierboven verkregen netwerkmatrixen en kies de Gemiddelde (functionele) bewerking. Selecteer een uitgangsrichting om de gemiddelde netwerkmatrix te behouden; zie figuur 4 voor meer details.
- Verdeel het gemiddelde netwerk dat uit de bovenstaande stap is verkregen in 2, 3, 4 en 5 modules met behulp van de functie spectraalcluster in MATLAB.
- Bereken het aandeel knooppunten dat in dezelfde module is verdeeld in REST 1 en REST 2 na het uitlijnen van de module-indelingen met behulp van het script procrustes_alignment.m. Gebruik de verhouding als de index van de herhaalbaarheid van de modulepartitie.
- Selecteer het aantal modules met de hoogste herhaalbaarheid.
- Statistische analyses
OPMERKING: De volgende analyses zijn voornamelijk bedoeld voor validatie en zijn niet nodig bij het toepassen van dit protocol op individuele variatiestudies.- Onderzoek de gelijkenis van deze netwerkmetrieken tussen twee netwerken met verschillende soorten strategieën voor knooppuntdefinitie (d.w.z. de nieuwe die in de huidige studie wordt gegenereerd, aangeduid als NaviNet_AICHA en de eerdere van Kong et al., aangeduid als NaviNet_AAL)12. Bereken de Pearson-correlatie met behulp van de functie corrcoef in MATLAB en herhaal de analyses voor elke sparsiteitsdrempel.
OPMERKING: Na het extraheren van de netwerkstatistieken, kan men alle statistische analyses uitvoeren waarin ze geïnteresseerd zijn.
- Controleer de test-hertestbetrouwbaarheid van deze netwerkmetrieken met behulp van de functie ICC in MATLAB27,28, die de berekening van de intraklassecorrelatiecoëfficiënt implementeert.
OPMERKING: De oorspronkelijke niet-gecorrigeerde p-waarden werden gerapporteerd in de sectie representatieve resultaten. 0,2 < ICC < 0,4 wordt geïnterpreteerd als indicatief voor een redelijke test-hertestbetrouwbaarheid en ICC > 0,4 wordt geïnterpreteerd als matige tot goede test-hertestbetrouwbaarheid29,30. Negatieve ICC-scores werden op nul gezet, gezien het feit dat de aanwezigheid van negatieve ICC's zinloos en moeilijk te interpreteren is31.

Afbeelding 3: Analyse van netwerkstatistieken. Deze analyse definieert de gewogen positieve netwerken met 10 drempels. Bereken twee globale netwerkmaatstaven van klein woord en efficiëntie, vier knooppuntnetwerkmaatstaven van clustercoëfficiënt, kortste padlengte, efficiëntie en mate van centraliteit. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4: De berekening van de gemiddelde vaarnetwerken. De gemiddelde (functionele) operatie helpt bij het berekenen van de gemiddelde netwerken van alle deelnemers. Klik hier om een grotere versie van deze figuur te bekijken.