Een succesvolle synchronisatie wordt aangetoond door statistisch significante circadiane ritmiek van de geanalyseerde klokgenen onder de controleconditie (p < 0,05). In deze studie werden drie onafhankelijke biologische experimenten uitgevoerd, waarbij 0,1% dimethylsulfoxide diende als negatieve controle. Bmal1 en Per2 werden geselecteerd als representatieve circadiane markers omdat zij kerncomponenten zijn van de moleculaire circadiane klok en een goed gekarakteriseerd antiphasisch expressiepatroon vertonen. Representatieve ritmische expressieprofielen verkregen na synchronisatie worden getoond in Figuur 2. De overeenkomstige circadiane parameters, inclusief ritmiek, MESOR, amplitude en piektijd, zijn samengevat in Tabel 1, terwijl de paarsgewijze vergelijking tussen Bmal1 en Per2 wordt gepresenteerd in Tabel 2.

Figuur 2. Circadiane oscillaties van Bmal1 en Per2 na synchronisatie door serumshock. Representatieve circadiane expressieprofielen van Bmal1 (donkerblauw) en Per2 (lichtblauw) in gesynchroniseerde mHypoE-42-cellen, gemeten met kwantitatieve real-time polymerasekettenreactie (qRT-PCR) van 24 tot 48 u na synchronisatie. De curven representeren de gefitte circadiane modellen die zijn gegenereerd met CircaCompare.Klik hier om een grotere versie van deze figuur te bekijken.
| Gen | Ritmiek
(p-waarde) | MESOR | Amplitude | Piekperiode
(h) |
| Bmal1 | 0.0205 | −0.2948 | 0.4648 | 20.4239 |
| Per2 | 0.0071 | −0.318 | 0.4629 | 9.4883 |
| Bmal1 + Compound A | 0.1068 | — | — | — |
Tabel 1: Geschatte circadiane ritmeparameters voor Bmal1 en Per2 onder basale omstandigheden en voor Bmal1 na behandeling met Compound A. De circadiane ritmiek werd geëvalueerd met behulp van CircaCompare. De tabel vermeldt de p-waarde voor ritmiek, de midline estimating statistic of rhythm (MESOR), de amplitude en de piektijd voor elke conditie. Circadiane parameters worden niet vermeld voor Bmal1 na behandeling met Compound A omdat de ritmiek statistisch niet significant was (p = 0,1068).
Schatting van het verschil in MESOR
(p-waarde) | Schatting van het amplitudoverschil
(p-waarde) | Schatting van het faseverschil
(p-waarde) | Gedeelde periode
(h) |
| −0.0232 (p = 0.8904) | −0.0019 (p = 0.9936) | −10.9356 (p = 2.3155 × 10⁻6) | 24 |
Tabel 2: Paarsgewijze vergelijking van circadiane parameters tussen Bmal1 en Per2.Paarsgewijze vergelijking van de midline estimating statistic of rhythm (MESOR), amplitude en fase tussen de gefitte circadiane expressieprofielen van Bmal1 en Per2, geschat met behulp van CircaCompare. Waarden worden gepresenteerd als parameterverschillen met de bijbehorende p-waarden. De gemeenschappelijke periode werd tijdens het fitten van het model vastgelegd op 24 h.
Onder basale condities vertoonden zowel Bmal1 als Per2 een statistisch significante ritmiciteit (Bmal1: p = 0,0205; Per2: p = 0,0071). CircaCompare-analyse toonde verder aan dat de twee genen in tegenovergestelde fasen oscilleerden, met een faseverschil van −10,9356 h (p = 2,3155 × 10⁻6), wat consistent is met hun verwachte antiphasische relatie binnen het moleculaire circadiane kloknetwerk. De RNA-kwaliteitsmetrieken die ten grondslag liggen aan de daaropvolgende genexpressieanalyse zijn opgenomen in Aanvullende Tabel 1.
Het protocol kan ook worden toegepast om de circadiane ritmiek tussen experimentele condities te vergelijken, waaronder gezonde modellen versus ziektegevoelige modellen of controlegroepen versus behandelgroepen. Als representatief voorbeeld werd een opkomende milieuverontreiniging en endocriene disruptor, die eerder in de menselijke hypothalamus is aangetroffen, geëvalueerd; deze wordt hierna aangeduid als Compound A. Representatieve resultaten worden getoond in Figuur 3, en de overeenkomstige rhythmiciteitsanalyse is samengevat in Tabel 1.

