Ein erfolgreicher Synchronisierungserfolg zeigt sich durch eine statistisch signifikante circadiane Rhythmik der analysierten Clock-Gene unter Kontrollbedingungen (p < 0,05). In dieser Studie wurden drei unabhängige biologische Experimente durchgeführt, wobei 0,1 % Dimethylsulfoxid als negativer Kontrollstoff diente. Bmal1 und Per2 wurden als repräsentative circadiane Marker ausgewählt, da sie zentrale Bestandteile der molekularen circadianen Uhr darstellen und ein gut charakterisiertes antiphasisches Expressionsmuster aufweisen. Repräsentative rhythmische Expressionsprofile nach der Synchronisierung sind in Abbildung 2 dargestellt. Die entsprechenden circadianen Parameter, einschließlich Rhythmizität, MESOR, Amplitude und Zeitpunkt des Peaks, sind in Tabelle 1 zusammengefasst, während der paarweise Vergleich zwischen Bmal1 und Per2 in Tabelle 2 dargestellt ist.

Abbildung 2. Circadiane Oszillationen von Bmal1 und Per2 nach Synchronisierung mittels Serumschock. Repräsentative circadiane Expressionsprofile von Bmal1 (dunkelblau) und Per2 (hellblau) in synchronisierten mHypoE-42-Zellen, gemessen mittels quantitativer Echtzeit-Polymerase-Kettenreaktion (qRT-PCR) von 24 bis 48 h nach der Synchronisierung. Die Kurven zeigen die angepassten circadianen Modelle, die mit CircaCompare erstellt wurden.Klicken Sie bitte hier, um eine vergrößerte Version dieser Abbildung anzusehen.
| Gen | Rhythmizität
(p-Wert) | MESOR | Amplitude | Maximalzeit
(h) |
| Bmal1 | 0.0205 | −0.2948 | 0.4648 | 20.4239 |
| Per2 | 0.0071 | −0.318 | 0.4629 | 9.4883 |
| Bmal1 + Verbindung A | 0.1068 | — | — | — |
Tabelle 1: Unter basalen Bedingungen sowie für Bmal1 nach Behandlung mit Verbindung A geschätzte Parameter der circadianen Rhythmik für Bmal1 und Per2. Die circadiane Rhythmik wurde mittels CircaCompare bewertet. In der Tabelle sind der p-Wert für die Rhythmik, die MESOR (midline estimating statistic of rhythm), die Amplitude und die Zeit des Peaks für jede Bedingung angegeben. Für Bmal1 nach Behandlung mit Verbindung A werden keine Parameter der circadianen Rhythmik angegeben, da die Rhythmik nicht statistisch signifikant war (p = 0,1068).
Schätzung der MESOR-Differenz
(p-Wert) | Schätzung der Amplituden-Differenz
(p-Wert) | Schätzung der Phasen-Differenz
(p-Wert) | Gemeinsame Periode
(h) |
| −0.0232 (p = 0.8904) | −0.0019 (p = 0.9936) | −10.9356 (p = 2.3155 × 10⁻6) | 24 |
Tabelle 2: Pairwise-Vergleich der circadianen Parameter zwischen Bmal1 und Per2. Pairwise-Vergleich der MESOR (midline estimating statistic of rhythm), Amplitude und Phase zwischen den angepassten circadianen Expressionsprofilen von Bmal1 und Per2, wie mittels CircaCompare geschätzt. Die Werte werden als Parameterdifferenzen mit den entsprechenden p-Werten angegeben. Die gemeinsame Periode wurde während der Modellanpassung auf 24 h festgelegt.
Unter basalen Bedingungen zeigten sowohl Bmal1 als auch Per2 eine statistisch signifikante Rhythmik (Bmal1: p = 0,0205; Per2: p = 0,0071). Die CircaCompare-Analyse zeigte zudem, dass die beiden Gene in entgegengesetzten Phasen oszillierten, mit einer Phasendifferenz von −10,9356 h (p = 2,3155 × 10⁻6), was ihrer erwarteten antiphasischen Beziehung im molekularen circadianen Uhrwerk entspricht. Die RNA-Qualitätskennzahlen, die der nachgeschalteten Genexpressionsanalyse zugrunde liegen, sind in Supplementary Table 1 angegeben.
Das Protokoll kann auch verwendet werden, um die circadiane Rhythmik zwischen experimentellen Bedingungen zu vergleichen, einschließlich gesunder versus Krankheitsmodelle oder Kontroll- versus Behandlungsgruppen. Als repräsentatives Beispiel wurde ein neu auftretender Umweltkontaminant und endokriner Disruptor, der zuvor im menschlichen Hypothalamus nachgewiesen wurde, untersucht und wird hier als Verbindung A bezeichnet. Repräsentative Ergebnisse sind in Abbildung 3 dargestellt, und die entsprechende Analyse der Rhythmik ist in Tabelle 1 zusammengefasst.

