Multiscale structural analysis of mouse lingual myoarchitecture employing diffusion spectrum magnetic resonance imaging and multiphoton microscopy

Lipid and smooth muscle architectural pathology in the rabbit atherosclerotic vessel wall using Q-space cardiovascular magnetic resonance

BACKGROUND

Atherosclerosis is an arterial vessel wall disease characterized by slow, progressive lipid accumulation, smooth muscle disorganization, and inflammatory infiltration. Atherosclerosis often remains subclinical until extensive inflammatory injury promotes vulnerability of the atherosclerotic plaque to rupture with luminal thrombosis, which can cause the acute event of myocardial infarction or stroke. Current bioimaging techniques are unable to capture the pathognomonic distribution of cellular elements of the plaque and thus cannot accurately define its structural disorganization.

METHODS

We applied cardiovascular magnetic resonance spectroscopy (CMRS) and diffusion weighted CMR (DWI) with generalized Q-space imaging (GQI) analysis to architecturally define features of atheroma and correlated these to the microscopic distribution of vascular smooth muscle cells (SMC), immune cells, extracellular matrix (ECM) fibers, thrombus, and cholesteryl esters (CE). We compared rabbits with normal chow diet and cholesterol-fed rabbits with endothelial balloon injury, which accelerates atherosclerosis and produces advanced rupture-prone plaques, in a well-validated rabbit model of human atherosclerosis.

RESULTS

Our methods revealed new structural properties of advanced atherosclerosis incorporating SMC and lipid distributions. GQI with tractography portrayed the locations of these components across the atherosclerotic vessel wall and differentiated multi-level organization of normal, pro-inflammatory cellular phenotypes, or thrombus. Moreover, the locations of CE were differentiated from cellular constituents by their higher restrictive diffusion properties, which permitted chemical confirmation of CE by high field voxel-guided CMRS.

CONCLUSIONS

GQI with tractography is a new method for atherosclerosis imaging that defines a pathological architectural signature for the atheromatous plaque composed of distributed SMC, ECM, inflammatory cells, and thrombus and lipid. This provides a detailed transmural map of normal and inflamed vessel walls in the setting of atherosclerosis that has not been previously achieved using traditional CMR techniques. Although this is an ex-vivo study, detection of micro and mesoscale level vascular destabilization as enabled by GQI with tractography could increase the accuracy of diagnosis and assessment of treatment outcomes in individuals with atherosclerosis.

Validation of Muscle Fiber Architecture of the Human Tongue Revealed by Diffusion Magnetic Resonance Imaging With Histology Verification

PURPOSE

The goal of this study is to validate the muscle architecture derived from both ex vivo and in vivo diffusion-weighted magnetic resonance imaging (dMRI) of the human tongue with histology of an ex vivo tongue.

METHOD

dMRI was acquired with a 200-direction high angular resolution diffusion imaging (HARDI) diffusion scheme for both a postmortem head (imaged within 48 hr after death) and a healthy volunteer. After MRI, the postmortem head was fixed and the tongue excised for hematoxylin and eosin (H&E) staining and histology imaging. Structure tensor images were generated from the stained images to better demonstrate muscle fiber orientations. The tongue muscle fiber orientations, estimated from dMRI, were visualized using the tractogram, a novel representation of crossing fiber orientations, and compared against the histology images of the ex vivo tongue.

RESULTS

Muscle fibers identified in the tractograms showed good correspondence with those appearing in the histology images. We further demonstrated tongue muscle architecture in in vivo tractograms for the entire tongue.

CONCLUSION

The study demonstrates that dMRI can accurately reveal the complex muscle architecture of the human tongue and may potentially benefit planning and evaluation of oral surgery and research on speech and swallowing.

Two dimensional computational model coupling myoarchitecture-based lingual tissue mechanics with liquid bolus flow during oropharyngeal swallowing

