Генетика, 2024, № 12
научный журнал
Покупка
Новинка
Тематика:
Общая генетика. Общая цитология
Издательство:
Наименование: Генетика
Год издания: 2024
Кол-во страниц: 126
Дополнительно
Тематика:
ББК:
УДК:
ОКСО:
ГРНТИ:
Скопировать запись
Фрагмент текстового слоя документа размещен для индексирующих роботов
Российская академия наук ГЕНЕТИКА Том 60 № 12 2024 Декабрь Основан в апреле 1965 r. ISSN: 0016-6758 Ежемесячный журнал Журнал издается под руководством Отделения биологических наук РАН Главный редактор Н.К. Янковский Редакционная коллегия: А.П. Рысков (зам. главного редактора), С.К. Абилев (зам. главного редактора), С.А. Брускин (ответственный секретарь), А.М. Боронин, А.В. Васильев, В.А. Гвоздев, Е.К. Гинтер, Т.А. Ежова, И.А. Захаров-Гезехус, С.Г. Инге-Вечтомов, Н.А. Колчанов, А.М. Кудрявцев, Л.А. Лутова, А.С. Миронов, Н.С. Мюге, Д.В. Политов, В.П. Пузырев, А.Ю. Ржецкий (США), Н.Б. Рубцов, М.В. Холодова, Э.К. Хуснутдинова Редакционный совет: В.Г. Дебабов, А.В. Кильчевский (Беларусь), С.В. Костров, К. Крутовский (Германия), С.А. Лимборская, И.А. Тихонович, Д. Уотсон (США), С.В. Шестаков, В. К. Шумный Зав. редакцией Е.В. Тихомирова Адрес редакции: 119991, ГСП-1, Москва ул. Губкина, д. 3, тел.: 8-499-135-50-45 e-mail: genetika@vigg.ru Сайт журнала: http://www.vigg.ru/genetika/ Москва ФГБУ «Издательство «Наука» © Российская академия наук, 2024 © Редколлегия журнала «Генетика» (составитель), 2024
СОДЕРЖАНИЕ ТОМ 60, номер 12, 2024 Обзорные и теоретические статьи Номенклатура маркеров и контроль ошибок первого и второго рода при проведении генетических экспертиз достоверности происхождения крупного рогатого скота с использованием фрагментного анализа М. В. Модоров, И. В. Ткаченко, А. А. Клещева, М. Ю. Севостьянов 3 Индукция тренированного иммунитета с помощью БЦЖ: недавние данные и мнения И. В. Алексеенко, Л. Г. Кондратьева, И. П. Чернов, Е. Д. Свердлов 16 Генетика микроорганизмов Репаративная ДНК полимераза eta (Polη) играет ключевую роль в мутагенезе при низких дозах УФ-облучения у дрожжей Saccharomyces cerevisiae Е. А. Алексеева, Т. А. Евстюхина, И. И. Скобелева, В. Т. Пешехонов, В. Г. Королев 25 Генетика растений Транскрипционная активность CCA1 у растений Arabidopsis thaliana северной популяции в условиях измененного светового режима М. В. Зарецкая, О. М. Федоренко 36 Гены партеногенеза кукурузы: сравнительный анализ мутаций Е. М. Моисеева, В. В. Фадеев, Ю. В. Фадеева, С. И. Мазилов, А. Ю. Колесова, М. И. Чумаков 43 Генетика животных Генетическая изменчивость и филогенетические взаимоотношения популяций обыкновенной щуки (Esox lucius L.) из некоторых рек России Л. Т. Бачевская, В. В. Переверзева, А. А. Примак, Г. А. Агапова 54 Оценка полиморфизма генов иммунной системы TNF и TLR2 у крупного рогатого скота Л. А. Калашникова, И. Е. Багаль, Н. Е. Муругина, В. Е. Калашников, Р. Ю. Сенина 66 Характеристика по STR-маркерам серого украинского скота, разводимого в Российской Федерации А. С. Мокеев, Н. Н. Фурса, С. В. Бекетов, Г. Р. Свищева, А. А. Онохов, Ю. А. Столповский 75 Генетика человека Популяционно-генетическая структура населения Сибири по данным о частотах полиморфных вариантов генов транспорта и рецепции витамина D Б. А. Тхоренко, А. В. Мейер, Г. В. Вавин, Д. О. Имекина, М. В. Ульянова, Ф. А. Лузина, М. Б. Лавряшина 83
Математические модели и методы О стабилизации темпа дивергенции изонимии В. П. Пасеков 94 Краткие сообщения Паритетное соотношение самцов и самок в выводках красавки Anthropoides virgo в Калмыкии и Бурятии Е. А. Мудрик, Е. И. Ильяшенко, К. Д. Кондракова, А. А. Абушин, Л. Д. Базаров, Ц. З. Доржиев, Д. В. Политов 110 Влияние генотипа по гену алкогольдегидрогеназы ADH1B на уровень биомаркера злоупотребления алкоголем CDT А. А. Ким, А. С. Гуреев, А. В. Рубанович, С. А. Боринская 115 Персоналии Поздравление с 90-летним юбилеем академика РАН Шестакова Сергея Васильевича 121 Поправки 122
ГЕНЕТИКА, 2024, том 60, № 12, с. 3–15 ОБЗОРНЫЕ И ТЕОРЕТИЧЕСКИЕ СТАТЬИ УДК 636.082.21 НОМЕНКЛАТУРА МАРКЕРОВ И КОНТРОЛЬ ОШИБОК ПЕРВОГО И ВТОРОГО РОДА ПРИ ПРОВЕДЕНИИ ГЕНЕТИЧЕСКИХ ЭКСПЕРТИЗ ДОСТОВЕРНОСТИ ПРОИСХОЖДЕНИЯ КРУПНОГО РОГАТОГО СКОТА С ИСПОЛЬЗОВАНИЕМ ФРАГМЕНТНОГО АНАЛИЗА © 2024 г. М. В. Модоров1, *, И. В. Ткаченко1, А. А. Клещева1, М. Ю. Севостьянов1 1Уральский федеральный аграрный научно-исследовательский центр Уральского отделения Российской академии наук, Екатеринбург, 620142 Россия *e-mail: mmodorov@gmail.com Поступила в редакцию 13.05.2024 г. После доработки 18.06.2024 г. Принята к публикации 09.07.2024 г. Ошибочные записи родословных снижают качество племенной работы с крупным рогатым скотом, поэтому проведение генетических экспертиз достоверности происхождения стало неотъемлемой частью племенной работы. В течение многих лет на территории Российской Федерации экспертизы проводили с использованием иммуногенетических маркеров, однако совершенствование технологий и ужесточение требований регулятора привело к тому, что в последние годы на смену иммуногенетике приходят микросателлитные маркеры. В настоящее время отечественный протокол, четко регламентирующий процедуру проведения генетической экспертизы достоверности происхождения крупного рогатого скота с использованием микросателлитных локусов, отсутствует, что затрудняет работу испытательных лабораторий. В частности, требования к числу и номенклатуре генетических маркеров, которые должны быть использованы при проведении экспертиз, имеются только для скота, эмбрионов и спермопродукции, перемещаемых между государствами-членами Евразийского экономического союза. Отсутствуют требования к контролю ошибок первого (ложноположительные результаты) и второго (ложноотрицательные результаты) рода, учет которых необходим при формировании заключений экспертизы. В настоящей работе мы рассмотрим подходы к решению этих вопросов, предложенные Международным обществом генетики животных (ISAG), Международным комитетом регистрации животных (ICAR), Коллегией Евразийской экономической комиссии, а также отечественными нормативными документами, регламентирующими производство судебно-медицинских экспертиз, связанных с установлением родительства. По результатам обзора будут предложены номенклатура микросателлитных маркеров и протокол экспертизы, в который заложен контроль ошибок первого и второго рода. Особое внимание будет уделено описанию источников ошибок второго рода и необходимости их контроля при выдаче заключений. Ключевые слова: Bos taururs, микросателлиты, селекция, ISAG, STR, индексы исключения родителей. DOI: 10.31857/S0016675824120012 EDN: WAWXIF родословных. Вначале подобные тесты проводили с использованием иммуногенетических маркеров (групп крови). В настоящее время в качестве генетических маркеров используют STR-локусы (short tandem repeat, микросателлиты) либо SNP-локусы (single nucleotide polymorphism, однонуклеотидные полиморфизмы). На территории РФ работа лабораторий, проводящих подтверждение достоверности происхождения скота с использованием генетических маркеров, в настоящее время требует получения К середине XX в. искусственное осеменение крупного рогатого скота прочно вошло в практику животноводства. При этом возросло количество ошибочных записей “отец – потомок”, что вносит значительное искажение как в оценку индексов наследования хозяйственно ценных признаков, так и в оценку качества быков по потомству [1, 2]. Для преодоления возникших трудностей в практику работы с племенным скотом вошли генетические тесты, использование которых позволило значительно снизить число ошибочных записей в 3
МОДОРОВ и др. Свидетельства о регистрации в государственном племенном регистре. Ситуация изменилась с вступлением в силу Приказа Минсельхоза России № 336 от 2 июня 2022 г. “Об утверждении требований к видам племенных хозяйств”. Согласно пункту 8 данного Приказа к лабораториям молекулярно-генетической экспертизы в качестве дополнительного требования устанавливается “наличие документа об аккредитации в национальной системе аккредитации”. К настоящему моменту отечественный ГОСТ либо методические указания, содержащие обязательные условия к процедуре подтверждения достоверности происхождения скота с использованием фрагментного анализа, отсутствуют. Это вносит значительную неопределенность при подготовке лабораторий к аккредитации, в связи с чем интерес представляет анализ требований, предъявляемых к данному типу анализов за рубежом, а также положений отечественных нормативных документов, регламентирующих производство судебно-медицинских экспертиз, связанных с установлением родительства. В настоящей работе мы рассмотрим вопросы, касающиеся научной составляющей проведения исследований достоверности происхождения крупного рогатого скота (обсуждение требований к лабораторным помещениям и организации менеджмента качества испытательных лабораторий в наши задачи не входит), а именно: 1) генетическая идентификация (выбор генетических маркеров, их числа и номенклатуры); 2) достоверность происхождения (сопоставление микросателлитных профилей с учетом контроля ошибок первого и второго рода, т. е. ложноположительных и ложноотрицательных результатов теста). Цель настоящей работы – разработка протокола сопоставления генетических профилей предполагаемых родителей и потомка, в котором учтены интересы основных участников и выгодоприобретателей проводимых экспертиз, а именно племенной службы, сельхозпроизводителей, поставщиков спермопродукции и лабораторий, проводящих экспертизы. Обзор требований, предъявляемых к проведению фрагментного анализа, связанного с подтверждением родительства ISAG. Наиболее влиятельной международной организацией, определяющей тип, число и номенклатуру генетических маркеров, используемых в тестах для подтверждения достоверности происхождения сельскохозяйственных животных, является Международное общество генетики животных (International Society for Animal Genetics, ISAG). В частности, ISAG разрабатывает руководства по управлению генетическими ресурсами сельскохозяйственных животных в интересах FAO (Food and Agriculture Organization, Продовольственная и сельскохозяйственная организация ООН) [3, 4]. К 2006 г. для проведения тестов на достоверность происхождения крупного рогатого скота ISAG было рекомендовано девять микросателлитных локусов: BM1824, INRA23, BM2113, SPS115, ETH10, TGLA122, ETH225, TGLA126 и TGLA227 [5]. В 2008 г. к ним было добавлено еще три локуса, а именно: BM1818, ETH3 и TGLA53 [6]. По состоянию на 2021-й год минимальное число STR-локусов, рекомендованных ISAG для подтверждения достоверности происхождения крупного рогатого скота, не изменилось, оно остается равным 12 [7]. Раз в два года ISAG разрабатывает правила и проводит сравнительные испытания (comparison tests), в ходе которых лаборатории получают пробирки с выделенной ДНК нескольких пород скота. Задачей лабораторий, добровольно участвующих в данном испытании, является корректная расшифровка неизвестных генотипов по 12 локусам, т. е. генетическая идентификация скота. ICAR. Помимо ISAG, сертификацию лабораторий, проводящих подтверждение достоверности происхождения крупного рогатого скота, проводит ICAR (International Committee for Animal Recording, Международный комитет регистрации животных). ICAR рекомендует использовать тот же набор STR-маркеров, что и ISAG. Для аккредитации в ICAR лаборатории необходимо соответствовать более жестким требованиям, чем в ISAG, в частности, иметь сертификат аккредитации по ISO 17025. Кроме того, лабораториям необходимо проводить расчеты индексов исключения родителей (еxclusion probability) как для каждого тестируемого локуса, так и для всей панели маркеров. Популяция, для которой проводят расчет индексов, должна быть охарактеризована по породной принадлежности и числу проанализированных голов, которых должно быть не менее 150 [8]. Решение Коллегии Евразийской экономической комиссии. Русскоязычным документом, определяющим номенклатуру локусов, используемых для подтверждения достоверности происхождения крупного рогатого скота, и требования к проводящим тесты лабораториям, является Решение Коллегии Евразийской экономической комиссии от 2 июня 2020 г. № 74 “Об утверждении Положения о проведении молекулярной генетической экспертизы племенной продукции государств – членов Евразийского экономического союза” [9]. Согласно данному документу, проведение молекулярной генетической экспертизы осуществляется лабораториями, аккредитованными в государственных (национальных) системах аккредитации ГЕНЕТИКА том 60 № 12 2024
НОМЕНКЛАТУРА МАРКЕРОВ И КОНТРОЛЬ ОШИБОК ПЕРВОГО И ВТОРОГО РОДА 5 государств – членов Евразийского экономического союза (далее государства – члены) либо ICAR (пункт 3). Обязательной молекулярной генетической экспертизе подлежат племенные производители сельскохозяйственных животных (крупный рогатый скот, лошади, овцы, козы, свиньи, олени, верблюды), перемещаемые между государствами-членами, а также племенные производители и доноры эмбрионов сельскохозяйственных животных, спермопродукция и эмбрионы которых перемещаются между государствами-членами (п. 6). Молекулярная генетическая экспертиза племенной продукции осуществляется методом ДНК-типирования с использованием методик, разработанных с учетом рекомендаций ISAG, в соответствии с областью аккредитации лаборатории (п. 7). Проведение молекулярной генетической экспертизы с целью подтверждения достоверности происхождения племенной продукции осуществляется методом генотипирования STR-маркеров или SNP-маркеров (п. 9). Перечень STR-маркеров должен составлять не менее 12 STR-маркеров, включенных в базовую STR-панель ISAG. В Приказе Минсельхоза России № 336 от 2 июня 2022 г. “Об утверждении требований к видам племенных хозяйств” номенклатура генетических маркеров, требуемых для проведения генетических тестов на достоверность происхождения скота, отсутствует [10]. Однако этот документ определяет группы крупного рогатого скота, для которых генетическая экспертиза достоверности происхождения должна быть проведена. Это скот, содержащийся в племенных заводах, племенных репродукторах и генофондных хозяйствах; молочные породы скота (100% ремонтного молодняка для реализации, 100% матерей быков-производителей), мясные породы скота (100% быков-производителей и коров быкопроизводящей группы). Нормативная база, регламентирующая производство судебно-медицинских экспертиз, связанных с установлением родительства. В настоящее время (весна 2024 г.) на территории РФ действует Приказ Министерства здравоохранения и социального развития Российской Федерации от 12 мая 2010 г. № 346н "Об утверждении Порядка организации и производства судебно-медицинских экспертиз в государственных судебно-экспертных учреждениях Российской Федерации" [11]. Особенности порядка производства генетической экспертизы описаны в 84-м разделе данного документа. Согласно ему генетическая экспертиза по поводу спорного происхождения детей отвечает на вопросы: 1) исключается или не исключается отцовство, материнство данного индивидуума в отношении данного ребенка (плода); 2) если отцовство, материнство не исключается, то какова вероятность того, что полученный результат не является следствием случайного совпадения индивидуализирующих признаков неродственных лиц. Данный документ не регламентирует число и номенклатуру используемых локусов, однако предъявляет требования к вероятностной оценке неисключенного родства (пп. 84.12.5 и 84.12.6). Принимается, что совпадение условно отцовских (нематеринских) аллелей в генотипе ребенка с аллелями, присутствующими в геноме предполагаемого отца, не означает доказанного отцовства. Вероятностная оценка неисключенного родства строго обязательна. Искомое условное значение вероятности отцовства соответствует байесовской постериорной вероятности при 50%-ной априорной вероятности отцовства. Уровень доказательности экспертного исследования в случае неисключения отцовства, материнства должен составлять для полного трио (мать – ребенок – предполагаемый отец) при условии, что истинность другого родителя считается бесспорной, не ниже 99.90% (рассчитываемый как байесовская вероятность отцовства/материнства); не ниже 1000 (рассчитываемый как индекс отцовства PI). Для дуэта (ребенок – предполагаемый отец) в отсутствие другого родителя – не ниже 99.75% (рассчитываемый как байесовская вероятность отцовства/материнства); не ниже 400 (рассчитываемый как индекс отцовства PI). Согласно п. 84.12.7 при количественном анализе результата для расчета вероятности отцовства (материнства) следует учитывать этническую принадлежность обследуемых лиц. Кроме этого, в Приказе предъявляются требования к отрицательному заключению (п. 84.12.4). Для обоснованного вывода о безусловном исключении отцовства, материнства аллели ребенка, не свойственные ни одному из указанных родителей, должны быть зарегистрированы как минимум в двух несцепленных локусах. С 1 сентября 2024 г. на смену Приказу Министерства здравоохранения и социального развития Российской Федерации от 12 мая 2010 г. № 346н приходит Приказ № 491н от 25 сентября 2023 г. “Об утверждении Порядка проведения судебно-медицинской экспертизы” [12]. Приложение № 11 нового приказа посвящено “Правилам организации деятельности отделения генетической экспертизы и проведения генетической экспертизы”. Согласно этому документу, для вывода об исключении отцовства или материнства аллели ребенка, не свойственные ни одному из предполагаемых родителей, регистрируются сразу в нескольких несцепленных локусах. Для вывода о том, какое лицо является родителем данного ребенка, необходима вероятностная оценка наблюдаемого ГЕНЕТИКА том 60 № 12 2024
МОДОРОВ и др. В исследованиях крупного рогатого скота для контроля ошибки первого рода ICAR требует рассчитывать индексы исключения родителей, однако в просмотренных документах нам не удалось найти порогового значения индексов, превышение которого необходимо для оформления заключений. Контроль ошибки второго рода у крупного рогатого скота при использовании панели STR-маркеров, по-видимому, не проводится. Однако при использовании для анализа SNP-маркеров допускается, что подтверждение родительства происходит при отсутствии общих аллелей в двух из 186 протестированных локусов, при несовпадении аллелей в трех – пяти локусах предполагаемый родитель получает статус “сомнительный” (doubtful) и только при несовпадении более пяти локусов родитель исключается [8]. Индексы исключения родителей (контроль ошибки первого рода или выдача ложноположительных заключений) совпадения/несовпадения индивидуальных профилей полиморфизма ДНК сравниваемых биологических объектов в рамках интерпретационной модели противоположных версий: по умолчанию – версии несовпадения аллелей по причине нарушения базовых закономерностей их совпадения и версии аллельного несовпадения у не состоящих в родстве лиц. Расчетная вероятностная оценка рассматриваемых версий спорного родительства (исключение/неисключение отцовства или материнства) соответствует условной байесовской апостериорной вероятности при 50%-ной априорной вероятности родительства. При количественном анализе результата для расчета вероятности отцовства (материнства) учитывается этническая принадлежность обследуемых лиц. Методическое содержание и объем исполнения завершенного исследования спорного происхождения детей предусматривают в типовом случае анализа генотипирование у каждого из обследуемых лиц не менее 16 несцепленных полиморфных тандемных локусов-маркеров аутосомной ДНК. Проанализированные материалы позволяют говорить о том, что к настоящему времени существует консенсус в номенклатуре микросателлитных маркеров, используемых для проведения экспертиз на достоверность происхождения крупного рогатого скота, а именно использование не менее 12 локусов базовой панели ISAG. В то же время в отличие от требований, предъявляемых к проведению судебно-медицинских экспертиз, контроль ошибок первого и второго рода у крупного рогатого скота не регламентирован. Источники ошибок при сравнении генотипов предполагаемых родителей и потомка При проведении экспертизы достоверности происхождения крупного рогатого скота необходимо контролировать ошибки первого и второго рода. Ошибка первого рода (ложноположительное заключение) – ситуация, в которой подтверждается родительство особи, которая родителем не является. Ошибка второго рода (ложноотрицательное заключение) – выдача отрицательного заключения на настоящего родителя. При проведении судебно-медицинских экспертиз для контроля ошибки первого рода вводятся пороговые значения индексов байесовской вероятности отцовства/материнства либо индекса отцовства (PI). Кроме этого, в Приказе № 491н от 25 сентября 2023 г. вводится требование к минимальному числу локусов, используемых для анализа. Для контроля ошибки второго рода вывод об исключении отцовства или материнства может быть сделан в случае наличия у ребенка аллелей, не свойственных ни одному из предполагаемых родителей, сразу в нескольких несцепленных локусах (минимум двух). Наличие общих аллелей по каждому исследованному локусу у предполагаемого родителя и потомка может быть объяснено двумя способами: 1) потомок происходит от предполагаемого родителя; 2) аллели у анализируемых индивидов совпали случайно. Высокая вероятность случайного совпадения аллелей у предполагаемого родителя и потомка приведет к высокой частоте ложноположительных экспертных заключений, подтверждающих происхождение крупного рогатого скота (т. е. к статистической ошибке первого рода). Для контроля этой ошибки при проведении судебно-медицинских экспертиз требуется рассчитать индекс байесовской вероятности отцовства/материнства либо индекс отцовства (PI). Эти индексы рассчитываются с учетом частоты встречаемости каждого обнаруженного аллеля в популяции (этносе) [13–15]. Таким образом, при использовании единой панели генетических маркеров в каждой проведенной экспертизе эта вероятность требует отдельного расчета. ICAR следует иному подходу и для контроля ошибки первого рода предлагает рассчитывать индексы исключения родителей [8]. Формулы, позволяющие оценить эти индексы, приведены на сайте ISAG [16], более подробную информацию об этих формулах можно найти в работе Джеймесона и Тэйлора [17]. Эти индексы характеризуют локус (а не аллель), поэтому могут быть рассчитаны для определенной популяции единожды. По сути индексы исключения родителей показывают вероятность того, что у случайно отобранных в ГЕНЕТИКА том 60 № 12 2024
НОМЕНКЛАТУРА МАРКЕРОВ И КОНТРОЛЬ ОШИБОК ПЕРВОГО И ВТОРОГО РОДА 7 однако при проведении судебно-медицинских экспертиз этот индекс может быть использован: 2 n n n n n n n 2 2 3 4 5 3 3 2 = − + + − − (3) P1 1 2 2 3 2 i p p i i p p p p p + 1 i 1 i i 1 i i 1 i i 1 i i 1 i i 1 = = = = = ∑ ∑ ∑ ∑ ∑ = = ∑∑ , где pi – частоты аллелей в тестируемом локусе. С ростом значения индекса растет его информативность для решения задач по проведению генетических экспертиз на достоверность происхождения крупного рогатого скота. Если в исследовании использовано несколько (n) локусов, то для оценки индексов исключения родителей с использованием всей панели маркеров может быть применена формула: популяции особей в тестируемом локусе не будет сходного аллеля [18]. Отметим, что значительную путаницу при использовании индексов исключения родителей вносит то, что в публикации Джеймесона и Тэйлора все они называются просто “P”, без какого-либо дополнительного буквенного или цифрового обозначения. В настоящей работе вслед за разработчиками программы GenAlEx 6.5 [19] мы будем называть эти индексы P1, P2, P3, что соответствует последовательности написания формул в работе [17]. 1. Формула, используемая для расчета индекса исключения родителей, определяется тем, какие особи были выбраны для сравнения генотипов. В наиболее простом случае проводят сравнение генотипов потомка и одного предполагаемого родителя (данные о втором родителе в этом тесте не важны, формула 1, индекс P2). (4) n n n n 2 2 2 3 , (1) P2 1 4 2 4 3 4 = − + i p p p p i i i + − 1 i 1 i 1 i 1 = = = = ∑ ∑ ∑ ∑ i где pi – частоты аллелей в тестируемом локусе. 2. Во второй ситуации тестируется происхождение потомка от двух предполагаемых родителей (индекс P3, формула 2), т. е. вероятность отсутствия у потомка хотя бы одного аллеля, один из которых сходен с аллелем предполагаемого отца, а другой с аллелем предполагаемой матери: n n n n 4 5 6 2 2 P3 1 4 4 3 8 p p p p = + − − − i + ∑ ∑ ∑ ∑ i 1 i i 1 i i 1 i i 1 = = = = n n n (2) 3 3 2 8 + i 2 ∑ ∑ ∑ + 2 ( ) ( ) , 1 i i 1 i 1 = = p p p i i = где pi – частоты аллелей в тестируемом локусе. В практической работе лабораторий, специализирующихся на проведении генетических экспертиз на достоверность происхождения крупного рогатого скота, абсолютное большинство тестов описывается двумя вышеприведенными формулами. 3. Формула 3 (индекс P1) применима в ситуации установления отцовства, когда проверке подлежит не два родителя, а только один (например, в происхождении ребенка от матери сомнений нет, а вот отец сомнения вызывает). При этом известны генотипы всей триады: потомок, известный родитель и предполагаемый родитель. В практической работе по проведению генетических экспертиз на достоверность происхождения крупного рогатого скота с подобной ситуацией мы не сталкивались, , где Pk – значение индекса в каждом тестируемом локусе. В табл. 1 приведены рассчитанные значения индексов исключения родителей для популяции племенного голштинизированного черно-пестрого скота Свердловской области. Для расчетов использована база данных, включающая генотипы 24544 голов скота, расшифровка которых была проведена в лаборатории молекулярно-генетической экспертизы Уральского НИИСХ – филиала ФГБНУ УрФАНИЦ УрО РАН в период с 2017 по 2021 г. Методика генотипирования и описание популяции приведены нами ранее [20]. Индекс P3 для панели, включающей 12 рекомендованных ISAG локусов, составляет 0.999999705. Это значение позволяет сформировать ожидания относительно результатов 10 миллионов сравнений генотипов, случайно взятых из популяции отцов, матерей и потомка. Результаты 9 999 997 подобных сравнений покажут, что хотя бы в одном локусе триады отец – мать – потомок для потомка нельзя будет подобрать аллели, наследуемые от предполагаемого отца и предполагаемой матери. В трех оставшихся тестах в каждом локусе потомка будет обнаружен аллель, сходный с аллелем предполагаемого отца и предполагаемой матери (ложноположительное сходство). Значение индекса P2 для панели, включающей 12 рекомендованных ISAG локусов, составляет 0.9946. Это значение показывает, что при сравнении генотипов, случайно взятых из популяции матери и потомка, в 9946 тестах из 10000 будет получен результат, в котором хотя бы в одном локусе не будет обнаружен общий аллель, а в 54 оставшихся тестах общий аллель будет обнаружен во всех ГЕНЕТИКА том 60 № 12 2024
МОДОРОВ и др. Контроль ошибки второго рода (ложноотрицательные заключения) Таблица 1. Индексы исключения родителей в популяции голштинизированного черно-пестрого скота Свердловской области по данным изменчивости 24544 голов. Значения индексов приведены для каждого локуса, для панели из 12 локусов, рекомендованных ISAG, а также панели из 15 локусов, входящих в набор COrDIS Cattle Локус Индексы исключения родителей P1 P2 P3 BM1824* 0.47 0.30 0.65 BM2113* 0.55 0.37 0.74 ETH3* 0.44 0.27 0.62 ETH10* 0.58 0.40 0.77 ETH225* 0.44 0.27 0.63 INRA23* 0.53 0.35 0.71 SPS115* 0.34 0.18 0.53 TGLA53* 0.71 0.55 0.88 TGLA122* 0.68 0.52 0.86 TGLA126* 0.33 0.19 0.50 TGLA227* 0.66 0.48 0.84 BM1818* 0.35 0.21 0.52 Csrm60 0.42 0.25 0.61 Cssm66 0.52 0.34 0.72 Ilsts006 0.40 0.24 0.57 12 локусов 0.999868 0.9946 0.999999705 15 локусов 0.999978 0.9979 0.999999986 Примечание. * – локус включен в базовую панель ISAG. локусах. Это же касается пары предполагаемый отец – потомок. Высокие значения индексов P3, полученные для панели, включающей как 12, так и 15 локусов, на наш взгляд, позволяют проводить не только подтверждение достоверности происхождения, но и восстановление пары родителей (подбор родителей на основании данных о генотипах скота в стаде). Уровни значений индексов P2 позволяют проводить подтверждение достоверности происхождения в популяции, но недостаточны для восстановления родителя. Отметим, что индексы исключения родителей у локальных пород скота во многих случаях ниже, чем у голштинской породы либо голштинизированного скота [21]. Причинами отсутствия общих аллелей у родителя и потомка могут быть de novo мутации, главным образом изменение в числе мотивов микросателлита, мутации в последовательности ДНК, комплементарной последовательности праймера, ошибки генотипирования, путаница при отборе проб и составлении описей. Мутационный процесс. Частота мутаций в динуклеотидных микросателлитных локусах у человека составляет 7.88 × 10−5 на поколение или одно событие на 12690 проведенных локус-тестов. Наиболее часто мутационное событие заключается в изменении числа мотивов микросателлита на один, при этом четкой линейной связи числа повторов с числом мутационных событий для динуклеотидных повторов не отмечено [22]. Если принять сходный уровень скорости мутационного процесса у человека и крупного рогатого скота, то при использовании панели, включающей 12 микросателлитных локусов, мы вправе ожидать около одного мутационного события примерно на 1000 (1 × 10-3) проведенных тестов (точнее на 1058 тестов, 12690 событий / 12 локусов в тесте = 1058). Таким образом, при выставлении требований к 100%-ному совпадению аллелей во всех 12 локусах крупного рогатого скота ложноотрицательные результаты будут получены в каждом тысячном сравнении родителя и потомка либо в каждом 500-м тесте при сравнении двух родителей и потомка. Отметим, что в одном из миллиона проведенных тестов (1 × 10-3 × 1 × 10-3 = 1 × 10-6) вследствие мутационного процесса будет наблюдаться несовпадение аллелей в двух локусах. Мутационный процесс может привести к изменениям в последовательности ДНК, комплементарной последовательности праймера. Следствием мутации может стать либо снижение высоты пика, либо его выпадение (так называемые нуль-аллели). В этом случае гетерозиготный генотип будет считан как гомозигота. Маркером нуль-аллелей может быть отклонение числа наблюдаемых генотипов от их числа, рассчитанного с использованием уравнения Харди–Вайнберга, в сторону увеличения доли гомозигот либо, если нуль-аллели встретятся в гомозиготе, отсутствие пиков в данном локусе. В случае наследования нуль-аллеля генотип потомка будет выглядеть как гомозигота по аллелю, унаследованному от другого родителя. На рис. 1 рассмотрены гипотетические примеры ошибок второго рода, связанных с наличием нуль-аллелей. Отметим, что часть этих ошибок может быть связана с использованием наборов реагентов различных производителей. При этом как ISAG, так и ICAR вводят номенклатуру праймеров [4, 8], однако последовательности праймеров в производимых в РФ ГЕНЕТИКА том 60 № 12 2024
НОМЕНКЛАТУРА МАРКЕРОВ И КОНТРОЛЬ ОШИБОК ПЕРВОГО И ВТОРОГО РОДА 9 Рис. 1. Гипотетический пример ошибок второго рода (ложноотрицательное заключение о происхождении потомка от заявленных родителей), связанных с присутствием мутаций в последовательности ДНК, комплементарной последовательности праймера. В первом тесте при определении генотипа быка (♂), коровы (♀) и потомка (F1) был использован набор реагентов “a”, в состав которого входит праймер, в месте комплементарного связывания которого с матрицей ДНК у быка произошла мутация (показана стрелкой). В результате мутации фрагмент размером 66 пн не амплифицируется (нуль-аллель). Генотип быка, обозначенный числом мотивов во фрагменте, будет прочитан как “10 / 10” (гомозигота), хотя фактически его генотип “8 / 10” (гетерозигота). Генотип коровы с использованием этого набора реагентов читается без ошибок “4 / 12”. Потомок наследует аллель “8” быка и аллель “4” коровы. При расшифровке генотипа потомка этим же набором реагентов считывается генотип “4 / 4”, так как аллель “8” не амплифицируется из-за мутации в месте комплементарного связывания праймера с матрицей ДНК (вариант F1a). Получившиеся данные позволяют сделать заключение об отсутствии в данном локусе общих аллелей у быка и потомка. В другой лаборатории использовали генетические паспорта родителей, а генотип потомка анализировали с использованием другого набора реагентов “b”. В наборе реагентов “b” для анализа данного локуса разработчики включили другие праймеры, комплементарные иным последовательностям матричной ДНК, в результате чего у потомка читаются оба аллеля как “4”, так и “8” (F1b). Однако и в этой ситуации будет сделан вывод об отсутствии в данном локусе общих аллелей у быка и потомка. Ошибки расшифровки. При анализе больших массивов данных ошибки неизбежны. В случае анализа изменчивости микросателлитов они могут быть связаны с малым числом копий выделенной ДНК, присутствием загрязнения в ходе ПЦР, низким качеством электрофореза, плохой калибровкой генетического анализатора и т. п. Оценить коммерческих наборах реагентов, предназначенных для мультиплексного анализа микросателлитных маркеров крупного рогатого скота, а именно: COrDIS Cattle (ООО “Гордиз” [23]) и Gene Profile Cattle (НПК Синтол [24]), являются коммерческой тайной разработчиков. ГЕНЕТИКА том 60 № 12 2024