$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Een studie naar de C. trachomatis LGV-uitbraak in Buenos Aires, Argentinië, toonde enkele nieuwe kenmerken in de monsters, namelijk 10 nucleotidesubstituties in het ompA-gen in vergelijking met ompA-genotype L1, en een ongebruikelijke Multi-Locus Sequence Typing (MLST) sequentie type 2 17,18,19. We wilden het genoom van deze klinische monsters onderzoeken en ze in de context van wereldwijde gegevens plaatsen. Om dit te bereiken, vergeleken we de prestaties van een doelverrijkingsbenadering met directe metagenomische sequencing voor WGS van bacteriële soa's uit klinische monsters, in een poging de uitdagingen van lage bacteriële belasting en overvloedig menselijk en microbieel DNA te overwinnen.
Het gebruikte DNA was ofwel afkomstig van rectale uitstrijkjes die tussen 2017 en 2023 in Buenos Aires, Argentinië, waren verzameld, geplaatst in 2-SP-medium (sucrosefosfaatmedium) aangevuld met 2% foetaal kalfsserum (FCS) en antibiotica (gentamicine 25 μg/ml, vancomycine 0,5 mg/ml en amfotericine B 2,5 μg/ml) en bewaard bij -80 °C, of gepoold urethraal, anorectaal, en faryngeale uitstrijkjes verzameld in het Universitair Ziekenhuis Zürich (USZ) in NAAT's (Nucleic Acid Amplification Tests) buffer. DNA uit de monsters van Buenos Aires werd geëxtraheerd met behulp van een handmatige extractiekit op basis van silica, terwijl DNA uit USZ-monsters werd geëxtraheerd met behulp van een geautomatiseerd systeem op basis van magnetische kralen.
Alle klinische monsters die in deze studie zijn opgenomen, werden positief bevestigd door een intern ontwikkelde multiplex real-time PCR-assay (zie materiaaltabel) met qPCR-cyclusdrempelwaarden (Ct) die als volgt variëren: C. trachomatis: 17-34, N. gonorrhoeae: 24-31, T. pallidum: 27-36 en M. genitalium: 30-37. De DNA-concentratie werd gekwantificeerd met behulp van een zeer gevoelige fluorometrische test voorafgaand aan sequencing om het inputmateriaal van de klinische monsters te karakteriseren (zie Tabel met materialen).
Om directe metagenomische sequencing te vergelijken met doelverrijking, gebruikten we 32 monsters, waarvan er 29 positief waren voor C. trachomatis, negen voor N. gonorrhoeae, vijf voor T. pallidum en drie voor M. genitalium. Directe metagenomische sequencing resulteerde niet in volledige genomen en bereikte inderdaad nauwelijks doelgegevens boven de basislijn die werd gezien in negatieve monsters. Doelverrijking met hybridisatie bij 65 °C resulteerde in veel betere resultaten voor elke ziekteverwekker, en in veel gevallen werden succesvolle genomen verkregen uit deze methode (Figuur 1).
Om de optimale steekproefselectie voor doelverrijking te beoordelen, evalueerden we de relatie tussen diagnostische qPCR Ct-waarden en de efficiëntie van genoomherstel, zoals gemeten aan de hand van het percentage on-target reads (%OTR). Deze metriek weerspiegelt het aandeel van de sequencing-lezingen die overeenkomen met het genoom van het doelorganisme en wordt vaak gebruikt om de verrijkingsprestaties te beoordelen. Er werd een significante correlatie gevonden tussen de ziekteverwekkerbelasting en het succes van de sequencing voor C. trachomatis, met een Ct-cut-off voor succes rond Ct30. Hoewel de gegevens voor de andere doelpathogenen beperkt waren, werd een vergelijkbare trend waargenomen, waarbij succes bij genoomsequencing voornamelijk optrad in monsters met een pathogene belasting onder Ct30 (Figuur 2).
Dit experiment met een hybridisatietemperatuur van 65 °C gaf genoomsequencing succespercentages met doelverrijking van 75%-80% voor C. trachomatis, N. gonorrhoeae en T. pallidum, op basis van monsters met Ct-waarden onder de 30. De drie M. genitalium-positieve monsters hadden Ct-waarden van 30-37 en bereikten een lagere genoomdekking en gemiddelde leesdiepte. Om deze resultaten te verbeteren, hebben we geëxperimenteerd met het verlagen van de hybridisatietemperatuur, wat de hybridisatie verbetert vanwege het lage %G+C-gehalte in het genoom van M. genitalium . Met behulp van lagere hybridisatietemperaturen, zoals de nu aanbevolen temperatuur van 62,5 °C, kunnen we het herstel van M. genitalium-genomen verbeteren zonder het herstel van andere genomen in gevaar te brengen (Figuur 3). Deze resultaten tonen aan dat doelverrijking een effectieve benadering is voor het verkrijgen van waardevolle genomische informatie uit klinische monsters en kan worden aangepast aan een breed scala aan genomen, ook in een panelgebaseerd ontwerp.
Fylogenetische analyse van de succesvol verrijkte C. trachomatis-genomen toonde aan dat de meerderheid van de monsters uit Buenos Aires clusteren binnen de ompA-genotype L2b-clade, wat wijst op een succesvolle afstamming met wereldwijde circulatie. De analyse onthulde ook de aanwezigheid van een nieuwe, aparte afstamming in Argentinië, waardoor LGV ontstond. De monsters van deze clade verschillen ongeveer 600 SNP's van de andere LGV-lijnen, en deze afstamming wordt voorgesteld als ompA-genotype L4 (Figuur 4). Een gedetailleerde beschrijving van de monsters die zijn gebruikt in figuur 1, figuur 2 en figuur 4 wordt gegeven in aanvullende tabel 1.
Deze bevindingen onderstrepen de waarde van doelverrijking bij het bestuderen van de diversiteit van pathogenen en tonen een hoge mate van succes bij genoomsequencing aan van klinische monsters.