Biomechanical relationships involving lingual myoanatomy, contractility, and bolus movement are fundamental properties of human swallowing. To portray the relationship between lingual deformation and bolus flow during swallowing, a weakly one-way solid-fluid finite element model (FEM) was derived employing an elemental mesh aligned to magnetic resonance diffusional tractography (Q-space MRI, QSI) of the human tongue, an arbitrary Lagrangian-Eulerian (ALE) formulation with remeshing to account for the effects of lingual surface (boundary) deformation, an implementation of patterned fiber shortening, and a computational visualization of liquid bolus flow. Representing lingual tissue deformation in terms of its 2D principal Lagrangian strain in the mid-sagittal plane, we demonstrated that the swallow sequence was characterized by initial superior-anterior expansion directed towards the hard palate, followed by sequential, radially directed, contractions of the genioglossus and verticalis to promote lingual rotation (lateral perspective) and propulsive displacement. We specifically assessed local bolus velocity as a function of viscosity (perfect slip conditions) and observed that a low viscosity bolus (5 cP) exhibited maximal displacement, surface spreading and local velocity compared to medium (110 cP, 300 cP) and high (525 cP) viscosity boluses. Analysis of local nodal velocity revealed that all bolus viscosities exhibited a bi-phasic progression, with the low viscosity bolus being the most heterogeneous and fragmented and the high viscosity bolus being the most homogenous and cohesive. Intraoral bolus cohesion was depicted in terms of the distributed velocity gradient, with higher gradients being associated with increased shear rate and bolus fragmentation. Lastly, we made a sensitivity analysis on tongue stiffness and contractility by varying the degree of extracellular matrix (ECM) stiffness through effects on the Mooney-Rivlin derived passive matrix and by varying maximum tetanized isometric stress, and observed that a graded increase of ECM stiffness was associated with reduced bolus spreading, posterior displacement, and surface velocity gradients, whereas a reduction of global contractility resulted in a graded reduction of obtainable accommodation volume, absent bolus spreading, and loss of posterior displacement. We portray a unidirectionally coupled solid-liquid FEM which associates myoarchitecture-based lingual deformation with intra-oral bolus flow, and deduce that local elevation of the velocity gradient correlates with bolus fragmentation, a precondition believed to be associated with aspiration vulnerability during oropharyngeal swallowing.

Crossing muscle fibers of the human tongue resolved in vivo using constrained spherical deconvolution

BACKGROUND

Surgical resection of tongue cancer may impair swallowing and speech. Knowledge of tongue muscle architecture affected by the resection could aid in patient counseling. Diffusion tensor imaging (DTI) enables reconstructions of muscle architecture in vivo. Reconstructing crossing fibers in the tongue requires a higher-order diffusion model.

PURPOSE

To develop a clinically feasible diffusion imaging protocol, which facilitates both DTI and constrained spherical deconvolution (CSD) reconstructions of tongue muscle architecture in vivo.

STUDY TYPE

Cross-sectional study.

SUBJECTS/SPECIMEN

One ex vivo bovine tongue resected en bloc from mandible to hyoid bone. Ten healthy volunteers (mean age 25.5 years; range 21-34 years; four female).

FIELD STRENGTH/SEQUENCE

Diffusion-weighted echo planar imaging at 3 T using a high-angular resolution diffusion imaging scheme acquired twice with opposing phase-encoding for B₀ -field inhomogeneity correction. The scan of the healthy volunteers was divided into four parts, in between which the volunteers were allowed to swallow, resulting in a total acquisition time of 10 minutes.

ASSESSMENT

The ability of resolving crossing muscle fibers using CSD was determined on the bovine tongue specimen. A reproducible response function was estimated and the optimal peak threshold was determined for the in vivo tongue. The quality of tractography of the in vivo tongue was graded by three experts.

STATISTICAL TESTS

The within-subject coefficient of variance was calculated for the response function. The qualitative results of the grading of DTI and CSD tractography were analyzed using a multilevel proportional odds model.

RESULTS

Fiber orientation distributions in the bovine tongue specimen showed that CSD was able to resolve crossing muscle fibers. The response function could be determined reproducibly in vivo. CSD tractography displayed significantly improved tractography compared with DTI tractography (P = 0.015).

DATA CONCLUSION

The 10-minute diffusion imaging protocol facilitates CSD fiber tracking with improved reconstructions of crossing tongue muscle fibers compared with DTI.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:96-105.

Association between tumor architecture derived from generalized Q-space MRI and survival in glioblastoma

While it is recognized that the overall resistance of glioblastoma to treatment may be related to intra-tumor patterns of structural heterogeneity, imaging methods to assess such patterns remain rudimentary.

METHODS

We utilized a generalized Q-space imaging (GQI) algorithm to analyze magnetic resonance imaging (MRI) derived from a rodent model of glioblastoma and 2 clinical datasets to correlate GQI, histology, and survival.

RESULTS

