Nonlinear Gradient Field Estimation in Diffusion MRI Tensor Simulation.

Gradient nonlinearities not only induce spatial distortion in magnetic resonance imaging (MRI), but also introduce discrepancies between intended and acquired diffusion sensitization in diffusion weighted (DW) MRI. Advances in scanner performance have increased the importance of correcting gradient nonlinearities. The most common approaches for gradient nonlinear field estimations rely on phantom calibration field maps which are not always feasible, especially on retrospective data. Here, we derive a quadratic minimization problem for the complete gradient nonlinear field (L(r)). This approach starts with corrupt diffusion signal and estimates the L(r) in two scenarios: (1) the true diffusion tensor known and (2) the true diffusion tensor unknown (i.e., diffusion tensor is estimated). We show the validity of this mathematical approach, both theoretically and through tensor simulation. The estimated field is assessed through diffusion tensor metrics: mean diffusivity (MD), fractional anisotropy (FA), and principal eigenvector (V1). In simulation with 300 diffusion tensors, the study shows the mathematical model is not ill-posed and remains stable. We find when the true diffusion tensor is known (1) the change in determinant of the estimated L(r) field and the true field is near zero and (2) the median difference in estimated L(r) corrected diffusion metrics to true values is near zero. We find the results of L(r) estimation are dependent on the level of L(r) corruption. This work provides an approach to estimate gradient field without the need for additional calibration scans. To the best of our knowledge, the mathematical derivation presented here is novel.

Assessment of Subject Head Motion in Diffusion MRI.

Subject head motion during the acquisition of diffusion-weighted imaging (DWI) of the brain induces artifacts and affects image quality. Information about the frequency and extent of motion could reveal which aspects of motion correction are most necessary. Therefore, we investigate the extent of translation and rotation among participants, and how the motion changes during the scan acquisition. We analyze 5,380 DWI scans from 1,034 participants. We measure the rotations and translations in the sagittal, coronal and transverse planes needed to align the volumes to the first and previous volumes, as well as the displacement. The different types of motion are compared with each other and compared over time. The largest rotation (per minute) is around the right - left axis (median 0.378 °/min, range 0.000 - 11.466°) and the largest translation (per minute) is along the anterior - posterior axis (median 1.867 mm/min, range 0.000 - 10.944 mm). We additionally observe that spikes in movement occur at the beginning of the scan, particularly in anterior - posterior translation. The results show that all scans are affected by subtle head motion, which may impact subsequent image analysis.

Tractography with T1-weighted MRI and associated anatomical constraints on clinical quality diffusion MRI.

Diffusion MRI (dMRI) streamline tractography, the gold-standard for in vivo estimation of white matter (WM) pathways in the brain, has long been considered as a product of WM microstructure. However, recent advances in tractography demonstrated that convolutional recurrent neural networks (CoRNN) trained with a teacher-student framework have the ability to learn to propagate streamlines directly from T1 and anatomical context. Training for this network has previously relied on high resolution dMRI. In this paper, we generalize the training mechanism to traditional clinical resolution data, which allows generalizability across sensitive and susceptible study populations. We train CoRNN on a small subset of the Baltimore Longitudinal Study of Aging (BLSA), which better resembles clinical scans. We define a metric, termed the epsilon ball seeding method, to compare T1 tractography and traditional diffusion tractography at the streamline level. We show that under this metric T1 tractography generated by CoRNN reproduces diffusion tractography with approximately three millimeters of error.

Functional correlation tensors in brain white matter and the effects of normal aging.

Resting state correlations between blood oxygenation level dependent (BOLD) MRI signals from voxels in white matter (WM) are demonstrably anisotropic, so that functional correlation tensors (FCT) may be used to quantify the underlying microstructure of BOLD effects in WM tracts. However, the overall spatial distribution of FCTs and their metrics in specific populations has not yet been established, and the factors that affect their precise arrangements remain unclear. Changes in WM occur with normal aging, and these may be expected to affect FCTs. We hypothesized that FCTs exhibit a characteristic spatial pattern and may show systematic changes with aging or other factors. Here we report our analyses of the FCT characteristics of fMRI images of a large cohort of 461 cognitively normal subjects (190 females, 271 males) sourced from the Open Access Series of Imaging Studies (OASIS), with age distributions of 42 y/o - 95 y/o. Group averages and statistics of FCT indices, including axial functional correlations, radial functional correlations, mean functional correlations and fractional anisotropy, were quantified in WM bundles defined by the JHU ICBM-DTI-81 WM atlas. In addition, their variations with normal aging were examined. The results revealed a dimorphic distribution of changes in FCT metrics with age, with decreases of the functional correlations in some regions and increases in others. Supplementary analysis revealed that females exhibited significant age effects on a greater number of WM areas, but the interaction between age and sex was not significant. The findings demonstrate the reproducibility of the spatial distribution of FCT metrics and reveal subtle regional changes with age.

