The Mind of a Mouse.

Large scientific projects in genomics and astronomy are influential not because they answer any single question but because they enable investigation of continuously arising new questions from the same data-rich sources. Advances in automated mapping of the brain's synaptic connections (connectomics) suggest that the complicated circuits underlying brain function are ripe for analysis. We discuss benefits of mapping a mouse brain at the level of synapses.

Biomechanical models and mechanisms of cellular morphogenesis and cerebral cortical expansion and folding.

Morphogenesis of the nervous system involves a highly complex spatio-temporal pattern of physical forces (mainly tension and pressure) acting on cells and tissues that are pliable but have an intricately organized cytoskeletal infrastructure. This review begins by covering basic principles of biomechanics and the core cytoskeletal toolkit used to regulate the shapes of cells and tissues during embryogenesis and neural development. It illustrates how the principle of 'tensegrity' provides a useful conceptual framework for understanding how cells dynamically respond to forces that are generated internally or applied externally. The latter part of the review builds on this foundation in considering the development of mammalian cerebral cortex. The main focus is on cortical expansion and folding - processes that take place over an extended period of prenatal and postnatal development. Cortical expansion and folding are likely to involve many complementary mechanisms, some related to regulating cell proliferation and migration and others related to specific types and patterns of mechanical tension and pressure. Three distinct multi-mechanism models are evaluated in relation to a set of 18 key experimental observations and findings. The Composite Tension Plus (CT+) model is introduced as an updated version of a previous multi-component Differential Expansion Sandwich Plus (DES+) model (Van Essen, 2020); the new CT+ model includes 10 distinct mechanisms and has the greatest explanatory power among published models to date. Much needs to be done in order to validate specific mechanistic components and to assess their relative importance in different species, and important directions for future research are suggested.

Age-related differences in resting-state functional connectivity from childhood to adolescence.

The human brain is active at rest, and spontaneous fluctuations in functional MRI BOLD signals reveal an intrinsic functional architecture. During childhood and adolescence, functional networks undergo varying patterns of maturation, and measures of functional connectivity within and between networks differ as a function of age. However, many aspects of these developmental patterns (e.g. trajectory shape and directionality) remain unresolved. In the present study, we characterised age-related differences in within- and between-network resting-state functional connectivity (rsFC) and integration (i.e. participation coefficient, PC) in a large cross-sectional sample of children and adolescents (n = 628) aged 8-21 years from the Lifespan Human Connectome Project in Development. We found evidence for both linear and non-linear differences in cortical, subcortical, and cerebellar rsFC, as well as integration, that varied by age. Additionally, we found that sex moderated the relationship between age and putamen integration where males displayed significant age-related increases in putamen PC compared with females. Taken together, these results provide evidence for complex, non-linear differences in some brain systems during development.

Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin.

Adolescence is characterized by the maturation of cortical microstructure and connectivity supporting complex cognition and behavior. Axonal myelination influences brain connectivity during development by enhancing neural signaling speed and inhibiting plasticity. However, the maturational timing of cortical myelination during human adolescence remains poorly understood. Here, we take advantage of recent advances in high-resolution cortical T1w/T2w mapping methods, including principled correction of B1+ transmit field effects, using data from the Human Connectome Project in Development (HCP-D; N = 628, ages 8-21). We characterize microstructural changes relevant to myelination by estimating age-related differences in T1w/T2w throughout the cerebral neocortex from childhood to early adulthood. We apply Bayesian spline models and clustering analysis to demonstrate graded variation in age-dependent cortical T1w/T2w differences that are correlated with the sensorimotor-association (S-A) axis of cortical organization reported by others. In sensorimotor areas, T1w/T2w ratio measures start at high levels at early ages, increase at a fast pace, and decelerate at later ages (18-21). In intermediate multimodal areas along the S-A axis, T1w/T2w starts at intermediate levels and increases linearly at an intermediate pace. In transmodal/paralimbic association areas, T1w/T2w starts at low levels and increases linearly at the slowest pace. These data provide evidence for graded variation of the T1w/T2w ratio along the S-A axis that may reflect cortical myelination changes during adolescence underlying the development of complex information processing and psychological functioning. We discuss the implications of these results as well as caveats in interpreting magnetic resonance imaging (MRI)-based estimates of myelination.SIGNIFICANCE STATEMENT Myelin is a lipid membrane that is essential to healthy brain function. Myelin wraps axons to increase neural signaling speed, enabling complex neuronal functioning underlying learning and cognition. Here, we characterize the developmental timing of myelination across the cerebral cortex during adolescence using a noninvasive proxy measure, T1w/T2w mapping. Our results provide new evidence demonstrating graded variation across the cortex in the timing of T1w/T2w changes during adolescence, with rapid T1w/T2w increases in lower-order sensory areas and gradual T1w/T2w increases in higher-order association areas. This spatial pattern of microstructural brain development closely parallels the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny.

