A Generalized Framework of Antisymmetric Polyspectral Indices for Identifying High-Order Neural Interactions

Technische Samenvatting: Een Gegeneraliseerd Kader van Antisymmetrische Polyspectrale Indices voor het Identificeren van Hoog-orde Neuronale Interacties

Probleemstelling
Cross-frequency coupling (CFC) is een fundamenteel mechanisme voor het integreren van informatie over temporale schalen in de hersenen. Het nauwkeurig identificeren van echte hoog-orde neuronale interacties wordt echter gehinderd door twee primaire uitdagingen: complexe multi-frequentie niet-lineariteiten en schijnbare, nul-lag artefacten veroorzaakt door volumgeleiding en veldverspreiding in EEG/MEG-opnames. Conventionele connectiviteitsmetrieken (bijv. PLV, ImCoh, wPLI) gaan doorgaans uit van interacties op dezelfde frequentie en falen in het karakteriseren van afhankelijkheden waarbij een frequentie van belang $f_N$ voortkomt uit de combinatie van $N-1$ componenten ($f_N = \sum_{i=1}^{N-1} f_i$). Bovendien ontbreekt het bestaande cross-frequency maatstaven vaak aan intrinsieke robuustheid tegen lineaire instantane menging, generaliseren ze niet op natuurlijke wijze tot interacties tussen meer dan twee regio's, en scheiden ze fase en amplitude op manieren die het signaal-ruisverhouding kunnen verminderen. Deze beperkingen zijn bijzonder kritiek voor opkomende multi-locus transcraniële magnetische stimulatie (mTMS) protocollen, die emergente aandrijvingen induceren bij samengestelde frequenties in overlappende corticale regio's; het monitoren van alleen connectiviteit op dezelfde frequentie in dergelijke scenario's riskeert het missen van de daadwerkelijke dynamische kenmerken die stimulatie-effecten bemiddelen.

Methodologie
De auteurs introduceren een algemene familie van antisymmetrische cross-polyspectrale indices die zijn ontworpen om harmonische afhankelijkheden (specifiek $1:n$ interacties) te kwantificeren, terwijl ze robuust blijven tegen instantane lineaire menging.

1. Theoretische Formulering:

   • Het kader begint met de definitie van een trispectrum voor twee tijdreeksen $x(t)$ en $y(t)$, met focus op kubische afhankelijkheden ($1:3$ koppeling waarbij $f_1=f_2=f_3=f$ en $f_4=3f$).
   • Onder een lineair mengmodel ($x = \sum a_i s_i$, $y = \sum b_i s_i$) bevat het standaard trispectrum bijdragen van zowel echte bronnen als lineaire mengartefacten.
   • Om echte koppeling te isoleren, definiëren de auteurs een antisymmetrische cross-trispectrum door de gepermuteerde versie van de schatter af te trekken: $T[x|xx|y] = T_{xxxy} - T_{yxxx}$.
   • Een theoretische afleiding toont aan dat onder de nulhypothese van onafhankelijke bronnen die lineair gemengd zijn op de sensoren, de symmetrische mengbijdragen exact opheffen, waardoor alleen termen overblijven die geassocieerd zijn met echte niet-lineaire koppeling vanuit dezelfde bron.
   • Deze constructie wordt genormaliseerd met behulp van amplitude-normen van de vierde orde om de Antisymmetrische Cross-Tricoherentie (ACT) te verkrijgen, begrensd binnen $[0, 1]$.
   • Het kader wordt gegeneraliseerd tot een willekeurige orde $m$ (Antisymmetrische Cross-Polyspectrum), waarbij de Antisymmetrische Cross-Polycoherentie (ACP) wordt gedefinieerd.
   • De auteurs tonen aan dat voor $m=2$ deze formulering reduceert tot het imaginaire deel van de coherentie (ImCoh), en voor $m=3$ correspondeert het met de antisymmetrische cross-bicoherentie.