Identification of functional white matter networks in BOLD fMRI.

White matter signals in resting state blood oxygen level dependent functional magnetic resonance (BOLD-fMRI) have been largely discounted, yet there is growing evidence that these signals are indicative of brain activity. Understanding how these white matter signals capture function can provide insight into brain physiology. Moreover, functional signals could potentially be used as early markers for neurological changes, such as in Alzheimer's Disease. To investigate white matter brain networks, we leveraged the OASIS-3 dataset to extract white matter signals from resting state BOLD-FMRI data on 711 subjects. The imaging was longitudinal with a total of 2,026 images. Hierarchical clustering was performed to investigate clusters of voxel-level correlations on the timeseries data. The stability of clusters was measured with the average Dice coefficients on two different cross fold validations. The first validated the stability between scans, and the second validated the stability between populations. Functional clusters at hierarchical levels 4, 9, 13, 18, and 24 had local maximum stability, suggesting better clustered white matter. In comparison with JHU-DTI-SS Type-I Atlas defined regions, clusters at lower hierarchical levels identified well-defined anatomical lobes. At higher hierarchical levels, functional clusters mapped motor and memory functional regions, identifying 50.00%, 20.00%, 27.27%, and 35.14% of the frontal, occipital, parietal, and temporal lobe regions respectively.

Harmonized connectome resampling for variance in voxel sizes.

To date, there has been no comprehensive study characterizing the effect of diffusion-weighted magnetic resonance imaging voxel resolution on the resulting connectome for high resolution subject data. Similarity in results improved with higher resolution, even after initial down-sampling. To ensure robust tractography and connectomes, resample data to 1 mm isotropic resolution.

Characterizing patterns of diffusion tensor imaging variance in aging brains.

Purpose: As large analyses merge data across sites, a deeper understanding of variance in statistical assessment across the sources of data becomes critical for valid analyses. Diffusion tensor imaging (DTI) exhibits spatially varying and correlated noise, so care must be taken with distributional assumptions. Here, we characterize the role of physiology, subject compliance, and the interaction of the subject with the scanner in the understanding of DTI variability, as modeled in the spatial variance of derived metrics in homogeneous regions. Approach: We analyze DTI data from 1035 subjects in the Baltimore Longitudinal Study of Aging, with ages ranging from 22.4 to 103 years old. For each subject, up to 12 longitudinal sessions were conducted. We assess the variance of DTI scalars within regions of interest (ROIs) defined by four segmentation methods and investigate the relationships between the variance and covariates, including baseline age, time from the baseline (referred to as "interval"), motion, sex, and whether it is the first scan or the second scan in the session. Results: Covariate effects are heterogeneous and bilaterally symmetric across ROIs. Inter-session interval is positively related ( p ≪ 0.001 ) to FA variance in the cuneus and occipital gyrus, but negatively ( p ≪ 0.001 ) in the caudate nucleus. Males show significantly ( p ≪ 0.001 ) higher FA variance in the right putamen, thalamus, body of the corpus callosum, and cingulate gyrus. In 62 out of 176 ROIs defined by the Eve type-1 atlas, an increase in motion is associated ( p < 0.05 ) with a decrease in FA variance. Head motion increases during the rescan of DTI ( Δ µ = 0.045 mm per volume). Conclusions: The effects of each covariate on DTI variance and their relationships across ROIs are complex. Ultimately, we encourage researchers to include estimates of variance when sharing data and consider models of heteroscedasticity in analysis. This work provides a foundation for study planning to account for regional variations in metric variance.

Empirical assessment of the assumptions of ComBat with diffusion tensor imaging.

