Fick's law adapted to cationic diffusion in a semi-infinite solid: application to spinel megacrystals

In alkaline lavas, the chemical zoning of megacrystals of spinel is due to the cationic exchange between the latter and the host lava. The application of Fick's law to cationic diffusion profiles allows to calculate the time these crystals have stayed in the lava. Those which are in a chemical equilibrium were in contact with the lava during 20 to 30 days, whereas megacrystals lacking this equilibrium were in contact only for 3 or 4 days. The duration of the rise of an ultrabasic nodule in the volcanic chimney was calculated by applying Stokes' law.

The application of Fick's law to evaluate the period of stay of mega-olivines in alkaline lavas

The integration of the differential equation of the second law of Fick applied to the diffusion of chemical elements in a semi-infinite solid made it easier to estimate the time of stay of olivine mega-cristals in contact with the host lava The results of this research show the existence of two groups of olivine. The first remained in contact with the magmatic liquid during 19 to 22 days, while the second remained so during only 5 to 9 days. This distinction is correlative to that based on the qualitative observation.

Semi-supervised segmentation of ultrasound images based on patch representation and continuous min cut.

Ultrasound segmentation is a challenging problem due to the inherent speckle and some artifacts like shadows, attenuation and signal dropout. Existing methods need to include strong priors like shape priors or analytical intensity models to succeed in the segmentation. However, such priors tend to limit these methods to a specific target or imaging settings, and they are not always applicable to pathological cases. This work introduces a semi-supervised segmentation framework for ultrasound imaging that alleviates the limitation of fully automatic segmentation, that is, it is applicable to any kind of target and imaging settings. Our methodology uses a graph of image patches to represent the ultrasound image and user-assisted initialization with labels, which acts as soft priors. The segmentation problem is formulated as a continuous minimum cut problem and solved with an efficient optimization algorithm. We validate our segmentation framework on clinical ultrasound imaging (prostate, fetus, and tumors of the liver and eye). We obtain high similarity agreement with the ground truth provided by medical expert delineations in all applications (94% DICE values in average) and the proposed algorithm performs favorably with the literature.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

A musculoskeletal shoulder model based on pseudo-inverse and null-space optimization.

The goal of the present work was assess the feasibility of using a pseudo-inverse and null-space optimization approach in the modeling of the shoulder biomechanics. The method was applied to a simplified musculoskeletal shoulder model. The mechanical system consisted in the arm, and the external forces were the arm weight, 6 scapulo-humeral muscles and the reaction at the glenohumeral joint, which was considered as a spherical joint. The muscle wrapping was considered around the humeral head assumed spherical. The dynamical equations were solved in a Lagrangian approach. The mathematical redundancy of the mechanical system was solved in two steps: a pseudo-inverse optimization to minimize the square of the muscle stress and a null-space optimization to restrict the muscle force to physiological limits. Several movements were simulated. The mathematical and numerical aspects of the constrained redundancy problem were efficiently solved by the proposed method. The prediction of muscle moment arms was consistent with cadaveric measurements and the joint reaction force was consistent with in vivo measurements. This preliminary work demonstrated that the developed algorithm has a great potential for more complex musculoskeletal modeling of the shoulder joint. In particular it could be further applied to a non-spherical joint model, allowing for the natural translation of the humeral head in the glenoid fossa.

Diffusion-weighted MR imaging of the spine at 3-T: feasibility, optimization of b-value and utility to differentiate benign from pathological vertebral compression fractures

