Distrofia LGMD R2: 2 traiectorii de progresie prezise prin inteligență artificială

Context și relevanță clinică

Distrofia musculară de centură de tip R2, legată de deficitul de disferlină (LGMD R2, cunoscută anterior ca LGMD 2B), este o boală musculară progresivă, ereditară, determinată de variante patogene în gena DYSF. Disferlina este o proteină esențială pentru repararea membranei fibrei musculare; absența ei duce la degenerare musculară lentă, dar inexorabilă. Pacienții își pierd treptat forța la nivelul centurilor scapulară și pelvină, ajung să meargă greu, apoi își pierd capacitatea de a se ridica de pe podea, de a urca scări și, în final, de a mai merge independent.

Marea problemă clinică nu este atât diagnosticul, cât imposibilitatea de a prezice cum va evolua un pacient anume. Variabilitatea cursului bolii este izbitoare: doi pacienți cu aceeași mutație și aceeași vârstă pot avea ritmuri de declin complet diferite. Această heterogenitate face dificilă consilierea pacienților, planificarea îngrijirii și, mai ales, proiectarea studiilor clinice. Atunci când testezi un medicament nou într-o boală rară, cu progresie lentă și inegală, ai nevoie să știi din start cine declină repede și cine declină încet — altfel efectul real al tratamentului se pierde în zgomotul variabilității naturale.

Evaluarea clinică standard se bazează pe scale funcționale și pe teste de forță, dar acestea surprind târziu modificările și au o sensibilitate limitată la schimbări fine. Rezonanța magnetică cantitativă (RMN cantitativ) a apărut ca un instrument complementar capabil să detecteze înlocuirea grăsoasă a mușchiului mai devreme și mai precis decât examenul clinic. Problema este că integrarea unor surse de date atât de diferite — măsurători clinice, teste cronometrate, valori de laborator și hărți imagistice de fracție grăsoasă — depășește metodele statistice convenționale. Aici intervine inteligența artificială, capabilă să combine modalități heterogene de date și să extragă tipare predictive pe care analiza clasică le ratează.

Design și metodologie

Studiul a folosit date provenind din Cohorta de Studii Clinice 1 (COS 1) dedicată disferlinopatiilor, o colecție longitudinală internațională de pacienți urmăriți sistematic. Au fost analizați 188 de participanți. Progresia bolii a fost cuantificată cu ajutorul scalei North Star Assessment for Limb Girdle type Muscular Dystrophies (NSAD), un instrument funcțional validat specific pentru distrofiile de centură. Au fost incluși doar pacienții ambulatori, cu un scor NSAD maxim de cel puțin 20, pentru a evita efectele de plafon și planșeu care distorsionează măsurarea progresiei.

Subgrupul folosit pentru modelarea traiectoriilor a cuprins 101 pacienți care aveau cel puțin 5 vizite consecutive, urmăriți în medie 1.280,8 zile, adică aproximativ trei ani și jumătate. Dintre aceștia, 58 dispuneau și de RMN cantitativ al coapsei și gambei. Obiectivul primar a fost evoluția scorului NSAD pe parcursul a trei ani.

• Selecția traiectoriilor: clustering ierarhic cu distanțe euclidiene și metoda Ward, aplicat pacienților cu serii temporale suficient de lungi.
• Selecția variabilelor: dintr-un set inițial de 307 variabile, după excluderi (funcția pulmonară, capacitatea vitală forțată, variantele genetice individuale, scorul EK) au rămas 126 de variabile; măsurătorile bilaterale au fost mediate, iar mușchii au fost grupați pe compartimente anatomice.
• Optimizarea trăsăturilor: o procedură de Recursive Feature Optimization a izolat cele mai discriminante 30 de trăsături; un clasificator Random Forest a fost antrenat prin validare încrucișată stratificată.
• Importanța trăsăturilor: cadrul SHAP a fost folosit pentru a determina contribuția fiecărei variabile la predicție.
• Construcția modelelor: regresie logistică pas cu pas, cu selecție înainte, a dus la două modele — unul exclusiv clinic și unul combinat clinic plus RMN.