Purpose: Diffusion tensor imaging (DTI) is a magnetic resonance imaging technique that provides unique information about white matter microstructure in the brain but is susceptible to confounding effects introduced by scanner or acquisition differences. ComBat is a leading approach for addressing these site biases. However, despite its frequent use for harmonization, ComBat's robustness toward site dissimilarities and overall cohort size have not yet been evaluated in terms of DTI. Approach: As a baseline, we match N=358 participants from two sites to create a "silver standard" that simulates a cohort for multi-site harmonization. Across sites, we harmonize mean fractional anisotropy and mean diffusivity, calculated using participant DTI data, for the regions of interest defined by the JHU EVE-Type III atlas. We bootstrap 10 iterations at 19 levels of total sample size, 10 levels of sample size imbalance between sites, and 6 levels of mean age difference between sites to quantify (i) ßAGE, the linear regression coefficient of the relationship between FA and age; (ii) Î³/f*, the ComBat-estimated site-shift; and (iii) Î´/f*, the ComBat-estimated site-scaling. We characterize the reliability of ComBat by evaluating the root mean squared error in these three metrics and examine if there is a correlation between the reliability of ComBat and a violation of assumptions. Results: ComBat remains well behaved for ßAGE when N>162 and when the mean age difference is less than 4 years. The assumptions of the ComBat model regarding the normality of residual distributions are not violated as the model becomes unstable. Conclusion: Prior to harmonization of DTI data with ComBat, the input cohort should be examined for size and covariate distributions of each site. Direct assessment of residual distributions is less informative on stability than bootstrap analysis. We caution use ComBat of in situations that do not conform to the above thresholds.

Neuronal microstructural changes in the human brain are associated with neurocognitive aging.

Gray matter (GM) alterations play a role in aging-related disorders like Alzheimer's disease and related dementias, yet MRI studies mainly focus on macroscopic changes. Although reliable indicators of atrophy, morphological metrics like cortical thickness lack the sensitivity to detect early changes preceding visible atrophy. Our study aimed at exploring the potential of diffusion MRI in unveiling sensitive markers of cortical and subcortical age-related microstructural changes and assessing their associations with cognitive and behavioral deficits. We leveraged the Human Connectome Project-Aging cohort that included 707 participants (394 female; median age = 58, range = 36-90 years) and applied the powerful mean apparent diffusion propagator model to measure microstructural parameters, along with comprehensive behavioral and cognitive test scores. Both macro- and microstructural GM characteristics were strongly associated with age, with widespread significant microstructural correlations reflective of cellular morphological changes, reduced cellular density, increased extracellular volume, and increased membrane permeability. Importantly, when correlating MRI and cognitive test scores, our findings revealed no link between macrostructural volumetric changes and neurobehavioral performance. However, we found that cellular and extracellular alterations in cortical and subcortical GM regions were associated with neurobehavioral performance. Based on these findings, it is hypothesized that increased microstructural heterogeneity and decreased neurite orientation dispersion precede macrostructural changes, and that they play an important role in subsequent cognitive decline. These alterations are suggested to be early markers of neurocognitive performance that may distinctly aid in identifying the mechanisms underlying phenotypic aging and subsequent age-related functional decline.

Quantification of mediation effects of white matter functional characteristics on cognitive decline in aging.

Cognitive decline with aging involves multifactorial processes, including changes in brain structure and function. This study focuses on the role of white matter functional characteristics, as reflected in blood oxygenation level-dependent signals, in age-related cognitive deterioration. Building on previous research confirming the reproducibility and age-dependence of blood oxygenation level-dependent signals acquired via functional magnetic resonance imaging, we here employ mediation analysis to test if aging affects cognition through white matter blood oxygenation level-dependent signal changes, impacting various cognitive domains and specific white matter regions. We used independent component analysis of resting-state blood oxygenation level-dependent signals to segment white matter into coherent hubs, offering a data-driven view of white matter's functional architecture. Through correlation analysis, we constructed a graph network and derived metrics to quantitatively assess regional functional properties based on resting-state blood oxygenation level-dependent fluctuations. Our analysis identified significant mediators in the age-cognition relationship, indicating that aging differentially influences cognitive functions by altering the functional characteristics of distinct white matter regions. These findings enhance our understanding of the neurobiological basis of cognitive aging, highlighting the critical role of white matter in maintaining cognitive integrity and proposing new approaches to assess interventions targeting cognitive decline in older populations.

DeepN4: Learning N4ITK Bias Field Correction for T1-weighted Images.

