$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Badanie nad wybuchem epidemii Lv C. trachomatis w Buenos Aires w Argentynie wykazało pewne nowe cechy w próbkach, a mianowicie 10 podstawień nukleotydów w genie ompA w porównaniu z genotypem L1 ompA oraz niezwykłą sekwencję typu 2 17,18,19 (Multi-Locus Sequence Typing) (MLST). Naszym celem było zbadanie genomów z tych próbek klinicznych, umieszczając je w kontekście danych globalnych. Aby to osiągnąć, porównaliśmy wydajność podejścia do wzbogacania docelowego z bezpośrednim sekwencjonowaniem metagenomicznym dla WGS bakteryjnych chorób przenoszonych drogą płciową z próbek klinicznych, próbując przezwyciężyć wyzwania związane z niskim obciążeniem bakteryjnym oraz obfitym ludzkim i mikrobiologicznym DNA.
Wykorzystane DNA pochodziło z próbek wymazów z odbytu pobranych w latach 2017-2023 w Buenos Aires w Argentynie, umieszczonych w pożywce 2-SP (pożywka sacharozowo-fosforanowa) uzupełnionej 2% płodową surowicą cielęcą (FCS) i antybiotykami (gentamycyna 25 μg/ml, wankomycyna 0,5 mg/ml i amfoterycyna B 2,5 μg/ml) i przechowywane w temperaturze -80 °C, lub zbiorcze wymazy z cewki moczowej, odbytnicy i gardła pobrane w Szpitalu Uniwersyteckim w Zurychu (USZ) w buforze NAAT (Nucleic Acid Amplification Tests). DNA z próbek z Buenos Aires zostało wyekstrahowane przy użyciu zestawu do ekstrakcji ręcznej na bazie krzemionki, podczas gdy DNA z próbek USZ zostało wyekstrahowane za pomocą zautomatyzowanego systemu opartego na kulkach magnetycznych.
Wszystkie próbki kliniczne uwzględnione w tym badaniu zostały potwierdzone dodatnim wynikiem za pomocą opracowanego przez nas multipleksowego testu PCR w czasie rzeczywistym (patrz tabela materiałów) z wartościami progowymi cyklu qPCR (Ct) w następującym zakresie: C. trachomatis: 17-34, N. gonorrhoeae: 24-31, T. pallidum: 27-36 i M. genitalium: 30-37. Stężenie DNA określono ilościowo za pomocą testu fluorometrycznego o wysokiej czułości przed sekwencjonowaniem w celu scharakteryzowania materiału wejściowego próbek klinicznych (patrz Tabela materiałów).
Aby porównać bezpośrednie sekwencjonowanie metagenomiczne ze wzbogaceniem docelowym, użyliśmy 32 próbek, z których 29 było dodatnich dla C. trachomatis, dziewięć dla N. gonorrhoeae, pięć dla T. pallidum i trzy dla M. genitalium. Bezpośrednie sekwencjonowanie metagenomiczne nie dało żadnych kompletnych genomów, a w rzeczywistości ledwo osiągnęło dane docelowe powyżej wartości wyjściowej obserwowanej w próbkach ujemnych. Wzbogacenie celu za pomocą hybrydyzacji w temperaturze 65 °C dało znacznie lepsze wyniki dla każdego patogenu, a w wielu przypadkach uzyskano udane genomy dzięki tej metodzie (ryc. 1).
Aby ocenić optymalny wybór próbki do wzbogacenia docelowego, oceniliśmy związek między diagnostycznymi wartościami qPCR Ct a skutecznością odzyskiwania genomu, mierzoną procentem odczytów na miejscu (%OTR). Wskaźnik ten odzwierciedla proporcję odczytów sekwencjonowania, które odwzorowują genom organizmu docelowego i jest powszechnie stosowany do oceny wydajności wzbogacania. Stwierdzono istotną korelację między obciążeniem patogenem a sukcesem sekwencjonowania C. trachomatis, z punktem odcięcia Ct dla powodzenia około Ct30. Chociaż dane dotyczące innych docelowych patogenów były ograniczone, zaobserwowano podobny trend, w którym sukces sekwencjonowania genomu wystąpił głównie w próbkach z ładunkiem patogenu poniżej Ct30 (ryc. 2).
Ten eksperyment z temperaturą hybrydyzacji 65 °C dał wskaźniki powodzenia sekwencjonowania genomu ze wzbogaceniem docelowym o 75%-80% dla C. trachomatis, N. gonorrhoeae i T. pallidum, na podstawie próbek o wartościach Ct poniżej 30. Trzy próbki dodatnie M. genitalium miały wartości Ct 30-37 i osiągnęły niższe pokrycie genomu oraz średnią głębokość odczytu. Aby poprawić te wyniki, eksperymentowaliśmy z obniżeniem temperatury hybrydyzacji, co poprawia hybrydyzację ze względu na niską zawartość %G+C w genomie M. genitalium . Stosując niższe temperatury hybrydyzacji, takie jak obecnie zalecana temperatura 62,5 °C, możemy poprawić odzyskiwanie genomów M. genitalium bez uszczerbku dla odzyskiwania innych genomów (ryc. 3). Wyniki te pokazują, że wzbogacanie docelowe jest skutecznym podejściem do uzyskiwania cennych informacji genomowych z próbek klinicznych i może być dostosowane do szerokiego zakresu genomów, również w ramach projektu panelowego.
Analiza filogenetyczna pomyślnie wzbogaconych genomów C. trachomatis wykazała, że większość próbek z Buenos Aires skupia się w obrębie kladu ompA-genotyp L2b, co sugeruje udaną linię o globalnym obiegu. Analiza ujawniła również obecność nowej, odrębnej linii rodowej w Argentynie, powodującej LGV. Próbki z tego kladu różnią się o około 600 SNP od innych linii LGV, a linia ta jest proponowana jako genotyp L4 ompA (ryc. 4). Szczegółowy opis próbek użytych na rys. 1, rys. 2 i rys. 4 znajduje się w tabeli uzupełniającej 1.
Odkrycia te podkreślają wartość wzbogacania celu w badaniu różnorodności patogenów i wykazują wysoki stopień sukcesu sekwencjonowania genomu na podstawie próbek klinicznych.

