Simplifying Daily Cortisol Cycle Analysis: Validation and Benchmarking of the Cortisol Sine Score Against Cosinor and JTK_CYCLE models

Dit artikel introduceert en valideert de Cortisol Sine Score (CSS), een nieuwe, modelvrije scalair maatstaf die de dagelijkse cortisolcyclus effectief samenvat tot één interpreteerbare waarde, waarmee complexe circadiane analyses worden vereenvoudigd voor gebruik in klinische workflows, machine learning en grote epidemiologische studies.

De kern

Het Cortisol Sine Score: Een Simpele "Ritme-Check" voor je Lichaam

Stel je voor dat je lichaam een orkest is en je hormonen de muzikanten. Een van de belangrijkste muzikanten is cortisol, het stresshormoon. Op een gezonde dag speelt deze muzikant een heel duidelijk ritme: hij begint de ochtend met een krachtige, hoge noot (wanneer je wakker wordt) en zakt gedurende de dag langzaam af tot een zachte, lage noot in de avond. Dit is je natuurlijke "dagritme".

Maar wat als dit ritme verstoord raakt? Misschien speelt de muzikant te laat, te vroeg, of speelt hij de hele dag even hard (een platte lijn zonder pieken en dalen). Dit gebeurt vaak bij mensen met veel stress, slaapproblemen of bepaalde ziektes.

Het Probleem: Te Muzieknoten om te tellen
Vroeger was het heel moeilijk om dit ritme te meten. Wetenschappers moesten de hele dag door cortisol meten (bijvoorbeeld om 2:00, 6:00, 10:00 uur, etc.) en toen met ingewikkelde wiskundige formules proberen te berekenen of het ritme goed was. Het was alsof je een heel orkeststuk moet analyseren door elke noot apart op te schrijven. Dit was te lastig voor artsen in de spreekkamer of voor grote computerprogramma's (kunstmatige intelligentie) die snel een oordeel moeten vellen.

De Oplossing: De "Cortisol Sine Score" (CSS)
De auteurs van dit onderzoek hebben een slimme, nieuwe manier bedacht: de Cortisol Sine Score.

Stel je voor dat je in plaats van alle noten apart te schrijven, gewoon kijkt naar de hoofdrichting van de muziek:

• Een positief getal (+): De muziek begint hoog en zakt af. Dit is een gezond ritme (je bent 's ochtends wakker en 's avonds rustig).
• Een negatief getal (-): De muziek begint laag en wordt 's avonds juist harder. Dit is verward (je bent 's avonds nog piekerig en 's ochtends moe).
• Een getal dicht bij 0: De muziek is de hele dag even hard. Geen piek, geen dal. Dit is een "plat" ritme, vaak een teken van chronische stress of uitputting.

Deze score is als een thermometer voor je dagritme. In plaats van een heel ingewikkeld rapport te lezen, krijg je één simpel getal dat je direct vertelt: "Is je ritme gezond, verward of plat?"

Waarom is dit zo slim?

