$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Eine Studie über den Ausbruch von C. trachomatis LGV in Buenos Aires, Argentinien, zeigte einige neue Merkmale in den Proben, nämlich 10 Nukleotidsubstitutionen im ompA-Gen im Vergleich zum ompA-Genotyp L1 und eine ungewöhnliche Multi-Locus Sequence Typing (MLST) Sequenz Typ 2 17,18,19. Unser Ziel war es, die Genome dieser klinischen Proben zu untersuchen und sie in den Kontext globaler Daten zu stellen. Um dies zu erreichen, verglichen wir die Leistung eines Target-Enrichment-Ansatzes mit der direkten metagenomischen Sequenzierung für WGS von bakteriellen STIs aus klinischen Proben, um die Herausforderungen zu überwinden, die sich aus der geringen bakteriellen Belastung und der reichlich vorhandenen menschlichen und mikrobiellen DNA ergeben.
Die verwendete DNA stammte entweder aus rektalen Abstrichproben, die zwischen 2017 und 2023 in Buenos Aires, Argentinien, entnommen und in 2-SP-Medium (Saccharose-Phosphat-Medium) gegeben wurden, das mit 2 % fötalem Kälberserum (FCS) und Antibiotika (Gentamicin 25 μg/ml, Vancomycin 0,5 mg/ml und Amphotericin B 2,5 μg/ml) ergänzt und bei -80 °C gelagert wurde, oder gepoolte Harnröhren-, anorektale, und Rachenabstriche, die am Universitätsspital Zürich (USZ) in NAATs (Nucleic Acid Amplification Tests) gesammelt wurden. Die DNA aus den Proben von Buenos Aires wurde mit einem manuellen Extraktionskit auf Siliziumdioxidbasis extrahiert, während die DNA aus den USZ-Proben mit einem automatisierten System auf Basis von magnetischen Beads extrahiert wurde.
Alle klinischen Proben, die in diese Studie eingeschlossen waren, wurden durch einen intern entwickelten Multiplex-Real-Time-PCR-Assay (siehe Materialtabelle) mit qPCR-Zyklus-Schwellenwerten (Ct) wie folgt bestätigt: C. trachomatis: 17-34, N. gonorrhoeae: 24-31, T. pallidum: 27-36 und M. genitalium: 30-37. Die DNA-Konzentration wurde vor der Sequenzierung mit einem hochempfindlichen fluorometrischen Assay quantifiziert, um das Eingangsmaterial der klinischen Proben zu charakterisieren (siehe Materialtabelle).
Um die direkte metagenomische Sequenzierung mit der Zielanreicherung zu vergleichen, verwendeten wir 32 Proben, von denen 29 positiv für C. trachomatis, neun für N. gonorrhoeae, fünf für T. pallidum und drei für M. genitalium waren. Die direkte metagenomische Sequenzierung führte zu keinem vollständigen Genom und erreichte in der Tat kaum Zieldaten, die über dem Ausgangswert der negativen Proben lagen. Die Zielanreicherung mit Hybridisierung bei 65 °C führte zu deutlich besseren Ergebnissen für jeden Erreger, und in vielen Fällen wurden mit dieser Methode erfolgreiche Genome erzielt (Abbildung 1).
Um die optimale Probenauswahl für die Zielanreicherung zu bewerten, untersuchten wir die Beziehung zwischen diagnostischen qPCR-Ct-Werten und der Effizienz der Genomwiederherstellung, gemessen am Prozentsatz der On-Target-Reads (%OTR). Diese Metrik spiegelt den Anteil der Sequenzierungs-Reads wider, die dem Genom des Zielorganismus zugeordnet sind, und wird häufig zur Bewertung der Anreicherungsleistung verwendet. Es wurde eine signifikante Korrelation zwischen der Erregerlast und dem Sequenzierungserfolg für C. trachomatis gefunden, mit einem Ct-Cut-off für den Erfolg um Ct30. Obwohl die Daten für die anderen Zielerreger begrenzt waren, wurde ein ähnlicher Trend beobachtet, bei dem der Erfolg der Genomsequenzierung überwiegend in Proben mit einer Erregerlast unter Ct30 auftrat (Abbildung 2).
Dieses Experiment mit einer Hybridisierungstemperatur von 65 °C ergab Erfolgsraten bei der Genomsequenzierung mit einer angestrebten Anreicherung von 75 % bis 80 % für C. trachomatis, N. gonorrhoeae und T. pallidum, basierend auf Proben mit Ct-Werten unter 30. Die drei M . genitalium-positiven Proben wiesen Ct-Werte von 30-37 auf und erreichten eine geringere Genomabdeckung und mittlere Lesetiefe. Um diese Ergebnisse zu verbessern, experimentierten wir mit der Reduzierung der Hybridisierungstemperatur, was die Hybridisierung aufgrund des niedrigen %G+C-Gehalts im M . genitalium-Genom verbessert. Durch die Verwendung niedrigerer Hybridisierungstemperaturen, wie z. B. der jetzt empfohlenen Temperatur von 62,5 °C, können wir die Wiederherstellung von M. genitalium-Genomen verbessern, ohne die Wiederherstellung anderer Genome zu beeinträchtigen (Abbildung 3). Diese Ergebnisse zeigen, dass die Target-Anreicherung ein effektiver Ansatz ist, um wertvolle genomische Informationen aus klinischen Proben zu gewinnen, und dass sie an eine Vielzahl von Genomen angepasst werden können, auch über ein panelbasiertes Design hinweg.
Die phylogenetische Analyse der erfolgreich angereicherten C. trachomatis-Genome zeigte, dass sich die Mehrzahl der Proben aus Buenos Aires innerhalb der ompA-Genotyp-L2b-Klade gruppiert, was auf eine erfolgreiche Abstammungslinie mit globaler Zirkulation hindeutet. Die Analyse ergab auch das Vorhandensein einer neuen, eigenständigen Abstammungslinie in Argentinien, die LGV verursacht. Die Proben aus dieser Gruppe unterscheiden sich um etwa 600 SNPs von den anderen LGV-Linien, und diese Linie wird als ompA-Genotyp L4 vorgeschlagen (Abbildung 4). Eine detaillierte Beschreibung der in Abbildung 1, Abbildung 2 und Abbildung 4 verwendeten Beispiele finden Sie in der ergänzenden Tabelle 1.
Diese Ergebnisse unterstreichen den Wert der Zielanreicherung bei der Untersuchung der Erregerdiversität und zeigen einen hohen Grad an Erfolg bei der Genomsequenzierung aus klinischen Proben.