Rysunek 1: Wzbogacanie genomów za pomocą wzbogacania docelowego w porównaniu z bezpośrednim sekwencjonowaniem metagenomicznym. Odsetek odczytów po osiągnięciu celu (%OTR), zdefiniowany jako odsetek odczytów sekwencjonowania mapujących genom organizmu docelowego, pokazuje stopień wzbogacenia osiągnięty dzięki wzbogaceniu docelowemu. Linia reprezentuje medianę wszystkich punktów danych, a pole wskazuje 50% przedział międzykwartylowy (IQR). Kolory punktów danych odpowiadają źródłu próbki: bladoniebieski dla Argentyny i czerwony dla Szwajcarii. Sukces sekwencjonowania genomu, zdefiniowany jako pokrycie >95% i średnia głębokość odczytu >10, jest reprezentowany przez duże punkty, podczas gdy niepowodzenie jest wskazywane przez małe punkty. Zwróć uwagę na różne skale na osi y. Skróty: CT = C. trachomatis, NG = N. gonorrhoeae, TP = T. pallidum, MG = M. genitalium. Liczba ta została zmodyfikowana z10. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Korelacja %OTR ze wzbogacenia celu z diagnostycznymi wartościami qPCR Ct. Kolory punktów danych reprezentują źródło próbki: bladoniebieski dla Argentyny i czerwony dla Szwajcarii. Powodzenie sekwencjonowania genomu, zdefiniowane jako pokrycie >95% i średnia głębokość odczytu >10, jest oznaczone dużymi punktami, podczas gdy niepowodzenie jest reprezentowane przez małe punkty. Zwróć uwagę na różne skale na osi y. Skróty: CT = C. trachomatis, NG = N. gonorrhoeae, TP = T. pallidum, MG = M. genitalium. Liczba ta została zmodyfikowana z10. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rycina 3: Wpływ temperatury hybrydyzacji na odzyskiwanie kompletnych genomów wszystkich patogenów bakteryjnych chorób przenoszonych drogą płciową w dziesięciu próbkach. Ludzkie DNA w zakresie próbek klinicznych zostało wzbogacone znanymi stężeniami docelowego bakteryjnego DNA osobno (pojedyncze), a także razem, aby reprezentować próbkę z wieloma infekcjami (połączone). Przeprowadzono dwa eksperymenty, każdy z pięcioma patogenami/kombinacjami: 0,01% i 0,03%, aby przedstawić przybliżony procent DNA z każdego mikroorganizmu w stosunku do całkowitego DNA w próbce (pozostałe DNA odpowiada ludzkiemu DNA). Regeneracja M. genitalium poprawia się w niższych temperaturach. C. trachomatis (C. trachomatis) DNA zostało wzbogacone w niskim stężeniu; Dalsze eksperymenty pokazują, że temperatury nie wpłynęły na odtworzenie genomów C. trachomatis (dane nie pokazane). Skróty: CT = C. trachomatis, MG = M. genitalium, NG = N. gonorrhoeae, TP = T. pallidum. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rycina 4: filogeneza skorygowana o rekombinację szczepów LGV z Argentyny. Analiza obejmuje wszystkie dostępne genomy LGV wraz z opublikowanymi danymi. Metadane wyświetlane po prawej stronie filogenezy obejmują zduplikowane próbki pacjentów, bayesowską analizę struktury populacji (BAPS; klaster = klad, z polifiletycznym klastrem BAPS 6 reprezentującym różnorodność genetyczną populacji), metodę sekwencjonowania genomu (dane z tego badania), genotyp ompA , kraj i rok, jak wskazano w kluczu w lewym górnym rogu. Brakujące dane są wyświetlane w kolorze białym. Nazwy genomów z tego badania są wyróżnione kolorem jasnoniebieskim dla Argentyny i ciemnoniebieskim dla Finlandii (nie uzyskano ich przez docelowe wzbogacenie). Kolory te są używane wyłącznie do rozróżnienia genomów wygenerowanych w tym badaniu. Wartości bootstrap z 1000 ultraszybkich bootstrapów są wyświetlane jako wartości procentowe w głównych gałęziach. Liczba ta została zmodyfikowana z13. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Tabela uzupełniająca 1: Opis materiałów, odczynników, przyrządów i sprzętu wykorzystanych w badaniu. Tabela zawiera szczegółowe informacje na temat nazwy handlowej, producenta i numeru katalogowego składników używanych do ekstrakcji DNA, przygotowania biblioteki, wzbogacania celu i sekwencjonowania. Obejmuje również pozycje używane do wewnętrznych testów qPCR i ogólnych procedur biologii molekularnej. Równoważne alternatywy mogą być stosowane w przypadku produktów oznaczonych gwiazdką (*). Kliknij tutaj, aby pobrać tę tabelę.