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.

The correlation between inbreeding and fitness: does allele size matter?

The availability of highly polymorphic genetic markers, in particular microsatellites, has made it possible to test the effect of inbreeding on fitness in the field and in the absence of pedigree information. It has been suggested that the squared difference in allele size at a locus (d(2)) might be a better indicator of the level of inbreeding than is heterozygosity. Using an elegant new analytical model, Tsitrone et al. now put this idea to the test, and to rest.

Universality of rain event size distributions

We compare rain event size distributions derived from measurements in climatically different regions, which we find to be well approximated by power laws of similar exponents over broad ranges. Differences can be seen in the large-scale cutoffs of the distributions. Event duration distributions suggest that the scale-free aspects are related to the absence of characteristic scales in the meteorological mesoscale.

Active mitral ring for post-surgical remote correction of residual mitral regurgitation on the beating heart.

OBJECTIVES: Residual mitral regurgitation after valve repair worsens patients' clinical outcome. Postimplant adjustable mitral rings potentially address this issue, allowing the reshaping of the annulus on the beating heart under echocardiography control. We developed an original mitral ring allowing valve geometry remodelling after the implantation and designed an animal study to assess device effectiveness in correcting residual mitral regurgitation. METHODS: The device consists of two concentric rings: one internal and flexible, sutured to the mitral annulus and a second external and rigid. A third conic element slides between the two rings, modifying the shape of the flexible ring. This sliding element is remotely activated with a rotating tool. Animal model: in adult swine, under cardio pulmonary bypass and cardiac arrest, we shortened the primary chordae of P2 segment to reproduce Type III regurgitation and implanted the active ring. We used intracardiac ultrasound to assess mitral regurgitation and the efficacy of the active ring to correct it. RESULTS: Severe mitral regurgitation (3+ and 4+) was induced in eight animals, 54 ± 6 kg in weight. Vena contracta width decreased from 0.8 ± 0.2 to 0.1 cm; proximal isovelocity surface area radius decreased from 0.8 ± 0.2 to 0.1 cm and effective regurgitant orifice area decreased from 0.50 ± 0.1 to 0.1 ± 0.1 cm(2). Six animals had a reversal of systolic pulmonary flow that normalized following the activation of the device. All corrections were reversible. CONCLUSIONS: Postimplant adjustable mitral ring corrects severe mitral regurgitation through the reversible modification of the annulus geometry on the beating heart. It addresses the frequent and morbid issue of recurrent mitral valve regurgitation.

On dome-shaped norms of reaction for size-to-age at maturity in fishes

The optimal size-to-age at maturity depends on growth and mortality rates, which vary with environment. Therefore, organisms in spatially or temporaly changing environments should develop adaptative phenotypic plasticity for this trait. Experimental work by Alm (1959) on several fish species shows a dome-shape norm of reaction for size-to-age at maturity: size at maturity is smaller in both fast-growing and slow-growing fishes, than it is in fish with a medium growth rate. Using computer simulations, we show that such a dome-shaped norm of reaction is optimal when assuming a finite life span and a negative relationship between production and survival rates. This latter assumption is supported by empirical data, as well as by physiological and emographic arguments.

Some solutions to obtain very efficient separations in isocratic and gradient modes using small particles size and ultra-high pressure.

The UHPLC strategy which combines sub-2 microm porous particles and ultra-high pressure (>1000 bar) was investigated considering very high resolution criteria in both isocratic and gradient modes, with mobile phase temperatures between 30 and 90 degrees C. In isocratic mode, experimental conditions to reach the maximal efficiency were determined using the kinetic plot representation for DeltaP(max)=1000 bar. It has been first confirmed that the molecular weight of the compounds (MW) was a critical parameter which should be considered in the construction of such curves. With a MW around 1000 g mol(-1), efficiencies as high as 300,000 plates could be theoretically attained using UHPLC at 30 degrees C. By limiting the column length to 450 mm, the maximal plate count was around 100,000. In gradient mode, the longest column does not provide the maximal peak capacity for a given analysis time in UHPLC. This was attributed to the fact that peak capacity is not only related to the plate number but also to column dead time. Therefore, a compromise should be found and a 150 mm column should be preferentially selected for gradient lengths up to 60 min at 30 degrees C, while the columns coupled in series (3x 150 mm) were attractive only for t(grad)>250 min. Compared to 30 degrees C, peak capacities were increased by about 20-30% for a constant gradient length at 90 degrees C and gradient time decreased by 2-fold for an identical peak capacity.

Workers cooperation within the firm: an analysis using small and medium size firms