Compararea grupurilor s-a făcut cu testul Mann-Whitney U, cu corecție Benjamini-Hochberg pentru testări multiple, pentru a controla rata erorilor de tip fals pozitiv inevitabile când se compară zeci de variabile simultan.

Populația studiată

La momentul inițial, cohorta avea o vârstă medie de 35,2 ani (deviație standard 12,6), cu vârsta medie de debut a simptomelor la 23,1 ani (deviație standard 10,2) și o durată medie a bolii de 12,1 ani (deviație standard 7,7). Scorul NSAD mediu de start a fost de 35 de puncte (deviație standard 12), confirmând o populație ambulatorie cu boală activă, dar nu în stadiu terminal. Repartiția pe sexe a fost echilibrată — 53 de femei, reprezentând 52,5% din lot. Originea etnică predominantă a fost cea caucaziană, în proporție de 71,3%.

Criteriile de includere au cerut confirmarea genetică sau biochimică a bolii: fie două variante patogene în gena DYSF, fie o singură variantă patogenă însoțită de absența demonstrată a expresiei disferlinei la nivel tisular. Astfel, lotul a fost omogen din punct de vedere diagnostic, eliminând cazurile incerte care ar fi putut introduce zgomot în model. Condiția de ambulație și pragul NSAD de cel puțin 20 au asigurat că toți participanții se aflau într-o fereastră a bolii în care progresia funcțională poate fi măsurată semnificativ, fereastra cea mai relevantă pentru studiile clinice de intervenție.

Această selecție atentă reflectă o realitate a cercetării în bolile rare: numărul mic de pacienți disponibili obligă la maximizarea valorii fiecărui caz. Cohorta COS 1 reprezintă una dintre cele mai mari colecții longitudinale de pacienți cu disferlinopatie din lume, ceea ce conferă rezultatelor o greutate deosebită într-un domeniu unde majoritatea studiilor lucrează cu câteva zeci de pacienți.

Intervenție și comparator

Studiul nu a fost unul de tratament, ci de modelare prognostică. „Intervenția" a constat în construirea și validarea a două modele predictive concurente, comparate direct pentru capacitatea de a clasifica pacienții în progresie rapidă versus progresie moderată.

• Modelul exclusiv clinic: bazat numai pe date accesibile la patul pacientului — durata bolii, creatinfosfokinaza (CK) și viteza la testul de mers/alergare pe 10 metri (10MWR).
• Modelul combinat clinic plus RMN: adaugă biomarkeri imagistici cantitativi — fracția de grăsime și aria secțiunii transversale contractile în mușchii coapsei anterioare, vastul intermediar și gracilis.

Logica comparației este eminamente practică: modelul clinic se poate aplica oriunde, fără echipament costisitor, în timp ce modelul combinat necesită RMN cantitativ, disponibil doar în centre specializate. Întrebarea pe care studiul o adresează direct este dacă informația imagistică suplimentară justifică efortul și costul, aducând un câștig real de acuratețe predictivă.

Rezultate principale

Cele două traiectorii de progresie

Clusteringul ierarhic a separat clar populația în două stadii de progresie. Stadiul rapid a cuprins 30 de pacienți, cu o pierdere medie de 14,4 puncte NSAD în trei ani (deviație standard 4,7). Stadiul moderat a inclus 71 de pacienți, cu o pierdere medie de doar 3,8 puncte NSAD în aceeași perioadă (deviație standard 2,8). Diferența de aproape patru ori în ritmul de declin funcțional demonstrează că „progresia LGMD R2 poate fi grupată în traiectorii clinice distincte", nu este un continuum uniform.