A 2020 view of tension-based cortical morphogenesis.

Mechanical tension along the length of axons, dendrites, and glial processes has been proposed as a major contributor to morphogenesis throughout the nervous system [D. C. Van Essen, Nature 385, 313-318 (1997)]. Tension-based morphogenesis (TBM) is a conceptually simple and general hypothesis based on physical forces that help shape all living things. Moreover, if each axon and dendrite strive to shorten while preserving connectivity, aggregate wiring length would remain low. TBM can explain key aspects of how the cerebral and cerebellar cortices remain thin, expand in surface area, and acquire their distinctive folds. This article reviews progress since 1997 relevant to TBM and other candidate morphogenetic mechanisms. At a cellular level, studies of diverse cell types in vitro and in vivo demonstrate that tension plays a major role in many developmental events. At a tissue level, I propose a differential expansion sandwich plus (DES+) revision to the original TBM model for cerebral cortical expansion and folding. It invokes tangential tension and "sulcal zipping" forces along the outer cortical margin as well as tension in the white matter core, together competing against radially biased tension in the cortical gray matter. Evidence for and against the DES+ model is discussed, and experiments are proposed to address key tenets of the DES+ model. For cerebellar cortex, a cerebellar multilayer sandwich (CMS) model is proposed that can account for many distinctive features, including its unique, accordion-like folding in the adult, and experiments are proposed to address its specific tenets.

Empirical transmit field bias correction of T1w/T2w myelin maps.

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Anatomical variability, multi-modal coordinate systems, and precision targeting in the marmoset brain.

Localising accurate brain regions needs careful evaluation in each experimental species due to their individual variability. However, the function and connectivity of brain areas is commonly studied using a single-subject cranial landmark-based stereotactic atlas in animal neuroscience. Here, we address this issue in a small primate, the common marmoset, which is increasingly widely used in systems neuroscience. We developed a non-invasive multi-modal neuroimaging-based targeting pipeline, which accounts for intersubject anatomical variability in cranial and cortical landmarks in marmosets. This methodology allowed creation of multi-modal templates (MarmosetRIKEN20) including head CT and brain MR images, embedded in coordinate systems of anterior and posterior commissures (AC-PC) and CIFTI grayordinates. We found that the horizontal plane of the stereotactic coordinate was significantly rotated in pitch relative to the AC-PC coordinate system (10 degrees, frontal downwards), and had a significant bias and uncertainty due to positioning procedures. We also found that many common cranial and brain landmarks (e.g., bregma, intraparietal sulcus) vary in location across subjects and are substantial relative to average marmoset cortical area dimensions. Combining the neuroimaging-based targeting pipeline with robot-guided surgery enabled proof-of-concept targeting of deep brain structures with an accuracy of 0.2 mm. Altogether, our findings demonstrate substantial intersubject variability in marmoset brain and cranial landmarks, implying that subject-specific neuroimaging-based localization is needed for precision targeting in marmosets. The population-based templates and atlases in grayordinates, created for the first time in marmoset monkeys, should help bridging between macroscale and microscale analyses.

A multi-modal parcellation of human cerebral cortex.