T1-weighted (T1w) MRI has low frequency intensity artifacts due to magnetic field inhomogeneities. Removal of these biases in T1w MRI images is a critical preprocessing step to ensure spatially consistent image interpretation. N4ITK bias field correction, the current state-of-the-art, is implemented in such a way that makes it difficult to port between different pipelines and workflows, thus making it hard to reimplement and reproduce results across local, cloud, and edge platforms. Moreover, N4ITK is opaque to optimization before and after its application, meaning that methodological development must work around the inhomogeneity correction step. Given the importance of bias fields correction in structural preprocessing and flexible implementation, we pursue a deep learning approximation / reinterpretation of the N4ITK bias fields correction to create a method which is portable, flexible, and fully differentiable. In this paper, we trained a deep learning network "DeepN4" on eight independent cohorts from 72 different scanners and age ranges with N4ITK-corrected T1w MRI and bias field for supervision in log space. We found that we can closely approximate N4ITK bias fields correction with naïve networks. We evaluate the peak signal to noise ratio (PSNR) in test dataset against the N4ITK corrected images. The median PSNR of corrected images between N4ITK and DeepN4 was 47.96 dB. In addition, we assess the DeepN4 model on eight additional external datasets and show the generalizability of the approach. This study establishes that incompatible N4ITK preprocessing steps can be closely approximated by naïve deep neural networks, facilitating more flexibility. All code and models are released at https://github.com/MASILab/DeepN4 .

Evaluation of Mean Shift, ComBat, and CycleGAN for Harmonizing Brain Connectivity Matrices Across Sites.

Connectivity matrices derived from diffusion MRI (dMRI) provide an interpretable and generalizable way of understanding the human brain connectome. However, dMRI suffers from inter-site and between-scanner variation, which impedes analysis across datasets to improve robustness and reproducibility of results. To evaluate different harmonization approaches on connectivity matrices, we compared graph measures derived from these matrices before and after applying three harmonization techniques: mean shift, ComBat, and CycleGAN. The sample comprises 168 age-matched, sex-matched normal subjects from two studies: the Vanderbilt Memory and Aging Project (VMAP) and the Biomarkers of Cognitive Decline Among Normal Individuals (BIOCARD). First, we plotted the graph measures and used coefficient of variation (CoV) and the Mann-Whitney U test to evaluate different methods' effectiveness in removing site effects on the matrices and the derived graph measures. ComBat effectively eliminated site effects for global efficiency and modularity and outperformed the other two methods. However, all methods exhibited poor performance when harmonizing average betweenness centrality. Second, we tested whether our harmonization methods preserved correlations between age and graph measures. All methods except for CycleGAN in one direction improved correlations between age and global efficiency and between age and modularity from insignificant to significant with p-values less than 0.05.

Identification of functional white matter networks in BOLD fMRI.

White matter signals in resting state blood oxygen level dependent functional magnetic resonance (BOLD-fMRI) have been largely discounted, yet there is growing evidence that these signals are indicative of brain activity. Understanding how these white matter signals capture function can provide insight into brain physiology. Moreover, functional signals could potentially be used as early markers for neurological changes, such as in Alzheimer's Disease. To investigate white matter brain networks, we leveraged the OASIS-3 dataset to extract white matter signals from resting state BOLD-FMRI data on 711 subjects. The imaging was longitudinal with a total of 2,026 images. Hierarchical clustering was performed to investigate clusters of voxel-level correlations on the timeseries data. The stability of clusters was measured with the average Dice coefficients on two different cross fold validations. The first validated the stability between scans, and the second validated the stability between subject populations. Functional clusters at hierarchical levels 4, 9, 13, 18, and 24 had local maximum stability, suggesting better clustered white matter. In comparison with JHU-DTI-SS Type-I Atlas defined regions, clusters at lower hierarchical levels identified well defined anatomical lobes. At higher hierarchical levels, functional clusters mapped motor and memory functional regions, identifying 50.00%, 20.00%, 27.27%, and 35.14% of the frontal, occipital, parietal, and temporal lobe regions respectively.

Influence of preprocessing, distortion correction and cardiac triggering on the quality of diffusion MR images of spinal cord.

Diffusion MRI of the spinal cord (SC) is susceptible to geometric distortion caused by field inhomogeneities, and prone to misalignment across time series and signal dropout caused by biological motion. Several modifications of image acquisition and image processing techniques have been introduced to overcome these artifacts, but their specific benefits are largely unproven and warrant further investigations. We aim to evaluate two specific aspects of image acquisition and processing that address image quality in diffusion studies of the spinal cord: susceptibility corrections to reduce geometric distortions, and cardiac triggering to minimize motion artifacts. First, we evaluate 4 distortion preprocessing strategies on 7 datasets of the cervical and lumbar SC and find that while distortion correction techniques increase geometric similarity to structural images, they are largely driven by the high-contrast cerebrospinal fluid, and do not consistently improve the geometry within the cord nor improve white-to-gray matter contrast. We recommend at a minimum to perform bulk-motion correction in preprocessing and posit that improvements/adaptations are needed for spinal cord distortion preprocessing algorithms, which are currently optimized and designed for brain imaging. Second, we design experiments to evaluate the impact of removing cardiac triggering. We show that when triggering is foregone, images are qualitatively similar to triggered sequences, do not have increased prevalence of artifacts, and result in similar diffusion tensor indices with similar reproducibility to triggered acquisitions. When triggering is removed, much shorter acquisitions are possible, which are also qualitatively and quantitatively similar to triggered sequences. We suggest that removing cardiac triggering for cervical SC diffusion can be a reasonable option to save time with minimal sacrifice to image quality.