We investigate the determinants of teamwork and workers cooperation within the firm. Up to now the literature has almost exclusively focused on workers incentives as the main determinants for workers cooperation. We take a broader look at the firm's organizational design and analyze the impact that different aspects of it might have on cooperation. In particular, we consider the way in which the degree of decentralization of decisions and the use of complementary HRM practices (what we call the .rm.s vertical organizational design) can affect workers'collaboration with each other. We test the model's predictions on a unique dataset on Spanish small and medium size firms containing a rich set of variables that allows us to use sensible proxies for workers cooperation. We find that the decentralization of labor decisions (and to a less extent that of task planning) has a positive impact on workers cooperation. Likewise, cooperation is positively correlated to many of the HRM practices that seem to favor workers'interaction the most. We also confirm the previous finding that collaborative efforts respond positively to pay incentives, and particularly, to group or company incentives.

Bone mineral density (BMD), micro-architecture estimation (TBS) and vertebral fracture assessment (VFA) extracted from a single DXA device in combination with clinical risk factors improve significantly the identification of women at high risk of fracture

Introduction: Osteoporosis (OP) is a systemic skeletal disease characterized by a low bone mineral density (BMD) and a micro-architectural (MA) deterioration. Clinical risk factors (CRF) are often used as a MA approximation. MA is yet evaluable in daily practice by the Trabecular Bone Score (TBS) measure. TBS is a novel grey-level texture measurement reflecting bone micro-architecture based on the use of experimental variograms of 2D projection images. TBS is very simple to obtain, by reanalyzing a lumbar DXA-scan. TBS has proven to have diagnosis and prognosis value, partially independent of CRF and BMD. The aim of the OsteoLaus cohort is to combine in daily practice the CRF and the information given by DXA (BMD, TBS and vertebral fracture assessment (VFA)) to better identify women at high fracture risk. Method: The OsteoLaus cohort (1400 women 50 to 80 years living in Lausanne, Switzerland) started in 2010. This study is derived from the cohort COLAUS who started in Lausanne in 2003. The main goals of COLAUS is to obtain information on the epidemiology and genetic determinants of cardiovascular risk in 6700 men and women. CRF for OP, bone ultrasound of the heel, lumbar spine and hip BMD, VFA by DXA and MA evaluation by TBS are recorded in OsteoLaus. Preliminary results are reported. Results: We included 631 women: mean age 67.4±6.7 y, BMI 26.1±4.6, mean lumbar spine BMD 0.943±0.168 (T-score -1.4 SD), TBS 1.271±0.103. As expected, correlation between BMD and site matched TBS is low (r2=0.16). Prevalence of VFx grade 2/3, major OP Fx and all OP Fx is 8.4%, 17.0% and 26.0% respectively. Age- and BMI-adjusted ORs (per SD decrease) are 1.8 (1.2- 2.5), 1.6 (1.2-2.1), 1.3 (1.1-1.6) for BMD for the different categories of fractures and 2.0 (1.4-3.0), 1.9 (1.4-2.5), 1.4 (1.1-1.7) for TBS respectively. Only 32 to 37% of women with OP Fx have a BMD < -2.5 SD or a TBS < 1.200. If we combine a BMD < -2.5 SD or a TBS < 1.200, 54 to 60% of women with an osteoporotic Fx are identified. Conclusion: As in the already published studies, these preliminary results confirm the partial independence between BMD and TBS. More importantly, a combination of TBS subsequent to BMD increases significantly the identification of women with prevalent OP Fx which would have been miss-classified by BMD alone. For the first time we are able to have complementary information about fracture (VFA), density (BMD), micro- and macro architecture (TBS & HAS) from a simple, low ionizing radiation and cheap device: DXA. Such complementary information is very useful for the patient in the daily practice and moreover will likely have an impact on cost effectiveness analysis.

A gamma camera count rate saturation correction method for whole-body planar imaging.

Whole-body (WB) planar imaging has long been one of the staple methods of dosimetry, and its quantification has been formalized by the MIRD Committee in pamphlet no 16. One of the issues not specifically addressed in the formalism occurs when the count rates reaching the detector are sufficiently high to result in camera count saturation. Camera dead-time effects have been extensively studied, but all of the developed correction methods assume static acquisitions. However, during WB planar (sweep) imaging, a variable amount of imaged activity exists in the detector's field of view as a function of time and therefore the camera saturation is time dependent. A new time-dependent algorithm was developed to correct for dead-time effects during WB planar acquisitions that accounts for relative motion between detector heads and imaged object. Static camera dead-time parameters were acquired by imaging decaying activity in a phantom and obtaining a saturation curve. Using these parameters, an iterative algorithm akin to Newton's method was developed, which takes into account the variable count rate seen by the detector as a function of time. The algorithm was tested on simulated data as well as on a whole-body scan of high activity Samarium-153 in an ellipsoid phantom. A complete set of parameters from unsaturated phantom data necessary for count rate to activity conversion was also obtained, including build-up and attenuation coefficients, in order to convert corrected count rate values to activity. The algorithm proved successful in accounting for motion- and time-dependent saturation effects in both the simulated and measured data and converged to any desired degree of precision. The clearance half-life calculated from the ellipsoid phantom data was calculated to be 45.1 h after dead-time correction and 51.4 h with no correction; the physical decay half-life of Samarium-153 is 46.3 h. Accurate WB planar dosimetry of high activities relies on successfully compensating for camera saturation which takes into account the variable activity in the field of view, i.e. time-dependent dead-time effects. The algorithm presented here accomplishes this task.

