Diese einzentrische retrospektive Kohortenstudie wurde gemäß der STROBE-Erklärung14 konzipiert und berichtet. Die Studie verwendete zuvor erhobene ambulante Krankenakten, intraorale Fotografien und digitale Bilddaten, ohne zusätzliche Untersuchungen oder Behandlungen erforderlich zu machen. Das Studienprotokoll wurde von dem Ethikkomitee des Sinopharm Dongfeng Stomatological Hospital, Hubei University of Medicine (Genehmigungsnummer LW-2026-01) genehmigt. Da alle Daten vor der Analyse anonymisiert worden waren, verzichtete die Ethikkommission auf die Einholung einer informierten Einwilligung.
Fallquelle und Studienpopulation
Die Studie wurde in der Abteilung für Präventive Zahnheilkunde am Sinopharm Dongfeng Stomatological Hospital der Medizinischen Universität Hubei durchgeführt. Der retrospektive Studienzeitraum umfasste aufeinanderfolgende Patienten, die die Indexbehandlung zwischen dem 1. Januar 2022 und dem 31. Dezember 2023 erhielten. Die Nachbeobachtungsdaten wurden bis zum 28. Februar 2025 abgerufen, und die Datenbank wurde am 31. März 2025 gesperrt. Die Daten wurden aus dem elektronischen Gesundheitsakten-System, dem Behandlungsrekorde-System, dem Abrechnungscode-System, dem digitalen Bildarchiv sowie dem intraoralen Fotografiesystem gewonnen. Als Indexdatum wurde der Tag definiert, an dem die qualifizierende minimalinvasive Behandlung abgeschlossen wurde. Die Analyseeinheit war die Indexokklusionsläsion des betroffenen Zahns und nicht der Patient.
Die Therapiezuweisung erfolgte nicht randomisiert und wurde nicht zu Forschungszwecken festgelegt. Vielmehr wurde die Therapiewahl vom behandelnden Zahnarzt im Rahmen der routinemäßigen klinischen Praxis getroffen und nicht systematisch als unabhängige Entscheidungsvariable dokumentiert. Läsionsmerkmale, Durchbruchstatus, die Machbarkeit der Feuchtigkeitskontrolle, die Präferenz des Behandlers sowie andere klinische Überlegungen konnten alle die Therapiewahl beeinflussen. Daher konnten unterschiedliche in Frage kommende Zähne desselben Patienten unterschiedliche Behandlungsstrategien erhalten, und jeder Zahn wurde entsprechend der tatsächlich durchgeführten Behandlung klassifiziert. Die ungleichen Gruppengrößen spiegelten die relative Häufigkeit wider, mit der die beiden Behandlungsstrategien in der routinemäßigen klinischen Praxis angewandt wurden, und nicht eine geplante 1:1-Zuordnung. Bei den statistischen Analysen wurde berücksichtigt, dass die Daten auf Patientenebene geclustert waren.
Die Studienpopulation bestand aus bleibenden ersten oder zweiten Molaren bei Patienten im Alter von 6 bis 18 Jahren. Zahnbezüge erfüllten alle folgenden Einschlusskriterien: Bei der Indexuntersuchung wurde die Zielokklusionsfläche nach Reinigung, Isolation und kontinuierlichem Trocknen mit Luftstrom für 5 s als International Caries Detection and Assessment System (ICDAS) Code 1 oder 2 diagnostiziert; die digitale Röntgenaufnahme im Seitenzahnbiss zeigte keine dentinäre Durchleuchtung an der Okklusionsfläche im Basisbefund; und entweder eine Harzinfiltration in Kombination mit einer pit-and-fissure-Versiegelung mittels fließfähigem Kunststoffkomposit oder eine alleinige pit-and-fissure-Versiegelung mittels fließfähigem Kunststoffkomposit wurde während desselben Behandlungstermins abgeschlossen.
Bei Zähnen ohne ein Ergebnisereignis innerhalb von 12 Monaten war eine Nachuntersuchung im selben Krankenhaus zwischen 10 und 14 Monaten nach der Indexbehandlung erforderlich, wobei interpretierbare klinische Unterlagen, intraorale Fotografien und digitale Bissflügel-Aufnahmen vorliegen mussten. Bei Zähnen, die innerhalb von 12 Monaten aufgrund einer Progression der Karies an derselben Okklusionsoberfläche einer operativen Wiederherstellung unterzogen wurden, diente die Behandlungsdokumentation der Restauration als Nachweis des Ergebnisses, und es war keine zusätzliche bildgebende Endpunktuntersuchung erforderlich.
Die Ausschlusskriterien waren wie folgt: vorherige Infiltration mit Harz, Fissurenversiegelung, Restauration, Kronenrestauration oder Wurzelkanalbehandlung an der Index-Occlusalfläche; lokalisierte Zahnschmelzdefekte, ein eindeutiger Kariesbefund, ein darunterliegender Dentinschatten oder radiologische Dentintransluzenz im Ausgangszustand; begleitende Enamelhypoplasie, von mittelschwer bis schwer ausgeprägte strukturelle Defekte im Rahmen einer Fluorose oder traumatische Kronendefekte; orthodontische Attachments, die die Index-Occlusalfläche bedeckten; Unfähigkeit, die primäre Exposition zu bestätigen; oder nicht auswertbare Nachuntersuchungsdaten. Wurden während des Studienzeitraums mehrere in Frage kommende Behandlungsakten für denselben Zahn identifiziert, so wurde lediglich die früheste geeignete Behandlung als Indexbehandlung berücksichtigt.
Diagnose und Klassifizierungskriterien von Läsionen
Die Diagnose der Basisläsion basierte auf einer standardisierten klinischen Untersuchung, intraoralen Fotografien und digitalen Bissflügel-Aufnahmen, die am Tag des Indexbesuchs angefertigt wurden. Digitale Bissflügel-Aufnahmen zur Basisuntersuchung und zum Follow-up wurden mit dem gleichen Modell eines komplementären metall-oxid-halbleiterbasierten (CMOS) intraoralen Sensors gemäß einem standardisierten klinischen Protokoll (70 kVp, 7 mA und einer Belichtungszeit von 0,10–0,14 s, abhängig von der Zahngröße und dem Bogen) angefertigt. Zur Standardisierung der Sensorpositionierung und der Röntgenstrahlenausrichtung kam ein Bissflügel-Positionierungssystem bestehend aus einem Sensorhalter, Bissblock, Ausrichtungsarm und Zielring zum Einsatz. Der Sensor wurde parallel zu den Kronen der Indexmolaren positioniert, und der zentrale Röntgenstrahl wurde durch die interproximalen Kontaktpunkte geleitet. Bei der Nachuntersuchung wurden, soweit technisch möglich, derselbe Sensortyp, dasselbe Positionierungssystem, dasselbe Belichtungsprotokoll und dieselben Software-Anzeigeeinstellungen verwendet. Aufnahmen mit erheblichem proximalem Kontaktoverlap, geometrischer Verzerrung oder unvollständiger Darstellung der Indexokklusionsfläche galten als nicht auswertbar. Während des gesamten Studienzeitraums wurden routinemäßige tägliche Systemkontrollen sowie monatliche Qualitätskontrollverfahren durchgeführt.
Vor der klinischen Untersuchung wurde die Zahnoberfläche mit einer fluoridfreien Prophylaxepaste gereinigt. Nach der Isolation wurde der Zahn kontinuierlich 5 s lang mit Luft getrocknet, und die Läsionsbeurteilung erfolgte unter standardisiertem Behandlungslicht durch Zahnärzte, die im International Caries Detection and Assessment System (ICDAS) geschult waren. Gegebenenfalls wurde eine Kugelendsonde lediglich verwendet, um Ablagerungen zu entfernen und die Kontinuität der Oberfläche zu überprüfen; eine scharfe Sonde wurde nicht eingesetzt. ICDAS-Kode 1 wurde definiert als die erste sichtbare Veränderung des Schmelzes, die erst nach dem Trocknen mit Luft erkennbar wurde, einschließlich einer kalkig-weißen oder verfärbten Läsion. ICDAS-Kode 2 wurde definiert als eine deutliche sichtbare Veränderung der Schmelzfarbe oder Durchsichtigkeit, die bereits im feuchten Zustand sichtbar war und sich nach dem Trocknen mit Luft verstärkte15. Beide ICDAS-Kodes setzten eine intakte Schmelzoberfläche voraus, ohne lokalisierte Schmelzdefekte oder sichere Kavitation.
Da Bitewing-Aufnahmen für okkludale Läsionen nicht die gleiche anatomische Präzision wie für proximale Läsionen bieten, wurden die radiographischen Kategorien E0, E1, E2, D1, D2 und D3 nur als vorgegebene, studienspezifische ordinale Beschreibungen des radiographischen Erscheinungsbilds und nicht als direkte Maße der histologischen Läsionstiefe verwendet. E0 bedeutete keine sichtbare okkludale Durchsichtigkeit; E1 und E2 deuteten auf eine Durchsichtigkeit hin, die offenbar auf die äußere bzw. innere Hälfte des Schmelzes beschränkt war; D1, D2 und D3 deuteten auf eine Durchsichtigkeit hin, die offenbar in das äußere, mittlere bzw. innere Drittel des Dentins reichte16. Infrage kommende Läsionen wurden klinisch als ICDAS-Code 1 oder 2 klassifiziert und zeigten keinen radiographischen Hinweis auf eine Beteiligung des Dentins, was den studienspezifischen radiographischen Kategorien E0–E2 entsprach.
Wenn Unstimmigkeiten zwischen den klinischen Befunden, intraoralen Fotografien und den radiologischen Befunden bestanden, überprüften zwei Beurteiler den Fall unabhängig voneinander. Bestand nach der unabhängigen Überprüfung weiterhin Uneinigkeit, wurde der Fall gemäß dem im Folgenden beschriebenen vorab festgelegten Schlichtungsverfahren einem dritten erfahrenen Zahnarzt zur Entscheidung vorgelegt. Fälle, die nach der Schlichtung weiterhin nicht interpretierbar waren, wurden von der Studie ausgeschlossen.
Expositions- und Behandlungsprotokolle
Die interessierende Exposition war die mikroinvasive Behandlungsstrategie, die beim Erstbesuch angewendet wurde. Alle Eingriffe wurden gemäß eines standardisierten Abteilungsprotokolls von drei Zahnärzten mit mehr als fünf Jahren Erfahrung in präventiver und restaurativer Zahnheilkunde durchgeführt: Operator 1 (X.Q.), Operator 2 (Z.Z.) und Operator 3 (J.Z.). Vor Studienbeginn erhielten alle Operateure eine standardisierte Schulung zu den klinischen Verfahren und Behandlungsprotokollen, die in dieser Studie verwendet wurden.
Vor der Behandlung wurde die Zahnoberfläche mit einer fluoridfreien Prophylaxepaste gereinigt. Die Isolation erfolgte entweder mittels eines Kofferdams oder mit Watte-Rollen in Kombination mit einer Hochvakuumsaugung. Die Kofferdam-Isolation wurde bevorzugt, wenn der Durchbruchstatus und die Zahnposition eine stabile Klammerfixierung und vollständige Feldabdeckung ermöglichten. Watte-Rollen in Kombination mit Hochvakuumsaugung wurden eingesetzt, wenn ein teilweiser Durchbruch, die Zahnposition oder eine unzureichende Klammerstabilität eine zuverlässige Kofferdam-Platzierung ausschlossen. Die Isolationsmethode wurde für jeden betroffenen Zahn dokumentiert und als vorgegebene Kovariable berücksichtigt, da die Feuchtigkeitskontrolle sowohl die Behandlungswahl als auch die Materialeigenschaften beeinflussen könnte. Die Verteilung der Isolationsmethoden zwischen den Behandlungsgruppen ist in Tabelle 1 dargestellt. Alle Lichtaushärtungsverfahren wurden mit einem LED-Härtungsgerät mit einer Ausgangsintensität von mindestens 1000 mW/cm2 durchgeführt.
In der kombinierten Behandlungsgruppe wurde ein 15%iges Vorbehandlungs-Gel aus Salzsäure 120 s lang auf die okklusale Punkt-und-Fissur-Läsion aufgetragen, anschließend 30 s lang mit Wasser abgespült und luftgetrocknet. Danach wurde die Läsion 30 s lang mit 99%igem Ethanol dehydriert. Anschließend wurde ein niedrigviskoses, lichtgehärtetes Harzinfiltrat 3 Minuten lang aufgetragen und 40 s lang belichtet. Die Infiltration wurde erneut 1 Minute lang durchgeführt, gefolgt von einem zweiten 40-s-Lichthärtungszyklus. Nach Abschluss der Harzinfiltration wurden der Fissureingang und das angrenzende Schmelzgewebe 20 s lang mit 37%iger Phosphorsäure geätzt, gründlich gespült und luftgetrocknet. Danach wurde ein universeller Haftvermittler aufgetragen und 20 s lang belichtet, woraufhin das Fissursystem mit einem fließfähigen Kompositversiegelungsharz verschlossen wurde. Ein zahnärztlicher Prüfer wurde verwendet, um das Eindringen des Materials in die Fissuren zu unterstützen und eingeschlossene Luftblasen vor der 40-s-Belichtung zu entfernen. Abschließend wurden die Randanpassung und okklusalen Kontakte überprüft und gegebenenfalls angepasst.
In der nur-Versiegelung-Gruppe wurde keine Harzinfiltration durchgeführt. Stattdessen wurden die Fissuren und das angrenzende Schmelzmaterial 20 s lang direkt mit 37%iger Phosphorsäure geätzt, anschließend gespült und luftgetrocknet. Der gleiche universelle Haftvermittler wurde aufgetragen und 20 s lang lichtgehärtet, und derselbe fließfähige Kompositharz wurde zur vollständigen Fissurenversiegelung verwendet. Die endgültige Lichtaushärtung, die Beurteilung der Randanpassung sowie die okklusale Anpassung wurden mit denselben Verfahren wie in der kombinierten Behandlungsgruppe durchgeführt.
Nach der Behandlung erhielten die Patienten beider Gruppen identische Anweisungen zur Mundhygiene und wurden angewiesen, zweimal täglich mit einer fluoridhaltigen Zahnpasta mit 1450 ppm Fluorid zu putzen. Die Einhaltung dieser Mundhygieneanweisungen wurde aufgrund des retrospektiven Studiendesigns nicht systematisch überwacht. Ein zusätzlicher Versiegelungsschritt, der innerhalb von 30 Tagen nach der Indexbehandlung durchgeführt wurde, um marginale Bläschen oder lokal begrenzte Unterfüllungen zu korrigieren, galt als Teil des ursprünglichen Behandlungsverfahrens und wurde nicht als Ergebnisereignis gewertet.
Basisdaten und mögliche Störfaktoren
Mögliche Störfaktoren wurden vor der Datenauswertung im Voraus festgelegt. Patientenspezifische Variablen umfassten das Alter, das Geschlecht, die Gesamt-Dmft+DMFT-Skala und den sichtbaren Plaqueindex. Die Gesamt-Dmft+DMFT-Skala wurde anhand der vollständigen Mundaufzeichnung über kariöse, fehlende und gefüllte Zähne am Tag des Indexbesuchs ermittelt, und der sichtbare Plaqueindex wurde anhand der Ganzmunduntersuchung erhoben, die während desselben Besuchs durchgeführt wurde.
Zahnspezifische Variablen umfassten den Kieferbogen (maxillär oder mandibulär), den Molartyp (erster oder zweiter bleibender Molar), den Durchbruchstatus, die ICDAS-Klassifikation zum Studienbeginn sowie die zu Studienbeginn spezifische radiologische Erscheinungskategorie. Operator-spezifische Variablen enthielten die Identität des Operators und die Isolationsmethode.
Vollständige Durchbruch wurde definiert als die vollständige okklusale Fläche oberhalb des gingivalen Randes, wodurch eine vollständige Isolation ermöglicht wurde. Teilweiser Durchbruch wurde definiert als das Bestehen einer gingivalen Operkulum, das einen Teil der okklusalen Fläche bedeckt. Der sichtbare Plaque-Index wurde als Anteil der untersuchten Zahnoberflächen berechnet, die beim Indexbesuch sichtbaren Plaque aufwiesen. Alle Kovariaten wurden zum Zeitpunkt der Indexbehandlung ermittelt, und keine Nachverfolgungsdaten wurden zur Definition der Basislinienvariablen herangezogen.
Nachuntersuchung und Outcome-Bewertung
Die routinemäßige ambulante Nachsorge wurde von den behandelnden Zahnärzten etwa 6 und 12 Monate nach der Behandlung im Rahmen der standardmäßigen klinischen Versorgung durchgeführt. Der Beobachtungszeitraum dieser Studie betrug 12 Monate nach der Indexbehandlung. Bei betroffenen Zähnen ohne Anzeichen einer Läsionsprogression während des Beobachtungszeitraums wurde der Nachsorgetermin, der dem 12-Monats-Zeitpunkt innerhalb eines vorgegebenen Fensters von 10–14 Monaten am nächsten lag, als Endpunkttermin festgelegt. Bei betroffenen Zähnen, die aufgrund einer Progression kariöser Läsionen an derselben okklusalen Fläche innerhalb von 12 Monaten einer operativen Wiederherstellung unterzogen wurden, diente der erste Termin, bei dem die Läsionsprogression dokumentiert wurde, als Zeitpunkt zur Erfassung des Ergebnisses. Betroffene Zähne ohne dokumentiertes Ergebnis innerhalb von 12 Monaten und ohne qualifizierenden Nachsorgetermin zwischen 10 und 14 Monaten wurden als verloren für die Nachsorge eingestuft. Der 6-Monats-Termin diente hauptsächlich einer routinemäßigen klinischen Überprüfung. Die Befunde dieses Termins flossen in die anschließende klinische Behandlung ein und wurden zur Ergebniserfassung herangezogen, falls eine Läsionsprogression festgestellt wurde; andernfalls wurde das primäre Ergebnis anhand des Nachsorgetermins ermittelt, der dem 12-Monats-Zeitpunkt am nächsten lag.
Das primäre Ergebnis war das Fortschreiten von Läsionen innerhalb von 12 Monaten. Ein Fortschreiten der Läsion wurde definiert als das Eintreten eines der folgenden Ereignisse während des Beobachtungszeitraums: (1) Entwicklung einer eindeutigen Kavität oder eines lokalisierten Zahnschmelzdefekts am Rand der Versiegelung oder auf der entsprechenden okklusalen Fläche unterhalb der Versiegelung, wobei die Läsion im klinischen Befund als kariös bestätigt wurde; (2) Auftreten einer neuen Röntgendurchsichtigkeit, die sich von der Schmelzschicht bis in das Dentin erstreckt, auf einer digitalen Bissflügel-Aufnahme der Index-Okklusionsfläche, die für Studienzwecke als D1, D2 oder D3 kategorisiert wurde; oder (3) Anbringen einer operativen Füllung auf der okklusalen Fläche aufgrund des Fortschreitens der kariösen Läsion.
Ein teilweiser oder vollständiger Verlust des Versiegelungsmaterials ohne Vorliegen einer der oben genannten Befunde wurde nicht als Fortschreiten der Läsion angesehen. Ebenso wurden Füllungen, die aufgrund von Trauma, okklusaler Störung oder frakturiertem nicht-kariösem Material eingebracht wurden, nicht als Ergebnisereignisse klassifiziert.
Datenextraktion und Qualitätskontrolle
Die Datenauswertung erfolgte mithilfe eines vorab festgelegten Datenwörterbuchs und standardisierter Erfassungsformulare. Zwei Untersucher (X.Q. und Z.Z.) extrahierten die Daten unabhängig voneinander aus dem elektronischen Krankenakten-System, dem Behandlungs-Aufzeichnungssystem, dem Abrechnungs-Codesystem, dem digitalen Bildarchiv und dem intraoralen Fotografiesystem. Die Studiendatenbank wurde mittels doppelter Dateneingabe angelegt, und alle Unstimmigkeiten wurden durch Überprüfung der ursprünglichen Quelldokumente behoben.
Die Zuordnung der Behandlung musste sowohl durch die Behandlungsunterlagen als auch durch die entsprechenden Abrechnungscodes bestätigt werden. Bei Unstimmigkeiten zwischen diesen Quellen wurden die Verfahrensunterlagen und intraoperativen Fotografien überprüft. Fälle, bei denen die Behandlungszuordnung nach der Überprüfung nicht bestätigt werden konnte, wurden ausgeschlossen.
Vor der Beurteilung wurden alle Bilddateien anonymisiert, mit Studienkennnummern neu beschriftet und von Behandlungsakten sowie Abrechnungsinformationen getrennt. Die erste Interpretation der Bilder erfolgte ohne Zugriff auf Patientenidentifikatoren oder dokumentierte Behandlungszuweisungen. Wenn zur Klärung des Ergebnisses ein Einblick in die klinischen Unterlagen erforderlich war, erfolgte dieser erst nach Abschluss der initialen Bildinterpretation.
Zwei gleichmäßig geschulte Zahnärzte (X.Q. und Z.Z.) bewerteten unabhängig voneinander die basale ICDAS-Klassifikation, die studienspezifische radiologische Erscheinungskategorie sowie die Ergebnisse der Nachuntersuchung. Beide Beurteiler hatten zudem als behandelnde Zahnärzte in der klinischen Kohorte fungiert. Daher konnte, obwohl die Behandlungsinformationen dokumentiert waren, nicht sichergestellt werden, dass während der ersten Bildbeurteilung eine vollständige Maskierung der Behandlungszuordnung gewährleistet war.
Vor der eigentlichen Bildinterpretation wurde eine Bewertereichung anhand von 50 Nicht-Studienfällen durchgeführt. Die formale Bewertung begann erst, nachdem die gewichteten κ-Koeffizienten sowohl für die ICDAS-Klassifikation als auch für die studienspezifische radiographische Erscheinungskategorie 0,80 erreicht hatten. Im Kalibrierungssatz lagen die endgültigen gewichteten κ-Koeffizienten bei 0,84 für die ICDAS-Klassifikation und bei 0,88 für die studienspezifische radiographische Erscheinungskategorie.
Wenn während der formalen Beurteilung Uneinigkeiten auftraten, führten die beiden Beurteiler zunächst eine gemeinsame Neubeurteilung durch. Falls kein Konsens erzielt werden konnte, wurde der Fall von einem dritten, leitenden Zahnarzt (J.L.) entschieden, der nicht an der ursprünglichen paargleichen Beurteilung teilgenommen hatte. Eine Neubeurteilung durch einen dritten Gutachter war bei 58 von 1.117 Basisbeurteilungen (5,19 %) und bei 46 von 989 Nachuntersuchungen (4,65 %) erforderlich.
Fälle mit fehlenden Angaben zur primären Exposition oder zum primären Endpunkt wurden von der Analyse ausgeschlossen. Fehlende Daten betrafen lediglich die Gesamtwerte des dmft+DMFT-Scores und des sichtbaren Plaque-Indexes und wurden mittels multipler Imputation durch verkettete Gleichungen mit 20 Imputationen behandelt. Das Imputationsmodell umfasste die Behandlungsgruppe, das Fortschreiten der Läsion nach 12 Monaten, das Alter, das Geschlecht, den Gesamtwert des dmft+DMFT-Scores, den sichtbaren Plaque-Index, den Zahnkieferbogen, den Molartyp, den Durchbruchstatus, die ICDAS-Klassifikation zum Studienbeginn, die studienspezifische radiologische Erscheinungskategorie zum Studienbeginn, den Operateur und die Isolationsmethode.
Statistische Analyse
Es wurde keine a priori-Stichprobengrößenberechnung durchgeführt. Stattdessen wurden alle während des Studienzeitraums identifizierten aufeinanderfolgenden in Frage kommenden Fälle in die Analyse eingeschlossen. Die betroffene Zahnzelle diente als Analyseeinheit, und alle statistischen Modelle berücksichtigten die Clusterung auf Patientenebene, um die Korrelation zwischen mehreren betroffenen Zähnen innerhalb desselben Individuums zu berücksichtigen.
Kontinuierliche Variablen wurden zunächst hinsichtlich ihrer Verteilung überprüft. Normalverteilte Variablen werden als Mittelwert ± Standardabweichung angegeben, nicht normalverteilte Variablen hingegen als Median und Interquartilsabstand. Kategoriale Variablen werden als Anzahlen und Prozentanteile dargestellt. Die Gleichgewichtsbeurteilung der basislinienbezogenen Kovariaten erfolgte mittels standardisierter Mittelwertdifferenzen (SMD), wobei eine absolute SMD <0,10 auf eine ausreichende Balance hinweist. Zur Bestimmung der Kovariatenaufnahme in die Regressionsmodelle wurden keine Hypothesentests zur Basislinie verwendet.
Die vorab festgelegte Nullhypothese besagte, dass das 12-Monats-Risiko für Läsionsprogression zwischen den beiden Behandlungsgruppen nicht unterschiedlich sein würde, was einem relativen Risiko (RR) von 1,00 entspricht. Die primäre Analyse schätzte zunächst die kumulative Inzidenz und das rohe RR der 12-Monats-Läsionsprogression in den beiden Gruppen sowie 95 %-Konfidenzintervalle (CIs). Unadjustierte Analysen wurden mithilfe einer modifizierten Poisson-Regression mit generalisierten Schätzgleichungen (GEEs) durchgeführt, wobei die Behandlungsgruppe als einzige unabhängige Variable berücksichtigt wurde. Die adjustierte Analyse verwendete denselben Modellierungsansatz und fügte zusätzlich Alter, Geschlecht, Gesamtwert des dmft+DMFT-Scores, Sichtbarkeitsindex für Plaque, Molartyp, Kieferbogen, Durchbruchstatus, ICDAS-Klassifikation zum Studienbeginn, radiologische Erscheinungskategorie zum Studienbeginn nach studienspezifischer Definition, Operateur und Isolationsmethode hinzu. Zur Schätzung der robusten Varianz wurde eine austauschbare Arbeitskorrelationsstruktur verwendet, und adjustierte RRs mit entsprechenden 95 %-Konfidenzintervallen wurden angegeben. Kovariaten wurden a priori basierend auf klinischer Relevanz und zeitlicher Abfolge und nicht anhand univariable p-Werte ausgewählt.
Für die primäre Exposition oder das primäre Ergebnis wurde keine Imputation durchgeführt. Fehlende Kovariatendaten wurden mithilfe der multiplen Imputation durch verkettete Gleichungen behandelt, wobei 20 imputierte Datensätze erzeugt wurden; die Schätzungen wurden gemäß Rubins Regeln zusammengefasst. Vier vorab festgelegte Sensitivitätsanalysen wurden durchgeführt, um die Robustheit der primären Ergebnisse zu überprüfen. Erstens wurde die primäre Analyse mittels einer Complete-Case-Analyse wiederholt. Zweitens wurde die primäre Analyse mittels inverser Wahrscheinlichkeitsgewichtung (IPTW) basierend auf Propensity-Scores wiederholt. Das Propensity-Score-Modell enthielt dieselben basislinienbezogenen Kovariaten wie die primäre adjustierte Analyse, und die Kovariatenbalance nach der Gewichtung wurde erneut anhand standardisierter Mittelwertdifferenzen bewertet. Drittens wurde die Analyse wiederholt, nachdem die Studienpopulation auf Zähne mit ausschließlich ICDAS-Baselinestatus 2-Läsionen eingeschränkt worden war. Viertens wurde die Analyse wiederholt, nachdem teilweise durchgebrochene Molaren ausgeschlossen worden waren.
Alle statistischen Tests waren zweiseitig, und ein p-Wert < 0,05 galt als statistisch signifikant. Die statistischen Analysen wurden mit der Software R, Version 4.4.2, durchgeführt.