The missing third dimension-Functional correlations of BOLD signals incorporating white matter.

Correlations between magnetic resonance imaging (MRI) blood oxygenation level-dependent (BOLD) signals from pairs of gray matter areas are used to infer their functional connectivity, but they are unable to describe how white matter is engaged in brain networks. Recently, evidence that BOLD signals in white matter are robustly detectable and are modulated by neural activities has accumulated. We introduce a three-way correlation between BOLD signals from pairs of gray matter volumes (nodes) and white matter bundles (edges) to define the communication connectivity through each white matter bundle. Using MRI images from publicly available databases, we show, for example, that the three-way connectivity is influenced by age. By integrating functional MRI signals from white matter as a third component in network analyses, more comprehensive descriptions of brain function may be obtained.

Neuronal microstructural changes in the human brain are associated with neurocognitive aging.

Gray matter (GM) alterations play a role in aging-related disorders like Alzheimer's disease and related dementias, yet MRI studies mainly focus on macroscopic changes. Although reliable indicators of atrophy, morphological metrics like cortical thickness lack the sensitivity to detect early changes preceding visible atrophy. Our study aimed at exploring the potential of diffusion MRI in unveiling sensitive markers of cortical and subcortical age-related microstructural changes and assessing their associations with cognitive and behavioral deficits. We leveraged the Human Connectome Project-Aging cohort that included 707 unimpaired participants (394 female; median age = 58, range = 36-90 years) and applied the powerful mean apparent diffusion propagator model to measure microstructural parameters, along with comprehensive behavioral and cognitive test scores. Both macro- and microstructural GM characteristics were strongly associated with age, with widespread significant microstructural correlations reflective of cellular morphological changes, reduced cellular density, increased extracellular volume, and increased membrane permeability. Importantly, when correlating MRI and cognitive test scores, our findings revealed no link between macrostructural volumetric changes and neurobehavioral performance. However, we found that cellular and extracellular alterations in cortical and subcortical GM regions were associated with neurobehavioral performance. Based on these findings, it is hypothesized that increased microstructural heterogeneity and decreased neurite orientation dispersion precede macrostructural changes, and that they play an important role in subsequent cognitive decline. These alterations are suggested to be early markers of neurocognitive performance that may distinctly aid in identifying the mechanisms underlying phenotypic aging and subsequent age-related functional decline.

Assessing White Matter Engagement in Brain Networks through Functional and Structural Connectivity Mapping.

Understanding the intricate interplay between gray matter (GM) and white matter (WM) is crucial for deciphering the complex activities of the brain. While diffusion tensor imaging (DTI) has advanced the mapping of these structural pathways, the relationship between structural connectivity (SC) and functional connectivity (FC) remains inadequately understood. This study addresses the need for a more integrative approach by mapping the importance of the inter-GM functional link to its structural counterparts in WM. This mapping yields a spatial distribution of engagement that is not only highly reproducible but also aligns with direct structural, functional, and bioenergetic measures within WM, illustrating a notable interdependence between the function of GM and the characteristics of WM. Additionally, our research has uncovered a set of unique engagement modes through a clustering analysis of window-wise engagement maps, highlighting the dyanmic nature of the engagement. The engagement along with their temporal variations revealed significant differences across genders and age groups. These findings suggest the potential of WM engagement as a biomarker for neurological and cognitive conditions, offering a more nuanced understanding of individualized brain activity and connectivity patterns.

Characterizing patterns of DTI variance in aging brains.

