Papillary thyroid carcinoma associated with papillary renal neoplasia: genetic linkage analysis of a distinct heritable tumor syndrome

Papillary thyroid carcinoma usually is sporadic, but may occur in a familial form. The complete clinical and pathological phenotype of familial papillary thyroid carcinoma (fPTC) has not been determined, and the susceptibility gene(s) is unknown. We investigated the clinical and pathological characteristics of an unusually large three-generation fPTC kindred to characterize more fully the clinical phenotype. We performed linkage analysis to determine the chromosomal location of a fPTC susceptibility gene. In addition to the known association of fPTC with nodular thyroid disease, we observed the otherwise rare entity of papillary renal neoplasia (PRN) in two kindred members, one affected with PTC and the other an obligate carrier. The multifocality of PRN in one subject adds weight to the likelihood of a true genetic predisposition to PRN. Both genetic linkage and sequence analysis excluded MET, the protooncogene of isolated familial PRN, as the cause of the fPTC/PRN phenotype. A genome-wide screening and an investigation of specific candidate genes demonstrated that the fPTC/PRN phenotype was linked to 1q21. A maximum three-point log of likelihood ratio score of 3.58 was observed for markers D1S2343 and D1S2345 and for markers D1S2343 and D1S305. Critical recombination events limited the region of linkage to approximately 20 cM. A distinct inherited tumor syndrome has been characterized as the familial association of papillary thyroid cancer, nodular thyroid disease, and papillary renal neoplasia. The predisposing gene in a large kindred with this syndrome has been mapped to 1q21.

Methods for quantitative susceptibility and R2* mapping in whole post-mortem brains at 7T applied to amyotrophic lateral sclerosis

Susceptibility weighted magnetic resonance imaging (MRI) is sensitive to the local concentration of iron and myelin. Here, we describe a robust image processing pipeline for quantitative susceptibility mapping (QSM) and R2* mapping of fixed post-mortem, whole-brain data. Using this pipeline, we compare the resulting quantitative maps in brains from patients with amyotrophic lateral sclerosis (ALS) and controls, with validation against iron and myelin histology. Twelve post-mortem brains were scanned with a multi-echo gradient echo sequence at 7T, from which susceptibility and R2* maps were generated. Semi-quantitative histological analysis for ferritin (the principal iron storage protein) and myelin proteolipid protein was performed in the primary motor, anterior cingulate and visual cortices. Magnetic susceptibility and R2* values in primary motor cortex were higher in ALS compared to control brains. Magnetic susceptibility and R2* showed positive correlations with both myelin and ferritin estimates from histology. Four out of nine ALS brains exhibited clearly visible hyperintense susceptibility and R2* values in the primary motor cortex. Our results demonstrate the potential for MRI-histology studies in whole, fixed post-mortem brains to investigate the biophysical source of susceptibility weighted MRI signals in neurodegenerative diseases like ALS.

Dissecting the pathobiology of altered MRI signal in amyotrophic lateral sclerosis: A post mortem whole brain sampling strategy for the integration of ultra-high-field MRI and quantitative neuropathology

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a clinically and histopathologically heterogeneous neurodegenerative disorder, in which therapy is hindered by the rapid progression of disease and lack of biomarkers. Magnetic resonance imaging (MRI) has demonstrated its potential for detecting the pathological signature and tracking disease progression in ALS. However, the microstructural and molecular pathological substrate is poorly understood and generally defined histologically. One route to understanding and validating the pathophysiological correlates of MRI signal changes in ALS is to directly compare MRI to histology in post mortem human brains.

RESULTS