1. Minder gedoe voor de proefpersonen:
   Normaal moet je 6 keer per dag een speekselmonster geven (ook 's nachts om 2 uur!). Dat is oncomfortabel. De onderzoekers ontdekten dat je met deze nieuwe score al heel nauwkeurig werkt met slechts 3 of 4 metingen per dag. Je hoeft niet meer wakker te worden om 2 uur 's nachts; je kunt gewoon meten bij het opstaan, 's middags en 's avonds. Het is alsof je de temperatuur van een kamer goed kunt schatten met één meting in plaats van tien.

2. Werkt in de computerwereld:
   Omdat het maar één getal is, kunnen computers (AI) dit heel makkelijk gebruiken. De onderzoekers hebben getoond dat dit getal zelfs kan voorspellen hoe de darmbacteriën van iemand eruitzien! Het is alsof je door één simpel getal te kijken, de gezondheid van een heel complex ecosysteem in je buik kunt begrijpen.

3. Koppeling aan stress:
   Ze hebben bewezen dat mensen met meer sociale problemen (zoals lage inkomens of moeilijke woonomstandigheden) vaker een "plat" of verward ritme hebben (een score dicht bij 0 of negatief). De nieuwe score pakt dit probleem heel goed op, beter dan de oude, ingewikkelde methoden.

Conclusie
Deze nieuwe "Cortisol Sine Score" is als een simpel dashboard in een auto. Vroeger moest je als bestuurder (of arts) naar tientallen aparte meters kijken om te weten of de motor goed liep. Nu heb je één lampje dat je direct vertelt: "Alles goed" (groen/positief), "Let op" (geel/0) of "Probleem" (rood/negatief).

Het maakt het meten van stress en ritme makkelijker, sneller en toegankelijker voor iedereen, van de huisarts tot de grote wetenschappelijke studies.

Technische samenvatting

1. Het Probleem

De dagelijkse cortisolcyclus is een cruciale indicator voor de functie van de hypothalamus-hypofyse-bijnias (HPA-as). Veranderingen in dit ritme zijn geassocieerd met chronische stress, psychiatrische aandoeningen, metabole disfunctie en verstoorde slaap-waakcycli.

Huidige geavanceerde analysemethoden, zoals Cosinor-regressie (parametrisch) en JTK_CYCLE (niet-parametrisch), genereren echter meerdere uitvoerparameters (amplitude, acrofasie, mesor, p-waarden voor ritmiciteit). Deze multi-output resultaten zijn moeilijk te integreren in eenvoudige statistische modellen, klinische workflows en Machine Learning/AI-pijplijnen die vaak een enkele numerieke invoer vereisen. Er is een kritiek tekort aan een unificerende, scalair meetwaarde die de timing, uitlijning en richting van de cortisolcyclus samenvat zonder de interpretatiebaarheid of toepasbaarheid te verliezen.

2. Methodologie

De auteurs ontwikkelden en valideerden een nieuwe, modelvrije, scalair meetwaarde: de Cortisol Sine Score (CSS).

• Wiskundige Basis: De CSS is geïnspireerd op Fourier-decompositie. Het is een gewogen som van cortisolmetingen over de dag, waarbij de gewichten worden afgeleid van de sinus van het tijdstip van de dag (omgezet in radialen op een 24-uurscirkel).
  • Formule: $CSS = \sum_{i=1}^{n} (Cortisol_i \times \sin(\frac{2\pi t_i}{24}))$
  • Interpretatie:
    • Positieve waarden: Duiden op een ochtend-dominant patroon (gezond circadiaans ritme).
    • Negatieve waarden: Duiden op een vertraagd of omgekeerd ritme (avond-dominant).
    • Waarden dicht bij nul: Duiden op een afgevlakt ritme (chronische stress of circadiaanse disruptie).
• Validatie-Data: De studie gebruikte data van 501 zwangere vrouwen uit de March of Dimes cohortstudie, met in totaal 3.006 speekselmonsters (6 metingen per dag: 02:00, 06:00, 10:00, 14:00, 18:00, 22:00).
• Validatie-aanpak:
  • Vergelijking: De CSS werd vergeleken met waargenomen en geschatte amplitude en acrofasie via Cosinor-regressie en JTK_CYCLE.
  • Statistische Tests: Lineaire gemengde modellen (LMM) en Generalized Additive Models (GAM) met cyclische splines werden gebruikt om correlaties te testen.
  • Benchmarking: De prestaties van de CSS werden vergeleken met andere metrics (AUCg, MESOR, atypiciteit) in associaties met Sociaal Nadeel (Social Disadvantage) en Darmmicrobioom-samenstelling (via Random Forest regressie).
  • Optimalisatie: Er werden 63 mogelijke combinaties van steekproeftijdstippen geanalyseerd om het minimale aantal benodigde metingen te bepalen.

3. Belangrijkste Bijdragen

• Ontwikkeling van de CSS: Een nieuwe, modelvrije scalair meetwaarde die complexe tijdsreeksdata reduceert tot één interpreteerbare waarde.
• Open Source Software: Het ontwikkelen van het R-pakket CortSineScore om de berekening toegankelijk te maken.
• Protocol-Optimalisatie: Identificatie van de minimale steekproefgrootte die nodig is voor betrouwbare CSS-berekening, wat de last voor deelnemers aanzienlijk vermindert.
• Toepassing in Multi-omics: Demonstratie van de CSS als effectieve invoer voor ML/AI-modellen in complexe studies (bijv. microbioom), waar traditionele methoden faalden.

4. Resultaten

• Validatie: De CSS toonde sterke correlaties met amplitude en acrofasie zoals geschat door Cosinor en JTK_CYCLE.
  • Classificatieprestatie (AUC) voor het onderscheiden van hoge vs. lage amplitude was uitstekend (AUC = 0,89).
  • Er was een significante associatie tussen CSS en de timing van de piek (acrofasie), met een verschuiving van ongeveer 5 uur tussen de hoogste en laagste CSS-tertiles.
• Minimale Steekproefgrootte:
  • De CSS kan betrouwbaar worden geschat met slechts 3 tot 4 metingen per dag, in plaats van de volledige 6.
  • Een 4-metingen protocol (06:00, 10:00, 14:00, 18:00) bereikte een correlatie van r = 0,953 en een MAE van 0,023 ten opzichte van het volledige protocol.
  • Een 3-metingen protocol (02:00, 06:00, 18:00) leverde ook hoge nauwkeurigheid op (r = 0,952).
  • De 06:00-meting (ochtendpiek) werd geïdentificeerd als essentieel; zonder deze meting daalde de nauwkeurigheid drastisch.
• Benchmarking:
  • Sociaal Nadeel: De CSS toonde de sterkste associatie met sociaal nadeel (Spearman's ρ = -0,37), beter dan of gelijk aan geavanceerde metrics zoals Cosinor-amplitude.
  • Darmmicrobioom: In een Random Forest-model voorspelde de CSS de samenstelling van het darmmicrobioom het best, terwijl traditionele amplitude-metrics minder goed presteerden. Dit suggereert dat de CSS (die richting en timing combineert) biologisch relevanter is voor circadiaus-gevoelige systemen dan amplitude alleen.

5. Betekenis en Conclusie

De Cortisol Sine Score (CSS) overbrugt de kloof tussen gedetailleerde circadiaanse biologie en grote schaal datawetenschap.

• Praktisch: Het biedt een snelle, interpreteerbare screeningtool voor HPA-as-functie in klinische en epidemiologische studies.
• Efficiëntie: Door de mogelijkheid om met 3-4 metingen te werken, wordt de last voor deelnemers (vooral bij kwetsbare groepen zoals zwangere vrouwen of kinderen) aanzienlijk verminderd.
• Toekomst: De eenvoudige wiskundige aard van de CSS maakt het ideaal voor implementatie in draagbare biosensoren en real-time monitoring, zonder zware rekenkracht te vereisen.

De CSS vervangt niet de gedetailleerde analysemethoden (zoals Cosinor) wanneer volledige parameterisatie nodig is, maar fungeert als een krachtige, complementaire tool die de integratie van cortisoldata in moderne statistische en AI-frameworks mogelijk maakt.