Understanding the amazingly complex human cerebral cortex requires a map (or parcellation) of its major subdivisions, known as cortical areas. Making an accurate areal map has been a century-old objective in neuroscience. Using multi-modal magnetic resonance images from the Human Connectome Project (HCP) and an objective semi-automated neuroanatomical approach, we delineated 180 areas per hemisphere bounded by sharp changes in cortical architecture, function, connectivity, and/or topography in a precisely aligned group average of 210 healthy young adults. We characterized 97 new areas and 83 areas previously reported using post-mortem microscopy or other specialized study-specific approaches. To enable automated delineation and identification of these areas in new HCP subjects and in future studies, we trained a machine-learning classifier to recognize the multi-modal 'fingerprint' of each cortical area. This classifier detected the presence of 96.6% of the cortical areas in new subjects, replicated the group parcellation, and could correctly locate areas in individuals with atypical parcellations. The freely available parcellation and classifier will enable substantially improved neuroanatomical precision for studies of the structural and functional organization of human cerebral cortex and its variation across individuals and in development, aging, and disease.

Cerebral cortical folding, parcellation, and connectivity in humans, nonhuman primates, and mice.

Advances in neuroimaging and neuroanatomy have yielded major insights concerning fundamental principles of cortical organization and evolution, thus speaking to how well different species serve as models for human brain function in health and disease. Here, we focus on cortical folding, parcellation, and connectivity in mice, marmosets, macaques, and humans. Cortical folding patterns vary dramatically across species, and individual variability in cortical folding increases with cortical surface area. Such issues are best analyzed using surface-based approaches that respect the topology of the cortical sheet. Many aspects of cortical organization can be revealed using 1 type of information (modality) at a time, such as maps of cortical myelin content. However, accurate delineation of the entire mosaic of cortical areas requires a multimodal approach using information about function, architecture, connectivity, and topographic organization. Comparisons across the 4 aforementioned species reveal dramatic differences in the total number and arrangement of cortical areas, particularly between rodents and primates. Hemispheric variability and bilateral asymmetry are most pronounced in humans, which we evaluated using a high-quality multimodal parcellation of hundreds of individuals. Asymmetries include modest differences in areal size but not in areal identity. Analyses of cortical connectivity using anatomical tracers reveal highly distributed connectivity and a wide range of connection weights in monkeys and mice; indirect measures using functional MRI suggest a similar pattern in humans. Altogether, a multifaceted but integrated approach to exploring cortical organization in primate and nonprimate species provides complementary advantages and perspectives.

Modelling white matter in gyral blades as a continuous vector field.

Many brain imaging studies aim to measure structural connectivity with diffusion tractography. However, biases in tractography data, particularly near the boundary between white matter and cortical grey matter can limit the accuracy of such studies. When seeding from the white matter, streamlines tend to travel parallel to the convoluted cortical surface, largely avoiding sulcal fundi and terminating preferentially on gyral crowns. When seeding from the cortical grey matter, streamlines generally run near the cortical surface until reaching deep white matter. These so-called "gyral biases" limit the accuracy and effective resolution of cortical structural connectivity profiles estimated by tractography algorithms, and they do not reflect the expected distributions of axonal densities seen in invasive tracer studies or stains of myelinated fibres. We propose an algorithm that concurrently models fibre density and orientation using a divergence-free vector field within gyral blades to encourage an anatomically-justified streamline density distribution along the cortical white/grey-matter boundary while maintaining alignment with the diffusion MRI estimated fibre orientations. Using in vivo data from the Human Connectome Project, we show that this algorithm reduces tractography biases. We compare the structural connectomes to functional connectomes from resting-state fMRI, showing that our model improves cross-modal agreement. Finally, we find that after parcellation the changes in the structural connectome are very minor with slightly improved interhemispheric connections (i.e, more homotopic connectivity) and slightly worse intrahemispheric connections when compared to tracers.

Comparative connectomics of the primate social brain.

