Name: Analiza czynników ekspresji genów nowotworowych za pomocą portalu internetowego CorExplorer (Wideo) | JoVE Journal
Uploaded: 2019-10-11T15:00:00.000+00:00
Duration: 8 min
Description: Odkryj analizę czynników ekspresji genów nowotworowych za pomocą portalu CorExplorer z algorytmem CorEx. Ta metoda pozwala badaczom identyfikować klastry w danych RNA-seq i oceniać powiązania z przeżyciem pacjentów.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

GV był wspierany przez DARPA award W911NF-16-0575.

1. Badanie czynników zawierających interesujący nas gen
<ol>
	<li>Otwórz przeglądarkę internetową i przejdź do http://corex.isi.edu, strony głównej CorExplorer.</li>
	<li>Po prawej stronie, w sekcji Szybkie linki, kliknij przycisk + rozwiń obok pozycji Jajnik (TCGA-OV), aby wyświetlić podsumowanie wykresu czynnika CorEx, który został wytrenowany na danych dotyczących raka jajnika TCGA (pokazanych w <silna klasa="xfig">Rysunek 1). Opcjonalnie kliknij inne, aby porównać.</li>
	<li>Po zakończeniu sprawdzania wykresów czynników kliknij Płuca (TCGA-LUAD), aby uzyskać dostęp do strony CorExplorer dla sekwencji RNA raka płuc.
	<ol>
		<li>Zapoznaj się z wykresem czynnika CorEx dla interesującego genu za pomocą okna "Wykres czynnika" CorExplorer.
		<ol>
			<li>Najedź kursorem myszy na okno wyświetlania wykresu czynnikowego. Powiększ wykres czynników za pomocą kółka przewijania myszy lub gładzika, aby zobaczyć szczegóły wykresu, takie jak najważniejsze geny w każdym czynniku i połączenia między węzłami na różnych warstwach. Możesz też kliknąć i przeciągnąć, aby przesunąć obszar obserwacji lub dowolny węzeł.</li>
			<li>Aby znaleźć gen docelowy (tutaj użyjemy BRCA1), kliknij menu rozwijane Gen w górnej części okna wykresu czynnikowego. Wpisz "BRCA1", aby wybrać go z listy rozwijanej i naciśnij Return, aby powiększyć widok do czynnika 26, współczynnika, z którym BRCA1 jest najsilniej skorelowany.</li>
			<li>Umieść kursor myszy na ekranie wykresu i przewiń, aby pomniejszyć, aby wyświetlić węzeł poziomu 2, L2_8 i skojarzone z nim czynniki, które są sąsiadami czynnika 26. Należy pamiętać, że wyświetlane są tylko geny o wadze większej niż próg wskazany na suwaku Min link weight.</li>
			<li>Aby zobaczyć wszystkie geny powiązane z tym czynnikiem, kliknij węzeł L1_26 i wybierz opcję Załaduj dodatkowe geny w wyskakującym okienku. Gdy pojawi się słowo "Gotowe", zamknij wyskakujące okienko.</li>
			<li>Teraz wróć do sekcji nagłówka nad oknem wykresu czynnikowego i chwyć i przeciągnij modyfikator Min. Teraz, gdy suwak wagi ogniwa zostanie przesunięty w dół do 0,05, inne geny czynnika L1_26, w tym BRCA2, pojawią się w kolejności wagowej. Opcjonalnie zmień położenie węzłów, chwytając i przeciągając, aby poprawić układ.</li>
		</ol>
		</li>
		<li>Określ, w jaki sposób stratyfikacja pacjentów w odniesieniu do czynnika wpływa na przeżycie, pytając w oknie przeżycia.
		<ol>
			<li>W oknie przeżycia odznacz opcję Sortuj według p-wartość, a następnie wybierz czynnik 26 z menu rozwijanego Pojedynczy czynnik, aby wyświetlić krzywe przeżycia dla czynnika 26.</li>
			<li>Przewiń w dół wykres przeżycia, aby pokazać liczbę zagrożonych pacjentów wzdłuż osi x.</li>
		</ol>
		</li>
		<li>Znajdź powiązania z funkcją biologiczną, wysyłając zapytanie w oknie Adnotacja.
		<ol>
			<li>W oknie adnotacji, aby posortować menu rozwijane Czynnik według numeru czynnika, a nie wskaźnika fałszywych wykryć (FDR), usuń zaznaczenie opcji Sortowanie FDR.</li>
			<li>Przewiń i kliknij, aby wybrać czynnik 26 z listy rozwijanej okna adnotacji, aby wyświetlić adnotacje wzbogacenia dla współczynnika.</li>
			<li>Przewiń listę adnotacji w dół, aż widoczna będzie naprawa DNA, a następnie kliknij ją, aby natychmiast zobaczyć powiązane geny podświetlone na żółto na ekranie wykresu. Zobacz środkowy panel Rysunek 2.</li>
			<li>Należy zauważyć, że czynniki znikają lub pojawiają się w miarę wybierania różnych terminów GO, w zależności od tego, czy są one wzbogacone o geny z wybraną adnotacją, np. "wewnętrzny apoptotyczny szlak sygnałowy w odpowiedzi na uszkodzenie DNA".</li>
		</ol>
		</li>
		<li>Zapoznaj się z tymi czynnikami, dodając okna o różnych funkcjach.
		<ol>
			<li>Na górnym pasku menu dodaj okno sieci interakcji białko-białko (PPI), wybierając PPI z listy rozwijanej Dodaj okno, a następnie kliknij przycisk Dodaj, aby dodać okno wykresu PPI do obszaru wyświetlania. W oknie wykresu PPI wybierz współczynnik "Warstwa1: 26", aby pokazać interakcje białko-białko. Zwróć uwagę na gęstość połączeń.</li>
			<li>Na górnym pasku menu, zamiast PPI, wybierz Mapa ciepła z listy rozwijanej Dodaj okno, a następnie kliknij przycisk Dodaj, aby dodać okno mapy ciepła do obszaru wyświetlania. W oknie mapy cieplnej wybierz współczynnik "Warstwa1: 26", aby wyświetlić wzorce ekspresji genów.</li>
			<li>Chwyć i zmień położenie okna mapy cieplnej tak, aby okno przetrwania było również widoczne. Wzdłuż górnej części mapy cieplnej obserwuj, jak pomarańczowy/niebieski/szary pasek odpowiada warstwom ryzyka pacjenta na wykresie przeżycia. Wyniki są pokazane na dole Rysunek 2.</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