The article delineates a universal whole brain sampling strategy of pathologically relevant grey matter (cortical and subcortical) and white matter tracts of interest suitable for histological evaluation and direct correlation with MRI. A standardised systematic sampling strategy that was compatible with co-registration of images across modalities was established for regions representing phosphorylated 43-kDa TAR DNA-binding protein (pTDP-43) patterns that were topographically recognisable with defined neuroanatomical landmarks. Moreover, tractography-guided sampling facilitated accurate delineation of white matter tracts of interest. A digital photography pipeline at various stages of sampling and histological processing was established to account for structural deformations that might impact alignment and registration of histological images to MRI volumes. Combined with quantitative digital histology image analysis, the proposed sampling strategy is suitable for routine implementation in a high-throughput manner for acquisition of large-scale histology datasets. Proof of concept was determined in the spinal cord of an ALS patient where multiple MRI modalities (T1, T2, FA and MD) demonstrated sensitivity to axonal degeneration and associated heightened inflammatory changes in the lateral corticospinal tract. Furthermore, qualitative comparison of R2* and susceptibility maps in the motor cortex of 2 ALS patients demonstrated varying degrees of hyperintense signal changes compared to a control. Upon histological evaluation of the same region, intensity of signal changes in both modalities appeared to correspond primarily to the degree of microglial activation.

CONCLUSION

The proposed post mortem whole brain sampling methodology enables the accurate intraindividual study of pathological propagation and comparison with quantitative MRI data, to more fully understand the relationship of imaging signal changes with underlying pathophysiology in ALS.

Longitudinal connections and the organization of the temporal cortex in macaques, great apes, and humans

The temporal association cortex is considered a primate specialization and is involved in complex behaviors, with some, such as language, particularly characteristic of humans. The emergence of these behaviors has been linked to major differences in temporal lobe white matter in humans compared with monkeys. It is unknown, however, how the organization of the temporal lobe differs across several anthropoid primates. Therefore, we systematically compared the organization of the major temporal lobe white matter tracts in the human, gorilla, and chimpanzee great apes and in the macaque monkey. We show that humans and great apes, in particular the chimpanzee, exhibit an expanded and more complex occipital-temporal white matter system; additionally, in humans, the invasion of dorsal tracts into the temporal lobe provides a further specialization. We demonstrate the reorganization of different tracts along the primate evolutionary tree, including distinctive connectivity of human temporal gray matter.

Phenotypic and genetic associations of quantitative magnetic susceptibility in UK Biobank brain imaging

A key aim in epidemiological neuroscience is identification of markers to assess brain health and monitor therapeutic interventions. Quantitative susceptibility mapping (QSM) is an emerging magnetic resonance imaging technique that measures tissue magnetic susceptibility and has been shown to detect pathological changes in tissue iron, myelin and calcification. We present an open resource of QSM-based imaging measures of multiple brain structures in 35,273 individuals from the UK Biobank prospective epidemiological study. We identify statistically significant associations of 251 phenotypes with magnetic susceptibility that include body iron, disease, diet and alcohol consumption. Genome-wide associations relate magnetic susceptibility to 76 replicating clusters of genetic variants with biological functions involving iron, calcium, myelin and extracellular matrix. These patterns of associations include relationships that are unique to QSM, in particular being complementary to T2* signal decay time measures. These new imaging phenotypes are being integrated into the core UK Biobank measures provided to researchers worldwide, creating the potential to discover new, non-invasive markers of brain health.

Familial papillary thyroid carcinoma is genetically distinct from familial adenomatous polyposis coli

Familial papillary thyroid carcinoma (fPTC) is an inherited tumor syndrome characterized by isolated papillary thyroid carcinoma (PTC) in affected subjects. Its etiology is unknown. Large multigeneration families with PTC are very rare, and therefore, modern genetic linkage studies have not been applied extensively to this disorder. Familial adenomatous polyposis coli (FAP) is an inherited tumor syndrome enriched in PTC. FAP is caused by germline mutations of the adenomatous polyposis coli (APC) gene that is located in the 5q21 region. It is not known if fPTC is a phenotypic variant of FAP, or if it is a genetically distinct disorder. We report a large 3-generation fPTC kindred and use linkage analysis to test the hypothesis that fPTC and FAP are genetically distinct. In this kindred there are 25 living informative subjects; 5 have PTC, and 1 is an obligate carrier. Inheritance is autosomal dominant with incomplete penetrance. There is vertical transmission, multifocal disease, an average age of onset of 36 years, and 1 subject has colon cancer. The probability is approximately 1 in 2 billion against the clustering of 5 sporadic PTC cases in this kindred. To test for linkage to the APC gene we used 2 highly polymorphic markers, D5S656 and D5S421, which are located within a maximum distance of 1.7 megabase (Mb) of the APC gene and within an estimated genetic region of less than 1 centimorgan (cM) from each other. After polymerase chain reaction (PCR) amplification 18 family members were genotyped. Construction and inspection of haplotypes showed that the affected subjects do not share the same allelic composition. Using a penetrance ratio of 75%, linkage was excluded at 2 cM and 3 cM on both sides of D5S656 and D5S421, respectively. The combined haplotype of these 2 markers provided an exclusion region of 4 cM. We conclude that fPTC is genetically distinct from FAP.