Profilul pacienților din stadiul rapid a fost paradoxal la prima vedere, dar coerent cu istoria naturală a bolii: aceștia erau mai tineri (vârstă mediană 29 ani față de 35 ani, p<0,001) și aveau o durată a bolii mai scurtă (8 ani față de 14 ani, p<0,001). Cu alte cuvinte, declinul cel mai accelerat survine devreme în evoluția bolii, înainte ca mușchiul să fie deja înlocuit masiv de grăsime. Markerii de leziune musculară activă au fost mult mai mari în acest grup: CK de 6.630 U/L față de 3.486 U/L (p<0,001), reflectând o degenerare musculară în plină desfășurare.

Forța musculară a fost, de asemenea, mai afectată funcțional în grupul rapid: forța de extensie a șoldului (testare musculară manuală) de 3,8 față de 6,0 (p<0,001), forța de prehensiune de 81,0 livre față de 90,0 livre (p<0,01) și timpul de ridicare de pe podea de 5,7 secunde față de 4,6 secunde (p<0,01). Aparenta contradicție — pacienți mai tineri, dar mai slabi — se explică prin faptul că ei se află în faza de declin cel mai abrupt.

Biomarkerii imagistici cantitativi

RMN cantitativ a dezvăluit un tipar revelator: pacienții cu progresie rapidă aveau, paradoxal, mușchi relativ mai bine păstrați, cu fracție de grăsime mai mică, tocmai pentru că boala nu apucase încă să distrugă țesutul. Fracția grăsoasă în vastul intermediar a fost de 17,8% în grupul rapid față de 39,5% în grupul moderat (p<0,001) — diferența cea mai marcantă. Tipare similare au apărut la gracilis (7,9% față de 11,3%, p=0,048), sartorius (10,8% față de 13,5%, p=0,044) și vastul medial (22,3% față de 38,2%, p=0,036).

Interpretarea este esențială pentru practică: un mușchi încă „plin", neînlocuit de grăsime, dar la un pacient tânăr cu CK ridicată, este un mușchi aflat la începutul unei prăbușiri rapide. Conservarea relativă a vastului intermediar și a gracilisului devine astfel un semnal de alarmă, nu de liniștire.

Performanța modelelor predictive

Modelul combinat clinic plus RMN a atins o acuratețe echilibrată de 83,7%, cu sensibilitate de 78,7% și specificitate de 88,7%, o arie sub curbă (AUC) de 0,91 (interval de încredere 95%: 0,86–0,95) și un R² Nagelkerke de 0,587, explicând aproape 59% din varianța progresiei. Predictorii reținuți au fost fracția grăsoasă în vastul intermediar, aria contractilă a gracilisului, fracția grăsoasă și aria contractilă în coapsa anterioară, plus durata bolii.

Modelul exclusiv clinic a atins o acuratețe echilibrată de 78,7%, cu sensibilitate de 74,3%, specificitate de 83,1%, AUC de 0,84 (interval de încredere 95%: 0,79–0,89) și R² Nagelkerke de 0,443. Cei trei predictori clinici au fost CK (cu un raport de șanse de 1,42 per 1.000 U/L creștere, p<0,001), durata bolii (raport de șanse 0,88, p<0,001) și viteza la testul de 10 metri (raport de șanse 0,18, p<0,001).

Analize secundare și de sensibilitate

Concordanța dintre cele două modele a fost de 77,9%. Printre predicțiile concordante, 71,0% au fost corecte și doar 6,9% incorecte, ceea ce arată o robustețe ridicată atunci când ambele modele „sunt de acord". Mai interesant, atunci când modelele divergeau, modelul combinat avea dreptate în 15,2% din cazuri, față de doar 6,9% pentru modelul clinic — dovada directă că biomarkerii imagistici aduc informație suplimentară reală, nu redundantă, exact în cazurile dificile, ambigue.

