Экспрессия FOXD1 повышена в тканях LUSC человека
Запрос в базе данных интерактивного анализа профилей экспрессии генов (GEPIA) показал более высокое содержание транскриптов FOXD1 в образцах LUSC, чем в образцах нормальной ткани легкого, в то время как при сравнении с аденокарциномой легкого значимых различий выявлено не было (Рисунок 1A). Анализ 12 пар образцов LUSC и прилегающих нетрансформированных тканей с помощью qRT-PCR и вестерн-блоттинга показал изменение в том же направлении (Рисунок 1B,C). Иммуногистохимическое окрашивание также продемонстрировало более интенсивный сигнал FOXD1 в срезах опухоли, чем в соответствующих образцах прилегающей ткани (Рисунок 1D).
Нокдаун FOXD1 снижает жизнеспособность и миграцию, а также увеличивает апоптоз в клетках LUSC
Уровень экспрессии транскрипта FOXD1 был выше в клетках SK-MES-1 и NCI-H226, чем в клетках BEAS-2B (Рисунок 2A). Был проведен скрининг трех малых интерферирующих РНК, нацеленных на FOXD1; siRNA-1 обеспечила наибольшее снижение экспрессии, выявленное с помощью qRT-PCR и вестерн-блоттинга, и была выбрана для последующих анализов (Рисунок 2B,C). По сравнению с отрицательным контролем, нокдаун FOXD1 привел к снижению сигнала Cell Counting Kit-8 (Рисунок 2D), увеличению количества Annexin V-положительных клеток (Рисунок 2E) и снижению миграции в системе Transwell (Рисунок 2F) в обеих клеточных линиях LUSC. Для дальнейших исследований был использован один выбранный реагент, поэтому функциональные результаты не являются доказательством исключения нецелевых эффектов, независимым от последовательности.
Экспрессия METTL3 коррелирует с FOXD1 в образцах LUSC
Метилирование m6A представляет собой динамическую и обратимую химическую модификацию, широко распространенную у эукариот, которая регулируется метилтрансферазами и деметилазами, а также распознается и обрабатывается m6A-связывающими белками18. Общий сигнал m6A был выше в 12 образцах тканей LUSC, чем в соответствующих прилежащих тканях, а также выше в обеих клеточных линиях LUSC, чем в клетках BEAS-2B (Рисунок 3A). Экспрессия транскриптов, кодирующих METTL3, METTL14, WTAP и RBM15, была повышена в парных образцах опухолей (Рисунок 3B). Согласно данным GEPIA и анализу набора из 12 опухолей, FOXD1 продемонстрировал наиболее выраженную положительную корреляцию с METTL3; связь с RBM15 была слабее (Рисунок 3C,D). Вестерн-блоттинг дополнительно подтвердил более высокое содержание METTL3 в опухолевых тканях (Рисунок 3E). Данные анализы устанавливают ковариацию, а не прямую регуляцию.
Нокдаун METTL3 изменяет уровень экспрессии FOXD1 и фенотип клеток LUSC
Три малые интерферирующие РНК, нацеленные на METTL3, были протестированы в клетках SK-MES-1 и NCI-H226 с помощью qRT-PCR и вестерн-блоттинга (Рисунок 4A,B). Нокдаун METTL3 привел к снижению общего уровня сигнала m6A (Рисунок 4C). Оверэкспрессия FOXD1 повысила уровень FOXD1 без увеличения количества METTL3 (Рисунок 4D), а эксперимент по рескью с четырьмя группами подтвердил снижение содержания FOXD1 после нокдауна METTL3 и его восстановление с помощью конструкта FOXD1 (Рисунок 4E,F). Нокдаун METTL3 снизил жизнеспособность и миграцию клеток, а также увеличил апоптоз; реэкспрессия FOXD1 частично противодействовала этим изменениям (Рисунок 4G–I). Таким образом, данные по рескью согласуются с тем, что FOXD1 действует ниже METTL3 в культивируемых клетках.
Нарушение функции METTL3 связано со стабильностью РНК и белка FOXD1
Растущее число данных свидетельствует о том, что модификация мРНК m6A модулирует различные аспекты цикла мРНК, включая сплайсинг пре-мРНК, трансляцию, нуклеоцитоплазматический экспорт, деградацию и сплайсинг19,20. С помощью SRAMP были выявлены потенциальные мотивы m6A в РНК FOXD1, а в базе m6A2Target белок METTL3 был указан как предсказанный writer, связанный с FOXD1 в клетках H1299 (Рисунок 5A,B). Результаты этих баз данных позволяют формулировать гипотезы, но не демонстрируют связывание или сайт-специфическое метилирование в клетках LUSC. После блокады транскрипции актиномицином D уровень РНК FOXD1 снижался быстрее после нокдауна METTL3, чем в контрольных клетках (Рисунок 5C). Содержание FOXD1 было выше в клетках с гиперэкспрессией METTL3, обработанных MG132, чем в контрольных клетках с вектором, обработанных MG132, в то время как такая закономерность не наблюдалась после обработки CHX (Рисунок 5D,E). Временной ход обработки CHX показал более быструю потерю белка FOXD1 после нокдауна METTL3 (Рисунок 5F). В совокупности эти анализы согласуются с изменениями стабильности РНК и белка FOXD1, связанными с METTL3, однако они не идентифицируют модифицированный аденозин, не демонстрируют прямого связывания и не определяют reader m6A.
Нокдаун FOXD1 подавляет рост ксенографтов LUSC
Для проверки функции FOXD1 in vivo контрольные клетки SK-MES-1 и клетки SK-MES-1 с нокдауном FOXD1 имплантировали подкожно голым мышам (n = 8 в группе). Начиная с первой недели после имплантации, размеры опухолей измеряли каждые 4 дня в течение 20 дней; сбор опухолей проводили на 27-й день после инокуляции. Удаленные опухоли в группе с нокдауном FOXD1 были меньше (Рисунок 6A), кривые их роста поднимались медленнее (Рисунок 6B), а их конечная масса была ниже (Рисунок 6C). Окрашивание на FOXD1 и Ki-67 в опухолях с нокдауном также было снижено (Рисунок 6D). Данные эксперименты оценивают только нокдаун FOXD1; они не проверяют ингибирование METTL3 in vivo.
В совокупности было установлено, что FOXD1 гиперэкспрессируется в тканях и клетках плоскоклеточного рака легкого и способствует пролиферации, миграции и росту ксенографтов. Нокдаун METTL3 приводил к снижению общего уровня m6A, уменьшению экспрессии FOXD1 и подавлению злокачественного фенотипа клеток, в то время как повторная экспрессия FOXD1 частично обращала эти эффекты. Анализы стабильности подтверждают посттранскрипционную регуляцию FOXD1, связанную с METTL3, но не устанавливают прямое сайт-специфическое метилирование или восстановление всего пути in vivo. Дизайн исследования обобщен на Рисунке 7.
ДОСТУПНОСТЬ ДАННЫХ:
Все данные, проанализированные в этом исследовании, представлены в статье и на ее рисунках. Обезличенные числовые данные, лежащие в основе представленных сводных показателей, могут быть предоставлены автором, ответственным за переписку, по обоснованному запросу при условии соблюдения требований институционального управления и защиты участников. Файлы моделирования или специализированный код для анализа не использовались.