Clinical experience with transdermal clonidine in African-American and Hispanic-American patients with hypertension: evaluation from a 12-week prospective, open-label clinical trial in community-based clinics

The objective of this study was to assess the efficacy and tolerability of transdermal clonidine in inner-city African-American and Hispanic-American patients with essential hypertension. A multiclinic open-label, prospective trial for 12 weeks was used. Dose titration was based on office blood pressure (BP) measurements of > 140/90 mm Hg. Clinical sites were community-based primary care centers. Untreated and treated hypertensive patients whose diastolic BP exceeded 90 mm Hg were administered transdermal clonidine at 0.1 mg or 0.2 mg delivery daily. The drug was titrated after 1 month if diastolic BP was greater than 90 mm Hg. At 12 weeks of treatment, change in blood pressure from baseline as well as adverse effects and patient satisfaction were assessed. A total of 357 patients entered the treatment phase of the study, and 315 patients (244 African-Americans, 67 Hispanic-Americans) had evaluable data. Transdermal clonidine significantly (P <.001) lowered BP in all patients by 15.7/12.8 +/- 18.1/9.6 mm Hg, and heart rate was reduced by 3 +/- 9 beats/min (P <.001). There were no differences in BP reduction according to race and ethnicity, gender, or age. The most common adverse effects were pruritus or discomfort at the patch site, dizziness, dry mouth, and fatigue. Eleven percent of the patients discontinued treatment because of one of these adverse effects. A large proportion of patients (67%) reported that transdermal clonidine was more convenient to use than oral therapy. Transdermal clonidine, alone or in combination with other antihypertensive therapies, significantly lowered BP and heart rate in inner-city hypertensive patients. The drug was generally well tolerated, with 89% of the patients remaining in the trial. Patient acceptability was high with the once-weekly treatment, which is an important feature for this particular hypertensive population.

Preserved extrastriate visual network in a monkey with substantial, naturally occurring damage to primary visual cortex

Lesions of primary visual cortex (V1) lead to loss of conscious visual perception with significant impact on human patients. Understanding the neural consequences of such damage may aid the development of rehabilitation methods. In this rare case of a Rhesus macaque (monkey S), likely born without V1, the animal's in-group behaviour was unremarkable, but visual task training was impaired. With multi-modal magnetic resonance imaging, visual structures outside of the lesion appeared normal. Visual stimulation under anaesthesia with checkerboards activated lateral geniculate nucleus of monkey S, while full-field moving dots activated pulvinar. Visual cortical activation was sparse but included face patches. Consistently across lesion and control monkeys, functional connectivity analysis revealed an intact network of bilateral dorsal visual areas temporally correlated with V5/MT activation, even without V1. Despite robust subcortical responses to visual stimulation, we found little evidence for strengthened subcortical input to V5/MT supporting residual visual function or blindsight-like phenomena.

The Digital Brain Bank, an open access platform for post-mortem imaging datasets