2. Filtrowanie i interpretacja czynników CorEx na podstawie danych o wadze genów, przeżyciu i adnotacjach
<ol>
	<li>Filtrowanie pod kątem czynników zainteresowania przy użyciu przeżycia i jakości klastra.
	<ol>
		<li>Z menu rozwijanego Dataset (Zestaw danych) u góry wybierz TCGA_OVCA, aby przejść do strony CorExplorer dla sekwencji RNA raka jajnika TCGA.</li>
		<li>Po załadowaniu strony należy zauważyć w oknie przeżycia, że współczynnik o największej różnicy przeżycia dla różnych warstw wynosi 114.</li>
		<li>W górnej części okna wykresu czynników wybierz "Warstwa1: 114" z listy rozwijanej Czynnik.</li>
		<li>Chwyć suwak wagi linku za pomocą myszy i przesuń go w górę do 0.5. Należy zauważyć, że duża liczba genów w czynniku 114 (1609), z których żaden nie ma masy >0,35, wskazuje na stosunkowo słabe grupowanie.</li>
		<li>Następnie należy rozwinąć listę czynników w oknie przeżycia i wybrać następny najlepszy czynnik z listy rozwijanej okna przeżycia, czynnik 39, aby wyświetlić powiązane z nim krzywe przeżycia.</li>
		<li>Wybierz czynnik 39 w oknie adnotacji, klikając na niego. Pokazane są znaczące adnotacje GO i KEGG.</li>
	</ol>
	</li>
	<li>Aby lepiej zrozumieć biologiczną rolę genów w czynniku 39, zinterpretuj czynniki za pomocą informacji o adnotacjach sąsiedztwa w następujący sposób.
	<ol>
		<li>W górnej części okna wykresu czynnika wybierz czynnik "Warstwa1: 39" z listy rozwijanej współczynnika. Następnie najedź myszką na okno wykresu czynnikowego i pomniejsz je, aby odsłonić cały klaster L2_14 z 6 czynnikami: 14, 32, 39, 42, 52 i 82 (pokazanymi w <silna klasa="xfig">Rysunek 3).</li>
		<li>Aby zrozumieć względne znaczenie czynników związanych z węzłem L2_14, zacznij od wyświetlenia różnic w przeżywalności dla każdego z L2_14 czynników. Usuń zaznaczenie opcji Sortuj według p-val w oknie przetrwania, a następnie kliknij kolejno na każdą z liczb czynników. Robiąc to, zwróć uwagę, że tylko czynniki 14, 32 i 39 wykazują związek z przeżyciem.</li>
		<li>Teraz z górnego paska menu ponownie wybierz PPI z listy rozwijanej Dodaj okno. Naciśnij przycisk Dodaj, aby dodać okno wykresu PPI do obszaru wyświetlania. W oknie wykresu PPI wybierz współczynnik "Warstwa1: 52", aby pokazać interakcje białko-białko, które są istotne. Przykładowy układ okien w tym miejscu pokazany jest w Rysunek 3.</li>
		<li>Kliknij łącze View at StringDB (Wyświetl w bazie danych StringDB) u dołu okna PPI, aby połączyć się z internetową bazą danych StringDB. Kliknij przycisk Kontynuuj na pierwszym ekranie, a następnie wybierz kartę Analiza poniżej wykresu sieci, tak jak poprzednio, aby uzyskać analizę GO online dla genów sieci PPI. Najważniejszym składnikiem komórkowym jest "kompleks białkowy MHC klasy II".</li>
		<li>Wróć do zakładki CorExplorer i okna PPI i wybierz współczynnik 32, tym razem z listy rozwijanej współczynnika. Kliknij łącze Wyświetl w bazie danych StringDB do analizy bazy danych StringDB. Najważniejszym składnikiem komórkowym jest "kompleks białkowy MHC klasy I", w przeciwieństwie do klasy II dla czynnika 52 w poprzednim kroku!</li>
		<li>Na koniec wróć do okna PPI i wybierz "Warstwa 1: 39" z menu rozwijanego czynnika u góry. Kliknij łącze Wyświetl w bazie danych StringDB, aby utworzyć łącze do analizy bazy danych StringDB.</li>
		<li>Kliknij przycisk Kontynuuj na pierwszym ekranie, a następnie wybierz kartę Analiza poniżej wykresu sieci, aby uzyskać analizę GO online dla genów sieci PPI. Zauważ, że główną funkcją molekularną jest "wiązanie receptora chemokiny CXCR3".</li>
	</ol>
	</li>
</ol>
3. Korzystanie z adnotacji o przeżyciu i bazie danych w celu poszukiwania obiecujących kombinacji terapeutycznych
<ol>
	<li>Przełącz się na CorExplorer czerniaka TCGA, wybierając TCGA_SKCM z menu rozwijanego Dataset (Zestaw danych).</li>
	<li>Należy zauważyć, że czynnikiem o największej różnicy przeżycia jest czynnik 171. Zbadaj adnotacje czynnika 171, przewijając i zauważ, że "odpowiedź immunologiczna" i "szlak sygnałowy za pośrednictwem cytokin" znajdują się na górze (tak jak w przypadku najwyższego czynnika jajnika).</li>
	<li>Aby znaleźć czynnik komplementarny, należy zapoznać się z najważniejszymi czynnikami związanymi z przeżyciem wraz z ich najważniejszymi terminami adnotacji. Aby to zrobić, kliknij link Przegląd zbioru danych na górnym pasku menu, aby otworzyć osobną zakładkę zawierającą tabelę ze szczegółami przetwarzania zbioru danych, a także podsumowanie najważniejszych czynników według wartości p różnicy przeżycia. Należy pamiętać, że pierwszym czynnikiem nieimmunologicznym jest 88.</li>
	<li>Wróć do karty przeglądarki TCGA_SKCM.</li>
	<li>Wybierz współczynnik 88 w oknach przeżycia, adnotacji i wykresu. Kilka najpopularniejszych terminów GO jest związanych z "przetwarzaniem rRNA" i "organizacją mitochondrium", co potwierdza, że różni się od czynników związanych z odpornością.</li>
	<li>W oknie przeżycia, na liście rozwijanej sparowanych czynników, wybierz "88_171", aby zobaczyć, jak poprawia się przeżycie pacjentów z warstwy środkowej dla połączonych czynników ekspresji 171 i 88. Adnotacje i porównania przeżycia są zilustrowane na <silna klasa="xfig">Rysunek 4.</li>
</ol>
4. Znajdowanie podobieństw i różnic w zmienności ekspresji genów w różnych typach nowotworów za pomocą strony Szukaj
<ol>
	<li>Kliknij nagłówek CorExplorer, aby powrócić do strony głównej.</li>
	<li>Kliknij Szukaj na górnym pasku menu, aby przejść do strony umożliwiającej przeszukiwanie wszystkich zestawów danych w witrynie CorExplorer.</li>
	<li>W polu wyszukiwania genów wpisz "FLT1" (VEGFR1) i naciśnij Return lub naciśnij Szukaj. FLT1 występuje ze stosunkowo dużą masą w następujących czynnikach: OVCA - 76, LUAD - 162, SKCM - 195 i SKCM - 184, a także COAD - 112 i COAD - 74.</li>
	<li>Możesz też wyszukać powiązany termin języka GO we wszystkich zestawach danych. Spróbuj tego w polu "GO Search", wpisując "angiogeneza" i naciskając Return lub naciskając Szukaj. Wszystkie czynniki FLT1, z wyjątkiem SKCM-195, są wymienione jako statystycznie wzbogacone dla genów "angiogenezy" – czynnik 195 w rzeczywistości ma adnotację, ale poniżej domyślnego progu 10-8. Wyniki wyszukiwania dla tego i poprzedniego kroku są pokazane w Rysunek 5.</li>
	<li>Jako dalsze przykłady, w polu wyszukiwania GO najpierw wpisz "receptor naskórkowego czynnika wzrostu". Tylko LUAD jest wzbogacony o ten termin, dobrze znany czynnik stratyfikacji raka płuc. Następnie wpisz "mezenchymalny" w polu wyszukiwania. Termin ten jest wzbogacony w grupy ekspresji genów dla OVCA, gdzie jest dobrze zbadanym czynnikiem stratyfikacji.</li>
</ol>

