Ekspresja TMOD3 w raku jajnika
Najpierw, baza danych GEO wykazała, że poziomy ekspresji mRNA TMOD3 były podwyższone w zbiorach danych z mikromacierzy GSE51088 i GSE66957 (Ryc. 1A,B). TMOD3 był również wysoce ekspresyjny w raku jajnika w porównaniu z prawidłowymi tkankami jajnika, co wykazano za pomocą narzędzia internetowego TNMplot (Ryc. 1C). Analiza danych CPTAC za pomocą narzędzia internetowego UALCAN wykazała, że poziom białka TMOD3 był również istotnie wyższy w raku jajnika niż w prawidłowych tkankach jajnika (Ryc. 1D). Ponadto zbiór danych GEO GSE53963 wykazał, że ekspresja TMOD3 była istotnie wyższa w rakach jajnika stopnia 4 niż w rakach stopnia 2 i 3 (Ryc. 1E), a zbiór danych GSE23554 wykazał, że poziomy TMOD3 były wyższe u pacjentek z rakiem jajnika w stadium 3 niż w stadiach 1 i 2 (Ryc. 1F). Dodatkowo zbiór danych GSE131978 wykazał, że ekspresja TMOD3 była również istotnie wyższa w tkankach raka jajnika z przerzutami niż w tkankach bez przerzutów (Ryc. 1G). Wyniki te sugerują, że TMOD3 może wpływać na progresję raka jajnika.
Wpływ platyny na ekspresję TMOD3 w raku jajnika
Chemioterapia jest istotnym czynnikiem wpływającym na progresję raka jajnika33. Aby zbadać, czy TMOD3 wpływa na progresję raka jajnika poprzez oddziaływanie na chemioterapię, w niniejszym badaniu najpierw sprawdzono, czy ekspresja TMOD3 ulega zmianie w raku jajnika po zastosowaniu leków chemioterapeutycznych. Zbiór danych GSE47856 wykazał, że poziom ekspresji TMOD3 spadł po leczeniu cisplatyną w linii komórkowej raka jajnika CH1 (Rysunek 2A). Zbiór danych GSE49577 wykazał, że ekspresja mRNA TMOD3 była niższa u myszy z ksenoprzeszczepami raka jajnika leczonych karboplatyną niż u myszy nieleczonych (Rysunek 2B). Ponadto zbiór danych GSE63885 wykazał również, że ekspresja TMOD3 była znacząco obniżona po chemioterapii opartej na platynie u pacjentek z rakiem jajnika wykazujących wysoką wrażliwość na platynę w porównaniu z pacjentkami z rakiem jajnika, które nie były leczone chemioterapią; nie stwierdzono jednak różnic u pacjentek z rakiem jajnika o umiarkowanej wrażliwości na platynę (Rysunek 2C). Powyższe wyniki sugerują, że TMOD3 może być potencjalnym celem działania platyny.
Ekspresja TMOD3 w raku jajnika wrażliwym lub opornym na platynę
Na podstawie powyższych wyników w badaniu postawiono hipotezę, że wysoka ekspresja TMOD3 w raku jajnika może wpływać na progresję tej choroby poprzez regulację oporności na platynę. W związku z tym w badaniu najpierw przeanalizowano ekspresję TMOD3 w liniach komórkowych opornych na leki lub w tkankach nowotworowych pacjentów, korzystając z wielu zbiorów danych GEO. Zarówno zbiory danych GSE15372, jak i GSE33482 wykazały, że ekspresja TMOD3 w linii komórkowej raka jajnika a2780 opornej na cisplatynę była znacznie wyższa niż w grupie wrażliwej na cisplatynę (Rysunek 3A,B). Zbiór danych GSE45553 wykazał, że ekspresja TMOD3 w sferoidach raka jajnika opornych na cisplatynę była znacznie wyższa niż w grupie wrażliwej na cisplatynę (Rysunek 3C). Wreszcie kompleksowa analiza pacjentów z rakiem jajnika leczonych chemioterapią z baz TCGA i GEO przy użyciu programu ROC plotter wykazała, że ekspresja TMOD3 była znacznie wyższa u pacjentek z rakiem jajnika z opornością na platynę niż u pacjentek wrażliwych na platynę (Rysunek 3D). Ponadto krzywa ROC wskazała TMOD3 jako marker pozwalający określić oporność na platynę w raku jajnika (Rysunek 3E). Powyższe wyniki sugerują, że wysoka ekspresja TMOD3 może odgrywać kluczową rolę w oporności na platynę w raku jajnika.
Wpływ ekspresji TMOD3 na przeżywalność pacjentek z rakiem jajnika
Analiza prognostyczna pacjentek z rakiem jajnika z bazy TCGA, przeprowadzona za pomocą narzędzia KM-plotter, wykazała, że wysoka ekspresja mRNA TMOD3 była istotnie skorelowana z niskim OS i PFS w całej grupie pacjentek oraz u pacjentek leczonych platyną (Rysunek 4). Co ciekawe, PFS u pacjentek leczonych platyną był bardziej znaczący niż w całej grupie pacjentek (Rysunek 4D). Ponieważ PFS jest lepszym wskaźnikiem wyników leczenia niż OS, analiza prognostyczna dodatkowo wykazuje krytyczną rolę TMOD3 u pacjentek z rakiem jajnika poddanych chemioterapii opartej na platynie.
Analiza mRNA-miRNA dla TMOD3 w raku jajnika
Regulacja potranskrypcyjna zapośredniczona przez mikroRNA jest istotnym mechanizmem wpływającym na ekspresję genów. Funkcję genu można przewidzieć na podstawie funkcji mikroRNA, które go celują34. Aby zbadać możliwy mechanizm podwyższonej ekspresji TMOD3 w raku jajnika, za pomocą predykcji TargetScan zidentyfikowano dwadzieścia jeden mikroRNA celujących w TMOD3. Ich korelacja z TMOD3 została również potwierdzona w próbkach raka jajnika z bazy TCGA (Rysunek 5A). Wykazano, że wszystkie one pośredniczą w chemorezystencji nowotworów, a większość (17/21) powiązano z opornością na cisplatynę w raku jajnika35,36,37,38,39,40,41,42,43,44,45,46,47,48, z wyjątkiem hsa-miR-301a-3p, hsa-miR-301b-3p, hsa-miR-365a-5p oraz hsa-miR-365b-5p. Ponadto analiza próbek raka jajnika z TCGA wykazała, że zarówno hsa-miR-454-3p, jak i hsa-miR-133a-3p najsilniej ujemnie korelują z ekspresją TMOD3 (Rysunek 5A-C). Co więcej, analiza Linkedomics wykazała, że hsa-miR-133a-3p wykazuje niską ekspresję w rakach jajnika opornych na cisplatynę (Rysunek 5E), podczas gdy w przypadku hsa-miR-454-3p nie stwierdzono istotności (Rysunek 5D). Zatem możliwe jest, że regulacja potranskrypcyjna zapośredniczona przez mikroRNA, a w szczególności przez hsa-miR-133a-3p, odpowiada za wysoką ekspresję TMOD3 w raku jajnika opornym na cisplatynę.
Związek między ekspresją TMOD3 a infiltracją immunologiczną w raku jajnika
Immunoterapia stanowi nowatorską i obiecującą strategię terapeutyczną w porównaniu z chemioterapią i chirurgią. Około 50% pacjentek z rakiem jajnika wykazuje spontaniczną infiltrację wieloma typami komórek odpornościowych49,50. Ponadto istotne jest zastosowanie immunoterapii u pacjentek z rakiem jajnika wykazującym chemooporność. Na przykład inhibitor liganda 1 programowanej śmierci komórki (PD-L1) Durvalumab, w połączeniu z inhibitorem PARP Olaparib, wykazał korzystne efekty terapeutyczne u pacjentek z rakiem jajnika z opornością na związki platyny51. Zbiór danych HPA wykazał, że TMOD3 jest szeroko rozpowszechniony w różnych komórkach odpornościowych (Rycina 6A,B), co sugeruje, że TMOD3 odgrywa istotną rolę w tych komórkach. W związku z tym w niniejszym badaniu za pomocą narzędzia TIMER zbadano korelację między TMOD3 a infiltracją immunologiczną w raku jajnika. TMOD3 był ujemnie skorelowany z czystością guza (Rycina 6C). Sugeruje to, że wiele innych komórek infiltruje tkankę raka jajnika i że TMOD3 może potencjalnie wpływać na raka jajnika poprzez te infiltrujące komórki. Rzeczywiście, dalsze wyniki wykazały, że ekspresja TMOD3 była dodatnio skorelowana z limfocytami B, komórkami T CD8+, neutrofilami oraz komórkami dendrytycznymi, ale nie z komórkami T CD4+ i makrofagami (Rycina 6D-I). W badaniu tym przeanalizowano ponadto związek między ekspresją TMOD3 a immunomodulatorami lub chemokinami w raku jajnika. Stwierdzono, że ekspresja TMOD3 była istotnie skorelowana z różnymi immunomodulatorami (Rycina 7A), w tym z uznanymi punktami kontrolnymi układu odpornościowego, takimi jak programowana śmierć 1 (PD1), PD-L1 (CD274) oraz cytotoksyczny antygen 4 limfocytów T (CTLA4) (Rycina 7B-D). Wyniki te sugerują, że TMOD3 może odgrywać istotną rolę w infiltracji immunologicznej w raku jajnika.
Sieci koekspresji TMOD3 w raku jajnika
W niniejszym badaniu zgłębiono potencjalne mechanizmy udziału TMOD3 w rozwoju raka jajnika. W pierwszej kolejności skonstruowano sieć koekspresji TMOD3 w raku jajnika. Na podstawie danych z LinkedOmics w próbkach raka jajnika z TCGA zidentyfikowano łącznie 4774 genów istotnie powiązanych z TMOD3. Mapa ciepła przedstawiła 50 najważniejszych genów (TOP50) wykazujących pozytywną lub negatywną korelację z TMOD3 (Rysunek 8A,B). Następnie, przy użyciu narzędzia STRING, wyselekcjonowano 295 genów oddziałujących z TMOD3. Po nałożeniu zbioru 4774 genów istotnie powiązanych z TMOD3 oraz 295 genów z nim oddziałujących, wyłoniono 76 genów (Rysunek 8C). Analiza wzbogacenia funkcjonalnego tych 76 genów wykazała, że najsilniej wzbogacone były procesy związane z cytoszkieletem aktynowym oraz adhezją ogniskową (Rysunek 8D). Analiza wzbogacenia szlaków KEGG wykazała, że najbardziej wzbogacona była ścieżka proteoglikanów w nowotworach (Rysunek 8E).
Terapeutyki małocząsteczkowe
Korelacje pomiędzy TMOD3 a potencjalnymi lekami lub małymi cząsteczkami przeanalizowano przy użyciu bazy danych CTD. Zidentyfikowano łącznie osiem leków wykazujących potencjalny efekt hamujący wobec TMOD3 (Rycina 9A), w tym paracetamol, kwas arystolochowy I, bufalinę, doksorubicynę, kwas foliowy, haloperidol, toksynę T-2 oraz tetrachlorodibenzo-p-dioksynę (Rycina 9B-I). Struktury cząsteczkowe zidentyfikowanych leków określono na podstawie bazy danych PubChem. Leki te stanowią wskazówki do przełamania procesu oporności na leki w raku jajnika pośredniczonego przez nadekspresję TMOD3.