Post-mortem magnetic resonance imaging (MRI) provides the opportunity to acquire high-resolution datasets to investigate neuroanatomy and validate the origins of image contrast through microscopy comparisons. We introduce the Digital Brain Bank (open.win.ox.ac.uk/DigitalBrainBank), a data release platform providing open access to curated, multimodal post-mortem neuroimaging datasets. Datasets span three themes-Digital Neuroanatomist: datasets for detailed neuroanatomical investigations; Digital Brain Zoo: datasets for comparative neuroanatomy; and Digital Pathologist: datasets for neuropathology investigations. The first Digital Brain Bank data release includes 21 distinctive whole-brain diffusion MRI datasets for structural connectivity investigations, alongside microscopy and complementary MRI modalities. This includes one of the highest-resolution whole-brain human diffusion MRI datasets ever acquired, whole-brain diffusion MRI in fourteen nonhuman primate species, and one of the largest post-mortem whole-brain cohort imaging studies in neurodegeneration. The Digital Brain Bank is the culmination of our lab's investment into post-mortem MRI methodology and MRI-microscopy analysis techniques. This manuscript provides a detailed overview of our work with post-mortem imaging to date, including the development of diffusion MRI methods to image large post-mortem samples, including whole, human brains. Taken together, the Digital Brain Bank provides cross-scale, cross-species datasets facilitating the incorporation of post-mortem data into neuroimaging studies.

Diffusion MRI data, sulcal anatomy, and tractography for eight species from the Primate Brain Bank

Large-scale comparative neuroscience requires data from many species and, ideally, at multiple levels of description. Here, we contribute to this endeavor by presenting diffusion and structural MRI data from eight primate species that have not or rarely been described in the literature. The selected samples from the Primate Brain Bank cover a prosimian, New and Old World monkeys, and a great ape. We present preliminary labelling of the cortical sulci and tractography of the optic radiation, dorsal part of the cingulum bundle, and dorsal parietal-frontal and ventral temporal-frontal longitudinal white matter tracts. Both dorsal and ventral association fiber systems could be observed in all samples, with the dorsal tracts occupying much less relative volume in the prosimian than in other species. We discuss the results in the context of known primate specializations and present hypotheses for further research. All data and results presented here are available online as a resource for the scientific community.

Frequency difference mapping applied to the corpus callosum at 7T

Purpose

Frequency difference mapping (FDM) is a phase processing technique which characterizes the nonlinear temporal evolution of the phase of gradient echo (GE) signals. Here, a novel FDM‐processing algorithm is introduced, which is shown to reveal information about white matter microstructure. Unlike some other phase‐processing techniques, the FDM algorithm presented here does not require the use of phase unwrapping or sophisticated image processing. It uses a series of scaled complex divisions to unwrap phase and remove background fields.

Methods

Ten healthy subjects underwent a series of single‐slice, sagittal multi‐echo GE scans at 7T with the slice positioned at the midline. Phase data were processed with the novel FDM algorithm, and the temporal evolution of the magnitude signal and frequency difference was examined in 5 regions of the corpus callosum (CC; genu, anterior body, middle body, posterior body, and splenium).

Results

Consistent frequency difference contrast relative to surrounding tissue was observed in all subjects in the CC and in other white matter regions where the nerve fibers run perpendicular to B0, such as the superior cerebellar peduncle. Examination of the frequency difference curves shows distinct variations over the CC, with the genu and splenium displaying larger frequency differences than the other regions (in addition to a faster decay of signal magnitude).

Conclusion

The novel FDM algorithm presented here yields images sensitive to tissue microstructure and microstructural differences over the CC in a simple manner, without the requirement for phase unwrapping or sophisticated image processing.

Brain gyrification in wild and domestic canids: Has domestication changed the gyrification index in domestic dogs?

Over the last 15 years, research on canid cognition has revealed that domestic dogs possess a surprising array of complex sociocognitive skills pointing to the possibility that the domestication process might have uniquely altered their brains; however, we know very little about how evolutionary processes (natural or artificial) might have modified underlying neural structure to support species-specific behaviors. Evaluating the degree of cortical folding (i.e., gyrification) within canids may prove useful, as this parameter is linked to functional variation of the cerebral cortex. Using quantitative magnetic resonance imaging to investigate the impact of domestication on the canine cortical surface, we compared the gyrification index (GI) in 19 carnivore species, including six wild canid and 13 domestic dog individuals. We also explored correlations between global and local GI with brain mass, cortical thickness, white and gray matter volume and surface area. Our results indicated that GI values for domestic dogs are largely consistent with what would be expected for a canid of their given brain mass, although more variable than that observed in wild canids. We also found that GI in canids is positively correlated with cortical surface area, cortical thickness and total cortical gray matter volumes. While we found no evidence of global differences in GI between domestic and wild canids, certain regional differences in gyrification were observed.

