Wartość badań transkryptomicznych zależy przede wszystkim od jakości wyjściowego materiału biologicznego. Jeśli ekstrakcja RNA zostanie przeprowadzona w optymalnych warunkach, wskaźnik integralności RNA (RIN) zazwyczaj wynosi 7 lub więcej (Rysunek 4A). Konieczność hybrydyzacji 2 µg cDNA na chipie Affymetrix HERV-V2 wiąże się z zastosowaniem procesu amplifikacji. Pomyślny etap amplifikacji prowadzi do uzyskania rozkładu w kształcie dzwonu (Rysunek 4B). Następnie przeprowadza się fragmentację za pomocą DNAse1 w celu ujednolicenia rozkładu wielkości cDNA wokół 100 nukleotydów przed hybrydyzacją (Rysunek 4C). Po hybrydyzacji i skanowaniu (Rysunek 4D) wizualna inspekcja obrazu pozwala sprawdzić, czy siatka jest prawidłowo dopasowana do plamek (Rysunek 4E) oraz czy kontrole hybrydyzacji są spójne (Rysunek 4F). Etap ten jest również przydatny do wykluczenia mikrotablic, w których podczas eksperymentu pojawiły się pęcherzyki powietrza lub błędy.
Po przejściu chipów przez kontrolę jakości QC (Rycina 5) i po normalizacji, analiza statystyczna 5 par dopasowanych próbek RNA z guza i zdrowej tkanki prostaty z Lyon-Sud Hospital doprowadziła do zidentyfikowania 207 zestawów sond (probesets) HERV o wartościach ekspresji różnicowej (p.val <0,05) (Rycina 6A). Aby wesprzeć te wyniki i uzyskać informacje specyficzne dla prostaty, do analizy dodano 35 dodatkowych par dopasowanych próbek (jelita grubego, jajnika, jądra, piersi, płuc i prostaty), a procedura SAM-FDR (FDR = 20%) ostatecznie zidentyfikowała 44 specyficzne dla prostaty zestawy sond HERV. Wśród nich opisano 10 najistotniejszych struktur HERV (Rycina 6B). Dalsze badania kliniczne będą wymagane w celu oceny czułości i swoistości tych potencjalnych biomarkerów.

Rycina 1. Schemat całej procedury od etapu klinicznego (1: prostatektomia wykonana przez lekarza i przygotowanie tkanki przez patomorfologa) po etap laboratoryjny (2-6: przygotowanie próbek, przygotowanie celu, procesy mikroarray), prowadzący do identyfikacji potencjalnych biomarkerów (7: analiza bioinformatyczna mikroarray HERV). Kwasy nukleinowe pochodzące z tkanki prawidłowej przedstawiono na pomarańczowo; kwasy nukleinowe z obszaru guza składają się z mieszaniny kwasów nukleinowych prawidłowych (pomarańczowych) i specyficznych dla nowotworu (czarnych). Kliknij tutaj, aby wyświetlić powiększony obraz.

Rysunek 2. Koncepcja i zawartość chipa HERV-V2: Sekwencje HERV pobrane z ludzkiego genomu są przechowywane w bazie danych zwanej HERV-gDB3, następnie 25-merowe sondy kandydackie przechodzą przez dedykowaną procedurę modelowania hybrydyzacji (EDA+), zanim zostaną ostatecznie zsyntetyzowane na macierzy (wynikowe ukierunkowane podregiony przedstawiono dla każdej rodziny). Kliknij tutaj, aby wyświetlić większy obraz.

Rysunek 3. Przygotowanie prostaty przez patomorfologa. (A) Świeży preparat po radykalnej prostatektomii zostaje przekazany do laboratorium. (B-C) Prostata zostaje oznakowana (kolor zielony po prawej stronie, czarny po lewej stronie). (D) Duży przekrój poprzeczny gruczołu od strony tylnej. (E) Pozostawiając marginesy nienaruszone, wycina się fragmenty tkanek z różnych obszarów gruczołu przedstawczego. (F) Rdzenie tkankowe zostają umieszczone w probówce Eppendorf. (G) Nitka chirurgiczna służy do zamknięcia prostaty w celu zapobiegania zniekształceniu gruczołu i minimalizacji naruszenia marginesu chirurgicznego. Następnie preparat po radykalnej prostatektomii jest gotowy do utrwalenia w formalinie zgodnie ze standardową procedurą analizy histologicznej. Kliknij tutaj, aby wyświetlić powiększony obraz.

