$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Чтобы одновременно записывать до 12 областей мозга SBN, мы создали настраиваемый мультиволоконный массив по вышеуказанному протоколу. После восстановления после операции мыши были зафиксированы во время социального взаимодействия (рисунок 5). Во всех экспериментах самцам мышей с CD1 разрешалось взаимодействовать с самцами мышей Balb/c в их домашней клетке. Мы использовали мышей Vgat-Cre для экспрессии jRCaMP1b40 в ГАМКергических нейронах в зависимости от Cre в 11 областях мозга SBN, GCaMP6f39 в не-ГАМКергических клетках— 45 в тех же областях мозга, и GRAB-DA1m33 в Nucleus accumbens42 (Таблица 1). Гистологическое подтверждение таргетирования клетчаткой можно найти в наших предыдущихпубликациях 42,43. Здесь мы стремились проверить эффективность различных методов коррекции артефактов движения. В частности, мы сравнили эффективность использования: (1) только медианного фильтра; (2) Регрессия лассо с регуляризацией L1 с использованием только одноволокнистых изосбестических сигналов (метод с одним предиктором)11 и (3) регрессия лассо с использованием изосбестических сигналов всех волокон в качестве предикторов (мультипредикторный подход)41,42 (Рисунок 4A). Критически важно, что методы на основе лассо корректируют исходный сигнал, используя независимый от активности флуоресцентный сигнал, генерируемый в ответ на светодиод с изосбестической частотой (415 нм), тогда как сглаживание медианного фильтра не учитывает независимый от активности сигнал.
Мы оценивали коррекцию движения с помощью двух разных показателей: абсолютного высокочастотного шума и джиттера. Высокочастотный шум определялся как стандартное отклонение первой производной сигнала, отражающее величину резких колебаний сигнала. Джиттер был определен как коэффициент вариации производной сигнала, который отражает величину изменчивости сигнала от кадра к кадру относительно среднего значения производной сигнала. Мы обнаружили, что оба подхода, основанные на Лассо (однопредиктор и мультипредиктор) значительно превосходят медианный фильтр по снижению как высокочастотного шума, так и джеттера, при этом коррекция с мультипредикторами немного превосходит коррекцию одного предиктора (рисунок 4C–D; Высокочастотный шум, n=414 для 23 сигналов в 18 сессиях записи: тест Крускаля-Уоллиса, H = 501.9420, p < 0.001; пост-hoc тесты Mann-Whitney U с коррекцией Бонферрони: сглаженная против односайтовой регрессии ***p < 0,001, сглаженная против многосайтовой регрессии***p < 0,001, односайтовая против многосайтовой **p = 0,0073; Джиттер: тест Крускал-Уоллис, H = 495.6685, стр. < 0.001; пост-hoc тесты Mann-Whitney U с коррекцией Бонферрони: сглаженная против односайтовой регрессии ***p < 0,001, сглаженная против многоточковой регрессии***p < 0,001, одноместная против многосайтовой ***p = 0,0006). В частности, мультипредикторный подход (с использованием всех волоконных изосбестиков) показал преимущества для сигналов красного канала (RCaMP): 5 из 12 красных сигналов показали значительно меньший остаточный шум по сравнению с коррекцией с одним предиктором (n = 18 записей/регион). Для джиттера BNST (t-тест, t = 5,42, ***p < 0,001), POA (t-тест, t = 2,68, *p = 0,011), AH (t-тест, t = 2,66, *p = 0,012) и PAG (тест Mann-Whitney U, U = 279, ***p < 0,001) показали превосходство мультипредикторов. Для высокочастотного шума BNST (тест Mann-Whitney U, U = 285, ***p < 0,001), POA (тест Mann-Whitney U, U = 292, ***p < 0,001), AH (t-тест, t = 7,10, ***p < 0,001), PAG (t-тест, t = 13,72, ***p < 0,001) и LHb (t-тест, t = 5,47, ***p < 0,001) показали превосходство по мультипредикторам. Сигналы зелёного канала (GCaMP) не показали значимых различий между двумя подходами Lasso во всех 12 регионах по метрикам джиттера или высокочастотного шума (все p>0.05, полный список статистических тестов см. Таблицу 3). Это говорит о том, что для двухцветной фотометрии включение кросс-канальной изосбестовой информации улучшает коррекцию движения для красных сигналов, вероятно, потому что изосбестики зелёного канала обеспечивают более чистую информацию о движении. Эти данные показывают, что коррекция артефактов движения с использованием методов регрессии на основе изосбеста значительно превосходит простые методы фильтрации, при этом многоканальная регрессия даёт особые преимущества для записи с двумя цветами.