Purpose: To evaluate the feasibility, determine the optimal b-value, and assess the utility of 3-T diffusion-weighted MR imaging (DWI) of the spine in differentiating benign from pathologic vertebral compression fractures.Methods and Materials: Twenty patients with 38 vertebral compression fractures (24 benign, 14 pathologic) and 20 controls (total: 23 men, 17 women, mean age 56.2years) were included from December 2010 to May 2011 in this IRB-approved prospective study. MR imaging of the spine was performed on a 3-T unit with T1-w, fat-suppressed T2-w, gadolinium-enhanced fat-suppressed T1-w and zoomed-EPI (2D RF excitation pulse combined with reduced field-of-view single-shot echo-planar readout) diffusion-w (b-values: 0, 300, 500 and 700s/mm2) sequences. Two radiologists independently assessed zoomed-EPI image quality in random order using a 4-point scale: 1=excellent to 4=poor. They subsequently measured apparent diffusion coefficients (ADCs) in normal vertebral bodies and compression fractures, in consensus.Results: Lower b-values correlated with better image quality scores, with significant differences between b=300 (mean±SD=2.6±0.8), b=500 (3.0±0.7) and b=700 (3.6±0.6) (all p<0.001). Mean ADCs of normal vertebral bodies (n=162) were 0.23, 0.17 and 0.11×10-3mm2/s with b=300, 500 and 700s/mm2, respectively. In contrast, mean ADCs were 0.89, 0.70 and 0.59×10-3mm2/s for benign vertebral compression fractures and 0.79, 0.66 and 0.51×10-3mm2/s for pathologic fractures with b=300, 500 and 700s/mm2, respectively. No significant difference was found between ADCs of benign and pathologic fractures.Conclusion: 3-T DWI of the spine is feasible and lower b-values (300s/mm2) are recommended. However, our preliminary results show no advantage of DWI in differentiating benign from pathologic vertebral compression fractures.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Artificial Snow Optimization in Winter Sport Destinations Using a Multi-agent Simulation

This paper presents the Juste-Neige system for predicting the snow height on the ski runs of a resort using a multi-agent simulation software. Its aim is to facilitate snow cover management in order to i) reduce the production cost of artificial snow and to improve the profit margin for the companies managing the ski resorts; and ii) to reduce the water and energy consumption, and thus to reduce the environmental impact, by producing only the snow needed for a good skiing experience. The software provides maps with the predicted snow heights for up to 13 days. On these maps, the areas most exposed to snow erosion are highlighted. The software proceeds in three steps: i) interpolation of snow height measurements with a neural network; ii) local meteorological forecasts for every ski resort; iii) simulation of the impact caused by skiers using a multi-agent system. The software has been evaluated in the Swiss ski resort of Verbier and provides useful predictions.

Optimization of spatial resolution for peripheral magnetic resonance angiography.

RATIONALE AND OBJECTIVES: To determine optimum spatial resolution when imaging peripheral arteries with magnetic resonance angiography (MRA). MATERIALS AND METHODS: Eight vessel diameters ranging from 1.0 to 8.0 mm were simulated in a vascular phantom. A total of 40 three-dimensional flash MRA sequences were acquired with incremental variations of fields of view, matrix size, and slice thickness. The accurately known eight diameters were combined pairwise to generate 22 "exact" degrees of stenosis ranging from 42% to 87%. Then, the diameters were measured in the MRA images by three independent observers and with quantitative angiography (QA) software and used to compute the degrees of stenosis corresponding to the 22 "exact" ones. The accuracy and reproducibility of vessel diameter measurements and stenosis calculations were assessed for vessel size ranging from 6 to 8 mm (iliac artery), 4 to 5 mm (femoro-popliteal arteries), and 1 to 3 mm (infrapopliteal arteries). Maximum pixel dimension and slice thickness to obtain a mean error in stenosis evaluation of less than 10% were determined by linear regression analysis. RESULTS: Mean errors on stenosis quantification were 8.8% +/- 6.3% for 6- to 8-mm vessels, 15.5% +/- 8.2% for 4- to 5-mm vessels, and 18.9% +/- 7.5% for 1- to 3-mm vessels. Mean errors on stenosis calculation were 12.3% +/- 8.2% for observers and 11.4% +/- 15.1% for QA software (P = .0342). To evaluate stenosis with a mean error of less than 10%, maximum pixel surface, the pixel size in the phase direction, and the slice thickness should be less than 1.56 mm2, 1.34 mm, 1.70 mm, respectively (voxel size 2.65 mm3) for 6- to 8-mm vessels; 1.31 mm2, 1.10 mm, 1.34 mm (voxel size 1.76 mm3), for 4- to 5-mm vessels; and 1.17 mm2, 0.90 mm, 0.9 mm (voxel size 1.05 mm3) for 1- to 3-mm vessels. CONCLUSION: Higher spatial resolution than currently used should be selected for imaging peripheral vessels.

Cutting edge: influence of the TCR Vbeta domain on the selection of semi-invariant NKT cells by endogenous ligands.

