1
|
Liu W, Guo J, Luo J, Ren Q, Chen Z, Qu Z, Wu Z, Ni J, Xu X, Rashid M, Luo J, Yin H, Yang Z, Liu G. Analysis of microRNA expression profiles dynamic in different life stages of Haemaphysalis longicornis ticks by deep sequencing of small RNA libraries. Ticks Tick Borne Dis 2020; 11:101427. [PMID: 32370927 DOI: 10.1016/j.ttbdis.2020.101427] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2019] [Revised: 03/18/2020] [Accepted: 03/23/2020] [Indexed: 10/24/2022]
Abstract
The three-host tick Haemaphysalis longicornis is an obligate blood-sucking ectoparasite. In life-stage transitions, microRNAs (miRNAs) show a variety of expression changes. To investigate these changes, deep sequencing technology was applied to identify the conserved and potentially novel miRNAs expressed during the different life stages of H. longicornis. Total RNA from eggs, unfed larvae, unfed nymphs and unfed adults was extracted for deep sequence analysis. Deep sequencing on a Hiseq 4000 generated a total of 111,192,069 reads, grouped into four small RNA (sRNA) libraries, one for each of the four developmental stages of H. longicornis. Among these sequences, 78 conserved and 55 potentially novel miRNAs were identified, including stage-specific and differentially expressed miRNAs. Gene ontology (GO) analysis indicated significantly enriched GO terms related to cell proliferation and differentiation, including specific terms for the processes of development, growth, metabolism, regulation of biological functions, reproduction, and membrane enzyme regular activity. Kyoto Encyclopedia of Gene and Genomes (KEGG) analysis revealed a significant enrichment of the insulin, notch, Hippo, and Wnt signaling pathways for growth and development. Our data highlight the abundance of miRNA changes (conserved and potentially novel) in the different life stages of H. longicornis. In particular, stage-specific miRNAs, as observed, are essential regulators for the development of H. longicornis.
Collapse
Affiliation(s)
- Wenge Liu
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China; College of Animal Veterinary Medicine, Northwest A & F University, Yangling, China.
| | - Junhui Guo
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Jin Luo
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Qiaoyun Ren
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Ze Chen
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Zhiqiang Qu
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Zegong Wu
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Jun Ni
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Xiaofeng Xu
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Muhammad Rashid
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Jianxun Luo
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| | - Hong Yin
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China; Jiangsu Co-Innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou, China.
| | - Zengqi Yang
- College of Animal Veterinary Medicine, Northwest A & F University, Yangling, China.
| | - Guangyuan Liu
- State Key Laboratory of Veterinary Etiological Biology, Key Laboratory of Veterinary Parasitology of Gansu Province, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Science, Lanzhou, China.
| |
Collapse
|
2
|
Hernandez M. PDE-constrained LDDMM via geodesic shooting and inexact Gauss-Newton-Krylov optimization using the incremental adjoint Jacobi equations. Phys Med Biol 2019; 64:025002. [PMID: 30523830 DOI: 10.1088/1361-6560/aaf598] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
The class of non-rigid registration methods proposed in the framework of PDE-constrained large deformation diffeomorphic metric mapping is a particularly interesting family of physically meaningful diffeomorphic registration methods. Inexact Gauss-Newton-Krylov optimization has shown an excellent numerical accuracy and an extraordinarily fast convergence rate in this framework. However, the Galerkin representation of the non-stationary velocity fields does not provide proper geodesic paths. In this work, we propose a method for PDE-constrained LDDMM parameterized in the space of initial velocity fields under the EPDiff equation. The derivation of the gradient and the Hessian-vector products are performed on the final velocity field and transported backward using the adjoint and the incremental adjoint Jacobi equations. This way, we avoid the complex dependence on the initial velocity field in the computations. We also avoid the computation of the adjoint equation and its incremental counterpart that has been recently identified as a subtle problem in PDE-constrained LDDMM. The proposed method provides geodesics in the framework of PDE-constrained LDDMM, and it shows performance competing with benchmark PDE-constrained LDDMM and EPDiff-LDDMM methods.
Collapse
Affiliation(s)
- Monica Hernandez
- Department of Computer Science, Aragon Institute on Engineering Research (I3A), University of Zaragoza, Zaragoza, Spain
| |
Collapse
|
3
|
Abstract
The class of registration methods proposed in the framework of Stokes large deformation diffeomorphic metric mapping (LDDMM) is a particularly interesting family of physically meaningful diffeomorphic registration methods. Stokes-LDDMM methods are formulated as constrained variational problems, where the different physical models are imposed using the associated partial differential equations as hard constraints. The most significant limitation of Stokes-LDDMM framework is its huge computational complexity. The objective of this paper is to promote the use of Stokes-LDDMM in computational anatomy applications with an efficient approximation of the original variational problem. Thus, we propose a novel method for efficient Stokes-LDDMM diffeomorphic registration. Our method poses the constrained variational problem in the space of band-limited vector fields and it is implemented in the GPU. The performance of band-limited Stokes-LDDMM has been compared and evaluated with original Stokes-LDDMM, EPDiff-LDDMM, and band-limited EPDiff-LDDMM. The evaluation has been conducted in 3-D with the nonrigid image registration evaluation project database. Since the update equation in Stokes-LDDMM involves the action of low-pass filters, the computational complexity has been greatly alleviated with a modest accuracy lose. We have obtained a competitive performance for some method configurations. Overall, our proposed method may make feasible the extensive use of novel physically meaningful Stokes-LDDMM methods in different computational anatomy applications. In addition, our results reinforce the usefulness of band-limited vector fields in diffeomorphic registration methods involving the action of low-pass filters in the optimization, even in algorithmically challenging environments such as Stokes-LDDMM.
Collapse
|
4
|
Holm DD, Tyranowski TM. Variational principles for stochastic soliton dynamics. Proc Math Phys Eng Sci 2016; 472:20150827. [PMID: 27118922 DOI: 10.1098/rspa.2015.0827] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
We develop a variational method of deriving stochastic partial differential equations whose solutions follow the flow of a stochastic vector field. As an example in one spatial dimension, we numerically simulate singular solutions (peakons) of the stochastically perturbed Camassa-Holm (CH) equation derived using this method. These numerical simulations show that peakon soliton solutions of the stochastically perturbed CH equation persist and provide an interesting laboratory for investigating the sensitivity and accuracy of adding stochasticity to finite dimensional solutions of stochastic partial differential equations. In particular, some choices of stochastic perturbations of the peakon dynamics by Wiener noise (canonical Hamiltonian stochastic deformations, CH-SD) allow peakons to interpenetrate and exchange order on the real line in overtaking collisions, although this behaviour does not occur for other choices of stochastic perturbations which preserve the Euler-Poincaré structure of the CH equation (parametric stochastic deformations, P-SD), and it also does not occur for peakon solutions of the unperturbed deterministic CH equation. The discussion raises issues about the science of stochastic deformations of finite-dimensional approximations of evolutionary partial differential equation and the sensitivity of the resulting solutions to the choices made in stochastic modelling.
Collapse
Affiliation(s)
- Darryl D Holm
- Mathematics Department , Imperial College , London, UK
| | | |
Collapse
|
5
|
Diffeomorphic Metric Landmark Mapping Using Stationary Velocity Field Parameterization. Int J Comput Vis 2015. [DOI: 10.1007/s11263-015-0802-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
6
|
Reduction by Lie Group Symmetries in Diffeomorphic Image Registration and Deformation Modelling. Symmetry (Basel) 2015. [DOI: 10.3390/sym7020599] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
|
7
|
Bruveris M, Holm DD. Geometry of Image Registration: The Diffeomorphism Group and Momentum Maps. GEOMETRY, MECHANICS, AND DYNAMICS 2015. [DOI: 10.1007/978-1-4939-2441-7_2] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
8
|
Hernandez M, Bossa MN, Olmos S. Registration of Anatomical Images Using Paths of Diffeomorphisms Parameterized with Stationary Vector Field Flows. Int J Comput Vis 2009. [DOI: 10.1007/s11263-009-0219-z] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
9
|
Holm DD, O Náraigh L, Tronci C. Singular solutions of a modified two-component Camassa-Holm equation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 79:016601. [PMID: 19257154 DOI: 10.1103/physreve.79.016601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/13/2008] [Indexed: 05/27/2023]
Abstract
The Camassa-Holm (CH) equation is a well-known integrable equation describing the velocity dynamics of shallow water waves. This equation exhibits spontaneous emergence of singular solutions (peakons) from smooth initial conditions. The CH equation has been recently extended to a two-component integrable system (CH2), which includes both velocity and density variables in the dynamics. Although possessing peakon solutions in the velocity, the CH2 equation does not admit singular solutions in the density profile. We modify the CH2 system to allow a dependence on the average density as well as the pointwise density. The modified CH2 system (MCH2) does admit peakon solutions in the velocity and average density. We analytically identify the steepening mechanism that allows the singular solutions to emerge from smooth spatially confined initial data. Numerical results for the MCH2 system are given and compared with the pure CH2 case. These numerics show that the modification in the MCH2 system to introduce the average density has little short-time effect on the emergent dynamical properties. However, an analytical and numerical study of pairwise peakon interactions for the MCH2 system shows a different asymptotic feature. Namely, besides the expected soliton scattering behavior seen in overtaking and head-on peakon collisions, the MCH2 system also allows the phase shift of the peakon collision to diverge in certain parameter regimes.
Collapse
Affiliation(s)
- Darryl D Holm
- Department of Mathematics, Imperial College London, 180 Queen's Gate, London SW7 2AZ, United Kingdom
| | | | | |
Collapse
|
10
|
Younes L, Arrate F, Miller MI. Evolutions equations in computational anatomy. Neuroimage 2008; 45:S40-50. [PMID: 19059343 DOI: 10.1016/j.neuroimage.2008.10.050] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2008] [Accepted: 10/15/2008] [Indexed: 10/21/2022] Open
Abstract
One of the main purposes in computational anatomy is the measurement and statistical study of anatomical variations in organs, notably in the brain or the heart. Over the last decade, our group has progressively developed several approaches for this problem, all related to the Riemannian geometry of groups of diffeomorphisms and the shape spaces on which these groups act. Several important shape evolution equations that are now used routinely in applications have emerged over time. Our goal in this paper is to provide an overview of these equations, placing them in their theoretical context, and giving examples of applications in which they can be used. We introduce the required theoretical background before discussing several classes of equations of increasingly complexity. These equations include energy minimizing evolutions deriving from Riemannian gradient descent, geodesics, parallel transport and Jacobi fields.
Collapse
Affiliation(s)
- Laurent Younes
- Center for Imaging Science, The Johns Hopkins University, 3400N Charles St., Baltimore, MD 21218, USA.
| | | | | |
Collapse
|
11
|
Holm DD, Tronci C. Geodesic flows on semidirect-product Lie groups: geometry of singular measure-valued solutions. Proc Math Phys Eng Sci 2008. [DOI: 10.1098/rspa.2008.0263] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The EPDiff equation (or the dispersionless Camassa–Holm equation in one dimension) is a well-known example of geodesic motion on the Diff group of smooth invertible maps (diffeomorphisms). Its recent two-component extension governs geodesic motion on the semidirect product DiffⓈ
, where
denotes the space of scalar functions. This paper generalizes the second construction to consider geodesic motion on DiffⓈ
, where
denotes the space of scalar functions that take values on a certain Lie algebra (e.g.
=
⊗
(3)). Measure-valued delta-like solutions are shown to be momentum maps possessing a dual pair structure, thereby extending previous results for the EPDiff equation. The collective Hamiltonians are shown to fit into the Kaluza–Klein theory of particles in a Yang–Mills field and these formulations are shown to apply also at the continuum partial differential equation level. In the continuum description, the Kaluza–Klein approach produces the Kelvin circulation theorem.
Collapse
Affiliation(s)
- Darryl D Holm
- Department of Mathematics, Imperial College London180 Queen's Gate, London SW7 2AZ, UK
- Institute for Mathematical Sciences, Imperial College London53 Prince's Gate, London SW7 2PG, UK
| | - Cesare Tronci
- Department of Mathematics, Imperial College London180 Queen's Gate, London SW7 2AZ, UK
- TERA Foundation for Oncological Hadrontherapy11 Via Puccini, Novara 28100, Italy
| |
Collapse
|
12
|
Younes L, Qiu A, Winslow RL, Miller MI. Transport of Relational Structures in Groups of Diffeomorphisms. JOURNAL OF MATHEMATICAL IMAGING AND VISION 2008; 32:41-56. [PMID: 19809583 PMCID: PMC2756699 DOI: 10.1007/s10851-008-0074-5] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
This paper focuses on the issue of translating the relative variation of one shape with respect to another in a template centered representation. The context is the theory of Diffeomorphic Pattern Matching which provides a representation of the space of shapes of objects, including images and point sets, as an infinite dimensional Riemannian manifold which is acted upon by groups of diffeomorphisms. We discuss two main options for achieving our goal; the first one is the parallel translation, based on the Riemannian metric; the second one, based on the group action, is the coadjoint transport. These methods are illustrated with 3D experiments.
Collapse
Affiliation(s)
- Laurent Younes
- Center for Imaging Science, Johns Hopkins University, 3400 N. Charles St., Baltimore, MD 21218, USA
| | - Anqi Qiu
- Center for Imaging Science, Johns Hopkins University, 3400 N. Charles St., Baltimore, MD 21218, USA
| | - Raimond L. Winslow
- Department of Biomedical Engineering, Johns Hopkins University, 3400 N. Charles St., Baltimore, MD 21218, USA
| | - Michael I. Miller
- Center for Imaging Science, Johns Hopkins University, 3400 N. Charles St., Baltimore, MD 21218, USA
| |
Collapse
|
13
|
Lepore N, Brun C, Chou YY, Chiang MC, Dutton RA, Hayashi KM, Luders E, Lopez OL, Aizenstein HJ, Toga AW, Becker JT, Thompson PM. Generalized tensor-based morphometry of HIV/AIDS using multivariate statistics on deformation tensors. IEEE TRANSACTIONS ON MEDICAL IMAGING 2008; 27:129-41. [PMID: 18270068 PMCID: PMC2832297 DOI: 10.1109/tmi.2007.906091] [Citation(s) in RCA: 103] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
This paper investigates the performance of a new multivariate method for tensor-based morphometry (TBM). Statistics on Riemannian manifolds are developed that exploit the full information in deformation tensor fields. In TBM, multiple brain images are warped to a common neuroanatomical template via 3-D nonlinear registration; the resulting deformation fields are analyzed statistically to identify group differences in anatomy. Rather than study the Jacobian determinant (volume expansion factor) of these deformations, as is common, we retain the full deformation tensors and apply a manifold version of Hotelling's $T(2) test to them, in a Log-Euclidean domain. In 2-D and 3-D magnetic resonance imaging (MRI) data from 26 HIV/AIDS patients and 14 matched healthy subjects, we compared multivariate tensor analysis versus univariate tests of simpler tensor-derived indices: the Jacobian determinant, the trace, geodesic anisotropy, and eigenvalues of the deformation tensor, and the angle of rotation of its eigenvectors. We detected consistent, but more extensive patterns of structural abnormalities, with multivariate tests on the full tensor manifold. Their improved power was established by analyzing cumulative p-value plots using false discovery rate (FDR) methods, appropriately controlling for false positives. This increased detection sensitivity may empower drug trials and large-scale studies of disease that use tensor-based morphometry.
Collapse
Affiliation(s)
- N Lepore
- Department of Neurology, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA.
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
14
|
Cotter C, Holm D, Hydon P. Multisymplectic formulation of fluid dynamics using the inverse map. Proc Math Phys Eng Sci 2007. [DOI: 10.1098/rspa.2007.1892] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
We construct multisymplectic formulations of fluid dynamics using the inverse of the Lagrangian path map. This inverse map, the ‘back-to-labels’ map, gives the initial Lagrangian label of the fluid particle that currently occupies each Eulerian position. Explicitly enforcing the condition that the fluid particles carry their labels with the flow in Hamilton's principle leads to our multisymplectic formulation. We use the multisymplectic one-form to obtain conservation laws for energy, momentum and an infinite set of conservation laws arising from the particle relabelling symmetry and leading to Kelvin's circulation theorem. We discuss how multisymplectic numerical integrators naturally arise in this approach.
Collapse
Affiliation(s)
- C.J Cotter
- Department of Aeronautics, Imperial College LondonLondon SW7 2AZ, UK
| | - D.D Holm
- Department of Mathematics, Imperial College LondonLondon SW7 2AZ, UK
| | - P.E Hydon
- Department of Mathematics, University of SurreyGuildford, Surrey GU2 7XH, UK
| |
Collapse
|
15
|
Marsland S, McLachlan R. A hamiltonian particle method for diffeomorphic image registration. INFORMATION PROCESSING IN MEDICAL IMAGING : PROCEEDINGS OF THE ... CONFERENCE 2007; 20:396-407. [PMID: 17633716 DOI: 10.1007/978-3-540-73273-0_33] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Diffeomorphic image registration, where images are aligned using diffeomorphic warps, is a popular subject for research in medical image analysis. We introduce a novel algorithm for computing diffeomorphic warps that solves the Euler equations on the diffeomorphism group explicitly, based on a discretisation of the Hamiltonian, rather than using an optimiser. The result is an algorithm that is many times faster than those considered previously.
Collapse
Affiliation(s)
- Stephen Marsland
- Massey University, Private Bag 11-222, Palmerston North, New Zealand.
| | | |
Collapse
|
16
|
Leow AD, Yanovsky I, Chiang MC, Lee AD, Klunder AD, Lu A, Becker JT, Davis SW, Toga AW, Thompson PM. Statistical properties of Jacobian maps and the realization of unbiased large-deformation nonlinear image registration. IEEE TRANSACTIONS ON MEDICAL IMAGING 2007; 26:822-32. [PMID: 17679333 DOI: 10.1109/tmi.2007.892646] [Citation(s) in RCA: 137] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Maps of local tissue compression or expansion are often computed by comparing magnetic resonance imaging (MRI) scans using nonlinear image registration. The resulting changes are commonly analyzed using tensor-based morphometry to make inferences about anatomical differences, often based on the Jacobian map, which estimates local tissue gain or loss. Here, we provide rigorous mathematical analyses of the Jacobian maps, and use themto motivate a new numerical method to construct unbiased nonlinear image registration. First, we argue that logarithmic transformation is crucial for analyzing Jacobian values representing morphometric differences. We then examine the statistical distributions of log-Jacobian maps by defining the Kullback-Leibler (KL) distance on material density functions arising in continuum-mechanical models. With this framework, unbiased image registration can be constructed by quantifying the symmetric KL-distance between the identity map and the resulting deformation. Implementation details, addressing the proposed unbiased registration as well as the minimization of symmetric image matching functionals, are then discussed and shown to be applicable to other registration methods, such as inverse consistent registration. In the results section, we test the proposed framework, as well as present an illustrative application mapping detailed 3-D brain changes in sequential magnetic resonance imaging scans of a patient diagnosed with semantic dementia. Using permutation tests, we show that the symmetrization of image registration statistically reduces skewness in the log-Jacobian map.
Collapse
Affiliation(s)
- Alex D Leow
- Neuropsychiatric Hospital and the Laboratory of Neuro Imaging, Department of Neurology, UCLA David Geffen School of Medicine, Los Angeles, CA 90095, USA.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
17
|
Shi Y, Thompson PM, de Zubicaray GI, Rose SE, Tu Z, Dinov I, Toga AW. Direct mapping of hippocampal surfaces with intrinsic shape context. Neuroimage 2007; 37:792-807. [PMID: 17625918 PMCID: PMC2227952 DOI: 10.1016/j.neuroimage.2007.05.016] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2006] [Revised: 05/05/2007] [Accepted: 05/10/2007] [Indexed: 10/23/2022] Open
Abstract
We propose in this paper a new method for the mapping of hippocampal (HC) surfaces to establish correspondences between points on HC surfaces and enable localized HC shape analysis. A novel geometric feature, the intrinsic shape context, is defined to capture the global characteristics of the HC shapes. Based on this intrinsic feature, an automatic algorithm is developed to detect a set of landmark curves that are stable across population. The direct map between a source and target HC surface is then solved as the minimizer of a harmonic energy function defined on the source surface with landmark constraints. For numerical solutions, we compute the map with the approach of solving partial differential equations on implicit surfaces. The direct mapping method has the following properties: (1) it has the advantage of being automatic; (2) it is invariant to the pose of HC shapes. In our experiments, we apply the direct mapping method to study temporal changes of HC asymmetry in Alzheimer's disease (AD) using HC surfaces from 12 AD patients and 14 normal controls. Our results show that the AD group has a different trend in temporal changes of HC asymmetry than the group of normal controls. We also demonstrate the flexibility of the direct mapping method by applying it to construct spherical maps of HC surfaces. Spherical harmonics (SPHARM) analysis is then applied and it confirms our results on temporal changes of HC asymmetry in AD.
Collapse
Affiliation(s)
- Yonggang Shi
- Laboratory of Neuro Imaging, Department of Neurology, UCLA School of Medicine, Los Angeles, CA 90095, USA 1
| | - Paul M. Thompson
- Laboratory of Neuro Imaging, Department of Neurology, UCLA School of Medicine, Los Angeles, CA 90095, USA 1
| | - Greig I. de Zubicaray
- Centre for Magnetic Resonance, University of Queensland, Brisbane, QLD 4072, Australia
| | - Stephen E. Rose
- Centre for Magnetic Resonance, University of Queensland, Brisbane, QLD 4072, Australia
| | - Zhuowen Tu
- Laboratory of Neuro Imaging, Department of Neurology, UCLA School of Medicine, Los Angeles, CA 90095, USA 1
| | - Ivo Dinov
- Laboratory of Neuro Imaging, Department of Neurology, UCLA School of Medicine, Los Angeles, CA 90095, USA 1
| | - Arthur W. Toga
- Laboratory of Neuro Imaging, Department of Neurology, UCLA School of Medicine, Los Angeles, CA 90095, USA 1
- Laboratory of Neuro Imaging, Department of Neurology, UCLA School of Medicine, Los Angeles, CA 90095, USA, Email address: (Arthur W. Toga)
| |
Collapse
|
18
|
Avants BB, Schoenemann PT, Gee JC. Lagrangian frame diffeomorphic image registration: Morphometric comparison of human and chimpanzee cortex. Med Image Anal 2006; 10:397-412. [PMID: 15948659 DOI: 10.1016/j.media.2005.03.005] [Citation(s) in RCA: 121] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2004] [Revised: 02/09/2005] [Accepted: 03/04/2005] [Indexed: 11/23/2022]
Abstract
We develop a novel Lagrangian reference frame diffeomorphic image and landmark registration method. The algorithm uses the fixed Langrangian reference frame to define the map between coordinate systems, but also generates and stores the inverse map from the Eulerian to the Lagrangian frame. Computing both maps allows facile computation of both Eulerian and Langrangian quantities. We apply this algorithm to estimating a putative evolutionary change of coordinates between a population of chimpanzee and human cortices. Inter-species functional homologues fix the map explicitly, where they are known, while image similarities guide the alignment elsewhere. This map allows detailed study of the volumetric change between chimp and human cortex. Instead of basing the inter-species study on a single species atlas, we diffeomorphically connect the mean shape and intensity templates for each group. The human statistics then map diffeomorphically into the space of the chimpanzee cortex providing a comparison between species. The population statistics show a significant doubling of the relative prefrontal lobe size in humans, as compared to chimpanzees.
Collapse
Affiliation(s)
- Brian B Avants
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA 19104-6389, USA.
| | | | | |
Collapse
|
19
|
Vaillant M, Miller MI, Younes L, Trouvé A. Statistics on diffeomorphisms via tangent space representations. Neuroimage 2004; 23 Suppl 1:S161-9. [PMID: 15501085 DOI: 10.1016/j.neuroimage.2004.07.023] [Citation(s) in RCA: 102] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2004] [Accepted: 07/01/2004] [Indexed: 11/27/2022] Open
Abstract
In this paper, we present a linear setting for statistical analysis of shape and an optimization approach based on a recent derivation of a conservation of momentum law for the geodesics of diffeomorphic flow. Once a template is fixed, the space of initial momentum becomes an appropriate space for studying shape via geodesic flow since the flow at any point along the geodesic is completely determined by the momentum at the origin through geodesic shooting equations. The space of initial momentum provides a linear representation of the nonlinear diffeomorphic shape space in which linear statistical analysis can be applied. Specializing to the landmark matching problem of Computational Anatomy, we derive an algorithm for solving the variational problem with respect to the initial momentum and demonstrate principal component analysis (PCA) in this setting with three-dimensional face and hippocampus databases.
Collapse
Affiliation(s)
- M Vaillant
- Center for Imaging Science, The Johns Hopkins University, Baltimore, MD 21218, USA.
| | | | | | | |
Collapse
|