In a rodent glioblastoma model, GQI demonstrated a poorly coherent core region, consisting of diffusion tracts <5 mm, surrounded by a shell of highly coherent diffusion tracts, 6-25 mm. Histologically, the core region possessed a high degree of necrosis, whereas the shell consisted of organized sheets of anaplastic cells with elevated mitotic index. These attributes define tumor architecture as the macroscopic organization of variably aligned tumor cells. Applied to MRI data from The Cancer Imaging Atlas (TCGA), the core-shell diffusion tract-length ratio (c/s ratio) correlated linearly with necrosis, which, in turn, was inversely associated with survival (p = 0.00002). We confirmed in an independent cohort of patients (n = 62) that the c/s ratio correlated inversely with survival (p = 0.0004).

CONCLUSIONS

The analysis of MR images by GQI affords insight into tumor architectural patterns in glioblastoma that correlate with biological heterogeneity and clinical outcome.

Assessing the multiscale architecture of muscular tissue with Q-space magnetic resonance imaging: Review

Contraction of muscular tissue requires the synchronized shortening of myofibers arrayed in complex geometrical patterns. Imaging such myofiber patterns with diffusion-weighted MRI reveals architectural ensembles that underlie force generation at the organ scale. Restricted proton diffusion is a stochastic process resulting from random translational motion that may be used to probe the directionality of myofibers in whole tissue. During diffusion-weighted MRI, magnetic field gradients are applied to determine the directional dependence of proton diffusion through the analysis of a diffusional probability distribution function (PDF). The directions of principal (maximal) diffusion within the PDF are associated with similarly aligned diffusion maxima in adjacent voxels to derive multivoxel tracts. Diffusion-weighted MRI with tractography thus constitutes a multiscale method for depicting patterns of cellular organization within biological tissues. We provide in this review, details of the method by which generalized Q-space imaging is used to interrogate multidimensional diffusion space, and thereby to infer the organization of muscular tissue. Q-space imaging derives the lowest possible angular separation of diffusion maxima by optimizing the conditions by which magnetic field gradients are applied to a given tissue. To illustrate, we present the methods and applications associated with Q-space imaging of the multiscale myoarchitecture associated with the human and rodent tongues. These representations emphasize the intricate and continuous nature of muscle fiber organization and suggest a method to depict structural "blueprints" for skeletal and cardiac muscle tissue.

Alterations in Multi-Scale Cardiac Architecture in Association With Phosphorylation of Myosin Binding Protein-C

Background

The geometric organization of myocytes in the ventricular wall comprises the structural underpinnings of cardiac mechanical function. Cardiac myosin binding protein‐C (MYBPC3) is a sarcomeric protein, for which phosphorylation modulates myofilament binding, sarcomere morphology, and myocyte alignment in the ventricular wall. To elucidate the mechanisms by which MYBPC3 phospho‐regulation affects cardiac tissue organization, we studied ventricular myoarchitecture using generalized Q‐space imaging (GQI). GQI assessed geometric phenotype in excised hearts that had undergone transgenic (TG) modification of phospho‐regulatory serine sites to nonphosphorylatable alanines (MYBPC3^{AllP−/(t/t)}) or phospho‐mimetic aspartic acids (MYBPC3^AllP+/(t/t)).

Methods and Results

Myoarchitecture in the wild‐type (MYBPC3^WT) left‐ventricle (LV) varied with transmural position, with helix angles ranging from −90/+90 degrees and contiguous circular orientation from the LV mid‐myocardium to the right ventricle (RV). Whereas MYBPC3^AllP+/(t/t) hearts were not architecturally distinct from MYBPC3^WT, MYBPC3^{AllP−/(t/t)} hearts demonstrated a significant reduction in LV transmural helicity. Null MYBPC3^(t/t) hearts, as constituted by a truncated MYBPC3 protein, demonstrated global architectural disarray and loss in helicity. Electron microscopy was performed to correlate the observed macroscopic architectural changes with sarcomere ultrastructure and demonstrated that impaired phosphorylation of MYBPC3 resulted in modifications of the sarcomere aspect ratio and shear angle. The mechanical effect of helicity loss was assessed through a geometric model relating cardiac work to ejection fraction, confirming the mechanical impairments observed with echocardiography.

Conclusions

We conclude that phosphorylation of MYBPC3 contributes to the genesis of ventricular wall geometry, linking myofilament biology with multiscale cardiac mechanics and myoarchitecture.

Patterns of intersecting fiber arrays revealed in whole muscle with generalized Q-space imaging