Background: As large analyses merge data across sites, a deeper understanding of variance in statistical assessment across the sources of data becomes critical for valid analyses. Diffusion tensor imaging (DTI) exhibits spatially varying and correlated noise, so care must be taken with distributional assumptions. Purpose: We characterize the role of physiology, subject compliance, and the interaction of subject with the scanner in the understanding of DTI variability, as modeled in spatial variance of derived metrics in homogeneous regions. Methods: We analyze DTI data from 1035 subjects in the Baltimore Longitudinal Study of Aging (BLSA), with ages ranging from 22.4 to 103 years old. For each subject, up to 12 longitudinal sessions were conducted. We assess variance of DTI scalars within regions of interest (ROIs) defined by four segmentation methods and investigate the relationships between the variance and covariates, including baseline age, time from the baseline (referred to as "interval"), motion, sex, and whether it is the first scan or the second scan in the session. Results: Covariate effects are heterogeneous and bilaterally symmetric across ROIs. Inter-session interval is positively related (p ≪ 0.001) to FA variance in the cuneus and occipital gyrus, but negatively (p ≪ 0.001) in the caudate nucleus. Males show significantly (p ≪ 0.001) higher FA variance in the right putamen, thalamus, body of the corpus callosum, and cingulate gyrus. In 62 out of 176 ROIs defined by the Eve type-1 atlas, an increase in motion is associated (p < 0.05) with a decrease in FA variance. Head motion increases during the rescan of DTI (Δµ = 0.045 millimeters per volume). Conclusions: The effects of each covariate on DTI variance, and their relationships across ROIs are complex. Ultimately, we encourage researchers to include estimates of variance when sharing data and consider models of heteroscedasticity in analysis. This work provides a foundation for study planning to account for regional variations in metric variance.

Predicting Age from White Matter Diffusivity with Residual Learning.

Imaging findings inconsistent with those expected at specific chronological age ranges may serve as early indicators of neurological disorders and increased mortality risk. Estimation of chronological age, and deviations from expected results, from structural magnetic resonance imaging (MRI) data has become an important proxy task for developing biomarkers that are sensitive to such deviations. Complementary to structural analysis, diffusion tensor imaging (DTI) has proven effective in identifying age-related microstructural changes within the brain white matter, thereby presenting itself as a promising additional modality for brain age prediction. Although early studies have sought to harness DTI's advantages for age estimation, there is no evidence that the success of this prediction is owed to the unique microstructural and diffusivity features that DTI provides, rather than the macrostructural features that are also available in DTI data. Therefore, we seek to develop white-matter-specific age estimation to capture deviations from normal white matter aging. Specifically, we deliberately disregard the macrostructural information when predicting age from DTI scalar images, using two distinct methods. The first method relies on extracting only microstructural features from regions of interest (ROIs). The second applies 3D residual neural networks (ResNets) to learn features directly from the images, which are non-linearly registered and warped to a template to minimize macrostructural variations. When tested on unseen data, the first method yields mean absolute error (MAE) of 6.11 ± 0.19 years for cognitively normal participants and MAE of 6.62 ± 0.30 years for cognitively impaired participants, while the second method achieves MAE of 4.69 ± 0.23 years for cognitively normal participants and MAE of 4.96 ± 0.28 years for cognitively impaired participants. We find that the ResNet model captures subtler, non-macrostructural features for brain age prediction.

DeepN4: Learning N4ITK Bias Field Correction for T1-weighted Images.

T1-weighted (T1w) MRI has low frequency intensity artifacts due to magnetic field inhomogeneities. Removal of these biases in T1w MRI images is a critical preprocessing step to ensure spatially consistent image interpretation. N4ITK bias field correction, the current state-of-the-art, is implemented in such a way that makes it difficult to port between different pipelines and workflows, thus making it hard to reimplement and reproduce results across local, cloud, and edge platforms. Moreover, N4ITK is opaque to optimization before and after its application, meaning that methodological development must work around the inhomogeneity correction step. Given the importance of bias fields correction in structural preprocessing and flexible implementation, we pursue a deep learning approximation / reinterpretation of the N4ITK bias fields correction to create a method which is portable, flexible, and fully differentiable. In this paper, we trained a deep learning network "DeepN4" on eight independent cohorts from 72 different scanners and age ranges with N4ITK-corrected T1w MRI and bias field for supervision in log space. We found that we can closely approximate N4ITK bias fields correction with naïve networks. We evaluate the peak signal to noise ratio (PSNR) in test dataset against the N4ITK corrected images. The median PSNR of corrected images between N4ITK and DeepN4 was 47.96 dB. In addition, we assess the DeepN4 model on eight additional external datasets and show the generalizability of the approach. This study establishes that incompatible N4ITK preprocessing steps can be closely approximated by naïve deep neural networks, facilitating more flexibility. All code and models are released at https://github.com/MASILab/DeepN4.