<ol><li>Petryszak, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+RNASeq-er+API-a+gateway+to+systematically+updated+analysis+of+public+RNA-seq+data%5BTitle%5D)+AND+%22Bioinformatics%22%5BJournal%5D)">The RNASeq-er API-a gateway to systematically updated analysis of public RNA-seq data.</a> Bioinformatics. 33, 2218-2220 (2017).</li>
<li>Steeg, G. V., Galstyan, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Maximally+Informative+Hierarchical+Representations+of+High-Dimensional+Data%5BTitle%5D)+AND+%22Proceedings+of+the+Eighteenth+International+Conference+on+Artificial+Intelligence+and+Statistics+%28AISTATS%29%22%5BJournal%5D)">Maximally Informative Hierarchical Representations of High-Dimensional Data.</a> Proceedings of the Eighteenth International Conference on Artificial Intelligence and Statistics (AISTATS). , San Diego, CA. (2015).</li>
<li>Ver Steeg, G., Galstyan, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Discovering+structure+in+high-dimensional+data+through+correlation+explanation%5BTitle%5D)+AND+%22Advances+in+Neural+Information+Processing+Systems%22%5BJournal%5D)">Discovering structure in high-dimensional data through correlation explanation.</a> Advances in Neural Information Processing Systems. , Montreal, Canada. (2014).</li>
<li>Pepke, S., Ver Steeg, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+discovery+of+subsample+gene+expression+components+by+information+explanation%3A+therapeutic+implications+in+cancer%5BTitle%5D)+AND+%22BMC+medical+Genomics%22%5BJournal%5D)">Comprehensive discovery of subsample gene expression components by information explanation: therapeutic implications in cancer.</a> BMC medical Genomics. 10, 12(2017).</li>
<li>Byron, S. A., Van Keuren-Jensen, K. R., Engelthaler, D. M., Carpten, J. D., Craig, D. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Translating+RNA+sequencing+into+clinical+diagnostics%3A+opportunities+and+challenges%5BTitle%5D)+AND+%22Nature+Reviews+Genetics%22%5BJournal%5D)">Translating RNA sequencing into clinical diagnostics: opportunities and challenges.</a> Nature Reviews Genetics. 17, 257(2016).</li>
<li>Cancer Genome Atlas Research Network. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+molecular+profiling+of+lung+adenocarcinoma%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Comprehensive molecular profiling of lung adenocarcinoma.</a> Nature. 511, 543(2014).</li>
<li>Cancer Genome Atlas Network. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comprehensive+molecular+characterization+of+human+colon+and+rectal+cancer%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Comprehensive molecular characterization of human colon and rectal cancer.</a> Nature. 487, 330(2012).</li>
<li>Akbani, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genomic+classification+of+cutaneous+melanoma%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Genomic classification of cutaneous melanoma.</a> Cell. 161, 1681-1696 (2015).</li>
<li>Cancer Genome Atlas Research Network. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrated+genomic+analyses+of+ovarian+carcinoma%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Integrated genomic analyses of ovarian carcinoma.</a> Nature. 474, 609(2011).</li>
<li>Grossman, R. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Toward+a+shared+vision+for+cancer+genomic+data%5BTitle%5D)+AND+%22New+England+Journal+of+Medicine%22%5BJournal%5D)">Toward a shared vision for cancer genomic data.</a> New England Journal of Medicine. 375, 1109-1112 (2016).</li>
<li>Moynahan, M. E., Chiu, J. W., Koller, B. H., Jasin, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Brca1+controls+homology-directed+DNA+repair%5BTitle%5D)+AND+%22Molecular+Cell%22%5BJournal%5D)">Brca1 controls homology-directed DNA repair.</a> Molecular Cell. 4, 511-518 (1999).</li>
<li>Szklarczyk, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((STRING+v11%3A+protein%E2%80%93protein+association+networks+with+increased+coverage%2C+supporting+functional+discovery+in+genome-wide+experimental+datasets%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets.</a> Nucleic Acids Research. 47, 607-613 (2018).</li>
<li>Durgeau, A., Virk, Y., Corgnac, S., Mami-Chouaib, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Recent+advances+in+targeting+CD8+T-cell+immunity+for+more+effective+cancer+immunotherapy%5BTitle%5D)+AND+%22Frontiers+in+Immunology%22%5BJournal%5D)">Recent advances in targeting CD8 T-cell immunity for more effective cancer immunotherapy.</a> Frontiers in Immunology. 9, 14(2018).</li>
<li>Sato, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Intraepithelial+CD8%2B+tumor-infiltrating+lymphocytes+and+a+high+CD8%2B%2Fregulatory+T+cell+ratio+are+associated+with+favorable+prognosis+in+ovarian+cancer%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Intraepithelial CD8+ tumor-infiltrating lymphocytes and a high CD8+/regulatory T cell ratio are associated with favorable prognosis in ovarian cancer.</a> Proceedings of the National Academy of Sciences of the United States of America. 102, 18538-18543 (2005).</li>
<li>De Moura, M. B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mitochondrial+respiration-an+important+therapeutic+target+in+melanoma%5BTitle%5D)+AND+%22PLoS+One%22%5BJournal%5D)">Mitochondrial respiration-an important therapeutic target in melanoma.</a> PLoS One. 7, 40690(2012).</li>
<li>Folkman, J., Merler, E., Abernathy, C., Williams, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Isolation+of+a+tumor+factor+responsible+for+angiogenesis%5BTitle%5D)+AND+%22Journal+of+Experimental+Medicine%22%5BJournal%5D)">Isolation of a tumor factor responsible for angiogenesis.</a> Journal of Experimental Medicine. 133, 275-288 (1971).</li>
<li>Takahashi, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Vascular+endothelial+growth+factor+%28VEGF%29%2C+VEGF+receptors+and+their+inhibitors+for+antiangiogenic+tumor+therapy%5BTitle%5D)+AND+%22Biological+and+Pharmaceutical+Bulletin%22%5BJournal%5D)">Vascular endothelial growth factor (VEGF), VEGF receptors and their inhibitors for antiangiogenic tumor therapy.</a> Biological and Pharmaceutical Bulletin. 34, 1785-1788 (2011).</li>
<li>Subramanian, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Gene+set+enrichment+analysis%3A+a+knowledge-based+approach+for+interpreting+genome-wide+expression+profiles%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.</a> Proceedings of the National Academy of Sciences of the United States of America. 102, 15545-15550 (2005).</li>
<li>Cerami, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+cBio+Cancer+Genomics+Portal%3A+An+Open+Platform+for+Exploring+Multidimensional+Cancer+Genomics+Data%5BTitle%5D)+AND+%22Cancer+Discovery%22%5BJournal%5D)">The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data.</a> Cancer Discovery. 2, 401-404 (2012).</li>
<li>Gao, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrative+Analysis+of+Complex+Cancer+Genomics+and+Clinical+Profiles+Using+the+cBioPortal%5BTitle%5D)+AND+%22Science+Signalling%22%5BJournal%5D)">Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal.</a> Science Signalling. 6, 1(2013).</li>
<li>Reich, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((GenePattern+2.0%5BTitle%5D)+AND+%22Nature+Genetics%22%5BJournal%5D)">GenePattern 2.0.</a> Nature Genetics. 38, 500(2006).</li>
<li>Wang, Y. E., Kutnetsov, L., Partensky, A., Farid, J., Quackenbush, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((WebMeV%3A+A+Cloud+Platform+for+Analyzing+and+Visualizing+Cancer+Genomic+Data%5BTitle%5D)+AND+%22Cancer+Research%22%5BJournal%5D)">WebMeV: A Cloud Platform for Analyzing and Visualizing Cancer Genomic Data.</a> Cancer Research. 77, 11-14 (2017).</li>
<li>Morpheus. , Available from: 
 <a target="_blank" href="https://software.broadinstitute.org/morpheus">https://software.broadinstitute.org/morpheus</a> (2019).</li>