Modeling an equivalent b-value in diffusion-weighted steady-state free precession

PURPOSE

Diffusion-weighted steady-state free precession (DW-SSFP) is shown to provide a means to probe non-Gaussian diffusion through manipulation of the flip angle. A framework is presented to define an effective b-value in DW-SSFP.

THEORY

The DW-SSFP signal is a summation of coherence pathways with different b-values. The relative contribution of each pathway is dictated by the flip angle. This leads to an apparent diffusion coefficient (ADC) estimate that depends on the flip angle in non-Gaussian diffusion regimes. By acquiring DW-SSFP data at multiple flip angles and modeling the variation in ADC for a given form of non-Gaussianity, the ADC can be estimated at a well-defined effective b-value.

METHODS

A gamma distribution is used to model non-Gaussian diffusion, embedded in the Buxton signal model for DW-SSFP. Monte-Carlo simulations of non-Gaussian diffusion in DW-SSFP and diffusion-weighted spin-echo sequences are used to verify the proposed framework. Dependence of ADC on flip angle in DW-SSFP is verified with experimental measurements in a whole, human postmortem brain.

RESULTS

Monte-Carlo simulations reveal excellent agreement between ADCs estimated with diffusion-weighted spin-echo and the proposed framework. Experimental ADC estimates vary as a function of flip angle over the corpus callosum of the postmortem brain, estimating the mean and standard deviation of the gamma distribution as 1.50 · 10 - 4  mm2 /s and 2.10 · 10 - 4  mm2 /s.

CONCLUSION

DW-SSFP can be used to investigate non-Gaussian diffusion by varying the flip angle. By fitting a model of non-Gaussian diffusion, the ADC in DW-SSFP can be estimated at an effective b-value, comparable to more conventional diffusion sequences.

Tensor image registration library: Deformable registration of stand-alone histology images to whole-brain post-mortem MRI data

BACKGROUND

Accurate registration between microscopy and MRI data is necessary for validating imaging biomarkers against neuropathology, and to disentangle complex signal dependencies in microstructural MRI. Existing registration methods often rely on serial histological sampling or significant manual input, providing limited scope to work with a large number of stand-alone histology sections. Here we present a customisable pipeline to assist the registration of stand-alone histology sections to whole-brain MRI data.

METHODS

Our pipeline registers stained histology sections to whole-brain post-mortem MRI in 4 stages, with the help of two photographic intermediaries: a block face image (to undistort histology sections) and coronal brain slab photographs (to insert them into MRI space). Each registration stage is implemented as a configurable stand-alone Python script using our novel platform, Tensor Image Registration Library (TIRL), which provides flexibility for wider adaptation. We report our experience of registering 87 PLP-stained histology sections from 14 subjects and perform various experiments to assess the accuracy and robustness of each stage of the pipeline.

RESULTS

All 87 histology sections were successfully registered to MRI. Histology-to-block registration (Stage 1) achieved 0.2-0.4 mm accuracy, better than commonly used existing methods. Block-to-slice matching (Stage 2) showed great robustness in automatically identifying and inserting small tissue blocks into whole brain slices with 0.2 mm accuracy. Simulations demonstrated sub-voxel level accuracy (0.13 mm) of the slice-to-volume registration (Stage 3) algorithm, which was observed in over 200 actual brain slice registrations, compensating 3D slice deformations up to 6.5 mm. Stage 4 combined the previous stages and generated refined pixelwise aligned multi-modal histology-MRI stacks.

CONCLUSIONS

Our open-source pipeline provides robust automation tools for registering stand-alone histology sections to MRI data with sub-voxel level precision, and the underlying framework makes it readily adaptable to a diverse range of microscopy-MRI studies.

A method to remove the influence of fixative concentration on postmortem T2 maps using a kinetic tensor model

