$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
В этой статье мы предлагаем метод надежного обнаружения косоглазия с высоким временным разрешением с помощью DeepLabCut. Мы оптимизировали параметры обучения и дали оценку сильных и слабых сторон этого метода (рис. 1).
После обучения наших моделей мы убедились, что они способны правильно оценить верхнюю и нижнюю точки века (рис. 2), которые служат координатными точками для евклидовой меры расстояния. Евклидово расстояние определяется как средняя длина расстояний между двумя верхними и нижними точками глаза. Наша модель смогла обнаружить случаи отсутствия косоглазия (рисунок 2A) и косоглазия (рисунок 2B). Синими точками обозначены точки, используемые для определения евклидова расстояния для каждого кадра. Зеленая, желтая, оранжевая и фиолетовая точки были использованы, чтобы помочь модели правильно оценить евклидово расстояние и уменьшить значение вероятности, когда голова находится в неоптимальном положении (т.е. с учетом движения головы и изменений положения между сеансами). Затем мы проверили точность модели с помощью ряда различных методов.
Чтобы проверить идеальное количество кадров, используемых для модели, мы обучили и протестировали четыре модели с различным размером выборки кадров (рис. 3). Сначала мы сравнили значения среднеквадратичной ошибки (RMSE) между тестовыми и обучающими данными, чтобы проверить, насколько хорошо модели могут точно прогнозировать тестовые данные, на которых они не были обучены. Это сравнение показало, что вариативность между точками, помеченными вручную, и точками, помеченными моделью, выровнялась после 300 кадров. Эта тенденция коррелировала со средними значениями вероятности, которые также выровнялись после 300 помеченных кадров. Мы использовали эти значения вероятности для фильтрации точек, которые были меньше 0,92. Эти значения вероятности показывают, насколько модель уверена в том, что данная точка была правильно помечена на основе обучающих данных. Мы усреднили эти значения для точек, которые вносят вклад в метрику евклидова расстояния, чтобы проверить, насколько хорошо модели работают относительно друг друга. Несмотря на то, что существенной разницы между 300 и 400 кадрами не было, мы использовали 400 кадров, потому что их среднее значение превышало значение правдоподобия 0,95, которое приближается к нашему пороговому значению для ручной фильтрации и соответствует пороговому значению, используемому в аналогичных моделях для оценки позы16.
Еще одним способом проверки точности модели была матрица несоответствий, сравнивающая вручную аннотированные кадры с кадрами, помеченными DLC. Двое слепых вручную аннотировали 300 кадров одного и того же глаза в восьми видеороликах. Мы использовали эти данные для построения матрицы неточностей для оценки истинных и ложных положительных и отрицательных результатов (рис. 4), где данные, полученные вручную, использовались в качестве исходной истины. Для DLC положительное значение косоглазия регистрировалось, когда евклидово расстояние было записано как менее 75 пикселей (т.е. животное косоглазо), и отрицательное значение регистрировалось при значениях больше 75 пикселей (т.е. животное не щурится). Мы обнаружили положительное прогностическое значение, равное 96,96%, что представляет собой процент времени, в течение которого модель точно предсказывает косоглазие относительно косоглазия, аннотированного вручную. Мы обнаружили отрицательную прогностическую ценность в 99,66%, которая представляет собой процент времени, в течение которого модель точно предсказывает отсутствие косоглазия по сравнению с аннотированным вручную косоглазием. Они показывают соотношение отрицательных и положительных значений, которые были правильно помечены. Мы также обнаружили истинно положительный уровень 98,1% и истинно отрицательный уровень 99,46%, которые представляют собой точное предсказание модели положительных и отрицательных значений относительно всех положительных и отрицательных значений соответственно. Наш коэффициент корреляции Мэтьюса, или MCC, составил 93,8%, что указывает на коэффициент корреляции между наблюдаемыми и прогнозируемыми значениями.
Убедившись, что наша модель надежно отслеживает косоглазие, мы сравнили этот метод DLC с ранее опубликованным методом отслеживания косоглазия с использованием доклинического набора данных о мигрени14. Мы будем называть этот другой метод «моделью площадного косоглазия (ASM)», потому что он был разработан с использованием области открытого глаза в качестве непрерывной переменной измерения косоглазия14. Модель косоглазия использует обученное программное обеспечение для распознавания лиц в сочетании с пользовательским скриптом MATLAB для анализа средней площади глаза в пикселях, исключая кадры с коэффициентом ошибок отслеживания >15%14. Одним из основных ограничений является то, что «ASM» не является открытым исходным кодом и, следовательно, не широко доступен. DLC обеспечивает повышенную оптимизацию и адаптивность, не требуя значительной покупки программного и аппаратного обеспечения.
Мы использовали набор данных из 10 самок и 10 самцов мышей CD1. Экспериментально все животные были акклиматизированы в щадящих фиксаторах в течение 30 минут в общей сложности за 3 дня до начала записи. Каждое животное записывалось в течение 5 минут исходного уровня, а затем 5 минут для записи лечения. Во время сеансов лечения животных вводили либо PBS (носитель), либо 0,1 мг/кг CGRP (лечение) внутрибрюшинно, чтобы вызвать мигреньподобное состояние. Данные собирались в хорошо освещенном помещении с помощью камер, оснащенных инфракрасным светом для освещения лица, обеспечивая точное обнаружение ориентиров. Инфракрасная камера включала в себя объектив Kowa LM35JC 2/3" 35 mm F1.6 с ручной диафрагмой C-mount и фокусным расстоянием 254 мм и соответствующим образом отрегулированной диафрагмой. После того, как мы собрали данные, мы использовали ASM и DLC для их анализа. Поскольку ручная оценка традиционно используется в полевых условиях для количественной оценки гримасы лица, а косоглазие является одним из компонентов гримасы лица14, мы также сравнили наши данные с данными, полученными вручную.
Основываясь на предыдущих выводах о том, что периферическая инъекция CGRP вызывает косоглазие у мышей, мы ожидали наблюдать значительные различия в реакции косоглазия между носителем и лечением CGRP 6,14. Мы сравнили методы ASM, ручные методы и методы DLC и обнаружили, что наша модель надежно обнаруживает фенотип косоглазия, как и ручные методы и методы ASM (рис. 5). Важно отметить, что модель ASM использовалась для оценки боли и косоглазия, вызванных CGRP. В этом исследовании Rea et al. сравнили реакцию косоглазия после CGRP с реакцией на косоглазие после инъекции формалина в заднюю лапу в качестве «более традиционного» метода индукции боли14. Кроме того, хорошо известно, что CGRP вызывает гиперчувствительность к прикосновениям у мышей при использовании метода фон Фрея 3,17. В соответствии с полевыми условиями, мы нормализовали среднее косоглазие во время сеанса лечения до 5 минут до начала обработки для каждого животного и сравнили PBS (n = 10) с животными, получавшими CGRP (n = 11). Статистический анализ групп PBS в сравнении с группами, получавшими CGRP, выглядит следующим образом. Мы обнаружили, что животные, получавшие CGRP, демонстрировали уменьшение средней площади пикселя при использовании метода отслеживания с использованием площадного косоглазия (p = 0,012, рисунок 5A) и уменьшение евклидова расстояния при ручной оценке (p = 0,0007, рисунок 5B) и использовании нашей модели DLC (p = 0,007, рисунок 5C). Когда мы сравнивали каждый метод с течением времени на одном репрезентативном животном, наблюдалась та же закономерность (рис. 5). У этого животного наблюдался очень четкий фенотип косоглазия в ответ на лечение CGRP, но не на PBS. Все модели смогли обнаружить эти различия, но наиболее четко данные были представлены в нашей модели DLC (рис. 5). Точные и достоверные метрики особенно важны, когда данные должны анализироваться с более высоким разрешением, где усреднение не указывает на полное поведенческое считывание (например, активность мозга). Метод DLC для обнаружения косоглазия у мышей позволяет нам собирать данные в миллисекундном масштабе времени и привязывать их к измерениям активности мозга (например, потенциалов локального поля), которая происходит в миллисекундном масштабе времени. Затем мы можем использовать эту технику для построения более надежного профиля состояния мозга, свидетельствующего о спонтанной боли в контексте мигрени и других сложных расстройств мозга.