Рисунок 1: Экспрессия FOXD1 в образцах рака легкого человека. (A) Сравнение уровней транскриптов FOXD1 с помощью GEPIA при аденокарциноме легкого (LUAD) и плоскоклеточном раке легкого (LUSC) по сравнению с образцами нормальной ткани легкого. (B) Уровень транскриптов FOXD1 в 12 парах образцов LUSC и прилегающих нетуморных тканей. (C) Уровень белка FOXD1 в тех же парах тканей. (D) Репрезентативные изображения иммуногистохимического окрашивания на FOXD1. Ткани человека n = 12 пар; технические повторы в стиле in vitro не изменяют значение n на уровне пациентов. Масштабная линейка на панели D = 50 µm. ***P < 0.001. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 2: Эффекты нокдауна FOXD1 в клетках SK-MES-1 и NCI-H226. (A) Относительное количество транскриптов FOXD1 в клетках BEAS-2B, SK-MES-1 и NCI-H226. (B) Скрининг трех малых интерферирующих РНК, нацеленных на FOXD1, методом qRT-PCR через 48 ч после трансфекции. (C) Скрининг тех же трех малых интерферирующих РНК, нацеленных на FOXD1, методом вестерн-блоттинга. (D) Анализ с использованием набора для подсчета клеток Cell counting kit-8 после нокдауна FOXD1. (E) Анализ апоптоза с помощью аннексина V, меченного флуоресцеинизотиоцианатом, и пропидия иодида через 48 ч. (F) Миграция в системе Transwell, измеренная через 24 ч после посева трансфецированных клеток; архивные изображения полей зрения представлены без ретроспективного масштабного отрезка, так как метаданные калибровки были недоступны. Проведено три независимых биологических эксперимента. *P < 0.05; **P < 0.01; ***P < 0.001. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 3: Общий уровень N6-метиладенозина и связь между METTL3 и FOXD1. (A) Общий сигнал m6A в парных тканях, а также в клетках BEAS-2B, SK-MES-1 и NCI-H226. (B) Транскрипты METTL3, METTL14, WTAP и RBM15 в 12 парах тканей LUSC и прилежащих тканей. (C) Корреляции GEPIA между FOXD1 и каждым из белков-райтеров. (D) Корреляции между FOXD1 и транскриптами райтеров в 12 опухолях LUSC. (E) Белок METTL3 в парных тканях LUSC и прилежащих тканях. Ткани человека n = 12 пар; клеточные анализы представляют три независимых биологических эксперимента. *P < 0.05; **P < 0.01; ***P < 0.001. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 4: Нокдаун METTL3, спасение экспрессии FOXD1 и фенотипы клеток LUSC. (A) Скрининг трех малых интерферирующих РНК, нацеленных на METTL3, методом qRT-PCR через 48 ч после трансфекции. (B) Скрининг тех же трех малых интерферирующих РНК, нацеленных на METTL3, методом вестерн-блоттинга. (C) Общий сигнал m6A после нокдауна METTL3. (D) Содержание транскрипта и белка FOXD1 после гиперэкспрессии FOXD1. (E) Содержание транскриптов METTL3 и FOXD1 в контрольных условиях, при si-METTL3, ov-FOXD1 и комбинированном воздействии si-METTL3 + ov-FOXD1, измеренное методом qRT-PCR. (F) Содержание белков METTL3 и FOXD1 в тех же группах спасения, измеренное методом вестерн-блоттинга. (G) Анализ с использованием набора Cell Counting Kit-8. (H) Апоптоз, определенный методом проточной цитометрии. (I) Миграция в системе Transwell. Было проведено три независимых биологических эксперимента. Числа под полосами вестерн-блоттинга указывают на нормированные соотношения целевого белка/GAPDH. ns — статистически не значимо; *P < 0.05; **P < 0.01; ***P < 0.001. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 5: METTL3 связан с регуляцией транскрипта и белка FOXD1. (A) Сайты N6-метиладенозина в транскрипте FOXD1, предсказанные с помощью SRAMP. (B) Прогноз m6A2Target, идентифицирующий METTL3 как потенциальный белок-райтер, связанный с FOXD1. (C) Относительное содержание мРНК FOXD1 после обработки актиномицином D в концентрации 5 µg/mL в клетках отрицательного контроля и клетках с нокдауном METTL3 через 0 h, 2 h, 4 h и 6 h. (D) Содержание белка FOXD1 в необработанных клетках с векторным контролем, клетках с векторным контролем, обработанных MG132, и клетках с гиперэкспрессией METTL3, обработанных MG132 (10 µM MG132). (E) Содержание белка FOXD1 в необработанных клетках с векторным контролем, клетках с векторным контролем, обработанных CHX, и клетках с гиперэкспрессией METTL3, обработанных CHX (20 µM CHX). (F) Экспрессия белка FOXD1 после обработки CHX в концентрации 20 µM в течение 0 h, 3 h, 6 h и 9 h в клетках отрицательного контроля и клетках с нокдауном METTL3. GAPDH использовался в качестве контроля нагрузки. Эксперименты на клетках проводились в трех независимых биологических повторностях. Представленные блоты являются репрезентативными; числа под полосами указывают на нормализованные соотношения FOXD1/GAPDH. *P < 0.05; **P < 0.01; ***P < 0.001. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 6: Нокдаун FOXD1 подавляет рост ксенографтов SK-MES-1. (A) Извлеченные опухоли из групп отрицательного контроля и нокдауна FOXD1; масштабная линейка = 1 cm. (B) Рост опухоли, измеряемый каждые 4 дня в течение 20 дней, начиная через 1 неделю после имплантации; сбор материала на 27-й день после инокуляции. (C) Вес опухоли при сборе. (D) Репрезентативные результаты иммуногистохимического окрашивания на FOXD1 и Ki-67 в ткани ксенографта; масштабная линейка = 50 µm. n = 8 мышей в группе; **P < 0.01. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 7: Общая схема исследования и границы доказательной базы. Пары тканей LUSC (n = 12 пар) и клеточные линии использовались для профилирования FOXD1, METTL3 и глобального m6A; за генетическим воздействием и восстановлением экспрессии следовали анализы жизнеспособности, апоптоза, миграции и стабильности. Модель ксенотрансплантата с нокдауном FOXD1 (n = 8 мышей в группе) наблюдалась в течение 20 дней. На диаграмме разграничены данные по METTL3, полученные на клеточных моделях, и специфические для FOXD1 данные in vivo, а также указаны неисследованные аспекты: прямое связывание, сайт-специфический m6A, идентификация ридера и активность METTL3 in vivo. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
| Ген | Направление | Последовательность (5′–3′) |
| FOXD1 | Прямой | TCTGAGTAACGCGTTATGCTT |
| FOXD1 | Обратный | CGGAATCGGACGCAAGTCTT |
| GAPDH | Прямой | GGAGCGAGATCCCTCCAAAAT |
| GAPDH | Обратный | GGCTGTTGTCATACTTCTCATGG |
| METTL14 | Прямой | CTGAAAGTGCCGACAGCATTGG |
| METTL14 | Обратный | CTCTCCTTCATCCAGATACTTACG |
| METTL3 | Прямой | CCAGCACAGCTTCAGCAGTTCC |
| METTL3 | Обратный | GCGTGGAGATGGCAAGACAGATG |
| RBM15 | Прямой | GTTGTGGCTTATGTGGAGTTTAC |
| RBM15 | Обратный | CACTTAAAACACCGGCATTGG |
| WTAP | Прямой | GCAACAACAGCAGGAGTCTGCA |
| WTAP | Обратный | CTGCTGGACTTGCTTGAGGTAC |
Таблица 1: Последовательности праймеров для количественной полимеразной цепной реакции с обратной транскрипцией. Последовательности указаны в направлении 5′–3′.