DXA-Derived Skeletal Phenotypes and Hip Fracture Risk: A Backdoor-Adjusted Causal Analysis

Résumé technique : Phénotypes squelettiques dérivés de la DXA et risque de fracture de la hanche

Énoncé du problème
La fracture de la hanche représente un fardeau de santé important avec des taux de mortalité et de morbidité élevés. L'évaluation clinique actuelle repose largement sur l'outil d'évaluation du risque de fracture (FRAX), qui intègre des facteurs de risque cliniques et, facultativement, la densité minérale osseuse (DMO) du col du fémur, ainsi que sur les mesures de synthèse standard de l'absorptiométrie à rayons X à double énergie (DXA), telles que la DMO de surface et les T-scores. Cependant, les examens DXA contiennent une multitude de mesures squelettiques spécifiques à chaque région (par exemple, à travers le col du fémur, le fémur total, la région trochantérique, la région de Ward, la diaphyse fémorale et le bassin) qui ne sont pas systématiquement utilisées dans la stratification du risque. De plus, les études observationnelles existantes se concentrent souvent sur la précision de la prédiction ou sur l'inférence causale informée génétiquement (randomisation mendélienne) plutôt que de comparer explicitement plusieurs phénotypes dérivés de la DXA au sein d'une cohorte observationnelle sous un cadre rigoureux d'ajustement des facteurs de confusion. Il existe une lacune critique dans la compréhension de quels phénotypes spécifiques dérivés de la DXA présentent les relations causales les plus fortes avec le risque de fracture de la hanche après ajustement pour les facteurs de confusion, et si le classement de ces phénotypes par taille d'effet causal améliore la stratification du risque en aval.

Méthodologie
L'étude a utilisé les données de 21 098 participants de la UK Biobank avec des dossiers de santé liés et des mesures DXA de la hanche. L'analyse a suivi un cadre causal pré-spécifié :

1. Sélection des phénotypes : Seize phénotypes squelettiques dérivés de la DXA ont été évalués, couvrant trois familles de mesures (Contenu Minéral Osseux [CMO], Densité Minérale Osseuse [DMO] et T-score) à travers six régions anatomiques (col du fémur, fémur total, région trochantérique, région de Ward, diaphyse fémorale et bassin).
2. **Cadre causal : la construction d'un graphe acyclique dirigé (DAG) informé par la matière a guidé la sélection des facteurs de confusion. L'ensemble d'ajustement comprenait des facteurs démographiques/de taille corporelle (âge, sexe, taille, poids), des facteurs socio-économiques (revenu du ménage), des facteurs de mode de vie/nutritionnels (tabagisme, alcool, régime alimentaire), des antécédents cliniques (polyarthrite rhumatoïde, chutes, fractures antérieures) et des indicateurs fonctionnels (vitesse de marche, force de préhension manuelle).
3. Estimation de l'ATE par ajustement de la porte dérobée (Backdoor-Adjusted) : Pour chaque l'un des 16 phénotypes, un modèle de probabilité linéaire spécifique au phénotype a été ajusté pour estimer l'Effet de Traitement Moyen (ATE) ajusté par la porte dérobée sur le risque de fracture de la hanche. Les effets ont été rapportés en tant que différences de risque absolu par écart-type (SD) d'augmentation du phénotype, moyennées sur la distribution des facteurs de confusion pré-spécifiés. La robustesse a été vérifiée via des analyses de placebo d'exposition et de réfutation de cause commune aléatoire.
4. Analyse de l'hétérogénéité : Les Effets de Traitement Moyens Conditionnels (CATE) ont été estimés pour le phénotype de tête (DMO du fémur total) en utilisant des modèles de forêt causale (causal forest) pour évaluer la variation de l'effet entre les sous-groupes (âge, sexe, IMC, etc.).
5. Évaluation prédictive : Une analyse prédictive supplémentaire a évalué si le classement des phénotypes par magnitude de leur ATE ajusté par la porte dérobée améliorait la stratification du risque. Des modèles de régression logistique ont été entraînés en utilisant uniquement des variables cliniques, diverses combinaisons de phénotypes squelettiques classés par ATE, et l'ensemble complet des caractéristiques. La performance a été comparée à FRAX en utilisant l'aire sous la courbe (AUC), la sensibilité et la spécificité.