Рисунок 1. Запись из сети социального поведения с использованием многопрофильной двухцветной фотометрии волоконного соединения. (A) Схематическое представление того, как построить мультиволоконный массив (multifiber). (B) Схематическое представление рабочего процесса хирургии. (C) Схематическое представление двухцветных экспериментов по кальциевой визуализации. (D) Пример использования изосбестического сигнала для вычисления сигнала, скорректированного с изменением движения. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 2. Изображения ключевых процедур изготовления и хирургии. (A) Изготовленный патчкорд, готовый к записи. (B) Мультиволоконная решетка, помещённая в индивидуальный стереотаксический держатель с сверлом размером 0,25 мм x 5,5 мм, расположенной над сквозным отверстием для волокна. (C) Многоволоконный массив с очищенными сквозными отверстиями и приклеенными шестигранными гайками. (D) Разрезание волокна с помощью волокнистого приписчика. (E) Многоволоконный массив со всеми волокнами, приклеенными суперклеем и размещёнными на индивидуальную полировочную шайбу. (F) Увеличенный вид того, как выравнивать и прикручивать разъём патчкорда к мультиволоконному массиву. (G) Полностью собранный мультиволоконный массив с приклеенными штифтами. (H) Мышь, прикреплённая к стереотаксической раме с разрезом на коже головы, чтобы обнажить череп. (I) Изготовленная многоволоконная решетка, прикреплённая к стереотактическому манипулятору для управления мультифибром во время операции. (J) Увеличенный вид многоволоконного массива, опускаемого в мозг. (K) Многоволокнистый массив полностью опущен и прикреплен к черепу зубным цементом. (L) Мультиволоконный массив на животном, надёжно подключённый к патчкорду перед записью. (M) Изображение животного в клетке с подключением патчкорда до и/или во время записи. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 3. Спектральное разделение. Примеры примеров демонстрируют минимальную спектральную перекрёстную связь между кальциевыми индикаторами. Сигналы GCaMP6f (синий) демонстрируют устойчивую активацию во время возбуждения 470 нм, но незначительную реакцию на возбуждение 560 нм. Сигналы jRCaMP1b (красный) демонстрируют обратный рисунок с сильной активацией от 560 нм, но минимальным откликом на освещение 470 нм. (A) Одновременно записывались сигналы GCaMP6f и jRCaMP1b с освещением только 470 нм. (B) То же, что и (a), но с освещением только 560 нм. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 4. Сравнение различных методов мультиволоконной обработки сигналов. (A) Схема, изображающая многосайтовую регрессионную коррекцию фотометрического сигнала. Этот подход включает последовательные этапы обработки, при которых все доступные 415-нм изосбестические сигналы линейно соотносятся с каждым кальциезависимым сигналом. Предсказанные компоненты движения вычитаются из этого сигнала для получения итогового показания активности. z-баллированные сигналы ΔF/F от двух сигналов с выраженным наличием артефактов движения. (B) Сверху вниз: изосбестовый сигнал из того же волоконного/изображенного канала, отображающий артефакты движения, сглаженный только RCaMP (слева) или GCaMP (справа) (медианно-фильтрованный, зелёный), тот же сигнал с регрессионной корректировкой по оптоволокном/каналу изображения (циан), многоволоконный регрессионный сигнал (синий). (C) Слева: средний балл высокочастотного шума для каждого записанного сигнала (n = 21 сессия). Высокочастотный шум вычисляется как стандартное отклонение первой производной каждого сигнала. Оценки шума показаны для 1) сигнала с корректированной базовой линией, но некорректированного сигнала (серый), 2) сигнала с коррекцией базовой линии и медианной фильтрации (розовый), 3) сигнала с корректированной регрессией базовой и одноволокнистой линии с сглаживанием Савицкого-Голея (циан), и 4) сигнала с коррекцией регрессии по базовой и многоволоконной линии с сглаживанием Савицкого-Голея (синий). Центр: средний процент шума, удалённый только с помощью сглаживания, одноволоконной или многоволоконной регрессии, агрегированных по всем 23 сигналам. Процент удалённого шума рассчитывается как доля высокочастотного шума, возникающего в результате коррекции относительно общего уровня шума в сыром сигнале, скорректированном на базовую линию. Данные отображаются как среднее ± SEM (тест Крускал-Уоллис, H = 501.9420, p < 0.001; пост-хок тесты Mann-Whitney U с коррекцией Бонферрони: сглаженная против односайтовой регрессии ***p < 0.001, сглаженная против многосайтовой регрессии***p < 0.001, одноместная против многоточочной **p = 0.0073). Справа: Процент сигналов, при которых высокочастотный шум статистически отличается после коррекции на основе регрессии с одним волоконным и мультиволоконным режимом, разделённый по типу сигнала. Для каждого из 23 сигналов уровни шума после коррекции на одном участке сравнивались с уровнями шума после коррекции на несколько точек с помощью независимых t-тестов (если они распределены нормально), либо тестов Mann-Whitney U (непараметрическая альтернатива). Синие сегменты указывают на сигналы, где коррекция на несколько точек приводила к значительно меньшему шуму, чем коррекция на одном сайте (p < 0,05). Серые сегменты указывают на сигналы без значительной разницы между методами коррекции. Циановые сегменты указывают на сигналы, где коррекция на одном месте приводила к значительно снижению шума. E: Vgat- GCaMP6f. ДА: ХВАТАЙ-ДА. I: Vgat+ jRCaMP1b. (D) То же самое, что и (C), но с изображением вариативности кадра к кадру (джиттер), рассчитанной как производная коэффициента вариации каждого сигнала (тест Крускаля-Уоллиса, H = 495.6685, p < 0.001; пост-хок тесты Mann-Whitney U с коррекцией Бонферрони: сглаженная против односайтовой регрессии ***p < 0.001, сглаженная против многоточочной регрессии***p < 0.001, одноточковая против многоточочной ***p = 0.0006). Полный список статистических тестов и результатов см. Таблицу 3. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.