Formalin fixation has been shown to substantially reduce T₂ estimates, primarily driven by the presence of fixative in tissue. Prior to scanning, post‐mortem samples are often placed into a fluid that has more favourable imaging properties. This study investigates whether there is evidence for a change in T₂ in regions close to the tissue surface due to fixative outflux into this surrounding fluid. Furthermore, we investigate whether a simulated spatial map of fixative concentration can be used as a confound regressor to reduce T₂ inhomogeneity. To achieve this, T₂ maps and diffusion tensor estimates were obtained in 14 whole, formalin‐fixed post‐mortem brains placed in Fluorinert approximately 48 hr prior to scanning. Seven brains were fixed with 10% formalin and seven brains were fixed with 10% neutral buffered formalin (NBF). Fixative outflux was modelled using a proposed kinetic tensor (KT) model, which incorporates voxelwise diffusion tensor estimates to account for diffusion anisotropy and tissue‐specific diffusion coefficients. Brains fixed with 10% NBF revealed a spatial T₂ pattern consistent with modelled fixative outflux. Confound regression of fixative concentration reduced T₂ inhomogeneity across both white and grey matter, with the greatest reduction attributed to the KT model versus simpler models of fixative outflux. No such effect was observed in brains fixed with 10% formalin. Correlations between the transverse relaxation rate R₂ and ferritin/myelin proteolipid protein (PLP) histology lead to an increased similarity for the relationship between R₂ and PLP for the two fixative types after KT correction.

Use of multi-flip angle measurements to account for transmit inhomogeneity and non-Gaussian diffusion in DW-SSFP

Diffusion-weighted steady-state free precession (DW-SSFP) is an SNR-efficient diffusion imaging method. The improved SNR and resolution available at ultra-high field has motivated its use at 7T. However, these data tend to have severe B1 inhomogeneity, leading not only to spatially varying SNR, but also to spatially varying diffusivity estimates, confounding comparisons both between and within datasets. This study proposes the acquisition of DW-SSFP data at two-flip angles in combination with explicit modelling of non-Gaussian diffusion to address B1 inhomogeneity at 7T. Data were acquired from five fixed whole human post-mortem brains with a pair of flip angles that jointly optimize the diffusion contrast-to-noise (CNR) across the brain. We compared one- and two-flip angle DW-SSFP data using a tensor model that incorporates the full DW-SSFP Buxton signal, in addition to tractography performed over the cingulum bundle and pre-frontal cortex using a ball & sticks model. The two-flip angle DW-SSFP data produced angular uncertainty and tractography estimates close to the CNR optimal regions in the single-flip angle datasets. The two-flip angle tensor estimates were subsequently fitted using a modified DW-SSFP signal model that incorporates a gamma distribution of diffusivities. This allowed us to generate tensor maps at a single effective b-value yielding more consistent SNR across tissue, in addition to eliminating the B1 dependence on diffusion coefficients and orientation maps. Our proposed approach will allow the use of DW-SSFP at 7T to derive diffusivity estimates that have greater interpretability, both within a single dataset and between experiments.

Why every lab needs a handbook

A lab handbook is a flexible document that outlines the ethos of a research lab or group. A good handbook will outline the different roles within the lab, explain what is expected of all lab members, provide an overview of the culture the lab aims to create, and describe how the lab supports its members so that they can develop as researchers. Here we describe how we wrote a lab handbook for a large research group, and provide resources to help other labs write their own handbooks.

Considerations and recommendations from the ISMRM diffusion study group for preclinical diffusion MRI: Part 1: In vivo small-animal imaging

Small-animal diffusion MRI (dMRI) has been used for methodological development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. The steps from animal setup and monitoring, to acquisition, analysis, and interpretation are complex, with many decisions that may ultimately affect what questions can be answered using the resultant data. This work aims to present selected considerations and recommendations from the diffusion community on best practices for preclinical dMRI of in vivo animals. We describe the general considerations and foundational knowledge that must be considered when designing experiments. We briefly describe differences in animal species and disease models and discuss why some may be more or less appropriate for different studies. We, then, give recommendations for in vivo acquisition protocols, including decisions on hardware, animal preparation, and imaging sequences, followed by advice for data processing including preprocessing, model-fitting, and tractography. Finally, we provide an online resource that lists publicly available preclinical dMRI datasets and software packages to promote responsible and reproducible research. In each section, we attempt to provide guides and recommendations, but also highlight areas for which no guidelines exist (and why), and where future work should focus. Although we mainly cover the central nervous system (on which most preclinical dMRI studies are focused), we also provide, where possible and applicable, recommendations for other organs of interest. An overarching goal is to enhance the rigor and reproducibility of small animal dMRI acquisitions and analyses, and thereby advance biomedical knowledge.