The multiscale attributes of mammalian muscle confer significant challenges for structural imaging in vivo. To achieve this, we employed a magnetic resonance method, termed "generalized Q-space imaging", that considers the effect of spatially distributed diffusion-weighted magnetic field gradients and diffusion sensitivities on the morphology of Q-space. This approach results in a subvoxel scaled probability distribution function whose shape correlates with local fiber orientation. The principal fiber populations identified within these probability distribution functions can then be associated by streamline methods to create multivoxel tractlike constructs that depict the macroscale orientation of myofiber arrays. We performed a simulation of Q-space input parameters, including magnetic field gradient strength and direction, diffusion sensitivity, and diffusional sampling to determine the optimal achievable fiber angle separation in the minimum scan time. We applied this approach to resolve intravoxel crossing myofiber arrays in the setting of the human tongue, an organ with anatomic complexity based on the presence of hierarchical arrays of intersecting myocytes. Using parameters defined by simulation, we imaged at 3T the fanlike configuration of the human genioglossus and the laterally positioned merging fibers of the styloglossus, inferior longitudinalis, chondroglossus, and verticalis. Comparative scans of the excised mouse tongue at 7T demonstrated similar midline and lateral crossing fiber patterns, whereas histological analysis confirmed the presence and distribution of these myofiber arrays at the microscopic scale. Our results demonstrate a magnetic resonance method for acquiring and displaying diffusional data that defines highly ordered myofiber patterns in architecturally complex tissue. Such patterns suggest inherent multiscale fiber organization and provide a basis for structure-function analyses in vivo and in model tissues.

Resolving intravoxel fiber architecture using nonconvex regularized blind compressed sensing

In diffusion magnetic resonance imaging, accurate and reliable estimation of intravoxel fiber architectures is a major prerequisite for tractography algorithms or any other derived statistical analysis. Several methods have been proposed that estimate intravoxel fiber architectures using low angular resolution acquisitions owing to their shorter acquisition time and relatively low b-values. But these methods are highly sensitive to noise. In this work, we propose a nonconvex regularized blind compressed sensing approach to estimate intravoxel fiber architectures in low angular resolution acquisitions. The method models diffusion-weighted (DW) signals as a sparse linear combination of unfixed reconstruction basis functions and introduces a nonconvex regularizer to enhance the noise immunity. We present a general solving framework to simultaneously estimate the sparse coefficients and the reconstruction basis. Experiments on synthetic, phantom, and real human brain DW images demonstrate the superiority of the proposed approach.

Purse-string morphology of external anal sphincter revealed by novel imaging techniques

The external anal sphincter (EAS) may be injured in 25-35% of women during the first and subsequent vaginal childbirths and is likely the most common cause of anal incontinence. Since its first description almost 300 years ago, the EAS was believed to be a circular or a "donut-shaped" structure. Using three-dimensional transperineal ultrasound imaging, MRI, diffusion tensor imaging, and muscle fiber tracking, we delineated various components of the EAS and their muscle fiber directions. These novel imaging techniques suggest "purse-string" morphology, with "EAS muscles" crossing contralaterally in the perineal body to the contralateral transverse perineal (TP) and bulbospongiosus (BS) muscles, thus attaching the EAS to the pubic rami. Spin-tag MRI demonstrated purse-string action of the EAS muscle. Electromyography of TP/BS and EAS muscles revealed their simultaneous contraction and relaxation. Lidocaine injection into the TP/BS muscle significantly reduced anal canal pressure. These studies support purse-string morphology of the EAS to constrict/close the anal canal opening. Our findings have implications for the effect of episiotomy on anal closure function and the currently used surgical technique (overlapping sphincteroplasty) for EAS reconstructive surgery to treat anal incontinence.

In vivo 3 T MR diffusion tensor imaging for detection of the fibre architecture of the human uterus: a feasibility and quantitative study

OBJECTIVE

The aim of this study was to investigate the feasibility of depicting fibre architecture of human uteri in vivo using 3 T MR diffusion tensor imaging (MR-DTI) with a three-dimensional (3D) tractography approach. Quantitative results were provided.

METHODS

In vivo 3 T MR-DTI was performed on 30 volunteers (9 Caesarean delivery). Main diffusion directions reflecting the fibre orientation were determined using sensitivity-encoding single-shot echo planar imaging with diffusion-sensitised gradients (b=600 mm(2) s(-1)) along 32 directions. A deterministic fibre-tracking algorithm was used to show in vivo fibre architecture, compared with ex vivo histological slides of cadaveric uteri. The number of fibres, the fibre density, the fractional anisotropy (FA) and the apparent diffusion coefficient (ADC) were measured in 13 volunteers.

RESULTS

