Профиль экспрессии клеточных микроРНК и их потенциальная роль при раке шейки матки

Обложка


Цитировать

Полный текст

Открытый доступ Открытый доступ
Доступ закрыт Доступ предоставлен
Доступ закрыт Доступ платный или только для подписчиков

Аннотация

Обоснование. Рак шейки матки (РШМ) занимает 4-е место в структуре заболеваемости и смертности от онкологических заболеваний у женщин. Несмотря на установленный этиологический фактор РШМ, а именно инфекцию вирусами папиллом человека (HPV, human papillomavirus) высокого онкогенного риска и имеющуюся вакцинопрофилактику, существует проблема как в понимании механизма HPV-ассоциированного канцерогенеза, так и в разработке новых диагностических и терапевтических подходов для РШМ. Важную роль в патогенезе РШМ играет эпигенетическая регуляция экспрессии генов с помощью микроРНК, поэтому актуален поиск новых дифференциально экспрессирующихся микроРНК, потенциально задействованных в механизмах HPV-ассоциированной трансформации опухолевых клеток.

Цель. Поиск микроРНК, дифференциально экспрессирующихся в опухолевой ткани РШМ, и in silico определение их потенциальной функциональной роли.

Материалы и методы. Спектр микроРНК, экспрессия которых изменяется в опухолевой ткани плоскоклеточных HPV16-положительных карцином шейки матки, определяли с использованием NGS-секвенирования (Next Generation Sequencing) опухолевой ткани и условно-нормального прилежащего эпителия, полученных с помощью микродиссекции. Потенциальные гены-мишени исследуемых микроРНК выявляли с помощью сервисов MirTargetLink, Tarbase v9.0, LinkedOmics. Функциональное обогащение генов выполняли с помощью Metascape. Поиск ассоциаций уровня экспрессии микроРНК с клиническими параметрами при РШМ (по данным TCGA, выборка CESC) вели с использованием сервиса USCS Xena.

Результаты. По данным NGS-секвенирования парных HPV16-положительных образцов опухолевой и нормальной ткани шейки матки мы определили 42 дифференциально экспрессируемые микроРНК, среди которых уровень 22 повышен и 20 снижен в опухолевой ткани по сравнению с условно-нормальным прилежащим эпителием. Анализ потенциальных мишеней наиболее значимых микроРНК выявил множество функциональных категорий генов, потенциально вовлечённых в процессы канцерогенеза, а также ассоциацию с клиническими характеристиками. Мы показали, что высокий уровень экспрессии микроРНК-20b в опухолевой ткани коррелирует с риском возникновения отдалённых метастазов, в то время как сниженный уровень микроРНК-218-1 и микроРНК-218-2 связан с неблагоприятным прогнозом заболевания. Для микроРНК-363, -615, -769 впервые описано их увеличение при РШМ и определены потенциальные мишени и сигнальные пути, ассоциированные с уровнем их экспрессии.

Заключение. Поиск дифференциально экспрессирующихся микроРНК при РШМ выявил спектр микроРНК, которые потенциально играют важную роль в процессах онкотрансформации и поддержании опухолевого фенотипа. Помимо микроРНК, функциональное значение которых описано в мировой литературе, мы обнаружили микроРНК, роль которых при РШМ неизвестна и для которых определены потенциально значимые мишени, важные для понимания механизмов канцерогенеза. Полученные данные могут лечь в основу разработки новых подходов к диагностике и терапии РШМ.

Полный текст

Доступ закрыт

Об авторах

Данила Сергеевич Елкин

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Автор, ответственный за переписку.
Email: d.elkin@ronc.ru
ORCID iD: 0000-0002-4793-6063
SPIN-код: 9946-6863
Россия, Москва

Радик Сяитович Фасхутдинов

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: r.faskhutdinov@ronc.ru
ORCID iD: 0000-0002-0050-7798
Россия, Москва

Ольга Владимировна Зверева

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: zvereva.owl@mail.ru
ORCID iD: 0009-0007-7691-3235
Россия, Москва

Мария Дмитриевна Федорова

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: m.d.fedorova@ronc.ru
ORCID iD: 0000-0002-8813-7516
SPIN-код: 4943-5931

канд. биол. наук

Россия, Москва

Алексей Николаевич Катаргин

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: a.katargin@ronc.ru
ORCID iD: 0000-0002-7405-0671

