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.

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.

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.

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.

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.

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.

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.

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 T2 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 T2 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 T2 inhomogeneity. To achieve this, T2 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 T2 pattern consistent with modelled fixative outflux. Confound regression of fixative concentration reduced T2 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 R2 and ferritin/myelin proteolipid protein (PLP) histology lead to an increased similarity for the relationship between R2 and PLP for the two fixative types after KT correction.

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.

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.

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.

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.

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.

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.

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.

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.

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.