Рисунок 1: Обзор процедуры генерации обученной сети с помощью DLC. Общая схема процесса, с помощью которого отслеживаются особенности глаз животного, а затем анализируются с помощью машинного обучения. Аббревиатура: DLC = DeepLabCut. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 2: Пример автоматического отслеживания косоглазия у репрезентативной мыши CD1. (A) Пример кадра, показывающего косоглазие DLC (цветные точки) на контуре глаза в день лечения, когда мышь не щурит. (B) Пример кадра, показывающего автоматическое обнаружение косоглазия в день лечения с использованием нашей модели DLC. Евклидово расстояние измерялось с использованием среднего расстояния между B и C, синими точками, в верхней и нижней части глаза. Синие наборы точек в верхней и нижней части глаза используются при отслеживании евклидова расстояния. Другие точки (зеленая, желтая, оранжевая, фиолетовая) являются обрамляющими ориентирами, используемыми как для оценки моделью точек евклидова расстояния, так и для фильтрации неоптимального положения головы после сбора данных. Аббревиатура: DLC = DeepLabCut. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 3: Обоснование количества кадров, используемых для обучения модели. (A) Анализ среднеквадратичной ошибки показывает среднее расстояние между прогнозируемыми и наблюдаемыми значениями для тестовых и обучающих наборов данных. Набор обучающих данных представляет кадры, отобранные при обучении модели, а набор тестовых данных представляет необучающие кадры, используемые для проверки того, насколько хорошо модель может идентифицировать похожие, но разные изображения. Мы использовали пять наборов обучающих и тестовых данных и обнаружили, что значения RMSE выровнялись примерно на 300 кадров для тестовой группы. (B) Вероятность того, что данная точка правильно обозначена (среднее + SEM). Это показало, что 400 помеченных вручную кадров были идеальными, потому что исходные наборы данных в среднем имели вероятность выше 0,95, при этом показатель RMSE был наиболее близок к индексу обучающих данных. Это означало, что модель могла точно приблизиться к точкам, на которых она была обучена, а также сообщать о большинстве кадров с высокой вероятностью. Сокращение: RMSE = среднеквадратичная ошибка. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 4: Матрица нерезкости для измерений косоглазия DLC. Мы сэмплировали 300 секунд из восьми видео (пять CGRP и три PBS) и сравнили эти точки с вручную помеченной двоичной оценкой «да» или «нет» для косоглазия. Мы количественно оценили прогнозируемые значения как те, которые были идентифицированы DLC, а фактические значения — как те, которые были оценены вручную человеком. Затем мы сравнили их с данными, полученными вручную, чтобы увидеть, как часто косоглазие было правильно идентифицировано по сравнению с двоичным значением «да» или «нет» косоглазия, оцененным вручную. Сокращения: DLC = DeepLabCut; CGRP = пептид, связанный с геном кальцитонина; PBS = фосфатно-солевой буфер; TP = истинные положительные результаты; FP = ложные срабатывания; FN = ложноотрицательные результаты; TN = истинные отрицательные значения; PPV = положительная прогностическая ценность; NPV = отрицательная прогностическая ценность; TPR = истинно положительный процент; TNR = истинно отрицательный показатель; MCC = коэффициент корреляции Мэтью. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.

Рисунок 5: Фенотип косоглазия в трех различных моделях для обнаружения косоглазия. Верхние две строки содержат одно и то же репрезентативное животное с каждым условием (PBS или CGRP) в трех разных моделях обнаружения косоглазия. Нижняя строка отражает средние значения по всем животным. (A) Наблюдалось уменьшение средней площади пикселя (средняя общая площадь пикселя/исходный уровень) у мышей, получавших CGRP, по сравнению с мышами, получавшими PBS (t(18) = 2,805, p = 0,012) после обработки всех данных с использованием ранее опубликованной и валидированной модели косоглазия14. (B) Аналогичная реакция наблюдалась в данных, полученных вручную (t(18) = 4,064, p = 0,0007). (C) У мышей, получавших CGRP, наблюдалось снижение среднего расстояния от века до века (евклидово расстояние/евклидово расстояние до обработки, исходный уровень) по сравнению с мышами, получавшими PBS (t(18) = 3,040, p = 0,007 при использовании DLC для обработки всех данных. N = 20 (10 самок, 10 самцов). Полосы погрешностей указывают среднее значение ± SEM. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.