Investigating time-independent and time-dependent diffusion phenomena using steady-state diffusion MRI

Diffusion MRI is a leading method to non-invasively characterise brain tissue microstructure across multiple domains and scales. Diffusion-weighted steady-state free precession (DW-SSFP) is an established imaging sequence for post-mortem MRI, addressing the challenging imaging environment of fixed tissue with short T2 and low diffusivities. However, a current limitation of DW-SSFP is signal interpretation: it is not clear what diffusion 'regime' the sequence probes and therefore its potential to characterise tissue microstructure. Building on Extended Phase Graphs (EPG), I establish two alternative representations of the DW-SSFP signal in terms of (1) conventional b-values (time-independent diffusion) and (2) encoding power-spectra (time-dependent diffusion). The proposed representations provide insights into how different parameter regimes and gradient waveforms impact the diffusion sensitivity of DW-SSFP. I subsequently introduce an approach to incorporate existing biophysical models into DW-SSFP without the requirement of extensive derivations, with time dependence estimated via a Gaussian phase approximation representation of the DW-SSFP signal. Investigations incorporating free-diffusion and tissue-relevant microscopic restrictions (cylinder of varying radius) give excellent agreement to complementary analytical models and Monte Carlo simulations. Experimentally, the time-independent representation is used to derive Tensor and proof-of-principle NODDI estimates in a whole human post-mortem brain. A final SNR-efficiency investigation demonstrates the theoretical potential of DW-SSFP for ultra-high field microstructural imaging.

Considerations and recommendations from the ISMRM Diffusion Study Group for preclinical diffusion MRI: Part 3-Ex vivo imaging: Data processing, comparisons with microscopy, and tractography

Preclinical diffusion MRI (dMRI) has proven value in methods development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. While dMRI enables in vivo non-invasive characterization of tissue, ex vivo dMRI is increasingly being used to probe tissue microstructure and brain connectivity. Ex vivo dMRI has several experimental advantages that facilitate high spatial resolution and high SNR images, cutting-edge diffusion contrasts, and direct comparison with histological data as a methodological validation. However, there are a number of considerations that must be made when performing ex vivo experiments. The steps from tissue preparation, image acquisition and processing, and interpretation of results are complex, with many decisions that not only differ dramatically from in vivo imaging of small animals, but ultimately affect what questions can be answered using the data. This work concludes a three-part series of recommendations and considerations for preclinical dMRI. Herein, we describe best practices for dMRI of ex vivo tissue, with a focus on image pre-processing, data processing, and comparisons with microscopy. In each section, we attempt to provide guidelines and recommendations but also highlight areas for which no guidelines exist (and why), and where future work should lie. We end by providing guidelines on code sharing and data sharing and point toward open-source software and databases specific to small animal and ex vivo imaging.

Both noise-floor and tissue compartment difference in diffusivity contribute to FA dependence on b-value in diffusion MRI

Noninvasive diffusion magnetic resonance imaging (dMRI) has been widely employed in both clinical and research settings to investigate brain tissue microstructure. Despite the evidence that dMRI-derived fractional anisotropy (FA) correlates with white matter properties, the metric is not specific. Recent studies have reported that FA is dependent on the b-value, and its origin has primarily been attributed to either the influence of microstructure or the noise-floor effect. A systematic investigation into the inter-relationship of these two effects is however still lacking. This study aims to quantify contributions of the reported differences in intra- and extra-neurite diffusivity to the observed changes in FA, in addition to the noise in measurements. We used in-vivo and post-mortem human brain imaging, as well as numerical simulations and histological validation, for this purpose. Our investigations reveal that the percentage difference of FA between b-values (pdFA) has significant positive associations with neurite density index (NDI), which is derived from in-vivo neurite orientation dispersion and density imaging (NODDI), or Bielschowsky's silver impregnation (BIEL) staining sections of fixed post-mortem human brain samples. Furthermore, such an association is found to be varied with Signal-to-Noise Ratio (SNR) level, indicating a nonlinear interaction effect between tissue microstructure and noise. Finally, a multicompartment model simulation revealed that these findings can be driven by differing diffusivities of intra- and extra-neurite compartments in tissue, with the noise-floor further amplifying the effect. In conclusion, both the differences in intra- and extra-neurite diffusivity and noise-floor effects significantly contribute to the FA difference associated with the b-value.

