Najpierw wykreśliliśmy średnie potencjały wywołane (ERP) dla każdego bodźca słuchowego u marmosetów (Rysunek 2). Potencjał wywołany słuchowo (AEP) był wyraźny w warunku Noise, co odzwierciedlało wyraźny początek bodźców (patrz Rysunek 1D). Aby porównać uśrednione ERP pomiędzy typami wołań a bodźcami szumu, zastosowaliśmy jednoczynnikową analizę wariancji (ANOVA), przyjmując bodźce jako czynnik międzygrupowy dla odpowiedzi w punkcie Cz. Stwierdzono istotny główny efekt bodźców na aktywność w punkcie Cz odpowiednio w przedziałach 13-18 ms, 28-36 ms oraz 45-88 ms po wystąpieniu bodźca (p < 0.05). Analiza porównań wielokrotnych post hoc metodą Tukeya wykazała, że różnica występowała pomiędzy szumem a pozostałymi wołaniami (p < 0.05), natomiast nie stwierdzono różnic pomiędzy typami wołań marmosetów. Wynik ten sugeruje, że różnice w aktywności mózgu w zależności od typu wołania nie mogły zostać zaobserwowane wyłącznie na podstawie potencjałów wywołanych.
Następnie przeprowadzono analizę czasowo-częstotliwościową dla każdego badanego. Rysunek 3 przedstawia przykład map czasowo-częstotliwościowych dla sygnału Tsik-string uzyskanego od badanego R (starszego, Rysunek 3A) oraz badanego Y (młodszego, Rysunek 3B). Stwierdzono, że moc widmowa związana z bodźcem wzrosła przy niższej częstotliwości wynoszącej około 20-50 Hz bezpośrednio po wystąpieniu bodźca. Reakcje te były wyraźnie widoczne w elektrodzie Cz. W przeciwieństwie do tego, moc w zakresie gamma (powyżej 30 Hz) spadła po wystąpieniu bodźca w porównaniu z okresem bazowym. Ponadto spadek ten trwał przez 1 s. U osób starszych powyżej 8 lat nie zaobserwowano spadku mocy związanej z bodźcem. W tych przykładach osoby starsze wykazały silniejszą pierwotną odpowiedź na sygnał w obszarze korony (Cz), podczas gdy u osób młodszych zaobserwowano trwały spadek aktywności w pasmie γ podczas prezentacji sygnału. Wyniki sugerują istnienie różnic w odpowiedziach pierwotnych i trwałych w zależności od wieku badanego.
Na koniec zbadano związek między wiekiem badanych a mocą perturbacji widmowej związanej z zdarzeniem (ERSP) w początkowej odpowiedzi przejściowej (Rysunek 4A) i odpowiedzi podtrzymywanej (Rysunek 4B). Przeprowadzono dwuczynnikową analizę wariancji (ANOVA) z typem bodźca jako czynnikiem wewnątrzgrupowym oraz wiekiem jako czynnikiem międzygrupowym, aby określić wpływ rodzaju bodźca słuchowego i wieku badanych na aktywność EEG. W początkowych odpowiedziach przejściowych w punkcie Fz stwierdzono istotne efekty główne dla typu bodźca (F (3,24) = 9.020, p < 0.001) oraz wieku (F (8,24) = 3.934, p = 0.004). Nie zaobserwowano jednak istotnej interakcji między typem bodźca a wiekiem (p = 0.144). W odpowiedziach przejściowych w punkcie Cz wystąpił istotny efekt główny typu bodźca (F (3,24) = 8.533, p < 0.001), lecz nie stwierdzono wpływu wieku (F (8,24) = 2.215, p = 0.073) ani interakcji (p = 0.228). W przypadku odpowiedzi podtrzymywanych w punkcie Fz odnotowano istotne efekty główne zarówno typu bodźca, jak i wieku (odpowiednio F (3,24) = 9.020, p < 0.001; F (8,24) = 3.934, p = 0.004). Nie zaobserwowano istotnej interakcji (P = 0.144). Odpowiedzi podtrzymywane w punkcie Cz wykazały istotny efekt główny typu bodźca (F (3,24) = 8.533, p < 0.001), ale nie stwierdzono efektu głównego wieku (F(8,24) = 2.215, p = 0.073) ani interakcji (p = 0.228). Wyniki te sugerują, że w środkowym obszarze czołowym (Fz) początkowe odpowiedzi przejściowe na bodźce w postaci wołania i szumu znacznie różniły się wraz z wiekiem, a odpowiedzi podtrzymywane były tłumione w młodszych grupach wiekowych. Może to odzwierciedlać funkcjonalne dojrzewanie obszaru czołowego.
Poprzednie badania neurofizjologiczne zgłaszały odpowiedzi neuronalne w pierwotnej korze słuchowej podczas wokalizacji u marmosetów9,38. Ponadto ponad połowa tych neuronów wykazuje odpowiedź hamującą, która utrzymuje się podczas wokalizacji9,38. Co więcej, wcześniejsze badania elektrofizjologiczne u naczelnych i ludzi wykazały, że aktywność wysokiego pasma gamma w potencjale polowym (LFP) i ECoG dobrze koreluje z częstotliwością wyładowań neuronów39,40,41,42. EEG skóry głowy jest spatiotempelalnie wygładzoną wersją LFP, zintegrowaną nad obszarem o powierzchni 1 cm2 lub większej43. Chociaż komponent high-gamma w EEG wykazuje niższą korelację z częstotliwością wyładowań niż LFP i ECoG, uważa się, że koduje on sygnał wyjściowy jako zintegrowany zakres kilku centymetrów. Zatem trwały spadek mocy pasma gamma zaobserwowany w punkcie Cz w naszych eksperymentach może odzwierciedlać aktywność klastrów neuronalnych wykazujących stłumioną aktywność podczas emisji wokalizacji, co występuje w korach słuchowych. W przeciwieństwie do tego, poprzednie badania elektrofizjologiczne zgłosiły, że więcej neuronów w korze czołowej, głównie w korze przedmotorycznej, wykazuje odpowiedzi pobudzające podczas wokalizacji sygnałów13. Co ciekawe, trwałą aktywność hamującą zaobserwowano w naszych eksperymentach nawet wtedy, gdy Fz znajdował się w obszarze czołowym, mimo że zaobserwowano odrębne mechanizmy neuronalne pomiędzy percepcją a produkcją wokalną. Niedawne badanie fMRI zidentyfikowało ponadto kilka podregionów w korze czołowej, w tym przednią korę obręczy oraz korę przedmotoryczną, jako „plamy wokalne” (vocal patches), które reagują na gatunkowe sygnały u marmosetów19. Nasze wyniki odzwierciedlają ogólną aktywność mózgu w tych obszarach.
W niniejszym eksperymencie do badania wykorzystano wyłącznie elektrody w pozycjach linii środkowej (Fz, Cz, Pz oraz Oz); w związku z tym nie możemy odnieść się do różnic między prawostronną a lewostronną aktywnością EEG w procesach słuchowych. W przyszłości konieczne będzie zbadanie lateralizacji aktywności neuronalnej leżącej u podstaw przetwarzania wokalnego.
Wcześniejsze badania wykazały, że niska aktywność gamma w LFP i ECoG jest generowana przez wejścia synaptyczne do komórek piramidalnych. Zatem wysoka aktywność gamma odzwierciedla komponent sygnału bliższy wyjściu, podczas gdy niska aktywność gamma odzwierciedla komponenty bliższe wejściu39. W niniejszym badaniu zaobserwowaliśmy przejściową aktywność w pasmach beta i niskiej gamma bezpośrednio po ekspozycji na wokalizację. Reakcje te mogą odzwierciedlać sygnały wejścia sensorycznego do kory. Zaletą naszej metody jest możliwość rejestracji aktywności mózgu z różnych populacji neuronalnych jako dynamicznych zmian z wysoką rozdzielczością czasową. Z naszej wiedzy jest to pierwsze badanie ujawniające, w jaki sposób sygnały EEG ze skóry głowy zmieniają się podczas percepcji gatunkowych wokalizacji u marmosetów. Przedstawione wyniki dostarczają nowych informacji na temat integracji reprezentacji neuronalnych dzięki rejestracji z szerokiego zakresu obszarów mózgu.

