Клональный гемопоэз долго определяли по наличию известных драйверных мутаций вроде DNMT3A или TET2. Но недавние полногеномные исследования одиночных клеток показали, что большинство клональных разрастаний в стареющей крови происходит без этих канонических драйверов. Вокруг данной работы в архиве лежат исследования, которые используют машинное обучение для предсказания болезней по обычному анализу крови, реконструируют генеалогию опухолевых клонов и описывают эволюцию как лоскутное одеяло. Наша статья стоит на отшибе: она не предсказывает болезнь и не строит дерево, а вскрывает наследуемую архитектуру самого процесса накопления мутаций. Это шаг назад от клиники к фундаменту — но шаг, который объясняет, почему одни люди накапливают клоны быстрее других.
Авторы взяли полногеномные последовательности крови 791 067 человек из двух биобанков: 481 079 из UK Biobank и 309 988 из All of Us. Вместо поиска драйверных мутаций они подсчитали «пассажирские» — одиночные соматические варианты, которые встречаются только у одного человека и почти наверняка нейтральны. Чтобы отсеять технический шум, оставили только варианты с поддержкой 4–8 прочтений из 30-кратного покрытия, что соответствует клону размером в четверть или половину клеток крови. Затем с помощью регрессии убрали влияние возраста, пола и курения, получив количественный признак «груза пассажирских мутаций». Для этого признака провели полногеномный поиск ассоциаций (GWAS) у 481 079 участников UK Biobank, редковариантный анализ (RVAS) и стратифицированные по полу и происхождению анализы. Гены-кандидаты приоритизировали с помощью эмпирического байесовского подхода, который объединяет патогенность вариантов, данные об открытом хроматине в стволовых клетках крови и колокализацию с экспрессионными локусами.
Мета-анализ двух когорт выявил 81 локус, ассоциированный с грузом пассажирских мутаций на уровне значимости P < 5×10⁻⁸. Из них 39 совпали с известными локусами клонального гемопоэза, включая TERT, TCL1A, ATM, CHEK2 и BCL2, а 42 оказались новыми. Среди новых — NFE2, главный фактор транскрипции эритроидного и мегакариоцитарного созревания; MSANTD2, регулятор адгезии стволовых клеток; PTPN7, фосфатаза пути MAPK; и CD163, макрофагальный рецептор, связанный с противовоспалительной поляризацией. Редковариантный анализ выявил 11 генов-кандидатов, среди которых MBD2 — белок, читающий метильные метки ДНК. Вредные варианты в MBD2 повышали груз мутаций (β = 0.353, P = 10⁻¹⁸·⁵¹). Это замыкает ось метилирования: DNMT3A пишет метки, TET2 стирает, а MBD2 интерпретирует. Стратификация по полу дала 14 новых локусов, включая женский RMND5B и мужской OLIG3. Анализ по происхождению в All of Us выявил африканский локус KAZN/TMEM51 и восточноазиатский PTPRN2 с большим эффектом (β = 0.538). Наблюдаемый груз мутаций сильнее ассоциировался с гематологическими злокачественными новообразованиями и смертностью, чем полигенный риск, а пространственный анализ обнаружил географическую кластеризацию, не объяснимую доходом, качеством воздуха или химиотерапией.
Работа отличается от соседних исследований тем, что использует количественный признак пассажирских мутаций вместо бинарного наличия драйверов. Это позволило найти больше локусов при меньшем размере выборки, чем в предыдущих драйвер-центричных GWAS. Впервые показана конвергенция общих и редких герминативных вариантов с соматическими драйверами на одних и тех же биологических осях: метилирование ДНК, сплайсинг РНК, репарация ДНК и сигналинг mTORC1. Обнаружение MBD2 как редкого герминативного фактора расширяет ось метилирования с «пишущих» и «стирающих» ферментов на «читающие» белки. Пространственная гетерогенность, сохраняющаяся после учёта известных факторов, указывает на неизмеренные средовые воздействия, которые могут быть столь же важны, как и генетика.
Тема может развиваться в сторону интеграции более богатых данных об экспозиции: радиация, мелкодисперсные частицы, доступ к медицинской помощи. По мере роста неевропейских выборок в биобанках станет возможным проверить, насколько найденные локусы переносятся на другие популяции. Одноклеточные технологии, такие как ReDeeM, позволят связать герминативные варианты с конкретными субпопуляциями стволовых клеток и их клональной динамикой. Также перспективно объединение полигенного риска и наблюдаемого груза мутаций в единый инструмент стратификации риска острого миелоидного лейкоза.
Результаты повлияют на генетическую эпидемиологию старения крови и на разработку инструментов стратификации риска гематологических злокачественных новообразований. Они также дают новые мишени для экспериментальной проверки в биологии стволовых клеток.
Следующие шаги включают функциональную проверку новых генов-кандидатов в клеточных моделях гемопоэза и расширение феномных ассоциаций на неевропейские группы. Также необходимо идентифицировать конкретные средовые факторы, ответственные за географическую кластеризацию.
В области клонального гемопоэза остаётся нерешённым вопрос о том, какие механизмы обеспечивают селективное преимущество клонов без канонических драйверных мутаций. Работа напрямую связывает эту проблему с герминативной изменчивостью в эпигенетических регуляторах и сигнальных путях, предлагая конкретные гены для проверки. Другой нерешённый вопрос — почему клональный гемопоэз по-разному проявляется у мужчин и женщин и в разных популяциях; стратифицированные анализы дают первые зацепки, но требуют больших выборок.
Ген MBD2, найденный в этой работе, кодирует белок, который читает метильные метки на ДНК. Если DNMT3A — это писатель, а TET2 — ластик, то MBD2 — читатель, который интерпретирует написанное. Вредные варианты в этом «читателе» повышают груз мутаций так же, как соматические мутации в «писателе» или «ластике».
Идея о том, что нейтральные «пассажирские» мутации могут служить часами, отсчитывающими возраст клона, перекликается с концепцией молекулярных часов из романа Грега Игана «Город перестановок», где возраст клеток определяется накоплением случайных изменений.
Работа опубликована не на arXiv. Здесь её страница, PDF и формализация: проверить можно самим, а не только прочитать наш разбор.