Figuur 1: Verrijking van genomen met behulp van doelverrijking in vergelijking met directe metagenomische sequencing. Het percentage on-target reads (%OTR), gedefinieerd als het percentage sequentie-reads dat in kaart wordt gebracht op het genoom van het doelorganisme, geeft de mate van verrijking weer die is bereikt met doelverrijking. De lijn geeft de mediaan van alle gegevenspunten weer en het vak geeft het interkwartielbereik (IQR) van 50% aan. De kleuren van de datapunten komen overeen met de steekproefbron: lichtblauw voor Argentinië en rood voor Zwitserland. Het succes van genoomsequencing, gedefinieerd als dekking >95% en gemiddelde leesdiepte >10, wordt weergegeven door grote punten, terwijl falen wordt aangegeven door kleine punten. Let op de verschillende schalen op de y-as. Afkortingen: CT = C. trachomatis, NG = N. gonorrhoeae, TP = T. pallidum, MG = M. genitalium. Dit cijfer is gewijzigd van10. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2: Correlatie van %OTR van doelverrijking met diagnostische qPCR Ct-waarden. De kleuren van het gegevenspunt geven de steekproefbron weer: lichtblauw voor Argentinië en rood voor Zwitserland. Het succes van genoomsequencing, gedefinieerd als dekking >95% en gemiddelde leesdiepte >10, wordt aangegeven door grote punten, terwijl falen wordt weergegeven door kleine punten. Let op de verschillende schalen op de y-as. Afkortingen: CT = C. trachomatis, NG = N. gonorrhoeae, TP = T. pallidum, MG = M. genitalium. Dit cijfer is gewijzigd van10. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3: Effect van hybridisatietemperatuur op het herstel van volledige genomen van alle soa-bacteriële pathogenen in tien monsters. Menselijk DNA in het bereik van dat in klinische monsters werd verrijkt met bekende concentraties van bacterieel doel-DNA afzonderlijk (enkelvoudig) en ook in totaal om een monster met meerdere infecties te vertegenwoordigen (gepoold). Er werden twee experimenten uitgevoerd, elk met de vijf pathogenen/combinaties: 0,01% en 0,03% om het geschatte percentage DNA van elk micro-organisme weer te geven ten opzichte van het totale DNA in het monster (het resterende DNA komt overeen met menselijk DNA). Het herstel van M. genitalium wordt verbeterd bij lagere temperaturen. C. trachomatis DNA was in een lage concentratie geprikt; verdere experimenten tonen aan dat de temperaturen geen invloed hadden op het herstel van het genoom van C. trachomatis (gegevens niet getoond). Afkortingen: CT = C. trachomatis, MG = M. genitalium, NG = N. gonorrhoeae, TP = T. pallidum. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4: Voor recombinatie gecorrigeerde fylogenie van LGV-stammen uit Argentinië. De analyse omvat alle beschikbare LGV-genomen met gepubliceerde gegevens. Metagegevens die rechts van de fylogenie worden weergegeven, omvatten dubbele patiëntmonsters, Bayesiaanse analyse van de populatiestructuur (BAPS; cluster = clade, waarbij de polyfyletische BAPS-cluster 6 de genetische diversiteit van de populatie vertegenwoordigt), genoomsequencingmethode (gegevens uit deze studie), ompA-genotype , land en jaar, zoals aangegeven in de sleutel linksboven. Ontbrekende gegevens worden in het wit weergegeven. Genoomnamen uit deze studie zijn gemarkeerd in lichtblauw voor Argentinië en donkerblauw voor Finland (niet verkregen door doelverrijking). Deze kleuren worden uitsluitend gebruikt om de genomen te onderscheiden die in deze studie zijn gegenereerd. Bootstrap-waarden van 1000 ultrasnelle bootstraps worden weergegeven als percentages op de belangrijkste vestigingen. Dit cijfer is gewijzigd van13. Klik hier om een grotere versie van deze figuur te bekijken.
Aanvullende tabel 1: Beschrijving van de materialen, reagentia, instrumenten en apparatuur die in het onderzoek zijn gebruikt. De tabel bevat gedetailleerde informatie over de handelsnaam, de fabrikant en het catalogusnummer van de componenten die worden gebruikt voor DNA-extractie, bibliotheekvoorbereiding, doelverrijking en sequencing. Het omvat ook items die worden gebruikt voor interne qPCR-assays en algemene moleculair-biologische procedures. Gelijkwaardige alternatieven kunnen worden gebruikt voor items die zijn gemarkeerd met een asterisk (*). Klik hier om deze tabel te downloaden.