Analiza importanței trăsăturilor prin cadrul SHAP a confirmat că durata bolii a fost singura variabilă temporală reținută, surprinzând faptul că stadiul bolii, mai degrabă decât vârsta cronologică, guvernează ritmul de progresie. CK a captat informație despre stadiul bolii dincolo de simplele măsuri de severitate clinică, comportându-se ca un marker dinamic al degenerării active. Implicarea heterogenă a mușchilor coapsei anterioare s-a corelat cu tiparele de progresie rapidă, reflectând binecunoscutul semn al „diamantului" cvadricepsului — o bombare asimetrică a coapsei anterioare care trădează înlocuirea grăsoasă neuniformă.

Robustețea modelelor a fost asigurată prin validare încrucișată stratificată, esențială într-un eșantion mic unde supraînvățarea este un risc constant. Cercetătorii au dezvoltat și un instrument bazat pe web pentru a face modelul accesibil în practica clinică, transformând un exercițiu statistic complex într-un calculator de risc utilizabil la patul pacientului.

Discuție și limitări

Studiul demonstrează că tehnicile analitice bazate pe inteligență artificială permit o analiză mai profundă și mai sistematică a seturilor de date neuromusculare complexe decât metodele convenționale. Rezultatele întăresc cunoștințe deja stabilite — pacienții cu progresie mai rapidă sunt mai tineri și în stadii mai precoce ale bolii, așa cum se raportase anterior — dar identifică și tipare noi. CK captează informație despre stadiul bolii care depășește măsurile de severitate clinică, iar implicarea heterogenă a mușchiului coapsei anterioare se corelează cu tiparul de progresie rapidă.

Mesajul clinic central este contraintuitiv, dar puternic: tinerețea și mușchiul relativ păstrat nu sunt semne liniștitoare la un pacient cu LGMD R2 și CK ridicată, ci, dimpotrivă, semnalează fereastra de declin cel mai abrupt. Aceasta răstoarnă intuiția clasică și are implicații directe asupra momentului în care ar trebui intervenit terapeutic — cât mai devreme, înainte ca degenerarea rapidă să consume rezerva musculară.

Limitările sunt cele inerente bolilor rare. Eșantionul, deși mare pentru o disferlinopatie, rămâne modest în termeni absoluți, iar subgrupul cu RMN cantitativ a fost de numai 58 de pacienți, ceea ce limitează puterea statistică a modelului combinat. Caracterul observațional al cohortei nu permite concluzii de cauzalitate. Disponibilitatea RMN cantitativ doar în centre specializate restrânge aplicabilitatea imediată a modelului combinat. În plus, validarea s-a făcut intern, prin validare încrucișată, dar o validare externă pe o cohortă independentă rămâne necesară înainte de adoptarea pe scară largă. Generalizarea la pacienții non-ambulatori sau la cei cu boală foarte avansată nu este posibilă, întrucât aceștia au fost excluși prin design.

Concluzii și implicații practice

Această cercetare stabilește existența a două traiectorii distincte de progresie în LGMD R2, identificabile printr-o analiză integrată asistată de inteligență artificială. Modelul combinat, care încorporează biomarkeri de rezonanță magnetică cantitativă, oferă o acuratețe prognostică superioară pentru predicția progresiei pe trei ani, cu aplicații practice în stratificarea pacienților pentru studiile clinice și în consilierea individuală.

Pentru clinician, mesajul operațional este dublu. Întâi, un pacient tânăr, cu durată scurtă a bolii și CK marcat crescută, trebuie considerat la risc înalt de declin rapid, indiferent de cât de bine arată mușchiul la examenul clinic. Al doilea, acolo unde RMN cantitativ este disponibil, conservarea relativă a vastului intermediar la un astfel de pacient confirmă, nu infirmă, riscul de progresie accelerată. Pentru cercetarea terapeutică, capacitatea de a separa din start pacienții cu declin rapid de cei cu declin lent înseamnă studii clinice mai mici, mai scurte și mai sensibile la efectul real al unui tratament — un câștig esențial într-o boală rară unde fiecare pacient contează.