Iliac vein morphology and wall shear stress: a statistical shape modelling and CFD analysis of patient-specific geometries

Deze studie toont aan dat de keuze voor de mate van anatomische detaillering (van 2D-projecties tot volledige 3D-reconstructies) de statistische vormvariatie en de voorspellingen van wandschuifspanning in de iliacale aderen sterk beïnvloedt, waarbij vereenvoudigde modellen de gebieden met lage schuifspanning systematisch overschatten.

De kern

De Diepe Aders: Waarom de Vorm van je Ader Belangrijker is dan je denkt

Stel je voor dat je bloedstroom door je aderen niet als water in een rechte tuinslang stroomt, maar meer als een rivier die door een landschap met bochten, heuvels en smalle doorgangen stroomt. In dit onderzoek kijken wetenschappers naar de heupaderen (de grote aders in je bekken) en proberen ze te begrijpen waarom sommige mensen sneller bloedstolsels krijgen dan anderen.

Het probleem is diepe veneuze trombose (DVT). Dit is een gevaarlijke bloedstolsel dat zich vormt in de aderen. Wetenschappers weten dat stolsels vaak ontstaan waar het bloed 'traag' stroomt of waar de wand van de ader een lage 'schuifkracht' (wall shear stress) heeft. Maar hoe ziet die ader eruit? En hoe beïnvloedt de vorm precies de stroming?

De onderzoekers hebben een slimme manier bedacht om dit uit te zoeken, door twee dingen te combineren: 3D-schetsen van aderen en computersimulaties van waterstroming.

1. De Drie Manieren om een Ader te Kijken

Om te zien hoe belangrijk de details zijn, hebben ze dezelfde patiëntdata op drie verschillende manieren bekeken, alsof je een foto van een berg bekijkt:

• Manier A (De Volledige 3D-Scopie): Dit is als een gedetailleerde 3D-film van de berg. Je ziet elke steen, elke bocht en elke kuiltje. Dit is de meest nauwkeurige weergave van de echte ader van de patiënt (gemaakt van MRI/CT-scans).
• Manier B (De 2D-Schets): Dit is als een platte tekening van de berg, alsof je er tegenop kijkt. Je ziet de vorm, maar je mist de diepte. Dit is wat artsen vaak zien op een standaard röntgenfoto.
• Manier C (De Vereenvoudigde 3D-Blok): Dit is alsof je de platte tekening (Manier B) uitrekt tot een 3D-blok, maar dan zonder de echte krommingen. Het is een 'idealiserende' versie, alsof je de berg maakt van perfect gladde plastic buizen.

2. De Simulatie: Het Water in de Ader

Vervolgens lieten ze een computer (CFD) zien hoe bloed door deze drie versies stroomt. Ze keken specifiek naar gebieden waar het bloed heel traag stroomt (lage schuifkracht). Waarom? Omdat daar stolsels het makkelijkst ontstaan.

Het verrassende resultaat:
De 'idealiserende' versie (Manier C) gaf een veel grootere kans op stolsels dan de echte, gedetailleerde versie (Manier A).

• De Analogie: Stel je voor dat je een rivier simuleert. Als je de rivier tekent als een perfect rechte, gladde betonnen kanaal (vereenvoudigd model), denk je dat het water overal even traag stroomt. Maar in de echte natuur (de echte 3D-ader) zijn er kleine rotsjes en bochten die het water juist weer in beweging houden. De vereenvoudigde modellen waren dus te pessimistisch: ze dachten dat er meer gevaar was dan er echt was.

3. De Vorm-analyse: Het 'Statistische Vormmodel'

De onderzoekers wilden ook weten: Welke specifieke vormtrekjes zorgen voor deze trage stroming?
Ze gebruikten een slimme techniek (Statistical Shape Modelling) om de variatie tussen de 12 patiënten te meten.

• Bij de 2D-tekeningen: Ze zagen één grote, dominante vormtrek. Als de ader in de tekening een bepaalde bocht had, veranderde de stroming enorm. Het was alsof er één grote 'knop' was die de stroming regelde.
• Bij de echte 3D-modellen: Het was veel ingewikkelder. Er was geen enkele 'knop'. De variatie zat verspreid over heel veel kleine details (kleine krommingen, lichte dikteverschillen). De stroming werd beïnvloed door een combinatie van heel veel kleine factoren, niet door één grote vorm.

Waarom is dit belangrijk?

Dit onderzoek leert ons twee belangrijke dingen voor de toekomst:

1. Kies je model wijs: Als artsen in de toekomst computersimulaties gebruiken om te voorspellen of een patiënt een stolsel krijgt, moeten ze oppassen met te simpele modellen. Een vereenvoudigde 3D-versie kan je laten denken dat een patiënt een groot risico heeft, terwijl de echte, complexe vorm van de ader dat risico juist verkleint.
2. Details tellen: De relatie tussen de vorm van een ader en de bloedstroom is niet simpel. In 3D is het een complex dansje van veel kleine bewegingen, terwijl 2D-schetsen soms een te simpel verhaal vertellen.

Conclusie in één zin

Net zoals je niet kunt voorspellen hoe een rivier stroomt door alleen naar een platte kaart te kijken, kun je het risico op bloedstolsels niet goed inschatten zonder de volledige, complexe 3D-vorm van de ader te begrijpen; te simpele modellen kunnen je namelijk laten denken dat er meer gevaar is dan er echt is.

De onderzoekers hopen dat deze methode in de toekomst helpt om behandelingen (zoals het plaatsen van een stent) persoonlijker en veiliger te maken voor elke individuele patiënt.

Technische samenvatting

Probleemstelling

Diepe veneuze trombose (DVT) is een veelvoorkomende vasculaire aandoening waarbij bloedstolsels zich vormen in de diepe aderen van het been. Hoewel anatomische variaties en stromingsverstoringen bijdragen aan het risico op trombose, zijn de mechanistische koppelingen tussen de vorm van het vat en de hemodynamiek (bloedstroom) slecht gekwantificeerd.

• Huidige beperking: Behandelingen worden voornamelijk gebaseerd op anatomische beelden (MRI/CT), maar stromingskenmerken zoals lage wandschuifspanning (Wall Shear Stress - WSS), die geassocieerd worden met bloedstolling, worden zelden kwantitatief beoordeeld.
• Uitdaging: Het is onduidelijk in hoeverre de keuze voor de geometrische representatie (bijv. vereenvoudigde 2D-projecties versus volledige 3D-reconstructies) de voorspellingen van CFD-simulaties (Computational Fluid Dynamics) beïnvloedt. Bestaande studies gebruiken vaak vereenvoudigde modellen, maar de impact van anatomische precisie op de relatie tussen vorm en stroming is niet systematisch onderzocht.

Methodologie

De auteurs hebben een geïntegreerde werkstroom ontwikkeld die Statistische Vormmodelleren (SSM) combineert met CFD-simulaties om de relatie tussen veneuze geometrie en hemodynamische risico's te analyseren.

1. Dataverzameling:

   • Gebruik van patiëntspecifieke MRI- en CT-scans van 12 gevallen van de gemeenschappelijke iliacale aderen (common iliac veins).
   • Drie niveaus van geometrische fideliteit werden geconstrueerd:
     • (A) Volledige 3D-reconstructies: Direct afgeleid van volumetrische MRI/CT-data.
     • (B) 2D-projecties: Projecties van de 3D-geometrieën op het vlak van een standaard angiogram.
     • (C) 3D-extrusies: Vereenvoudigde 3D-modellen gebaseerd op de 2D-projecties (zonder uit-vlak-kromming), simulerend een workflow die alleen 2D-angiografie gebruikt.

2. Statistische Vormmodelleren (SSM):

   • Gebruik van de bibliotheek Deformetrica gebaseerd op Large Deformation Diffeomorphic Metric Mapping (LDDMM).
   • Deterministische Atlas: Een populatie-template werd gegenereerd om individuele vormen te aligneren.
   • Principal Geodesic Analysis (PGA): Identificeerde de dominante modi van vormvariatie.
   • Alignement-strategie:
     • 2D: Anatomisch gebaseerde alignement (behoud van in-flow hoek en kromming).
     • 3D: Centraal-lijn gebaseerde alignement (verwijdering van globale pose om numerieke stabiliteit te garanderen).

3. CFD-simulaties:

   • Uitgevoerd in ANSYS Fluent 2024 R2.
   • Standaardisatie: Stationaire stroming, laminair, Newtonse vloeistof (bloed), constante uitlaatdruk en parabolische inlaatprofiel.
   • Metriek: Het oppervlak van de vaatwand met lage wandschuifspanning (Low-WSS), gemeten onder drie drempels: ≤0.05, 0.10 en 0.15 Pa.
   • Geautomatiseerde workflow via PyFluent op een HPC-cluster.