Invariant Valpha14 (Valpha14i) NKT cells are a murine CD1d-dependent regulatory T cell subset characterized by a Valpha14-Jalpha18 rearrangement and expression of mostly Vbeta8.2 and Vbeta7. Whereas the TCR Vbeta domain influences the binding avidity of the Valpha14i TCR for CD1d-alpha-galactosylceramide complexes, with Vbeta8.2 conferring higher avidity binding than Vbeta7, a possible impact of the TCR Vbeta domain on Valpha14i NKT cell selection by endogenous ligands has not been studied. In this study, we show that thymic selection of Vbeta7(+), but not Vbeta8.2(+), Valpha14i NKT cells is favored in situations where endogenous ligand concentration or TCRalpha-chain avidity are suboptimal. Furthermore, thymic Vbeta7(+) Valpha14i NKT cells were preferentially selected in vitro in response to CD1d-dependent presentation of endogenous ligands or exogenously added self ligand isoglobotrihexosylceramide. Collectively, our data demonstrate that the TCR Vbeta domain influences the selection of Valpha14i NKT cells by endogenous ligands, presumably because Vbeta7 confers higher avidity binding.

Optimization of the piggyBac Transposon Using mRNA and Insulators: Toward a More Reliable Gene Delivery System.

Integrating and expressing stably a transgene into the cellular genome remain major challenges for gene-based therapies and for bioproduction purposes. While transposon vectors mediate efficient transgene integration, expression may be limited by epigenetic silencing, and persistent transposase expression may mediate multiple transposition cycles. Here, we evaluated the delivery of the piggyBac transposase messenger RNA combined with genetically insulated transposons to isolate the transgene from neighboring regulatory elements and stabilize expression. A comparison of piggyBac transposase expression from messenger RNA and DNA vectors was carried out in terms of expression levels, transposition efficiency, transgene expression and genotoxic effects, in order to calibrate and secure the transposition-based delivery system. Messenger RNA reduced the persistence of the transposase to a narrow window, thus decreasing side effects such as superfluous genomic DNA cleavage. Both the CTF/NF1 and the D4Z4 insulators were found to mediate more efficient expression from a few transposition events. We conclude that the use of engineered piggyBac transposase mRNA and insulated transposons offer promising ways of improving the quality of the integration process and sustaining the expression of transposon vectors.

Characterization of tumor antigen specific CD8+ T cells and semi-invariant Vα24/Vβ11 NKT cells in tumor invaded liver