Рисунок 5. Многопрофильная демонстрация двухцветной фотометрии волокон. Репрезентативные сигналы от нескольких типов индикаторов, фиксируемые одновременно во время социального поведения (тест резидент-незваный гость, эпоха входа незваного гость). Следы показывают ненормализованные и исходно-вычитаемые, но не корректированные артефактами движения GCaMP6f (Cre-отрицательные нейроны), jRCaMP1b (ГАМКергические Vgat-Cre+ нейроны) и GRAB-DA1m (дофаминовый сенсор в ядре accumbens) на сайтах социальных поведенческих сетей у одного и того же животного. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этой фигуры.
Таблица 1. Координаты вируса и длина волокон для каждой целевой области мозга социальной поведенческой сети. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Таблица 2. Сравнение методов коррекции артефактов движения для многопрофильной фотометрии волокон . Основные принципы, оптимальные сценарии использования, преимущества и ограничения следующих методов коррекции артефактов движения: простая линейная регрессия, регрессия Лассо с регрессией L1 (с использованием мультипредикторного подхода), байесовское генеративное моделирование, частотное масштабирование и итеративно перевзвешенные наименьшие квадраты. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Таблица 3. Детали и результаты статистических тестов, использованных на рисунке 4. Пожалуйста, нажмите здесь, чтобы скачать эту таблицу.
Дополнительный файл 1. Файл 3D-печати для разъёма патчкорда.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 2. 3D печать для разъёма патчкорда без сквозных отверстий.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 3. 3D печать для стереотактического сверлящего аксессуара.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 4. Индивидуальный чертёж шайбы для полировки.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 5. 3D печатный файл для многосайтового фотометрического массива, ориентированного на 12 сайтов SBN.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 6. 3D печатный файл для многосайтового фотометрического массива без сквозных отверстий, направленных на сайты.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 7. Конфигурационный файл DAQ для программного обеспечения DAQ.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 8. Файл визуального кодирования для фотометрической системы.Пожалуйста, нажмите здесь, чтобы скачать этот файл.
Дополнительный файл 9. Код на Python, используемый в Протоколе Step 4 и для сравнений, корректирующих артефакты движения. Пожалуйста, нажмите здесь, чтобы скачать этот файл.