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

Обложка

Цитировать

Полный текст

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

Аннотация

Обоснование. Рак шейки матки (РШМ) занимает 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


 


Согласие на обработку персональных данных

 

Используя сайт https://journals.rcsi.science, я (далее – «Пользователь» или «Субъект персональных данных») даю согласие на обработку персональных данных на этом сайте (текст Согласия) и на обработку персональных данных с помощью сервиса «Яндекс.Метрика» (текст Согласия).