4. Statistische Analyse:

   • Berekening van Spearman-rangcorrelaties tussen de SSM-modi (vormvariabele) en de Low-WSS-burden.
   • Toepassing van permutatietests (omdat n=12 klein is) en Benjamini-Hochberg FDR-correctie voor meervoudige toetsing.

Belangrijkste Bijdragen

• Vergelijking van Fideliteit: Het is de eerste studie die systematisch vergelijkt hoe verschillende niveaus van anatomische representatie (2D, 3D-extrusie, volledige 3D) de statistische structuur van vormvariatie en de afgeleide hemodynamische risico's beïnvloeden.
• Integratie SSM-CFD: Een robuust framework dat SSM-modi koppelt aan CFD-metrieken om te bepalen welke vormvariaties het meest relevant zijn voor trombose-risico.
• Kritische Analyse van Alignement: Het toont aan dat de keuze van de alignement-strategie (anatomisch vs. centraal-lijn gebaseerd) fundamenteel de latent ruimte en de interpretatie van vorm-stromingsrelaties verandert.

Resultaten

1. Invloed van Geometrische Fideliteit op CFD:

   • Vereenvoudigde geometrieën (3D-extrusies gebaseerd op 2D) produceerden consistent grotere gebieden met lage WSS dan de volledige 3D-reconstructies.
   • De gemiddelde toename in Low-WSS-oppervlak varieerde tussen 118% en 136% afhankelijk van de drempelwaarde. Dit suggereert dat vereenvoudigde modellen het risico op trombose systematisch overschatten.

2. Statistische Structuur van Vormvariatie:

   • 2D SSM: Toonde een sterk hiërarchisch spectrum. De eerste twee modi verklaarden respectievelijk 33% en 21% van de variantie. Deze modi correleerden sterk met globale eigenschappen zoals buiging en afname in diameter.
   • 3D SSM: Toonde een veel vlakker spectrum. De eerste modi verklaarde slechts 17,4% van de variantie, met de resterende modi die elk ongeveer 7-9% verklaarden. De alignement-strategie (centraal-lijn) had globale pose-variatie verwijderd, waardoor alleen subtiele, lokale morfologische verschillen overbleven.

3. Correlatie tussen Vorm en Stroming:

   • 2D SSM: Er was een sterke, significante negatieve correlatie tussen Modus 1 (globale buiging en diameterverandering) en de Low-WSS-burden. Dit betekent dat specifieke globale vormveranderingen in 2D sterk voorspellen waar lage stroming optreedt.
   • 3D SSM: Correlaties waren zwakker, heterogeen en verspreid over meerdere modi. Geen enkele univariate associatie overleefde de FDR-correctie. Dit komt omdat de vormvariatie in 3D over veel modi is verdeeld, waardoor de invloed van één enkele modus minder dominant is in een kleine steekproef.

4. Sensitiviteit:

   • De rangorde van patiënten op basis van hun Low-WSS-risico bleef stabiel bij variaties in de inlaat-snelheid (±20%), wat aangeeft dat de geometrie de primaire drijver is van hemodynamische variatie, niet de randvoorwaarden binnen een fysiologisch bereik.

Betekenis en Conclusie

De studie concludeert dat de keuze voor model fideliteit en alignement-strategie de interpretatie van hemodynamische risico's fundamenteel verandert:

• Voorzichtigheid bij vereenvoudiging: Het gebruik van vereenvoudigde 3D-extrusies of 2D-projecties leidt tot een systematische overschatting van gebieden met lage wandschuifspanning, wat kan leiden tot valse positieve risico-inschattingen voor trombose.
• Rol van Alignement: De manier waarop vormen worden uitgelijnd (behoud van anatomische oriëntatie vs. verwijdering daarvan) bepaalt welke vormvariabelen als "dominant" worden geïdentificeerd en hoe sterk deze correleren met stromingsparameters.
• Toekomstperspectief: Hoewel de huidige steekproefkleur (n=12) multivariate analyse beperkt, biedt deze studie een reproduceerbaar workflow voor het koppelen van vormbiomarkers aan hemodynamische risico's. Voor klinische toepassing zijn grotere cohorten en geharmoniseerde alignement-strategieën nodig om robuuste, patiëntspecifieke risicomodellen te ontwikkelen.

Kortom, de paper waarschuwt dat in DVT-onderzoek de keuze van het model niet slechts een technische detail is, maar een kritische factor die de conclusies over trombose-risico kan verdraaien.