Social interaction is thought to provide a selection pressure for human intelligence, yet little is known about its neurobiological basis and evolution throughout the primate lineage. Recent advances in neuroimaging have enabled whole brain investigation of brain structure, function, and connectivity in humans and non-human primates (NHPs), leading to a nascent field of comparative connectomics. However, linking social behavior to brain organization across the primates remains challenging. Here, we review the current understanding of the macroscale neural mechanisms of social behaviors from the viewpoint of system neuroscience. We first demonstrate an association between the number of cortical neurons and the size of social groups across primates, suggesting a link between neural information-processing capacity and social capabilities. Moreover, by capitalizing on recent advances in species-harmonized functional MRI, we demonstrate that portions of the mirror neuron system and default-mode networks, which are thought to be important for representation of the other's actions and sense of self, respectively, exhibit similarities in functional organization in macaque monkeys and humans, suggesting possible homologies. With respect to these two networks, we describe recent developments in the neurobiology of social perception, joint attention, personality and social complexity. Together, the Human Connectome Project (HCP)-style comparative neuroimaging, hyperscanning, behavioral, and other multi-modal investigations are expected to yield important insights into the evolutionary foundations of human social behavior.

The nonhuman primate neuroimaging and neuroanatomy project.

Multi-modal neuroimaging projects such as the Human Connectome Project (HCP) and UK Biobank are advancing our understanding of human brain architecture, function, connectivity, and their variability across individuals using high-quality non-invasive data from many subjects. Such efforts depend upon the accuracy of non-invasive brain imaging measures. However, 'ground truth' validation of connectivity using invasive tracers is not feasible in humans. Studies using nonhuman primates (NHPs) enable comparisons between invasive and non-invasive measures, including exploration of how "functional connectivity" from fMRI and "tractographic connectivity" from diffusion MRI compare with long-distance connections measured using tract tracing. Our NonHuman Primate Neuroimaging & Neuroanatomy Project (NHP_NNP) is an international effort (6 laboratories in 5 countries) to: (i) acquire and analyze high-quality multi-modal brain imaging data of macaque and marmoset monkeys using protocols and methods adapted from the HCP; (ii) acquire quantitative invasive tract-tracing data for cortical and subcortical projections to cortical areas; and (iii) map the distributions of different brain cell types with immunocytochemical stains to better define brain areal boundaries. We are acquiring high-resolution structural, functional, and diffusion MRI data together with behavioral measures from over 100 individual macaques and marmosets in order to generate non-invasive measures of brain architecture such as myelin and cortical thickness maps, as well as functional and diffusion tractography-based connectomes. We are using classical and next-generation anatomical tracers to generate quantitative connectivity maps based on brain-wide counting of labeled cortical and subcortical neurons, providing ground truth measures of connectivity. Advanced statistical modeling techniques address the consistency of both kinds of data across individuals, allowing comparison of tracer-based and non-invasive MRI-based connectivity measures. We aim to develop improved cortical and subcortical areal atlases by combining histological and imaging methods. Finally, we are collecting genetic and sociality-associated behavioral data in all animals in an effort to understand how genetic variation shapes the connectome and behavior.

Minimal specifications for non-human primate MRI: Challenges in standardizing and harmonizing data collection.

Recent methodological advances in MRI have enabled substantial growth in neuroimaging studies of non-human primates (NHPs), while open data-sharing through the PRIME-DE initiative has increased the availability of NHP MRI data and the need for robust multi-subject multi-center analyses. Streamlined acquisition and analysis protocols would accelerate and improve these efforts. However, consensus on minimal standards for data acquisition protocols and analysis pipelines for NHP imaging remains to be established, particularly for multi-center studies. Here, we draw parallels between NHP and human neuroimaging and provide minimal guidelines for harmonizing and standardizing data acquisition. We advocate robust translation of widely used open-access toolkits that are well established for analyzing human data. We also encourage the use of validated, automated pre-processing tools for analyzing NHP data sets. These guidelines aim to refine methodological and analytical strategies for small and large-scale NHP neuroimaging data. This will improve reproducibility of results, and accelerate the convergence between NHP and human neuroimaging strategies which will ultimately benefit fundamental and translational brain science.

The Human Connectome Project: A retrospective.

