[go: nahoru, domu]

Jump to content

Tractography: Difference between revisions

From Wikipedia, the free encyclopedia
Content deleted Content added
m Removed duplicate link
OAbot (talk | contribs)
m Open access bot: doi updated in citation with #oabot.
(12 intermediate revisions by 11 users not shown)
Line 1: Line 1:
{{short description|3D visualization of nerve tracts via diffusion MRI}}


{{Infobox diagnostic
{{Infobox diagnostic
Line 22: Line 23:
| LOINC = <!--{{LOINC|code}}-->
| LOINC = <!--{{LOINC|code}}-->
}}
}}
In [[neuroscience]], '''tractography''' is a [[3D modeling]] technique used to visually represent [[nerve tract]]s using data collected by [[diffusion MRI]].<ref name=":0" /> It uses special techniques of [[magnetic resonance imaging]] (MRI) and computer-based diffusion MRI. The results are presented in two- and three-dimensional images called '''tractograms'''.
In [[neuroscience]], '''tractography''' is a [[3D modeling]] technique used to visually represent [[nerve tract]]s using data collected by [[diffusion MRI]].<ref name=":0" /> It uses special techniques of [[magnetic resonance imaging]] (MRI) and computer-based diffusion MRI. The results are presented in two- and three-dimensional images called '''tractograms'''.<ref>{{cite book |doi=10.1093/med/9780199541164.001.0001 |title=Atlas of Human Brain Connections |date=2012 |last1=Catani |first1=Marco |last2=Thiebaut De Schotten |first2=Michel |isbn=978-0-19-954116-4 }}{{pn|date=September 2023}}</ref>


In addition to the long tracts that connect the [[brain]] to the rest of the body, there are complicated [[neural circuit]]s formed by short connections among different [[Cerebral cortex|cortical]] and [[subcortical]] regions. The existence of these tracts and circuits has been revealed by [[histochemistry]] and [[biology|biological]] techniques on [[post-mortem]] specimens. Nerve tracts are not identifiable by direct exam, [[computed tomography|CT]], or [[MRI]] scans. This difficulty explains the paucity of their description in [[neuroanatomy]] atlases and the poor understanding of their functions.
In addition to the long tracts that connect the [[brain]] to the rest of the body, there are complicated [[neural circuit]]s formed by short connections among different [[Cerebral cortex|cortical]] and [[subcortical]] regions. The existence of these tracts and circuits has been revealed by [[histochemistry]] and [[biology|biological]] techniques on [[post-mortem]] specimens. Nerve tracts are not identifiable by direct exam, [[computed tomography|CT]], or [[MRI]] scans. This difficulty explains the paucity of their description in [[neuroanatomy]] atlases and the poor understanding of their functions.