канд. биол. наук

Россия, Москва

Лариса Сергеевна Павлова

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: l.pavlova@ronc.ru
ORCID iD: 0000-0003-3993-4823
Россия, Москва

Кирилл Иосифович Жорданиа

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: k.zhordania@ronc.ru
ORCID iD: 0000-0003-1380-3710
SPIN-код: 6271-8954

д-р мед. наук, профессор

Россия, Москва

Екатерина Александровна Мустафина

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: e.mustafina@ronc.ru
ORCID iD: 0000-0002-1009-0383

канд. мед. наук

Россия, Москва

Светлана Владимировна Винокурова

Национальный медицинский исследовательский центр онкологии им. Н.Н. Блохина

Email: s.vinokurova@ronc.ru
ORCID iD: 0000-0003-1615-3928
SPIN-код: 3453-4502

канд. мед. наук

Россия, Москва

Список литературы

  1. Bray F., Laversanne M., Sung H., et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries // CA Cancer J Clin. 2024. Vol. 74, N 3. P. 229–263. doi: 10.3322/caac.21834
  2. Egawa N., Egawa K., Griffin H., Doorbar J. Human Papillomaviruses; Epithelial Tropisms, and the Development of Neoplasia // Viruses. 2015. Vol. 7, N 7. P. 3863–3890. doi: 10.3390/v7072802
  3. Scarth J., Patterson M., Morgan E., Macdonald A. The human papillomavirus oncoproteins: a review of the host pathways targeted on the road to transformation // J Gen Virol. 2021. Vol. 102, N 3. P. 001540. doi: 10.1099/jgv.0.001540
  4. Thierry F. Transcriptional regulation of the papillomavirus oncogenes by cellular and viral transcription factors in cervical carcinoma // Virology. 2009. Vol. 384, N 2. P. 375–379. doi: 10.1016/j.virol.2008.11.014
  5. Pett M., Coleman N. Integration of high-risk human papillomavirus: a key event in cervical carcinogenesis? // J Pathol. 2007. Vol. 212, N 4. P. 356–367. doi: 10.1002/path.2192
  6. Vinokurova S., Wentzensen N., Kraus I., et al. Type-dependent integration frequency of human papillomavirus genomes in cervical lesions // Cancer Res. 2008. Vol. 68, N 1. P. 307–313. doi: 10.1158/0008-5472.CAN-07-2754
  7. Arias-Pulido H., Peyton C., Joste N., et al. Human papillomavirus type 16 integration in cervical carcinoma in situ and in invasive cervical cancer // J Clin Microbiol. 2006. Vol. 44, N 5. P. 1755–1762. doi: 10.1128/JCM.44.5.1755-1762.2006
  8. Klaes R., Friedrich T., Spitkovsky D., et al. Overexpression of p16(INK4A) as a specific marker for dysplastic and neoplastic epithelial cells of the cervix uteri // Int J Cancer. 2001. Vol. 92, N 2. P. 276–284. doi: 10.1002/ijc.1174
  9. Sano T., Oyama T., Kashiwabara K., et al. Expression status of p16 protein is associated with human papillomavirus oncogenic potential in cervical and genital lesions // Am J Pathol. 1998. Vol. 153, N 6. P. 1741–1748. doi: 10.1016/S0002-9440(10)65689-1
  10. Shang R., Lee S., Senavirathne G., Lai E. microRNAs in action: biogenesis, function and regulation // Nat Rev Genet. 2023. Vol. 24, N 12. P. 816–833. doi: 10.1038/s41576-023-00611-y
  11. Chauhan P., Pramodh S., Hussain A., et al. Understanding the role of miRNAs in cervical cancer pathogenesis and therapeutic responses // Front Cell Dev Biol. 2024. Vol. 12. P. 1397945. doi: 10.3389/fcell.2024.1397945
  12. Robinson M., McCarthy D., Smyth G. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data // Bioinformatics. 2010. Vol. 26, N 1. P. 139–140. doi: 10.1093/bioinformatics/btp616
  13. Goedhart J., Luijsterburg M. VolcaNoseR is a web app for creating, exploring, labeling and sharing volcano plots // Sci Rep. 2020. Vol. 10, N 1. P. 20560. doi: 10.1038/s41598-020-76603-3
  14. Xu F., Wang Y., Ling Y., et al. dbDEMC 3.0: Functional Exploration of Differentially Expressed miRNAs in Cancers of Human and Model Organisms // Genomics Proteomics Bioinformatics. 2022. Vol. 20, N 3. P. 446–454. doi: 10.1016/j.gpb.2022.04.006
  15. Tang D., Chen M., Huang X., et al. SRplot: A free online platform for data visualization and graphing // PLoS One. 2023. Vol. 18, N 11. P. e0294236. doi: 10.1371/journal.pone.0294236
  16. Kern F., Aparicio-Puerta E., Li Y., et al. miRTargetLink 2.0-interactive miRNA target gene and target pathway networks // Nucleic Acids Res. 2021. Vol. 49, N W1. P. W409–W416. doi: 10.1093/nar/gkab297
  17. 17. Skoufos G., Kakoulidis P., Tastsoglou S., et al. TarBase-v9.0 extends experimentally supported miRNA-gene interactions to cell-types and virally encoded miRNAs // Nucleic Acids Res. 2024. Vol. 52, N D1. P. D304–D310. doi: 10.1093/nar/gkad1071
  18. Vasaikar S., Straub P., Wang J., Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types // Nucleic Acids Res. 2018. Vol. 46, N D1. P. D956–D963. doi: 10.1093/nar/gkx1090
  19. Zhou Y., Zhou B., Pache L., et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets // Nat Commun. 2019. Vol. 10, N 1. P. 1523. doi: 10.1038/s41467-019-09234-6
  20. Goldman M., Craft B., Hastie M., et al. Visualizing and interpreting cancer genomics data via the Xena platform // Nat Biotechnol. 2020. Vol. 38, N 6. P. 675–678. doi: 10.1038/s41587-020-0546-8
  21. Cruz-De la Rosa M., Jimenez-Wences H., Alarcon-Millan J., et al. miR-218-5p/RUNX2 Axis Positively Regulates Proliferation and Is Associated with Poor Prognosis in Cervical Cancer // Int J Mol Sci. 2022. Vol. 23, N 13. doi: 10.3390/ijms23136993
  22. Zhu L., Tu H., Liang Y., Tang D. MiR-218 produces anti-tumor effects on cervical cancer cells in vitro // World J Surg Oncol. 2018. Vol. 16, N 1. P. 204. doi: 10.1186/s12957-018-1506-3
  23. Wang P., Zhai G., Bai Y. Values of miR-34a and miR-218 expression in the diagnosis of cervical cancer and the prediction of prognosis // Oncol Lett. 2018. Vol. 15, N 3. P. 3580–3585. doi: 10.3892/ol.2018.7791
  24. Kogo R., How C., Chaudary N., et al. The microRNA-218~Survivin axis regulates migration, invasion, and lymph node metastasis in cervical cancer // Oncotarget. 2015. Vol. 6, N 2. P. 1090–1100. doi: 10.18632/oncotarget.2836
  25. Yamamoto N., Kinoshita T., Nohata N., et al. Tumor suppressive microRNA-218 inhibits cancer cell migration and invasion by targeting focal adhesion pathways in cervical squamous cell carcinoma // Int J Oncol. 2013. Vol. 42, N 5. P. 1523–1532. doi: 10.3892/ijo.2013.1851
  26. Szekerczes T., Galamb A., Varga N., et al. Increased miR-20b Level in High Grade Cervical Intraepithelial Neoplasia // Pathol Oncol Res. 2020. Vol. 26, N 4. P. 2633–2640. doi: 10.1007/s12253-020-00852-w
  27. Cheng Y., Geng L., Zhao L., et al. Human papillomavirus E6-regulated microRNA-20b promotes invasion in cervical cancer by targeting tissue inhibitor of metalloproteinase 2 // Mol Med Rep. 2017. Vol. 16, N 4. P. 5464–5470. doi: 10.3892/mmr.2017.7231
  28. Lajer C., Garnaes E., Friis-Hansen L., et al. The role of miRNAs in human papilloma virus (HPV)-associated cancers: bridging between HPV-related head and neck cancer and cervical cancer // Br J Cancer. 2012. Vol. 106, N 9. P. 1526–1534. doi: 10.1038/bjc.2012.109
  29. Wald A., Hoskins E., Wells S., et al. Alteration of microRNA profiles in squamous cell carcinoma of the head and neck cell lines by human papillomavirus // Head Neck. 2011. Vol. 33, N 4. P. 504–512. doi: 10.1002/hed.21475
  30. Haga R., Ridley A. Rho GTPases: Regulation and roles in cancer cell biology // Small GTPases. 2016. Vol. 7, N 4. P. 207–221. doi: 10.1080/21541248.2016.1232583
  31. Castro-Munoz L., Rocha-Zavaleta L., Lizano M., et al. Alteration of the IFN-Pathway by Human Papillomavirus Proteins: Antiviral Immune Response Evasion Mechanism // Biomedicines. 2022. Vol. 10, N 11. P. 2965. doi: 10.3390/biomedicines10112965