Rycina 1: Poziomy ekspresji mRNA i białka TMOD3 w wielu zbiorach danych dotyczących raka jajnika. (A) Analiza ekspresji mRNA TMOD3 w prawidłowych tkankach jajnika (n = 20) i raku jajnika (n = 152) na podstawie zbioru danych GSE51088. (B) Analiza ekspresji mRNA TMOD3 w prawidłowych tkankach jajnika (n = 12) i raku jajnika (n = 57) na podstawie zbioru danych GSE66957. (C) Analiza ekspresji mRNA TMOD3 w prawidłowej tkance jajnika (n = 133) i raku jajnika (n = 374) za pomocą narzędzia TNM plotter. (D) Analiza ekspresji białka TMOD3 w prawidłowej tkance jajnika (n = 25) i raku jajnika (n = 100) w danych CPTAC za pomocą narzędzia UALCAN. (E) Analiza ekspresji mRNA TMOD3 w raku jajnika stopnia 2 i 3 (n = 94) oraz w raku jajnika stopnia 4 (n = 80) na podstawie zbioru danych GSE53963. (F) Analiza ekspresji mRNA TMOD3 w raku jajnika stopnia 1 i 2 (n = 10) oraz w raku jajnika stopnia 3 (n = 18) na podstawie zbioru danych GSE23554. Enter text here. (G) Analiza ekspresji mRNA TMOD3 w raku jajnika z przerzutami i bez przerzutów na podstawie zbioru danych GSE131978. Dane przedstawiono jako średnia ± SEM. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Poziomy ekspresji mRNA TMOD3 w liniach komórkowych raka jajnika, przeszczepionych guzach i u pacjentów po leczeniu lekami na bazie platyny. (A) Analiza ekspresji mRNA TMOD3 w liniach komórkowych raka jajnika CH1 leczonych cisplatiną (n = 3) i nieleczonych (n = 3) na podstawie zbioru danych GSE47856. (B) Analiza ekspresji mRNA TMOD3 w przeszczepach raka jajnika CH1 leczonych cisplatiną (n = 12) i nieleczonych (n = 19). (C) Analiza ekspresji mRNA TMOD3 u pacjentek z rakiem jajnika leczonych cisplatiną, z wysoką wrażliwością (HS, n = 13) i umiarkowaną wrażliwością na cisplatinę (n = 28), a także u pacjentek nieleczonych cisplatiną (n = 26). Dane przedstawiono jako średnia ± SEM. *w porównaniu z grupą nieleczoną, P < 0.05. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 3: Poziomy ekspresji mRNA TMOD3 w liniach komórkowych lub u pacjentów z rakiem jajnika wrażliwym i opornym na platynę. (A,B) Poziomy ekspresji mRNA TMOD3 w liniach komórkowych A2780 wrażliwych i opornych na cisplatynę przeanalizowano za pomocą zbioru danych GSE15372 (A) oraz zbioru danych GSE33482 (B). (C) Poziomy ekspresji mRNA TMOD3 w sferoidach raka jajnika wrażliwych i opornych na cisplatynę przeanalizowano za pomocą zbioru danych GSE45553. (D,E) Poziomy ekspresji mRNA TMOD3 u pacjentek z rakiem jajnika wrażliwym i opornym na platynę przeanalizowano za pomocą narzędzia ROC plotter (D), a TMOD3 zidentyfikowano jako marker określający oporność na platynę na podstawie krzywych ROC (E). Dane przedstawiono jako średnia ± SEM. Aby zobaczyć powiększoną wersję tego rysunku, kliknij tutaj.

