Модельный расчёт
Совместная модель ПГТ-мт и ПГТ-А
Сценарии совместного отбора и чувствительность к допущениям
Версия модели: 2026-10-12
Что означает результат
P(≥1) — модельная вероятность того, что среди N бластоцист окажется хотя бы одна, удовлетворяющая обоим критериям модели: эуплоидная по ПГТ-А и с гетероплазмией не выше выбранной границы. Это не вероятность беременности, рождения ребёнка или рождения здорового ребёнка и не прогноз для конкретной пациентки.
В исходном наборе 187 записей из 12 публикаций: 79 явно отмечены как бластоцисты, у 57 стадия не указана. У 43 записей «не обнаружено» отсутствует предел обнаружения. Это смешанный набор опубликованных наблюдений, а не полные независимые циклы.
Дополнительный знаменатель для m.14487T>C без 3PN-эмбриона: ≤18% — 5/25. Основной набор сохраняет все опубликованные эмбриональные значения; 3PN не является кандидатом на перенос.
H — доля значений не выше выбранной границы; E — сценарная доля эуплоидных бластоцист; ◆ модельный расчёт.
При 1–2 носительницах доля описывает этих женщин, а не вариант. Неизвестные связи записей оставлены «не установлено». Числа по публикациям приведены отдельно; одна женщина может быть описана в нескольких статьях, поэтому складывать их нельзя.
Носительницы по публикациям (не суммировать)
- Ji et al., 2021 doi:10.1007/s10815-021-02354-3: 1; Clinical case, pp. 3252–3253; Table 3
- Tofilo et al., 2026 doi:10.1007/s10815-026-04012-y: 1; Clinical case; Table 1; Results
- Hu et al., 2026 doi:10.1007/s10815-026-04041-7: 1; Case presentation; Fig. 2; Supplementary Table S1
- Spath et al., 2021 doi:10.1016/j.fertnstert.2021.01.026: 4; Supplemental Table 3, patients A–D
- Treff et al., 2012 doi:10.1016/j.fertnstert.2012.07.1119: 1; Materials and Methods, clinical case; Table 1
- Heindryckx et al., 2014 doi:10.1016/j.mito.2014.08.005: 1; Clinical case; Table 1
- Monnot et al., 2011 doi:10.1002/humu.21417: 2; Table 3, patients 1–2 (embryo series)
- Hyslop et al., 2025 doi:10.1056/NEJMoa2415539: 39; Results, PGT programme: 39 patients; Table 1 is a live-birth subset
- Sallevelt et al., 2017 doi:10.1093/humrep/dew356: 1; Clinical case; Fig. 2; Supplementary Table SI
ρ относится к бинарному событию прохождения обоих тестов. Это заданная корреляция в бета-биномиальной модели, а не измеренная корреляция процентов гетероплазмии. Все варианты используют одно выбранное сценарное значение.
Учебный пример из текущих данных: при H = 41/88, E = 0,5, N = 6 вероятность составляет 0,8 при ρ = 0 и 0,64 при ρ = 0,2. Это сценарный расчёт, а не индивидуальный прогноз.
Сценарий E = 0,7
Сценарий E = 0,5
Сценарий E = 0,3
- H — доля опубликованных значений гетероплазмии не выше выбранной границы; она описывает включённые серии, а не конкретную пациентку.
- E — сценарная доля эуплоидных бластоцист. Значения не привязаны к возрасту. Ориентир: Volodyaev et al., 2025 doi:10.1016/j.xfnr.2024.100082 — Volodyaev et al. 2025, F S Rev 6:100082, doi:10.1016/j.xfnr.2024.100082: с. 2, раздел «Incidence, parental origin, and age relation» (анеуплоидии целых хромосом 20–30% в 25–35 лет, почти линейный рост ≈6% в год после 35 лет, 80–90% в 42–44 года); с. 3, рис. 2 (доля эуплоидных бластоцист по возрасту, данные Franasiak et al.); значения E — сценарные, с точками рисунка не сопоставлены
P(≥1) для всех вариантов набора
H = k/n рассчитывается по включённым записям варианта. В общем наборе смешаны стадии дробления, 3PN и значения только перенесённых эмбрионов (Hyslop 2025, Thorburn 2009); отбор по низкой гетероплазмии смещает H вверх. Фильтр оставляет только явно отмеченные биоптаты ТЭ бластоцист. Интервалы засчитываются по верхней границе; «не обнаружено» условно ниже границы. Интервалы Уилсона не учитывают кластеризацию.
| Вариант мтДНК | H ≤18% | H ≤30% | Состав входной выборки | Сценарий E = 0,7 | Сценарий E = 0,5 | Сценарий E = 0,3 |
|---|---|---|---|---|---|---|
| m.3243A>G | 0,47 (41/88) 95% ДИ Уилсона (без учёта кластеризации): 0,37–0,57 | 0,60 (53/88) 95% ДИ Уилсона (без учёта кластеризации): 0,50–0,70 | Публикаций: 5 Бластоцист: 28; стадия не указана: 49 Не обнаружено, предел обнаружения не указан: 14 | ≤18%: 0,69 ≤30%: 0,81 | ≤18%: 0,55 ≤30%: 0,66 | ≤18%: 0,36 ≤30%: 0,45 |
| m.3697G>A | 0,05 (1/20) 95% ДИ Уилсона (без учёта кластеризации): 0,01–0,24 | 0,15 (3/20) 95% ДИ Уилсона (без учёта кластеризации): 0,05–0,36 | Публикаций: 1 Бластоцист: 15; стадия не указана: 0 | ≤18%: 0,10 ≤30%: 0,28 | ≤18%: 0,07 ≤30%: 0,21 | ≤18%: 0,04 ≤30%: 0,13 |
| m.14487T>C | 0,21 (6/28) 95% ДИ Уилсона (без учёта кластеризации): 0,10–0,40 | 0,29 (8/28) 95% ДИ Уилсона (без учёта кластеризации): 0,15–0,47 | Публикаций: 2 Бластоцист: 6; стадия не указана: 0 | ≤18%: 0,39 ≤30%: 0,49 | ≤18%: 0,29 ≤30%: 0,37 | ≤18%: 0,18 ≤30%: 0,24 |
| m.10197G>A малое n — оценка доли неустойчива | 0,56 (5/9) 95% ДИ Уилсона (без учёта кластеризации): 0,27–0,81 | 0,56 (5/9) 95% ДИ Уилсона (без учёта кластеризации): 0,27–0,81 Пересекающих границу интервалов: 1 | Публикаций: 1 Бластоцист: 9; стадия не указана: 0 | ≤18%: 0,77 ≤30%: 0,77 | ≤18%: 0,62 ≤30%: 0,62 | ≤18%: 0,42 ≤30%: 0,42 |
| m.10191T>C малое n — оценка доли неустойчива | 1,00 (8/8) 95% ДИ Уилсона (без учёта кластеризации): 0,68–1,00 | 1,00 (8/8) 95% ДИ Уилсона (без учёта кластеризации): 0,68–1,00 | Публикаций: 1 Бластоцист: 0; стадия не указана: 3 Не обнаружено, предел обнаружения не указан: 8 | ≤18%: 0,97 ≤30%: 0,97 | ≤18%: 0,88 ≤30%: 0,88 | ≤18%: 0,66 ≤30%: 0,66 |
| m.8993T>G | 0,92 (24/26) 95% ДИ Уилсона (без учёта кластеризации): 0,76–0,98 | 0,92 (24/26) 95% ДИ Уилсона (без учёта кластеризации): 0,76–0,98 | Публикаций: 4 Бластоцист: 15; стадия не указана: 5 Не обнаружено, предел обнаружения не указан: 21 | ≤18%: 0,96 ≤30%: 0,96 | ≤18%: 0,84 ≤30%: 0,84 | ≤18%: 0,62 ≤30%: 0,62 |
| m.3250T>C Вариант неопределённого значения (VUS) малое n — оценка доли неустойчива | 1,00 (8/8) 95% ДИ Уилсона (без учёта кластеризации): 0,68–1,00 | 1,00 (8/8) 95% ДИ Уилсона (без учёта кластеризации): 0,68–1,00 | Публикаций: 1 Бластоцист: 6; стадия не указана: 0 | ≤18%: 0,97 ≤30%: 0,97 | ≤18%: 0,88 ≤30%: 0,88 | ≤18%: 0,66 ≤30%: 0,66 |
Формулы
q = H × E · P(≥1) = 1 − (1 − q)^N
При 0 < ρ < 1 общая вероятность Q имеет бета-распределение со средним q. Эмбрионы независимы условно на Q; после интегрирования возникает положительная зависимость. При ρ = 0 используется биномиальный предел.
Q ~ Beta(α, β); α = q(1/ρ − 1), β = (1 − q)(1/ρ − 1)
P(≥1) = 1 − B(α, β + N) / B(α, β)
Допущения (не доказаны)
- q = H × E: гетероплазмия и хромосомный статус считаются независимыми;
- при ρ = 0 бластоцисты считаются независимыми испытаниями с одинаковой вероятностью; N — заданное число бластоцист, а не число ооцитов или средний выход цикла;
- при ρ > 0 используется общая случайная вероятность совместного прохождения. Значение ρ выбрано для сценария и не оценено по данным Атласа;
- H описывает объединённые опубликованные данные разных стадий и методов; E — сценарная доля эуплоидных бластоцист.
Для оценки семейных и вариантных эффектов нужны полные связанные циклы нескольких независимых семей. Сценарные доли эуплоидии не калиброваны по возрасту.
Если независимость двух тестов не установлена
Маргинальные H и E определяют только границы возможного совместного q для одного эмбриона. Таблица показывает эти границы, а не доверительные интервалы. Они условны на входные доли и не описывают зависимость между эмбрионами.
max(0, H + E − 1) ≤ q ≤ min(H, E)
| Вариант мтДНК | Граница | E = 0,7 | E = 0,5 | E = 0,3 |
|---|---|---|---|---|
| m.8993T>G | ≤18% | 0,62–0,70 | 0,42–0,50 | 0,22–0,30 |
| m.8993T>G | ≤30% | 0,62–0,70 | 0,42–0,50 | 0,22–0,30 |
| m.3243A>G | ≤18% | 0,17–0,47 | 0,00–0,47 | 0,00–0,30 |
| m.3243A>G | ≤30% | 0,30–0,60 | 0,10–0,50 | 0,00–0,30 |
| m.14487T>C | ≤18% | 0,00–0,21 | 0,00–0,21 | 0,00–0,21 |
| m.14487T>C | ≤30% | 0,00–0,29 | 0,00–0,29 | 0,00–0,29 |
| m.3697G>A | ≤18% | 0,00–0,05 | 0,00–0,05 | 0,00–0,05 |
| m.3697G>A | ≤30% | 0,00–0,15 | 0,00–0,15 | 0,00–0,15 |
Совместные результаты первичных серий
Фиксированный критерий здесь — H≤18% и строгая эуплоидия. Знаменатель — ТЭ-бластоцисты каждой опубликованной серии. Эти малые серии описывают конкретных носительниц; их результаты не подставляются в сценарные кривые как параметры новых семей.
| Источник | Вариант мтДНК | H≤18% | Подтверждённо эуплоидные | Прошли оба критерия |
|---|---|---|---|---|
| Ji et al., 2021 doi:10.1007/s10815-021-02354-3 Table 3, p. 3256 Одна пациентка, два цикла; ПГТ-А информативен у 14/15. У оставшейся бластоцисты H=99,4%, поэтому совместный критерий не пройден независимо от неизвестной плоидности. | m.3697G>A | 1/15 | 3/15 | 0/15 |
| Spath et al., 2021 doi:10.1016/j.fertnstert.2021.01.026 Supplemental Table 3, pp. 1–2 Пациентка B. Во всех биопсиях вариант не обнаружен; предел обнаружения не указан. Доли условны на уровень ниже 18%, если вариант не обнаружен. При постоянном статусе H проверить независимость тестов нельзя. | m.8993T>G | 14/14 | 10/14 | 10/14 |
| Spath et al., 2021 doi:10.1016/j.fertnstert.2021.01.026 Supplemental Table 3, pp. 2–3 Пациентка D. Доли с ND условны на предел обнаружения ниже 18%. Совместные 3/9 отличаются от произведения маргинальных долей 8/27; малый объём не позволяет установить ассоциацию. | m.3243A>G | 8/9 | 3/9 | 3/9 |
| Hu et al., 2026 doi:10.1007/s10815-026-04041-7 Fig. 2; embryo results E3, E6, E8 Одна носительница варианта неопределённого значения (VUS); результат не характеризует патогенные варианты. | m.3250T>C (VUS) | 6/6 | 3/6 | 3/6 |
| Tofilo et al., 2026 doi:10.1007/s10815-026-04012-y Table 1; Results Одна носительница; один из шести биоптатов проходит оба критерия. | m.14487T>C | 3/6 | 3/6 | 1/6 |
Sallevelt 2017: 3/22 числовых записей ниже 18% включают низкое значение у 3PN-эмбриона. Для 2PN подтверждённый диагностический выход — 2/23. Четыре анализа неубедительны из-за положительного отрицательного контроля; при принятой классификации остальных 19 записей доля ниже 18% ограничена 2/23–6/23. Это границы неполных данных, не доверительный интервал. Sallevelt et al., 2017 doi:10.1093/humrep/dew356 — pp. 700–701; Supplementary Table SI.