Automated detection of cerebral microbleeds on MR images using knowledge distillation framework

Introduction

Cerebral microbleeds (CMBs) are associated with white matter damage, and various neurodegenerative and cerebrovascular diseases. CMBs occur as small, circular hypointense lesions on T2*-weighted gradient recalled echo (GRE) and susceptibility-weighted imaging (SWI) images, and hyperintense on quantitative susceptibility mapping (QSM) images due to their paramagnetic nature. Accurate automated detection of CMBs would help to determine quantitative imaging biomarkers (e.g., CMB count) on large datasets. In this work, we propose a fully automated, deep learning-based, 3-step algorithm, using structural and anatomical properties of CMBs from any single input image modality (e.g., GRE/SWI/QSM) for their accurate detections.

Methods

In our method, the first step consists of an initial candidate detection step that detects CMBs with high sensitivity. In the second step, candidate discrimination step is performed using a knowledge distillation framework, with a multi-tasking teacher network that guides the student network to classify CMB and non-CMB instances in an offline manner. Finally, a morphological clean-up step further reduces false positives using anatomical constraints. We used four datasets consisting of different modalities specified above, acquired using various protocols and with a variety of pathological and demographic characteristics.

Results

On cross-validation within datasets, our method achieved a cluster-wise true positive rate (TPR) of over 90% with an average of <2 false positives per subject. The knowledge distillation framework improves the cluster-wise TPR of the student model by 15%. Our method is flexible in terms of the input modality and provides comparable cluster-wise TPR and better cluster-wise precision compared to existing state-of-the-art methods. When evaluating across different datasets, our method showed good generalizability with a cluster-wise TPR >80 % with different modalities. The python implementation of the proposed method is openly available.

Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations

PURPOSE

Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination.

THEORY AND METHODS

One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation.

RESULTS

The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain).

CONCLUSION

This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.

An automated pipeline for extracting histological stain area fraction for voxelwise quantitative MRI-histology comparisons

The acquisition of MRI and histology in the same post-mortem tissue sample enables direct correlation between MRI and histologically-derived parameters. However, there still lacks a standardised automated pipeline to process histology data, with most studies relying on manual intervention. Here, we introduce an automated pipeline to extract a quantitative histological measure for staining density (stain area fraction, SAF) from multiple immunohistochemical (IHC) stains. The pipeline is designed to directly address key IHC artefacts related to tissue staining and slide digitisation. Here, the pipeline was applied to post-mortem human brain data from multiple subjects, relating MRI parameters (FA, MD, RD, AD, R2*, R1) to IHC slides stained for myelin, neurofilaments, microglia and activated microglia. Utilising high-quality MRI-histology co-registrations, we then performed whole-slide voxelwise comparisons (simple correlations, partial correlations and multiple regression analyses) between multimodal MRI- and IHC-derived parameters. The pipeline was found to be reproducible, robust to artefacts and generalisable across multiple IHC stains. Our partial correlation results suggest that some simple MRI-SAF correlations should be interpreted with caution, due to the co-localisation of other tissue features (e.g., myelin and neurofilaments). Further, we find activated microglia-a generic biomarker of inflammation-to consistently be the strongest predictor of high DTI FA and low RD, which may suggest sensitivity of diffusion MRI to aspects of neuroinflammation related to microglial activation, even after accounting for other microstructural changes (demyelination, axonal loss and general microglia infiltration). Together, these results show the utility of this approach in carefully curating IHC data and performing multimodal analyses to better understand microstructural relationships with MRI.