Данное исследование было одобрено Институциональным наблюдательным советом Второй аффилированной больницы Нанчанского университета (CDEFYYLK 3-05). До начала сбора образцов от всех пациентов было получено письменное информированное согласие в соответствии с Хельсинкской декларацией. Необработанные данные были экспортированы и проанализированы с использованием программного обеспечения и веб-ресурсов, перечисленных в Таблица материалов.
Сбор клинических образцов тканей
Парные образцы опухоли и соответствующей отдаленной нормальной ткани легкого были собраны у пациентов, перенесших хирургическую резекцию в период с января по июнь 2024 года с патологически подтвержденным диагнозом «инвазивная аденокарцинома легкого». Образцы извлекались оперирующим хирургом, обрабатывались в течение 5 min после резекции, немедленно подвергались шоковой заморозке в жидком азоте и хранились до начала анализа. Обезличенные клинико-патологические характеристики представлены в Дополнительной таблице 1.
Внутренняя экстрактивная масс-спектрометрия с электрораспылением (iEESI-MS) для дифференциации LUAD и нормальных тканей легких
Образцы тканей (1 mm3) анализировали методом iEESI-MS с использованием гибридного масс-спектрометра с линейной ионной ловушкой и электростатической орбитальной ловушкой в режиме положительных ионов. Растворитель для экстракции (метанол 100%) подавали со скоростью 3,0 µL/min при смещении +5 kV. Диапазон сканирования iEESI-MS1 составлял m/z 50–2,000. Для iEESI-MS2 использовали коллизионную диссоциацию с нормированной энергией соударений 25–35%, окном изоляции m/z 2,0 и временем динамического исключения 30 s. Ежедневная внешняя калибровка в соответствии со стандартной процедурой прибора обеспечивала погрешность определения массы <5 ppm. После медианной нормализации, log₂-преобразования и автоскейлинга проводили OPLS-DA с 200 пермутационными тестами. Аннотирование метаболитов осуществляли с помощью HMDB, а анализ обогащения метаболических путей проводили с использованием аннотаций KEGG.
Обработка данных масс-спектрометрии и идентификация метаболитов
Необработанные данные масс-спектрометрии (МС) были нормализованы методом медианной нормализации с последующим логарифмическим преобразованием (log₂) и автоскейлингом (центрирование по среднему и деление на стандартное отклонение каждой переменной). Для идентификации дискриминантных признаков m/z использовали ортогональный анализ дискриминантных функций на основе метода частичных наименьших квадратов (OPLS-DA). Признаки считались дифференциальными при VIP > 1.0, |log₂(FC)| > 0.58 и p < 0.05. Дифференциальные признаки подвергали анализу МС2 и аннотировали с помощью базы данных Human Metabolome Database. Затем для выявления измененных путей и приоритизации соответствующих им ферментов использовали анализ обогащения метаболических путей на основе KEGG. Аннотированные дифференциальные метаболиты представлены в Дополнительной таблице 2.
Анализ экспрессии SHMT2, показателей выживаемости и клинико-патологических ассоциаций на основе данных TCGA
Профили экспрессии RNA-seq и соответствующие клинические аннотации для TCGA-LUAD были загружены из Genomic Data Commons 18 декабря 2025 года; опухоли (n = 542) и нормальные ткани (n = 59). Значения количества фрагментов на килобазу транскрипта на миллион картированных прочтений (FPKM) были преобразованы как log2(FPKM + 1). В прогностический анализ были включены только образцы с полными данными о выживаемости. Значения экспрессии SHMT2 в опухолевых и нормальных тканях сравнивали с помощью рангового суммарного критерия Вилкоксона. Пациенты были разделены на группы с высокой и низкой экспрессией SHMT2 по медианному значению экспрессии. Общую выживаемость оценивали с помощью анализа Каплана-Мейера и лог-рангового критерия. Связь со стадией опухоли по системе TNM (tumor-node-metastasis), возрастом и полом оценивали с помощью критерия хи-квадрат или критерия Вилкоксона, в зависимости от ситуации. Зависимые от времени ROC-кривые (receiver operating characteristic) были построены для 1-, 3- и 5-летней выживаемости. Все статистические тесты были двухсторонними.
Анализ дифференциальной экспрессии и функциональное обогащение (GO/KEGG)
Дифференциально экспрессируемые гены (ДЭГ) между группами с высоким и низким уровнем экспрессии SHMT2 были определены с использованием ложноположительной частоты (FDR) < 0.05 и |log₂(FC)| ≥ 1. Результаты были визуализированы в виде графика «вулкан» и тепловой карты 40 наиболее значимых ДЭГ с Z-нормализацией по строкам. Для анализа обогащения Gene Ontology (GO) и KEGG использовалось скорректированное значение p < 0.05. Аннотации генов человека были получены из базы данных Ensembl в декабре 2025 года.
Иммунная инфильтрация, прогнозирование ответа на иммунотерапию и вывод о чувствительности к препаратам
Показатели микроокружения опухоли, включая ImmuneScore, StromalScore и ESTIMATEScore, были рассчитаны с помощью ESTIMATE. Состав иммунных клеток определяли с помощью CIBERSORT с использованием эталонного набора LM22 (загружен 5 декабря 2025 г.), 1000 перестановок и сохранением образцов с p < 0.05 для деконволюции. Погрупповые сравнения проводились между опухолями с высоким и низким уровнем экспрессии SHMT2. Корреляции Спирмена между SHMT2 и генами, связанными с иммунными контрольными точками, были визуализированы в виде корреляционной матрицы. Иммунофеноскоры были получены из TCIA, а TIDE использовался для дополнительного прогнозирования ответа на иммунотерапию. Чувствительность к препаратам определяли с помощью oncoPredict, обученного на выпуске GDSC от декабря 2023 г., для оценки значений IC50. Для анализа межгрупповых различий использовали критерий Вилкоксона, а для оценки взаимосвязей непрерывных переменных — корреляцию Спирмена.
Анализ РНК-секвенирования единичных клеток при LUAD
Данные транскриптомики единичных клеток19 были загружены из базы данных Gene Expression Omnibus 20 декабря 2025 г. и проанализированы. Были отобраны клетки с количеством выявленных признаков от 200 до 6 000 и содержанием митохондриальных транскриптов < 10%. Нормализация количества прочтений проводилась с помощью LogNormalize с коэффициентом масштабирования 10 000; были выбраны 2 000 наиболее вариабельных генов, после чего все гены были масштабированы. Для анализа главных компонент использовались первые 20 компонент, для UMAP — измерения 1–20, а для кластеризации — разрешение 0,5. Аннотация типов клеток выполнялась с использованием референсного атласа Human Primary Cell Atlas. Пропорции типов клеток сравнивались между опухолевыми и нормальными образцами, также были построены матрицы корреляции Спирмена. Распределение SHMT2 изучалось в различных аннотированных типах клеток. Состояния клеточного цикла определялись с помощью cyclone и визуализировались с помощью ridge-графиков.
Культура клеток
Клетки BEAS-2B, H1299 и A549 культивировали при 37 °C в увлажненном инкубаторе при содержании 5% CO22 и относительной влажности воздуха приблизительно 95%, используя среды и добавки, указанные в Список материалов.
Количественная полимеразная цепная реакция с обратной транскрипцией в реальном времени (ОТ-кПЦР)
Общую РНК экстрагировали с использованием реагента для экстракции РНК на основе фенола. Экстрагированную РНК подвергали обратной транскрипции в кДНК с помощью мастер-микса для обратной транскрипции. Обратную транскрипцию проводили в 20 µL система с удалением геномной ДНК в 42 °C в течение 2 мин, обратная транскрипция при 50 °C в течение 15 мин, и прекращение реакции при 85 °C в течение 5 с. Количественную ПЦР проводили с использованием готовой смеси для кПЦР на основе ДНК-связывающего красителя на системе ПЦР в реальном времени с использованием 95 °C в течение 3 мин, затем 40 циклов: 95 °C в течение 10 с и 60 °C в течение 30 с, а затем этап построения кривой плавления с 95 °C в течение 15 с, 60 °C в течение 60 с, и 95 °C в течение 15 с. Все реакции проводили в трех повторностях. Относительную экспрессию рассчитывали с использованием 2−ΔΔCt метод с использованием GAPDH в качестве референсного гена. Стоковые растворы праймеров готовили с концентрацией 10 µM в дважды дистиллированной воде и хранится при -20 °C в темноте. Последовательности праймеров приведены в Дополнительная таблица 3.
Вестерн-блоттинг
Общий клеточный белок экстрагировали с использованием буфера для радиоиммунопреципитационного анализа на льду в течение 30 min, после чего центрифугировали при 12,000 x g в течение 15 min при 4 °C. Белки разделяли методом 12% электрофореза в полиакриламидном геле с додецилсульфатом натрия и переносили на поливинилиденовые фторидные мембраны с размером пор 0.22 µm. Мембраны блокировали 5% (w/v) раствором обезжиренного молока в трис-солевом буфере с 0.1% Tween 20 (TBST) в течение 2 h при комнатной температуре, затем инкубировали в течение ночи при 4 °C с первичными антителами к anti-SHMT2 (1:1,000) и anti-GAPDH (1:2,000), разведенными в 5% обезжиренном молоке в TBST. После трех промывок по 8 min каждой в TBST (20 mM Tris-HCl, 150 mM NaCl, 0.1% Tween-20) мембраны инкубировали с вторичными антителами, конъюгированными с пероксидазой хрена (1:5,000), в течение 1 h при комнатной температуре. Сигналы визуализировали с использованием субстрата для усиленной хемилюминесценции, приготовленного путем смешивания равных объемов растворов A и B; инкубацию проводили в течение 2 min при комнатной температуре в темноте, после чего изображения фиксировали с помощью системы хемилюминесцентной визуализации.
Валидация иммуногистохимического анализа HPA
Профили экспрессии белка SHMT2 в нормальной ткани легкого и при LUAD оценивали с помощью HPA. В базе HPA был выполнен поиск по запросу «SHMT2», после чего для сравнения были загружены репрезентативные изображения иммуногистохимического (ИГХ) анализа и соответствующие патоморфологические аннотации. Чтобы минимизировать систематическую ошибку, связанную со специфичностью антител, приоритет отдавался антителам с более высоким уровнем валидации в HPA. При наличии данных результаты ИГХ с использованием нескольких антител подвергались перекрестной проверке (HPA020543 и HPA020549).
Нокдаун SHMT2 в клетках LUAD
Для изучения внутренней роли SHMT2 в опухолевых клетках при LUAD был выполнен нокдаун SHMT2 в клетках H1299 и A549 с использованием shRNA. shRNA переносился лентивирусным вектором с промотором H1 и селективным маркером пуромицина; при производстве вируса совместно трансфецировали упаковочные плазмиды. sh-SHMT2 был направлен на кодирующую область человеческого гена SHMT2, в то время как sh-NC служил скрамблированным отрицательным контролем. Подготовленные вирусные стоки имели титр ≥ 1 x 108 трансдуцирующих единиц/мл и хранились при -80 °C. Клетки заражали подготовленными лентивирусными частицами sh-SHMT2 или sh-NC, описанными в таблице материалов. Через 48 ч эффективность нокдауна оценивали с помощью RT-qPCR и вестерн-блоттинга. Для анализов пролиферации, миграции и колониеобразования использовали только клетки с подтвержденным сайленсингом SHMT2.
Анализ пролиферации клеток
Пролиферацию клеток оценивали с помощью анализа CCK-8. Трансфицированные клетки H1299 и A549 высевали в 96-луночные планшеты в количестве 5 x 103 клеток на лунку в шестикратной повторности. На 1, 2, 3, 4 и 5 дни в каждую лунку добавляли 10 µL реагента CCK-8 (10% от объема культуральной среды) и инкубировали в течение 2 ч при 37 °C. Оптическую плотность измеряли при 450 nm и использовали в качестве показателя жизнеспособности клеток.
Анализ заживления ран
Миграцию клеток оценивали с помощью теста на заживление раны (wound-healing assay). Трансфицированные клетки H1299 и A549 высевали в 6-луночные планшеты и культивировали до достижения 90–95% конфлюэнтности. С помощью стерильного наконечника пипетки объемом 200 µL создавали линейную «рану», после чего отделившиеся клетки удаляли осторожным промыванием фосфатно-солевым буферным раствором. Затем клетки содержали в безсывороточной или низкосывороточной среде, а область раны фотографировали через 0 h и 24 h с помощью инвертированного микроскопа. Степень закрытия раны количественно определяли с использованием программного обеспечения для анализа изображений, указанного в таблице материалов, и рассчитывали как процент сокращения площади относительно начальной площади раны.
Анализ колониеобразующей способности
Для анализа колониеобразования трансфицированные клетки H1299 и A549 высевали в 6-луночные планшеты в количестве 800 клеток на лунку и культивировали в течение примерно 10–14 дней до образования видимых колоний, осуществляя микроскопический мониторинг в процессе культивирования. Колонии фиксировали 4% параформальдегидом в течение 30 мин при комнатной температуре и окрашивали 0,1% кристаллическим фиолетовым в течение 15 мин при комнатной температуре. Под световым микроскопом подсчитывали колонии, содержащие более 50 клеток.
Статистический анализ для экспериментов in vitro
Все эксперименты in vitro проводились независимо как минимум трижды. Данные представлены в виде среднего значения ± стандартного отклонения. Для сравнения двух групп использовался t-критерий Стьюдента; значение p < 0,05 считалось статистически значимым.