Rycina 4. Kontrola jakości preparatyki kwasów nukleinowych i wydajności hybrydyzacji. (A) integralność RNA, (B) amplifikowane cele cDNA oraz (C) pofragmentowane cele użyte na etapie hybrydyzacji. Te trzy kontrole jakości uzyskano za pomocą Bioanalyzer, stosując RNA nano chips oraz test Eukaryote Nano Serie II. (D) Ogólny obraz obszaru hybrydyzacji mikrociągu HERV-V2 po skanowaniu, (E) powiększenie lewego górnego rogu pokazujące kontrole wyrównania siatki oraz (F) powiększenie obszaru centralnego pokazujące kontrole hybrydyzacji naniesionych punktów. Kliknij tutaj, aby wyświetlić powiększony obraz.

Rysunek 5. Przetwarzanie sygnałów. (A) Kontrole amplifikacji polyA spike-in firmy Affymetrix. Kontrole polyA dla transkryptów Dap, Thr, Phe i Lys z genów B. subtilis są dodawane do próbki RNA i służą do oceny ogólnego sukcesu etapów przygotowania celu. Intensywność powinna być wykrywana przy malejących wartościach wśród tych kontroli spike-in, aby upewnić się, że podczas amplifikacji WT-Ovation nie wystąpiło obciążenie (bias) między genami o wysokiej i niskiej ekspresji. (B) Kontrole hybrydyzacji spike-in firmy Affymetrix. Cele te, wyizolowane z E. coli oraz bakteriofaga P1, są dodawane przed procedurą znakowania. Rosnące wartości dla BioB, BioC, BioD i Cre wskazują na ogólny sukces hybrydyzacji. (C) Rozkład intensywności sygnałów z chipa po normalizacji RMA. Większość zestawów sond (probesets) wykazuje sygnały o wartościach niższych niż 26 (tło), co wskazuje na ogólną ekspresję ograniczoną głównie do niektórych specyficznych loci HERV. Kliknij tutaj, aby wyświetlić większy obraz.

Rycina 6. Analiza danych. (A) Analiza hierarchicznego grupowania próbek prawidłowych i nowotworowych. W odniesieniu do znormalizowanych wartości ekspresji zastosowano algorytm grupowania partycyjnego z funkcją odległości euklidesowej, co pozwoliło na pogrupowanie zestawów sond w zależności od nadekspresji (kolor czerwony) i niedoekspresji (kolor niebieski) w próbkach prawidłowych i nowotworowych. (B) Wybór 10 najważniejszych struktur HERV zidentyfikowanych jako potencjalne biomarkery raka prostaty. Dla każdego elementu HERV podano powiązaną rodzinę HERV, współrzędne genomiczne (NCBI 36/hg18) oraz krótki opis struktury HERV. Kliknij tutaj, aby wyświetlić powiększony obraz.

Rysunek 7. Repertuar HERV. (A) Sekwencjonowanie ludzkiego genomu ujawniło 25 000 genów kodujących białka (egzony, 2%) oraz ogromną liczbę elementów transpozycyjnych, w tym 200 000 retrotranspozonów z długimi powtórzeniami końcowymi (LTR) (HERV, 8%). (B) Ekstrapolacja z zawartości chipa HERV-V2 i powiązanych danych o ekspresji (79 próbek pochodzących z 8 typów tkanek prawidłowych w porównaniu z nowotworowymi) sugeruje, że jedna trzecia repertuaru HERV jest aktywna transkrypcyjnie. Kliknij tutaj, aby wyświetlić powiększony obraz.

Rysunek 8. Interpretacja funkcjonalna sygnałów z chipa. (A) Identyfikacja promotora i kontrola epigenetyczna: ujemny sygnał U3 (czerwona sonda, 5'LTR) versus dodatni sygnał R-U5 (niebieska sonda, 5'LTR) sugerują transkrypcję sterowaną przez U3, co potwierdza różny poziom metylacji CpG (pełne czarne kółka) w obrębie U3 w tkankach prawidłowych okołoguzowych versus tkankach nowotworowych. (B) Strategia splicingu: domniemany mRNA kodujący otoczkę o długości 3,1 kb, ekspresowany wyłącznie w guzie, zidentyfikowano za pomocą sondy nakładającej się na złącze splicingowe SD1/SA2. *Wnioskowane na podstawie porównania z innymi tkankami pozapłodowymi. Kliknij tutaj, aby wyświetlić powiększony obraz.