A post-processing pipeline to reconstruct the developing fetal brain using low-resolution MRI

Introduction. Development of the fetal brain surfacewith concomitant gyrification is one of the majormaturational processes of the human brain. Firstdelineated by postmortem studies or by ultrasound, MRIhas recently become a powerful tool for studying in vivothe structural correlates of brain maturation. However,the quantitative measurement of fetal brain developmentis a major challenge because of the movement of the fetusinside the amniotic cavity, the poor spatial resolution,the partial volume effect and the changing appearance ofthe developing brain. Today extensive efforts are made todeal with the âeurooepost-acquisitionâeuro reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution (Rousseau, F., 2006;Jiang, S., 2007). We here propose a framework devoted tothe segmentation of the basal ganglia, the gray-whitetissue segmentation, and in turn the 3D corticalreconstruction of the fetal brain. Method. Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences in fetuses aged from 29 to 32gestational weeks (slice thickness 5.4mm, in planespatial resolution 1.09mm). For each fetus, 6 axialvolumes shifted by 1 mm were acquired (about 1 min pervolume). First, each volume is manually segmented toextract fetal brain from surrounding fetal and maternaltissues. Inhomogeneity intensity correction and linearintensity normalization are then performed. A highspatial resolution image of isotropic voxel size of 1.09mm is created for each fetus as previously published byothers (Rousseau, F., 2006). B-splines are used for thescattered data interpolation (Lee, 1997). Then, basalganglia segmentation is performed on this superreconstructed volume using active contour framework witha Level Set implementation (Bach Cuadra, M., 2010). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed (Bach Cuadra, M., 2009). Theresulting white matter image is then binarized andfurther given as an input in the Freesurfer software(http://surfer.nmr.mgh.harvard.edu/) to provide accuratethree-dimensional reconstructions of the fetal brain.Results. High-resolution images of the cerebral fetalbrain, as obtained from the low-resolution acquired MRI,are presented for 4 subjects of age ranging from 29 to 32GA. An example is depicted in Figure 1. Accuracy in theautomated basal ganglia segmentation is compared withmanual segmentation using measurement of Dice similarity(DSI), with values above 0.7 considering to be a verygood agreement. In our sample we observed DSI valuesbetween 0.785 and 0.856. We further show the results ofgray-white matter segmentation overlaid on thehigh-resolution gray-scale images. The results arevisually checked for accuracy using the same principlesas commonly accepted in adult neuroimaging. Preliminary3D cortical reconstructions of the fetal brain are shownin Figure 2. Conclusion. We hereby present a completepipeline for the automated extraction of accuratethree-dimensional cortical surface of the fetal brain.These results are preliminary but promising, with theultimate goal to provide âeurooemovieâeuro of the normal gyraldevelopment. In turn, a precise knowledge of the normalfetal brain development will allow the quantification ofsubtle and early but clinically relevant deviations.Moreover, a precise understanding of the gyraldevelopment process may help to build hypotheses tounderstand the pathogenesis of several neurodevelopmentalconditions in which gyrification have been shown to bealtered (e.g. schizophrenia, autismâeuro¦). References.Rousseau, F. (2006), 'Registration-Based Approach forReconstruction of High-Resolution In Utero Fetal MR Brainimages', IEEE Transactions on Medical Imaging, vol. 13,no. 9, pp. 1072-1081. Jiang, S. (2007), 'MRI of MovingSubjects Using Multislice Snapshot Images With VolumeReconstruction (SVR): Application to Fetal, Neonatal, andAdult Brain Studies', IEEE Transactions on MedicalImaging, vol. 26, no. 7, pp. 967-980. Lee, S. (1997),'Scattered data interpolation with multilevel B-splines',IEEE Transactions on Visualization and Computer Graphics,vol. 3, no. 3, pp. 228-244. Bach Cuadra, M. (2010),'Central and Cortical Gray Mater Segmentation of MagneticResonance Images of the Fetal Brain', ISMRM Conference.Bach Cuadra, M. (2009), 'Brain tissue segmentation offetal MR images', MICCAI.

Regional and cellular gene expression changes in human Huntington's disease brain.

Huntington's disease (HD) pathology is well understood at a histological level but a comprehensive molecular analysis of the effect of the disease in the human brain has not previously been available. To elucidate the molecular phenotype of HD on a genome-wide scale, we compared mRNA profiles from 44 human HD brains with those from 36 unaffected controls using microarray analysis. Four brain regions were analyzed: caudate nucleus, cerebellum, prefrontal association cortex [Brodmann's area 9 (BA9)] and motor cortex [Brodmann's area 4 (BA4)]. The greatest number and magnitude of differentially expressed mRNAs were detected in the caudate nucleus, followed by motor cortex, then cerebellum. Thus, the molecular phenotype of HD generally parallels established neuropathology. Surprisingly, no mRNA changes were detected in prefrontal association cortex, thereby revealing subtleties of pathology not previously disclosed by histological methods. To establish that the observed changes were not simply the result of cell loss, we examined mRNA levels in laser-capture microdissected neurons from Grade 1 HD caudate compared to control. These analyses confirmed changes in expression seen in tissue homogenates; we thus conclude that mRNA changes are not attributable to cell loss alone. These data from bona fide HD brains comprise an important reference for hypotheses related to HD and other neurodegenerative diseases.

Expression of dual angiogenic/neurogenic growth factors in human primary brain tumors.

Brain tumors, benign or malignant, are characterized by a very high degree of vascularization. Recent accumulating evidence suggests that during development the neuronal wiring follows the same routes as the vasculature and that these two systems may share some of the same factors for guidance. Thus, expression of dual angiogenic/neurogenic growth factors was evaluated by in situ hybridization in human primary brain tumors of three different types, i.e., astrocytomas, oligodendrogliomas, and ependymomas, of increasing grades, in relation with the grade and type of the tumor. For this evaluation we selected vascular endothelial growth factor (VEGF-A) and its receptors VEGF-R1 and VEGF-R2 and the neuropilins 1 and 2 (NRP-1 and NRP-2), which have proangiogenic properties, platelet-derived growth factor (PDGF) receptor-beta (PDGF-Rβ), which is required for the functional maturation of blood vessels, the ephrins and their Eph receptors, angiotensinogen (AGT) and thrombospondin-2 (TSP-2), which have potential antiangiogenic properties, and netrin-1 (Net-1), which regulates vascular architecture. We show that the expression of the VEGF-NRP system, PDGF-Rβ, TSP-2, AGT, and Net-1 are differentially regulated, either increased or decreased, in relation with the type and grade of the tumor, whereas regulation of the ephrinB system does not seem to be relevant in these human brain tumors.

Postmortem measurement of human chorionic gonadotropin in vitreous humor and bile.

Postmortem human chorionic gonadotrophin (HCG) blood assay can confirm postmortem diagnosis of pregnancy or document situations in which HCG levels are elevated. In some cases, however, blood sampling is not possible at autopsy. In this study, HCG was quantified by enzyme-linked fluorescent assay (ELFA) in the bile (n = 5), vitreous humor (n = 4), and postmortem blood (n = 4) of five pregnant women. There were no false negatives in the pregnant subjects (n = 5) or false positives in controls (n = 34), enabling this test to be recommended for routine use in forensic contexts in which the detection of elevated HCG levels could be of interest.

Decreased Brain Levels of Vitamin B12 in Aging, Autism and Schizophrenia.

Many studies indicate a crucial role for the vitamin B12 and folate-dependent enzyme methionine synthase (MS) in brain development and function, but vitamin B12 status in the brain across the lifespan has not been previously investigated. Vitamin B12 (cobalamin, Cbl) exists in multiple forms, including methylcobalamin (MeCbl) and adenosylcobalamin (AdoCbl), serving as cofactors for MS and methylmalonylCoA mutase, respectively. We measured levels of five Cbl species in postmortem human frontal cortex of 43 control subjects, from 19 weeks of fetal development through 80 years of age, and 12 autistic and 9 schizophrenic subjects. Total Cbl was significantly lower in older control subjects (> 60 yrs of age), primarily reflecting a >10-fold age-dependent decline in the level of MeCbl. Levels of inactive cyanocobalamin (CNCbl) were remarkably higher in fetal brain samples. In both autistic and schizophrenic subjects MeCbl and AdoCbl levels were more than 3-fold lower than age-matched controls. In autistic subjects lower MeCbl was associated with decreased MS activity and elevated levels of its substrate homocysteine (HCY). Low levels of the antioxidant glutathione (GSH) have been linked to both autism and schizophrenia, and both total Cbl and MeCbl levels were decreased in glutamate-cysteine ligase modulatory subunit knockout (GCLM-KO) mice, which exhibit low GSH levels. Thus our findings reveal a previously unrecognized decrease in brain vitamin B12 status across the lifespan that may reflect an adaptation to increasing antioxidant demand, while accelerated deficits due to GSH deficiency may contribute to neurodevelopmental and neuropsychiatric disorders.

EEG-based functional brain networks: does the network size matter?

Functional connectivity in human brain can be represented as a network using electroencephalography (EEG) signals. These networks--whose nodes can vary from tens to hundreds--are characterized by neurobiologically meaningful graph theory metrics. This study investigates the degree to which various graph metrics depend upon the network size. To this end, EEGs from 32 normal subjects were recorded and functional networks of three different sizes were extracted. A state-space based method was used to calculate cross-correlation matrices between different brain regions. These correlation matrices were used to construct binary adjacency connectomes, which were assessed with regards to a number of graph metrics such as clustering coefficient, modularity, efficiency, economic efficiency, and assortativity. We showed that the estimates of these metrics significantly differ depending on the network size. Larger networks had higher efficiency, higher assortativity and lower modularity compared to those with smaller size and the same density. These findings indicate that the network size should be considered in any comparison of networks across studies.

Distinct expression pattern of microtubule-associated protein-2 in human oligodendrogliomas and glial precursor cells.

Microtubule-associated protein 2 (MAP2), a protein linked to the neuronal cytoskeleton in the mature central nervous system (CNS), has recently been identified in glial precursors indicating a potential role during glial development. In the present study, we systematically analyzed the expression of MAP2 in a series of 237 human neuroepithelial tumors including paraffin-embedded specimens and tumor tissue microarrays from oligodendrogliomas, mixed gliomas, astrocytomas, glioblastomas, ependymomas, as well as dysembryoplastic neuroepithelial tumors (DNT), and central neurocytomas. In addition, MAP2-immunoreactive precursor cells were studied in the developing human brain. Three monoclonal antibodies generated against MAP2A-B or MAP2A-D isoforms were used. Variable immunoreactivity for MAP2 could be observed in all gliomas with the exception of ependymomas. Oligodendrogliomas exhibited a consistently strong and distinct pattern of expression characterized by perinuclear cytoplasmic staining without significant process labeling. Tumor cells with immunoreactive bi- or multi-polar processes were mostly encountered in astroglial neoplasms, whereas the small cell component in neurocytomas and DNT was not labeled. These features render MAP2 immunoreactivity a helpful diagnostic tool for the distinction of oligodendrogliomas and other neuroepithelial neoplasms. RT-PCR, Western blot analysis, and in situ hybridization confirmed the expression of MAP2A-C (including the novel MAP2+ 13 transcript) in both oligodendrogliomas and astrocytomas. Double fluorescent laser scanning microscopy showed that GFAP and MAP2 labeled different tumor cell populations. In embryonic human brains, MAP2-immunoreactive glial precursor cells were identified within the subventricular or intermediate zones. These precursors exhibit morphology closely resembling the immunolabeled neoplastic cells observed in glial tumors. Our findings demonstrate MAP2 expression in astrocytic and oligodendroglial neoplasms. The distinct pattern of immunoreactivity in oligodendrogliomas may be useful as a diagnostic tool. Since MAP2 expression occurs transiently in migrating immature glial cells, our findings are in line with an assumed origin of diffuse gliomas from glial precursors.

Quantitative Measurement of Brain Perfusion with Intravoxel Incoherent Motion MR Imaging.

Purpose: To evaluate the sensitivity of the perfusion parameters derived from Intravoxel Incoherent Motion (IVIM) MR imaging to hypercapnia-induced vasodilatation and hyperoxygenation-induced vasoconstriction in the human brain. Materials and Methods: This study was approved by the local ethics committee and informed consent was obtained from all participants. Images were acquired with a standard pulsed-gradient spin-echo sequence (Stejskal-Tanner) in a clinical 3-T system by using 16 b values ranging from 0 to 900 sec/mm(2). Seven healthy volunteers were examined while they inhaled four different gas mixtures known to modify brain perfusion (pure oxygen, ambient air, 5% CO(2) in ambient air, and 8% CO(2) in ambient air). Diffusion coefficient (D), pseudodiffusion coefficient (D*), perfusion fraction (f), and blood flow-related parameter (fD*) maps were calculated on the basis of the IVIM biexponential model, and the parametric maps were compared among the four different gas mixtures. Paired, one-tailed Student t tests were performed to assess for statistically significant differences. Results: Signal decay curves were biexponential in the brain parenchyma of all volunteers. When compared with inhaled ambient air, the IVIM perfusion parameters D*, f, and fD* increased as the concentration of inhaled CO(2) was increased (for the entire brain, P = .01 for f, D*, and fD* for CO(2) 5%; P = .02 for f, and P = .01 for D* and fD* for CO(2) 8%), and a trend toward a reduction was observed when participants inhaled pure oxygen (although P > .05). D remained globally stable. Conclusion: The IVIM perfusion parameters were reactive to hyperoxygenation-induced vasoconstriction and hypercapnia-induced vasodilatation. Accordingly, IVIM imaging was found to be a valid and promising method to quantify brain perfusion in humans. © RSNA, 2012.

Central and Cortical Gray Mater Segmentation of Magnetic Resonance Images of the Fetal Brain

Motivation. The study of human brain development in itsearly stage is today possible thanks to in vivo fetalmagnetic resonance imaging (MRI) techniques. Aquantitative analysis of fetal cortical surfacerepresents a new approach which can be used as a markerof the cerebral maturation (as gyration) and also forstudying central nervous system pathologies [1]. However,this quantitative approach is a major challenge forseveral reasons. First, movement of the fetus inside theamniotic cavity requires very fast MRI sequences tominimize motion artifacts, resulting in a poor spatialresolution and/or lower SNR. Second, due to the ongoingmyelination and cortical maturation, the appearance ofthe developing brain differs very much from thehomogenous tissue types found in adults. Third, due tolow resolution, fetal MR images considerably suffer ofpartial volume (PV) effect, sometimes in large areas.Today extensive efforts are made to deal with thereconstruction of high resolution 3D fetal volumes[2,3,4] to cope with intra-volume motion and low SNR.However, few studies exist related to the automatedsegmentation of MR fetal imaging. [5] and [6] work on thesegmentation of specific areas of the fetal brain such asposterior fossa, brainstem or germinal matrix. Firstattempt for automated brain tissue segmentation has beenpresented in [7] and in our previous work [8]. Bothmethods apply the Expectation-Maximization Markov RandomField (EM-MRF) framework but contrary to [7] we do notneed from any anatomical atlas prior. Data set &Methods. Prenatal MR imaging was performed with a 1-Tsystem (GE Medical Systems, Milwaukee) using single shotfast spin echo (ssFSE) sequences (TR 7000 ms, TE 180 ms,FOV 40 x 40 cm, slice thickness 5.4mm, in plane spatialresolution 1.09mm). Each fetus has 6 axial volumes(around 15 slices per volume), each of them acquired inabout 1 min. Each volume is shifted by 1 mm with respectto the previous one. Gestational age (GA) ranges from 29to 32 weeks. Mother is under sedation. Each volume ismanually segmented to extract fetal brain fromsurrounding maternal tissues. Then, in-homogeneityintensity correction is performed using [9] and linearintensity normalization is performed to have intensityvalues that range from 0 to 255. Note that due tointra-tissue variability of developing brain someintensity variability still remains. For each fetus, ahigh spatial resolution image of isotropic voxel size of1.09 mm is created applying [2] and using B-splines forthe scattered data interpolation [10] (see Fig. 1). Then,basal ganglia (BS) segmentation is performed on thissuper reconstructed volume. Active contour framework witha Level Set (LS) implementation is used. Our LS follows aslightly different formulation from well-known Chan-Vese[11] formulation. In our case, the LS evolves forcing themean of the inside of the curve to be the mean intensityof basal ganglia. Moreover, we add local spatial priorthrough a probabilistic map created by fitting anellipsoid onto the basal ganglia region. Some userinteraction is needed to set the mean intensity of BG(green dots in Fig. 2) and the initial fitting points forthe probabilistic prior map (blue points in Fig. 2). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed as described in [8]. Results.The case study presented here has 29 weeks of GA. Thehigh resolution reconstructed volume is presented in Fig.1. The steps of BG segmentation are shown in Fig. 2.Overlap in comparison with manual segmentation isquantified by the Dice similarity index (DSI) equal to0.829 (values above 0.7 are considered a very goodagreement). Such BG segmentation has been applied on 3other subjects ranging for 29 to 32 GA and the DSI hasbeen of 0.856, 0.794 and 0.785. Our segmentation of theinner (red and blue contours) and outer cortical surface(green contour) is presented in Fig. 3. Finally, torefine the results we include our WM segmentation in theFreesurfer software [12] and some manual corrections toobtain Fig.4. Discussion. Precise cortical surfaceextraction of fetal brain is needed for quantitativestudies of early human brain development. Our workcombines the well known statistical classificationframework with the active contour segmentation forcentral gray mater extraction. A main advantage of thepresented procedure for fetal brain surface extraction isthat we do not include any spatial prior coming fromanatomical atlases. The results presented here arepreliminary but promising. Our efforts are now in testingsuch approach on a wider range of gestational ages thatwe will include in the final version of this work andstudying as well its generalization to different scannersand different type of MRI sequences. References. [1]Guibaud, Prenatal Diagnosis 29(4) (2009). [2] Rousseau,Acad. Rad. 13(9), 2006, [3] Jiang, IEEE TMI 2007. [4]Warfield IADB, MICCAI 2009. [5] Claude, IEEE Trans. Bio.Eng. 51(4) (2004). [6] Habas, MICCAI (Pt. 1) 2008. [7]Bertelsen, ISMRM 2009 [8] Bach Cuadra, IADB, MICCAI 2009.[9] Styner, IEEE TMI 19(39 (2000). [10] Lee, IEEE Trans.Visual. And Comp. Graph. 3(3), 1997, [11] Chan, IEEETrans. Img. Proc, 10(2), 2001 [12] Freesurfer,http://surfer.nmr.mgh.harvard.edu.

Cytotoxic effects of combination of oxidosqualene cyclase inhibitors with atorvastatin in human cancer cells.

Ten oxidosqualene cyclase inhibitors with high efficacy as cholesterol-lowering agents and of different chemical structure classes were evaluated as potential anticancer agents against human cancer cells from various tissue origins and nontumoral human-brain-derived endothelial cells. Inhibition of cancer cell growth was demonstrated at micromolar concentrations, comparable to the concentrations of statins necessary for antitumor effect. Human glioblastoma cells were among the most sensitive cells. These compounds were also able to decrease the proliferation of angiogenic brain-derived endothelial cells, as a model of tumor-induced neovasculation. Additive effects in human glioblastoma cells were also demonstrated for oxidosqualene cyclase inhibitors in combination with atorvastatin while maintaining selectivity against endothelial cells. Thus, not only statins targeting the 3-hydroxy-3-methylglutaryl coenzyme A reductase but also inhibitors of oxidosqualene cyclase decrease tumor growth, suggesting new therapeutic opportunities of combined anti-cholesterol agents for dual treatment of glioblastoma.

Dependence of brain intravoxel incoherent motion perfusion parameters on the cardiac cycle.

Measurement of microvascular perfusion with Intravoxel Incoherent Motion (IVIM) MRI is gaining interest. Yet, the physiological influences on the IVIM perfusion parameters ("pseudo-diffusion" coefficient D*, perfusion fraction f, and flow related parameter fD*) remain insufficiently characterized. In this article, we hypothesize that D* and fD*, which depend on blood speed, should vary during the cardiac cycle. We extended the IVIM model to include time dependence of D* = D*(t), and demonstrate in the healthy human brain that both parameters D* and fD* are significantly larger during systole than diastole, while the diffusion coefficient D and f do not vary significantly. The results non-invasively demonstrate the pulsatility of the brain's microvasculature.

Distribution of the monocarboxylate transporter MCT2 in human cerebral cortex: an immunohistochemical study

The monocarboxylate transporter MCT2 belongs to a large family of membrane proteins involved in the transport of lactate, pyruvate and ketone bodies. Although its expression in rodent brain has been well documented, the presence of MCT2 in the human brain has been questioned on the basis of low mRNA abundance. In this study, the distribution of the monocarboxylate transporter MCT2 has been investigated in the cortex of normal adult human brain using an immunohistochemical approach. Widespread neuropil staining in all cortical layers was observed by light microscopy. Such a distribution was very similar in three different cortical areas investigated. At the cellular level, the expression of MCT2 could be observed in a large number of neurons, in fibers both in grey and white matter, as well as in some astrocytes, mostly localized in layer I and in the white matter. Double staining experiments combined with confocal microscopy confirmed the neuronal expression but also suggested a preferential postsynaptic localization of synaptic MCT2 expression. A few astrocytes in the grey matter appeared to exhibit MCT2 labelling but at low levels. Electron microscopy revealed strong MCT2 expression at asymmetric synapses in the postsynaptic density and also within the spine head but not in the presynaptic terminal. These data not only demonstrate neuronal MCT2 expression in human, but since a portion of it exhibits a distinct synaptic localization, it further supports a putative role for MCT2 in adjustment of energy supply to levels of activity.

Cerebral extracellular lactate increase is predominantly nonischemic in patients with severe traumatic brain injury.

Growing evidence suggests that endogenous lactate is an important substrate for neurons. This study aimed to examine cerebral lactate metabolism and its relationship with brain perfusion in patients with severe traumatic brain injury (TBI). A prospective cohort of 24 patients with severe TBI monitored with cerebral microdialysis (CMD) and brain tissue oxygen tension (PbtO2) was studied. Brain lactate metabolism was assessed by quantification of elevated CMD lactate samples (>4 mmol/L); these were matched to CMD pyruvate and PbtO2 values and dichotomized as glycolytic (CMD pyruvate >119 μmol/L vs. low pyruvate) and hypoxic (PbtO2 <20 mm Hg vs. nonhypoxic). Using perfusion computed tomography (CT), brain perfusion was categorized as oligemic, normal, or hyperemic, and was compared with CMD and PbtO2 data. Samples with elevated CMD lactate were frequently observed (41±8%), and we found that brain lactate elevations were predominantly associated with glycolysis and normal PbtO2 (73±8%) rather than brain hypoxia (14±6%). Furthermore, glycolytic lactate was always associated with normal or hyperemic brain perfusion, whereas all episodes with hypoxic lactate were associated with diffuse oligemia. Our findings suggest predominant nonischemic cerebral extracellular lactate release after TBI and support the concept that lactate may be used as an energy substrate by the injured human brain.

Anatomical correlates for category-specific naming of objects and actions: A brain stimulation mapping study.

The production of object and action words can be dissociated in aphasics, yet their anatomical correlates have been difficult to distinguish in functional imaging studies. To investigate the extent to which the cortical neural networks underlying object- and action-naming processing overlap, we performed electrostimulation mapping (ESM), which is a neurosurgical mapping technique routinely used to examine language function during brain-tumor resections. Forty-one right-handed patients who had surgery for a brain tumor were asked to perform overt naming of object and action pictures under stimulation. Overall, 73 out of the 633 stimulated cortical sites (11.5%) were associated with stimulation-induced language interferences. These interference sites were very much localized (<1 cm(2) ), and showed substantial variability across individuals in their exact localization. Stimulation interfered with both object and action naming over 44 sites, whereas it specifically interfered with object naming over 19 sites and with action naming over 10 sites. Specific object-naming sites were mainly identified in Broca's area (Brodmann area 44/45) and the temporal cortex, whereas action-naming specific sites were mainly identified in the posterior midfrontal gyrus (Brodmann area 6/9) and Broca's area (P = 0.003 by the Fisher's exact test). The anatomical loci we emphasized are in line with a cortical distinction between objects and actions based on conceptual/semantic features, so the prefrontal/premotor cortex would preferentially support sensorimotor contingencies associated with actions, whereas the temporal cortex would preferentially underpin (functional) properties of objects. Hum Brain Mapp 35:429-443, 2014. © 2012 Wiley Periodicals, Inc.