Abbildung 3. Wirkung von Verbindung A auf die circadiane Rhythmik der Bmal1-Expression. Repräsentatives Expressionsprofil von Bmal1 in synchronisierten mHypoE-42-Zellen nach Behandlung mit Verbindung A, gemessen mittels quantitativer Echtzeit-Polymerase-Kettenreaktion (qRT-PCR). Die Kurven stellen die mithilfe von CircaCompare erstellten angepassten circadianen Modelle dar. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Nach der Exposition gegenüber Verbindung A zeigte Bmal1 keine statistisch signifikante Rhythmizität mehr (p = 0,1068). Unter diesen Bedingungen erfüllte das angepasste zirkadiane Modell die statistischen Kriterien für Rhythmizität nicht. Daher werden zirkadiane Parameter wie MESOR, Amplitude und Zeitpunkt des Peaks nicht angegeben, da sie biologisch nicht interpretierbar gelten, wenn keine statistisch signifikante Rhythmizität vorliegt. Dieses Beispiel verdeutlicht, wie der Arbeitsablauf rhythmische von nicht-rhythmischen Expressionsprofilen unterscheiden und zur Beurteilung möglicher Veränderungen in der zirkadianen Genexpression unter verschiedenen experimentellen Bedingungen herangezogen werden kann. Um potenzielle technische Variabilität im Zusammenhang mit der gestaffelten Probennahmestrategie zu minimieren, wurden alle Kulturplatten aus derselben Zellsuspension hergestellt, unter identischen Kulturbedingungen gehalten und mit demselben experimentellen Arbeitsablauf synchronisiert, wobei lediglich der Zeitpunkt der Synchronisierung variierte.
Abbildung 1, Zusatzmaterial. Wöchentlicher Zeitplan des gestaffelten Synchronisationsprotokolls. Repräsentativer wöchentlicher Ablauf, der das Zellansäen, die Serumverarmung, die Synchronisation durch Serumstoss, die Arrestierung der Synchronisation und die zirkadiane Probennahme für die beiden gestaffelten Experimentalsätze (Satz 1, lila; Satz 2, grün) veranschaulicht. Durch den gestaffelten Versuchsansatz ist es möglich, zirkadiane Zeitpunkte (CT) von CT0 bis CT48 innerhalb der regulären Arbeitszeit zu erfassen und gleichzeitig die Probennahme über Nacht zu vermeiden. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle 1. RNA-Qualitätskennzahlen für alle biologischen Proben, die in das circadiane Zeitverlaufs-Experiment einbezogen wurden. Die RNA-Reinheit wurde anhand der Absorptionsverhältnisse A260/A280 und A260/A230 bewertet, und die RNA-Konzentration wurde vor der Synthese von komplementärer DNA (cDNA) spektrophotometrisch bestimmt. Die Proben sind durch Sammelzeitpunkt, Behandlungsgruppe und biologische Wiederholung gekennzeichnet. Diese Messwerte dienten der Überprüfung der RNA-Qualität vor der nachgeschalteten quantitativen Echtzeit-Polymerase-Kettenreaktion (qRT-PCR). Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzungstabelle 2. Primersequenzen und Amplifikationsbedingungen, die für die quantitative Echtzeit-Polymerase-Kettenreaktion (qRT-PCR) verwendet wurden. Die Vorwärts- und Rückwärts-Oligonukleotid-Primersequenzen sind in 5′-3′-Orientierung angegeben, zusammen mit den entsprechenden NCBI-Zugangsnummer(n), der Annealing-Temperatur und der verwendeten Primerkonzentration für die Amplifikation jedes Zielgens. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzungstabelle 3. Beispiel-Eingabedatensatz zur Analyse von circadianen Rhythmen mit CircaCompare. Beispielhafte durch Kommas getrennte Werte (CSV) im Format für die Analyse mit CircaCompare. Der Datensatz enthält drei erforderliche Variablen: Zeit (h), Gruppe und Ergebnis (ΔΔCq), entsprechend dem im Protokoll zur Datenanalyse beschriebenen Format, das als Beispiel-Eingabedatei für die Schätzung zirkadianer Parameter verwendet wird. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Zusatzdatei 1. Beispiel-R-Skript zur Analyse zirkadianer Rhythmen mit CircaCompare. Beispiel-R-Skript zur Analyse zirkadianer Genexpressionsdaten unter Verwendung des CircaCompare-Pakets (Version 0.2.0) in R (Version 4.5.1). Das Skript importiert den Eingabedatensatz, führt eine Schätzung der zirkadianen Parameter mit einer Periode von 24 Stunden durch, exportiert die statistische Zusammenfassung und erzeugt eine publikationsfähige Abbildung. Vor der Verwendung müssen die Platzhalter-Dateinamen (XXX.csv, YYY.csv und ZZZ.svg) durch die gewünschten Namen für Ein- und Ausgabedateien ersetzt und die Plot-Parameter (z. B. ylim) entsprechend dem Datensatz angepasst werden. Bitte klicken Sie hier, um diese Datei herunterzuladen.