The Human Connectome Project (HCP) was launched in 2010 as an ambitious effort to accelerate advances in human neuroimaging, particularly for measures of brain connectivity; apply these advances to study a large number of healthy young adults; and freely share the data and tools with the scientific community. NIH awarded grants to two consortia; this retrospective focuses on the "WU-Minn-Ox" HCP consortium centered at Washington University, the University of Minnesota, and University of Oxford. In just over 6 years, the WU-Minn-Ox consortium succeeded in its core objectives by: 1) improving MR scanner hardware, pulse sequence design, and image reconstruction methods, 2) acquiring and analyzing multimodal MRI and MEG data of unprecedented quality together with behavioral measures from more than 1100 HCP participants, and 3) freely sharing the data (via the ConnectomeDB database) and associated analysis and visualization tools. To date, more than 27 Petabytes of data have been shared, and 1538 papers acknowledging HCP data use have been published. The "HCP-style" neuroimaging paradigm has emerged as a set of best-practice strategies for optimizing data acquisition and analysis. This article reviews the history of the HCP, including comments on key events and decisions associated with major project components. We discuss several scientific advances using HCP data, including improved cortical parcellations, analyses of connectivity based on functional and diffusion MRI, and analyses of brain-behavior relationships. We also touch upon our efforts to develop and share a variety of associated data processing and analysis tools along with detailed documentation, tutorials, and an educational course to train the next generation of neuroimagers. We conclude with a look forward at opportunities and challenges facing the human neuroimaging field from the perspective of the HCP consortium.

A Domain-General Cognitive Core Defined in Multimodally Parcellated Human Cortex.

Numerous brain imaging studies identified a domain-general or "multiple-demand" (MD) activation pattern accompanying many tasks and may play a core role in cognitive control. Though this finding is well established, the limited spatial localization provided by traditional imaging methods precluded a consensus regarding the precise anatomy, functional differentiation, and connectivity of the MD system. To address these limitations, we used data from 449 subjects from the Human Connectome Project, with the cortex of each individual parcellated using neurobiologically grounded multimodal MRI features. The conjunction of three cognitive contrasts reveals a core of 10 widely distributed MD parcels per hemisphere that are most strongly activated and functionally interconnected, surrounded by a penumbra of 17 additional areas. Outside cerebral cortex, MD activation is most prominent in the caudate and cerebellum. Comparison with canonical resting-state networks shows MD regions concentrated in the fronto-parietal network but also engaging three other networks. MD activations show modest relative task preferences accompanying strong co-recruitment. With distributed anatomical organization, mosaic functional preferences, and strong interconnectivity, we suggest MD regions are well positioned to integrate and assemble the diverse components of cognitive operations. Our precise delineation of MD regions provides a basis for refined analyses of their functions.

The impact of traditional neuroimaging methods on the spatial localization of cortical areas.

Localizing human brain functions is a long-standing goal in systems neuroscience. Toward this goal, neuroimaging studies have traditionally used volume-based smoothing, registered data to volume-based standard spaces, and reported results relative to volume-based parcellations. A novel 360-area surface-based cortical parcellation was recently generated using multimodal data from the Human Connectome Project, and a volume-based version of this parcellation has frequently been requested for use with traditional volume-based analyses. However, given the major methodological differences between traditional volumetric and Human Connectome Project-style processing, the utility and interpretability of such an altered parcellation must first be established. By starting from automatically generated individual-subject parcellations and processing them with different methodological approaches, we show that traditional processing steps, especially volume-based smoothing and registration, substantially degrade cortical area localization compared with surface-based approaches. We also show that surface-based registration using features closely tied to cortical areas, rather than to folding patterns alone, improves the alignment of areas, and that the benefits of high-resolution acquisitions are largely unexploited by traditional volume-based methods. Quantitatively, we show that the most common version of the traditional approach has spatial localization that is only 35% as good as the best surface-based method as assessed using two objective measures (peak areal probabilities and "captured area fraction" for maximum probability maps). Finally, we demonstrate that substantial challenges exist when attempting to accurately represent volume-based group analysis results on the surface, which has important implications for the interpretability of studies, both past and future, that use these volume-based methods.

Quantitative assessment of prefrontal cortex in humans relative to nonhuman primates.