Rysunek 4: Analiza korelacji poziomów ekspresji mRNA TMOD3 z OS i PFS u pacjentek z rakiem jajnika. (A,B) Analiza korelacji poziomów ekspresji mRNA TMOD3 z OS u wszystkich pacjentek z rakiem jajnika w bazie TCGA (A) lub u pacjentek otrzymujących chemioterapię opartą na platynie (B). (C,D) Analiza korelacji poziomów ekspresji mRNA TMOD3 z PFS u wszystkich pacjentek z rakiem jajnika w bazie TCGA (C) lub u pacjentek otrzymujących chemioterapię opartą na platynie (D). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5: Określanie zależności między TMOD3 a opornością na leki za pomocą analizy oddziaływań mRNA-microRNA. (A) Przewidywanie miRNA celujących w TMOD3. Kolor linii reprezentuje percentyl punktacji TargetScan Context++, a kolor węzła reprezentuje siłę korelacji między danym miRNA a TMOD3 w zbiorze danych TCGA dla raka jajnika. (B,C) Korelacja między poziomami ekspresji mRNA TMOD3 a hsa-miR454-3p (B) oraz hsa-miR133a-3p (C) w próbkach raka jajnika z bazy TCGA. (D,E) Ekspresja hsa-miR454-3p (D) oraz hsa-miR133a-3p (E) w wrażliwym i opornym na cisplatynę raku jajnika w próbkach z bazy TCGA. Dane przedstawiono jako średnia ± SEM. Aby zobaczyć powiększoną wersję tej ryciny, kliknij tutaj.