The most advanced tractography algorithm can produce 90% of the ground truth bundles, but it still contains a substantial amount of invalid results.<ref name="ChallengeNatComm2017">{{cite journal|display-authors=6|vauthors=Maier-Hein KH, Neher PF, Houde JC, Côté MA, Garyfallidis E, Zhong J, Chamberland M, Yeh FC, Lin YC, Ji Q, Reddick WE, Glass JO, Chen DQ, Feng Y, Gao C, Wu Y, Ma J, Renjie H, Li Q, Westin CF, Deslauriers-Gauthier S, González JO, Paquette M, St-Jean S, Girard G, Rheault F, Sidhu J, Tax CM, Guo F, Mesri HY, Dávid S, Froeling M, Heemskerk AM, Leemans A, Boré A, Pinsard B, Bedetti C, Desrosiers M, Brambati S, Doyon J, Sarica A, Vasta R, Cerasa A, Quattrone A, Yeatman J, Khan AR, Hodges W, Alexander S, Romascano D, Barakovic M, Auría A, Esteban O, Lemkaddem A, Thiran JP, Cetingul HE, Odry BL, Mailhe B, Nadar MS, Pizzagalli F, Prasad G, Villalon-Reina JE, Galvis J, Thompson PM, Requejo FS, Laguna PL, Lacerda LM, Barrett R, Dell'Acqua F, Catani M, Petit L, Caruyer E, Daducci A, Dyrby TB, Holland-Letz T, Hilgetag CC, Stieltjes B, Descoteaux M|date=November 2017|title=The challenge of mapping the human connectome based on diffusion tractography|journal=Nature Communications|volume=8|issue=1|pages=1349|doi=10.1038/s41467-017-01285-x|pmc=5677006|pmid=29116093}}</ref>
The most advanced tractography algorithm can produce 90% of the ground truth bundles, but it still contains a substantial amount of invalid results.<ref name="ChallengeNatComm2017">{{cite journal|display-authors=6|vauthors=Maier-Hein KH, Neher PF, Houde JC, Côté MA, Garyfallidis E, Zhong J, Chamberland M, Yeh FC, Lin YC, Ji Q, Reddick WE, Glass JO, Chen DQ, Feng Y, Gao C, Wu Y, Ma J, Renjie H, Li Q, Westin CF, Deslauriers-Gauthier S, González JO, Paquette M, St-Jean S, Girard G, Rheault F, Sidhu J, Tax CM, Guo F, Mesri HY, Dávid S, Froeling M, Heemskerk AM, Leemans A, Boré A, Pinsard B, Bedetti C, Desrosiers M, Brambati S, Doyon J, Sarica A, Vasta R, Cerasa A, Quattrone A, Yeatman J, Khan AR, Hodges W, Alexander S, Romascano D, Barakovic M, Auría A, Esteban O, Lemkaddem A, Thiran JP, Cetingul HE, Odry BL, Mailhe B, Nadar MS, Pizzagalli F, Prasad G, Villalon-Reina JE, Galvis J, Thompson PM, Requejo FS, Laguna PL, Lacerda LM, Barrett R, Dell'Acqua F, Catani M, Petit L, Caruyer E, Daducci A, Dyrby TB, Holland-Letz T, Hilgetag CC, Stieltjes B, Descoteaux M|date=November 2017|title=The challenge of mapping the human connectome based on diffusion tractography|journal=Nature Communications|volume=8|issue=1|pages=1349|doi=10.1038/s41467-017-01285-x|pmc=5677006|pmid=29116093|bibcode=2017NatCo...8.1349M }}</ref>


== MRI technique ==
== MRI technique ==
Line 36: Line 37:
Tractography is performed using data from [[diffusion MRI]]. The free water diffusion is termed "[[isotropic]]" diffusion. If the water diffuses in a medium with barriers, the diffusion will be uneven, which is termed [[anisotropic]] diffusion. In such a case, the relative mobility of the [[molecules]] from the origin has a shape different from a [[sphere]]. This shape is often modeled as an [[ellipsoid]], and the technique is then called [[diffusion tensor imaging]]. Barriers can be many things: cell membranes, axons, myelin, etc.; but in [[white matter]] the principal barrier is the [[myelin]] sheath of [[axons]]. Bundles of axons provide a barrier to perpendicular diffusion and a path for parallel diffusion along the orientation of the fibers.
Tractography is performed using data from [[diffusion MRI]]. The free water diffusion is termed "[[isotropic]]" diffusion. If the water diffuses in a medium with barriers, the diffusion will be uneven, which is termed [[anisotropic]] diffusion. In such a case, the relative mobility of the [[molecules]] from the origin has a shape different from a [[sphere]]. This shape is often modeled as an [[ellipsoid]], and the technique is then called [[diffusion tensor imaging]]. Barriers can be many things: cell membranes, axons, myelin, etc.; but in [[white matter]] the principal barrier is the [[myelin]] sheath of [[axons]]. Bundles of axons provide a barrier to perpendicular diffusion and a path for parallel diffusion along the orientation of the fibers.


Anisotropic diffusion is expected to be increased in areas of high mature axonal order. Conditions where the [[myelin]] or the structure of the axon are disrupted, such as [[Physical trauma|trauma]]<ref>{{cite journal |last1=Wade |first1=Ryckie G. |last2=Tanner |first2=Steven F. |last3=Teh |first3=Irvin |last4=Ridgway |first4=John P. |last5=Shelley |first5=David |last6=Chaka |first6=Brian |last7=Rankine |first7=James J. |last8=Andersson |first8=Gustav |last9=Wiberg |first9=Mikael |last10=Bourke |first10=Grainne |title=Diffusion Tensor Imaging for Diagnosing Root Avulsions in Traumatic Adult Brachial Plexus Injuries: A Proof-of-Concept Study |journal=Frontiers in Surgery |date=16 April 2020 |volume=7 |doi=10.3389/fsurg.2020.00019|doi-access=free }}</ref>, [[tumors]], and [[inflammation]] reduce anisotropy, as the barriers are affected by destruction or disorganization.
Anisotropic diffusion is expected to be increased in areas of high mature axonal order. Conditions where the [[myelin]] or the structure of the axon are disrupted, such as [[Physical trauma|trauma]],<ref>{{cite journal |last1=Wade |first1=Ryckie G. |last2=Tanner |first2=Steven F. |last3=Teh |first3=Irvin |last4=Ridgway |first4=John P. |last5=Shelley |first5=David |last6=Chaka |first6=Brian |last7=Rankine |first7=James J. |last8=Andersson |first8=Gustav |last9=Wiberg |first9=Mikael |last10=Bourke |first10=Grainne |title=Diffusion Tensor Imaging for Diagnosing Root Avulsions in Traumatic Adult Brachial Plexus Injuries: A Proof-of-Concept Study |journal=Frontiers in Surgery |date=16 April 2020 |volume=7 |page=19 |doi=10.3389/fsurg.2020.00019|pmid=32373625 |pmc=7177010|doi-access=free }}</ref> [[tumors]], and [[inflammation]] reduce anisotropy, as the barriers are affected by destruction or disorganization.


Anisotropy is measured in several ways. One way is by a ratio called [[fractional anisotropy]] (FA). An FA of 0 corresponds to a perfect sphere, whereas 1 is an ideal linear diffusion. Few regions have FA larger than 0.90. The number gives information about how aspherical the diffusion is but says nothing of the direction.
Anisotropy is measured in several ways. One way is by a ratio called [[fractional anisotropy]] (FA). An FA of 0 corresponds to a perfect sphere, whereas 1 is an ideal linear diffusion. Few regions have FA larger than 0.90. The number gives information about how aspherical the diffusion is but says nothing of the direction.
Line 51: Line 52:
== Mathematics ==
== Mathematics ==


Using [[Diffusion MRI|diffusion tensor MRI]], one can measure the [[Diffusion MRI#ADC|apparent diffusion coefficient]] at each [[voxel]] in the image, and after [[Multi-linear regression|multilinear regression]] across multiple images, the whole diffusion tensor can be reconstructed.<ref name=":0">{{cite journal | vauthors = Basser PJ, Pajevic S, Pierpaoli C, Duda J, Aldroubi A | title = In vivo fiber tractography using DT-MRI data | journal = Magnetic Resonance in Medicine | volume = 44 | issue = 4 | pages = 625–32 | date = October 2000 | pmid = 11025519 | doi = 10.1002/1522-2594(200010)44:4<625::AID-MRM17>3.0.CO;2-O }}</ref>
Using [[Diffusion MRI|diffusion tensor MRI]], one can measure the [[Diffusion MRI#ADC|apparent diffusion coefficient]] at each [[voxel]] in the image, and after [[Multi-linear regression|multilinear regression]] across multiple images, the whole diffusion tensor can be reconstructed.<ref name=":0">{{cite journal | vauthors = Basser PJ, Pajevic S, Pierpaoli C, Duda J, Aldroubi A | title = In vivo fiber tractography using DT-MRI data | journal = Magnetic Resonance in Medicine | volume = 44 | issue = 4 | pages = 625–32 | date = October 2000 | pmid = 11025519 | doi = 10.1002/1522-2594(200010)44:4<625::AID-MRM17>3.0.CO;2-O | doi-access = }}</ref>


Suppose there is a fiber tract of interest in the sample. Following the [[Frenet–Serret formulas]], we can formulate the space-path of the fiber tract as a parameterized curve:
Suppose there is a fiber tract of interest in the sample. Following the [[Frenet–Serret formulas]], we can formulate the space-path of the fiber tract as a parameterized curve:
Line 57: Line 58:
:<math> \frac{d\mathbf{r}(s)}{ds} = \mathbf{T}(s), </math>
:<math> \frac{d\mathbf{r}(s)}{ds} = \mathbf{T}(s), </math>


where <math>\mathbf{T}(s)</math> is the tangent vector of the curve. The reconstructed diffusion tensor <math> D </math> can be treated as a matrix, and we can easily compute its [[eigenvalues]] <math> \lambda_1, \lambda_2, \lambda_3 </math> and [[eigenvectors]] <math> \mathbf{u}_1, \mathbf{u}_2, \mathbf{u}_3 </math>. By equating the eigenvector corresponding to the largest eigenvalue with the direction of the curve:
where <math>\mathbf{T}(s)</math> is the tangent vector of the curve. The reconstructed diffusion tensor <math> D </math> can be treated as a matrix, and we can compute its [[eigenvalues]] <math> \lambda_1, \lambda_2, \lambda_3 </math> and [[eigenvectors]] <math> \mathbf{u}_1, \mathbf{u}_2, \mathbf{u}_3 </math>. By equating the eigenvector corresponding to the largest eigenvalue with the direction of the curve:


:<math> \frac{d\mathbf{r}(s)}{ds} = \mathbf{u}_1(\mathbf{r}(s)) </math>
:<math> \frac{d\mathbf{r}(s)}{ds} = \mathbf{u}_1(\mathbf{r}(s)) </math>

Revision as of 18:34, 2 December 2023

Tractography
Tractography of human brain
Purposeused to visually represent nerve tracts

In neuroscience, tractography is a 3D modeling technique used to visually represent nerve tracts using data collected by diffusion MRI.[1] It uses special techniques of magnetic resonance imaging (MRI) and computer-based diffusion MRI. The results are presented in two- and three-dimensional images called tractograms.[2]

In addition to the long tracts that connect the brain to the rest of the body, there are complicated neural circuits formed by short connections among different cortical and subcortical regions. The existence of these tracts and circuits has been revealed by histochemistry and biological techniques on post-mortem specimens. Nerve tracts are not identifiable by direct exam, CT, or MRI scans. This difficulty explains the paucity of their description in neuroanatomy atlases and the poor understanding of their functions.

The most advanced tractography algorithm can produce 90% of the ground truth bundles, but it still contains a substantial amount of invalid results.[3]

MRI technique

DTI of the brachial plexus - see https://doi.org/10.3389/fsurg.2020.00019 for more information
Tractographic reconstruction of neural connections by diffusion tensor imaging (DTI)
MRI tractography of the human subthalamic nucleus

Tractography is performed using data from diffusion MRI. The free water diffusion is termed "isotropic" diffusion. If the water diffuses in a medium with barriers, the diffusion will be uneven, which is termed anisotropic diffusion. In such a case, the relative mobility of the molecules from the origin has a shape different from a sphere. This shape is often modeled as an ellipsoid, and the technique is then called diffusion tensor imaging. Barriers can be many things: cell membranes, axons, myelin, etc.; but in white matter the principal barrier is the myelin sheath of axons. Bundles of axons provide a barrier to perpendicular diffusion and a path for parallel diffusion along the orientation of the fibers.

Anisotropic diffusion is expected to be increased in areas of high mature axonal order. Conditions where the myelin or the structure of the axon are disrupted, such as trauma,[4] tumors, and inflammation reduce anisotropy, as the barriers are affected by destruction or disorganization.

Anisotropy is measured in several ways. One way is by a ratio called fractional anisotropy (FA). An FA of 0 corresponds to a perfect sphere, whereas 1 is an ideal linear diffusion. Few regions have FA larger than 0.90. The number gives information about how aspherical the diffusion is but says nothing of the direction.

Each anisotropy is linked to an orientation of the predominant axis (predominant direction of the diffusion). Post-processing programs are able to extract this directional information.

This additional information is difficult to represent on 2D grey-scaled images. To overcome this problem, a color code is introduced. Basic colors can tell the observer how the fibers are oriented in a 3D coordinate system, this is termed an "anisotropic map". The software could encode the colors in this way:

  • Red indicates directions in the X axis: right to left or left to right.
  • Green indicates directions in the Y axis: posterior to anterior or from anterior to posterior.
  • Blue indicates directions in the Z axis: foot-to-head direction or vice versa.

The technique is unable to discriminate the "positive" or "negative" direction in the same axis.

Mathematics

Using diffusion tensor MRI, one can measure the apparent diffusion coefficient at each voxel in the image, and after multilinear regression across multiple images, the whole diffusion tensor can be reconstructed.[1]

Suppose there is a fiber tract of interest in the sample. Following the Frenet–Serret formulas, we can formulate the space-path of the fiber tract as a parameterized curve:

where is the tangent vector of the curve. The reconstructed diffusion tensor can be treated as a matrix, and we can compute its eigenvalues and eigenvectors . By equating the eigenvector corresponding to the largest eigenvalue with the direction of the curve:

we can solve for given the data for . This can be done using numerical integration, e.g., using Runge–Kutta, and by interpolating the principal eigenvectors.

See also

References

  1. ^ a b Basser PJ, Pajevic S, Pierpaoli C, Duda J, Aldroubi A (October 2000). "In vivo fiber tractography using DT-MRI data". Magnetic Resonance in Medicine. 44 (4): 625–32. doi:10.1002/1522-2594(200010)44:4<625::AID-MRM17>3.0.CO;2-O. PMID 11025519.
  2. ^ Catani, Marco; Thiebaut De Schotten, Michel (2012). Atlas of Human Brain Connections. doi:10.1093/med/9780199541164.001.0001. ISBN 978-0-19-954116-4.[page needed]
  3. ^ Maier-Hein KH, Neher PF, Houde JC, Côté MA, Garyfallidis E, Zhong J, et al. (November 2017). "The challenge of mapping the human connectome based on diffusion tractography". Nature Communications. 8 (1): 1349. Bibcode:2017NatCo...8.1349M. doi:10.1038/s41467-017-01285-x. PMC 5677006. PMID 29116093.
  4. ^ Wade, Ryckie G.; Tanner, Steven F.; Teh, Irvin; Ridgway, John P.; Shelley, David; Chaka, Brian; Rankine, James J.; Andersson, Gustav; Wiberg, Mikael; Bourke, Grainne (16 April 2020). "Diffusion Tensor Imaging for Diagnosing Root Avulsions in Traumatic Adult Brachial Plexus Injuries: A Proof-of-Concept Study". Frontiers in Surgery. 7: 19. doi:10.3389/fsurg.2020.00019. PMC 7177010. PMID 32373625.