Summary: Detailed knowledge on tumor antigen expression and specific immune cells is required for a rational design of immunotherapy for patients with tumor invaded liver. In this study, we confirmed that Cancer/Testis (CT) tumor-associated antigens are frequently expressed in hepatocellular carcinoma (HCC) and searched for the presence of CD8+ T cells specific for these antigens. In 2/10 HLA-A2+ patients with HCC, we found that MAGE-A10 and/or SSX-2 specific CD8+ T cells naturally responded to the disease, since they were enriched in tumor lesions but not in non-tumoral liver. Isolated T cells specifically and strongly killed tumor cells in vitro, suggesting that these CTL were selected in vivo for high avidity antigen recognition, providing the rational for specific immunotherapy of HCC, based on immunization with CT antigens such as MAGE-Al 0 and SSX-2. Type 1 NKT cells express an invariant TCR α chain (Vα24.1α18, paired with Vβ11 in human) and share a specific reactivity to αGalactosylceramide (αGC) presented by CD1d. These cells can display paradoxical immuno-regulatory properties including strong anti-tumor effects upon αGC administration in murine models. To understand why NKT cells were not sufficiently protective against tumor development in patients with tumor invaded liver, we characterized the diversity of Vα24/Vβ11 NKT cells in healthy donors (HD) and cancer patients: NKT cells from HD and patients were generally diverse in terms of TCR β chain (Vβ11) variability and NKT cells from HD showed a variable recognition of αGC loaded CD 1 d multimers. Vα24/ Vβ11 NKT cells can be divided in 3 populations, the CD4, DN (CD4-/CD8-) and CD8 NKT cell subsets that show distinct ability of cytokine production. In addition, our functional analysis revealed that DN and CD8 subsets displayed a higher cytolytic potential and a weaker IFNγ release than the CD4 NKT cell subset. NKT cell subsets were variably represented in the blood of HD and cancer patients. However, HD with high NKT cell frequencies displayed an enrichment of the DN and CD8 subsets, and few of them were suggestive of an oligoclonal expansion in vivo. Comparable NKT cell frequencies were found between blood, non-tumoral liver and tumor of patients. In contrast, we identified a gradual enrichment of CD4 NKT cells from blood to the liver and to the tumor, together with a decrease of DN and CD8 NKT cell subsets. Most patient derived NKT cells were unresponsive upon αGalactosylceramide stimulation ex vivo; NKT cells from few patients displayed a weak responsiveness with different cytokine polarization. The NKT cell repertoire was thus different in tumor tissue, suggesting that CD4 NKT cells infiltrating tumors may be detrimental for protection against tumors and instead may favour the tumor growth/recurrence as recently reported in mice. Résumé en français scientifique : Afin de développer le traitement des patients porteurs d'une tumeur dans le foie par immunothérapie, de nouvelles connaissances sont requises concernant l'expression d'antigènes par les tumeurs et les cellules immunitaires spécifiques de ces antigènes. Nous avons vérifié que des antigènes associés aux tumeurs, tels que les antigènes « Cancer-Testis » (CT), sont fréquemment exprimés par le carcinome hepatocéllulaire (CHC). La recherche de lymphocytes T CD8+ spécifiques (CTL) de ces antigènes a révélé que des CTL spécifiques de MAGE-A10 et/ou SSX-2 ont répondu naturellement à la tumeur chez 2/10 patients étudiés. Ces cellules étaient présentes dans les lésions tumorales mais pas dans le foie adjacent. De plus, ces CTL ont démontré une activité cytolytique forte et spécifique contre les cellules tumorales in vitro, ce qui suggère que ces CTL ont été sélectionnés pour une haute avidité de reconnaissance de l'antigène in vivo. Ces données fournissent une base pour l'immunothérapie spécifique du CHC, en proposant de cibler les antigènes CT tels que MAGE-A10 ou SSX-2. Les cellules NKT de type 1 ont une chaîne α de TCR qui est invariante (chez l'homme, Vα24Jα18, apparié avec Vβ11) et reconnaissent spécifiquement l'αGalactosylceramide (αGC) présenté par CD1d. Ces cellules ont des propriétés immuno¬régulatrices qui peuvent être parfois contradictoires et leur activation par l'αGC induit une forte protection anti-tumorale chez la souris: Afin de comprendre pourquoi ces cellules ne sont pas assez protectrices contre le développement des tumeurs dans le foie chez l'homme, nous avons étudié la diversité des cellules NKT Vα24/Vβ11 d'individus sains (IS) et de patients cancéreux. Les cellules NKT peuvent être sous-divisées en 3 populations : Les CD4, DN (CD4- /CD8-) ou CDS, qui ont la capacité de produire des cytokines différentes. Nos analyses fonctionnelles ont aussi révélé que les sous-populations DN et CD8 ont un potentiel cytolytique plus élevé et une production d'IFNγ plus faible que la sous-population CD4. Ces sous-populations sont représentées de manière variable dans le sang des IS ou des patients. Cependant, les IS avec un taux élevé de cellules NKT ont un enrichissement des sous- populations DN ou CDS, et certains suggèrent qu'il s'agit d'une expansion oligo-clonale in vivo. Les patients avaient des fréquences comparables de cellules NKT entre le sang, le foie et la tumeur. Par contre, la sous-population CD4 était progressivement enrichie du sang vers le foie et la tumeur, tandis que les sous-populations DN ou CD8 était perdues. La plupart des cellules NKT des patients ne réagissaient pas lors de stimulation avec l'αGC ex vivo et les cellules NKT de quelques patients répondaient faiblement et avec des polarisations de cytokines différentes. Ces données suggèrent que les cellules NKT CD4, prédominantes dans les tumeurs, sont inefficaces pour la lutte anti-tumorale et pourraient même favoriser la croissance ou la récurrence tumorale. Donc, une mobilisation spécifique des cellules NKT CD4 négatives par immunothérapie pourrait favoriser l'immunité contre des tumeurs chez l'homme. Résumé en français pour un large public Au sein des globules blancs, les lymphocytes T expriment un récepteur (le TCR), qui est propre à chacun d'entre eux et leur permet d'accrocher de manière très spécifique une molécule appelée antigène. Ce TCR est employé par les lymphocytes pour inspecter les antigènes associés avec des molécules présentatrices à la surface des autres cellules. Les lymphocytes T CD8 reconnaissent un fragment de protéine (ou peptide), qui est présenté par une des molécules du Complexe Majeur d'Histocompatibilité de classe I et tuent la cellule qui présente ce peptide. Ils sont ainsi bien adaptés pour éliminer les cellules qui présentent un peptide issu d'un virus quand la cellule est infectée. D'autres cellules T CD8 reconnaissent des peptides comme les antigènes CT, qui sont produits anormalement par les cellules cancéreuses. Nous avons confirmé que les antigènes CT sont fréquemment exprimés par le cancer du foie. Nous avons également identifié des cellules T CD8 spécifiques d'antigènes CT dans la tumeur, mais pas dans le foie normal de 2 patients sur 10. Cela signifie que ces lymphocytes peuvent être naturellement activés contre la tumeur et sont capables de la trouver. De plus les lymphocytes issus d'un patient ont démontré une forte sensibilité pour reconnaître l'antigène et tuent spécifiquement les cellules tumorales. Les antigènes CT représentent donc des cibles intéressantes qui pourront être intégrés dans des vaccins thérapeutiques du cancer du foie. De cette manière, les cellules T CD8 du patient lui-même pourront être induites à détruire de manière spécifique les cellules cancéreuses. Un nouveau type de lymphocytes T a été récemment découvert: les lymphocytes NKT. Quand ils reconnaissent un glycolipide présenté par la molécule CD1d, ils sont capables, de manière encore incomprise, d'initier, d'augmenter, ou à l'inverse d'inhiber la défense immunitaire. Ces cellules NKT ont démontré qu'elles jouent un rôle important dans la défense contre les tumeurs et particulièrement dans le foie des souris. Nous avons étudié les cellules NKT de patients atteints d'une tumeur dans le foie, afin de comprendre pourquoi elles ne sont pas assez protectrice chez l'homme. Les lymphocytes NKT peuvent être sous-divisés en 3 populations: Les CD4, les DN (CD4-/CD8-) et les CD8. Ces 3 classes de NKT peuvent produire différents signaux chimiques appelés cytokines. Contrairement aux cellules NKT DN ou CDS, seules les cellules NKT CD4 sont capables de produire des cytokines qui sont défavorables pour la défense anti-tumorale. Par ailleurs nous avons trouvé que les cellules NKT CD4 tuent moins bien les cellules cancéreuses que les cellules NKT DN ou CD8. L'analyse des cellules NKT, fraîchement extraites du sang, du foie et de la tumeur de patients a révélé que les cellules NKT CD4 sont progressivement enrichies du sang vers le foie et la tumeur. La large prédominance des NKT CD4 à l'intérieur des tumeurs suggère que, chez l'homme, ces cellules sont inappropriées pour la lutte anti-tumorale. Par ailleurs, la plupart des cellules NKT de patients n'étaient pas capables de produire des cytokines après stimulation avec un antigène. Cela explique également pourquoi ces cellules ne protègent pas contre les tumeurs dans le foie.