Rysunek 6: Korelacja poziomów ekspresji mRNA TMOD3 z infiltracją komórek odpornościowych w raku jajnika. (A,B) Poziomy ekspresji mRNA TMOD3 w komórkach odpornościowych w zbiorze danych HPA (A) oraz zbiorze danych Monaco (B). (C-I) Korelacja poziomów ekspresji mRNA TMOD3 z czystością guza w TCGA raka jajnika (C), limfocytami B (D), limfocytami T CD4+ (E), limfocytami T CD8+ (F), makrofagami (G), neutrofilami (H) i infiltracją komórek dendrytycznych (I). Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 7: Analiza korelacji poziomów ekspresji mRNA TMOD3 z czynnikami immunoregulatornymi u pacjentek z rakiem jajnika. (A) Mapa ciepła korelacji między ekspresją mRNA TMOD3 a immunoregulatorami w TCGA raka jajnika, obejmująca inhibitory odporności, immunostymulatory, MHC, chemokiny i receptory. (B-D) Analiza korelacji poziomów ekspresji mRNA TMOD3 z punktami kontrolnymi odporności, takimi jak PD1 (B), PD-L1 (CD274) (C) oraz CTLA4 (D) w próbkach TCGA raka jajnika. Prosimy kliknąć tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 8: Sieci koekspresji TMOD3 w raku jajnika. (A,B) Mapy ciepła przedstawiające 50 genów o najsilniejszej korelacji dodatniej (A) i ujemnej (B) z TMOD3 w zbiorach danych TCGA dla raka jajnika z wykorzystaniem LinkedOmics. Kolory czerwony i niebieski reprezentują odpowiednio geny dodatnio i ujemnie skorelowane. (C) Diagram Venna dla 295 genów oddziałujących z TMOD3 oraz 4774 genów skorelowanych z TMOD3 w TCGA dla raka jajnika. (D) Analiza wzbogacenia funkcjonalnego 76 wspólnych genów. Wykres bąbelkowy przedstawia 5 najważniejszych elementów istotnie wzbogaconych w BP, CC i MF. (E) TOP10 szlaków z analizy wzbogacenia szlaków KEGG. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 9: Przewidywanie potencjalnych leków lub małych cząsteczek celujących w TMOD3. (A) Przewidywanie potencjalnych leków celujących w TMOD3 z wykorzystaniem bazy danych CTD, gdzie kolor czerwony oznacza leki, które mogą hamować TMOD3. (B-H) Baza danych PubChem przedstawia struktury cząsteczkowe ośmiu wybranych leków. (B) Paracetamol, (C) kwas aristolochowy I, (D) bufalina, (E) doksorubicyna, (F) kwas foliowy, (G) haloperidol, (H) toksyna T-2. (I) Tetrachlorodibenzo-p-dioksyna. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Rycina uzupełniająca 1: Metody analizy poziomów ekspresji mRNA i białka TMOD3 w wielu zbiorach danych dotyczących raka jajnika. (A) Analiza ekspresji TMOD3 z wykorzystaniem bazy danych GEO. (B) Dalsza analiza danych wyekstrahowanych z bazy danych GEO przy użyciu oprogramowania statystycznego i graficznego. (C) Analiza ekspresji TMOD3 z wykorzystaniem TNMplot. (D) Analiza ekspresji białka z wykorzystaniem UALCAN. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 2: Analiza prognostyczna z wykorzystaniem KM-plotter. Prosimy kliknąć tutaj, aby pobrać ten plik.
Rysunek uzupełniający 3: Analiza krzywych ROC z wykorzystaniem programu ROC plotter. Aby pobrać ten plik, kliknij tutaj.
Rysunek uzupełniający 4: Metody analizy mRNA-miRNA. (A) Przewidywanie miRNA celujących w TMOD3 przy użyciu TargetScan. (B) Analiza korelacji TMOD3 z miRNA przy użyciu cBioportal. (C) Wizualizacja analizy korelacji przy użyciu Cytoscape. (D) Analiza ekspresji miRNA przy użyciu LinkedOmics. Prosimy kliknąć tutaj, aby pobrać ten plik.
Rysunek uzupełniający 5: Metody analizy immunoinfiltracji. (A) Analiza dystrybucji TMOD3 w różnych komórkach odpornościowych. (B) Analiza korelacji TMOD3 z infiltracją odpornościową przy użyciu cBioportal. (C) Analiza korelacji TMOD3 z immunomodulatorami przy użyciu TISIDB. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 6: Przewidywanie genów wchodzących w interakcje z TMOD3 przy użyciu STRING. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 7: Adnotacja funkcjonalna z użyciem DAVID. Prosimy kliknąć tutaj, aby pobrać ten plik.
Rycina uzupełniająca 8: Przewidywanie leków celujących w TMOD3 przy użyciu bazy danych CTD. Kliknij tutaj, aby pobrać ten plik.
Rycina uzupełniająca 9: Analiza struktury molekularnej leków z wykorzystaniem bazy danych PubChem. Proszę kliknąć tutaj, aby pobrać ten plik.