Дополнительные файлы

Доп. файлы
Действие
1. JATS XML
2. Рис. 1. Результаты NGS-секвенирования. a. Диаграмма Вулкана, демонстрирующая паттерн экспрессии микроРНК в опухолевой ткани рака шейки матки. Подписями обозначены микроРНК, имеющие наиболее значимое изменение экспрессии как по степени выраженности, так и по уровню значимости. b, c. Диаграммы Венна, отображающие пересечение списков дифференциально экспрессирующихся микроРНК при раке шейки матки по данным TCGA и собственных данных в зависимости от того, снижен (b) или повышен (c) уровень микроРНК в опухолевой ткани.

Скачать (325KB)
3. Рис. 2. Схематическое представление алгоритма поиска генов-мишеней микроРНК и их обогащения по функциональным категориям. Жёлтым цветом на диаграмме Венна схематично представлены гены, используемые в нижележащем анализе. ДЭ-микроРНК — дифференциально экспрессируемые микроРНК.

Скачать (195KB)
4. Рис. 3. Хордовая диаграмма, отражающая функциональные категории, по которым обогащаются гены-мишени исследуемых микроРНК. На левой половине диаграммы отображены обозначения категорий, на правой половине различными цветами указаны микроРНК.

Скачать (880KB)
5. Рис. 4. Диаграмма Сэнкей, отображающая наиболее общие функциональные категории генов-мишеней (левая часть диаграммы), соответствующие регулирующим их микроРНК (правая часть диаграммы).