Anisotropy was found in most regions of normal uteri and the preferential order of uterine fibres depicted, consisting of two representative fibre directions: circular and longitudinal, as in ex vivo studies. Two-thirds of uteri with a Caesarean scar did not have the same orientation of fibres in the anterior isthmus when compared with non-scarred myometrium. Quantitative data were obtained from 13 volunteers: Caesarean-scarred uteri (n=5) showed lower fibre number and density in the scarred anterior isthmus than the nulliparous uteri (n=8). No significant differences were found in FA (0.42 ± 0.02, 0.41 ± 0.02; p=0.25) and ADC (1.82 ± 0.18 × 10(-3) mm(2) s(-1), 1.93 ± 0.25 × 10(-3) mm(2) s(-1); p=0.20).

CONCLUSION

Fibre architecture of the human uterus can be depicted in vivo using 3 T MR-DTI. Advances in knowledge 3 T MR-DTI can help to provide an in vivo insight of uterine anatomy non-invasively, especially in females with previous Caesarean surgery, in order to provide better management of subsequent deliveries.

Derivation of a finite-element model of lingual deformation during swallowing from the mechanics of mesoscale myofiber tracts obtained by MRI

To demonstrate the relationship between lingual myoarchitecture and mechanics during swallowing, we performed a finite-element (FE) simulation of lingual deformation employing mesh aligned with the vector coordinates of myofiber tracts obtained by diffusion tensor imaging with tractography in humans. Material properties of individual elements were depicted in terms of Hill's three-component phenomenological model, assuming that the FE mesh was composed of anisotropic muscle and isotropic connective tissue. Moreover, the mechanical model accounted for elastic constraints by passive and active elements from the superior and inferior directions and the effect of out-of-plane muscles and connective tissue. Passive bolus effects were negligible. Myofiber tract activation was simulated over 500 ms in 1-ms steps following lingual tip association with the hard palate and incorporated specifically the accommodative and propulsive phases of the swallow. Examining the displacement field, active and passive muscle stress, elemental stretch, and strain rate relative to changes of global shape, we demonstrate that lingual reconfiguration during these swallow phases is characterized by (in sequence) the following: 1) lingual tip elevation and shortening in the anterior-posterior direction; 2) inferior displacement related to hyoglossus contraction at its inferior-most position; and 3) dominant clockwise rotation related to regional contraction of the genioglossus and contraction of the hyoglossus following anterior displacement. These simulations demonstrate that lingual deformation during the indicated phases of swallowing requires temporally patterned activation of intrinsic and extrinsic muscles and delineate a method to ascertain the mechanics of normal and pathological swallowing.

Resolving myoarchitectural disarray in the mouse ventricular wall with diffusion spectrum magnetic resonance imaging

The myoarchitecture of the ventricular wall provides a structural template dictating tissue-scale patterns of mechanical function. We studied whether myofiber tract imaging performed with MR diffusion spectrum imaging (DSI) tractography has the capacity to resolve abnormalities of ventricular myoarchitecture in a model of congenital hypertrophic cardiomyopathy (HCM) associated with the ablation of myosin binding protein-C (MyBP-C). Homozygous MyBP-C knockout mice were generated by deletion of exons 3-10 from the endogenous MyBP-C gene. Fiber alignment in the left ventricular wall of wild type mice was depicted through DSI tractography (and confirmed by multi-slice two-photon microscopy) as a set of helical structures whose angles display a continuous transition from negative in the subepicardium to positive in the subendocardium. In contrast, the hearts obtained from the MyBP-C knockouts displayed substantial myoarchitectural disarray, characterized by a loss of voxel-to-voxel orientational coherence for fibers principally located in the mid-myocardium-subendocardium and impairment of the transmural progression of helix angles. These results substantiate the use of DSI tractography in determining myoarchitectural disarray in models of cardiomyopathy and suggest a biological association between myofilament expression, cardiac fiber alignment, and torsional rotation in the setting of congenital HCM.

Diffusion MR tractography of the heart

Histological studies have shown that the myocardium consists of an array of crossing helical fiber tracts. Changes in myocardial fiber architecture occur in ischemic heart disease and heart failure, and can be imaged non-destructively with diffusion-encoded MR. Several diffusion-encoding schemes have been developed, ranging from scalar measurements of mean diffusivity to a 6-dimensional imaging technique known as diffusion spectrum imaging or DSI. The properties of DSI make it particularly suited to the generation of 3-dimensional tractograms of myofiber architecture. In this article we review the physical basis of diffusion-tractography in the myocardium and the attributes of the available techniques, placing particular emphasis on DSI. The application of DSI in ischemic heart disease is reviewed, and the requisites for widespread clinical translation of diffusion MR tractography in the heart are discussed.