Humans have the largest cerebral cortex among primates. The question of whether association cortex, particularly prefrontal cortex (PFC), is disproportionately larger in humans compared with nonhuman primates is controversial: Some studies report that human PFC is relatively larger, whereas others report a more uniform PFC scaling. We address this controversy using MRI-derived cortical surfaces of many individual humans, chimpanzees, and macaques. We present two parcellation-based PFC delineations based on cytoarchitecture and function and show that a previously used morphological surrogate (cortex anterior to the genu of the corpus callosum) substantially underestimates PFC extent, especially in humans. We find that the proportion of cortical gray matter occupied by PFC in humans is up to 1.9-fold greater than in macaques and 1.2-fold greater than in chimpanzees. The disparity is even more prominent for the proportion of subcortical white matter underlying the PFC, which is 2.4-fold greater in humans than in macaques and 1.7-fold greater than in chimpanzees.

Dynamic patterns of cortical expansion during folding of the preterm human brain.

During the third trimester of human brain development, the cerebral cortex undergoes dramatic surface expansion and folding. Physical models suggest that relatively rapid growth of the cortical gray matter helps drive this folding, and structural data suggest that growth may vary in both space (by region on the cortical surface) and time. In this study, we propose a unique method to estimate local growth from sequential cortical reconstructions. Using anatomically constrained multimodal surface matching (aMSM), we obtain accurate, physically guided point correspondence between younger and older cortical reconstructions of the same individual. From each pair of surfaces, we calculate continuous, smooth maps of cortical expansion with unprecedented precision. By considering 30 preterm infants scanned two to four times during the period of rapid cortical expansion (28-38 wk postmenstrual age), we observe significant regional differences in growth across the cortical surface that are consistent with the emergence of new folds. Furthermore, these growth patterns shift over the course of development, with noninjured subjects following a highly consistent trajectory. This information provides a detailed picture of dynamic changes in cortical growth, connecting what is known about patterns of development at the microscopic (cellular) and macroscopic (folding) scales. Since our method provides specific growth maps for individual brains, we are also able to detect alterations due to injury. This fully automated surface analysis, based on tools freely available to the brain-mapping community, may also serve as a useful approach for future studies of abnormal growth due to genetic disorders, injury, or other environmental variables.

Towards HCP-Style macaque connectomes: 24-Channel 3T multi-array coil, MRI sequences and preprocessing.

Macaque monkeys are an important animal model where invasive investigations can lead to a better understanding of the cortical organization of primates including humans. However, the tools and methods for noninvasive image acquisition (e.g. MRI RF coils and pulse sequence protocols) and image data preprocessing have lagged behind those developed for humans. To resolve the structural and functional characteristics of the smaller macaque brain, high spatial, temporal, and angular resolutions combined with high signal-to-noise ratio are required to ensure good image quality. To address these challenges, we developed a macaque 24-channel receive coil for 3-T MRI with parallel imaging capabilities. This coil enables adaptation of the Human Connectome Project (HCP) image acquisition protocols to the in-vivo macaque brain. In addition, we adapted HCP preprocessing methods to the macaque brain, including spatial minimal preprocessing of structural, functional MRI (fMRI), and diffusion MRI (dMRI). The coil provides the necessary high signal-to-noise ratio and high efficiency in data acquisition, allowing four- and five-fold accelerations for dMRI and fMRI. Automated FreeSurfer segmentation of cortex, reconstruction of cortical surface, removal of artefacts and nuisance signals in fMRI, and distortion correction of dMRI all performed well, and the overall quality of basic neurobiological measures was comparable with those for the HCP. Analyses of functional connectivity in fMRI revealed high sensitivity as compared with those from publicly shared datasets. Tractography-based connectivity estimates correlated with tracer connectivity similarly to that achieved using ex-vivo dMRI. The resulting HCP-style in vivo macaque MRI data show considerable promise for analyzing cortical architecture and functional and structural connectivity using advanced methods that have previously only been available in studies of the human brain.

Classification of temporal ICA components for separating global noise from fMRI data: Reply to Power.

We respond to a critique of our temporal Independent Components Analysis (ICA) method for separating global noise from global signal in fMRI data that focuses on the signal versus noise classification of several components. While we agree with several of Power's comments, we provide evidence and analysis to rebut his major criticisms and to reassure readers that temporal ICA remains a powerful and promising denoising approach.