<li>Weitschek, E., Lauro, S. D., Cappelli, E., Bertolazzi, P., Felici, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((CamurWeb%3A+a+classification+software+and+a+large+knowledge+base+for+gene+expression+data+of+cancer%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">CamurWeb: a classification software and a large knowledge base for gene expression data of cancer.</a> BMC Bioinformatics. 19, 354(2018).</li>
<li>Chou, P. -H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((tACCo%2C+a+Database+Connecting+transcriptome+Alterations%2C+pathway+Alterations+and+Clinical+outcomes+in+Cancers%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">tACCo, a Database Connecting transcriptome Alterations, pathway Alterations and Clinical outcomes in Cancers.</a> Scientific Reports. 9, 3877(2019).</li>
</ol>

Autorzy oświadczają, że nie mają konkurencyjnych interesów finansowych.

Zaprezentowaliśmy stronę CorExplorer, publicznie dostępny serwer WWW do interaktywnej eksploracji maksymalnie skorelowanych czynników ekspresji genów poznanych z sekwencji RNA guza za pomocą algorytmu CorEx. Pokazaliśmy, w jaki sposób strona internetowa może być wykorzystywana do stratyfikacji pacjentów według ekspresji genów nowotworowych oraz jak taka stratyfikacja odpowiada funkcji biologicznej i przeżyciu.
Zbudowano inne serwery sieciowe do analizy sekwencyjnej RNA. Analiza różnicowa i koekspresja nowotworów może być badana i integrowana z innymi typami danych w cbioPortal 19,20. Serwery GenePattern21, Mev22 i Morpheus23 zawierają sprawdzone techniki grupowania, takie jak analiza głównych składowych (PCA), kmeans lub samoorganizujące się mapy (SOM). Bardziej innowacyjne działania obejmują CamurWeb24, oparty na automatycznym klasyfikatorze generującym reguły, oraz TACCO25, który implementuje losowe klasyfikatory lasu i lasso. Zastosowany tutaj algorytm CorEx optymalizuje informacje wielowymiarowe w celu znalezienia hierarchii czynników, które wyjaśniają wzorce w danych. Wydaje się, że nieliniowe i hierarchiczne uczenie się czynników zapewnia lepszą interpretowalność w stosunku do liniowych czynników globalnych znalezionych za pomocą PCA4. Dodatkowo, precyzyjna analiza sygnałów próbek w ramach tej techniki umożliwia precyzyjne porównania guzów z częściej stosowanymi szerokimi podtypami. To połączenie nakładającej się i hierarchicznej analizy czynnikowej odróżnia CorExplorer od większości innych podejść i wymaga nowych narzędzi do wizualizacji i podsumowania.
Kluczową częścią analizy czynnikowej CorExplorer jest możliwość zbadania nie tylko kilku, ale ponad 100 czynników z informacyjnymi wzorcami genów, które są umieszczone w nakładającej się hierarchii. CorExplorer ułatwia eksplorację tych niezliczonych czynników w celu uzyskania powiązań biologicznych i klinicznych oraz pozwala na wyjątkowo szczegółową charakterystykę poszczególnych nowotworów. Nienadzorowane uczenie się tak dużej liczby czynników oznacza, że nie wszystkie będą miały znaczenie dla biologii choroby. W takim przypadku konieczne jest użycie adnotacji lub znanych genów w celu wyodrębnienia interesujących czynników lub wyszukania czynników związanych z danymi klinicznymi, takimi jak przeżycie. W ten sposób CorExplorer pozwala użytkownikom zaimplementować ten bardzo ważny krok filtrowania. Obecność wzorców genów czynnikowych w guzie może nawet sugerować podejście do spersonalizowanego leczenia onkologicznego. Co więcej, mnogość wyników czynników dla każdego nowotworu, która pozwala na odkrycie potencjalnie użytecznych kombinacji terapeutycznych.
Czasami zdarza się, że nie pojawiają się żadne znaczące adnotacje GO dla czynników silnie skorelowanych z przeżyciem. Chociaż może się to zdarzyć z powodu zaszumionych lub niedostatecznie pobranych próbek, istnieją inne możliwe przyczyny, takie jak rozmiar klastra, który jest zbyt mały, aby zarejestrować znaczące wyniki wzbogacenia lub grupa będąca "koszykiem" pojedynczych genów z różnych szlaków bez spójnego związku biologicznego. Ponadto odpowiednia może być kategoria adnotacji inna niż w procesie biologicznym KEGG i GO, np. kompartment komórkowy. Dostęp do nich można uzyskać, łącząc się z bazą danych StringDB, jak pokazano w protokole. Analiza wzbogacania ontologii genów na stronie CorExplorer obecnie nie uwzględnia wagi genów w czynniku, chociaż prawdopodobnie zostanie to naprawione w najbliższej przyszłości. Zwróć uwagę, że w oknie "Dodaj dostępną jest opcja listy genów", która umożliwia pobranie pełnej listy genów czynnika w celu dalszej analizy za pomocą narzędzi zewnętrznych.
Na potrzeby strony internetowej CorEx został uruchomiony na każdym ze zbiorów danych pięć razy, a przebieg, który doprowadził do największej ogólnej korelacji, został zachowany. Posiadanie statystycznej reprezentacji wyników wielu serii może być bardziej pouczające i stanowi cel przyszłych prac. Ponadto zestaw typów nowotworów dostępnych na serwerze jest raczej niewielki, ale spodziewamy się, że z czasem będzie się rozszerzał w zależności od zainteresowania użytkowników.
Jak opisano powyżej, CorExplorer wizualizuje relacje między czynnikami sekwencyjnymi RNA CorEx wraz z informacjami klinicznymi i informacjami z bazy danych, umożliwiając w ten sposób wiele różnych trybów zapytania. Mamy nadzieję, że to narzędzie doprowadzi do dalszych prac nad wykorzystaniem mocy analizy sekwencyjnej RNA do odkryć i zastosowań klinicznych w onkologii.

Od czasu wprowadzenia nieco ponad dekadę temu, sekwencja RNA stała się wszechobecnym narzędziem do pomiaru ekspresji genów1. Dzieje się tak dlatego, że umożliwia szybkie i tanie profilowanie de novo całego transkryptomu próbki. Jednak dane dotyczące guza z sekwencją RNA odzwierciedlają podstawową biologię, która jest wewnętrznie złożona i często niedostatecznie próbkowana, podczas gdy same dane są wielowymiarowe i zaszumione. Stanowi to poważne wyzwanie dla ekstrakcji wiarygodnych sygnałów. Algorytm CorEx wykorzystuje wielowymiarowe informacje wzajemne, aby znaleźć subtelne wzorce w takich sytuacjach2,3 . Technika ta została wcześniej zaadaptowana do analizy próbek RNA guza jajnika z Atlasu Genomu Raka (TCGA) i w tym kontekście wydaje się, że ma znaczące zalety w porównaniu z częściej stosowanymi metodami analizy4.
Chociaż użycie sekwencjonowania RNA jest niezwykle rozpowszechnione w zastosowaniach badawczych, w tym w onkologii, wysiłki te nie doprowadziły do szerokiego zastosowania dla celów interwencji klinicznych5. Jednym z powodów takiego stanu rzeczy jest brak przyjaznych dla użytkownika algorytmów i oprogramowania ukierunkowanego na te konkretne problemy. Aby wypełnić tę lukę, zaprojektowaliśmy portal internetowy CorExplorer, aby umożliwić naukowcom z różnych środowisk badanie czynników ekspresji genów w próbkach sekwencyjnych RNA guza, znalezionych przez algorytm uczenia maszynowego CorEx. Portal CorExplorer obsługuje interaktywną wizualizację i odpytywanie czynników z kilku różnych typów nowotworów, w tym płuc, okrężnicy, czerniaka i jajnika6,7,8,9,10, z zamiarem pomocy badaczom w przesiewaniu korelacji danych i identyfikowaniu potencjalnych ścieżek stratyfikacji pacjentów do celów terapeutycznych.
Spodziewamy się, że portal CorExplorer może być przydatny dla kilku typów użytkowników. Portal został zaprojektowany z myślą o użytkownikach, którzy chcą zrozumieć szerokie czynniki wpływające na różnice w ekspresji genów nowotworowych w publicznych bazach danych, a być może także umieścić indywidualne profile ekspresji genów w kontekście nowotworów o podobnej charakterystyce. Oprócz przedstawionych tutaj reprezentatywnych protokołów, badania CorExplorer mogą służyć jako punkt wyjścia do sugerowania hipotez do dalszych testów, porównywania i kontrastowania wyników CorEx na zestawach danych poza CorExplorer oraz do łączenia patologicznych sygnatur ekspresji jednego lub kilku genów w pojedynczym guzie z większymi grupami, które mogą być skoordynowanie dotknięte. Wreszcie, może służyć jako przyjazne dla użytkownika wprowadzenie do zastosowania uczenia maszynowego do sekwencjonowania RNA dla osób rozpoczynających pracę w tej dziedzinie.

<table><tbody><tr><td>Serwer publiczny dla strony CorExplorer USC</td><td></td><td>http://corex.isi.edu</td><td>Intel Xeon E5-2690 4-rdzeniowy 2,6 GHz, 8GB RAM. Architektura backendowa to LAMP: Linux, Apache, MySQL, PHP.</td></tr><tr><td>Przeglądarka internetowa</td><td>: Google/Apple</td><td>Chrome/Safari</td><td>Zweryfikowane przeglądarki internetowe.</td></tr></tbody></table>

analiza czynników ekspresji genów nowotworowych za pomocą portalu internetowego corexplorer

Analyzing Tumor Gene Expression Factors with the CorExplorer Web Portal

Różnicowa analiza ekspresji genów jest ważną techniką dla zrozumienia stanów chorobowych. Algorytm uczenia maszynowego CorEx okazał się przydatny w analizie różnicowej ekspresji grup genów w sekwencyjnym RNA guza w sposób, który może być pomocny w rozwoju onkologii precyzyjnej. Jednak CorEx wytwarza wiele czynników, które mogą być trudne do przeanalizowania i połączenia z istniejącym zrozumieniem. Aby ułatwić takie połączenia, zbudowaliśmy stronę internetową CorExplorer, która umożliwia użytkownikom interaktywną eksplorację danych i odpowiadanie na często zadawane pytania związane z ich analizą. Przeszkoliliśmy CorEx na danych dotyczących ekspresji genów sekwencyjnych RNA dla czterech typów nowotworów: jajnika, płuc, czerniaka i jelita grubego. Następnie włączyliśmy do strony internetowej odpowiednie przeżycia, interakcje białko-białko, wzbogacenie szlaków Gene Ontology (GO) i Kyoto Encyclopedia of Genes and Genomes (KEGG) oraz mapy cieplne w celu powiązania z wizualizacją wykresu czynnikowego. W tym miejscu wykorzystujemy przykładowe protokoły, aby zilustrować wykorzystanie bazy danych do zrozumienia znaczenia poznanych czynników nowotworowych w kontekście tych danych zewnętrznych.

Wyszukiwanie genu „BRCA1” w zbiorze danych dotyczącym raka płuca wykazuje, że jest on najsilniej powiązany z czynnikiem CorEx 26 (Rysunek 2). Wzbogacenie terminów GO dla tego czynnika jest niezwykle wysokie, przy czym naprawa DNA wykazuje FDR na poziomie zaledwie 1 x 10-19. Wybór ten zwraca również uwagę na klaster drugiego poziomu L2_8, który posiada sześć blisko powiązanych czynników jako elementy potomne. Wybranie terminu „DNA repair” w adnotacjach GO lub w rozwijanym menu wzbogacenia GO na grafie czynników podświetla powiązane geny w każdym z czynników, przy czym czynnik 26 posiada ich zdecydowanie najwięcej, zgodnie z oczekiwaniami11. Sieć oddziaływań białko-białko jest silnie połączona, co dodatkowo potwierdza ściśle powiązaną funkcjonalność genów w czynniku 26. Powiązany wykres przeżycia sugeruje możliwą korelację z przeżywalnością pacjentów, jednak wymagałoby to potwierdzenia w większym zbiorze danych.
Rozpoczęcie analizy od przeżywalności pozwala na badanie przyczyn poprawy przeżywalności związanej z konkretnymi grupami ekspresji genów. Na przykład głównym czynnikiem wpływającym na przeżywalność w raku jajnika jest czynnik nr 39, który wykazuje silne wzbogacenie w genach związanych z układem odpornościowym (Rysunek 3). Pięć innych czynników powiązanych z tym samym węzłem poziomu 2 również jest wskazanych jako związane z odpornością, jednak wpływ na przeżywalność wydaje się być wśród nich bardzo zróżnicowany, przy czym czynnik 39 ma najwyższy, a czynnik 52 najniższy wpływ. Dodanie okna interakcji białko-białko dla danego czynnika pokazuje bezpośrednią sieć interakcji i umożliwia przejście do strony StringDB12 w celu sprawdzenia różnych wzbogaceń dla genów w sieci PPI. Przeprowadzając taką analizę dla każdego z czynników L2_14 po kolei, można stwierdzić, że wzbogacenia StringDB dla genów sieci PPI sugerują następujące możliwe wyjaśnienie powiązań z przeżywalnością. Czynnik 32 zawiera geny tworzące kompleks białkowy głównego układu zgodności tkankowej (MHC) klasy I, który jest rozpoznawany przez cytotoksyczne limfocyty T. Czynnik 39 odpowiada sygnalizacji cytokin i wiązaniu receptora CXCR3, co jest związane z limfocytami T CD8+. Oba te czynniki wydają się zapewniać znaczną przewagę w przeżywalności pacjentom wykazującym stosunkowo wysoką ekspresję odpowiadających im genów. Cytotoksyczne limfocyty T CD8+ są główno odpowiedzialne za odporność przeciwnowotworową. Z kolei czynnik 52 składa się z genów kodujących białka kompleksu MHC klasy II, które są rozpoznawane głównie przez pomocnicze limfocyty T CD4+, a nie bezpośrednio przez cytotoksyczne limfocyty T. Pozostałe czynniki L2_14 odzwierciedlają ogólną aktywację układu odpornościowego, która nie różnicuje tych dwóch populacji limfocytów. Powiązanie przeżywalności specyficzne dla rozpoznawania komórkowych antygenów MHC klasy I przez cytotoksyczne limfocyty T jest zgodne z naszą wiedzą o odporności przeciwnowotworowej w ujęciu ogólnym oraz z obserwacjami w innych nowotworach, takich jak czerniak13,14.
Portal internetowy wspiera odkrywanie par czynników o komplementarnych funkcjach, co może sugerować skuteczne, specyficzne dla guza terapie skojarzone. Przegląd zbioru danych pozwala na wyszukiwanie czynników, które korelują z przeżywalnością, a jednocześnie wykazują odmienne wzbogacenia GO. W przypadku czerniaka (TCGA_SKCM; Rycina 4) widać, że najważniejszy czynnik przeżywalności 171 jest związany z odpornością, podczas gdy czynnik 88, znajdujący się niżej na liście, wykazuje wzbogacenie w geny związane z organizacją mitochondriów. Rzeczywiście, zostało to zasugerowane jako cel terapeutyczny w czerniaku15. Dodanie okien przeżywalności na stronie CorExplorer umożliwia porównanie stratyfikacji z użyciem pary czynników do stratyfikacji dla każdego czynnika z osobna, co pokazuje, że korzystne wzorce ekspresji genów z obu grup wykazują trend przeżywalności lepszy niż w przypadku każdego z czynników osobno. Najwyższa warstwa nie wydaje się jednak ulec poprawie, co sugeruje, że sama immunoterapia może być najlepszą opcją dla niektórych pacjentów.
Podobieństwa i różnice między nowotworami można zaobserwować, przeszukując zestawy danych pod kątem genów lub terminów GO (Rysunek 5). Jako przykład, FLT1 (znany również jako VEGFR1) jest dobrze przebadanym markerem proangiogennym16,17. Po wpisaniu go w pasek wyszukiwania okazuje się, że wszystkie nowotwory posiadają czynniki, w których FLT1 odgrywa istotną rolę. Odwrotnie, po wprowadzeniu terminu GO „angiogenesis” na stronie wyszukiwania, 5 z 6 grup FLT1 pojawia się z tym wzbogaceniem. Wszystkie czynniki FLT1, z wyjątkiem SKCM-195, są wymienione jako statystycznie wzbogacone o geny związane z „angiogenesis”. Szósty czynnik faktycznie posiada tę adnotację, ale znajduje się on poniżej domyślnego progu 10-8. Gdy w alternatywnym kalkulatorze wzbogaceń, np. Gene Set Enrichment Analysis (GSEA)18, zostanie wykorzystane ważenie wewnątrz listy czynników, stwierdza się, że szósty czynnik jest również znacząco wzbogacony o geny „angiogenesis”.
Ważne jest sprawdzenie map ciepła, aby upewnić się, że wzorzec ekspresji genów jest odpowiedniej jakości, by wspierać interpretacje biologiczne. Mapy ciepła wykazujące wyraźną, silną zmienność mogą przedstawiać albo skoordynowaną ekspresję genów czynnikowych w zakresie od niskiej do wysokiej, albo bardziej złożone wzorce, w których niska ekspresja niektórych genów jest skorelowana z wysoką ekspresją innych (Rycina 6). Kluczowym wskaźnikiem wysokiej jakości grupy jest obecność kilku genów o płynnej zmienności ekspresji w funkcji wyniku czynnika. Mapy ciepła czynników przedstawiają próbki uporządkowane zgodnie z wynikiem czynnika, zatem powinien występować płynny gradient przesuwający się od lewej do prawej strony. Może to jednak nie nastąpić na co najmniej dwa różne sposoby. Najczęściej korelacje mogą być niezwykle zaszumione (Rycina 5C), co poddaje w wątpliwość wiarygodność i użyteczność wszelkich wniosków dotyczących przeżywalności i/lub funkcji biologicznej. Ponadto wzorce występujące tylko w niewielkiej mniejszości próbek mogą nie być zgodne z modelem trzech stanów ekspresji założonym przez algorytm CorEx, co prowadzi do błędnej klasyfikacji próbek (prawa strona Ryciny 5D).
<img alt="Interfejs strony internetowej analizy ekspresji genów w raku jajnika z diagramem uczenia maszynowego CorEx." class="xfigimg" src="/files/ftp_upload/60431/60431fig01.jpg" data-altupdatedat="1788282455000" data-alt="Rycina 1"> 
Rycina 1: Strona główna CorExplorer. Po kliknięciu znaku + obok Ovarian Cancer w sekcji Quick Links wyświetlane są szczegóły grafu czynnikowego. Hierarchiczny model CorEx składa się ze zmiennych wejściowych (w tym przypadku ekspresji genów) w dolnej warstwie oraz wywnioskowanych czynników utajonych w warstwach wyższych. <a href="https://www.jove.com/files/ftp_upload/60431/60431fig01large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>
<img alt="Teoria grafów i analiza sieci w danych biologicznych za pomocą CoR Explorer; interaktywne kroki procesu." class="xfigimg" src="/files/ftp_upload/60431/60431fig02.jpg" data-altupdatedat="1788282455000" data-alt="Rycina 2"> 
Rycina 2: Wykorzystanie nazwy genu do prowadzenia eksploracji. Rycina przedstawia serię zrzutów ekranu ilustrujących eksplorację czynników raka płuca w CorEx silnie powiązanych z BRCA1. Po pierwsze, wybranie „BRCA1” w rozwijanym menu Gene dla grafu czynników powoduje przybliżenie widoku grafu do czynnika, dla którego BRCA1 ma największą wagę. Lekkie oddalenie kadru ukazuje węzeł drugiej warstwy L2_8 łączący ten czynnik z innymi powiązanymi czynnikami. Można porównać przeżywalność i adnotacje: kliknięcie terminu GO DNA repair wyróżnia adnotowane geny. Dodano okno PPI, aby pokazać interakcje sieciowe dla genów w danym czynniku. Użycie przycisku Add Window w celu dodania mapy ciepła pokazuje powiązanie wzorców ekspresji z przeżywalnością, co sugeruje, że zwiększona ekspresja genów naprawy DNA może być powiązana ze zmniejszoną przeżywalnością. <a href="https://www.jove.com/files/ftp_upload/60431/60431fig02large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.</a>
<img alt="Sieć oddziaływań białkowych z analizą przeżywalności, schemat przedstawiający relacje genowe i krzywe przeżywalności." class="xfigimg" src="/files/ftp_upload/60431/60431fig03.jpg" data-altupdatedat="1788282455000" data-alt="Rysunek 3"> 
Rysunek 3: Wykorzystanie danych klinicznych (przeżywalności) do ukierunkowanej eksploracji. Analiza głównego czynnika powiązanego z przeżywalnością (39) w raku jajnika ujawnia interesujące relacje między sąsiednimi czynnikami. Po wybraniu czynnika 39 na grafie czynników i niewielkim oddaleniu widać, że czynnik drugiej warstwy powiązany z czynnikiem 39 ma pięć innych powiązanych czynników. Dodatkowe okno przeżywalności umożliwia bezpośrednie porównanie powiązanych różnic w przeżywalności. Czynniki 39 i 32 wykazują pozytywną korelację z przeżywalnością, w przeciwieństwie do czynnika 52, który jej nie wykazuje. Sieci oddziaływań białko-białko są dobrze zdefiniowane. Odnośniki do StringDB pozwalają na porównanie adnotacji GO (nie pokazano): czynnik 39 jest powiązany z siecią sygnalizacji cytokin związanych z aktywacją cytotoksycznych limfocytów T CD8+, a w czynniku 32 dominują białka prezentujące antygen MHC klasy I, które wyzwalają rozpoznanie przez takie limfocyty; sąsiednie czynniki są jednak zdominowane przez inne komponenty układu odpornościowego, takie jak pomocnicze limfocyty T CD4+, i nie wykazują korelacji z przeżywalnością. <a href="https://www.jove.com/files/ftp_upload/60431/60431fig03large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.</a>
<img alt="Wykres analizy ekspresji genów z krzywymi przeżycia, adnotacjami i wizualizacją danych." class="xfigimg" src="/files/ftp_upload/60431/60431fig04v3.jpg" data-altupdatedat="1788282455000" data-alt="Rycina 4"> 
Rycina 4: Analiza głównych czynników przeżycia sugeruje potencjalne kombinacje terapeutyczne. Link „Datasets” na pasku menu strony głównej prowadzi do zwięzłej tabeli czynników przeżycia uporządkowanych według wartości p, wraz z główną adnotacją GO (niewidoczną na rysunku). Wykorzystując te informacje dla czerniaka, kombinacja czynnika 171 dla funkcji odpornościowej z czynnikiem 88 dla organizacji mitochondriów wydaje się komplementarna. Na rysunku przedstawiono okna adnotacji dla każdego z czynników obok siebie w celu ich zestawienia. Krzywe przeżycia dla pacjentów stratyfikowanych według dwóch czynników osobno lub razem wskazują, że kombinacja ta zwiększa różnicę w przeżywalności w porównaniu do każdego z czynników stosowanego pojedynczo. <a href="https://www.jove.com/files/ftp_upload/60431/60431fig04large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>
<img alt="Wykres wyników wyszukiwania genów i GO; zbiór danych, czynnik, waga, MI dla analizy badawczej." class="xfigimg" src="/files/ftp_upload/60431/60431fig05.jpg" data-altupdatedat="1788282455000" data-alt="Rysunek 5"> 
Rysunek 5: Strona wyszukiwania umożliwia analizę panrakową. Geny lub terminy procesów biologicznych GO można wyszukiwać we wszystkich zbiorach danych, korzystając z łącza Search na stronie głównej. Rysunek przedstawia wyniki wyszukiwania dla genu FLT1 i terminu GO „angiogenesis”. Wyniki wskazują obecność FLT1 w czynnikach opisanych terminem „angiogenesis” w różnych nowotworach. <a href="https://www.jove.com/files/ftp_upload/60431/60431fig05large.jpg" target="_blank">Prosimy kliknąć tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Analiza map ciepła ekspresji genów za pomocą CoRExplorer; wizualna reprezentacja zbioru danych TCGA-LUAD." class="xfigimg" src="/files/ftp_upload/60431/60431fig06v2.jpg" data-altupdatedat="1788282455000" data-alt="Rysunek 6"> 
Rysunek 6: Mapy ciepła mogą być wykorzystywane do jakościowej oceny korelacji między genami a próbkami w zależności od wyniku czynnika. Relacje wysokiej jakości w ekspresji genów są widoczne jako płynne przejścia, gdy pacjenci w mapach ciepła są uporządkowani według wyniku czynnika. Przykładem jest najbardziej lewa mapa ciepła dla czynnika 18. Wzorce mogą obejmować również złożone sygnatury nadekspresji i niedoekspresji, jak w środkowej dużej mapie ciepła dla czynnika 11. Wzorce niższej jakości czasami wykazują gwałtowne zmiany ekspresji dla podgrupy pacjentów, jak w mapie ciepła dla czynnika 9 po prawej stronie, lub proste, bardzo szumiące korelacje, jak w mapie ciepła dla czynnika 161 w prawym dolnym rogu. <a href="https://www.jove.com/files/ftp_upload/60431/60431fig06large.jpg" target="_blank">Aby zobaczyć powiększoną wersję tego rysunku, kliknij tutaj.</a>

Przedstawiamy portal internetowy CorExplorer, zasób do eksploracji czynników sekwencjonowania RNA nowotworów znalezionych przez algorytm uczenia maszynowego CorEx (Correlation Explanation) i pokazujemy, jak można analizować czynniki związane z przeżyciem, adnotacjami w bazie danych, interakcjami białko-białko i sobą nawzajem, aby uzyskać wgląd w biologię nowotworu i interwencje terapeutyczne.

Analiza czynników ekspresji genów nowotworowych za pomocą portalu internetowego CorExplorer

Differential gene expression analysis is an important technique for understanding disease states. The machine learning algorithm CorEx has shown utility ...

CorExplorer znajduje czynniki ekspresji genów związane z nowotworem w sposób, który jest matematycznie zasadniczy, za pomocą interaktywnych wizualizacji, a także informacji klinicznych i baz danych do odkryć biologicznych. CorEx maksymalizuje informacje w sposób, który wymaga stosunkowo niewielu próbek do znalezienia klastrów w danych wielowymiarowych. Hierarchia czynników utajonych, które wytwarza, ułatwia zrozumienie biologiczne.Metoda ta okazała się już interesująca, gdy zostanie zastosowana w biologii raka. Jego podstawy teorii informacji sprawiają, że można go zastosować do każdego systemu z wieloma zmiennymi i niewielką wcześniejszą wiedzą na temat ich relacji. Zacznij od przejścia do strony głównej CorExplorer.W sekcji Szybkie linki kliknij przycisk rozwijania plus, aby wyświetlić podsumowanie wykresu czynnika CorEx, który został wytrenowany na interesujących danych dotyczących raka. Na przykład w tym miejscu pokazane są wykresy czynników dla danych dotyczących raka jajnika TCGA. Po sprawdzeniu wykresów czynników kliknij płuca TCGA-LUAD, aby uzyskać dostęp do strony CorExplorer dotyczącej sekwencjonowania RNA raka płuc i użyj okna Wykres czynnika CorExplorer, aby zapoznać się z wykresem czynnika CorEx dla interesującego nas genu.Przesuwając kursor myszy nad oknem wyświetlania wykresu czynników, powiększ wykres czynników za pomocą gładzika myszy, aby wyświetlić szczegóły wykresu, takie jak najważniejsze geny w każdym czynniku i połączenia między węzłami na różnych warstwach. Aby zlokalizować gen docelowy, kliknij menu Gen i wpisz nazwę genu, aby wybrać go z listy rozwijanej. Naciśnij Return na klawiaturze, aby powiększyć widok do czynnika, z którym gen zainteresowania jest najsilniej skorelowany.Umieść kursor myszy na ekranie wykresu i przewiń, aby pomniejszyć, aby zobaczyć poziom węzła i powiązanych z nim czynników, które sąsiadują z najściślej powiązanym czynnikiem genowym. Należy pamiętać, że wyświetlane są tylko geny o wadze większej niż próg wskazany na suwaku minimalna waga linku. Aby wyświetlić wszystkie geny powiązane z czynnikiem, kliknij odpowiedni węzeł i wybierz opcję Załaduj dodatkowe geny w wyskakującym okienku.Gdy pojawi się Gotowe, zamknij wyskakujące okienko. W sekcji nagłówka kliknij i przeciągnij modyfikator min Link Weight do 0.05, aby umożliwić pojawienie się genów w kolejności wagowej. Aby zidentyfikować powiązania z funkcją biologiczną, w oknie Adnotacja usuń zaznaczenie opcji Sortuj wskaźnik fałszywych odkryć, aby posortować menu rozwijane Czynnik według liczby czynnika, a nie współczynnika fałszywych odkryć.Przewiń i kliknij, aby wybrać interesujący Cię czynnik w menu rozwijanym okna Adnotacja, aby wyświetlić adnotacje wzbogacenia dla tego czynnika. Następnie kliknij czynnik wzbogacenia, aby natychmiast wyświetlić powiązane geny, podświetlone na żółto na wykresie. Należy zauważyć, że czynniki, które znikają lub pojawiają się jako różne terminy GO, są wybierane w zależności od tego, czy są wzbogacone o geny z wybraną adnotacją, czy nie.W przypadku filtrowania czynników zainteresowania, przy użyciu przeżycia i jakości klastra, z menu rozwijanego Zestaw danych wybierz TCGA_OVCA, aby przejść do strony CorExplorer sekwencjonowania RNA raka jajnika TCGA. Zanotować w oknie Przeżycie czynnik o największej różnicy przeżycia i wybrać ten czynnik w oknie Wykres czynnika z menu rozwijanego Czynnik. Kliknij i przeciągnij suwak Link Weight (Waga łącza) do pozycji 0,5 i zanotuj liczbę genów we współczynniku.Rozwiń listę czynników w oknie Przetrwanie i kliknij następny najlepszy czynnik w menu rozwijanym Okno przetrwania, aby wyświetlić powiązane z nim krzywe przeżycia. Zostaną wyświetlone znaczące adnotacje GO i Kegg. Aby lepiej zrozumieć biologiczną rolę genów w tym czynniku, wybierz warstwę czynnika w górnej części okna Wykres czynnika i przesuń kursor myszy nad okno, jednocześnie pomniejszając, aby wyświetlić cały klaster i powiązane czynniki.Aby zrozumieć względne znaczenie czynników powiązanych z węzłem klastra, odznacz opcję Sortuj według p-wartości w oknie Przeżycie i kliknij na każdy z numerów czynników po kolei, aby je wyświetlić, zwracając uwagę na czynniki, które wyświetlają powiązanie przeżycia. W menu Dodaj okno wybierz opcję PPI i kliknij przycisk Dodaj, aby dodać okno wykresu PPI do obszaru wyświetlania. W oknie wykresu PPI wybierz warstwę czynnika, która Cię interesuje, aby wyświetlić interakcje białko-białko, które są istotne.Kliknij łącze Wyświetl w bazie danych StringDB, aby połączyć się z internetową bazą danych STRINGdb, a następnie kliknij przycisk Kontynuuj. Następnie otwórz zakładkę Anaylsis, aby uzyskać analizę GO online dla genów sieci PPI. Zostanie wyświetlony górny komponent komórkowy.Wróć do zakładki CorExplorer i okna PPI i wybierz inny czynnik. Kliknij ponownie łącze Wyświetl w bazie danych StringDB. Zostanie wyświetlony inny górny komponent komórkowy.Następnie w oknie PPI wybierz inny czynnik do analizy bazy danych STRING. Aby znaleźć podobieństwa i różnice w zmienności ekspresji genów w różnych typach nowotworów, kliknij nagłówek CorExplorer, aby powrócić do strony głównej, a następnie kliknij przycisk Szukaj, aby uzyskać dostęp do strony umożliwiającej przeszukiwanie wszystkich zestawów danych w witrynie CorExplorer. W polu Wyszukiwanie genów wprowadź nazwę genu, który Cię interesuje, a następnie kliknij przycisk Wyszukaj.Na przykład, jak wykazano, FLT1 występuje ze stosunkowo dużą masą i wieloma różnymi czynnikami. Poszukiwanie genu BRCA1 w zbiorze danych dotyczących raka płuc ujawnia, że gen ten jest najsilniej związany z czynnikiem CorEx 26. Wzbogacenie terminu GO dla tego czynnika jest niezwykle wysokie, a naprawa DNA wykazuje wskaźnik fałszywych odkryć tylko jeden razy 10 do ujemnych 19.Wybór zwraca również uwagę na L2_8 klastra drugiego poziomu, który w dzieciństwie ma sześć ściśle powiązanych ze sobą czynników. Sieć interakcji białko-białko naprawy DNA jest silnie powiązana, co dodatkowo wspiera ściśle powiązaną funkcjonalność genów czynnika 26. Powiązane wykresy przeżycia sugerują możliwy związek z przeżyciem pacjenta, który musiałby zostać potwierdzony w większym zbiorze danych.Rozpoczęcie od oceny przeżycia może pozwolić na analizę czynników, które korelują z poprawą przeżycia w odniesieniu do określonych grup ekspresji genów. Dodanie okna interakcji białko-białko dla każdego czynnika ułatwia z kolei określenie możliwych wyjaśnień ich powiązań z przeżyciem. Ważne jest, aby sprawdzić mapy cieplne dla każdego czynnika, aby potwierdzić, że wzorzec ekspresji genów jest odpowiedniej jakości, aby wspierać interpretacje biologiczne.Mapy cieplne, które pokazują silne, wyraźne różnice we wzorcach ekspresji genów, mogą wykazywać albo skoordynowaną ekspresję genów czynnika, od wysokiej do niskiej, albo bardziej złożone wzorce, przy czym niektóre geny mają niską ekspresję skorelowaną z innymi genami o wysokiej ekspresji. Nadrzędnym celem tej procedury jest ustalenie spersonalizowanej terapii poprzez mapowanie guza poza próbką do czynników CorExplorer w celu zidentyfikowania potencjalnych terapii specyficznych dla nowotworu.

analyzing-tumor-gene-expression-factors-with-corexplorer-web

Research

JoVE Journal

Cancer Research

Badania nad rakiem

Degradacja kwasów nukleinowych w tkance archiwalnej, niejednorodność guzów i brak świeżo zamrożonych próbek tkanek mogą negatywnie wpłynąć na usługi diagnostyczne raka w laboratoriach patologicznych na całym świecie. Niniejszy manuskrypt opisuje optymalizację panelu biomarkerów przy użyciu multipleksowego testu kulek magnetycznych do klasyfikacji raka piersi.

Optymalizacja testu ekspresji opartego na multipleksowym RNA z wykorzystaniem materiału archiwalnego raka piersi

Tutaj prezentujemy protokół do badania biomarkera i predyktora przeżycia raka piersi na podstawie kompleksowej analizy zbiorczych zestawów danych klinicznych pochodzących z różnych publicznie dostępnych baz danych, wykorzystując strategię ekspresji, korelacji i analizy przeżycia krok po kroku.

Eksploracja danych i analiza integracyjna biomarkera raka piersi przy użyciu wielu publicznie dostępnych baz danych

Bioinformatyka to użyteczny sposób na przetwarzanie dużych zbiorów danych. Dzięki wdrożeniu podejść bioinformatycznych naukowcy mogą szybko, niezawodnie i skutecznie uzyskiwać wnikliwe aplikacje i odkrycia naukowe. W artykule przedstawiono wykorzystanie bioinformatyki w badaniach nad rakiem jajnika. Z powodzeniem weryfikuje również wyniki badań bioinformatycznych poprzez eksperymenty.

Integracja podejść bioinformatycznych i walidacji eksperymentalnych w celu zrozumienia roli sygnalizacji Notch w raku jajnika

Prezentowany tutaj jest protokół do odkrycia nadekspresji genów napędzających, które utrzymują ustalone rakopodobne komórki macierzyste pochodzące z komórek HT29 jelita grubego. Przeprowadzono sekwencję RNAseq z dostępnymi danymi bioinformatycznymi w celu zbadania i zbadania sieci ekspresji genów w celu wyjaśnienia potencjalnego mechanizmu związanego z przetrwaniem docelowych komórek nowotworowych.

Odkrycie genów sterujących w guzach przypominających łodygi raka jelita grubego pochodzących z HT29

Tropomodulina 3 (TMOD3) jest coraz częściej badana w nowotworach w ostatnich latach. To badanie jest pierwszym, które donosi, że TMOD3 ulega silnej ekspresji w raku jajnika i jest ściśle związany z opornością na platynę i infiltracją immunologiczną. Wyniki te mogą pomóc w poprawie wyników leczenia raka jajnika.

Nadekspresja tropomoduliny 3 jako marker oporności na platynę i infiltracji immunologicznej w raku jajnika

Protokół ten przedstawia metody hodowli, manipulowania i analizowania mysich CAF, w tym testy funkcjonalne in vitro/in vivo , transkryptomikę i analizę obliczeniową. Ponadto opisano w nim platformę internetową opracowaną w celu analizy publicznie dostępnych transkryptomów mikrowypreparowanych składników zrębu i nabłonka z kolekcji ludzkich nowotworów raka piersi.

Fibroblasty związane z rakiem z guzów sutka myszy jako narzędzia do badań molekularnych i obliczeniowych

Ten protokół pokazuje, jak generować fluorescencyjnie oznaczone, genetycznie zdefiniowane guzy klonalne w oku/czułkach Drosophila imaginal discs (EAD). Opisuje, jak rozdzielić EAD i mózg od larw trzeciego stadium rozwojowego i jak je przetworzyć w celu wizualizacji i ilościowego określenia zmian ekspresji genów i inwazyjności guza.

Model guza dysku imaginalnego Drosophila: wizualizacja i kwantyfikacja ekspresji genów i inwazyjności guza za pomocą mozaik genetycznych

Prezentujemy protokół demonstrujący nowatorski somatyczny system transferu genów wykorzystujący RIP-Tag; Mysi model RIP-tva do badania funkcji genów w przerzutach. Ptasie retrowirusy są dostarczane dosercowo, aby zapewnić transfer genów do przedzłośliwych, nieinwazyjnych zmian komórek β trzustki u dorosłych myszy.

Identyfikacja i charakterystyka czynników przerzutowych poprzez transfer genów do nowego RIP-Tag; Model mysi RIP-tva

Ten protokół służy jako schemat do tworzenia funkcjonalnego systemu Tet-ON w liniach komórek nowotworowych i jego późniejszego wykorzystania, w szczególności do badania roli białek pochodzących z komórek nowotworowych w rekrutacji monocytów/makrofagów do mikrośrodowiska guza.

Warunkowe obniżenie ekspresji genów w liniach komórek nowotworowych w celu zbadania rekrutacji monocytów/makrofagów do mikrośrodowiska guza

Opisujemy protokół przedklinicznego śledzenia in vivo przerzutów nowotworowych. Opiera się on na reporterze radionuklidowo-fluorescencyjnym łączącym symporter jodku sodu, wykrywany przez nieinwazyjny [18F]tetrafluoroboran-PET, oraz białko fluorescencyjne w celu usprawnionego potwierdzenia ex vivo. Metoda ma zastosowanie do przedklinicznego śledzenia komórek in vivo poza biologią guza.

Obrazowanie genu reporterowego radionuklidowo-fluorescencyjnego w celu śledzenia progresji guza w modelach guzów gryzoni

Odkryj analizę czynników ekspresji genów nowotworowych za pomocą portalu CorExplorer z algorytmem CorEx. Ta metoda pozwala badaczom identyfikować klastry w danych RNA-seq i oceniać powiązania z przeżyciem pacjentów.

Analiza czynników ekspresji genów nowotworowych za pomocą portalu internetowego CorExplorer

Przeglądaj więcej filmów