Résultats clés

• Effets causaux : Les 16 phénotypes dérivés de la DXA présentaient des ATE négatifs ajustés par la porte dérobée, indiquant qu'une valeur de phénotype plus élevée est associée à un risque de fracture de la hanche plus faible. Les effets les plus importants ont été observés pour le CMO du fémur total et la DMO du fémur total, montrant tous deux une différence de risque de −0,0047 par augmentation de 1 SD (soit environ 4,7 fractures de moins pour 1 000 participants par SD). En revanche, le CMO de la région de Ward présentait le plus faible effet (−0,0019).
• Hétérogénéité : L'analyse CATE pour la DMO du fémur total a révélé une hétérogénéité significative. L'effet protecteur était plus fort (CATE plus négatifs) chez les participants plus âgés et ceux ayant un indice de masse corporelle (IMC) plus faible. Le revenu du ménage, l'âge, le poids, la polyarthrite rhumatoïde et la force de préhension manuelle ont été identifiés comme des contributeurs clés de cette hétérogénéité.
• Performance prédictive : La combinaison des variables cliniques avec les 11 meilleurs phénotypes squelettiques classés par ATE a produit la performance prédictive la plus élevée (AUC = 0,842), surpassant FRAX (AUC = 0,709) et un modèle uniquement clinique (AUC = 0,754). Le modèle des 11 meilleures caractéristiques a démontré une sensibilité significativement plus élevée (0,748 contre 0,443) que FRAX tout en maintenant une spécificité similaire.
• Importance des caractéristiques : L'analyse SHAP a indiqué que le modèle prédictif final était principalement piloté par les phénotypes squelettiques dérivés de la DXA (spécifiquement le CMO du fémur total, le CMO de la région trochantérique et diverses mesures de DMO), les variables cliniques apportant des contributions additionnelles modestes.

Signification et affirmations
Les auteurs affirment que cette étude fournit une évaluation systématique des phénotypes squelettiques dérivés de la DXA dans un cadre explicite d'ajustement des facteurs de confusion, distinguant les variables qui améliorent la classification du risque de celles qui présentent des effets causaux plus importants. La principale signification réside dans la démonstration que :

1. Variation au niveau des phénotypes : Les phénotypes squelettiques dérivés de la DXA diffèrent par leurs effets causaux ajustés par la porte dérobée sur le risque de fracture de la hanche ; ils ne contribuent pas uniformément selon les régions ou les familles de mesures.
2. Utilité du classement causal : Le classement des phénotypes par leurs ATE ajustés par la porte dérobée aide à identifier un ensemble parcimonieux de caractéristiques (les 11 meilleures) qui, combinées aux variables cliniques, améliorent la stratification du risque au-delà des outils standards comme FRAX.
3. Aperçu clinique : Bien que les T-scores restent la norme pour la classification de l'ostéoporose, l'étude suggère que les phénotypes spécifiques à une région et à une famille de mesures (particulièrement le CMO et la DMO du fémur total) peuvent contenir des informations de risque non entièrement capturées par les T-scores seuls.
4. Risque individualisé : Les conclusions soutiennent la prise en compte de l'hétérogénéité de l'effet, suggérant que la pertinence des phénotypes squelettiques (par exemple, la DMO du fémur total) varie en fonction des profils des patients tels que l'âge et l'IMC.

L'article maintient un ton modeste concernant la mise en œuvre clinique, notant que les résultats prédictifs sont dérivés d'un échantillon spécifique de la UK Biobank et nécessitent une validation externe avant que le modèle proposé ne puisse remplacer les outils de risque clinique existants. L'étude positionne ses conclusions comme des preuves complémentaires pour affiner l'interprétation de la DXA et les stratégies de stratification du risque.