Figuur 3. Effect van Compound A op de circadiane ritmiek van Bmal1-expressie. Representatief expressieprofiel van Bmal1 in gesynchroniseerde mHypoE-42-cellen na behandeling met Compound A, gemeten via kwantitatieve real-time polymerasekettingreactie (qRT-PCR). De curven representeren de gefitte circadiane modellen gegenereerd met CircaCompare. Klik hier om een grotere versie van deze figuur te bekijken.
Na blootstelling aan Compound A vertoonde Bmal1 geen statistisch significante rhythmiciteit meer (p = 0,1068). Onder deze omstandigheden voldeed het gefitte circadiaanse model niet aan de statistische criteria voor rhythmiciteit. Bijgevolg worden circadiaanse parameters zoals MESOR, amplitude en piektijd niet gerapporteerd, omdat deze niet als biologisch interpreteerbaar worden beschouwd wanneer de rhythmiciteit niet statistisch significant is. Dit voorbeeld illustreert hoe de workflow onderscheid kan maken tussen rhythmische en niet-rhythmische expressieprofielen en kan worden gebruikt om potentiële wijzigingen in de circadiaanse genexpressie onder verschillende experimentele condities te beoordelen. Om potentiële technische variabiliteit in verband met de verschoofde bemonsteringsstrategie te minimaliseren, werden alle cultuurplaten bereid uit dezelfde celsuspensie, onder identieke cultuurcondities gehouden en gesynchroniseerd met behulp van dezelfde experimentele workflow, waarbij alleen de timing van de synchronisatie verschilde.
Aanvullende figuur 1. Wekelijks schema van het protocol voor versprongen synchronisatie. Representatieve wekelijkse tijdlijn die het uitplaten van cellen, serumstarvatie, serumshock-synchronisatie, synchronisatiestop en circadiane monstername voor de twee versprongen experimentele sets illustreert (Set 1, paars; Set 2, groen). Het versprongen ontwerp maakt de collectie van circadiane tijdspunten (CT) CT0–CT48 tijdens standaard werktijden mogelijk, waardoor monstername gedurende de nacht wordt vermeden. Klik hier om dit bestand te downloaden.
Aanvullende tabel 1. RNA-kwaliteitsmetrieken voor alle biologische monsters opgenomen in het circadiane tijdreeks-experiment. De RNA-zuiverheid werd beoordeeld aan de hand van de absorbentieverhoudingen A260/A280 and A260/A230, en de RNA-concentratie werd spectrofotometrisch bepaald vóór de synthese van complementair DNA (cDNA). Monsters zijn geïdentificeerd op basis van verzameltijd, behandelingsgroep en biologisch replicaat. Deze metingen werden gebruikt om de RNA-kwaliteit te verifiëren voorafgaand aan de daaropvolgende kwantitatieve real-time polymeraseketenreactie (qRT-PCR) analyse. Klik hier om dit bestand te downloaden.
Aanvullende Tabel 2. Primersequenties en amplificatiecondities gebruikt voor kwantitatieve real-time polymerasekettingreactie (qRT-PCR). Forward- en revers oligonucleotide primersequenties worden gepresenteerd in de 5′-3′ oriëntatie, samen met het bijbehorende NCBI-toegangsnummer(s), de annealingtemperatuur en de primerconcentratie gebruikt voor de amplificatie van elk doelgen. Klik hier om dit bestand te downloaden.
Aanvullende tabel 3. Voorbeeld van een input-dataset voor circadiane ritmeanalyse met CircaCompare. Voorbeeld van een dataset met kommagescheiden waarden (CSV), geformatteerd voor analyse met CircaCompare. De dataset bevat drie vereiste variabelen: Time (h), Group en Outcome (ΔΔCq), overeenkomstig het formaat dat wordt beschreven in het protocol voor data-analyse en dat is gebruikt als voorbeeld-inputbestand voor de schatting van circadiane parameters. Klik hier om dit bestand te downloaden.
Aanvullend bestand 1. Voorbeeld R-script voor circadiane ritmeanalyse met CircaCompare. Voorbeeld R-script voor het analyseren van circadiane genexpressiedata met het CircaCompare-pakket (versie 0.2.0) in R (versie 4.5.1). Het script importeert de inputdataset, voert de schatting van circadiane parameters uit met een periode van 24 u, exporteert de statistische samenvatting en genereert een plot van publicatiekwaliteit. Vervang vóór gebruik de tijdelijke bestandsnamen (XXX.csv, YYY.csv en ZZZ.svg) door de gewenste input- en outputbestandsnamen en pas de plotparameters (bijv. ylim) aan naar wens voor de dataset. Klik hier om dit bestand te downloaden.