Sensitivity of predictive species distribution models to change in grain size

Predictive species distribution modelling (SDM) has become an essential tool in biodiversity conservation and management. The choice of grain size (resolution) of environmental layers used in modelling is one important factor that may affect predictions. We applied 10 distinct modelling techniques to presence-only data for 50 species in five different regions, to test whether: (1) a 10-fold coarsening of resolution affects predictive performance of SDMs, and (2) any observed effects are dependent on the type of region, modelling technique, or species considered. Results show that a 10 times change in grain size does not severely affect predictions from species distribution models. The overall trend is towards degradation of model performance, but improvement can also be observed. Changing grain size does not equally affect models across regions, techniques, and species types. The strongest effect is on regions and species types, with tree species in the data sets (regions) with highest locational accuracy being most affected. Changing grain size had little influence on the ranking of techniques: boosted regression trees remain best at both resolutions. The number of occurrences used for model training had an important effect, with larger sample sizes resulting in better models, which tended to be more sensitive to grain. Effect of grain change was only noticeable for models reaching sufficient performance and/or with initial data that have an intrinsic error smaller than the coarser grain size.

Genetic and environmental components of phenotypic variation in immune response and body size of a colonial bird, Delichon urbica (the house martin).

Directional selection for parasite resistance is often intense in highly social host species. Using a partial cross-fostering experiment we studied environmental and genetic variation in immune response and morphology in a highly colonial bird species, the house martin (Delichon urbica). We manipulated intensity of infestation of house martin nests by the haematophagous parasitic house martin bug Oeciacus hirundinis either by spraying nests with a weak pesticide or by inoculating them with 50 bugs. Parasitism significantly affected tarsus length, T cell response, immunoglobulin and leucocyte concentrations. We found evidence of strong environmental effects on nestling body mass, body condition, wing length and tarsus length, and evidence of significant additive genetic variance for wing length and haematocrit. We found significant environmental variance, but no significant additive genetic variance in immune response parameters such as T cell response to the antigenic phytohemagglutinin, immunoglobulins, and relative and absolute numbers of leucocytes. Environmental variances were generally greater than additive genetic variances, and the low heritabilities of phenotypic traits were mainly a consequence of large environmental variances and small additive genetic variances. Hence, highly social bird species such as the house martin, which are subject to intense selection by parasites, have a limited scope for immediate microevolutionary response to selection because of low heritabilities, but also a limited scope for long-term response to selection because evolvability as indicated by small additive genetic coefficients of variation is weak.

Managing the size of the Olympic Games

The last two decades have seen dramatic increases in the size and scope of the Summer Olympic Games. In many ways, London 2012 reached even higher summits than the Beijing Games in 2008. This growth is a major challenge for the Olympics and its future organizers, as it is making the Games increasingly difficult to stage and has greatly reduced the number of cities capable of hosting them. This study shows how various participation and organization indicators have expanded over six Olympiads, from Barcelona 1992 to London 2012, and examines the reasons for this growth. It suggests ways of reducing the size of the Summer Olympic Games in order to make them more manageable and to encourage candidatures from smaller cities and countries.

Larval biometry of Simulium rubrithorax (Diptera: Simuliidae) and size comparison between populations in the states of Minas Gerais and Roraima, Brazil

The number of larval instars of Simulium (Hemicnetha) rubrithorax Lutz (Diptera: Nematocera) was determined using the lateral length of the head capsule. In this study 1,035 larvae, of different sizes, were measured (639 from the state of Roraima and 396 from the state of Minas Gerais). A frequency distribution analysis was carried out on the measurements of the lateral length of the head capsule to determine the number of larval instars. The limits of each instar were defined by the lower frequency of the measurements falling in a range of values, by the presence of the "egg burster" that characterizes the first larval instar, and by the developmental stage of the gill histoblast. The determination of the instar number was tested using a Student's t-test (p < 0.05), the Dyar rule and the Crosby growth rule. The results indicate the existence of 7 larval instars for this species, although this result was not in accordance to the Crosby rule. Last-instar larvae from two widely separated geographical populations (Roraima and Minas Gerais), collected in habitats with different water temperature were compared and no differences (p > 0.05) were observed between them.