Скачать (515KB)
6. Риc. 5. Схематичное отображение функциональных категорий генов-мишеней соответствующих микроРНК с повышенной экспрессией в опухолевой ткани рака шейки матки. Категории расположены по оси ординат в порядке убывания в соответствии с уровнем значимости участия генов-мишеней соответствующих микроРНК (–log10 <2 по оси абсцисс).

Скачать (775KB)
7. Рис. 6. Схематичное отображение функциональных категорий генов-мишеней соответствующих микроРНК с пониженной экспрессией в опухолевой ткани рака шейки матки. Категории расположены по оси ординат в порядке убывания в соответствии с уровнем значимости участия генов-мишеней соответствующих микроРНК (–log10 <2 по оси абсцисс).

Скачать (655KB)
8. Рис. 7. Зависимость уровня экспрессии микроРНК-20b (a), микроРНК-345 (b), микроРНК-363 (c) от наличия отдалённых метастазов. *p <0,05.

Скачать (127KB)
9. Рис. 8. Зависимость экспрессии микроРНК-106b со степенью злокачественности опухолевого процесса. * p <0,05, ** p <0,01.

Скачать (71KB)
10. Рис. 9. Корреляция уровня экспрессии микроРНК-218-1 с общей выживаемостью пациентов с раком шейки матки (a), стадией заболевания (b), размером первичного очага опухоли (с). Корреляция уровня экспрессии микроРНК-218-2 с общей выживаемостью пациентов с раком шейки матки (d), стадией заболевания (e), размером первичного очага опухоли (f). * p <0,05, ** p <0,01, *** p <0,001.

Скачать (90KB)

© Эко-Вектор, 2024



СМИ зарегистрировано Федеральной службой по надзору в сфере связи, информационных технологий и массовых коммуникаций (Роскомнадзор).
Регистрационный номер и дата принятия решения о регистрации СМИ: ПИ № ФС 77 - 86496 от 11.12.2023 г
СМИ зарегистрировано Федеральной службой по надзору в сфере связи, информационных технологий и массовых коммуникаций (Роскомнадзор).
Регистрационный номер и дата принятия решения о регистрации СМИ: ЭЛ № ФС 77 - 80673 от 23.03.2021 г
.