2. Validatie-experimenten:

   • Synthetische Data: Een generatief model werd gebruikt om tijdreeksen te simuleren die een mix bevatten van echte kubische koppeling en artefactuele instantane lineaire menging, gecontroleerd door een parameter $\gamma$. De prestaties van ACT werden vergeleken met niet-antisymmetrische trispectrale maatstaven (CT1 en CT2).
   • Echte EEG-data: De indices werden toegepast op een single-subject rust-EEG dataset (61 kanalen). Statistische significantie werd beoordeeld met behulp van een parametrische test op basis van surrogate data (willekeurig permuteren van temporele segmenten) om p-waarden te genereren.

Belangrijkste Resultaten

• Synthetische Data:
  • In regimes gedomineerd door echte kubische koppeling ($\gamma \approx 0$) toonden zowel ACT als standaard trispectrale maatstaven verhoogde waarden, wat de gevoeligheid voor echte niet-lineaire afhankelijkheden bevestigt.
  • In regimes gedomineerd door instantane lineaire menging ($\gamma \approx 1$) werden standaard maatstaven (CT1, CT2) systematisch opgeblazen door nul-lag correlaties. Daarentegen bleven ACT-waarden dicht bij de basislijn, wat robuuste immuniteit tegen schijnbare mengartefacten aantoont.
  • In intermediaire regimes vertoonden standaard maatstaven complexe, verwarde gedragingen (bijv. U-vormige afhankelijkheid van $\gamma$), terwijl ACT specificiteit behield voor echte koppeling.
• Echte EEG-data:
  • Standaard trispectrale kaarten (CT1, CT2) vertoonden uitgesproken, ruimtelijk beperkte hoog-waarde vlekken gecentreerd rond zaad-sensoren, een topografie die consistent is met volumgeleiding in plaats van lange-afstandsinteractie.
  • ACT-kaarten vertoonden deze korte-afstands "bloeiingen" niet. In plaats daarvan onthulden ze bescheiden, ruimtelijk verspreide toenames richting frontale kanalen (voor occipitale zaden) en duidelijke lange-afstands interhemisferische patronen (voor pariëtale zaden).
  • Statistische toetsing bevestigde dat hoewel sommige lokale pieken in standaard maatstaven de drempelwaarde overleefden (waarschijnlijk artefactueel), de lange-afstands patronen geïdentificeerd door ACT statistisch robuust waren na surrogate drempelwaarde.
  • Een directe vergelijking met ImCoh (lineaire connectiviteit) toonde extreem lage correlatiecoëfficiënten ($\rho \approx 0.04-0.06$), wat aangeeft dat ACT grotendeels onafhankelijke kenmerken van neuronale connectiviteit vastlegt die niet worden verklaard door lineaire fase-synchronisatie.

Betekenis en Aanspraken
Het artikel claimt het eerste robuuste kader te bieden voor het identificeren van echte hoog-orde cross-frequency interacties tussen meerdere tijdreeksen dat intrinsiek resistent is tegen lineaire mengartefacten. De auteurs stellen dat deze indices bijzonder actueel zijn voor de analyse van multi-locus TMS (mTMS), waarbij stimulatieprotocollen samengestelde spectrale aandrijvingen induceren ($f_{overlap} = \sum f_i$) die niet kunnen worden gemonitord via connectiviteit op dezelfde frequentie alleen.

De studie dient als bewijs van concept, en demonstreert dat:

1. Antisymmetrisatie effectief volumgeleiding artefacten opheft in spectrale analyse van hogere orde.
2. Echte afhankelijkheden van hogere orde op sensorniveau kunnen worden gedetecteerd in single-subject EEG-data.
3. Deze metrieken een principiële basis bieden voor het leiden van gepersonaliseerde, multi-site, multi-frequentie stimulatiestrategieën.

De auteurs erkennen dat hoewel de methode robuust is, toekomstig werk vereist is om deze indices te valideren over grotere subjectcohorten om te bepalen of de waargenomen koppelingspatronen algemene principes van corticale organisatie of individuele kenmerken vertegenwoordigen. Zij merken ook de noodzaak op van een rigoureuze statistische karakterisering om de datalengte en signaal-ruisvereisten te optimaliseren voor potentiële real-time toepassingen in adaptieve neuromodulatie.