COMMIT: Convex Optimization Modeling for Microstructure Informed Tractography.

Tractography is a class of algorithms aiming at in vivo mapping the major neuronal pathways in the white matter from diffusion magnetic resonance imaging (MRI) data. These techniques offer a powerful tool to noninvasively investigate at the macroscopic scale the architecture of the neuronal connections of the brain. However, unfortunately, the reconstructions recovered with existing tractography algorithms are not really quantitative even though diffusion MRI is a quantitative modality by nature. As a matter of fact, several techniques have been proposed in recent years to estimate, at the voxel level, intrinsic microstructural features of the tissue, such as axonal density and diameter, by using multicompartment models. In this paper, we present a novel framework to reestablish the link between tractography and tissue microstructure. Starting from an input set of candidate fiber-tracts, which are estimated from the data using standard fiber-tracking techniques, we model the diffusion MRI signal in each voxel of the image as a linear combination of the restricted and hindered contributions generated in every location of the brain by these candidate tracts. Then, we seek for the global weight of each of them, i.e., the effective contribution or volume, such that they globally fit the measured signal at best. We demonstrate that these weights can be easily recovered by solving a global convex optimization problem and using efficient algorithms. The effectiveness of our approach has been evaluated both on a realistic phantom with known ground-truth and in vivo brain data. Results clearly demonstrate the benefits of the proposed formulation, opening new perspectives for a more quantitative and biologically plausible assessment of the structural connectivity of the brain.