$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Dane spektralne MS/MS uzyskano z pozbawionych elementów identyfikacyjnych próbek PTF, które zostały zebrane przy użyciu procedur zgodnych z wytycznymi i przepisami zatwierdzonymi przez radę instytucji, jak wcześniej opisano21,29,30.
UWAGA: Rysunek 1 przedstawia przegląd całego przepływu pracy, który składa się z pięciu modułów. Wszystkie dane wejściowe, wyjściowe i narzędzia programowe podsumowano w tabeli uzupełniającej 1.

Rysunek 1: Podsumowanie modułów przepływu pracy metaproteomiki klinicznej w Galaxy. Kompletny przepływ pracy metaproteomiki klinicznej składa się z pięciu modułów: Generowanie bazy danych, Odkrywanie, Weryfikacja, Kwantyfikacja i Interpretacja danych. (A) Obszerna, obszerna baza danych zawiera sekwencje białek z gatunków drobnoustrojów, które uważa się za obecne w próbce, ludzi i powszechnych zanieczyszczeń. Narzędzie programowe MetaNovo bezpośrednio dopasowało dane spektralne MS/MS do peptydów i wywnioskowało białka i ich organizm źródłowy na podstawie surowych danych MS i dużej bazy danych sekwencji wejściowych białek, tworząc zredukowaną bazę danych33. Zredukowana baza danych z MetaNovo jest następnie łączona z białkami ludzkimi i białkami zanieczyszczeń w celu stworzenia bazy danych do odkrywania peptydów. orazDwa algorytmy identyfikacji peptydów, SearchGUI/PeptideShaker i MaxQuant, dopasowują sekwencje peptydów do widm MS/MS i bazy danych białek wabiących cel49. orazPeptydy zidentyfikowane przez SearchGUI/PeptideShaker i MaxQuant są następnie weryfikowane za pomocą PepQuery2. PepQuery2 rygorystycznie ponownie bada przypuszczalnie zidentyfikowane sekwencje peptydów mikrobiologicznych i ich dopasowane widma MS/MS w porównaniu z innymi potencjalnymi dopasowaniami do ludzkiego proteomu gospodarza i/lub zanieczyszczeń, weryfikując w ten sposób wysoce pewne dopasowania mikrobiologiczne40,41. Zweryfikowane peptydy są wykorzystywane do generowania zweryfikowanej bazy danych sekwencji białek, która zostanie wykorzystana do ilościowego oznaczania peptydów i białek. (D) MaxQuant42 przeszukuje dane MS/MS w stosunku do zweryfikowanej sekwencji białek i określa ilościowo peptydy mikrobiologiczne i wywnioskowane białka wraz z białkami ludzkimi. (E) Unipept45 i MSstatsTMT46 są używane w ostatnim kroku do opisywania białek za pomocą taksonomii i informacji funkcjonalnych (akcesje Komisji ds. Enzymów), a także do generowania wykresów wulkanów i porównań. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
1. Etykietowanie TMT i generowanie widm MS/MS
- Aby przygotować się do analizy stwardnienia rozsianego, należy przeprowadzić pobieranie próbek klinicznych zgodnie z wytycznymi i przepisami.
UWAGA: Ponieważ niniejszy protokół kładzie nacisk na bioinformatyczny przepływ pracy, procedury pobierania próbek klinicznych mogą różnić się od tych, które zostały użyte w tym manuskrypcie. W tym przypadku białka zostały tryptycznie strawione do mieszaniny peptydów, znakowane, frakcjonowane i analizowane za pomocą spektrometrii mas w celu wygenerowania danych spektralnych MS/MS do dalszej analizy przy użyciu platformy Galaxy. Szczegółowe instrukcje przetwarzania próbek zostały wcześniej opisane przez Boylan et al.29 oraz Afiuni-Zadel et al.30.
- Wyizoluj białka z próbek klinicznych i strawij je na peptydy za pomocą trypsyny29,30.
- Oznaczyć białka odczynnikiem Tandem Mass Tag (TMT)-11-plex. Ten odczynnik do znakowania pomoże w ilościowym określeniu peptydów i białek31,32.
- Podziel oznaczone próbki losowo i równomiernie na cztery grupy eksperymentalne oparte na TMT.
- Dla każdej grupy eksperymentalnej dołącz jedną zbiorczą próbkę referencyjną oznaczoną unikalnym znacznikiem TMT, która posłuży jako wspólne odniesienie do porównania z każdą pojedynczą próbką w czterech grupach eksperymentalnych31,32.
- Przeprowadzić frakcjonowanie offline na próbkach zbiorczych za pomocą chromatografii cieczowej z odwróconymi fazami o wysokim pH (RPLC)29,30.
- Analizuj frakcje metodą chromatografii cieczowej tandemowej MS (LC-MS/MS) za pomocą hybrydowego spektrometru mas kwadrupolowo-orbitrapowego29,30. Zapisz wygenerowane dane spektralne MS/MS w formacie Thermo Raw (thermo.raw).
UWAGA: W razie potrzeby pliki Thermo Raw są konwertowane do formatu Mascot Generic Format (.mgf), aby były kompatybilne z różnymi programami. W tym tekście skróty "RAW" i "MGF" oznaczają format pliku wejściowych zestawów danych MS/MS. Na rysunkach zestawy danych MS/MS są dla uproszczenia reprezentowane przez te same ikony RAW.
2. Konfiguracja modułu
UWAGA: Wybrane przyciski/menu są pogrubione. Przykładowe pliki, przepływy pracy i parametry narzędzi są dostępne za pośrednictwem tabel uzupełniających. Więcej informacji na temat korzystania z Galaxy można znaleźć na stronie GTN FAQs (https://training.galaxyproject.org/training-material/faqs/galaxy/).
- Serwer Galaxy Europe
- Uzyskaj dostęp do serwera Galaxy Europe (Galaxy EU; https://usegalaxy.eu/).
- Utwórz konto lub zaloguj się. Do utworzenia nowego konta wymagany jest prawidłowy adres e-mail. Zaloguj się jako użytkownik, aby korzystać z Galaxy.
- Przygotowanie historii galaktyki
- Jeśli użytkownik importuje przykładowe dane wejściowe z Tabeli uzupełniającej 2, wykonaj kroki 2.2.1.1-2.2.1.3.
- Otwórz przykładowe historie galaktyk, korzystając z linków podanych w Tabeli Uzupełniającej 2.
- Kliknij szary przycisk Importuj tę historię znajdujący się w lewym górnym rogu (środkowego) panelu. Zmień nazwę historii i kliknij przycisk Kopiuj historię. W razie potrzeby dodaj ich zestawy danych do tej historii, klikając przycisk Prześlij w skrajnym lewym panelu i dodaj pliki do przesłania.
- Kliknij przycisk Start > Close (Rozpocznij) Zamknij. Przesłane pliki pojawią się w panelu historii po prawej stronie. Poczekaj, aż kolor zestawu danych zmieni kolor na zielony przed użyciem.
UWAGA: W przypadku importowania (kopiowania) istniejącej historii nie należy tworzyć osobnej (nowej) historii.
- Jeśli użytkownik tworzy nową historię i przesyła swoje dane, wykonaj kroki 2.2.2.1.-2.2.2.2.
- W panelu Historia (po prawej stronie) kliknij raz ikonę + (plus), aby utworzyć nową historię o nazwie "Historia nienazwana". Kliknij ikonę ołówka obok historii i kliknij Zapisz. Te same kroki dodawania zestawów danych do istniejącej (przykładowej) historii mają zastosowanie do przesyłania danych.
- W skrajnym lewym panelu kliknij Prześlij i dodaj pliki do przesłania. Kliknij przycisk Start > Close (Rozpocznij) Zamknij. Przesłane pliki pojawią się w nowej historii. Poczekaj, aż kolor zestawów danych zmieni kolor na zielony.
- Jeśli użytkownik analizuje jednocześnie wiele plików MS/MS, wykonaj kroki 2.2.3.1.-2.2.3.3.
- Umieść je w kolekcji zestawów danych, aby wybrać je jako jedne dane wejściowe. Kliknij ikonę znacznika wyboru w panelu Historia i wybierz (sprawdź) zestawy danych.
- Kliknij przycisk z informacją o liczbie wybranych zestawów danych (np. 4 z 8 wybranych), a następnie w menu rozwijanym kliknij opcję Utwórz listę zestawów danych. W wyskakującym okienku wpisz nazwę kolekcji (np. MGF Data, RAW Data). W razie potrzeby wybierz, czy oryginalne zestawy danych zostaną ukryte po utworzeniu kolekcji.
- Kliknij niebieski przycisk Utwórz kolekcję w prawym dolnym rogu wyskakującego okienka. Kliknij ikonę znacznika wyboru w panelu Historia, aby usunąć zaznaczenie zestawów danych.
UWAGA: Każdy z pięciu modułów powinien być uruchomiony we własnej (zaimportowanej lub nowej) historii Galaxy, aby poprawić wrażenia użytkownika. Aby uniknąć nadmiarowości, późniejsze instrukcje modułu pominą konfigurację i skupią się na krokach przepływu pracy.
- Importowanie i uruchamianie przepływu pracy
UWAGA: Zdecydowanie zaleca się, aby wszyscy użytkownicy, niezależnie od tego, czy korzystają z danych przykładowych, czy swoich danych, korzystali i/lub dostosowywali modułowe przepływy pracy z ustalonymi parametrami (Tabela uzupełniająca 2). W ten sposób użytkownicy mogą uniknąć konieczności wyszukiwania i ustawiania parametrów dla każdego narzędzia. W razie potrzeby użytkownicy mogą wyszukać narzędzia, klikając przycisk Narzędzia w skrajnym lewym panelu i wpisując nazwę narzędzia (tak dokładnie, jak to możliwe) w pasku wyszukiwania w sąsiednim panelu. Pasujące narzędzia pojawią się automatycznie. Kliknij właściwy wynik wyszukiwania i ustaw odpowiednie parametry (patrz Plik uzupełniający 1). Przed uruchomieniem narzędzia użytkownicy mogą skonfigurować powiadomienia e-mail, aby powiadamiać ich o zakończeniu zadania, wybierając przycisk na końcu parametrów. Dla wygody dostępne są dwa przyciski Uruchom: jeden w prawym górnym rogu środkowego panelu, a drugi po polach parametrów. Tabela uzupełniająca 3 zawiera dodatkowe zasoby szkoleniowe. Wersje narzędzi i baz danych są aktualne i działają w momencie pisania tego tekstu (czerwiec 2024 r.), ale mogą ulec zmianie wraz z aktualizacją Galaxy oraz powiązanych narzędzi i baz danych.
- Otwórz przepływ pracy w nowej karcie, korzystając z łączy w Tabeli uzupełniającej 2.
- Kliknij przycisk Importuj w prawym górnym rogu panelu. Otworzy się nowa karta z zielonym polem, które potwierdza, że przepływ pracy został zaimportowany. Zielone pole będzie również zawierało dwie opcje: zacznij od razu korzystać z tego przepływu pracy lub wróć do poprzedniej strony.
- Kliknij pierwszy przycisk ("rozpocznij korzystanie z tego przepływu pracy..."), aby otworzyć kartę Przepływ pracy w środkowym panelu interfejsu, w której wyświetlane są wszystkie zapisane przepływy pracy. Znajdź przepływ pracy, który został właśnie zaimportowany, i kliknij niebieski przycisk odtwarzania (trójkąt). Spowoduje to wyświetlenie pól wejściowych.
UWAGA: Dla każdego podanego przepływu pracy pola wejściowe odpowiadają przykładowym danym wejściowym (tabela uzupełniająca 2). Jeśli użytkownik analizuje swoje dane, jego dane wejściowe powinny być odpowiednio nazwane, aby upewnić się, że dla każdego modułu są używane właściwe pliki.
- Jeśli użytkownik chce wyświetlić przepływy pracy na serwerze Galaxy EU, wykonaj kroki 2.3.2.1-2.3.2.4.
- Kliknij przycisk Przepływ pracy w górnym pasku witryny Galaxy. Na tej karcie kliknij kartę podrzędną Moje przepływy pracy, aby wyświetlić wszystkie zaimportowane przepływy pracy. Aby wyświetlić przepływ pracy, kliknij przycisk Edytuj z ikoną ołówka w celu otwarcia edytora przepływu pracy.
- W edytorze przepływu pracy możesz wchodzić w interakcje z przepływem pracy, na przykład klikać i przeciągać w celu reorganizacji, klikać narzędzia, aby je wyświetlić, zmieniać parametry itp. Po wprowadzeniu zmian zapisz edytowany przepływ pracy, klikając ikonę dysku w górnej części prawego panelu, a w razie potrzeby uruchom przepływ pracy, klikając ikonę odtwarzania (również w górnej części prawego panelu).
- Tworzenie przepływów pracy specyficznych dla użytkownika w celu analizowania niestandardowych danych wejściowych. W zależności od wiedzy użytkownika z zakresu metaproteomiki i doświadczenia z platformą Galaxy, zbuduj przepływ pracy, a następnie przeanalizuj dane.
- Jeśli użytkownik jest mniej doświadczony, przetestuj różne narzędzia w historii, a następnie wyodrębnij przepływ pracy z ich ukończonej analizy.
UWAGA: Ten wyodrębniony przepływ pracy można rozszerzać, poprawiać i ponownie wykorzystywać, co pozwala użytkownikom na dokładne odtworzenie swojej pracy. Bardziej szczegółowe instrukcje można znaleźć w sekcji FAQ GTN dotyczącej przepływów pracy (https://training.galaxyproject.org/training-material/faqs/galaxy/#workflows).
- Kliknij każde pole wejściowe i wybierz odpowiednie dane wejściowe. Sekcje od 3 do 7 opisują wejścia modułu. Sprawdź, czy wszystkie dane wejściowe są w akceptowanym formacie, aby uniknąć błędów. Kliknij akceptowane formaty pod każdym polem wejściowym, aby sprawdzić, czy wszystkie pliki są kompatybilne z narzędziami. Po zakończeniu kliknij Uruchom workflow.
UWAGA: Jeśli użytkownik woli skonfigurować narzędzia ręcznie, materiały instruktażowe dla każdego modułu tego przepływu pracy metaproteomiki klinicznej są dostępne na stronie internetowej GTN (https://gxy.io/GTN:P00019). Szacowane czasy wykonania kluczowych narzędzi zostały uwzględnione w tabeli uzupełniającej 2, ale środowiska wykonawcze są zależne od rozmiaru danych wejściowych, zależności narzędzi (takich jak wymagania dotyczące pamięci w porównaniu z przydzieloną pamięcią), zaplanowanych czasów konserwacji, błędów itp. Stany zadań są wskazywane kolorem zestawu danych, a po wybraniu (kliknięciu) zestawu danych pojawi się komunikat informujący, czy zadanie oczekuje na umieszczenie w kolejce (szary), uruchomiony (pomarańczowy) czy zakończony niepowodzeniem (czerwony). Po zakończeniu zadania zestaw danych zmieni kolor na zielony (brak komunikatu potwierdzającego). Użytkownicy mogą wyrazić zgodę na otrzymywanie powiadomień e-mail, które będą ich powiadamiać o zakończeniu zadań (patrz UWAGA na początku kroku 2.3). Poniższe instrukcje dotyczące modułu pominą wyraźne kroki konfiguracji, ponieważ są one takie same dla każdego modułu (patrz sekcja 2 i często zadawane pytania dotyczące GTN, jeśli to konieczne) i opiszą kluczowe narzędzia dla każdego modułu. Pełna lista używanych narzędzi znajduje się w Tabeli Uzupełniającej 1. Nazwy narzędzi zostały pogrubione. Dla porównania, wszystkie nazwy, wersje i opisy narzędzi są zawarte w Tabeli materiałów. Jeśli użytkownik korzysta z przykładowych przepływów pracy z tabeli uzupełniającej 2, zapoznaj się z przykładowymi nazwami plików znajdującymi się w nawiasach na końcu każdego kroku. Jeśli użytkownik uruchamia narzędzia niezależnie, przykładowe nazwy plików mogą zostać zignorowane. Aby zmienić nazwę zestawu danych, kliknij ikonę ołówka w prawym górnym rogu zestawu danych. W polu "Nazwa" wpisz nową nazwę i kliknij Zapisz.
3. Moduł 1: Generowanie bazy danych sekwencji białek
UWAGA: Jeśli użytkownik chce skorzystać z przykładowych danych wejściowych i przepływu pracy z Tabeli Uzupełniającej 2, powinien postępować zgodnie z instrukcjami w sekcji 2. W przypadku modułu 1 zaimportuj dane wejściowe i przepływ pracy dla GENEROWANIA BAZY DANYCH. Kolumna danych wyjściowych w tabeli uzupełniającej 2 zawiera przykłady ukończonych historii wydruków w celach informacyjnych. Dla wszystkich modułów odpowiedni samouczek dotyczący OWS można znaleźć w tabeli uzupełniającej 3.
- Należy sporządzić wykaz gatunków, które są powiązane z chorobą lub stanem będącym przedmiotem zainteresowania i/lub miejscem pobrania próbek.
- Listę tych gatunków można uzyskać z przeglądu literatury. Alternatywnie, jeśli próbki były wcześniej analizowane, należy uzyskać listę gatunków z 16S rRNA lub sekwencjonowania metagenomicznego.
- Zapisz tę listę gatunków jako plik tabelaryczny (np. Species.tabular).
UWAGA: Korzystając z listy gatunków, zostanie wygenerowana duża, obszerna baza danych sekwencji białek znanych mikroorganizmów chorobotwórczych, a przy użyciu MetaNovo, ta duża baza danych, która zawiera miliony sekwencji białek, zostanie następnie zredukowana do łatwiejszej w zarządzaniu bazy danych zawierającej białka obecne w próbkach. Krok redukcji bazy danych ma kluczowe znaczenie, ponieważ wiele narzędzi do przeszukiwania baz danych nie jest w stanie obsłużyć milionów sekwencji. Zmniejszona baza danych zostanie połączona z białkami ludzkimi i zanieczyszczeniami w celu wygenerowania kompaktowej bazy danych w celu wygenerowania kompaktowej bazy danych, która zostanie wykorzystana do identyfikacji peptydów w następnym module (sekcja 4).
- Użyj listy gatunków (Species.tabular) jako danych wejściowych dla UniProt (pobierz proteom jako fasta), aby wygenerować bazę danych sekwencji białek (Species UniProt FASTA.fasta).
- Uruchom narzędzie Protein Database Downloader, aby wygenerować dwie dodatkowe bazy danych sekwencji białek: Human SwissProt (tylko recenzja) i białka zanieczyszczeń (Human SwissProt Protein Database.fasta, Contaminants [cRAP] Protein Database.fasta). Białka zanieczyszczające są również określane jako powszechne repozytorium białek przybyszowych lub cRAP.
- Użyj trzech baz danych białek jako danych wejściowych dla FASTA Merge Files i Filter Unique Sequences, aby wykluczyć duplikaty i wygenerować dużą bazę danych sekwencji białek (Human UniProt Microbial Proteins cRAP for MetaNovo.fasta).
- Użyj dużej (kompleksowej) bazy danych (z kroku 3.4) i zestawów danych MS (MGF) jako danych wejściowych dla MetaNovo33, aby wygenerować zredukowaną bazę danych (MetaNovo Compact Database.fasta).
- Uruchom FASTA Merge Files i Filter Unique Sequences w bazie danych wygenerowanej przez MetaNovo, bazach danych Human SwissProt (tylko recenzowane) i cRAP, aby wygenerować zredukowaną (docelową) bazę danych sekwencji białek mikrobiologicznych, ludzkich i zanieczyszczających, które będą używane do wykrywania peptydów (Human UniProt Microbial Proteins [z MetaNovo] i cRAP.fasta).
4. Moduł 2: Odkrywanie peptydów poprzez przeszukiwanie bazy danych
UWAGA: Jeśli użytkownik chce skorzystać z przykładowych danych wejściowych i przepływu pracy z Tabeli Uzupełniającej 2, koniecznie postępuj zgodnie z instrukcjami w Sekcji 2. W przypadku modułu 2 zaimportuj dane wejściowe i przepływ pracy dla funkcji DISCOVERY. Dla wszystkich modułów odpowiedni samouczek dotyczący OWS można znaleźć w tabeli uzupełniającej 3. SearchGUI34,35,36 i PeptideShaker37 są oddzielnymi programami, ale będą traktowane jako jeden program do identyfikacji i przetwarzania peptydów, ponieważ są używane razem. W celu zapewnienia kompatybilności oprogramowania, zestawy danych MS/MS zostaną przekonwertowane z RAW na MGF dla SearchGUI/PeptideShaker za pomocą narzędzia msconvert (w dostarczonym przepływie pracy). MaxQuant38 może przetwarzać pliki RAW.
- Uruchom FastaCLI, aby dodać sekwencje białek wabiących do zredukowanej (docelowej) bazy danych w celu wygenerowania bazy danych sekwencji białek wabika (FastaCLI MetaNovo Human SwissProt cRAP with decoys.fasta).
UWAGA: FastCLI będzie musiał zostać uruchomiony tylko dla SearchGUI/PeptideShaker. MaxQuant może dodawać wabiki i zanieczyszczenia do bazy danych sekwencji białek. W tym przypadku zredukowana baza danych zawiera już zanieczyszczenia (cRAP), więc MaxQuant został ustawiony tak, aby dodawać tylko wabiki.
- Uruchom SearchGUI/PeptideShaker i MaxQuant, aby przeszukać zestawy danych MS w zmniejszonej bazie danych w celu zidentyfikowania peptydów i ostatecznie przypisania ich do sekwencji białkowych za pomocą przeszukiwania bazy danych sekwencji. Patrz Tabela Uzupełniająca 4, aby zapoznać się z parametrami narzędzia.
UWAGA: W tym przypadku zostaną użyte dwa programy do identyfikacji peptydów (SearchGUI/PeptideShaker i MaxQuant) do identyfikacji sekwencji peptydów i białek poprzez przeszukiwanie bazy danych sekwencji. Programy te identyfikują peptydy w widmach MS/MS i przeszukują bazę danych sekwencji białek, dopasowując obserwowane i teoretyczne dane peptydowe, w tym masy peptydów i widma. W następnym module zidentyfikowane peptydy zostaną zweryfikowane za pomocą PepQuery2 w celu sprawdzenia, czy uzyskano peptydy mikrobiologiczne (sekcja 5).
- Uruchom SearchGUI, aby wygenerować plik archiwum zawierający pliki PSM (Wyszukaj GUI na danych [#].searchgui_archive).
- Użyj pliku archiwum SearchGUI jako danych wejściowych dla PeptideShaker w celu wygenerowania raportu PSM, raportu peptydowego i raportu białkowego (Peptide Shaker on data [#]: [nazwa raportu].tabular).
- Uruchom MaxQuant, aby wygenerować pliki grup białek i peptydów (MaxQuant Protein Groups.tabular, MaxQuant Peptides.tabular).
UWAGA: MaxQuant wymaga eksperymentalnego pliku projektowego, który zawiera warunki eksperymentu, grupy próbek i relacje między próbkami (Experimental Design Discovery MaxQuant.tabular). Ten plik informuje MaxQuant o tym, jak organizować i analizować dane MS. Przykład znajduje się w tabeli uzupełniającej 5. Jeśli używasz danych użytkownika, użytkownicy muszą zmodyfikować ten plik, aby pasował do ich zestawów danych MS.
- Użyj narzędzi do manipulowania tekstem, aby zarządzać danymi wyjściowymi z obu programów. Zapoznaj się z przepływem pracy DISCOVERY w Tabeli uzupełniającej 2, aby zobaczyć, które narzędzia mają zastosowanie do SearchGUI/PeptideShaker i MaxQuant.
UWAGA: W Galaxy zaimplementowano następujące narzędzia do manipulacji tekstem. Najważniejsze narzędzia zostały wyróżnione poniżej, dlatego zdecydowanie zaleca się, aby użytkownicy zapoznali się z przepływem pracy DISCOVERY, aby zobaczyć dodatkowe narzędzia, które nie zostały tutaj omówione. Zapoznaj się z sekcją 2, aby uzyskać instrukcje dotyczące wyświetlania przepływu pracy.
- Wybierz dopasowania mikrobiologiczne (Wybierz drobnoustrojowe PSMs.tabular z SGPS, Wybierz peptydy drobnoustrojne (MQ).tabular).
- Użyj Filter and Query Tabular39, aby wybrać pewne PSM i zapytać o ich numery dostępu do białek (Filter confident microbial PSMs.tabular, zapytaj o wyniki na danych [# i #].tabular).
- Użyj opcji Wytnij, aby wyodrębnić sekwencje peptydowe jako nowy zestaw danych (Wytnij dane [#].tabular).
- Użyj opcji Grupa, aby uzyskać unikatowe wpisy (np. unikalne sekwencje peptydów) dla każdego programu (MQ Peptides.tabular, SGPS Distinct Peptides.tabular).
- Połącz dwie listy peptydów w jeden zestaw danych (SGPS-MQ Peptides.tabular).
- Grupa, aby usunąć zduplikowane sekwencje peptydowe. Ostateczna lista odrębnych peptydów mikrobiologicznych zostanie wykorzystana do weryfikacji PepQuery2 (Distinct Peptides.tabular).
5. Moduł 3: Weryfikacja peptydów mikrobiologicznych
UWAGA: Jeśli użytkownik chce skorzystać z przykładowych danych wejściowych i przepływu pracy z Tabeli Uzupełniającej 2, koniecznie postępuj zgodnie z instrukcjami w Sekcji 2. W przypadku modułu 2 zaimportuj dane wejściowe i przepływ pracy dla WERYFIKACJI. Dla wszystkich modułów odpowiedni samouczek dotyczący OWS można znaleźć w tabeli uzupełniającej 3.
- Użyj następujących danych wejściowych dla PepQuery240,41 Lista odrębnych peptydów drobnoustrojów (Odrębne peptydy dla PepQuery.tabular); Zestawy danych spektralnych MS (MGF); Human UniProt Reference (wraz z izoformami) (Human UniProt+Isoforms FASTA.fasta) i bazami danych sekwencji białek cRAP (cRAP.fasta). Patrz parametry w tabeli uzupełniającej 6.
UWAGA: Weryfikacja obecności peptydów i białek ma kluczowe znaczenie dla uzyskania dokładnych danych i istotnych wglądów w proteom systemu biologicznego. PepQuery2 umożliwia walidację nowych, specyficznych dla choroby peptydów będących przedmiotem zainteresowania z czułością i swoistością. Zidentyfikowane peptydy mikrobiologiczne (z modułu 2) zostaną przeszukane w stosunku do sekwencji białek ludzkich i zanieczyszczających w celu sprawdzenia, czy są one pochodzenia mikrobiologicznego (należy uniknąć błędnego przypisania peptydów ludzkich). Zweryfikowane peptydy zostaną wykorzystane do wygenerowania bazy danych sekwencji zweryfikowanych białek, co jest niezbędne do ograniczenia wprowadzania wyników fałszywie dodatnich podczas kwantyfikacji białek w następnym module (sekcja 6).
- Jeden plik rangi PSM zostanie wygenerowany dla każdego zestawu danych MS/MS używanego jako dane wejściowe (PepQuery2 przy kolekcji [#]: psm_rank.tabular). Uruchom polecenie Zwiń kolekcję na plikach rang PSM, aby utworzyć jeden połączony zestaw danych (Zwiń kolekcję na danych [#] .tabular) i Filtruj, aby zachować pewne PSM (Filtruj według [PSM rank collection].tabular).
- Uruchom polecenie Usuń początek, aby wykluczyć nagłówki kolumn, i Wytnij, aby wyodrębnić zweryfikowane sekwencje peptydów jako nowy zestaw danych.
- Uruchom polecenie Cut w raportach peptydowych z SearchGUI/PeptideShaker i MaxQuant (SGPS Peptide Report.tabular, MaxQuant Peptide Report.tabular), aby wyodrębnić sekwencje peptydów i wpisy białek jako nowy zestaw danych peptyd-białko (dla każdego programu) i usuń początek, aby wykluczyć nagłówki kolumn.
- Połącz sekwencje peptydowe i wpisy białek z obu programów, aby utworzyć nowy (połączony) zestaw danych peptyd-białko.
- Uruchom tabelaryczne zapytanie na połączonym zestawie danych peptyd-białko i zweryfikowane peptydy, aby przypisać zweryfikowane peptydy do powiązanych z nimi wpisów białkowych (Peptide and Protein from Peptide Reports.tabular). Wpisy białek są katalogowane według ich numerów dostępu do białek (znanych również jako identyfikatory UniProt ID).
- Grupa, aby zachować unikalne zweryfikowane peptydy i powiązane z nimi identyfikatory UniProt .
- Uruchom tabelaryczne zapytanie, aby wyodrębnić identyfikatory UniProt (UniProt-ID ze zweryfikowanych Peptides.tabular).
- Umieść identyfikatory UniProt w UniProto, aby uzyskać powiązane z nimi sekwencje białek jako nową bazę danych (UniProt.fasta).
- Uruchom FASTA Merge Files i Filter Unique Sequences w bazie danych sekwencji białek generowanych przez UniProt, bazie danych Human UniProt (wraz z izoformami) i bazach danych zanieczyszczeń, aby wygenerować zweryfikowaną bazę danych, która zostanie wykorzystana do kwantyfikacji peptydów (Baza danych ilościowych dla MaxQuant.fasta).
6. Moduł 4: Kwantyfikacja MaxQuant
UWAGA: Jeśli użytkownik chce skorzystać z przykładowych danych wejściowych i przepływu pracy z Tabeli Uzupełniającej 2, koniecznie postępuj zgodnie z instrukcjami w Sekcji 2. W przypadku modułu 2 zaimportuj dane wejściowe i przepływ pracy dla modułu QUANTIFICATION. Dla wszystkich modułów odpowiedni samouczek dotyczący OWS można znaleźć w tabeli uzupełniającej 3.
- Użyj zweryfikowanej bazy danych sekwencji białek i zestawów danych MS (RAW) jako danych wejściowych dla MaxQuant42.
UWAGA: Pamiętaj, że MaxQuant wymaga eksperymentalnego pliku projektowego i może to być ten sam plik, który jest używany do identyfikacji peptydów (krok 4.2). W razie potrzeby zmień nazwy plików. Zweryfikowana baza danych z poprzedniego modułu jest wymagana w celu zmniejszenia liczby wyników fałszywie dodatnich podczas oznaczania ilościowego białka. Kwantyfikacja białek umożliwia naukowcom pomiar i porównanie obfitości peptydów i białek w próbkach biologicznych. Ten krok jest niezbędny do zrozumienia zróżnicowanej ekspresji białek poprzez uzyskanie wglądu w zmiany ilościowe w różnych warunkach.
- Wygeneruj pliki dowodów, grup białek i peptydów (MaxQuant Evidence.tabular, MaxQuant Protein Groups.tabular, MaxQuant Peptides.tabular).
- Wybierz peptydy mikrobiologiczne z pliku MaxQuant Peptides (Wybierz peptydy drobnoustrojowe.tabular).
- Wytnij tylko sekwencje peptydów mikrobiologicznych (Wytnij dane [#].tabular).
- Grupa, aby uzyskać listę ilościowo oznaczonych peptydów mikrobiologicznych (Quantified Peptides.tabular).
7. Moduł 5: Interpretacja danych
UWAGA: Jeśli użytkownik chce skorzystać z przykładowych danych wejściowych i przepływu pracy z Tabeli Uzupełniającej 2, koniecznie postępuj zgodnie z instrukcjami w Sekcji 2. W przypadku modułu 2 zaimportuj dane wejściowe i przepływ pracy dla INTERPRETACJI DANYCH. Dla wszystkich modułów odpowiedni samouczek dotyczący OWS można znaleźć w tabeli uzupełniającej 3. Dane wyjściowe z kwantyfikacji MaxQuant w poprzednim module zostaną tutaj wykorzystane do adnotacji taksonomicznych i funkcjonalnych przy użyciu Unipept oraz analizy statystycznej przy użyciu MSstatsTMT. Unipept umożliwia naukowcom identyfikację i ilościowe określenie mikroorganizmów w różnych środowiskach i integruje się z publicznymi bazami danych (takimi jak UniProt) w celu pobierania zaktualizowanych adnotacji. MSstatsTMT został zaprojektowany z myślą o solidnej analizie statystycznej ilościowych danych proteomicznych opartych na spektrometrii mas przy użyciu znakowania TMT.
- Użyj listy ilościowo oznaczonych peptydów mikrobiologicznych (Quantified Peptides.tabular) jako danych wejściowych dla Unipept43,44,45, aby wykonać adnotacje taksonomiczne i funkcjonalne. Patrz Tabela uzupełniająca 7, aby zapoznać się z parametrami i listą wyników.
- Interesującymi rezultatami Unipept są drzewo taksonomii drobnoustrojów oraz drzewo białek Komisji ds. Enzymów Mikrobiologicznych (EC) (Microbial Taxonomy Tree.d3_hierarchy, Microbial EC Proteins Tree.d3_hierarchy).
- Aby wyświetlić drzewa, kliknij zestaw danych, aby otworzyć opcje. Kliknij Visualize (4opcja od lewej) > przeglądarkę taksonomii Unipept.
- Aby wyświetlić adnotacje taksonomiczne i funkcjonalne w tabeli (Unipept peptinfo.tabular): kliknij ikonę oka w prawym górnym rogu zestawu danych. Przewiń, aby zobaczyć każdy peptyd w osobnym wierszu i informacje w różnych kolumnach.
- Przed wykonaniem analizy statystycznej przy użyciu MSstatsTMT uruchom polecenie Select w pliku MaxQuant Protein Groups, aby utworzyć dwa nowe zestawy danych: białka drobnoustrojów i białka ludzkie (Microbial Proteins.tabular, Human Proteins.tabular). Białka mają znaczniki taksonomiczne, które wskazują ich pochodzenie.
- Wyklucz białka zanieczyszczeń za pomocą tagu "con_".
- Zachowaj białka mikrobiologiczne i ludzkie, które są oznaczone odpowiednio znacznikami mikrobiologicznymi (np. "_9LACO") i "_HUMAN" (Microbial-Proteins.tabular, Human-Proteins.tabular).
- MSstatsTMT42,46,47 będzie używany do przeprowadzania analizy statystycznej. Jako dane wejściowe należy użyć pliku MaxQuant Evidence (z modułu 4) oraz wybranych białek mikrobiologicznych (lub białek ludzkich) z poprzedniego kroku. Ten przepływ pracy nadaje priorytet białkom mikrobiologicznym, ale oferuje również opcję charakteryzowania białek ludzkich. Patrz Tabela uzupełniająca 8, aby zapoznać się z parametrami i listą wyników.
UWAGA: MSstatsTMT wymaga pliku adnotacji i macierzy porównawczej (znanej również jako macierz kontrastu). Plik adnotacji określi, w jaki sposób kwantyfikacje zostaną połączone, podczas gdy macierz porównawcza pomieści różne grupy próbek. Przykłady tych plików zostały zamieszczone w Tabeli Uzupełniającej 9 i Tabeli Uzupełniającej 10.
- Interesujące wyniki MSstatsTMT to wulkan i wykresy porównawcze dla białek mikrobiologicznych (Microbial Proteins Volcano Plot.pdf, Microbial Proteins Comparison.pdf). Wyświetl wykresy, klikając ikonę oka w prawym górnym rogu zestawu danych.