Abbildung 1: Anreicherung von Genomen mittels Zielanreicherung im Vergleich zur direkten metagenomischen Sequenzierung. Der Prozentsatz der On-Target-Reads (%OTR), definiert als der Prozentsatz der Sequenzierungs-Reads, die dem Genom des Zielorganismus zugeordnet werden, zeigt den Grad der Anreicherung, der mit der Zielanreicherung erreicht wird. Die Linie stellt den Median aller Datenpunkte dar, und das Feld gibt den Interquartilsabstand (IQR) von 50 % an. Die Farben der Datenpunkte entsprechen der Probenquelle: Hellblau für Argentinien und Rot für die Schweiz. Der Erfolg der Genomsequenzierung, definiert als Abdeckung >95 % und mittlere Lesetiefe >10, wird durch große Punkte dargestellt, während das Scheitern durch kleine Punkte angezeigt wird. Beachten Sie die unterschiedlichen Maßstäbe auf der y-Achse. Abkürzungen: CT = C. trachomatis, NG = N. gonorrhoeae, TP = T. pallidum, MG = M. genitalium. Diese Zahl wurde von10 geändert. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 2: Korrelation von %OTR aus der Zielanreicherung mit diagnostischen qPCR-Ct-Werten. Die Farben der Datenpunkte stellen die Stichprobenquelle dar: Hellblau für Argentinien und Rot für die Schweiz. Der Erfolg der Genomsequenzierung, definiert als Abdeckung >95 % und mittlere Lesetiefe >10, wird durch große Punkte angezeigt, während das Scheitern durch kleine Punkte dargestellt wird. Beachten Sie die unterschiedlichen Maßstäbe auf der y-Achse. Abkürzungen: CT = C. trachomatis, NG = N. gonorrhoeae, TP = T. pallidum, MG = M. genitalium. Diese Zahl wurde von10 geändert. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 3: Einfluss der Hybridisierungstemperatur auf die Wiederfindung vollständiger Genome über alle bakteriellen STI-Krankheitserreger in zehn Proben. Menschliche DNA im Bereich der in klinischen Proben wurde mit bekannten Konzentrationen von bakterieller Ziel-DNA separat (einzeln) und auch insgesamt versetzt, um eine Probe mit mehreren Infektionen (gepoolt) darzustellen. Es wurden zwei Experimente mit jeweils den fünf Erregern/Kombinationen durchgeführt: 0,01 % und 0,03 %, um den ungefähren Prozentsatz der DNA jedes Mikroorganismus im Verhältnis zur Gesamt-DNA in der Probe darzustellen (die verbleibende DNA entspricht der menschlichen DNA). Die Rückgewinnung von M. genitalium wird bei niedrigeren Temperaturen verbessert. C. trachomatis Die DNA wurde in einer niedrigen Konzentration gespickt; weitere Experimente zeigen, dass die Temperaturen die Wiederherstellung des Genoms von C. trachomatis nicht beeinflussten (Daten nicht gezeigt). Abkürzungen: CT = C. trachomatis, MG = M. genitalium, NG = N. gonorrhoeae, TP = T. pallidum. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 4: Rekombinationsadjustierte Phylogenie von LGV-Stämmen aus Argentinien. Die Analyse umfasst alle verfügbaren LGV-Genome mit veröffentlichten Daten. Zu den Metadaten, die rechts neben der Phylogenie angezeigt werden, gehören doppelte Patientenproben, die Bayes'sche Analyse der Populationsstruktur (BAPS; Cluster = Klade, wobei der polyphyletische BAPS-Cluster 6 die genetische Vielfalt der Population darstellt), die Genomsequenzierungsmethode (Daten aus dieser Studie), der ompA-Genotyp , das Land und das Jahr, wie im Schlüssel oben links angegeben. Fehlende Daten werden weiß dargestellt. Die Genomnamen aus dieser Studie sind für Argentinien hellblau und für Finnland dunkelblau hervorgehoben (nicht durch Zielanreicherung erhalten). Diese Farben werden ausschließlich zur Unterscheidung der in dieser Studie generierten Genome verwendet. Bootstrap-Werte von 1000 ultraschnellen Bootstraps werden als Prozentsätze auf den Hauptzweigen angezeigt. Diese Zahl wurde von13 geändert. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.
Ergänzende Tabelle 1: Beschreibung der in der Studie verwendeten Materialien, Reagenzien, Instrumente und Geräte. Die Tabelle enthält detaillierte Informationen über den Handelsnamen, den Hersteller und die Katalognummer der Komponenten, die für die DNA-Extraktion, die Bibliothekserstellung, die Zielanreicherung und die Sequenzierung verwendet werden. Es umfasst auch Artikel, die für hauseigene qPCR-Assays und allgemeine molekularbiologische Verfahren verwendet werden. Gleichwertige Alternativen können für Elemente verwendet werden, die mit einem Sternchen (*) gekennzeichnet sind. Bitte klicken Sie hier, um diese Tabelle herunterzuladen.