Rysunek 1: Układ eksperymentalny. (A) Przykładowy obraz zwierzęcia podczas rejestracji. Marmoseta siedzi w krześle, a jego głowa jest przymocowana do krzesła za pomocą maski. Jama ustna pozostaje otwarta, aby ułatwić oddychanie i picie płynów będących nagrodą, a elektrody zostały przymocowane do górnej części głowy. (B) Wyposażenie: Bodźce dźwiękowe są prezentowane przez głośnik. Zainstalowano również wzmacniacz, moduł wejściowy elektrod oraz kamerę monitorującą. (C) Lokalizacja elektrod: Elektrody rozmieszczono w punktach Fz, Cz, Pz, Oz, A1 i A2 zgodnie z międzynarodowym systemem 10-20. Lokalizację elektrod określono, przyjmując inion, nasion oraz obustronne punkty przeduszne jako punkty orientacyjne anatomiczne. Jako elektrodę masy wykorzystano elektrodę C3 lub C4. (D) Sonogramy (lewy panel) i spektrogramy (lewy panel) dla wszystkich bodźców dźwiękowych. Wołanie Phee to pojedyncze długie wołanie trwające mniej niż 2 s. Wołanie Tsik-Ek jest wołaniem złożonym, składającym się z Tsik, po którym następuje Ek. Dwa zestawy wołania złożonego prezentowano przez około 1 s. Wołanie Tsik-string to powtarzające się wołanie Tsik5, a w bodźcu zaprezentowano cztery wołania Tsik przez około 1 s. Jako bodziec niebędący wołaniem zastosowano sygnał szumu białego wygenerowany przez specjalny skrypt, trwający około 1 s. Opóźnienia początku dźwięku sprawdzono wizualnie w cyfrowym edytorze audio. Czerwone strzałki wskazują poszczególne opóźnienia: wołanie Phee 49 ms, wołanie Tsik-Ek 35 ms, wołanie Tsik-string 16 ms oraz szum biały 0 ms. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Średnie potencjały wywołane w odpowiedzi na wołania Phee, Tsik-Ek i Tsik-string oraz szum. (n = 9) Aktywność została dopasowana do początku każdego wołania lub białego szumu. Czarna pozioma strzałka wskazuje okres prezentacji bodźca. Wołanie Phee trwało około 2 s, pozostałe trwały 1 s. Czerwone poziome linie przy Cz wskazują okresy z istotnymi różnicami między bodźcami słuchowymi. Wszystkie one występowały między Szumem a pozostałymi wołaniami; nie stwierdzono różnic między samymi wołaniami. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3: Przykład map czasowo-częstotliwościowych dla wezwań. (A) Mapa ERSP dla wezwania Tsik-string u badanego R (wiek 145 miesięcy); (B) Mapa EPSP dla tego samego wezwania Tsik-string u badanego Y (wiek 23 miesiące). Panele lewy i prawy przedstawiają dane zarejestrowane odpowiednio z elektrod Fz i Cz. Czerwone linie pionowe wskazują moment rozpoczęcia pliku dźwiękowego, a nie początek samego wezwania. Skrót: EPSP = zaburzenie widma powiązane z zdarzeniem. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 4: Zależność między średnią mocą powiązaną z bodźcem a wiekiem badanych. (A) Zależność między średnią mocą powiązaną z bodźcem w pasmach α i β a wiekiem badanych. Średnią moc ERSP obliczono dla zakresu 8-29 Hz w okresie 1-150 ms od początku każdego bodźca w porównaniu do wartości w okresie bazowym (-200 do 0 ms przed początkiem bodźca). Panele lewy i prawy przedstawiają dane zarejestrowane odpowiednio z elektrod Fz i Cz. (B) Zależność między średnią mocą powiązaną z bodźcem w pasmie γ a wiekiem badanych. Średnią moc ERSP obliczono dla zakresu 30-100 Hz w okresie 151-950 ms od początku każdego bodźca w porównaniu do wartości w okresie bazowym (-200 do 0 ms przed początkiem bodźca). Zatem ujemne wartości ERSP wskazują na spadek mocy w porównaniu do stanu przed prezentacją bodźca. Skrót: ERSP = zaburzenie widma powiązane z bodźcem (event-related spectral perturbation). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Plik uzupełniający 1: Plik zip zawierający cztery pliki audio wykorzystane w eksperymentach.Plik PH.wav zawiera bodziec z jednym długim głosem Phee; TE.wav zawiera dwa głosy Tsik-Ek z odstępem czasowym między nimi; MB.wav zawiera bodziec w postaci głosu Tsik-string z czterema następującymi po sobie głosami Tsik (zwanymi również głosami mobbingowymi); NO.wav zawiera wygenerowany programowo biały szum gaussowski. Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Plik zip zawierający kod do postprocessingu. Kliknij tutaj, aby pobrać ten plik.